An Atlas of Genetic Variation Linking Pathogen-Induced Cellular Traits to Human Disease

Genetic influence on vascular smooth muscle cell apoptosis

Vascular smooth muscle cell (VSMC) proliferation, migration, and apoptosis play important roles in many physiological processes and pathological conditions. To identify genetic influences on VSMC behavior, we measured these traits and undertook genome-wide association studies in primary umbilical artery-derived VSMCs from >2000 individuals. Although there were no genome-wide significant associations for VSMC proliferation or migration, genetic variants at two genomic loci (7p15.3 and 7q32.3) showed highly significant associations with VSMC apoptosis (P = 1.95 × 10-13 and P = 7.47 × 10-9, respectively). The lead variant at the 7p51.3 locus was associated with increased expression of the GSDME and PALS2 genes in VSMCs. Knockdown of GSDME or PALS2 in VSMCs attenuated apoptotic cell death. A protein co-immunoprecipitation assay indicated that GSDME complexed with PALS2. PALS2 knockdown attenuated activated caspase-3 and GSDME fragmentation, whilst GSDME knockdown also reduced activated caspase-3. These findings provide new insights into the genetic regulation of VSMC apoptosis, with potential utility for therapeutic development.

Zinc Finger and BTB Domain-Containing 20: A Newly Emerging Player in Pathogenesis and Development of Human Cancers

Zinc finger and BTB domain-containing 20 (ZBTB20), which was initially identified in human dendritic cells, belongs to a family of transcription factors (TFs) with an N-terminal BTB domain and one or more C-terminal DNA-binding zinc finger domains. Under physiological conditions, ZBTB20 acts as a transcriptional repressor in cellular development and differentiation, metabolism, and innate immunity. Interestingly, multiple lines of evidence from mice and human systems have revealed the importance of ZBTB20 in the pathogenesis and development of cancers. ZBTB20 is not only a hotspot of genetic variation or fusion in many types of human cancers, but also a key TF or intermediator involving in the dysregulation of cancer cells. Given the diverse functions of ZBTB20 in both health and disease, we herein summarize the structure and physiological roles of ZBTB20, with an emphasis on the latest findings on tumorigenesis and cancer progression.

Baseline immune states (BIS) associated with vaccine responsiveness and factors that shape the BIS

Vaccines are among the greatest inventions in medicine, leading to the elimination or control of numerous diseases, including smallpox, polio, measles, rubella, and, most recently, COVID-19. Yet, the effectiveness of vaccines varies among individuals. In fact, while some recipients mount a robust response to vaccination that protects them from the disease, others fail to respond. Multiple clinical and epidemiological factors contribute to this heterogeneity in responsiveness. Systems immunology studies fueled by advances in single-cell biology have been instrumental in uncovering pre-vaccination immune cell types and genomic features (i.e., the baseline immune state, BIS) that have been associated with vaccine responsiveness. Here, we review clinical factors that shape the BIS, and the characteristics of the BIS associated with responsiveness to frequently studied vaccines (i.e., influenza, COVID-19, bacterial pneumonia, malaria). Finally, we discuss potential strategies to enhance vaccine responsiveness in high-risk groups, focusing specifically on older adults.

Pathogen-specific innate immune response patterns are distinctly affected by genetic diversity

Innate immune responses vary by pathogen and host genetics. We analyze quantitative trait loci (eQTLs) and transcriptomes of monocytes from 215 individuals stimulated by fungal, Gram-negative or Gram-positive bacterial pathogens. We identify conserved monocyte responses to bacterial pathogens and a distinct antifungal response. These include 745 response eQTLs (reQTLs) and corresponding genes with pathogen-specific effects, which we find first in samples of male donors and subsequently confirm for selected reQTLs in females. reQTLs affect predominantly upregulated genes that regulate immune response via e.g., NOD-like, C-type lectin, Toll-like and complement receptor-signaling pathways. Hence, reQTLs provide a functional explanation for individual differences in innate response patterns. Our identified reQTLs are also associated with cancer, autoimmunity, inflammatory and infectious diseases as shown by external genome-wide association studies. Thus, reQTLs help to explain interindividual variation in immune response to infection and provide candidate genes for variants associated with a range of diseases.

Human variation impacting MCOLN2 restricts Salmonella Typhi replication by magnesium deprivation

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 defense during bacterial infection is nutritional immunity: host cells attempt to restrict bacterial replication by denying bacteria access to key nutrients or supplying toxic metabolites. 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 availability-demonstrates that the divalent cation channel mucolipin-2 (MCOLN2 or TRPML2) restricts S. Typhi intracellular replication through magnesium deprivation. Mg2+ currents, conducted through MCOLN2 and out of endolysosomes, were measured directly using patch-clamping of the endolysosomal membrane. Our results reveal Mg2+ limitation as a key component of nutritional immunity against S. Typhi and as a source of variable host resistance.

Serine/Threonine Phosphatase Calcineurin Orchestrates the Intrinsic Resistance to Micafungin in the Human-Pathogenic Fungus Mucor circinelloides

Procedures such as solid-organ transplants and cancer treatments can leave many patients in an immunocompromised state. This leads to their increased susceptibility to opportunistic diseases such as fungal infections. Mucormycosis infections are continually emerging and pose a serious threat to immunocompromised patients. Recently there has been a sharp increase in mucormycosis cases as a secondary infection in patients battling severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections. Mucorales fungi are notorious for presenting resistance to most antifungal drugs. The absence of effective means to treat these infections results in mortality rates approaching 100% in cases of disseminated infection. One of the most effective antifungal drug classes currently available is the echinocandins. Echinocandins seem to be efficacious in the treatment of many other fungal infections. Unfortunately, susceptibility testing has found that echinocandins have little to no effect on Mucorales fungi. In this study, we found that the model Mucorales Mucor circinelloides genome carries three copies of the genes encoding the echinocandin target protein β-(1,3)-d-glucan synthase (fksA, fksB, and fksC). Interestingly, we found that exposing M. circinelloides to micafungin significantly increased the expression of the fksA and fksB genes, resulting in an increased accumulation of β-(1,3)-d-glucan on the cell walls. However, this overexpression of the fks genes is not directly connected to the intrinsic resistance. Subsequent investigation discovered that the serine/threonine phosphatase calcineurin regulates the expression of fksA and fksB, and the deletion of calcineurin results in a decrease in expression of all three fks genes. Deletion of calcineurin also results in a lower minimum effective concentration (MEC) of micafungin. In addition, we found that duplication of the fks gene is also responsible for the intrinsic resistance, in which lack of either fksA or fksB led a lower MEC of micafungin. Together, these findings demonstrate that calcineurin and fks gene duplication contribute to the intrinsic resistance to micafungin we observe in M. circinelloides.

GarD-ing the pathogen-containing vacuole from destruction

The human pathogen Chlamydia trachomatis evades killing by IFNγ-activated mechanisms, yet how this occurs remains unclear. In this issue of Cell Host & Microbe, Walsh et al. identify an IFNγ-dependent antimicrobial mechanism mediated by the host ubiquitin ligase RNF213 that is evaded by the Chlamydia effector GarD.

The bacterial effector GarD shields Chlamydia trachomatis inclusions from RNF213-mediated ubiquitylation and destruction

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 Awesome Power of Human Genetics of Infectious Disease

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.

Single-cell genome-wide association reveals that a nonsynonymous variant in ERAP1 confers increased susceptibility to influenza virus

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.

A human iPSC-array-based GWAS identifies a virus susceptibility locus in the NDUFA4 gene and functional variants

Population-based studies to identify disease-associated risk alleles typically require samples from a large number of individuals. Here, we report a human-induced pluripotent stem cell (hiPSC)-based screening strategy to link human genetics with viral infectivity. A genome-wide association study (GWAS) identified a cluster of single-nucleotide polymorphisms (SNPs) in a cis-regulatory region of the NDUFA4 gene, which was associated with susceptibility to Zika virus (ZIKV) infection. Loss of NDUFA4 led to decreased sensitivity to ZIKV, dengue virus, and SARS-CoV-2 infection. Isogenic hiPSC lines carrying non-risk alleles of SNPs or deletion of the cis-regulatory region lower sensitivity to viral infection. Mechanistic studies indicated that loss/reduction of NDUFA4 causes mitochondrial stress, which leads to the leakage of mtDNA and thereby upregulation of type I interferon signaling. This study provides proof-of-principle for the application of iPSC arrays in GWAS and identifies NDUFA4 as a previously unknown susceptibility locus for viral infection.

Integration of the Salmonella Typhimurium Methylome and Transcriptome Reveals That DNA Methylation and Transcriptional Regulation Are Largely Decoupled under Virulence-Related Conditions

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⁶A DNA methylome of two Salmonella enterica serovar Typhimurium strains: 14028s and a ΔmetJ mutant with derepressed methionine metabolism, grown in Luria broth or medium that simulates the intracellular environment. We found that the methylome is remarkably static: >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, examination of the transcriptome in ΔyhdJ bacteria lacking the m⁶A methylase with the most dynamic methylation pattern in our data set revealed little evidence of YhdJ-mediated gene regulation. Curiously, despite G(m⁶A)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.

Understanding the genetic basis of immune responses to fungal infection

INTRODUCTION

Fungal infections represent a global public health problem that affect millions of people. Despite remarkable advances achieved over the last decades, available diagnostic and therapeutic tools remain insufficient for the optimal management of these diseases. The clinical course of fungal infection is highly variable, and evidence accumulated from patients with rare mutations and cohort-based studies suggests that the trajectory of disease is largely defined by patient genetics and its impact on immune responses. Therefore, there is an urgent need to elucidate the precise mechanisms by which which genetic variants influence the risk, progression, and outcome of fungal infection.

AREAS COVERED

In this review, we highlight recent advances in our understanding of the genetic factors that influence antifungal immune responses based on candidate gene studies and genome-wide approaches performed in different experimental and clinical models.

EXPERT OPINION

Research on genetics of susceptibility to infection is expected to lead to a detailed knowledge framework for the pathogenesis of human fungal infections and unveil novel targets and pathways amenable to clinical intervention.

A functional genomics approach in Tanzanian population identifies distinct genetic regulators of cytokine production compared to European population

Humans exhibit remarkable interindividual and interpopulation immune response variability upon microbial challenges. Cytokines play a vital role in regulating inflammation and immune responses, but dysregulation of cytokine responses has been implicated in different disease states. Host genetic factors were previously shown to significantly impact cytokine response heterogeneity mainly in European-based studies, but it is unclear whether these findings are transferable to non-European individuals. Here, we aimed to identify genetic variants modulating cytokine responses in healthy adults of East African ancestry from Tanzania. We leveraged both cytokine and genetic data and performed genome-wide cytokine quantitative trait loci (cQTLs) mapping. The results were compared with another cohort of healthy adults of Western European ancestry via direct overlap and functional enrichment analyses. We also performed meta-analyses to identify cQTLs with congruent effect direction in both populations. In the Tanzanians, cQTL mapping identified 80 independent suggestive loci and one genome-wide significant locus (TBC1D22A) at chromosome 22; SNP rs12169244 was associated with IL-1b release after Salmonella enteritidis stimulation. Remarkably, the identified cQTLs varied significantly when compared to the European cohort, and there was a very limited percentage of overlap (1.6% to 1.9%). We further observed ancestry-specific pathways regulating induced cytokine responses, and there was significant enrichment of the interferon pathway specifically in the Tanzanians. Furthermore, contrary to the Europeans, genetic variants in the TLR10-TLR1-TLR6 locus showed no effect on cytokine response. Our data reveal both ancestry-specific effects of genetic variants and pathways on cytokine response heterogeneity, hence arguing for the importance of initiatives to include diverse populations into genomics research.

Balancing Selection of the Intracellular Pathogen Response in Natural Caenorhabditis elegans Populations

Genetic variation in host populations may lead to differential viral susceptibilities. Here, we investigate the role of natural genetic variation in the Intracellular Pathogen Response (IPR), an important antiviral pathway in the model organism Caenorhabditis elegans against Orsay virus (OrV). The IPR involves transcriptional activity of 80 genes including the pals-genes. We examine the genetic variation in the pals-family for traces of selection and explore the molecular and phenotypic effects of having distinct pals-gene alleles. Genetic analysis of 330 global C. elegans strains reveals that genetic diversity within the IPR-related pals-genes can be categorized in a few haplotypes worldwide. Importantly, two key IPR regulators, pals-22 and pals-25, are in a genomic region carrying signatures of balancing selection, suggesting that different evolutionary strategies exist in IPR regulation. We infected eleven C. elegans strains that represent three distinct pals-22 pals-25 haplotypes with Orsay virus to determine their susceptibility. For two of these strains, N2 and CB4856, the transcriptional response to infection was also measured. The results indicate that pals-22 pals-25 haplotype shapes the defense against OrV and host genetic variation can result in constitutive activation of IPR genes. Our work presents evidence for balancing genetic selection of immunity genes in C. elegans and provides a novel perspective on the functional diversity that can develop within a main antiviral response in natural host populations.

Human genetic diversity regulating the TLR10/TLR1/TLR6 locus confers increased cytokines in response to Chlamydia trachomatis

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.

Regulation of intestinal immunity and tissue repair by enteric glia

Tissue maintenance and repair depend on the integrated activity of multiple cell types¹. Whereas the contributions of epithelial^2,3, immune^4,5 and stromal cells^6,7 in intestinal tissue integrity are well understood, the role of intrinsic neuroglia networks remains largely unknown. Here we uncover important roles of enteric glial cells (EGCs) in intestinal homeostasis, immunity and tissue repair. We demonstrate that infection of mice with Heligmosomoides polygyrus leads to enteric gliosis and the upregulation of an interferon gamma (IFNγ) gene signature. IFNγ-dependent gene modules were also induced in EGCs from patients with inflammatory bowel disease⁸. Single-cell transcriptomics analysis of the tunica muscularis showed that glia-specific abrogation of IFNγ signalling leads to tissue-wide activation of pro-inflammatory transcriptional programs. Furthermore, disruption of the IFNγ-EGC signalling axis enhanced the inflammatory and granulomatous response of the tunica muscularis to helminths. Mechanistically, we show that the upregulation of Cxcl10 is an early immediate response of EGCs to IFNγ signalling and provide evidence that this chemokine and the downstream amplification of IFNγ signalling in the tunica muscularis are required for a measured inflammatory response to helminths and resolution of the granulomatous pathology. Our study demonstrates that IFNγ signalling in enteric glia is central to intestinal homeostasis and reveals critical roles of the IFNγ-EGC-CXCL10 axis in immune response and tissue repair after infectious challenge.

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

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

ARHGEF26 enhances Salmonella invasion and inflammation in cells and mice

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.

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.

Response to Oxidative Burst-Induced Hypoxia Is Associated With Macrophage Inflammatory Profiles as Revealed by Cellular Genome-Wide Association

Background

In mammalian species, hypoxia is a prominent feature of inflammation. The role of hypoxia in regulating macrophage responses via alteration in metabolic pathways is well established. Recently, oxidative burst-induced hypoxia has been shown in murine macrophages after phagocytosis. Despite the available detailed information on the regulation of macrophage function at transcriptomic and epigenomic levels, the association of genetic polymorphism and macrophage function has been less explored. Previously, we have shown that host genetics controls approximately 80% of the variation in an oxidative burst as measured by nitric oxide (NO^-). Further studies revealed two clusters of transcription factors (hypoxia-related and inflammatory-related) are under the genetic control that shapes macrophages’ pro-inflammatory characteristics.

Material and Methods

In the current study, the association between 43,066 autosomal Single Nucleic Polymorphism (SNPs) and the ability of MDMs in production of NO^- in response to E. coli was evaluated in 58 Holstein cows. The positional candidate genes near significant SNPs were selected to perform functional analysis. In addition, the interaction between the positional candidate genes and differentially expressed genes from our previous study was investigated.

Results

Sixty SNPs on 22 chromosomes of the bovine genome were found to be significantly associated with NO^- production of macrophages. The functional genomic analysis showed a significant interaction between positional candidate genes and mitochondria-related differentially expressed genes from the previous study. Further examination showed 7 SNPs located in the vicinity of genes with roles in response to hypoxia, shaping approximately 73% of the observed individual variation in NO^- production by MDM. Regarding the normoxic condition of macrophage culture in this study, it was hypothesized that oxidative burst is responsible for causing hypoxia at the cellular level.

Conclusion

The results suggest that the genetic polymorphism via regulation of response to hypoxia is a candidate step that perhaps shapes macrophage functional characteristics in the pathway of phagocytosis leading to oxidative burst, hypoxia, cellular response to hypoxia and finally the pro-inflammatory responses. Since all cells in one individual carry the same alleles, the effect of genetic predisposition of sensitivity to hypoxia will likely be notable on the clinical outcome to a broad range of host-pathogen interactions.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

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 .

These Are the Genes You're Looking For: Finding Host Resistance Genes

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.

COVID-19: The Effect of Host Genetic Variations on Host-Virus Interactions

Spurred into action by the COVID-19 pandemic, the global scientific community has, in a short of period of time, made astonishing progress in understanding and combating COVID-19. Given the known human protein machinery for (a) SARS-CoV-2 entry, (b) the host innate immune response, and (c) virus-host interactions (protein-protein and RNA-protein), the potential effects of human genetic variation in this machinery, which may contribute to clinical differences in SARS-CoV-2 pathogenesis and help determine individual risk for COVID-19 infection, are explored. The Genome Aggregation Database (gnomAD) was used to show that several rare germline exome variants of proteins in these pathways occur in the human population, suggesting that carriers of these rare variants (especially for proteins of innate immunity pathways) are at risk for severe symptoms (like the severe symptoms in patients who are known to be rare variant carriers), whereas carriers of other variants could have a protective advantage against infection. The occurrence of genetic variation is thus expected to motivate the experimental probing of natural variants to understand the mechanistic differences in SARS-CoV-2 pathogenesis from one individual to another.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

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. 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); http://cpag.oit.duke.edu) 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 with severe COVID-19 demonstrated colocalization of the GWAS signal of the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN), pointing to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. 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.

Host genetics and infectious disease: new tools, insights and translational opportunities

Understanding how human genetics influence infectious disease susceptibility offers the opportunity for new insights into pathogenesis, potential drug targets, risk stratification, response to therapy and vaccination. As new infectious diseases continue to emerge, together with growing levels of antimicrobial resistance and an increasing awareness of substantial differences between populations in genetic associations, the need for such work is expanding. In this Review, we illustrate how our understanding of the host–pathogen relationship is advancing through holistic approaches, describing current strategies to investigate the role of host genetic variation in established and emerging infections, including COVID-19, the need for wider application to diverse global populations mirroring the burden of disease, the impact of pathogen and vector genetic diversity and a broad array of immune and inflammation phenotypes that can be mapped as traits in health and disease. Insights from study of inborn errors of immunity and multi-omics profiling together with developments in analytical methods are further advancing our knowledge of this important area.

Infectious diseases are an ever-present global threat. In this Review, Kwok, Mentzer and Knight discuss our latest understanding of how human genetics influence susceptibility to disease. Furthermore, they discuss emerging progress in the interplay between host and pathogen genetics, molecular responses to infection and vaccination, and opportunities to bring these aspects together for rapid responses to emerging diseases such as COVID-19.

Modeling of variables in cellular infection reveals CXCL10 levels are regulated by human genetic variation and the Chlamydia-encoded CPAF protease

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.

Understanding Mucormycoses in the Age of "omics"

Mucormycoses are deadly invasive infections caused by several fungal species belonging to the subphylum Mucoromycotina, order Mucorales. Hallmarks of disease progression include angioinvasion and tissue necrosis that aid in fungal dissemination through the blood stream, causing deeper infections and resulting in poor penetration of antifungal agents to the site of infection. In the absence of surgical removal of the infected focus, antifungal therapy alone is rarely curative. Even when surgical debridement is combined with high-dose antifungal therapy, the mortality associated with mucormycoses is >50%. The unacceptably high mortality rate, limited options for therapy and the extreme morbidity of highly disfiguring surgical therapy provide a clear mandate to understand the molecular mechanisms that govern pathogenesis with the hopes of developing alternative strategies to treat and prevent mucormycoses. In the absence of robust forward and reverse genetic systems available for this taxonomic group of fungi, unbiased next generation sequence (NGS)-based approaches have provided much needed insights into our understanding of many aspects of Mucormycoses, including genome structure, drug resistance, diagnostic development, and fungus-host interactions. Here, we will discuss the specific contributions that NGS-based approaches have made to the field and discuss open questions that can be addressed using similar approaches.

Eliciting priors and relaxing the single causal variant assumption in colocalisation analyses

Horizontal integration of summary statistics from different GWAS traits can be used to evaluate evidence for their shared genetic causality. One popular method to do this is a Bayesian method, coloc, which is attractive in requiring only GWAS summary statistics and no linkage disequilibrium estimates and is now being used routinely to perform thousands of comparisons between traits. Here we show that while most users do not adjust default software values, misspecification of prior parameters can substantially alter posterior inference. We suggest data driven methods to derive sensible prior values, and demonstrate how sensitivity analysis can be used to assess robustness of posterior inference. The flexibility of coloc comes at the expense of an unrealistic assumption of a single causal variant per trait. This assumption can be relaxed by stepwise conditioning, but this requires external software and an LD matrix aligned to study alleles. We have now implemented conditioning within coloc, and propose a new alternative method, masking, that does not require LD and approximates conditioning when causal variants are independent. Importantly, masking can be used in combination with conditioning where allelically aligned LD estimates are available for only a single trait. We have implemented these developments in a new version of coloc which we hope will enable more informed choice of priors and overcome the restriction of the single causal variant assumptions in coloc analysis.

A Novel Resistance Pathway for Calcineurin Inhibitors in the Human-Pathogenic Mucorales Mucor circinelloides

Mucormycosis is an emerging lethal fungal infection in immunocompromised patients. Mucor circinelloides is a causal agent of mucormycosis and serves as a model system to understand genetics in Mucorales. Calcineurin is a conserved virulence factor in many pathogenic fungi, and calcineurin inhibition or deletion of the calcineurin regulatory subunit (CnbR) in Mucor results in a shift from hyphal to yeast growth. We analyzed 36 calcineurin inhibitor-resistant or bypass mutants that exhibited hyphal growth in the presence of calcineurin inhibitors or in the yeast-locked cnbRΔ mutant background without carrying any mutations in known calcineurin components. We found that a majority of the mutants had altered sequence in a gene, named here bycA (bypass of calcineurin). bycA encodes an amino acid permease. We verified that both the bycAΔ single mutant and the bycAΔ cnbRΔ double mutant are resistant to calcineurin inhibitor FK506, thereby demonstrating a novel mechanism of resistance against calcineurin inhibitors. We also found that the level of expression of bycA was significantly higher in the wild-type strain treated with FK506 and in the cnbRΔ mutants but was significantly lower in the wild-type strain without FK506 treatment. These findings suggest that bycA is a negative regulator of hyphal growth and/or a positive regulator of yeast growth in Mucor and that calcineurin suppresses expression of the bycA gene at the mRNA level to promote hyphal growth. BycA is involved in the Mucor hypha-yeast transition as our data demonstrate positive correlations among bycA expression, protein kinase A activity, and Mucor yeast growth. Also, calcineurin, independently of its role in morphogenesis, contributes to virulence traits, including phagosome maturation blockade, host cell damages, and proangiogenic growth factor induction during interactions with hosts.IMPORTANCEMucor is intrinsically resistant to most known antifungals, which makes mucormycosis treatment challenging. Calcineurin is a serine/threonine phosphatase that is widely conserved across eukaryotes. When calcineurin function is inhibited in Mucor, growth shifts to a less virulent yeast growth form, which makes calcineurin an attractive target for development of new antifungal drugs. Previously, we identified two distinct mechanisms through which Mucor can become resistant to calcineurin inhibitors involving Mendelian mutations in the gene for FKBP12, including mechanisms corresponding to calcineurin A or B subunits and epimutations silencing the FKBP12 gene. Here, we identified a third novel mechanism where loss-of-function mutations in the amino acid permease corresponding to the bycA gene contribute to resistance against calcineurin inhibitors. When calcineurin activity is absent, BycA can activate protein kinase A (PKA) to promote yeast growth via a cAMP-independent pathway. Our data also show that calcineurin activity contributes to host-pathogen interactions primarily in the pathogenesis of Mucor.

Transcriptomic Profiles of Monocyte-Derived Macrophages in Response to Escherichia coli is Associated with the Host Genetics

Reactive Nitrogen Species (RNS) are a group of bactericidal molecules produced by macrophages in response to pathogens in a process called oxidative burst. Nitric oxide (NO^-) is a member of RNS produced from arginine by inducible Nitric Oxide Synthase (iNOS) enzyme. The activity of iNOS and production of NO^- by macrophages following stimulation is one of the indicators of macrophage polarization towards M1/proinflammatory. Production of NO^- by bovine monocyte-derived macrophage (MDM) and mouse peritoneal macrophages has been shown to be strongly associated with host genetic with the heritability of 0.776 in bovine MDM and 0.8 in mouse peritoneal macrophages. However, the mechanism of genetic regulation of macrophage response has remained less explored. In the current study, the transcriptome of bovine MDMs was compared between two extreme phenotypes that had been classified as high and low responder based on NO^- production. The results showed that 179 and 392 genes were differentially expressed (DE) between high and low responder groups at 3 and 18 hours after exposure to Escherichia coli, respectively. A set of 11 Transcription Factors (TFs) (STAT1, IRF7, SPI1, STAT4, IRF1, HIF1A, FOXO3, REL, NFAT5, HIC1, and IRF4) at 3 hours and a set of 13 TFs (STAT1, IRF1, HIF1A, STAT4, ATF4, TP63, EGR1, CDKN2A, RBL1, E2F1, PRDM1, GATA3, and IRF4) at 18 hours after exposure to E. coli were identified to be differentially regulated between the high and low responder phenotypes. These TFs were found to be divided into two clusters of inflammatory- and hypoxia-related TFs. Functional analysis revealed that some key canonical pathways such as phagocytosis, chemotaxis, antigen presentation, and cell-to-cell signalling are enriched among the over-expressed genes by high responder phenotype. Based on the results of this study, it was inferred that the functional characteristics of bovine MDMs are associated with NO-based classification. Since NO^- production is strongly associated with host genetics, this study for the first time shows the distinct proinflammatory profiles of macrophages are controlled by the natural genetic polymorphism in an outbred population. In addition, the results suggest that genetics can be considered as a new dimension in the current model of macrophage polarization which is currently described by the combination of stimulants, only.

Pathogen Evasion of Chemokine Response Through Suppression of CXCL10

Clearance of intracellular pathogens, such as Leishmania (L.) major, depends on an immune response with well-regulated cytokine signaling. Here we describe a pathogen-mediated mechanism of evading CXCL10, a chemokine with diverse antimicrobial functions, including T cell recruitment. Infection with L. major in a human monocyte cell line induced robust CXCL10 transcription without increasing extracellular CXCL10 protein concentrations. We found that this transcriptionally independent suppression of CXCL10 is mediated by the virulence factor and protease, glycoprotein-63 (gp63). Specifically, GP63 cleaves CXCL10 after amino acid A81 at the base of a C-terminal alpha-helix. Cytokine cleavage by GP63 demonstrated specificity, as GP63 cleaved CXCL10 and its homologs, which all bind the CXCR3 receptor, but not distantly related chemokines, such as CXCL8 and CCL22. Further characterization demonstrated that CXCL10 cleavage activity by GP63 was produced by both extracellular promastigotes and intracellular amastigotes. Crucially, CXCL10 cleavage impaired T cell chemotaxis in vitro, indicating that cleaved CXCL10 cannot signal through CXCR3. Ultimately, we propose CXCL10 suppression is a convergent mechanism of immune evasion, as Salmonella enterica and Chlamydia trachomatis also suppress CXCL10. This commonality suggests that counteracting CXCL10 suppression may provide a generalizable therapeutic strategy against intracellular pathogens.

Importance

Leishmaniasis, an infectious disease that annually affects over one million people, is caused by intracellular parasites that have evolved to evade the host's attempts to eliminate the parasite. Cutaneous leishmaniasis results in disfiguring skin lesions if the host immune system does not appropriately respond to infection. A family of molecules called chemokines coordinate recruitment of the immune cells required to eliminate infection. Here, we demonstrate a novel mechanism that Leishmania (L.) spp. employ to suppress host chemokines: a Leishmania-encoded protease cleaves chemokines known to recruit T cells that fight off infection. We observe that other common human intracellular pathogens, including Chlamydia trachomatis and Salmonella enterica, reduce levels of 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.

The effect of host genetics on in vitro performance of bovine monocyte-derived macrophages

The dynamic interaction between the host and pathogens, along with environmental factors, influences the regulation of mammalian immune responses. Therefore, comprehensive in vivo immune-phenotyping during an active response to a pathogen can be complex and prone to confounding effects. Evaluating critical fundamental aspects of the immune system at a cellular level is an alternative approach to reduce this complexity. Therefore, the objective of the current study was to examine an in vitro model for functional phenotyping of bovine monocyte-derived macrophages (MDM), cells which play a crucial role at all phases of inflammation, as well influence downstream immune responses. As indicators of MDM function, phagocytosis and nitric oxide (NO^-) production were tested in MDM of 16 cows in response to 2 common bacterial pathogens of dairy cows, Escherichia coli and Staphylococcus aureus. Notable functional variations were observed among the individuals (coefficient of variation: 33% for phagocytosis and 70% in the production of NO^-). The rank correlation analysis revealed a significant, positive, and strong correlation (rho = 0.92) between NO^- production in response to E. coli and S. aureus, and a positive but moderate correlation (rho = 0.58) between phagocytosis of E. coli and S. aureus. To gain further insight into this trait, another 58 cows were evaluated solely for NO^- response against E. coli. The pedigree of the tested animals was added to the statistical model and the heritability was estimated to be 0.776. Overall, the finding of this study showed a strong effect of host genetics on the in vitro activities of MDM and the possibility of ranking Holstein cows based on the in vitro functional variation of MDM.

Natural variation in the roles of C. elegans autophagy components during microsporidia infection

Natural genetic variation can determine the outcome of an infection, and often reflects the co-evolutionary battle between hosts and pathogens. We previously found that a natural variant of the nematode Caenorhabditis elegans from Hawaii (HW) has increased resistance against natural microsporidian pathogens in the Nematocida genus, when compared to the standard laboratory strain of N2. In particular, HW animals can clear infection, while N2 animals cannot. In addition, HW animals have lower levels of initial colonization of Nematocida inside intestinal cells, compared to N2. Here we investigate how this natural variation in resistance relates to autophagy. We found that there is much better targeting of autophagy-related machinery to parasites under conditions where they are cleared. In particular, ubiquitin targeting to Nematocida cells correlates very well with their subsequent clearance in terms of timing, host strain and age, as well as species of Nematocida. Furthermore, clearance correlates with targeting of the LGG-2/LC3 autophagy protein to parasite cells, with HW animals having much more efficient targeting of LGG-2 to parasite cells than N2 animals. Surprisingly, however, we found that LGG-2 is not required to clear infection. Instead, we found that LGG-2/LC3 regulates Nematocida colonization inside intestinal cells. Interestingly, LGG-2/LC3 regulates intracellular colonization only in the HW strain, and not in N2. Altogether these results demonstrate that there is natural genetic variation in an LGG-2-dependent process that regulates microsporidia colonization inside intestinal cells, although not microsporidia clearance.

Candida albicans Morphology-Dependent Host FGF-2 Response as a Potential Therapeutic Target

Angiogenesis mediated by proteins such as Fibroblast Growth Factor-2 (FGF-2) is a vital component of normal physiological processes and has also been implicated in contributing to the disease state associated with various microbial infections. Previous studies by our group and others have shown that Candida albicans, a common agent of candidiasis, induces FGF-2 secretion in vitro and angiogenesis in brains and kidneys during systemic infections. However, the underlying mechanism(s) via which the fungus increases FGF-2 production and the role(s) that FGF-2/angiogenesis plays in C. albicans disease remain unknown. Here we show, for the first time, that C. albicans hyphae (and not yeast cells) increase the FGF-2 response in human endothelial cells. Moreover, Candidalysin, a toxin secreted exclusively by C. albicans in the hyphal state, is required to induce this response. Our in vivo studies show that in the systemic C. albicans infection model, mice treated with FGF-2 exhibit significantly higher mortality rates when compared to untreated mice not given the angiogenic growth factor. Even treatment with fluconazole could not fully rescue infected animals that were administered FGF-2. Our data suggest that the increase of FGF-2 production/angiogenesis induced by Candidalysin contributes to the pathogenicity of C. albicans.