Search

Since the identification of sickle cell trait as a heritable form of resistance to malaria, candidate gene studies, linkage analysis paired with sequencing, and genome-wide association (GWA) studies have revealed many examples of genetic resistance and susceptibility to infectious diseases. GWA studies enabled the identification of many common variants associated with small shifts in susceptibility to infectious diseases. This is exemplified by multiple loci associated with leprosy, malaria, HIV, tuberculosis, and coronavirus disease 2019 (COVID-19), which illuminate genetic architecture and implicate pathways underlying pathophysiology. Despite these successes, most of the heritability of infectious diseases remains to be explained. As the field advances, current limitations may be overcome by applying methodological innovations such as cellular GWA studies and phenome-wide association (PheWA) studies as well as by improving methodological rigor with more precise case definitions, deeper phenotyping, increased cohort diversity, and functional validation of candidate loci in the laboratory or human challenge studies.

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human in vitro susceptibility testing (scHi-HOST), a method for rapidly identifying genetic variants that confer resistance and susceptibility. We applied this method to influenza A virus (IAV), the cause of four pandemics since the start of the 20th century. scHi-HOST leverages single-cell RNA sequencing (scRNA-seq) to simultaneously assign genetic identity to cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral burden, and identify expression quantitative trait loci. Integration of scHi-HOST with human challenge and experimental validation demonstrated that a missense variant in endoplasmic reticulum aminopeptidase 1 (ERAP1; rs27895) increased IAV burden in cells and human volunteers. rs27895 exhibits population differentiation, likely contributing to greater permissivity of cells from African populations to IAV. scHi-HOST is a broadly applicable method and resource for decoding infectious-disease genetics.

Despite being in a golden age of bacterial epigenomics, little work has systematically examined the plasticity and functional impacts of the bacterial DNA methylome. Here, we leveraged single-molecule, real-time sequencing (SMRT-seq) to examine the m

Human genetic diversity can reveal critical factors in host-pathogen interactions. This is especially useful for human-restricted pathogens like Salmonella enterica serovar Typhi (S. Typhi), the cause of Typhoid fever. One key dynamic during infection is competition for nutrients: host cells attempt to restrict intracellular replication by depriving bacteria of key nutrients or delivering toxic metabolites in a process called nutritional immunity. Here, a cellular genome-wide association study of intracellular replication by S. Typhi in nearly a thousand cell lines from around the world—and extensive follow-up using intracellular S. Typhi transcriptomics and manipulation of magnesium concentrations—demonstrates that the divalent cation channel mucolipin-2 (MCOLN2) restricts S. Typhi intracellular replication through magnesium deprivation. Our results reveal natural diversity in Mg2+ limitation as a key component of nutritional immunity against S. Typhi.One-Sentence SummaryHuman immune cells genetically vary in their ability to use magnesium deprivation to restrict growth of the typhoid fever bacterium.

Humanity's ongoing struggle with new, re-emerging and endemic infectious diseases serves as a frequent reminder of the need to understand host-pathogen interactions. Recent advances in genomics have dramatically advanced our understanding of how genetics contributes to host resistance or susceptibility to bacterial infection. Here we discuss current trends in defining host-bacterial interactions at the genome-wide level, including screens that harness CRISPR/Cas9 genome editing, natural genetic variation, proteomics, and transcriptomics. We report on the merits, limitations, and findings of these innovative screens and discuss their complementary nature. Finally, we speculate on future innovation as we continue to progress through the postgenomic era and towards deeper mechanistic insight and clinical applications.

SummaryDiversity in the human genome is one factor that confers resistance and susceptibility to infectious diseases. This is observed most dramatically during pandemics, where individuals exhibit large differences in risk and clinical outcomes against a pathogen infecting large portions of the world’s populations. Here, we developed scHi-HOST (single cell High-throughput Human in vitrO Susceptibility Testing), a method for rapidly identifying genetic variants that confer resistance and susceptibility to pathogens. scHi-HOST leverages scRNA-seq (single-cell RNA-sequencing) to simultaneously assign genetic identity to individual cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral entry and replication, and identify expression quantitative trait loci (eQTLs). Applying scHi-HOST to influenza A virus (IAV), we identified eQTLs at baseline and in genes that are induced by IAV infection. Integration of scHi-HOST with a human IAV challenge study (Prometheus) revealed that a missense variant inERAP1(Endoplasmic reticulum aminopeptidase 1;rs27895) was associated with IAV burden in cells and human volunteers. Functional studies using RNA interference, ERAP1 inhibitor, and overexpression of alternative alleles demonstrated that ERAP1 is exploited by IAV to promote infection. Specifically, the nonsynonymous substitution, which results in a glycine to aspartate substitution at ERAP1 residue 348, would disrupt the substrate binding pocket of ERAP1, likely resulting in a significantly altered preference for substrates, poorer catalytic efficiency, or both. Finally, rs27895 exhibits substantial population differentiation, with the higher frequency of the minor T allele in two African populations likely contributing to the greater permissivity of cells from these populations to IAV infection. scHi-HOST is an important resource for understanding susceptibility to influenza and is a broadly applicable method for decoding human genetics of infectious disease.

SummaryGenes associated with increased susceptibility to multiple sclerosis (MS) have been identified, but their functions are incompletely understood. One of these genes codes for the RNA helicase DExD/H-Box Polypeptide 39B (DDX39B), which shows genetic and functional epistasis with interleukin-7 receptor-α gene (IL7R) in MS-risk. Based on evolutionary and functional arguments we postulated that DDX39B enhances immune tolerance decreasing MS risk. Consistent with such a role we show that DDX39B controls the expression of many MS susceptibility genes and important immune-related genes. Among these we identified ForkheadBox P3 (FOXP3), which codes for the master transcriptional factor in CD4+/CD25+ T regulatory cells. DDX39B knockdown led to loss of immune-regulatory and gain of immune-effector expression signatures. Splicing of FOXP3 introns, which belong to a previously unrecognized subclass of introns with C-rich polypyrimidine tracts, was exquisitely sensitive to DDX39B levels. Given the importance of FOXP3 in autoimmunity, this work cements DDX39B as an important guardian of immune tolerance.

Summary Human genetic diversity can have profound effects on health outcomes upon exposure to infectious agents. For infections with Chlamydia trachomatis (C. trachomatis), the wide range of genital and ocular disease manifestations are likely influenced by human genetic differences that regulate interactions between C. trachomatis and host cells. We leveraged this diversity in cellular responses to demonstrate the importance of variation at the Toll-like receptor 1 (TLR1), TLR6, and TLR10 locus to cytokine production in response to C. trachomatis. We determined that a single-nucleotide polymorphism (SNP) (rs1057807), located in a region that forms a loop with the TLR6 promoter, is associated with increased expression of TLR1, TLR6, and TLR10 and secreted levels of ten C. trachomatis-induced cytokines. Production of these C. trachomatis-induced cytokines is primarily dependent on MyD88 and TLR6 based on experiments using inhibitors, blocking antibodies, RNAi, and protein overexpression. Population genetic analyses further demonstrated that the mean IL-6 response of cells from two European populations were higher than the mean response of cells from three African populations and that this difference was partially attributable to variation in rs1057807 allele frequency. In contrast, a SNP associated with a different pro-inflammatory cytokine (rs2869462 associated with the chemokine CXCL10) exhibited an opposite response, underscoring the complexity of how different genetic variants contribute to an individual’s immune response. This multidisciplinary study has identified a long-range chromatin interaction and genetic variation that regulates TLR6 to broaden our understanding of how human genetic variation affects the C. trachomatis-induced immune response., Chlamydia infections have diverse outcomes, and the causes of these variations are poorly understood. Here, we leverage a cellular genome-wide association approach to characterize a locus near the TLR10/TLR1/TLR6 region that is associated with levels of Chlamydia-induced cytokines, interacts with the TLR6 promoter, and contributes to differences among populations.

Chlamydia trachomatis is the leading cause of sexually transmitted bacterial infections and a major threat to women's reproductive health in particular. This obligate intracellular pathogen resides and replicates within a cellular compartment termed an inclusion, where it is sheltered by unknown mechanisms from gamma-interferon (IFNγ)-induced cell-autonomous host immunity. Through a genetic screen, we uncovered the Chlamydia inclusion membrane protein gamma resistance determinant (GarD) as a bacterial factor protecting inclusions from cell-autonomous immunity. In IFNγ-primed human cells, inclusions formed by garD loss-of-function mutants become decorated with linear ubiquitin and are eliminated. Leveraging cellular genome-wide association data, we identified the ubiquitin E3 ligase RNF213 as a candidate anti-Chlamydia protein. We demonstrate that IFNγ-inducible RNF213 facilitates the ubiquitylation and destruction of GarD-deficient inclusions. Furthermore, we show that GarD operates as a cis-acting stealth factor barring RNF213 from targeting inclusions, thus functionally defining GarD as an RNF213 antagonist essential for chlamydial growth during IFNγ-stimulated immunity.

The mechanisms underlying tumor immunosurveillance and their association with the immune-related adverse events (irAEs) associated with checkpoint inhibitor immunotherapies remain poorly understood. We describe a metastatic melanoma patient exhibiting multiple episodes of spontaneous disease regression followed by the development of several irAEs during the course of anti-programmed cell death protein 1 antibody immunotherapy. Whole-exome next-generation sequencing studies revealed this patient to harbor a pyrin inflammasome variant previously described to be associated with an atypical presentation of familial Mediterranean fever. This work highlights a potential role for inflammasomes in the regulation of tumor immunosurveillance and the pathogenesis of irAEs.

Despite being in a golden age of bacterial epigenomics, little work has systematically examined the plasticity and functional impacts of the bacterial DNA methylome. Here, we leveraged SMRT sequencing to examine the m6A DNA methylome of two Salmonella enterica ser. Typhimurium strains: 14028s and a ΔmetJ mutant with derepressed methionine metabolism, grown in Luria Broth or a media that simulates the intracellular environment. We find that the methylome is remarkably static—over 95% of adenosine bases retain their methylation status across conditions. Integration of methylation with transcriptomic data revealed limited correlation between changes in methylation and gene expression. Further, examining the transcriptome in ΔyhdJ bacteria, lacking the m6A methylase with the most dynamic methylation pattern in our dataset, revealed little evidence of YhdJ-mediated gene regulation. Curiously, despite G(m6A)TC motifs being particularly resistant to change across conditions, incorporating dam mutants into our analyses revealed two examples where changes in methylation and transcription may be linked across conditions. This includes the novel finding that the ΔmetJ motility defect may be partially driven by hypermethylation of the chemotaxis gene tsr. Together, these data redefine the S. Typhimurium epigenome as a highly stable system that has rare, but important, roles in transcriptional regulation. Incorporating these lessons into future studies will be critical as we progress through the epigenomic era.ImportanceWhile recent breakthroughs have enabled intense study of bacterial DNA modifications, limitations in current work have potentiated a surprisingly untested narrative that DNA methylation is a common mechanism of the bacterial response to environmental conditions. Essentially, whether epigenetic regulation of bacterial transcription is a common, generalizable phenomenon is a critical unanswered question that we address here. We find that most DNA methylation is static in Salmonella enterica serovar Typhimurium, even when the bacteria are grown under dramatically different conditions that cause broad changes in the transcriptome. Further, even when the methylation of individual bases change, these changes generally do not correlate with changes in gene expression. Finally, we demonstrate methods by which data can be stratified in order to identify coupled changes in methylation and gene expression.

Leishmaniasis is a neglected tropical disease with diverse infection outcomes ranging from self-healing lesions, to progressive non-healing lesion, to metastatic spread and destruction of mucous membranes. Although resolution of cutaneous leishmaniasis is a classic example of type-1 immunity leading to well controlled self-healing lesions, an excess of type-1 related inflammation can contribute to immunopathology and metastatic spread of disease.Leishmaniagenetic diversity can contribute to variation in polarization and robustness of the immune response through differences in both pathogen sensing by the host and immune evasion by the parasite. In this study, we observed a difference in parasite chemokine suppression between theLeishmania (L.)subgenus and theViannia (V.)subgenus, which is associated with severe immune mediated pathology such as mucocutaneous leishmaniasis. WhileLeishmania (L.)subgenus parasites utilize the virulence factor and metalloprotease glycoprotein-63(gp63)to suppress the type-1 associated host chemokine CXCL10,L. (V.) panamensisdid not suppress CXCL10. To understand the molecular basis for the inter-species variation in chemokine suppression, we usedin silicomodeling of the primary amino acid sequence and protein crystal structures to identify a putative CXCL10-binding site on GP63. We found the putative CXCL10 binding site to be located in a region ofgp63under significant positive selection and that it varies from theL. majorwild-type sequence in allgp63alleles identified in theL. (V.) panamensisreference genome. We determined that the predicted binding site and adjacent positively selected amino acids are required for CXCL10 suppression by mutating wild-typeL. (L.) major gp63to theL. (V.) panamensisallele and demonstrating impaired cleavage of CXCL10 but not a non-specific protease substrate. Notably,Vianniaclinical isolates confirmed thatL. (V.) panemensisprimarily encodes non-CXCL10-cleavinggp63alleles. In contrast,L. (V.) braziliensishas an intermediate level of activity, consistent with this species having more equal proportions of both alleles at the CXCL10 binding site, possibly due to balancing selection. Our results demonstrate how parasite genetic diversity can contribute to variation in the host immune response toLeishmaniaspp. infection that may play critical roles in the outcome of infection.

Intestinal immunity is coordinated by specialized mononuclear phagocyte populations, constituted by a diversity of cell subsets. Although the cell subsets constituting the mononuclear phagocyte network are thought to be similar in both small and large intestine, these organs have distinct anatomy, microbial composition, and immunological demands. Whether these distinctions demand organ-specific mononuclear phagocyte populations with dedicated organ-specific roles in immunity are unknown. Here we implement a new strategy to subset murine intestinal mononuclear phagocytes and identify two novel subsets which are colon-specific: a macrophage subset and a Th17-inducing dendritic cell (DC) subset. Colon-specific DCs and macrophages co-expressed CD24 and CD14, and surprisingly, both were dependent on the transcription factor IRF4. Novel IRF4-dependent CD14+CD24+ macrophages were markedly distinct from conventional macrophages and failed to express classical markers including CX3CR1, CD64 and CD88, and surprisingly expressed little IL-10, which was otherwise robustly expressed by all other intestinal macrophages. We further found that colon-specific CD14+CD24+ mononuclear phagocytes were essential for Th17 immunity in the colon, and provide definitive evidence that colon and small intestine have distinct antigen presenting cell requirements for Th17 immunity. Our findings reveal unappreciated organ-specific diversity of intestine-resident mononuclear phagocytes and organ-specific requirements for Th17 immunity.

Background While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility. Results Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb)) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs associated with both severe COVID-19 and other human traits demonstrated colocalization of the GWAS signal at the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN). This finding points to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. Conclusions Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches. The iCPAGdb web portal is accessible at http://cpag.oit.duke.edu and the software code at https://github.com/tbalmat/iCPAGdb.

SUMMARY Infectious diseases (ID) represent a significant proportion of morbidity and mortality across the world. Host genetic variation is likely to contribute to ID risk and downstream clinical outcomes, but there is a need for a genetics-anchored framework to decipher molecular mechanisms of disease risk, infer causal effect on potential complications, and identify instruments for drug target discovery. Here we perform transcriptome-wide association studies (TWAS) of 35 clinical ID traits in a cohort of 23,294 individuals, identifying 70 gene-level associations with 26 ID traits. Replication in two large-scale biobanks provides additional support for the identified associations. A phenome-scale scan of the 70 gene-level associations across hematologic, respiratory, cardiovascular, and neurologic traits proposes a molecular basis for known complications of the ID traits. Using Mendelian Randomization, we then provide causal support for the effect of the ID traits on adverse outcomes. The rich resource of genetic information linked to serologic tests and pathogen cultures from bronchoalveolar lavage, sputum, sinus/nasopharyngeal, tracheal, and blood samples (up to 7,699 positive pathogen cultures across 92 unique genera) that we leverage provides a platform to interrogate the genetic basis of compartment-specific infection and colonization. To accelerate insights into cellular mechanisms, we develop a TWAS repository of gene-level associations in a broad collection of human tissues with 79 pathogen-exposure induced cellular phenotypes as a discovery and replication platform. Cellular phenotypes of infection by 8 pathogens included pathogen invasion, intercellular spread, cytokine production, and pyroptosis. These rich datasets will facilitate mechanistic insights into the role of host genetic variation on ID risk and pathophysiology, with important implications for our molecular understanding of potentially severe phenotypic outcomes. HIGHLIGHTS Atlas of genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) results for 35 clinical infectious disease (ID) phenotypes, with genome-wide and transcriptome-wide significant results for 13 and 26 clinical ID traits, respectively Phenome-scale scan of ID-associated genes across 197 hematologic, respiratory, cardiovascular, and neurologic traits, facilitating identification of genes associated with known complications of the ID traits Mendelian Randomization analysis, leveraging naturally occurring DNA sequence variation to perform “randomized controlled trials” to test the causal effect of ID traits on potential outcomes and complications A genomic resource of TWAS associations for 79 pathogen-induced cellular traits from High-throughput Human in vitrO Susceptibility Testing (Hi-HOST) across 44 tissues as a discovery and replication platform to enable in silico cellular microbiology and functional genomic experiments

Susceptibility to infectious diseases is determined by a complex interaction between host and pathogen. For infections with the obligate intracellular bacterium Chlamydia trachomatis, variation in immune activation and disease presentation are regulated by both host genetic diversity and pathogen immune evasion. Previously, we discovered a single nucleotide polymorphism (rs2869462) associated with absolute abundance of CXCL10, a pro-inflammatory T-cell chemokine. Here, we report that levels of CXCL10 change during C. trachomatis infection of cultured cells in a manner dependent on both host and pathogen. Linear modeling of cellular traits associated with CXCL10 levels identified a strong, negative correlation with bacterial burden, suggesting that C. trachomatis actively suppresses CXCL10. We identified the pathogen-encoded factor responsible for this suppression as the chlamydial protease- or proteasome-like activity factor, CPAF. Further, we applied our modeling approach to other host cytokines in response to C. trachomatis and found evidence that RANTES, another T-cell chemoattractant, is actively suppressed by Chlamydia. However, this observed suppression of RANTES is not mediated by CPAF. Overall, our results demonstrate that CPAF suppresses CXCL10 to evade the host cytokine response and that modeling of cellular infection parameters can reveal previously unrecognized facets of host–pathogen interactions.

Leishmaniasis is a neglected tropical disease with diverse outcomes ranging from self-healing lesions, to progressive non-healing lesions, to metastatic spread and destruction of mucous membranes. Although resolution of cutaneous leishmaniasis is a classic example of type-1 immunity leading to self-healing lesions, an excess of type-1 related inflammation can contribute to immunopathology and metastatic spread. Leishmania genetic diversity can contribute to variation in polarization and robustness of the immune response through differences in both pathogen sensing by the host and immune evasion by the parasite. In this study, we observed a difference in parasite chemokine suppression between the Leishmania (L.) subgenus and the Viannia (V.) subgenus, which is associated with severe immune-mediated pathology such as mucocutaneous leishmaniasis. While Leishmania (L.) subgenus parasites utilize the virulence factor and metalloprotease glycoprotein-63 (gp63) to suppress the type-1 associated host chemokine CXCL10, L. (V.) panamensis did not suppress CXCL10. To understand the molecular basis for the inter-species variation in chemokine suppression, we used in silico modeling to identify a putative CXCL10-binding site on GP63. The putative CXCL10 binding site is in a region of gp63 under significant positive selection, and it varies from the L. major wild-type sequence in all gp63 alleles identified in the L. (V.) panamensis reference genome. Mutating wild-type L. (L.) major gp63 to the L. (V.) panamensis sequence at the putative binding site impaired cleavage of CXCL10 but not a non-specific protease substrate. Notably, Viannia clinical isolates confirmed that L. (V.) panamensis primarily encodes non-CXCL10-cleaving gp63 alleles. In contrast, L. (V.) braziliensis has an intermediate level of activity, consistent with this species having more equal proportions of both alleles. Our results demonstrate how parasite genetic diversity can contribute to variation in immune responses to Leishmania spp. infection that may play critical roles in the outcome of infection., Author summary Leishmaniasis is a neglected tropical disease caused by Leishmania parasites and spread by the bites of infected sand flies. Most cases of leishmaniasis present as self-healing sores that are resolved by a balanced immune response. Other cases of leishmaniasis involve spread to sites distant from the original bite, including damage of the inner surfaces of the mouth and nose. These cases of leishmaniasis involve an excessive immune response. Leishmania parasites produce virulence factor proteins, such as GP63, to trick the immune system into mounting a weaker response. GP63 specifically degrades signaling proteins that attract and activate certain immune cells. Here, we demonstrate that Leishmania parasite species have evolved to differ in their ability to degrade signaling proteins. In Leishmania species known to cause more immune-mediated tissue damage, the GP63 virulence factor has evolved to not degrade specific immune signaling proteins, thus attracting, and activating more immune cells. Our results demonstrate how diversity among Leishmania parasite species can contribute to variation in immune responses that may play critical roles in the outcome of infection.

Salmonella hijack host machinery in order to invade cells and establish infection. While considerable work has described the role of host proteins in invasion, much less is known regarding how natural variation in these invasion-associated host proteins affects Salmonella pathogenesis. Here we leveraged a candidate cellular GWAS screen to identify natural genetic variation in the ARHGEF26 (Rho Guanine Nucleotide Exchange Factor 26) gene that renders lymphoblastoid cells susceptible to Salmonella Typhi and Typhimurium invasion. Experimental follow-up redefined ARHGEF26’s role in Salmonella epithelial cell infection. Specifically, we identified complex serovar-by-host interactions whereby ARHGEF26 stimulation of S. Typhi and S. Typhimurium invasion into host cells varied in magnitude and effector-dependence based on host cell type. While ARHGEF26 regulated SopB- and SopE-mediated S. Typhi (but not S. Typhimurium) infection of HeLa cells, the largest effect of ARHGEF26 was observed with S. Typhimurium in polarized MDCK cells through a SopB- and SopE2-independent mechanism. In both cell types, knockdown of the ARHGEF26-associated protein DLG1 resulted in a similar phenotype and serovar specificity. Importantly, we show that ARHGEF26 plays a critical role in S. Typhimurium pathogenesis by contributing to bacterial burden in the enteric fever murine model, as well as inflammation in the colitis infection model. In the enteric fever model, SopB and SopE2 are required for the effects of Arhgef26 deletion on bacterial burden, and the impact of sopB and sopE2 deletion in turn required ARHGEF26. In contrast, SopB and SopE2 were not required for the impacts of Arhgef26 deletion on colitis. A role for ARHGEF26 on inflammation was also seen in cells, as knockdown reduced IL-8 production in HeLa cells. Together, these data reveal pleiotropic roles for ARHGEF26 during infection and highlight that many of the interactions that occur during infection that are thought to be well understood likely have underappreciated complexity., Author summary During infection, Salmonella manipulates host cells into engulfing the bacteria and establishing an intracellular niche. While many studies have identified genes involved in different stages of this Salmonella invasion process, few studies have examined how differences between human hosts contribute to infection susceptibility. Here we leveraged a candidate genetic screen to identify natural genetic variation in the human ARHGEF26 gene that correlates with Salmonella invasion. Springboarding from this result, we experimentally tested and redefined ARHGEF26’s role in Salmonella invasion, discovered a new role for ARHGEF26 in regulating inflammation during Salmonella disease, and demonstrated the relevance of these findings in mouse models. Building on how ARHGEF26 functions in other contexts, we implicated two ARHGEF26-interacting host proteins as contributors to Salmonella pathobiology. Collectively, these results identify a potential source of inter-person diversity in susceptibility to Salmonella disease and expand our molecular understanding of Salmonella infection to include a multifaceted role for ARHGEF26. They further identify important future directions in understanding how Salmonella recruit and manipulate ARHGEF26 as well as how ARHGEF26 is able to drive Salmonella-beneficial processes.

Clearance of intracellular pathogens, such asLeishmania (L.) major, depends on a well-regulated adaptive T cell response. Here we describe a pathogen-encoded mechanism to alter T cell recruitment by suppressing CXCL10, a chemokine that recruits CD4+ and CD8+ T cells by signaling through the CXCR3 receptor.L. majorsuppresses CXCL10 through the virulence factor and protease, glycoprotein-63(gp63).GP63 cleaves CXCL10 after amino acid A81, impairing T-cell chemotaxisin vitro.GP63 from either extracellular promastigotes or intracellular amastigotes is capable of cleaving CXCL10. Consistent with CXCL10 cleavage during infection, we observed GP63-mediated impairment of activated CD8+ T-cells in the draining lymph nodes of C57BL/6JWTmice. Correspondingly, in C57BL/6JWTmice,gp63deletion resulted in slower lesion development and a smaller maximum lesion size. However, infection in C57BL/6Jcxcr3−/−mice revealed the delay to lesion development required CXCR3 signaling. Finally,Salmonella entericaandChlamydia trachomatisinfection also suppress CXCL10, demonstrating convergent evolution of this immune evasion strategy in diverse intracellular pathogens. Understanding mechanisms of CXCL10 evasion may facilitate the development of novel therapeutic strategies to treat intracellular pathogens.Author SummaryLeishmaniasis is an infectious disease that affects over one million people annually. It is caused by intracellular parasites that have evolved to evade the host’s attempts to eliminate the parasite. The most common form of disease, cutaneous leishmaniasis, causes disfiguring skin lesions if the host immune system does not appropriately respond to the infection. A family of molecules called chemokines are critical for the recruitment of specific subsets of immune cells required to eliminate infection. Here, we demonstrate a novel mechanism that the parasiteLeishmania (L.) majoremploys to disrupt the immune response by interfering with chemokine signaling. By this strategy,L. majoruses a protease to cleave chemokines known to recruit T-cells required for fighting off infection. We show thatL. majorcleavage of these immune signals changes the recruitment of CD8+ T-cells and alters disease progression in mice. Finally, we observe that multiple common human intracellular pathogens, includingChlamydia trachomatisandSalmonella enterica, interfere with the same chemokines, suggesting a strong selective pressure to avoid this component of the immune response. Our study provides new insights into how intracellular pathogens interact with the host immune response to enhance pathogen survival and exacerbate disease outcomes.

Bacteria masterfully co-opt and subvert host signal transduction. As a paradigmatic example, Salmonella uses two type-3 secretion systems to inject effector proteins that facilitate Salmonella entry, establishment of an intracellular niche, and modulation of immune responses. We previously demonstrated that the Salmonella anti-inflammatory response activator SarA (Stm2585, GogC, PagJ, SteE) activates the host transcription factor STAT3 to drive expression of immunomodulatory STAT3-targets. Here, we demonstrate-by sequence, function, and biochemical measurement-that SarA mimics the cytoplasmic domain of glycoprotein 130 (gp130, IL6ST). SarA is phosphorylated at a YxxQ motif, facilitating binding to STAT3 with greater affinity than gp130. Departing from canonical gp130 signaling, SarA function is JAK-independent but requires GSK-3, a key regulator of metabolism and development. Our results reveal that SarA undergoes host phosphorylation to recruit a STAT3-activating complex, circumventing cytokine receptor activation. Effector mimicry of gp130 suggests GSK-3 can regulate normal cytokine signaling, potentially enabling metabolic and immune crosstalk.

The cholesterol storage disorder Niemann-Pick type C (NPC) disease is caused by defects in either of two late endosomal/lysosomal proteins, NPC1 and NPC2. NPC2 is a 16-kDa soluble protein that binds cholesterol in a 1:1 stoichiometry and can transfer cholesterol between membranes by a mechanism that involves protein-membrane interactions. To examine the structural basis of NPC2 function in cholesterol trafficking, a series of point mutations were generated across the surface of the protein. Several NPC2 mutants exhibited deficient sterol transport properties in a set of fluorescence-based assays. Notably, these mutants were also unable to promote egress of accumulated intracellular cholesterol from npc2(-/-) fibroblasts. The mutations mapped to several regions on the protein surface, suggesting that NPC2 can bind to more than one membrane simultaneously. Indeed, we have previously demonstrated that WT NPC2 promotes vesicle-vesicle interactions. These interactions were abrogated, however, by mutations causing defective sterol transfer properties. Molecular modeling shows that NPC2 is highly plastic, with several intense positively charged regions across the surface that could interact favorably with negatively charged membrane phospholipids. The point mutations generated in this study caused changes in NPC2 surface charge distribution with minimal conformational changes. The plasticity, coupled with membrane flexibility, probably allows for multiple cholesterol transfer routes. Thus, we hypothesize that, in part, NPC2 rapidly traffics cholesterol between closely appositioned membranes within the multilamellar interior of late endosomal/lysosomal proteins, ultimately effecting cholesterol egress from this compartment.

Pathogens have been a strong driving force for natural selection. Therefore, understanding how human genetic differences impact infection-related cellular traits can mechanistically link genetic variation to disease susceptibility. Here we report the Hi-HOST Phenome Project (H2P2): a catalog of cellular genome-wide association studies (GWAS) comprised of 79 infection-related phenotypes in response to 8 pathogens in 528 lymphoblastoid cell lines. Seventeen loci surpass genome-wide significance for infection-associated phenotypes ranging from pathogen replication to cytokine production. We combined H2P2 with clinical association data from patients to identify a SNP near CXCL10 as a risk factor for inflammatory bowel disease. A SNP in the transcriptional repressor ZBTB20 demonstrated pleiotropy, likely through suppression of multiple target genes, and was associated with viral hepatitis. These data are available on a web portal to facilitate interpreting human genome variation through the lens of cell biology and should serve as a rich resource for the research community.

In order to deploy virulence factors at appropriate times and locations, microbes must rapidly sense and respond to various metabolite signals. Previously, we showed a transient elevation of the methionine-derived metabolite methylthioadenosine (MTA) concentration in serum during systemic Salmonella enterica serovar Typhimurium infection., In order to deploy virulence factors at appropriate times and locations, microbes must rapidly sense and respond to various metabolite signals. Previously, we showed a transient elevation of the methionine-derived metabolite methylthioadenosine (MTA) concentration in serum during systemic Salmonella enterica serovar Typhimurium infection. Here we explored the functional consequences of increased MTA concentrations on S. Typhimurium virulence. We found that MTA, but not other related metabolites involved in polyamine synthesis and methionine salvage, reduced motility, host cell pyroptosis, and cellular invasion. Further, we developed a genetic model of increased bacterial endogenous MTA production by knocking out the master repressor of the methionine regulon, metJ. Like MTA-treated S. Typhimurium, the ΔmetJ mutant displayed reduced motility, host cell pyroptosis, and invasion. These phenotypic effects of MTA correlated with suppression of flagellar and Salmonella pathogenicity island 1 (SPI-1) networks. S. Typhimurium ΔmetJ had reduced virulence in oral and intraperitoneal infection of C57BL/6J mice independently of the effects of MTA on SPI-1. Finally, ΔmetJ bacteria induced a less severe inflammatory cytokine response in a mouse sepsis model. Together, these data indicate that exposure of S. Typhimurium to MTA or disruption of the bacterial methionine metabolism pathway suppresses S. Typhimurium virulence.

Eukaryotic parasites in the genusLeishmaniaplace approximately 350 million people per year at risk of disease. In addition to their global health significance,Leishmaniaspp. have served as an important model for delineating basic concepts in immunology such as T-helper cell polarization. There have been many qPCR-based assays reported for measuring parasite burden in humans and animals. However, these are largely optimized for use in clinical diagnosis and not specifically for animal models. This has led several of these assays to have suboptimal characteristics for use in animal models. For example, multi-copy number genes have been frequently used to increase sensitivity but are subject to greater plasticity within the genome and thus may confound effects of experimental manipulations in animal models. In this study, we developed a sybr-green based quantitative touchdown PCR assay for a highly conserved and single-copy putative RNA-binding protein, DRBD3. With primers that share greater than 90% sequence identity across all sequencedLeishmaniaspp., we demonstrate that this assay has a lower limit of detection of 100 fg of parasite DNA forLeishmania major,L. donovani,L. venezuelensis, andL. panamensis. Using C57BL6/J mice, we used this assay to monitor parasite burden over 1 month of infection with two strains ofL. major(Seidman and Friedlin), andL. venezeuelensis.These characteristics rival the sensitivity of previously reported qPCR based methods of parasite quantitation while amplifying a stable, single copy gene. Use of this protocol in the future will lead to improved accuracy in animal based models and help to tease apart differences in biology of host-parasite interactions.

Our studies of the SNP rs8060947 determine that the A allele is associated with increased invasion of Salmonella enterica serovar Typhi and increased susceptibility to typhoid fever (1). Gilchrist et al. (2) now provide evidence that the A allele is also associated with increased risk for bacteremia, and the association was driven primarily by nontyphoidal Salmonella , Streptococcus pneumoniae , Escherichia coli , and Acinetobacter , but not by Staphylococcus aureus, β-hemolytic streptococci, or Haemophilus influenza type b. Future studies are necessary to further validate the association between VAC14 and infectious diseases, determine the molecular mechanisms of protection, and investigate the … [↵][1]1To whom correspondence should be addressed. Email: dennis.ko{at}duke.edu. [1]: #xref-corresp-1-1

Calcineurin plays essential roles in virulence and growth of pathogenic fungi and is a target of the natural products FK506 and Cyclosporine A. In the pathogenic mucoralean fungus Mucor circinelloides, calcineurin mutation or inhibition confers a yeast-locked phenotype indicating that calcineurin governs the dimorphic transition. Genetic analysis in this study reveals that two calcineurin A catalytic subunits (out of three) are functionally diverged. Homology modeling illustrates modes of resistance resulting from amino substitutions in the interface between each calcineurin subunit and the inhibitory drugs. In addition, we show how the dimorphic transition orchestrated by calcineurin programs different outcomes during host-pathogen interactions. For example, when macrophages phagocytose Mucor yeast, subsequent phagosomal maturation occurs, indicating host cells respond appropriately to control the pathogen. On the other hand, upon phagocytosis of spores, macrophages fail to form mature phagosomes. Cytokine production from immune cells differs following exposure to yeast versus spores (which germinate into hyphae). Thus, the morphogenic transition can be targeted as an efficient treatment option against Mucor infection. In addition, genetic analysis (including gene disruption and mutational studies) further strengthens the understanding of calcineurin and provides a foundation to develop antifungal agents targeting calcineurin to deploy against Mucor and other pathogenic fungi.

Summary Salmonella enterica is an important foodborne pathogen that uses secreted effector proteins to manipulate host pathways to facilitate survival and dissemination. Different S. enterica serovars cause disease syndromes ranging from gastroenteritis to typhoid fever and vary in their effector repertoire. We leveraged this natural diversity to identify stm2585 , here designated sarA ( Salmonella anti-inflammatory response activator ), as a Salmonella effector that induces production of the anti-inflammatory cytokine IL-10. RNA-seq of cells infected with either Δ sarA or wild-type S. Typhimurium revealed that SarA activates STAT3 transcriptional targets. Consistent with this, SarA is necessary and sufficient for STAT3 phosphorylation, STAT3 inhibition blocks IL-10 production, and SarA and STAT3 interact by co-immunoprecipitation. These effects of SarA contribute to intracellular replication in vitro and bacterial load at systemic sites in mice. Our results demonstrate the power of using comparative genomics for identifying effectors and that Salmonella has evolved mechanisms for activating an important anti-inflammatory pathway.

Significance Salmonella enterica serovar Typhi ( S . Typhi) causes ∼20 million cases of typhoid fever every year. We carried out a genome-wide association study to identify genetic differences that correlate with the susceptibility of cells from hundreds of individuals to S . Typhi invasion. A SNP in VAC14 was associated with susceptibility to S . Typhi invasion and VAC14 expression. Cells mutated for VAC14 displayed increased S . Typhi docking due to increased plasma membrane cholesterol levels. The same SNP was associated with risk of typhoid fever in a Vietnamese population. Furthermore, treating zebrafish with a cholesterol-lowering drug reduced their susceptibility to S . Typhi infection. Therefore, this work demonstrates the power of coupling multiple genetic association studies with mechanistic dissection for understanding infectious disease susceptibility.

Summary Multiple sclerosis (MS) is an autoimmune disorder where T cells attack neurons in the central nervous system (CNS) leading to demyelination and neurological deficits. A driver of increased MS risk is the soluble form of the interleukin-7 receptor alpha chain gene (sIL7R) produced by alternative splicing of IL7R exon 6. Here, we identified the RNA helicase DDX39B as a potent activator of this exon and consequently a repressor of sIL7R, and we found strong genetic association of DDX39B with MS risk. Indeed, we showed that a genetic variant in the 5′ UTR of DDX39B reduces translation of DDX39B mRNAs and increases MS risk. Importantly, this DDX39B variant showed strong genetic and functional epistasis with allelic variants in IL7R exon 6. This study establishes the occurrence of biological epistasis in humans and provides mechanistic insight into the regulation of IL7R exon 6 splicing and its impact on MS risk.

The substrate of the methionine salvage pathway, methylthioadenosine, is a predictor of death in sepsis patients., Sepsis is a deleterious inflammatory response to infection with high mortality. Reliable sepsis biomarkers could improve diagnosis, prognosis, and treatment. Integration of human genetics, patient metabolite and cytokine measurements, and testing in a mouse model demonstrate that the methionine salvage pathway is a regulator of sepsis that can accurately predict prognosis in patients. Pathway-based genome-wide association analysis of nontyphoidal Salmonella bacteremia showed a strong enrichment for single-nucleotide polymorphisms near the components of the methionine salvage pathway. Measurement of the pathway’s substrate, methylthioadenosine (MTA), in two cohorts of sepsis patients demonstrated increased plasma MTA in nonsurvivors. Plasma MTA was correlated with levels of inflammatory cytokines, indicating that elevated MTA marks a subset of patients with excessive inflammation. A machine-learning model combining MTA and other variables yielded approximately 80% accuracy (area under the curve) in predicting death. Furthermore, mice infected with Salmonella had prolonged survival when MTA was administered before infection, suggesting that manipulating MTA levels could regulate the severity of the inflammatory response. Our results demonstrate how combining genetic data, biomolecule measurements, and animal models can shape our understanding of disease and lead to new biomarkers for patient stratification and potential therapeutic targeting.

Pyroptosis is proinflammatory cell death that occurs in response to certain microbes. Activation of the protease caspase-1 by molecular platforms called inflammasomes is required for pyroptosis. We performed a cellular genome-wide association study (GWAS) using Salmonella typhimurium infection of human lymphoblastoid cell lines as a means of dissecting the genetic architecture of susceptibility to pyroptosis and identifying unknown regulatory mechanisms. Cellular GWAS revealed that a common human genetic difference that regulates pyroptosis also alters microtubule stability. An intergenic single-nucleotide polymorphism on chromosome 18 is associated with decreased pyroptosis and increased expression of TUBB6 (tubulin, β 6 class V). TUBB6 is unique among tubulin isoforms in that its overexpression can completely disrupt the microtubule network. Cells from individuals with higher levels of TUBB6 expression have lower microtubule stability and less pyroptosis. Reducing TUBB6 expression or stabilizing microtubules pharmacologically with paclitaxel (Taxol) increases pyroptosis without affecting the other major readout of caspase-1 activation, interleukin-1β secretion. The results reveal a new role for microtubules and possibly specific tubulin isoforms in the execution of pyroptosis. Furthermore, the finding that there is common diversity in TUBB6 expression and microtubule stability could have broad consequences for other microtubule-dependent phenotypes, diseases, and pharmacological responses.

By performing two high-content small molecule screens on dextran sodium sulfate- and trinitrobenzene sulfonic acid-induced zebrafish enterocolitis models of inflammatory bowel disease, we have identified novel anti-inflammatory drugs from the John Hopkins Clinical Compound Library that suppress neutrophilic inflammation. Live imaging of neutrophil distribution was used to assess the level of acute inflammation and concurrently screen for off-target drug effects. Supporting the validity of our screening strategy, most of the anti-inflammatory drug hits were known antibiotics or anti-inflammatory agents. Novel hits included cholecystokinin (CCK) and dopamine receptor agonists. Using a pharmacological approach, we show that while CCK and dopamine receptor agonists alleviate enterocolitis-associated inflammation, receptor antagonists exacerbate inflammation in zebrafish. This work highlights the utility of small molecule screening in zebrafish enterocolitis models as a tool to identify novel bioactive molecules capable of modulating acute inflammation.

SifA is a virulence protein required for assembly and tubulation of a modified phagosome that promotes Salmonella replication. We show that SifA expressed in yeast induces membrane invagination during peroxisome proliferation and requires functional Rho1p. This is consistent with SifA ability to interact with RhoA and the fact that it is a GEF structural homologue., The Salmonella typhimurium effector protein SifA regulates the assembly and tubulation of the Salmonella phagosome. SifA localizes to the phagosome and interacts with the membrane via its prenylated tail. SifA is a structural homologue of another bacterial effector that acts as a GTP-exchange factor for Rho family GTPases and can bind GDP-RhoA. When coexpressed with a bacterial lipase that is activated by RhoA, SifA can induce tubulation of mammalian endosomes. In an effort to develop a genetic system to study SifA function, we expressed SifA and characterized its activity in yeast. GFP-SifA predominantly localized to yeast peroxisomal membranes. Under peroxisome-inducing conditions, GFP-SifA reduced the number of free peroxisomes and promoted the formation of large peroxisomes with membrane invaginations. GFP-SifA activity depended on the recruitment to peroxisomes of wild-type Rho1p and Pex25p, a receptor for Rho1p. GFP-SifA could also rescue the actin organization defects in pex25Δ and rho1 mutants, suggesting that SifA may recruit and potentiate Rho1p activity. We reexamined the distribution of GFP-SifA in mammalian cells and found the majority colocalizing with LAMP1-positive compartment and not with the peroxisomal marker PMP70. Together, these data suggest that SifA may use a similar mode of action via Rho proteins to alter yeast peroxisomal and mammalian endosomal membranes. Further definition of SifA activity on yeast peroxisomes could provide more insight into its role in regulating host membrane dynamics and small GTPases.

Motivation: Genome-wide association studies are beginning to elucidate how our genetic differences contribute to susceptibility and severity of disease. While computational tools have previously been developed to support various aspects of genome-wide association studies, there is currently a need for informatics solutions that facilitate the integration of data from multiple sources. Results: Here we present GWAS Analyzer, a database driven web-based tool that integrates genotype and phenotype data, association analysis results and genomic annotations from multiple public resources. GWAS Analyzer contains features for browsing these interrelated data, exploring phenotypic values by family or genotype, and filtering association results based on multiple criteria. The utility of the tool has been demonstrated by a genome-wide association study of human in vitro susceptibility to bacterial infection. GWAS Analyzer facilitated management of large sets of phenotype and genotype data, analysis of phenotypic variation and heritability, and most importantly, generation of a refined set of candidate single nucleotide polymorphisms (SNPs). The tool revealed a SNP that was experimentally validated to be associated with increased cell death among Salmonella infected HapMap cell lines. Availability: http://www.nwrce.org/gwas-analyzer Contact: mbrittna@u.washington.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

Gene duplication has been proposed to serve as the engine of evolutionary innovation. It is well recognized that eukaryotic genomes contain a large number of duplicated genes that evolve new functions or expression patterns. However, in mollusks, the evolutionary mechanisms underlying the divergence and the functional maintenance of duplicate genes remain little understood. In the present study, we performed a comprehensive analysis of duplicate genes in the protein kinase superfamily using whole genome and transcriptome data for the Pacific oyster. A total of 64 duplicated gene pairs were identified based on a phylogenetic approach and the reciprocal best BLAST method. By analyzing gene expression from RNA-seq data from 69 different developmental and stimuli-induced conditions (nine tissues, 38 developmental stages, eight dry treatments, seven heat treatments, and seven salty treatments), we found that expression patterns were significantly correlated for a number of duplicate gene pairs, suggesting the conservation of regulatory mechanisms following divergence. Our analysis also identified a subset of duplicate gene pairs with very high expression divergence, indicating that these gene pairs may have been subjected to transcriptional subfunctionalization or neofunctionalization after the initial duplication events. Further analysis revealed a significant correlation between expression and sequence divergence (as revealed by synonymous or nonsynonymous substitution rates) under certain conditions. Taken together, these results provide evidence for duplicate gene sequence and expression divergence in the Pacific oyster, accompanying its adaptation to harsh environments. Our results provide new insights into the evolution of duplicate genes and their expression levels in the Pacific oyster.

Meta-analyses of genome-wide association studies (GWAS) have demonstrated that the same genetic variants can be associated with multiple diseases and other complex traits. We present software called CPAG (Cross-Phenotype Analysis of GWAS) to look for similarities between 700 traits, build trees with informative clusters, and highlight underlying pathways. Clusters are consistent with pre-defined groups and literature-based validation but also reveal novel connections. We report similarity between plasma palmitoleic acid and Crohn's disease and find that specific fatty acids exacerbate enterocolitis in zebrafish. CPAG will become increasingly powerful as more genetic variants are uncovered, leading to a deeper understanding of complex traits. CPAG is freely available at www.sourceforge.net/projects/CPAG/. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0722-1) contains supplementary material, which is available to authorized users.