Clonal evolution of a tumor ecosystem depends on different selection pressures that are principally immune and treatment mediated. We integrate RNA-seq, DNA sequencing, TCR-seq and SNP array data ...across multiple regions of liver cancer specimens to map spatio-temporal interactions between cancer and immune cells. We investigate how these interactions reflect intra-tumor heterogeneity (ITH) by correlating regional neo-epitope and viral antigen burden with the regional adaptive immune response. Regional expression of passenger mutations dominantly recruits adaptive responses as opposed to hepatitis B virus and cancer-testis antigens. We detect different clonal expansion of the adaptive immune system in distant regions of the same tumor. An ITH-based gene signature improves single-biopsy patient survival predictions and an expression survey of 38,553 single cells across 7 regions of 2 patients further reveals heterogeneity in liver cancer. These data quantify transcriptomic ITH and how the different components of the HCC ecosystem interact during cancer evolution.
Abstract Background and Aims Agents that induce an immune response against tumors by altering T-cell regulation have increased survival times of patients with advanced-stage tumors, such as melanoma ...or lung cancer. We aimed to characterize molecular features of immune cells that infiltrate hepatocellular carcinomas (HCCs) to determine whether these types of agents might be effective against liver tumors. Methods We analyzed HCC samples from 956 patients. We separated gene expression profiles from tumor, stromal, and immune cells using a non-negative matrix factorization algorithm. We then analyzed the gene expression pattern of inflammatory cells in HCC tumors samples. We correlated expression patterns with the presence of immune cell infiltrates and immune regulatory molecules, determined by pathology and immunohistochemical analyses, in a training set of 228 HCC samples. We validated the correlation in a validation set of 728 tumor samples. Using data from 190 tumors in the Cancer Genome Atlas, we correlated immune cell gene expression profiles with numbers of chromosomal aberrations (based on single-nucleotide polymorphism array) and mutations (exome sequence data). Results We found approximately 25% of HCCs to have markers of an inflammatory response, with high expression levels of the CD274 molecule (PD-L1) and programmed cell death 1 (PD-1), markers of cytolytic activity, and fewer chromosomal aberrations. We called this group of tumors the Immune class. It contained 2 subtypes, characterized by markers of an adaptive T-cell response or exhausted immune response. The exhausted immune response subclass expressed many genes regulated by transforming growth factor beta 1 (TGFB) that mediate immunosuppression. We did not observe any differences in numbers of mutations or expression of tumor antigens between the immune-specific class and other HCCs. Conclusions In an analysis of HCC samples from 956 patients, we found almost 25% to express markers of an inflammatory response. We identified 2 subclasses, characterized by adaptive or exhausted immune responses. These findings indicate that some HCCs might be susceptible to therapeutic agents designed to block the regulatory pathways in T cells, such as PD-L1, PD-1, or TGFB inhibitors.
A major challenge in inflammatory bowel disease (IBD) is the integration of diverse IBD data sets to construct predictive models of IBD. We present a predictive model of the immune component of IBD ...that informs causal relationships among loci previously linked to IBD through genome-wide association studies (GWAS) using functional and regulatory annotations that relate to the cells, tissues, and pathophysiology of IBD. Our model consists of individual networks constructed using molecular data generated from intestinal samples isolated from three populations of patients with IBD at different stages of disease. We performed key driver analysis to identify genes predicted to modulate network regulatory states associated with IBD, prioritizing and prospectively validating 12 of the top key drivers experimentally. This validated key driver set not only introduces new regulators of processes central to IBD but also provides the integrated circuits of genetic, molecular, and clinical traits that can be directly queried to interrogate and refine the regulatory framework defining IBD.
Genome-wide association studies (GWAS) have identified hundreds of cardiometabolic disease (CMD) risk loci. However, they contribute little to genetic variance, and most downstream gene-regulatory ...mechanisms are unknown. We genotyped and RNA-sequenced vascular and metabolic tissues from 600 coronary artery disease patients in the Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task study (STARNET). Gene expression traits associated with CMD risk single-nucleotide polymorphism (SNPs) identified by GWAS were more extensively found in STARNET than in tissue- and disease-unspecific gene-tissue expression studies, indicating sharing of downstream cis-/trans-gene regulation across tissues and CMDs. In contrast, the regulatory effects of other GWAS risk SNPs were tissue-specific; abdominal fat emerged as an important gene-regulatory site for blood lipids, such as for the low-density lipoprotein cholesterol and coronary artery disease risk gene PCSK9. STARNET provides insights into gene-regulatory mechanisms for CMD risk loci, facilitating their translation into opportunities for diagnosis, therapy, and prevention.
Recently proposed tumor fitness measures, based on profiling neoepitopes for reactive viral epitope similarity, have been proposed to predict response to immune checkpoint inhibitors in melanoma and ...small-cell lung cancer. Here we applied these checkpoint based fitness measures to the matched checkpoint treatment naive Cancer Genome Atlas (TCGA) samples where cytolytic activity (CYT) imparts a known survival benefit. We observed no significant survival predictive power beyond that of overall patient tumor mutation burden, and furthermore, found no association between checkpoint based fitness and tumor T-cell infiltration, cytolytic activity, and abundance (tumor infiltrating lymphocyte, TIL, burden). In addition, we investigated the key assumption of viral epitope similarity driving immune response in the hepatitis B virally infected liver cancer TCGA cohort, and uncovered suggestive evidence that tumor neoepitopes actually dominate viral epitopes in putative immunogenicity and plausibly drive immune response and recruitment.
To identify factors that regulate gut microbiota density and the impact of varied microbiota density on health, we assayed this fundamental ecosystem property in fecal samples across mammals, human ...disease, and therapeutic interventions. Physiologic features of the host (carrying capacity) and the fitness of the gut microbiota shape microbiota density. Therapeutic manipulation of microbiota density in mice altered host metabolic and immune homeostasis. In humans, gut microbiota density was reduced in Crohn's disease, ulcerative colitis, and ileal pouch-anal anastomosis. The gut microbiota in recurrent
infection had lower density and reduced fitness that were restored by fecal microbiota transplantation. Understanding the interplay between microbiota and disease in terms of microbiota density, host carrying capacity, and microbiota fitness provide new insights into microbiome structure and microbiome targeted therapeutics.
This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).
Although diabetes results in part from a deficiency of normal pancreatic beta cells, inducing human beta cells to regenerate is difficult. Reasoning that insulinomas hold the "genomic recipe" for ...beta cell expansion, we surveyed 38 human insulinomas to obtain insights into therapeutic pathways for beta cell regeneration. An integrative analysis of whole-exome and RNA-sequencing data was employed to extensively characterize the genomic and molecular landscape of insulinomas relative to normal beta cells. Here, we show at the pathway level that the majority of the insulinomas display mutations, copy number variants and/or dysregulation of epigenetic modifying genes, most prominently in the polycomb and trithorax families. Importantly, these processes are coupled to co-expression network modules associated with cell proliferation, revealing candidates for inducing beta cell regeneration. Validation of key computational predictions supports the concept that understanding the molecular complexity of insulinoma may be a valuable approach to diabetes drug discovery.Diabetes results in part from a deficiency of functional pancreatic beta cells. Here, the authors study the genomic and epigenetic landscapes of human insulinomas to gain insight into possible pathways for therapeutic beta cell regeneration, highlighting epigenetic genes and pathways.
Cancer testis antigens (CTAs) are an extensive gene family with a unique expression pattern restricted to germ cells, but aberrantly reactivated in cancer tissues. Studies indicate that the ...expression (or re-expression) of CTAs within the MAGE-A family is common in hepatocellular carcinoma (HCC). However, no systematic characterization has yet been reported. The aim of this study is to perform a comprehensive profile of CTA de-regulation in HCC and experimentally evaluate the role of MAGEA3 as a driver of HCC progression. The transcriptomic analysis of 44 multi-regionally sampled HCCs from 12 patients identified high intra-tumor heterogeneity of CTAs. In addition, a subset of CTAs was significantly overexpressed in histologically poorly differentiated regions. Further analysis of CTAs in larger patient cohorts revealed high CTA expression related to worse overall survival and several other markers of poor prognosis. Functional analysis of MAGEA3 was performed in human HCC cell lines by gene silencing and in a genetic mouse model by overexpression of MAGEA3 in the liver. Knockdown of MAGEA3 decreased cell proliferation, colony formation and increased apoptosis. MAGEA3 overexpression was associated with more aggressive tumors in vivo. In conclusion MAGEA3 enhances tumor progression and should be considered as a novel therapeutic target in HCC.