Integrative temporal multi-omics reveals uncoupling of transcriptome and proteome during human T cell activation.

Engagement of the T cell receptor (TCR) triggers molecular reprogramming leading to the acquisition of specialized effector functions by CD4 helper and CD8 cytotoxic T cells. While transcription factors, chemokines, and cytokines are known drivers in this process, the temporal proteomic and transcriptomic changes that regulate different stages of human primary T cell activation remain to be elucidated. Here, we report an integrative temporal proteomic and transcriptomic analysis of primary human CD4 and CD8 T cells following ex vivo stimulation with anti-CD3/CD28 beads, which revealed major transcriptome-proteome uncoupling. The early activation phase in both CD4 and CD8 T cells was associated with transient downregulation of the mRNA transcripts and protein of the central glucose transport GLUT1. In the proliferation phase, CD4 and CD8 T cells became transcriptionally more divergent while their proteome became more similar. In addition to the kinetics of proteome-transcriptome correlation, this study unveils selective transcriptional and translational metabolic reprogramming governing CD4 and CD8 T cell responses to TCR stimulation. This temporal transcriptome/proteome map of human T cell activation provides a reference map exploitable for future discovery of biomarkers and candidates targeting T cell responses.

C-reactive protein is a prognostic biomarker in pancreatic ductal adenocarcinoma patients.

AIM: The 5-year survival rate of pancreatic ductal adenocarcinoma (PDAC) is approximately 11% and has only improved marginally over the last three decades. For operable PDAC, resection and adjuvant FOLFIRINOX chemotherapy is standard of care. There is growing interest in perioperative regimens to improve outcomes. The non-randomized Phase II study "Gemcitabine and Abraxane for resectable Pancreatic cancer" (GAP) demonstrated the feasibility of perioperative gemcitabine/abraxane. Long-term survival in PDAC requires an effective immune response; hence, we undertook this translational study of the GAP trial cohort to identify immune-oncology biomarkers for clinical use. METHODS: We combined Nanostring nCounter technology with immunohistochemistry to investigate the correlation between gene expression and overall patient survival. Findings were investigated in samples from the International Cancer Genome Consortium (ICGC, n = 88) and the Australian Pancreatic Genome Initiative (APGI, n = 227). RESULTS: We confirmed that human equilibrative nucleoside transporter 1 (hENT1) expression was not a prognostic marker in PDAC but patients with high levels of hENT1 were more likely to live longer than 24 months post-surgery. Additionally, CD274 (PD-L1) and two novel biomarkers of survival, cathepsin W (CTSW) and C-reactive protein (CRP), were identified in the GAP cohort (n = 19). CRP expression was confirmed in data from the ICGC. Although PD-L1 and CTSW proteins were not significant across all three cohorts, results show that low CRP mRNA and protein expression are associated with longer overall survival in all three patient groups. CONCLUSION: PDAC patients with long survival have higher hENT1 expression levels. Furthermore, CRP expression is a biomarker of poor prognosis following perioperative chemotherapy and resection in PDAC patients and thus may be useful for identifying patients who could benefit from more aggressive adjuvant strategies.

Comprehensive genomic and tumour immune profiling reveals potential therapeutic targets in malignant pleural mesothelioma.

BACKGROUND: Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials. METHODS: We analysed somatic mutations from 229 MPM samples, including previously published data and 58 samples that had undergone WGS within this study. This was combined with RNA-seq analysis to characterize the tumour immune environment. RESULTS: The comprehensive genome analysis identified 12 driver genes, including new candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed SBS5/40 were dominant in 93% of samples, and defects in homologous recombination repair were infrequent in our cohort. The tumour immune environment contained high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immune suppressive environment. The expression of TGFB1 was associated with overall survival. A small subset of samples (less than 10%) had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a 'hot' immune environment independent of the somatic mutations. CONCLUSIONS: We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.

Combined Inhibition of G9a and EZH2 Suppresses Tumor Growth via Synergistic Induction of IL24-Mediated Apoptosis.

G9a and EZH2 are two histone methyltransferases commonly upregulated in several cancer types, yet the precise roles that these enzymes play cooperatively in cancer is unclear. We demonstrate here that frequent concurrent upregulation of both G9a and EZH2 occurs in several human tumors. These methyltransferases cooperatively repressed molecular pathways responsible for tumor cell death. In genetically distinct tumor subtypes, concomitant inhibition of G9a and EZH2 potently induced tumor cell death, highlighting the existence of tumor cell survival dependency at the epigenetic level. G9a and EZH2 synergistically repressed expression of genes involved in the induction of endoplasmic reticulum (ER) stress and the production of reactive oxygen species. IL24 was essential for the induction of tumor cell death and was identified as a common target of G9a and EZH2. Loss of function of G9a and EZH2 activated the IL24-ER stress axis and increased apoptosis in cancer cells while not affecting normal cells. These results indicate that G9a and EZH2 promotes the evasion of ER stress-mediated apoptosis by repressing IL24 transcription, therefore suggesting that their inhibition may represent a potential therapeutic strategy for solid cancers. SIGNIFICANCE: These findings demonstrate a novel role for G9a and EZH2 histone methyltransferases in suppressing apoptosis, which can be targeted with small molecule inhibitors as a potential approach to improve solid cancer treatment.

Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance.

We concurrently examine the whole genome, transcriptome, methylome, and immune cell infiltrates in baseline tumors from 77 patients with advanced cutaneous melanoma treated with anti-PD-1 with or without anti-CTLA-4. We show that high tumor mutation burden (TMB), neoantigen load, expression of IFNγ-related genes, programmed death ligand expression, low PSMB8 methylation (therefore high expression), and T cells in the tumor microenvironment are associated with response to immunotherapy. No specific mutation correlates with therapy response. A multivariable model combining the TMB and IFNγ-related gene expression robustly predicts response (89% sensitivity, 53% specificity, area under the curve [AUC], 0.84); tumors with high TMB and a high IFNγ signature show the best response to immunotherapy. This model validates in an independent cohort (80% sensitivity, 59% specificity, AUC, 0.79). Except for a JAK3 loss-of-function mutation, for patients who did not respond as predicted there is no obvious biological mechanism that clearly explained their outlier status, consistent with intratumor and intertumor heterogeneity in response to immunotherapy.

ERK and mTORC1 Inhibitors Enhance the Anti-Cancer Capacity of the Octpep-1 Venom-Derived Peptide in Melanoma BRAF(V600E) Mutations.

Melanoma is the main cause of skin cancer deaths, with special emphasis in those cases carrying BRAF mutations that trigger the mitogen-activated protein kinases (MAPK) signaling and unrestrained cell proliferation in the absence of mitogens. Current therapies targeting MAPK are hindered by drug resistance and relapse that rely on metabolic rewiring and Akt activation. To identify new drug candidates against melanoma, we investigated the molecular mechanism of action of the Octopus Kaurna-derived peptide, Octpep-1, in human BRAF(V600E) melanoma cells using proteomics and RNAseq coupled with metabolic analysis. Fluorescence microscopy verified that Octpep-1 tagged with fluorescein enters MM96L and NFF cells and distributes preferentially in the perinuclear area of MM96L cells. Proteomics and RNAseq revealed that Octpep-1 targets PI3K/AKT/mTOR signaling in MM96L cells. In addition, Octpep-1 combined with rapamycin (mTORC1 inhibitor) or LY3214996 (ERK1/2 inhibitor) augmented the cytotoxicity against BRAF(V600E) melanoma cells in comparison with the inhibitors or Octpep-1 alone. Octpep-1-treated MM96L cells displayed reduced glycolysis and mitochondrial respiration when combined with LY3214996. Altogether these data support Octpep-1 as an optimal candidate in combination therapies for melanoma BRAF(V600E) mutations.

Whole-genome sequencing of acral melanoma reveals genomic complexity and diversity.

To increase understanding of the genomic landscape of acral melanoma, a rare form of melanoma occurring on palms, soles or nail beds, whole genome sequencing of 87 tumors with matching transcriptome sequencing for 63 tumors was performed. Here we report that mutational signature analysis reveals a subset of tumors, mostly subungual, with an ultraviolet radiation signature. Significantly mutated genes are BRAF, NRAS, NF1, NOTCH2, PTEN and TYRP1. Mutations and amplification of KIT are also common. Structural rearrangement and copy number signatures show that whole genome duplication, aneuploidy and complex rearrangements are common. Complex rearrangements occur recurrently and are associated with amplification of TERT, CDK4, MDM2, CCND1, PAK1 and GAB2, indicating potential therapeutic options.

Hematopoietic stem and progenitor cell-restricted Cdx2 expression induces transformation to myelodysplasia and acute leukemia.

The caudal-related homeobox transcription factor CDX2 is expressed in leukemic cells but not during normal blood formation. Retroviral overexpression of Cdx2 induces AML in mice, however the developmental stage at which CDX2 exerts its effect is unknown. We developed a conditionally inducible Cdx2 mouse model to determine the effects of in vivo, inducible Cdx2 expression in hematopoietic stem and progenitor cells (HSPCs). Cdx2-transgenic mice develop myelodysplastic syndrome with progression to acute leukemia associated with acquisition of additional driver mutations. Cdx2-expressing HSPCs demonstrate enrichment of hematopoietic-specific enhancers associated with pro-differentiation transcription factors. Furthermore, treatment of Cdx2 AML with azacitidine decreases leukemic burden. Extended scheduling of low-dose azacitidine shows greater efficacy in comparison to intermittent higher-dose azacitidine, linked to more specific epigenetic modulation. Conditional Cdx2 expression in HSPCs is an inducible model of de novo leukemic transformation and can be used to optimize treatment in high-risk AML.

Chromatin interactome mapping at 139 independent breast cancer risk signals.

BACKGROUND: Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression. RESULTS: We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region. CONCLUSIONS: Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

Non-coding RNAs underlie genetic predisposition to breast cancer.

BACKGROUND: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. RESULTS: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. CONCLUSIONS: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits.

Exacerbation of Chikungunya Virus Rheumatic Immunopathology by a High Fiber Diet and Butyrate.

Chikungunya virus (CHIKV) is a mosquito transmitted alphavirus associated with a robust systemic infection and an acute inflammatory rheumatic disease. A high fiber diet has been widely promoted for its ability to ameliorate inflammatory diseases. Fiber is fermented in the gut into short chain fatty acids such as acetate, propionate, and butyrate, which enter the circulation providing systemic anti-inflammatory activities. Herein we show that mice fed a high fiber diet show a clear exacerbation of CHIKV arthropathy, with increased edema and neutrophil infiltrates. RNA-Seq analyses illustrated that a high fiber diet, in this setting, promoted a range of pro-neutrophil responses including Th17/IL-17. Gene Set Enrichment Analyses demonstrated significant similarities with mouse models of inflammatory psoriasis and significant depression of macrophage resolution phase signatures in the CHIKV arthritic lesions from mice fed a high fiber diet. Supplementation of the drinking water with butyrate also increased edema after CHIKV infection. However, the mechanisms involved were different, with modulation of AP-1 and NF-κB responses identified, potentially implicating deoptimization of endothelial barrier repair. Thus, neither fiber nor short chain fatty acids provided benefits in this acute infectious disease setting, which is characterized by widespread viral cytopathic effects and a need for tissue repair.

MHC Class II Antigen Presentation by the Intestinal Epithelium Initiates Graft-versus-Host Disease and Is Influenced by the Microbiota.

Graft-versus-host disease (GVHD) in the gastrointestinal (GI) tract is the principal determinant of lethality following allogeneic bone marrow transplantation (BMT). Here, we examined the mechanisms that initiate GVHD, including the relevant antigen-presenting cells. MHC class II was expressed on intestinal epithelial cells (IECs) within the ileum at steady state but was absent from the IECs of germ-free mice. IEC-specific deletion of MHC class II prevented the initiation of lethal GVHD in the GI tract. MHC class II expression on IECs was absent from mice deficient in the TLR adaptors MyD88 and TRIF and required IFNγ secretion by lamina propria lymphocytes. IFNγ responses are characteristically driven by IL-12 secretion from myeloid cells. Antibiotic-mediated depletion of the microbiota inhibited IL-12/23p40 production by ileal macrophages. IL-12/23p40 neutralization prevented MHC class II upregulation on IECs and initiation of lethal GVHD in the GI tract. Thus, MHC class II expression by IECs in the ileum initiates lethal GVHD, and blockade of IL-12/23p40 may represent a readily translatable therapeutic strategy.

Whole-genome landscape of mucosal melanoma reveals diverse drivers and therapeutic targets.

Knowledge of key drivers and therapeutic targets in mucosal melanoma is limited due to the paucity of comprehensive mutation data on this rare tumor type. To better understand the genomic landscape of mucosal melanoma, here we describe whole genome sequencing analysis of 67 tumors and validation of driver gene mutations by exome sequencing of 45 tumors. Tumors have a low point mutation burden and high numbers of structural variants, including recurrent structural rearrangements targeting TERT, CDK4 and MDM2. Significantly mutated genes are NRAS, BRAF, NF1, KIT, SF3B1, TP53, SPRED1, ATRX, HLA-A and CHD8. SF3B1 mutations occur more commonly in female genital and anorectal melanomas and CTNNB1 mutations implicate a role for WNT signaling defects in the genesis of some mucosal melanomas. TERT aberrations and ATRX mutations are associated with alterations in telomere length. Mutation profiles of the majority of mucosal melanomas suggest potential susceptibility to CDK4/6 and/or MEK inhibitors.

Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma.

Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that form intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to develop treatment resistance, ultimately leading to disease recurrence. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogeneous tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. Single nucleotide polymorphism arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors which were consistent with the heterogeneous molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

Different genetic mechanisms mediate spontaneous versus UVR-induced malignant melanoma.

Genetic variation conferring resistance and susceptibility to carcinogen-induced tumorigenesis is frequently studied in mice. We have now turned this idea to melanoma using the collaborative cross (CC), a resource of mouse strains designed to discover genes for complex diseases. We studied melanoma-prone transgenic progeny across seventy CC genetic backgrounds. We mapped a strong quantitative trait locus for rapid onset spontaneous melanoma onset to Prkdc, a gene involved in detection and repair of DNA damage. In contrast, rapid onset UVR-induced melanoma was linked to the ribosomal subunit gene Rrp15. Ribosome biogenesis was upregulated in skin shortly after UVR exposure. Mechanistically, variation in the 'usual suspects' by which UVR may exacerbate melanoma, defective DNA repair, melanocyte proliferation, or inflammatory cell infiltration, did not explain melanoma susceptibility or resistance across the CC. Instead, events occurring soon after exposure, such as dysregulation of ribosome function, which alters many aspects of cellular metabolism, may be important.

The immune checkpoint CD96 defines a distinct lymphocyte phenotype and is highly expressed on tumor-infiltrating T cells.

CD96 has recently been shown to be a potent immune checkpoint molecule in mice, but a similar role in humans is not known. In this study, we provide a detailed map of CD96 expression across human lymphocyte lineages, the kinetics of CD96 regulation on T-cell activation and co-expression with other conventional and emerging immune checkpoint molecules. We show that CD96 is predominantly expressed by T cells and has a unique lymphocyte expression profile. CD96high T cells exhibited distinct effector functions on activation. Of note, CD96 expression was highly correlated with T-cell markers in primary and metastatic human tumors and was elevated on antigen-experienced T cells and tumor-infiltrating lymphocytes. Collectively, these data demonstrate that CD96 may be a promising immune checkpoint to enhance T-cell function against human cancer and infectious disease.

Bone marrow transplantation generates T cell-dependent control of myeloma in mice.

Transplantation with autologous hematopoietic progenitors remains an important consolidation treatment for patients with multiple myeloma (MM) and is thought to prolong the disease plateau phase by providing intensive cytoreduction. However, transplantation induces inflammation in the context of profound lymphodepletion that may cause hitherto unexpected immunological effects. We developed preclinical models of bone marrow transplantation (BMT) for MM using Vk*MYC myeloma-bearing recipient mice and donor mice that were myeloma naive or myeloma experienced to simulate autologous transplantation. Surprisingly, we demonstrated broad induction of T cell-dependent myeloma control, most efficiently from memory T cells within myeloma-experienced grafts, but also through priming of naive T cells after BMT. CD8+ T cells from mice with controlled myeloma had a distinct T cell receptor (TCR) repertoire and higher clonotype overlap relative to myeloma-free BMT recipients. Furthermore, T cell-dependent myeloma control could be adoptively transferred to secondary recipients and was myeloma cell clone specific. Interestingly, donor-derived IL-17A acted directly on myeloma cells expressing the IL-17 receptor to induce a transcriptional landscape that promoted tumor growth and immune escape. Conversely, donor IFN-γ secretion and signaling were critical to protective immunity and were profoundly augmented by CD137 agonists. These data provide new insights into the mechanisms of action of transplantation in myeloma and provide rational approaches to improving clinical outcomes.

Keratinocyte Cytokine Networks Associated with Human Melanocytic Nevus Development.

Melanocytes can group together in nevi, commonly thought to form because of intrinsic somatic mutations involving MAPK pathway activation. However, the role of the microenvironment, in particular keratinocytes, in nevogenesis is rarely studied. Melanocytes proliferate during the hair follicle growth phase and in some basal cell carcinomas, allowing us to construct keratinocyte gene expression clusters correlated with melanocyte activation. We asked whether such correlations are evident in the more subtle context of regulation of melanocyte behavior in normal skin. We considered genes which, when mutated in keratinocytes in mice, lead to nevogenesis. Across the human GTEx normal skin database, their expression was correlated with that of keratinocyte cytokines KITLG, HGF, FGF2, EDN1, and melanocyte markers. These cytokines have pleiotropic effects on melanocyte-specific and pigmentation genes and also influence mast cell gene expression. We show five classes of keratinocyte genes that, via germline genetic variation, influence melanocyte activity. These include genes involved in SHH signaling, structural keratins, ribosomal biogenesis, and stem cell governance. In agreement with the finding of KITLG linked to nevogenesis in human genome-wide association studies, we provide evidence that specific keratinocyte cytokines are components of networks that may drive or exacerbate nevus development.

HLA and KIR Associations of Cervical Neoplasia.

Background: Cervical cancer is the fourth most common cancer in women, and we recently reported human leukocyte antigen (HLA) alleles showing strong associations with cervical neoplasia risk and protection. HLA ligands are recognized by killer immunoglobulin-like receptors (KIRs) expressed on a range of immune cell subsets, governing their proinflammatory activity. We hypothesized that the inheritance of particular HLA-KIR combinations would increase cervical neoplasia risk. Methods: Here, we used HLA and KIR dosages imputed from single-nucleotide polymorphism genotype data from 2143 cervical neoplasia cases and 13858 healthy controls of European decent. Results: The following 4 novel HLA alleles were identified in association with cervical neoplasia, owing to their linkage disequilibrium with known cervical neoplasia-associated HLA-DRB1 alleles: HLA-DRB3*9901 (odds ratio [OR], 1.24; P = 2.49 × 10-9), HLA-DRB5*0101 (OR, 1.29; P = 2.26 × 10-8), HLA-DRB5*9901 (OR, 0.77; P = 1.90 × 10-9), and HLA-DRB3*0301 (OR, 0.63; P = 4.06 × 10-5). We also found that homozygosity of HLA-C1 group alleles is a protective factor for human papillomavirus type 16 (HPV16)-related cervical neoplasia (C1/C1; OR, 0.79; P = .005). This protective association was restricted to carriers of either KIR2DL2 (OR, 0.67; P = .00045) or KIR2DS2 (OR, 0.69; P = .0006). Conclusions: Our findings suggest that HLA-C1 group alleles play a role in protecting against HPV16-related cervical neoplasia, mainly through a KIR-mediated mechanism.