PRMT1 acts as a suppressor of MHC-I and anti-tumor immunity

Cancer immunotherapies have demonstrated remarkable success; however, the majority of patients do not respond or develop resistance. Here, we conduct epigenetic gene-targeted CRISPR-Cas9 screens to identify epigenomic factors that limit CD8+ T cell-mediated anti-tumor immunity. We identify that PRMT1 suppresses interferon gamma (Ifnγ)-induced MHC-I expression, thus dampening CD8+ T cell-mediated killing. Indeed, PRMT1 knockout or pharmacological targeting of type I PRMT with the clinical inhibitor GSK3368715 enhances Ifnγ-induced MHC-I expression through elevated STAT1 expression and activation, while re-introduction of PRMT1 in PRMT1-deficient cells reverses this effect. Importantly, loss of PRMT1 enhances the efficacy of anti-PD-1 immunotherapy, and The Cancer Genome Atlas analysis reveals that PRMT1 expression in human melanoma is inversely correlated with expression of human leukocyte antigen molecules, infiltration of CD8+ T cells, and overall survival. Taken together, we identify PRMT1 as a negative regulator of anti-tumor immunity, unveiling clinical type I PRMT inhibitors as immunotherapeutic agents or as adjuncts to existing immunotherapies.

Abstract P1-04-04: Dna barcoding reveals ongoing immunoediting of clonal cancer populations during metastatic progression and in response to immunotherapy

As cancers develop and spread they must continually evade immune destruction. Understanding mechanisms of immune evasion in cancer is clinically significant as demonstrated with the successes of immune checkpoint inhibitors. Breast cancer is known to be highly immune evasive and responds poorly to the current immunotherapies, indicating alternative immune pathways must be targeted. We hypothesise that there are unidentified genetic mechanisms that enable immune evasion in breast cancer. We aim to uncover and target these mechanisms to sensitise immune evasive breast cancer cells to immune destruction in the context of immunotherapy treatment. DNA barcoding technology offers a new approach to understanding immune evasion. By stably integrating a unique DNA barcode sequence into each cell, we can study clonal immune evasion in vivo. Using this technology, we identified cancer cell clones from the 4T1 murine mammary carcinoma cell line that are highly enriched in lung metastases following treatment with combination immunotherapy (anti-CTLA-4 plus anti-PD-1). We isolated these specific immune evasive clones and established them as clonal cell lines. We have identified stark clonal differences in both PD-L1 and MHC I expression at both the RNA and protein level, and shown that MHC I expression is only partially controlled by epigenetic mechanisms. In addition, immune evasive subclones co-cultured with stimulated T cells resulted in less activated T cells than their less evasive counterparts. Furthermore, RNA sequencing of these clones has identified a gene signature that is strongly associated with decreased survival in both the METABRIC and TCGA cohorts. We have demonstrated ongoing immunoediting in the 4T1 model in vivo, both during metastasis and immunotherapy treatment. We have also identified subclonal populations of cells within a single tumour utilising different mechanisms of immune evasion. RNA sequencing has revealed a gene signature strongly associated with poor survival of basal-like breast cancer in two cohorts. Further pathway-level analysis of the resulting gene signature is required to elucidate the drivers of this aggressive and immune evasive phenotype. By targeting newly identified mechanisms of immune evasion in combination with current immunotherapies, we hope to improve the long-term survival of breast cancer patients.

Citation Format: Louise A Baldwin, Nenad Bartonicek, Jessica Yang, Sunny Z Wu, Niantao Deng, Daniel Roden, Chia-Ling Chan, Ghamdan Al-Eryani, Damien J Zanker, Belinda S Parker, Alexander Swarbrick, Simon Junankar. Dna barcoding reveals ongoing immunoediting of clonal cancer populations during metastatic progression and in response to immunotherapy [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-04-04.

Abstract 2761: CODEX highly multiplex image mapping to CITEseq datasets reveal the spatial dynamics of the TME during the development of acquired resistant in immunotherapy treated melanoma

Primary and acquired resistance to anti-PD-1 and CTLA-4 immunotherapies remains a major obstacle to improving the outcomes of cancer patients. This study utilized a combination of scRNAseq, CITE-seq, a method for CODEX, 40+ protein marker tissue imaging to characterize the molecular and spatial features of cells within melanoma biopsies of patients receiving anti-PD-1 and CTLA-4 immunotherapies. Analysis was conducted on melanoma biopsies taken pre-treatment, early during treatment (7-14 days) and after disease progression. Five patients underwent all analyses; three were responding, two had primary resistance and one had acquired resistance. Tissue dissociates underwent ssRNAseq along with CITE-seq, which incorporated 153+ antibodies to extensively define tumor and immune cell phenotypes. A separate portion of the same lesion underwent formalin fixation and paraffin embedding and was then analyzed using the CODEX high-multiplex imaging system to spatially localize 40+ markers in the tumor. The CODEX imaging system and processor (Akoya Biosciences) were used to generate multiplex images. HALO AI (Indica Labs) was used to perform cell segmentation and quantify single-cell marker expression of the 40+ markers. Tumor areas were classified into vasculature (CD34+/CD31+), melanoma (SOX10+) and stroma regions using the DenseNet deep learning classifier within HALO AI. The subsequent output contained single-cell expression, tissue classification and x and y location data, which were used to map to the CITE-seq data using the CamaraLab STvEA (Govek et al, biorix 2019) to identify CITE-seq defined populations within tissue sections. CODEX analysis of 2.2M cells revealed distinct immune phenotypes of tumor infiltrated cells compared to their peritumor counterparts, with a gradient in the upregulation of immunosuppressive and exhaustive signaling upon entering the tumor region. Changes in cellular populations during treatment were significantly related to the spatial location of the cells within the tumor. Treatment with immunotherapies causes a dynamic change in these profiles which relates to the development of tumor resistance.

Citation Format: Camelia Quek, Ines P. Silva, Ghamdan Al-Eryani, Aaron Mayer, Nenad Bartonicek, Kate Harvey, Chi Buckley, Oliver Braubach, John F. Thompson, Jonathan R. Stretch, Robyn P M Saw, Kerwin F. Shannon, Alexander M. Menzies, Richard A. Scolyer, Georgina V. Long, Alexander Swarbrick, James S. Wilmott. CODEX highly multiplex image mapping to CITEseq datasets reveal the spatial dynamics of the TME during the development of acquired resistant in immunotherapy treated melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2761.

Cryopreservation of human cancers conserves tumour heterogeneity for single-cell multi-omics analysis

BACKGROUND

High throughput single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for exploring cellular heterogeneity among complex human cancers. scRNA-Seq studies using fresh human surgical tissue are logistically difficult, preclude histopathological triage of samples, and limit the ability to perform batch processing. This hindrance can often introduce technical biases when integrating patient datasets and increase experimental costs. Although tissue preservation methods have been previously explored to address such issues, it is yet to be examined on complex human tissues, such as solid cancers and on high throughput scRNA-Seq platforms.

METHODS

Using the Chromium 10X platform, we sequenced a total of ~ 120,000 cells from fresh and cryopreserved replicates across three primary breast cancers, two primary prostate cancers and a cutaneous melanoma. We performed detailed analyses between cells from each condition to assess the effects of cryopreservation on cellular heterogeneity, cell quality, clustering and the identification of gene ontologies. In addition, we performed single-cell immunophenotyping using CITE-Seq on a single breast cancer sample cryopreserved as solid tissue fragments.

RESULTS

Tumour heterogeneity identified from fresh tissues was largely conserved in cryopreserved replicates. We show that sequencing of single cells prepared from cryopreserved tissue fragments or from cryopreserved cell suspensions is comparable to sequenced cells prepared from fresh tissue, with cryopreserved cell suspensions displaying higher correlations with fresh tissue in gene expression. We showed that cryopreservation had minimal impacts on the results of downstream analyses such as biological pathway enrichment. For some tumours, cryopreservation modestly increased cell stress signatures compared to freshly analysed tissue. Further, we demonstrate the advantage of cryopreserving whole-cells for detecting cell-surface proteins using CITE-Seq, which is impossible using other preservation methods such as single nuclei-sequencing.

CONCLUSIONS

We show that the viable cryopreservation of human cancers provides high-quality single-cells for multi-omics analysis. Our study guides new experimental designs for tissue biobanking for future clinical single-cell RNA sequencing studies.

Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

Leptin antagonism inhibits prostate cancer xenograft growth and progression

Hyperleptinaemia is a well-established therapeutic side effect of drugs inhibiting the androgen axis in prostate cancer (PCa), including main stay androgen deprivation therapy (ADT) and androgen targeted therapies (ATT). Given significant crossover between the adipokine hormone signalling of leptin and multiple cancer-promoting hallmark pathways, including growth, proliferation, migration, angiogenesis, metabolism and inflammation, targeting the leptin axis is therapeutically appealing, especially in advanced PCa where current therapies fail to be curative. In this study, we uncover leptin as a novel universal target in PCa and are the first to highlight increased intratumoural leptin and leptin receptor (LEPR) expression in PCa cells and patients' tumours exposed to androgen deprivation, as is observed in patients' tumours of metastatic and castrate resistant (CRPC) PCa. We also reveal the world-first preclinical evidence that demonstrates marked efficacy of targeted leptin-signalling blockade, using Allo-aca, a potent, specific, and safe LEPR peptide antagonist. Allo-aca-suppressed tumour growth and delayed progression to CRPC in mice bearing LNCaP xenografts, with reduced tumour vascularity and altered pathways of apoptosis, transcription/translation, and energetics in tumours determined as potential mechanisms underpinning anti-tumour efficacy. We highlight LEPR blockade in combination with androgen axis inhibition represents a promising new therapeutic strategy vital in advanced PCa treatment.

Adiponectin receptor activation inhibits prostate cancer xenograft growth

Adiponectin is an adipokine originally identified as dysregulated in obesity, with a key role in insulin sensitisation and in maintaining systemic energy balance. However, adiponectin is progressively emerging as having aberrant signalling in multiple disease states, including prostate cancer (PCa). Circulating adiponectin is lower in patients with PCa than in non-malignant disease, and inversely correlates with cancer severity. More severe hypoadiponectinemia is observed in advanced PCa than in organ-confined disease. Given the crossover between adiponectin signalling and several cancer hallmark pathways that influence PCa growth and progression, we hypothesised that targeting dysregulated adiponectin signalling may inhibit tumour growth and progression. We, therefore, aimed to test the efficacy of correcting the hypoadiponectinemia and dysregulated adiponectin signalling observed in PCa, a world-first PCa therapeutic approach, using peptide adiponectin receptor (ADIPOR) agonist ADP355 in mice bearing subcutaneous LNCaP xenografts. We demonstrate significant evidence for PCa growth inhibition by ADP355, which slowed tumour growth and delayed progression of serum PCa biomarker, prostate-specific antigen (PSA), compared to vehicle. ADP355 conferred a significant advantage by increasing time on treatment with a delayed ethical endpoint. mRNA sequencing and protein expression analyses of tumours revealed ADP355 PCa growth inhibition may be through altered cellular energetics, cellular stress and protein synthesis, which may culminate in apoptosis, as evidenced by the increased apoptotic marker in ADP355-treated tumours. Our findings highlight the efficacy of ADP355 in targeting classical adiponectin-associated signalling pathways in vivo and provide insights into the promising future for modulating adiponectin signalling through ADIPOR agonism as a novel anti-tumour treatment modality.

Stromal cell diversity associated with immune evasion in human triple-negative breast cancer

The tumour stroma regulates nearly all stages of carcinogenesis. Stromal heterogeneity in human triple-negative breast cancers (TNBCs) remains poorly understood, limiting the development of stromal-targeted therapies. Single-cell RNA sequencing of five TNBCs revealed two cancer-associated fibroblast (CAF) and two perivascular-like (PVL) subpopulations. CAFs clustered into two states: the first with features of myofibroblasts and the second characterised by high expression of growth factors and immunomodulatory molecules. PVL cells clustered into two states consistent with a differentiated and immature phenotype. We showed that these stromal states have distinct morphologies, spatial relationships and functional properties in regulating the extracellular matrix. Using cell signalling predictions, we provide evidence that stromal-immune crosstalk acts via a diverse array of immunoregulatory molecules. Importantly, the investigation of gene signatures from inflammatory-CAFs and differentiated-PVL cells in independent TNBC patient cohorts revealed strong associations with cytotoxic T-cell dysfunction and exclusion, respectively. Such insights present promising candidates to further investigate for new therapeutic strategies in the treatment of TNBCs.

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.

ELF5 modulates the estrogen receptor cistrome in breast cancer

Acquired resistance to endocrine therapy is responsible for half of the therapeutic failures in the treatment of breast cancer. Recent findings have implicated increased expression of the ETS transcription factor ELF5 as a potential modulator of estrogen action and driver of endocrine resistance, and here we provide the first insight into the mechanisms by which ELF5 modulates estrogen sensitivity. Using chromatin immunoprecipitation sequencing we found that ELF5 binding overlapped with FOXA1 and ER at super enhancers, enhancers and promoters, and when elevated, caused FOXA1 and ER to bind to new regions of the genome, in a pattern that replicated the alterations to the ER/FOXA1 cistrome caused by the acquisition of resistance to endocrine therapy. RNA sequencing demonstrated that these changes altered estrogen-driven patterns of gene expression, the expression of ER transcription-complex members, and 6 genes known to be involved in driving the acquisition of endocrine resistance. Using rapid immunoprecipitation mass spectrometry of endogenous proteins, and proximity ligation assays, we found that ELF5 interacted physically with members of the ER transcription complex, such as DNA-PKcs. We found 2 cases of endocrine-resistant brain metastases where ELF5 levels were greatly increased and ELF5 patterns of gene expression were enriched, compared to the matched primary tumour. Thus ELF5 alters ER-driven gene expression by modulating the ER/FOXA1 cistrome, by interacting with it, and by modulating the expression of members of the ER transcriptional complex, providing multiple mechanisms by which ELF5 can drive endocrine resistance.

The long noncoding RNA lncNB1 promotes tumorigenesis by interacting with ribosomal protein RPL35

The majority of patients with neuroblastoma due to MYCN oncogene amplification and consequent N-Myc oncoprotein over-expression die of the disease. Here our analyses of RNA sequencing data identify the long noncoding RNA lncNB1 as one of the transcripts most over-expressed in MYCN-amplified, compared with MYCN-non-amplified, human neuroblastoma cells and also the most over-expressed in neuroblastoma compared with all other cancers. lncNB1 binds to the ribosomal protein RPL35 to enhance E2F1 protein synthesis, leading to DEPDC1B gene transcription. The GTPase-activating protein DEPDC1B induces ERK protein phosphorylation and N-Myc protein stabilization. Importantly, lncNB1 knockdown abolishes neuroblastoma cell clonogenic capacity in vitro and leads to neuroblastoma tumor regression in mice, while high levels of lncNB1 and RPL35 in human neuroblastoma tissues predict poor patient prognosis. This study therefore identifies lncNB1 and its binding protein RPL35 as key factors for promoting E2F1 protein synthesis, N-Myc protein stability and N-Myc-driven oncogenesis, and as therapeutic targets.

MYCN amplification is common in neuroblastomas. Here, the authors identify a long noncoding RNA, lncNB1 in these cancers and show that it promotes tumorigenesis by binding to ribosomal protein, RPL35 to enhance E2F1 and DEPDC1B protein synthesis, which phosphorylates ERK to stabilise N-Myc.

Droplet-based single cell RNAseq tools: a practical guide

Droplet based scRNA-seq systems such as Drop-seq, inDrop and Chromium 10X have been the catalyst for the wide adoption of high-throughput scRNA-seq technologies in the research laboratory. In order to understand the capabilities of these systems to deeply interrogate biology; here we provide a practical guide through all the steps involved in a typical scRNA-seq experiment. Through comparing and contrasting these three main droplet based systems (and their derivatives), we provide an overview of all critical considerations in obtaining high quality and biologically relevant data. We also discuss the limitations of these systems and how they fit into the emerging field of Genomic Cytometry.

Targeted, High-Resolution RNA Sequencing of Non-coding Genomic Regions Associated With Neuropsychiatric Functions

The human brain is one of the last frontiers of biomedical research. Genome-wide association studies (GWAS) have succeeded in identifying thousands of haplotype blocks associated with a range of neuropsychiatric traits, including disorders such as schizophrenia, Alzheimer's and Parkinson's disease. However, the majority of single nucleotide polymorphisms (SNPs) that mark these haplotype blocks fall within non-coding regions of the genome, hindering their functional validation. While some of these GWAS loci may contain cis-acting regulatory DNA elements such as enhancers, we hypothesized that many are also transcribed into non-coding RNAs that are missing from publicly available transcriptome annotations. Here, we use targeted RNA capture ('RNA CaptureSeq') in combination with nanopore long-read cDNA sequencing to transcriptionally profile 1,023 haplotype blocks across the genome containing non-coding GWAS SNPs associated with neuropsychiatric traits, using post-mortem human brain tissue from three neurologically healthy donors. We find that the majority (62%) of targeted haplotype blocks, including 13% of intergenic blocks, are transcribed into novel, multi-exonic RNAs, most of which are not yet recorded in GENCODE annotations. We validated our findings with short-read RNA-seq, providing orthogonal confirmation of novel splice junctions and enabling a quantitative assessment of the long-read assemblies. Many novel transcripts are supported by independent evidence of transcription including cap analysis of gene expression (CAGE) data and epigenetic marks, and some show signs of potential functional roles. We present these transcriptomes as a preliminary atlas of non-coding transcription in human brain that can be used to connect neurological phenotypes with gene expression.

Lipid Uptake Is an Androgen-Enhanced Lipid Supply Pathway Associated with Prostate Cancer Disease Progression and Bone Metastasis

De novo lipogenesis is a well-described androgen receptor (AR)-regulated metabolic pathway that supports prostate cancer tumor growth by providing fuel, membrane material, and steroid hormone precursor. In contrast, our current understanding of lipid supply from uptake of exogenous lipids and its regulation by AR is limited, and exogenous lipids may play a much more significant role in prostate cancer and disease progression than previously thought. By applying advanced automated quantitative fluorescence microscopy, we provide the most comprehensive functional analysis of lipid uptake in cancer cells to date and demonstrate that treatment of AR-positive prostate cancer cell lines with androgens results in significantly increased cellular uptake of fatty acids, cholesterol, and low-density lipoprotein particles. Consistent with a direct, regulatory role of AR in this process, androgen-enhanced lipid uptake can be blocked by the AR-antagonist enzalutamide, but is independent of proliferation and cell-cycle progression. This work for the first time comprehensively delineates the lipid transporter landscape in prostate cancer cell lines and patient samples by analysis of transcriptomics and proteomics data, including the plasma membrane proteome. We show that androgen exposure or deprivation regulates the expression of multiple lipid transporters in prostate cancer cell lines and tumor xenografts and that mRNA and protein expression of lipid transporters is enhanced in bone metastatic disease when compared with primary, localized prostate cancer. Our findings provide a strong rationale to investigate lipid uptake as a therapeutic cotarget in the fight against advanced prostate cancer in combination with inhibitors of lipogenesis to delay disease progression and metastasis. IMPLICATIONS: Prostate cancer exhibits metabolic plasticity in acquiring lipids from uptake and lipogenesis at different disease stages, indicating potential therapeutic benefit by cotargeting lipid supply.

Abstract 2453: Eradication of neuroblastoma by suppressing the expression of a single noncoding RNA

N-Myc gene amplification occurs in one quarter of human neuroblastoma tissues, and is a marker for poor patient prognosis. We performed RNA sequencing experiments, and identified 5 transcripts, including RP1NB1, which were most considerably differentially expressed between N-Myc gene amplified and nonamplified human neuroblastoma cell lines. Affymetrix microarray studies revealed that DEPD was one of the few genes considerably downregulated in neuroblastoma cells after RP1NB1 depletion. Chromatin immunoprecipitation assays showed that knocking down RP1NB1 expression reduced histone H3 lysine 4 trimethylation, a marker for active gene transcription, at the DEPD gene promoter. Luciferase assays demonstrated that knocking down RP1NB1 decreased DEPD gene promoter activity. Depletion of RP1NB1 or DEPD with two independent siRNAs or shRNAs significantly reduced ERK protein phosphorylation, N-Myc protein phosphorylation at Serine 62, N-Myc protein stabilization, neuroblastoma cell proliferation and survival. Clonogenic assays showed that knocking down RP1NB1 with doxycycline completely abolished colony formation capacity of neuroblastoma cells stably transfected with doxycycline-inducible RP1NB1 shRNAs. Importantly, treatment with doxycycline in mice xenografted with neuroblastoma cells stably transfected with doxycycline-inducible RP1NB1 shRNA led to tumor eradication. In human neuroblastoma tissues from 600 neuroblastoma patients, high levels of RP1NB1 gene expression correlated with DEPD gene expression and poor patient prognosis. In conclusion, this study identifies the novel long noncoding RNA RP1NB1 as an important regulator of N-Myc protein stability and neuroblastoma tumorigenesis.

Citation Format: Andrew E. Tee, Pei Y. Liu, Giorgio Milazzo, Kate M. Hannan, Jesper Maag, Nenad Bartonicek, Renhua Song, Chen C. Jiang, Xu D. Zhang, Murray D. Norris, Michelle Haber, Glenn M. Marshall, Jinyan Li, Jo Vandesompele, John S. Mattick, Pieter Mestdagh, Giovanni Perini, Ross D. Hannan, Marcel E. Dinger, Tao Liu. Eradication of neuroblastoma by suppressing the expression of a single noncoding RNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2453.

Evidence that TLR4 Is Not a Receptor for Saturated Fatty Acids but Mediates Lipid-Induced Inflammation by Reprogramming Macrophage Metabolism

Chronic inflammation is a hallmark of obesity and is linked to the development of numerous diseases. The activation of toll-like receptor 4 (TLR4) by long-chain saturated fatty acids (lcSFAs) is an important process in understanding how obesity initiates inflammation. While experimental evidence supports an important role for TLR4 in obesity-induced inflammation in vivo, via a mechanism thought to involve direct binding to and activation of TLR4 by lcSFAs, several lines of evidence argue against lcSFAs being direct TLR4 agonists. Using multiple orthogonal approaches, we herein provide evidence that while loss-of-function models confirm that TLR4 does, indeed, regulate lcSFA-induced inflammation, TLR4 is not a receptor for lcSFAs. Rather, we show that TLR4-dependent priming alters cellular metabolism, gene expression, lipid metabolic pathways, and membrane lipid composition, changes that are necessary for lcSFA-induced inflammation. These results reconcile previous discordant observations and challenge the prevailing view of TLR4's role in initiating obesity-induced inflammation.

Intergenic disease-associated regions are abundant in novel transcripts

BACKGROUND

Genotyping of large populations through genome-wide association studies (GWAS) has successfully identified many genomic variants associated with traits or disease risk. Unexpectedly, a large proportion of GWAS single nucleotide polymorphisms (SNPs) and associated haplotype blocks are in intronic and intergenic regions, hindering their functional evaluation. While some of these risk-susceptibility regions encompass cis-regulatory sites, their transcriptional potential has never been systematically explored.

RESULTS

To detect rare tissue-specific expression, we employed the transcript-enrichment method CaptureSeq on 21 human tissues to identify 1775 multi-exonic transcripts from 561 intronic and intergenic haploblocks associated with 392 traits and diseases, covering 73.9 Mb (2.2%) of the human genome. We show that a large proportion (85%) of disease-associated haploblocks express novel multi-exonic non-coding transcripts that are tissue-specific and enriched for GWAS SNPs as well as epigenetic markers of active transcription and enhancer activity. Similarly, we captured transcriptomes from 13 melanomas, targeting nine melanoma-associated haploblocks, and characterized 31 novel melanoma-specific transcripts that include fusion proteins, novel exons and non-coding RNAs, one-third of which showed allelically imbalanced expression.

CONCLUSIONS

This resource of previously unreported transcripts in disease-associated regions ( http://gwas-captureseq.dingerlab.org ) should provide an important starting point for the translational community in search of novel biomarkers, disease mechanisms, and drug targets.

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.

Identification of a novel fusion transcript between human relaxin-1 (RLN1) and human relaxin-2 (RLN2) in prostate cancer

Simultaneous expression of highly homologous RLN1 and RLN2 genes in prostate impairs their accurate delineation. We used PacBio SMRT sequencing and RNA-Seq in LNCaP cells in order to dissect the expression of RLN1 and RLN2 variants. We identified a novel fusion transcript comprising the RLN1 and RLN2 genes and found evidence of its expression in the normal and prostate cancer tissues. The RLN1-RLN2 fusion putatively encodes RLN2 isoform with the deleted secretory signal peptide. The identification of the fusion transcript provided information to determine unique RLN1-RLN2 fusion and RLN1 regions. The RLN1-RLN2 fusion was co-expressed with RLN1 in LNCaP cells, but the two gene products were inversely regulated by androgens. We showed that RLN1 is underrepresented in common PCa cell lines in comparison to normal and PCa tissue. The current study brings a highly relevant update to the relaxin field, and will encourage further studies of RLN1 and RLN2 in PCa and broader.

lncRNAdb v2.0: expanding the reference database for functional long noncoding RNAs

Despite the prevalence of long noncoding RNA (lncRNA) genes in eukaryotic genomes, only a small proportion have been examined for biological function. lncRNAdb, available at http://lncrnadb.org, provides users with a comprehensive, manually curated reference database of 287 eukaryotic lncRNAs that have been described independently in the scientific literature. In addition to capturing a great proportion of the recent literature describing functions for individual lncRNAs, lncRNAdb now offers an improved user interface enabling greater accessibility to sequence information, expression data and the literature. The new features in lncRNAdb include the integration of Illumina Body Atlas expression profiles, nucleotide sequence information, a BLAST search tool and easy export of content via direct download or a REST API. lncRNAdb is now endorsed by RNAcentral and is in compliance with the International Nucleotide Sequence Database Collaboration.

Kraken: a set of tools for quality control and analysis of high-throughput sequence data

New sequencing technologies pose significant challenges in terms of data complexity and magnitude. It is essential that efficient software is developed with performance that scales with this growth in sequence information. Here we present a comprehensive and integrated set of tools for the analysis of data from large scale sequencing experiments. It supports adapter detection and removal, demultiplexing of barcodes, paired-end data, a range of read architectures and the efficient removal of sequence redundancy. Sequences can be trimmed and filtered based on length, quality and complexity. Quality control plots track sequence length, composition and summary statistics with respect to genomic annotation. Several use cases have been integrated into a single streamlined pipeline, including both mRNA and small RNA sequencing experiments. This pipeline interfaces with existing tools for genomic mapping and differential expression analysis.

miR-221 affects multiple cancer pathways by modulating the level of hundreds messenger RNAs

microRNA miR-221 is frequently over-expressed in a variety of human neoplasms. Aim of this study was to identify new miR-221 gene targets to improve our understanding on the molecular tumor-promoting mechanisms affected by miR-221. Gene expression profiling of miR-221-transfected-SNU-398 cells was analyzed by the Sylamer algorithm to verify the enrichment of miR-221 targets among down-modulated genes. This analysis revealed that enforced expression of miR-221 in SNU-398 cells caused the down-regulation of 602 mRNAs carrying sequences homologous to miR-221 seed sequence within their 3'UTRs. Pathways analysis performed on these genes revealed their prominent involvement in cell proliferation and apoptosis. Activation of E2F, MYC, NFkB, and β-catenin pathways was experimentally proven. Some of the new miR-221 target genes, including RB1, WEE1 (cell cycle inhibitors), APAF1 (pro-apoptotic), ANXA1, CTCF (transcriptional repressor), were individually validated as miR-221 targets in SNU-398, HepG2, and HEK293 cell lines. By identifying a large set of miR-221 gene targets, this study improves our knowledge about miR-221 molecular mechanisms involved in tumorigenesis. The modulation of mRNA level of 602 genes confirms the ability of miR-221 to promote cancer by affecting multiple oncogenic pathways.

Extent, causes, and consequences of small RNA expression variation in human adipose tissue

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

Genetic information is transmitted to the cell only through RNA molecules. A special class of RNAs is comprised of the small (up to 30 nucleotide) ones, known to be potent regulators of various cellular processes. At the same time, they have not been as widely studied as messenger RNAs—we do not know how much variation in their sequence and expression level occurs naturally in human populations or how this variability influences other traits. We measured small RNA levels and genetic variability in fat tissue from 131 individuals by high-throughput sequencing. We could associate the expression levels with genetic background of the individuals, as well as changes in metabolic traits. Surprisingly, we found no large scale influence of small RNA variation on mRNA levels, their main regulatory target. Overall, our study is the first to give a quantitative picture of the naturally occurring variation in these important regulatory molecules in human fat tissue.

Targeted deletion of microRNA-22 promotes stress-induced cardiac dilation and contractile dysfunction

BACKGROUND

Delineating the role of microRNAs (miRNAs) in the posttranscriptional gene regulation offers new insights into how the heart adapts to pathological stress. We developed a knockout of miR-22 in mice and investigated its function in the heart.

METHODS AND RESULTS

Here, we show that miR-22-deficient mice are impaired in inotropic and lusitropic response to acute stress by dobutamine. Furthermore, the absence of miR-22 sensitized mice to cardiac decompensation and left ventricular dilation after long-term stimulation by pressure overload. Calcium transient analysis revealed reduced sarcoplasmic reticulum Ca(2+) load in association with repressed sarcoplasmic reticulum Ca(2+) ATPase activity in mutant myocytes. Genetic ablation of miR-22 also led to a decrease in cardiac expression levels for Serca2a and muscle-restricted genes encoding proteins in the vicinity of the cardiac Z disk/titin cytoskeleton. These phenotypes were attributed in part to inappropriate repression of serum response factor activity in stressed hearts. Global analysis revealed increased expression of the transcriptional/translational repressor purine-rich element binding protein B, a highly conserved miR-22 target implicated in the negative control of muscle expression.

CONCLUSION

These data indicate that miR-22 functions as an integrator of Ca(2+) homeostasis and myofibrillar protein content during stress in the heart and shed light on the mechanisms that enhance propensity toward heart failure.

MiR-221 influences effector functions and actin cytoskeleton in mast cells

Mast cells have essential effector and immunoregulatory functions in IgE-associated allergic disorders and certain innate and adaptive immune responses, but the role of miRNAs in regulating mast cell functions is almost completely unexplored. To examine the role of the activation-induced miRNA miR-221 in mouse mast cells, we developed robust lentiviral systems for miRNA overexpression and depletion. While miR-221 favored mast cell adhesion and migration towards SCF or antigen in trans-well migration assays, as well as cytokine production and degranulation in response to IgE-antigen complexes, neither miR-221 overexpression, nor its ablation, interfered with mast cell differentiation. Transcriptional profiling of miR-221-overexpressing mast cells revealed modulation of many transcripts, including several associated with the cytoskeleton; indeed, miR-221 overexpression was associated with reproducible increases in cortical actin in mast cells, and with altered cellular shape and cell cycle in murine fibroblasts. Our bioinformatics analysis showed that this effect was likely mediated by the composite effect of miR-221 on many primary and secondary targets in resting cells. Indeed, miR-221-induced cellular alterations could not be recapitulated by knockdown of one of the major targets of miR-221. We propose a model in which miR-221 has two different roles in mast cells: in resting cells, basal levels of miR-221 contribute to the regulation of the cell cycle and cytoskeleton, a general mechanism probably common to other miR-221-expressing cell types, such as fibroblasts. Vice versa, upon induction in response to mast cell stimulation, miR-221 effects are mast cell-specific and activation-dependent, contributing to the regulation of degranulation, cytokine production and cell adherence. Our studies provide new insights into the roles of miR-221 in mast cell biology, and identify novel mechanisms that may contribute to mast cell-related pathological conditions, such as asthma, allergy and mastocytosis.

The miR-144/451 locus is required for erythroid homeostasis

The process of erythropoiesis must be efficient and robust to supply the organism with red bloods cells both under condition of homeostasis and stress. The microRNA (miRNA) pathway was recently shown to regulate erythroid development. Here, we show that expression of the locus encoding miR-144 and miR-451 is strictly dependent on Argonaute 2 and is required for erythroid homeostasis. Mice deficient for the miR-144/451 cluster display a cell autonomous impairment of late erythroblast maturation, resulting in erythroid hyperplasia, splenomegaly, and a mild anemia. Analysis of gene expression profiles from wild-type and miR-144/451–deficient erythroblasts revealed that the miR-144/451 cluster acts as a “tuner” of gene expression, influencing the expression of many genes. MiR-451 imparts a greater impact on target gene expression than miR-144. Accordingly, mice deficient in miR-451 alone exhibited a phenotype indistinguishable from miR-144/451–deficient mice. Thus, the miR-144/451 cluster tunes gene expression to impart a robustness to erythropoiesis that is critical under conditions of stress.

MADNet: microarray database network web server

MADNet is a user-friendly data mining and visualization tool for rapid analysis of diverse high-throughput biological data such as microarray, phage display or even metagenome experiments. It presents biological information in the context of metabolic and signalling pathways, transcription factors and drug targets through minimal user input, consisting only of the file with the experimental data. These data are integrated with information stored in various biological databases such as NCBI nucleotide and protein databases, metabolic and signalling pathway databases (KEGG), transcription regulation (TRANSFAC©) and drug target database (DrugBank). MADNet is freely available for academic use at http://www.bioinfo.hr/madnet.