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.

Correction: Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

[This corrects the article DOI: 10.1371/journal.pgen.1006866.].

Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

A small percentage of women with cervical HPV infection progress to cervical neoplasia, and the risk factors determining progression are incompletely understood. We sought to define the genetic loci involved in cervical neoplasia and to assess its heritability using unbiased unrelated case/control statistical approaches. We demonstrated strong association of cervical neoplasia with risk and protective HLA haplotypes that are determined by the amino-acids carried at positions 13 and 71 in pocket 4 of HLA-DRB1 and position 156 in HLA-B. Furthermore, 36% (standard error 2.4%) of liability of HPV-associated cervical pre-cancer and cancer is determined by common genetic variants. Women in the highest 10% of genetic risk scores have approximately >7.1% risk, and those in the highest 5% have approximately >21.6% risk, of developing cervical neoplasia. Future studies should examine genetic risk prediction in assessing the risk of cervical neoplasia further, in combination with other screening methods.

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.

Lower temperatures reduce type I interferon activity and promote alphaviral arthritis.

Chikungunya virus (CHIKV) belongs to a group of mosquito-borne alphaviruses associated with acute and chronic arthropathy, with peripheral and limb joints most commonly affected. Using a mouse model of CHIKV infection and arthritic disease, we show that CHIKV replication and the ensuing foot arthropathy were dramatically reduced when mice were housed at 30°C, rather than the conventional 22°C. The effect was not associated with a detectable fever, but was dependent on type I interferon responses. Bioinformatics analyses of RNA-Seq data after injection of poly(I:C)/jetPEI suggested the unfolded protein response and certain type I interferon responses are promoted when feet are slightly warmer. The ambient temperature thus appears able profoundly to effect anti-viral activity in the periphery, with clear consequences for alphaviral replication and the ensuing arthropathy. These observations may provide an explanation for why alphaviral arthropathies are largely restricted to joints of the limbs and the extremities.

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.

Further assessment of exome-wide UVR footprints in melanoma and their possible relevance.

C>T substitutions at dipyrimidine sites dominate the melanoma genome. We recently analyzed the exomes of spontaneous and neonatal UVR-induced murine melanomas, noting a dramatic change in the genomic footprint at C>T substitutions in the latter. Here we re-analyzed published exome-wide footprints in human melanomas stratified in terms of likely previous sun exposure. Acral and mucosal melanomas were heterogeneous in terms of base substitution types, but most C>Ts occurred in the context of 3'G, probably resulting from spontaneous deamination of the cytosine. C>Ts in sun-exposed melanomas were statistically different from acral/mucosal lesions only in preferring an adjacent 5'T and 3'C. Pyrimidine dimer adducts can form between any pyrimidine (TT, TC, CT, CC). Hence in melanoma C>Ts are overwhelmingly induced at TC or CC photoproducts, or, there are peculiarities in DNA repair that favor the mutation of cytosines with these two pyrimidines adjacent. If melanoma UVR footprints at C>Ts reflect a specific dimer type (eg, 6-4 photoproduct or cyclobutane pyrimidine dimer), these could be removed post UVR, for instance using photolyases, to potentially reduce melanoma risk. If specific modes of DNA repair and/or replication cause these footprints, methodically downregulating selected DNA polymerases in UVR-induced animal models of melanoma, combined with exome sequencing, could begin to assess this. Finally, a preponderance of TpCpC as opposed to NpCpG at C>Ts exome-wide is likely to be a good indicator of whether a melanoma has incurred even a small amount of sun damage. This information will assist epidemiological studies in predicting individual levels of sun exposure.

Interaction of c-Myb with p300 is required for the induction of acute myeloid leukemia (AML) by human AML oncogenes.

The MYB oncogene is widely expressed in acute leukemias and is important for the continued proliferation of leukemia cells, suggesting that MYB may be a therapeutic target in these diseases. However, realization of this potential requires a significant therapeutic window for MYB inhibition, given its essential role in normal hematopoiesis, and an approach for developing an effective therapeutic. We previously showed that the interaction of c-Myb with the coactivator CBP/p300 is essential for its transforming activity. Here, by using cells from Booreana mice which carry a mutant allele of c-Myb, we show that this interaction is essential for in vitro transformation by the myeloid leukemia oncogenes AML1-ETO, AML1-ETO9a, MLL-ENL, and MLL-AF9. We further show that unlike cells from wild-type mice, Booreana cells transduced with AML1-ETO9a or MLL-AF9 retroviruses fail to generate leukemia upon transplantation into irradiated recipients. Finally, we have begun to explore the molecular mechanisms underlying these observations by gene expression profiling. This identified several genes previously implicated in myeloid leukemogenesis and HSC function as being regulated in a c-Myb-p300-dependent manner. These data highlight the importance of the c-Myb-p300 interaction in myeloid leukemogenesis and suggest disruption of this interaction as a potential therapeutic strategy for acute myeloid leukemia.

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.

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.

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.

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.

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.

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.

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.