Human Epistatic Interaction Controls IL7R Splicing and Increases Multiple Sclerosis Risk.

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 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.

Host genotype-specific therapies can optimize the inflammatory response to mycobacterial infections.

Susceptibility to tuberculosis is historically ascribed to an inadequate immune response that fails to control infecting mycobacteria. In zebrafish, we find that susceptibility to Mycobacterium marinum can result from either inadequate or excessive acute inflammation. Modulation of the leukotriene A(4) hydrolase (LTA4H) locus, which controls the balance of pro- and anti-inflammatory eicosanoids, reveals two distinct molecular routes to mycobacterial susceptibility converging on dysregulated TNF levels: inadequate inflammation caused by excess lipoxins and hyperinflammation driven by excess leukotriene B(4). We identify therapies that specifically target each of these extremes. In humans, we identify a single nucleotide polymorphism in the LTA4H promoter that regulates its transcriptional activity. In tuberculous meningitis, the polymorphism is associated with inflammatory cell recruitment, patient survival and response to adjunctive anti-inflammatory therapy. Together, our findings suggest that host-directed therapies tailored to patient LTA4H genotypes may counter detrimental effects of either extreme of inflammation.

Targeting senescent hepatocytes using the thrombomodulin-PAR1 inhibitor vorapaxar ameliorates NAFLD progression.

BACKGROUND AND AIMS: Senescent hepatocytes accumulate in parallel with fibrosis progression during NASH. The mechanisms that enable progressive expansion of nonreplicating cell populations and the significance of that process in determining NASH outcomes are unclear. Senescing cells upregulate thrombomodulin-protease-activated receptor-1 (THBD-PAR1) signaling to remain viable. Vorapaxar blocks the activity of that pathway. We used vorapaxar to determine if and how THBD-PAR1 signaling promotes fibrosis progression in NASH. APPROACH AND RESULTS: We evaluated the THBD-PAR1 pathway in liver biopsies from patients with NAFLD. Chow-fed mice were treated with viral vectors to overexpress p16 in hepatocytes and induce replicative senescence. Effects on the THBD-PAR1 axis and regenerative capacity were assessed; the transcriptome of p16-overexpressing hepatocytes was characterized, and we examined how conditioned medium from senescent but viable (dubbed "undead") hepatocytes reprograms HSCs. Mouse models of NASH caused by genetic obesity or Western diet/CCl 4 were treated with vorapaxar to determine effects on hepatocyte senescence and liver damage. Inducing senescence upregulates the THBD-PAR1 signaling axis in hepatocytes and induces their expression of fibrogenic factors, including hedgehog ligands. Hepatocyte THBD-PAR1 signaling increases in NAFLD and supports sustained hepatocyte senescence that limits effective liver regeneration and promotes maladaptive repair. Inhibiting PAR1 signaling with vorapaxar interrupts this process, reduces the burden of 'undead' senescent cells, and safely improves NASH and fibrosis despite ongoing lipotoxic stress. CONCLUSION: The THBD-PAR1 signaling axis is a novel therapeutic target for NASH because blocking this pathway prevents accumulation of senescing but viable hepatocytes that generate factors that promote maladaptive liver repair.

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.

Human genetic variation in VAC14 regulates Salmonella invasion and typhoid fever through modulation of cholesterol.

Risk, severity, and outcome of infection depend on the interplay of pathogen virulence and host susceptibility. Systematic identification of genetic susceptibility to infection is being undertaken through genome-wide association studies, but how to expeditiously move from genetic differences to functional mechanisms is unclear. Here, we use genetic association of molecular, cellular, and human disease traits and experimental validation to demonstrate that genetic variation affects expression of VAC14, a phosphoinositide-regulating protein, to influence susceptibility to Salmonella enterica serovar Typhi (S Typhi) infection. Decreased VAC14 expression increased plasma membrane cholesterol, facilitating Salmonella docking and invasion. This increased susceptibility at the cellular level manifests as increased susceptibility to typhoid fever in a Vietnamese population. Furthermore, treating zebrafish with a cholesterol-lowering agent, ezetimibe, reduced susceptibility to S Typhi. Thus, coupling multiple genetic association studies with mechanistic dissection revealed how VAC14 regulates Salmonella invasion and typhoid fever susceptibility and may open doors to new prophylactic/therapeutic approaches.

Methylthioadenosine Suppresses Salmonella Virulence.

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.

Multiple Surface Regions on the Niemann-Pick C2 Protein Facilitate Intracellular Cholesterol Transport.

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.

Calcineurin orchestrates dimorphic transitions, antifungal drug responses and host-pathogen interactions of the pathogenic mucoralean fungus Mucor circinelloides.

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.

Functional genetic screen of human diversity reveals that a methionine salvage enzyme regulates inflammatory cell death.

Genome-wide association studies can identify common differences that contribute to human phenotypic diversity and disease. When genome-wide association studies are combined with approaches that test how variants alter physiology, biological insights can emerge. Here, we used such an approach to reveal regulation of cell death by the methionine salvage pathway. A common SNP associated with reduced expression of a putative methionine salvage pathway dehydratase, apoptotic protease activating factor 1 (APAF1)-interacting protein (APIP), was associated with increased caspase-1-mediated cell death in response to Salmonella. The role of APIP in methionine salvage was confirmed by growth assays with methionine-deficient media and quantitation of the methionine salvage substrate, 5'-methylthioadenosine. Reducing expression of APIP or exogenous addition of 5'-methylthioadenosine increased Salmonellae-induced cell death. Consistent with APIP originally being identified as an inhibitor of caspase-9-dependent apoptosis, the same allele was also associated with increased sensitivity to the chemotherapeutic agent carboplatin. Our results show that common human variation affecting expression of a single gene can alter susceptibility to two distinct cell death programs. Furthermore, the same allele that promotes cell death is associated with improved survival of individuals with systemic inflammatory response syndrome, suggesting a possible evolutionary pressure that may explain the geographic pattern observed for the frequency of this SNP. Our study shows that in vitro association screens of disease-related traits can not only reveal human genetic differences that contribute to disease but also provide unexpected insights into cell biology.

Biological sex affects functional variation across the human genome.

Humans display sexual dimorphism across many traits, but little is known about underlying genetic mechanisms and impacts on disease. We utilized single-cell RNA-seq of 480 lymphoblastoid cell lines to reveal that the vast majority (79%) of sex-biased genes are targets of transcription factors that display sex-biased expression. Further, we developed a two-step regression method that identified sex-biased expression quantitative trait loci (sb-eQTL) across the genome. In contrast to previous work, these sb-eQTL are abundant (n=10,754; FDR 5%) and reproducible (replication up to π1=0.56). These sb-eQTL are enriched in over 600 GWAS phenotypes, including 120 sb-eQTL associated with the female-biased autoimmune disease multiple sclerosis. Our results demonstrate widespread genetic impacts on sexual dimorphism and identify possible mechanisms and clinical targets for sex differences in diverse diseases.

Aging promotes metabolic dysfunction-associated steatotic liver disease by inducing ferroptotic stress.

Susceptibility to the biological consequences of aging varies among organs and individuals. We analyzed hepatocyte transcriptomes of healthy young and aged male mice to generate an aging hepatocyte gene signature, used it to deconvolute transcriptomic data from humans and mice with metabolic dysfunction-associated liver disease, validated findings with functional studies in mice and applied the signature to transcriptomic data from other organs to determine whether aging-sensitive degenerative mechanisms are conserved. We discovered that the signature enriches in diseased livers in parallel with degeneration. It is also enriched in failing human hearts, diseased kidneys and pancreatic islets from individuals with diabetes. The signature includes genes that control ferroptosis. Aged mice develop more hepatocyte ferroptosis and liver degeneration than young mice when fed diets that induce metabolic stress. Inhibiting ferroptosis shifts the liver transcriptome of old mice toward that of young mice and reverses aging-exacerbated liver damage, identifying ferroptosis as a tractable, conserved mechanism for aging-related tissue degeneration.

A single amino acid in the Salmonella effector SarA/SteE triggers supraphysiological activation of STAT3 for anti-inflammatory target gene expression.

Non-typhoidal Salmonella enterica cause an estimated 1 million cases of gastroenteritis annually in the United States. These serovars use secreted protein effectors to mimic and reprogram host cellular functions. We previously discovered that the secreted effector SarA (Salmonella anti-inflammatory response activator; also known as SteE) was required for increased intracellular replication of S. Typhimurium and production of the anti-inflammatory cytokine interleukin-10 (IL-10). SarA facilitates phosphorylation of STAT3 through a region of homology with the host cytokine receptor gp130. Here, we demonstrate that a single amino acid difference between SarA and gp130 is critical for the anti-inflammatory bias of SarA-STAT3 signaling. An isoleucine at the pY+1 position of the YxxQ motif in SarA (which binds the SH2 domain in STAT3) causes increased STAT3 phosphorylation and expression of anti-inflammatory target genes. This isoleucine, completely conserved in ~4000 Salmonella isolates, renders SarA a better substrate for tyrosine phosphorylation by GSK-3. GSK-3 is canonically a serine/threonine kinase that nonetheless undergoes tyrosine autophosphorylation at a motif that has an invariant isoleucine at the pY+1 position. Our results provide a molecular basis for how a Salmonella secreted effector achieves supraphysiological levels of STAT3 activation to control host genes during infection.

Nuclear factor kappa B-dependent persistence of Salmonella Typhi and Paratyphi in human macrophages.

Salmonella serovars Typhi and Paratyphi cause a prolonged illness known as enteric fever, whereas other serovars cause acute gastroenteritis. Mechanisms responsible for the divergent clinical manifestations of nontyphoidal and enteric fever Salmonella infections have remained elusive. Here, we show that S. Typhi and S. Paratyphi A can persist within human macrophages, whereas S. Typhimurium rapidly induces apoptotic macrophage cell death that is dependent on Salmonella pathogenicity island 2 (SPI2). S. Typhi and S. Paratyphi A lack 12 specific SPI2 effectors with pro-apoptotic functions, including nine that target nuclear factor κB (NF-κB). Pharmacologic inhibition of NF-κB or heterologous expression of the SPI2 effectors GogA or GtgA restores apoptosis of S. Typhi-infected macrophages. In addition, the absence of the SPI2 effector SarA results in deficient signal transducer and activator of transcription 1 (STAT1) activation and interleukin 12 production, leading to impaired TH1 responses in macrophages and humanized mice. The absence of specific nontyphoidal SPI2 effectors may allow S. Typhi and S. Paratyphi A to cause chronic infections. IMPORTANCE: Salmonella enterica is a common cause of gastrointestinal infections worldwide. The serovars Salmonella Typhi and Salmonella Paratyphi A cause a distinctive systemic illness called enteric fever, whose pathogenesis is incompletely understood. Here, we show that enteric fever Salmonella serovars lack 12 specific virulence factors possessed by nontyphoidal Salmonella serovars, which allow the enteric fever serovars to persist within human macrophages. We propose that this fundamental difference in the interaction of Salmonella with human macrophages is responsible for the chronicity of typhoid and paratyphoid fever, suggesting that targeting the nuclear factor κB (NF-κB) complex responsible for macrophage survival could facilitate the clearance of persistent bacterial infections.

Context-specific eQTLs reveal causal genes underlying shared genetic architecture of critically ill COVID-19 and idiopathic pulmonary fibrosis.

Most genetic variants identified through genome-wide association studies (GWAS) are suspected to be regulatory in nature, but only a small fraction colocalize with expression quantitative trait loci (eQTLs, variants associated with expression of a gene). Therefore, it is hypothesized but largely untested that integration of disease GWAS with context-specific eQTLs will reveal the underlying genes driving disease associations. We used colocalization and transcriptomic analyses to identify shared genetic variants and likely causal genes associated with critically ill COVID-19 and idiopathic pulmonary fibrosis. We first identified five genome-wide significant variants associated with both diseases. Four of the variants did not demonstrate clear colocalization between GWAS and healthy lung eQTL signals. Instead, two of the four variants colocalized only in cell-type and disease-specific eQTL datasets. These analyses pointed to higher ATP11A expression from the C allele of rs12585036, in monocytes and in lung tissue from primarily smokers, which increased risk of IPF and decreased risk of critically ill COVID-19. We also found lower DPP9 expression (and higher methylation at a specific CpG) from the G allele of rs12610495, acting in fibroblasts and in IPF lungs, and increased risk of IPF and critically ill COVID-19. We further found differential expression of the identified causal genes in diseased lungs when compared to non-diseased lungs, specifically in epithelial and immune cell types. These findings highlight the power of integrating GWAS, context-specific eQTLs, and transcriptomics of diseased tissue to harness human genetic variation to identify causal genes and where they function during multiple 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.

The RNA helicase DDX39B activates FOXP3 RNA splicing to control T regulatory cell fate.

Genes 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 thereby 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 Forkhead Box 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 type 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.

A genome-wide in vitro bacterial-infection screen reveals human variation in the host response associated with inflammatory disease.

Recent progress in cataloguing common genetic variation has made possible genome-wide studies that are beginning to elucidate the causes and consequences of our genetic differences. Approaches that provide a mechanistic understanding of how genetic variants function to alter disease susceptibility and why they were substrates of natural selection would complement other approaches to human-genome analysis. Here we use a novel cell-based screen of bacterial infection to identify human variation in Salmonella-induced cell death. A loss-of-function allele of CARD8, a reported inhibitor of the proinflammatory protease caspase-1, was associated with increased cell death in vitro (p = 0.013). The validity of this association was demonstrated through overexpression of alternative alleles and RNA interference in cells of varying genotype. Comparison of mammalian CARD8 orthologs and examination of variation among different human populations suggest that the increase in infectious-disease burden associated with larger animal groups (i.e., herds and colonies), and possibly human population expansion, may have naturally selected for loss of CARD8. We also find that the loss-of-function CARD8 allele shows a modest association with an increased risk of systemic inflammatory response syndrome in a small study (p = 0.05). Therefore, a by-product of the selected benefit of loss of CARD8 could be increased inflammatory diseases. These results demonstrate the utility of genome-wide cell-based association screens with microbes in the identification of naturally selected variants that can impact human health.