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

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

Human genetic and metabolite variation reveals that methylthioadenosine is a prognostic biomarker and an inflammatory regulator in sepsis.

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.

A cellular genome-wide association study reveals human variation in microtubule stability and a role in inflammatory cell death.

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.

Analysis of Epstein-Barr virus-regulated host gene expression changes through primary B-cell outgrowth reveals delayed kinetics of latent membrane protein 1-mediated NF-κB activation.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that dramatically reorganizes host gene expression to immortalize primary B cells. In this study, we analyzed EBV-regulated host gene expression changes following primary B-cell infection, both during initial proliferation and through transformation into lymphoblastoid cell lines (LCLs). While most EBV-regulated mRNAs were changed during the transition from resting, uninfected B cells through initial B-cell proliferation, a substantial number of mRNAs changed uniquely from early proliferation through LCL outgrowth. We identified constitutively and dynamically EBV-regulated biological processes, protein classes, and targets of specific transcription factors. Early after infection, genes associated with proliferation, stress responses, and the p53 pathway were highly enriched. However, the transition from early to long-term outgrowth was characterized by genes involved in the inhibition of apoptosis, the actin cytoskeleton, and NF-κB activity. It was previously thought that the major viral protein responsible for NF-κB activation, latent membrane protein 1 (LMP1), is expressed within 2 days after infection. Our data indicate that while this is true, LCL-level LMP1 expression and NF-κB activity are not evident until 3 weeks after primary B-cell infection. Furthermore, heterologous NF-κB activation during the first week after infection increased the transformation efficiency, while early NF-κB inhibition had no effect on transformation. Rather, inhibition of NF-κB was not toxic to EBV-infected cells until LMP1 levels and NF-κB activity were high. These data collectively highlight the dynamic nature of EBV-regulated host gene expression and support the notion that early EBV-infected proliferating B cells have a fundamentally distinct growth and survival phenotype from that of LCLs.

SplicerEX: a tool for the automated detection and classification of mRNA changes from conventional and splice-sensitive microarray expression data.

The key postulate that one gene encodes one protein has been overhauled with the discovery that one gene can generate multiple RNA transcripts through alternative mRNA processing. In this study, we describe SplicerEX, a novel and uniquely motivated algorithm designed for experimental biologists that (1) detects widespread changes in mRNA isoforms from both conventional and splice sensitive microarray data, (2) automatically categorizes mechanistic changes in mRNA processing, and (3) mitigates known technological artifacts of exon array-based detection of alternative splicing resulting from 5' and 3' signal attenuation, background detection limits, and saturation of probe set signal intensity. In this study, we used SplicerEX to compare conventional and exon-based Affymetrix microarray data in a model of EBV transformation of primary human B cells. We demonstrated superior detection of 3'-located changes in mRNA processing by the Affymetrix U133 GeneChip relative to the Human Exon Array. SplicerEX-identified exon-level changes in the EBV infection model were confirmed by RT-PCR and revealed a novel set of EBV-regulated mRNA isoform changes in caspases 6, 7, and 8. Finally, SplicerEX as compared with MiDAS analysis of publicly available microarray data provided more efficiently categorized mRNA isoform changes with a significantly higher proportion of hits supported by previously annotated alternative processing events. Therefore, SplicerEX provides an important tool for the biologist interested in studying changes in mRNA isoform usage from conventional or splice-sensitive microarray platforms, especially considering the expansive amount of archival microarray data generated over the past decade. SplicerEX is freely available upon request.

The Epstein-Barr virus (EBV)-induced tumor suppressor microRNA MiR-34a is growth promoting in EBV-infected B cells.

Epstein-Barr virus (EBV) infection of primary human B cells drives their indefinite proliferation into lymphoblastoid cell lines (LCLs). B cell immortalization depends on expression of viral latency genes, as well as the regulation of host genes. Given the important role of microRNAs (miRNAs) in regulating fundamental cellular processes, in this study, we assayed changes in host miRNA expression during primary B cell infection by EBV. We observed and validated dynamic changes in several miRNAs from early proliferation through immortalization; oncogenic miRNAs were induced, and tumor suppressor miRNAs were largely repressed. However, one miRNA described as a p53-targeted tumor suppressor, miR-34a, was strongly induced by EBV infection and expressed in many EBV and Kaposi's sarcoma-associated herpesvirus (KSHV)-infected lymphoma cell lines. EBV latent membrane protein 1 (LMP1) was sufficient to induce miR-34a requiring downstream NF-κB activation but independent of functional p53. Furthermore, overexpression of miR-34a was not toxic in several B lymphoma cell lines, and inhibition of miR-34a impaired the growth of EBV-transformed cells. This study identifies a progrowth role for a tumor-suppressive miRNA in oncogenic-virus-mediated transformation, highlighting the importance of studying miRNA function in different cellular contexts.