An atlas of human vector-borne microbe interactions reveals pathogenicity mechanisms.

Vector-borne diseases are a leading cause of death worldwide and pose a substantial unmet medical need. Pathogens binding to host extracellular proteins (the "exoproteome") represents a crucial interface in the etiology of vector-borne disease. Here, we used bacterial selection to elucidate host-microbe interactions in high throughput (BASEHIT)-a technique enabling interrogation of microbial interactions with 3,324 human exoproteins-to profile the interactomes of 82 human-pathogen samples, including 30 strains of arthropod-borne pathogens and 8 strains of related non-vector-borne pathogens. The resulting atlas revealed 1,303 putative interactions, including hundreds of pairings with potential roles in pathogenesis, including cell invasion, tissue colonization, immune evasion, and host sensing. Subsequent functional investigations uncovered that Lyme disease spirochetes recognize epidermal growth factor as an environmental cue of transcriptional regulation and that conserved interactions between intracellular pathogens and thioredoxins facilitate cell invasion. In summary, this interactome atlas provides molecular-level insights into microbial pathogenesis and reveals potential host-directed targets for next-generation therapeutics.

Joint host-pathogen genomic analysis identifies hepatitis B virus mutations associated with human NTCP and HLA class I variation.

Evolutionary changes in the hepatitis B virus (HBV) genome could reflect its adaptation to host-induced selective pressure. Leveraging paired human exome and ultra-deep HBV genome-sequencing data from 567 affected individuals with chronic hepatitis B, we comprehensively searched for the signatures of this evolutionary process by conducting "genome-to-genome" association tests between all human genetic variants and viral mutations. We identified significant associations between an East Asian-specific missense variant in the gene encoding the HBV entry receptor NTCP (rs2296651, NTCP S267F) and mutations within the receptor-binding region of HBV preS1. Through in silico modeling and in vitro preS1-NTCP binding assays, we observed that the associated HBV mutations are in proximity to the NTCP variant when bound and together partially increase binding affinity to NTCP S267F. Furthermore, we identified significant associations between HLA-A variation and viral mutations in HLA-A-restricted T cell epitopes. We used in silico binding prediction tools to evaluate the impact of the associated HBV mutations on HLA presentation and observed that mutations that result in weaker binding affinities to their cognate HLA alleles were enriched. Overall, our results suggest the emergence of HBV escape mutations that might alter the interaction between HBV PreS1 and its cellular receptor NTCP during viral entry into hepatocytes and confirm the role of HLA class I restriction in inducing HBV epitope variations.

The NRC0 gene cluster of sensor and helper NLR immune receptors is functionally conserved across asterid plants.

Nucleotide-binding domain and leucine-rich repeat-containing receptor (NLR) proteins can form complex receptor networks to confer innate immunity. NLR-REQUIRED FOR CELL DEATH (NRCs) are phylogenetically related nodes that function downstream of a massively expanded network of disease resistance proteins that protect against multiple plant pathogens. Here, we used phylogenomic methods to reconstruct the macroevolution of the NRC family. One of the NRCs, termed NRC0, is the only family member shared across asterid plants, leading us to investigate its evolutionary history and genetic organization. In several asterid species, NRC0 is genetically clustered with other NLRs that are phylogenetically related to NRC-dependent disease resistance genes. This prompted us to hypothesize that the ancestral state of the NRC network is an NLR helper-sensor gene cluster that was present early during asterid evolution. We provide support for this hypothesis by demonstrating that NRC0 is essential for the hypersensitive cell death that is induced by its genetically linked sensor NLR partners in four divergent asterid species: tomato (Solanum lycopersicum), wild sweet potato (Ipomoea trifida), coffee (Coffea canephora), and carrot (Daucus carota). In addition, activation of a sensor NLR leads to higher-order complex formation of its genetically linked NRC0, similar to other NRCs. Our findings map out contrasting evolutionary dynamics in the macroevolution of the NRC network over the last 125 million years, from a functionally conserved NLR gene cluster to a massive genetically dispersed network.

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.

An interstrand DNA crosslink glycosylase aids Acinetobacter baumannii pathogenesis.

Maintenance of DNA integrity is essential to all forms of life. DNA damage generated by reaction with genotoxic chemicals results in deleterious mutations, genome instability, and cell death. Pathogenic bacteria encounter several genotoxic agents during infection. In keeping with this, the loss of DNA repair networks results in virulence attenuation in several bacterial species. Interstrand DNA crosslinks (ICLs) are a type of DNA lesion formed by covalent linkage of opposing DNA strands and are particularly toxic as they interfere with replication and transcription. Bacteria have evolved specialized DNA glycosylases that unhook ICLs, thereby initiating their repair. In this study, we describe AlkX, a DNA glycosylase encoded by the multidrug resistant pathogen Acinetobacter baumannii. AlkX exhibits ICL unhooking activity similar to that of its Escherichia coli homolog YcaQ. Interrogation of the in vivo role of AlkX revealed that its loss sensitizes cells to DNA crosslinking and impairs A. baumannii colonization of the lungs and dissemination to distal tissues during pneumonia. These results suggest that AlkX participates in A. baumannii pathogenesis and protects the bacterium from stress conditions encountered in vivo. Consistent with this, we found that acidic pH, an environment encountered during host colonization, results in A. baumannii DNA damage and that alkX is induced by, and contributes to, defense against acidic conditions. Collectively, these studies reveal functions for a recently described class of proteins encoded in a broad range of pathogenic bacterial species.

Trichomonas vaginalis extracellular vesicles up-regulate and directly transfer adherence factors promoting host cell colonization.

Trichomonas vaginalis, a common sexually transmitted parasite that colonizes the human urogenital tract, secretes extracellular vesicles (TvEVs) that are taken up by human cells and are speculated to be taken up by parasites as well. While the crosstalk between TvEVs and human cells has led to insight into host:parasite interactions, roles for TvEVs in infection have largely been one-sided, with little known about the effect of TvEV uptake by T. vaginalis. Approximately 11% of infections are found to be coinfections of multiple T. vaginalis strains. Clinical isolates often differ in their adherence to and cytolysis of host cells, underscoring the importance of understanding the effects of TvEV uptake within the parasite population. To address this question, our lab tested the ability of a less adherent strain of T. vaginalis, G3, to take up fluorescently labeled TvEVs derived from both itself (G3-EVs) and TvEVs from a more adherent strain of the parasite (B7RC2-EVs). Here, we showed that TvEVs generated from the more adherent strain are internalized more efficiently compared to the less adherent strain. Additionally, preincubation of G3 parasites with B7RC2-EVs increases parasite aggregation and adherence to host cells. Transcriptomics revealed that TvEVs up-regulate expression of predicted parasite membrane proteins and identified an adherence factor, heteropolysaccharide binding protein (HPB2). Finally, using comparative proteomics and superresolution microscopy, we demonstrated direct transfer of an adherence factor, cadherin-like protein, from TvEVs to the recipient parasite's surface. This work identifies TvEVs as a mediator of parasite:parasite communication that may impact pathogenesis during mixed infections.

Vibrio MARTX toxin processing and degradation of cellular Rab GTPases by the cytotoxic effector Makes Caterpillars Floppy.

Vibrio vulnificus causes life-threatening wound and gastrointestinal infections, mediated primarily by the production of a Multifunctional-Autoprocessing Repeats-In-Toxin (MARTX) toxin. The most commonly present MARTX effector domain, the Makes Caterpillars Floppy-like (MCF) toxin, is a cysteine protease stimulated by host adenosine diphosphate (ADP) ribosylation factors (ARFs) to autoprocess. Here, we show processed MCF then binds and cleaves host Ras-related proteins in brain (Rab) guanosine triphosphatases within their C-terminal tails resulting in Rab degradation. We demonstrate MCF binds Rabs at the same interface occupied by ARFs. Moreover, we show MCF preferentially binds to ARF1 prior to autoprocessing and is active to cleave Rabs only subsequent to autoprocessing. We then use structure prediction algorithms to demonstrate that structural composition, rather than sequence, determines Rab target specificity. We further determine a crystal structure of aMCF as a swapped dimer, revealing an alternative conformation we suggest represents the open, activated state of MCF with reorganized active site residues. The cleavage of Rabs results in Rab1B dispersal within cells and loss of Rab1B density in the intestinal tissue of infected mice. Collectively, our work describes an extracellular bacterial mechanism whereby MCF is activated by ARFs and subsequently induces the degradation of another small host guanosine triphosphatase (GTPase), Rabs, to drive organelle damage, cell death, and promote pathogenesis of these rapidly fatal infections.

The germ theory revisited: A noncentric view on infection outcome.

The germ theory states that pathogenic microorganisms are responsible for causing infectious diseases. The theory is inherently microbe-centric and does not account for variability in disease severity among individuals and asymptomatic carriership-two phenomena indicating an important role for host variability in infection outcome. The basic tenet of the germ theory was recently challenged, and a radically host-centric paradigm referred to as the "full-blown host theory" was proposed. According to this view, the pathogen is reduced to a passive environmental trigger, and the development of disease is instead due to pre-existing immunodeficiencies of the host. Here, we consider the factors that determine disease severity using established knowledge concerning evolutionary biology, microbial pathogenesis, and host-pathogen interactions. We note that the available data support a noncentric view that recognizes key roles for both the causative microbe and the host in dictating infection outcome.

Mitochondrial recovery by the UPRmt: Insights from C. elegans.

Mitochondria are multifaceted organelles, with such functions as the production of cellular energy to the regulation of cell death. However, mitochondria incur various sources of damage from the accumulation of reactive oxygen species and DNA mutations that can impact the protein folding environment and impair their function. Since mitochondrial dysfunction is often associated with reductions in organismal fitness and possibly disease, cells must have safeguards in place to protect mitochondrial function and promote recovery during times of stress. The mitochondrial unfolded protein response (UPRmt) is a transcriptional adaptation that promotes mitochondrial repair to aid in cell survival during stress. While the earlier discoveries into the regulation of the UPRmt stemmed from studies using mammalian cell culture, much of our understanding about this stress response has been bestowed to us by the model organism Caenorhabditis elegans. Indeed, the facile but powerful genetics of this relatively simple nematode has uncovered multiple regulators of the UPRmt, as well as several physiological roles of this stress response. In this review, we will summarize these major advancements originating from studies using C. elegans.

Host caveolin-1 facilitates Zika virus infection by promoting viral RNA replication.

Zika virus (ZIKV) has gained notoriety in recent years because there are no targeted therapies or vaccines available so far. Caveolin-1 (Cav-1) in host cells plays crucial functions in the invasion of many viruses. However, its specific involvement in ZIKV infection has remained unclear. Here, we reveal that depleting Cav-1 leads to a substantial reduction in ZIKV RNA levels, protein expression and viral particle production, indicating that ZIKV exploits Cav-1 for its infection. By dissecting each stage of the viral life cycle, we unveil that, unlike its invasion role in many other viruses, Cav-1 depletion selectively impairs ZIKV replication, resulting in altered replication dynamics and reduced strand-specific RNA levels, but does not affect viral entry, maturation and release. These results reveal an unforeseen function of Cav-1 in facilitating ZIKV replication, which provides new insights into the intricate interaction between Cav-1 and ZIKV and underscores Cav-1 as a potential candidate for anti-ZIKV approaches.

Quantification of Adaptive Immune Responses Against Protein-Binding Interfaces in the Streptococcal M1 Protein.

Bacterial or viral antigens can contain subdominant protein regions that elicit weak antibody responses upon vaccination or infection although there is accumulating evidence that antibody responses against subdominant regions can enhance the protective immune response. One proposed mechanism for subdominant protein regions is the binding of host proteins that prevent antibody production against epitopes hidden within the protein binding interfaces. Here, we used affinity purification combined with quantitative mass spectrometry (AP-MS) to examine the level of competition between antigen-specific antibodies and host-pathogen protein interaction networks using the M1 protein from Streptococcus pyogenes as a model system. As most humans have circulating antibodies against the M1 protein, we first used AP-MS to show that the M1 protein interspecies protein network formed with human plasma proteins is largely conserved in naïve mice. Immunizing mice with the M1 protein generated a time-dependent increase of anti-M1 antibodies. AP-MS analysis comparing the composition of the M1-plasma protein network from naïve and immunized mice showed significant enrichment of 292 IgG peptides associated with 56 IgG chains in the immune mice. Despite the significant increase of bound IgGs, the levels of interacting plasma proteins were not significantly reduced in the immune mice. The results indicate that the antigen-specific polyclonal IgG against the M1 protein primarily targets epitopes outside the other plasma protein binding interfaces. In conclusion, this study demonstrates that AP-MS is a promising strategy to determine the relationship between antigen-specific antibodies and host-pathogen interaction networks that could be used to define subdominant protein regions of relevance for vaccine development.

Scs system links copper and redox homeostasis in bacterial pathogens.

The bacterial envelope is an essential compartment involved in metabolism and metabolites transport, virulence, and stress defense. Its roles become more evident when homeostasis is challenged during host-pathogen interactions. In particular, the presence of free radical groups and excess copper in the periplasm causes noxious reactions, such as sulfhydryl group oxidation leading to enzymatic inactivation and protein denaturation. In response to this, canonical and accessory oxidoreductase systems are induced, performing quality control of thiol groups, and therefore contributing to restoring homeostasis and preserving survival under these conditions. Here, we examine recent advances in the characterization of the Dsb-like, Salmonella-specific Scs system. This system includes the ScsC/ScsB pair of Cu+-binding proteins with thiol-oxidoreductase activity, an alternative ScsB-partner, the membrane-linked ScsD, and a likely associated protein, ScsA, with a role in peroxide resistance. We discuss the acquisition of the scsABCD locus and its integration into a global regulatory pathway directing envelope response to Cu stress during the evolution of pathogens that also harbor the canonical Dsb systems. The evidence suggests that the canonical Dsb systems cannot satisfy the extra demands that the host-pathogen interface imposes to preserve functional thiol groups. This resulted in the acquisition of the Scs system by Salmonella. We propose that the ScsABCD complex evolved to connect Cu and redox stress responses in this pathogen as well as in other bacterial pathogens.

Proximity labeling of host factor ANXA3 in HCV infection reveals a novel LARP1 function in viral entry.

Hepatitis C virus (HCV) infection is tightly connected to the lipid metabolism with lipid droplets (LDs) serving as assembly sites for progeny virions. A previous LD proteome analysis identified annexin A3 (ANXA3) as an important HCV host factor that is enriched at LDs in infected cells and required for HCV morphogenesis. To further characterize ANXA3 function in HCV, we performed proximity labeling using ANXA3-BioID2 as bait in HCV-infected cells. Two of the top proteins identified proximal to ANXA3 during HCV infection were the La-related protein 1 (LARP1) and the ADP ribosylation factor-like protein 8B (ARL8B), both of which have been previously described to act in HCV particle production. In follow-up experiments, ARL8B functioned as a pro-viral HCV host factor without localizing to LDs and thus likely independent of ANXA3. In contrast, LARP1 interacts with HCV core protein in an RNA-dependent manner and is translocated to LDs by core protein. Knockdown of LARP1 decreased HCV spreading without altering HCV RNA replication or viral titers. Unexpectedly, entry of HCV particles and E1/E2-pseudotyped lentiviral particles was reduced by LARP1 depletion, whereas particle production was not altered. Using a recombinant vesicular stomatitis virus (VSV)ΔG entry assay, we showed that LARP1 depletion also decreased entry of VSV with VSV, MERS, and CHIKV glycoproteins. Therefore, our data expand the role of LARP1 as an HCV host factor that is most prominently involved in the early steps of infection, likely contributing to endocytosis of viral particles through the pleiotropic effect LARP1 has on the cellular translatome.

Crystal structure of bacterial ubiquitin ADP-ribosyltransferase CteC reveals a substrate-recruiting insertion.

ADP-ribosylation is a post-translational modification involved in regulation of diverse cellular pathways. Interestingly, many pathogens have been identified to utilize ADP-ribosylation as a way for host manipulation. A recent study found that CteC, an effector from the bacterial pathogen Chromobacterium violaceum, hinders host ubiquitin (Ub) signaling pathways via installing mono-ADP-ribosylation on threonine 66 of Ub. However, the molecular basis of substrate recognition by CteC is not well understood. In this article, we probed the substrate specificity of this effector at protein and residue levels. We also determined the crystal structure of CteC in complex with NAD+, which revealed a canonical mono-ADP-ribosyltransferase fold with an additional insertion domain. The AlphaFold-predicted model differed significantly from the experimentally determined structure, even in regions not used in crystal packing. Biochemical and biophysical studies indicated unique features of the NAD+ binding pocket, while showing selectivity distinction between Ub and structurally close Ub-like modifiers and the role of the insertion domain in substrate recognition. Together, this study provides insights into the enzymatic specificities and the key structural features of a novel bacterial ADP-ribosyltransferase involved in host-pathogen interaction.

Epstein-Barr virus nuclear antigen EBNA3A modulates IRF3-dependent IFNß expression.

Epstein-Barr virus (EBV), the causative agent of infectious mononucleosis, persistently infects over 90% of the human adult population and is associated with several human cancers. To establish life-long infection, EBV tampers with the induction of type I interferon (IFN I)-dependent antiviral immunity in the host. How various EBV genes help orchestrate this crucial strategy is incompletely defined. Here, we reveal a mechanism by which the EBV nuclear antigen 3A (EBNA3A) may inhibit IFNß induction. Using proximity biotinylation we identify the histone acetyltransferase P300, a member of the IFNß transcriptional complex, as a binding partner of EBNA3A. We further show that EBNA3A also interacts with the activated IFN-inducing transcription factor IRF3 that collaborates with P300 in the nucleus. Both events are mediated by the N-terminal domain of EBNA3A. We propose that EBNA3A limits binding of IRF3 to the IFNß promoter, thereby hampering downstream IFN I signaling. Collectively, our findings suggest a new mechanism of immune evasion by EBV, affected by its latency gene EBNA3A.

Toxoplasma gondii mitochondrial association factor 1b interactome reveals novel binding partners including Ral GTPase accelerating protein α1.

The intracellular parasite, Toxoplasma gondii, has developed sophisticated molecular strategies to subvert host processes and promote growth and survival. During infection, T. gondii replicates in a parasitophorous vacuole (PV) and modulates host functions through a network of secreted proteins. Of these, Mitochondrial Association Factor 1b (MAF1b) recruits host mitochondria to the PV, a process that confers an in vivo growth advantage, though the precise mechanisms remain enigmatic. To address this knowledge gap, we mapped the MAF1b interactome in human fibroblasts using a commercial Yeast-2-hybrid (Y2H) screen, which revealed several previously unidentified binding partners including the GAP domain of Ral GTPase Accelerating Protein α1 (RalGAPα1(GAP)). Recombinantly produced MAF1b and RalGAPα1(GAP) formed as a stable binary complex as shown by size exclusion chromatography with a Kd of 334 nM as measured by isothermal titration calorimetry (ITC). Notably, no binding was detected between RalGAPα1(GAP) and the structurally conserved MAF1b homolog, MAF1a, which does not recruit host mitochondria. Next, we used hydrogen deuterium exchange mass spectrometry (HDX-MS) to map the RalGAPα1(GAP)-MAF1b interface, which led to identification of the "GAP-binding loop" on MAF1b that was confirmed by mutagenesis and ITC to be necessary for complex formation. A high-confidence Alphafold model predicts the GAP-binding loop to lie at the RalGAPα1(GAP)-MAF1b interface further supporting the HDX-MS data. Mechanistic implications of a RalGAPα1(GAP)-MAF1b complex are discussed in the context of T. gondii infection and indicates that MAF1b may have evolved multiple independent functions to increase T. gondii fitness.

Structural and functional characterization of the IpaD π-helix reveals critical roles in DOC interaction, T3SS apparatus maturation, and Shigella virulence.

Shigella spp. are highly pathogenic members of the Enterobacteriaceae family, causing ∼269 million cases of bacillary dysentery and >200,000 deaths each year. Like many Gram-negative pathogens, Shigella rely on their type three secretion system (T3SS) to inject effector proteins into eukaryotic host cells, driving both cellular invasion and evasion of host immune responses. Exposure to the bile salt deoxycholate (DOC) significantly enhances Shigella virulence and is proposed to serve as a critical environmental signal present in the small intestine that prepares Shigella's T3SS for efficient infection of the colonic epithelium. Here, we uncover critical mechanistic details of the Shigella-specific DOC signaling process by describing the role of a π-helix secondary structure element within the T3SS tip protein IpaD. Biophysical characterization and high-resolution structures of IpaD mutants lacking the π-helix show that it is not required for global protein structure, but that it defines the native DOC binding site and prevents off target interactions. Additionally, Shigella strains expressing the π-helix deletion mutants illustrate the pathogenic importance of its role in guiding DOC interaction as flow cytometry and gentamycin protection assays show that the IpaD π-helix is essential for DOC-mediated apparatus maturation and enhanced invasion of eukaryotic cells. Together, these findings add to our understanding of the complex Shigella pathogenesis pathway and its evolution to respond to environmental bile salts by identifying the π-helix in IpaD as a critical structural element required for translating DOC exposure to virulence enhancement.

Suppression of plant immunity by Verticillium dahliae effector Vd6317 through AtNAC53 association.

Verticillium dahliae, a soil-borne fungal pathogen, compromises host innate immunity by secreting a plethora of effectors, thereby facilitating host colonization and causing substantial yield and quality losses. The mechanisms underlying the modulation of cotton immunity by V. dahliae effectors are predominantly unexplored. In this study, we identified that the V. dahliae effector Vd6317 inhibits plant cell death triggered by Vd424Y and enhances PVX viral infection in Nicotiana benthamiana. Attenuation of Vd6317 significantly decreased the virulence of V. dahliae, whereas ectopic expression of Vd6317 in Arabidopsis and cotton enhanced susceptibility to V. dahliae infection, underscoring Vd6317's critical role in pathogenicity. We observed that Vd6317 targeted the Arabidopsis immune regulator AtNAC53, thereby impeding its transcriptional activity on the defense-associated gene AtUGT74E2. Arabidopsis nac53 and ugt74e2 mutants exhibited heightened sensitivity to V. dahliae compared to wild-type plants. A mutation at the conserved residue 193L of Vd6317 abrogated its interaction with AtNAC53 and reduced the virulence of V. dahliae, which was partially attributable to a reduction in Vd6317 protein stability. Our findings unveil a hitherto unrecognized regulatory mechanism by which the V. dahliae effector Vd6317 directly inhibits the plant transcription factor AtNAC53 activity to suppress the expression of AtUGT74E2 and plant defense.

Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat.

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

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 x genotype (GxG) interactions at the genetic level determining the outcome and severity of infection. While co-analyses of host and pathogen genomes (co-GWAs) 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 GxG interactions and to infer their underlying interaction matrix. In a combined SNP data set 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 (GxG) 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 GxG interactions with host and pathogen genomic data. It, therefore, paves the way for understanding the evolution of key GxG interactions underpinning HCV adaptation to the European human population after a recent expansion.