The NK cell granule protein NKG7 regulates cytotoxic granule exocytosis and inflammation.

Immune-modulating therapies have revolutionized the treatment of chronic diseases, particularly cancer. However, their success is restricted and there is a need to identify new therapeutic targets. Here, we show that natural killer cell granule protein 7 (NKG7) is a regulator of lymphocyte granule exocytosis and downstream inflammation in a broad range of diseases. NKG7 expressed by CD4+ and CD8+ T cells played key roles in promoting inflammation during visceral leishmaniasis and malaria-two important parasitic diseases. Additionally, NKG7 expressed by natural killer cells was critical for controlling cancer initiation, growth and metastasis. NKG7 function in natural killer and CD8+ T cells was linked with their ability to regulate the translocation of CD107a to the cell surface and kill cellular targets, while NKG7 also had a major impact on CD4+ T cell activation following infection. Thus, we report a novel therapeutic target expressed on a range of immune cells with functions in different immune responses.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Whole genomes redefine the mutational landscape of pancreatic cancer.

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Long Noncoding RNAs CUPID1 and CUPID2 Mediate Breast Cancer Risk at 11q13 by Modulating the Response to DNA Damage.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates CCND1. Here, we report that, in addition to regulating CCND1, this enhancer regulates two estrogen-regulated long noncoding RNAs, CUPID1 and CUPID2. We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the CUPID1 and CUPID2 bidirectional promoter. We further show that CUPID1 and CUPID2 are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.

Divergent reprogramming routes lead to alternative stem-cell states.

Pluripotency is defined by the ability of a cell to differentiate to the derivatives of all the three embryonic germ layers: ectoderm, mesoderm and endoderm. Pluripotent cells can be captured via the archetypal derivation of embryonic stem cells or via somatic cell reprogramming. Somatic cells are induced to acquire a pluripotent stem cell (iPSC) state through the forced expression of key transcription factors, and in the mouse these cells can fulfil the strictest of all developmental assays for pluripotent cells by generating completely iPSC-derived embryos and mice. However, it is not known whether there are additional classes of pluripotent cells, or what the spectrum of reprogrammed phenotypes encompasses. Here we explore alternative outcomes of somatic reprogramming by fully characterizing reprogrammed cells independent of preconceived definitions of iPSC states. We demonstrate that by maintaining elevated reprogramming factor expression levels, mouse embryonic fibroblasts go through unique epigenetic modifications to arrive at a stable, Nanog-positive, alternative pluripotent state. In doing so, we prove that the pluripotent spectrum can encompass multiple, unique cell states.

Genome-wide characterization of the routes to pluripotency.

Somatic cell reprogramming to a pluripotent state continues to challenge many of our assumptions about cellular specification, and despite major efforts, we lack a complete molecular characterization of the reprograming process. To address this gap in knowledge, we generated extensive transcriptomic, epigenomic and proteomic data sets describing the reprogramming routes leading from mouse embryonic fibroblasts to induced pluripotency. Through integrative analysis, we reveal that cells transition through distinct gene expression and epigenetic signatures and bifurcate towards reprogramming transgene-dependent and -independent stable pluripotent states. Early transcriptional events, driven by high levels of reprogramming transcription factor expression, are associated with widespread loss of histone H3 lysine 27 (H3K27me3) trimethylation, representing a general opening of the chromatin state. Maintenance of high transgene levels leads to re-acquisition of H3K27me3 and a stable pluripotent state that is alternative to the embryonic stem cell (ESC)-like fate. Lowering transgene levels at an intermediate phase, however, guides the process to the acquisition of ESC-like chromatin and DNA methylation signature. Our data provide a comprehensive molecular description of the reprogramming routes and is accessible through the Project Grandiose portal at http://www.stemformatics.org.

Ssb1 and Ssb2 cooperate to regulate mouse hematopoietic stem and progenitor cells by resolving replicative stress.

Hematopoietic stem and progenitor cells (HSPCs) are vulnerable to endogenous damage and defects in DNA repair can limit their function. The 2 single-stranded DNA (ssDNA) binding proteins SSB1 and SSB2 are crucial regulators of the DNA damage response; however, their overlapping roles during normal physiology are incompletely understood. We generated mice in which both Ssb1 and Ssb2 were constitutively or conditionally deleted. Constitutive Ssb1/Ssb2 double knockout (DKO) caused early embryonic lethality, whereas conditional Ssb1/Ssb2 double knockout (cDKO) in adult mice resulted in acute lethality due to bone marrow failure and intestinal atrophy featuring stem and progenitor cell depletion, a phenotype unexpected from the previously reported single knockout models of Ssb1 or Ssb2 Mechanistically, cDKO HSPCs showed altered replication fork dynamics, massive accumulation of DNA damage, genome-wide double-strand breaks enriched at Ssb-binding regions and CpG islands, together with the accumulation of R-loops and cytosolic ssDNA. Transcriptional profiling of cDKO HSPCs revealed the activation of p53 and interferon (IFN) pathways, which enforced cell cycling in quiescent HSPCs, resulting in their apoptotic death. The rapid cell death phenotype was reproducible in in vitro cultured cDKO-hematopoietic stem cells, which were significantly rescued by nucleotide supplementation or after depletion of p53. Collectively, Ssb1 and Ssb2 control crucial aspects of HSPC function, including proliferation and survival in vivo by resolving replicative stress to maintain genomic stability.

Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339.

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Re-thinking miRNA-mRNA interactions: intertwining issues confound target discovery.

Despite a library full of literature on miRNA biology, core issues relating to miRNA target detection, biological effect, and mode of action remain controversial. This essay proposes that the predominant mechanism of direct miRNA action is translational inhibition, whereas the bulk of miRNA effects are mRNA based. It explores several issues confounding miRNA target detection, and discusses their impact on the dominance of "miRNA seed" dogma and the exploration of non-canonical binding sites. Finally, it makes comparisons between miRNA target prediction and transcription factor binding prediction, and questions the value of characterizing miRNA binding sites based on which miRNA nucleotides are paired with an mRNA.

Minor class splicing shapes the zebrafish transcriptome during development.

Minor class or U12-type splicing is a highly conserved process required to remove a minute fraction of introns from human pre-mRNAs. Defects in this splicing pathway have recently been linked to human disease, including a severe developmental disorder encompassing brain and skeletal abnormalities known as Taybi-Linder syndrome or microcephalic osteodysplastic primordial dwarfism 1, and a hereditary intestinal polyposis condition, Peutz-Jeghers syndrome. Although a key mechanism for regulating gene expression, the impact of impaired U12-type splicing on the transcriptome is unknown. Here, we describe a unique zebrafish mutant, caliban (clbn), with arrested development of the digestive organs caused by an ethylnitrosourea-induced recessive lethal point mutation in the rnpc3 [RNA-binding region (RNP1, RRM) containing 3] gene. rnpc3 encodes the zebrafish ortholog of human RNPC3, also known as the U11/U12 di-snRNP 65-kDa protein, a unique component of the U12-type spliceosome. The biochemical impact of the mutation in clbn is the formation of aberrant U11- and U12-containing small nuclear ribonucleoproteins that impair the efficiency of U12-type splicing. Using RNA sequencing and microarrays, we show that multiple genes involved in various steps of mRNA processing, including transcription, splicing, and nuclear export are disrupted in clbn, either through intron retention or differential gene expression. Thus, clbn provides a useful and specific model of aberrant U12-type splicing in vivo. Analysis of its transcriptome reveals efficient mRNA processing as a critical process for the growth and proliferation of cells during vertebrate development.

MicroRNA-182-5p targets a network of genes involved in DNA repair.

MicroRNAs are noncoding regulators of gene expression, which act by repressing protein translation and/or degrading mRNA. Many have been shown to drive tumorigenesis in cancer, but functional studies to understand their mode of action are typically limited to single-target genes. In this study, we use synthetic biotinylated miRNA to pull down endogenous targets of miR-182-5p. We identified more than 1000 genes as potential targets of miR-182-5p, most of which have a known function in pathways underlying tumor biology. Specifically, functional enrichment analysis identified components of both the DNA damage response pathway and cell cycle to be highly represented in this target cohort. Experimental validation confirmed that miR-182-5p-mediated disruption of the homologous recombination (HR) pathway is a consequence of its ability to target multiple components in that pathway. Although there is a strong enrichment for the cell cycle ontology, we do not see primary proliferative defects as a consequence of miR-182-5p overexpression. We highlight targets that could be responsible for miR-182-5p-mediated disruption of other biological processes attributed in the literature so far. Finally, we show that miR-182-5p is highly expressed in a panel of human breast cancer samples, highlighting its role as a potential oncomir in breast cancer.

miR-139-5p is a regulator of metastatic pathways in breast cancer.

Metastasis is a complex, multistep process involved in the progression of cancer from a localized primary tissue to distant sites, often characteristic of the more aggressive forms of this disease. Despite being studied in great detail in recent years, the mechanisms that govern this process remain poorly understood. In this study, we identify a novel role for miR-139-5p in the inhibition of breast cancer progression. We highlight its clinical relevance by reviewing miR-139-5p expression across a wide variety of breast cancer subtypes using in-house generated and online data sets to show that it is most frequently lost in invasive tumors. A biotin pull-down approach was then used to identify the mRNA targets of miR-139-5p in the breast cancer cell line MCF7. Functional enrichment analysis of the pulled-down targets showed significant enrichment of genes in pathways previously implicated in breast cancer metastasis (P < 0.05). Further bioinformatic analysis revealed a predicted disruption to the TGFß, Wnt, Rho, and MAPK/PI3K signaling cascades, implying a potential role for miR-139-5p in regulating the ability of cells to invade and migrate. To corroborate this finding, using the MDA-MB-231 breast cancer cell line, we show that overexpression of miR-139-5p results in suppression of these cellular phenotypes. Furthermore, we validate the interaction between miR-139-5p and predicted targets involved in these pathways. Collectively, these results suggest a significant functional role for miR-139-5p in breast cancer cell motility and invasion and its potential to be used as a prognostic marker for the aggressive forms of breast cancer.

Addition of interleukin-6 inhibition with tocilizumab to standard graft-versus-host disease prophylaxis after allogeneic stem-cell transplantation: a phase 1/2 trial.

BACKGROUND: Interleukin 6 mediates graft-versus-host disease (GVHD) in experimental allogeneic stem-cell transplantation (allogeneic SCT) and represents an attractive therapeutic target. We aimed to assess whether the humanised anti-interleukin-6 receptor monoclonal antibody, tocilizumab, could attenuate the incidence of acute GVHD. METHODS: We undertook a single-group, single-institution phase 1/2 study at the Royal Brisbane and Women's Hospital Bone Marrow Transplantation unit, QLD, Australia. Eligible patients were 18-65 years old and underwent T-replete HLA-matched allogeneic SCT with either total body irradiation-based myeloablative or reduced-intensity conditioning from unrelated or sibling donors. One intravenous dose of tocilizumab (8 mg/kg, capped at 800 mg, over 60 mins' infusion) was given the day before allogeneic SCT along with standard GVHD prophylaxis (cyclosporin [5 mg/kg per day on days -1 to +1, then 3 mg/kg per day to maintain therapeutic levels (trough levels of 140-300 ng/mL) for 100 days plus methotrexate [15 mg/m(2) on day 1, then 10 mg/m(2) on days 3, 6, and 11]). The primary endpoint was incidence of grade 2-4 acute GVHD at day 100, assessed and graded as per the Seattle criteria. Immunological profiles were compared with a non-randomised group of patients receiving allogeneic SCT, but not treated with tocilizumab. This trial is registered with the Australian and New Zealand Clinical Trials Registry, number ACTRN12612000726853. FINDINGS: Between Jan 19, 2012, and Aug 27, 2013, 48 eligible patients receiving cyclosporin and methotrexate as GVHD prophylaxis were enrolled into the study. The incidence of grade 2-4 acute GVHD in patients treated with tocilizumab at day 100 was 12% (95% CI 5-24), and the incidence of grade 3-4 acute GVHD was 4% (1-13). Grade 2-4 acute GVHD involving the skin developed in five (10%) patients of 48 treated with tocilizumab, involving the gastrointestinal tract in four (8%) patients; there were no reported cases involving the liver. Low incidences of grade 2-4 acute GVHD were noted in patients receiving both myeloablative total body irradiation-based conditioning (12% [95% CI 2-34) and fludarabine and melphalan reduced-intensity conditioning (12% [4-27]). Immune reconstitution was preserved in recipients of interleukin-6 receptor inhibition, but qualitatively modified with suppression of known pathogenic STAT3-dependent pathways. INTERPRETATION: Interleukin 6 is the main detectable and dysregulated cytokine secreted after allogeneic SCT and its inhibition is a potential new and simple strategy to protect from acute GVHD despite robust immune reconstitution; a randomised, controlled trial assessing tocilizumab in addition to standard GVHD prophylaxis in these patients is warranted. FUNDING: National Health and Medical Research Council and Queensland Health.

Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells.

Recent RNA-sequencing studies have shown remarkable complexity in the mammalian transcriptome. The ultimate impact of this complexity on the predicted proteomic output is less well defined. We have undertaken strand-specific RNA sequencing of multiple cellular RNA fractions (>20 Gb) to uncover the transcriptional complexity of human embryonic stem cells (hESCs). We have shown that human embryonic stem (ES) cells display a high degree of transcriptional diversity, with more than half of active genes generating RNAs that differ from conventional gene models. We found evidence that more than 1000 genes express long 5' and/or extended 3'UTRs, which was confirmed by "virtual Northern" analysis. Exhaustive sequencing of the membrane-polysome and cytosolic/untranslated fractions of hESCs was used to identify RNAs encoding peptides destined for secretion and the extracellular space and to demonstrate preferential selection of transcription complexity for translation in vitro. The impact of this newly defined complexity on known gene-centric network models such as the Plurinet and the cell surface signaling machinery in human ES cells revealed a significant expansion of known transcript isoforms at play, many predicting possible alternative functions based on sequence alterations within key functional domains.

Integrated genome analysis suggests that most conserved non-coding sequences are regulatory factor binding sites.

More than 98% of a typical vertebrate genome does not code for proteins. Although non-coding regions are sprinkled with short (<200 bp) islands of evolutionarily conserved sequences, the function of most of these unannotated conserved islands remains unknown. One possibility is that unannotated conserved islands could encode non-coding RNAs (ncRNAs); alternatively, unannotated conserved islands could serve as promoter-distal regulatory factor binding sites (RFBSs) like enhancers. Here we assess these possibilities by comparing unannotated conserved islands in the human and mouse genomes to transcribed regions and to RFBSs, relying on a detailed case study of one human and one mouse cell type. We define transcribed regions by applying a novel transcript-calling algorithm to RNA-Seq data obtained from total cellular RNA, and we define RFBSs using ChIP-Seq and DNAse-hypersensitivity assays. We find that unannotated conserved islands are four times more likely to coincide with RFBSs than with unannotated ncRNAs. Thousands of conserved RFBSs can be categorized as insulators based on the presence of CTCF or as enhancers based on the presence of p300/CBP and H3K4me1. While many unannotated conserved RFBSs are transcriptionally active to some extent, the transcripts produced tend to be unspliced, non-polyadenylated and expressed at levels 10 to 100-fold lower than annotated coding or ncRNAs. Extending these findings across multiple cell types and tissues, we propose that most conserved non-coding genomic DNA in vertebrate genomes corresponds to promoter-distal regulatory elements.

A global role for KLF1 in erythropoiesis revealed by ChIP-seq in primary erythroid cells.

KLF1 regulates a diverse suite of genes to direct erythroid cell differentiation from bipotent progenitors. To determine the local cis-regulatory contexts and transcription factor networks in which KLF1 operates, we performed KLF1 ChIP-seq in the mouse. We found at least 945 sites in the genome of E14.5 fetal liver erythroid cells which are occupied by endogenous KLF1. Many of these recovered sites reside in erythroid gene promoters such as Hbb-b1, but the majority are distant to any known gene. Our data suggests KLF1 directly regulates most aspects of terminal erythroid differentiation including production of alpha- and beta-globin protein chains, heme biosynthesis, coordination of proliferation and anti-apoptotic pathways, and construction of the red cell membrane and cytoskeleton by functioning primarily as a transcriptional activator. Additionally, we suggest new mechanisms for KLF1 cooperation with other transcription factors, in particular the erythroid transcription factor GATA1, to maintain homeostasis in the erythroid compartment.