The retroelement Lx9 puts a brake on the immune response to virus infection.

The notion that mobile units of nucleic acid known as transposable elements can operate as genomic controlling elements was put forward over six decades ago1,2. However, it was not until the advancement of genomic sequencing technologies that the abundance and repertoire of transposable elements were revealed, and they are now known to constitute up to two-thirds of mammalian genomes3,4. The presence of DNA regulatory regions including promoters, enhancers and transcription-factor-binding sites within transposable elements5-8 has led to the hypothesis that transposable elements have been co-opted to regulate mammalian gene expression and cell phenotype8-14. Mammalian transposable elements include recent acquisitions and ancient transposable elements that have been maintained in the genome over evolutionary time. The presence of ancient conserved transposable elements correlates positively with the likelihood of a regulatory function, but functional validation remains an essential step to identify transposable element insertions that have a positive effect on fitness. Here we show that CRISPR-Cas9-mediated deletion of a transposable element-namely the LINE-1 retrotransposon Lx9c11-in mice results in an exaggerated and lethal immune response to virus infection. Lx9c11 is critical for the neogenesis of a non-coding RNA (Lx9c11-RegoS) that regulates genes of the Schlafen family, reduces the hyperinflammatory phenotype and rescues lethality in virus-infected Lx9c11-/- mice. These findings provide evidence that a transposable element can control the immune system to favour host survival during virus infection.

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.

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.

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.

The endonuclease activity of Mili fuels piRNA amplification that silences LINE1 elements.

Piwi proteins and Piwi-interacting RNAs (piRNAs) have conserved functions in transposon silencing. The murine Piwi proteins Mili and Miwi2 (also called Piwil2 and Piwil4, respectively) direct epigenetic LINE1 and intracisternal A particle transposon silencing during genome reprogramming in the embryonic male germ line. Piwi proteins are proposed to be piRNA-guided endonucleases that initiate secondary piRNA biogenesis; however, the actual contribution of their endonuclease activities to piRNA biogenesis and transposon silencing remain unknown. To investigate the role of Piwi-catalysed endonucleolytic activity, we engineered point mutations in mice that substitute the second aspartic acid to an alanine in the DDH catalytic triad of Mili and Miwi2, generating the Mili(DAH) and Miwi2(DAH) alleles, respectively. Analysis of Mili-bound piRNAs from homozygous Mili(DAH) fetal gonadocytes revealed a failure of transposon piRNA amplification, resulting in the marked reduction of piRNA bound within Miwi2 ribonuclear particles. We find that Mili-mediated piRNA amplification is selectively required for LINE1, but not intracisternal A particle, silencing. The defective piRNA pathway in Mili(DAH) mice results in spermatogenic failure and sterility. Surprisingly, homozygous Miwi2(DAH) mice are fertile, transposon silencing is established normally and no defects in secondary piRNA biogenesis are observed. In addition, the hallmarks of piRNA amplification are observed in Miwi2-deficient gonadocytes. We conclude that cycles of intra-Mili secondary piRNA biogenesis fuel piRNA amplification that is absolutely required for LINE1 silencing.

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.

Long noncoding RNAs in cancer: mechanisms of action and technological advancements.

The previous decade has seen long non-coding RNAs (lncRNAs) rise from obscurity to being defined as a category of genetic elements, leaving its mark on the field of cancer biology. With the current number of curated lncRNAs increasing by 10,000 in the last five years, the field is moving from annotation of lncRNA expression in various tumours to understanding their importance in the key cancer signalling networks and characteristic behaviours. Here, we summarize the previously identified as well as recently discovered mechanisms of lncRNA function and their roles in the hallmarks of cancer. Furthermore, we identify novel technologies for investigation of lncRNA properties and their function in carcinogenesis, which will be important for their translation to the clinic as novel biomarkers and therapeutic targets.

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.

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.

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.

Single-cell spatial multiomics reveals tumor microenvironment vulnerabilities in cancer resistance to immunotherapy.

Heterogeneous resistance to immunotherapy remains a major challenge in cancer treatment, often leading to disease progression and death. Using CITE-seq and matched 40-plex PhenoCycler tissue imaging, we performed longitudinal multimodal single-cell analysis of tumors from metastatic melanoma patients with innate resistance, acquired resistance, or response to immunotherapy. We established the multimodal integration toolkit to align transcriptomic features, cellular epitopes, and spatial information to provide deeper insights into the tumors. With longitudinal analysis, we identified an "immune-striving" tumor microenvironment marked by peri-tumor lymphoid aggregates and low infiltration of T cells in the tumor and the emergence of MITF+SPARCL1+ and CENPF+ melanoma subclones after therapy. The enrichment of B cell-associated signatures in the molecular composition of lymphoid aggregates was associated with better survival. These findings provide further insights into the establishment of microenvironmental cell interactions and molecular composition of spatial structures that could inform therapeutic intervention.

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.

Signaling Switching from Hedgehog-GLI to MAPK Signaling Potentially Serves as a Compensatory Mechanism in Melanoma Cell Lines Resistant to GANT-61.

BACKGROUND: Melanoma represents the deadliest skin cancer due to its cell plasticity which results in high metastatic potential and chemoresistance. Melanomas frequently develop resistance to targeted therapy; therefore, new combination therapy strategies are required. Non-canonical signaling interactions between HH-GLI and RAS/RAF/ERK signaling were identified as one of the drivers of melanoma pathogenesis. Therefore, we decided to investigate the importance of these non-canonical interactions in chemoresistance, and examine the potential for HH-GLI and RAS/RAF/ERK combined therapy. METHODS: We established two melanoma cell lines resistant to the GLI inhibitor, GANT-61, and characterized their response to other HH-GLI and RAS/RAF/ERK inhibitors. RESULTS: We successfully established two melanoma cell lines resistant to GANT-61. Both cell lines showed HH-GLI signaling downregulation and increased invasive cell properties like migration potential, colony forming capacity, and EMT. However, they differed in MAPK signaling activity, cell cycle regulation, and primary cilia formation, suggesting different potential mechanisms responsible for resistance occurrence. CONCLUSIONS: Our study provides the first ever insights into cell lines resistant to GANT-61 and shows potential mechanisms connected to HH-GLI and MAPK signaling which may represent new hot spots for noncanonical signaling interactions.

RNA-seq and ChIP-seq Identification of Unique and Overlapping Targets of GLI Transcription Factors in Melanoma Cell Lines.

BACKGROUND: Despite significant progress in therapy, melanoma still has a rising incidence worldwide, and novel treatment strategies are needed. Recently, researchers have recognized the involvement of the Hedgehog-GLI (HH-GLI) signaling pathway in melanoma and its consistent crosstalk with the MAPK pathway. In order to further investigate the link between the two pathways and to find new target genes that could be considered for combination therapy, we set out to find transcriptional targets of all three GLI proteins in melanoma. METHODS: We performed RNA sequencing on three melanoma cell lines (CHL-1, A375, and MEL224) with overexpressed GLI1, GLI2, and GLI3 and combined them with the results of ChIP-sequencing on endogenous GLI1, GLI2, and GLI3 proteins. After combining these results, 21 targets were selected for validation by qPCR. RESULTS: RNA-seq revealed a total of 808 differentially expressed genes (DEGs) for GLI1, 941 DEGs for GLI2, and 58 DEGs for GLI3. ChIP-seq identified 527 genes that contained GLI1 binding sites in their promoters, 1103 for GLI2 and 553 for GLI3. A total of 15 of these targets were validated in the tested cell lines, 6 of which were detected by both RNA-seq and ChIP-seq. CONCLUSIONS: Our study provides insight into the unique and overlapping transcriptional output of the GLI proteins in melanoma. We suggest that our findings could provide new potential targets to consider while designing melanoma-targeted therapy.

DNA barcoding reveals ongoing immunoediting of clonal cancer populations during metastatic progression and immunotherapy response.

Cancers evade the immune system through the process of cancer immunoediting. While immune checkpoint inhibitors are effective for reactivating tumour immunity in some cancer types, many other solid cancers, including breast cancer, remain largely non-responsive. Understanding how non-responsive cancers evade immunity and whether this occurs at the clonal level will improve immunotherapeutic design. Here we use DNA barcoding to track murine mammary cancer cell clones during immunoediting and determine clonal transcriptional profiles that allow immune evasion following anti-PD1 plus anti-CTLA4 immunotherapy. Clonal diversity is significantly restricted by immunotherapy treatment in both primary tumours and metastases, demonstrating selection for pre-existing breast cancer cell populations and ongoing immunoediting during metastasis and treatment. Immunotherapy resistant clones express a common gene signature associated with poor survival of basal-like breast cancer patient cohorts. At least one of these genes has an existing small molecule that can potentially be used to improve immunotherapy response.

SylArray: a web server for automated detection of miRNA effects from expression data.

UNLABELLED: A useful step for understanding the function of microRNAs (miRNA) or siRNAs is the detection of their effects on genome-wide expression profiles. Typically, approaches look for enrichment of words in the 3(')UTR sequences of the most deregulated genes. A number of tools are available for this purpose, but they require either in-depth computational knowledge, filtered 3(')UTR sequences for the genome of interest, or a set of genes acquired through an arbitrary expression cutoff. To this end, we have developed SylArray; a web-based resource designed for the analysis of large-scale expression datasets. It simply requires the user to submit a sorted list of genes from an expression experiment. SylArray utilizes curated sets of 3(')UTRs to attach sequences to these genes and then applies the Sylamer algorithm for detection of miRNA or siRNA signatures in those sequences. An intuitive system for visualization and interpretation of the small RNA signatures is included. AVAILABILITY: SylArray is written in Perl-CGI, Perl and Java and also uses the R statistical package. The source-code, database and web resource are freely available under GNU Public License (GPL). The web server is freely accessible at http://www.ebi.ac.uk/enright/sylarray.

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.

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.

A single-cell and spatially resolved atlas of human breast cancers.

Breast cancers are complex cellular ecosystems where heterotypic interactions play central roles in disease progression and response to therapy. However, our knowledge of their cellular composition and organization is limited. Here we present a single-cell and spatially resolved transcriptomics analysis of human breast cancers. We developed a single-cell method of intrinsic subtype classification (SCSubtype) to reveal recurrent neoplastic cell heterogeneity. Immunophenotyping using cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) provides high-resolution immune profiles, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcome. Mesenchymal cells displayed diverse functions and cell-surface protein expression through differentiation within three major lineages. Stromal-immune niches were spatially organized in tumors, offering insights into antitumor immune regulation. Using single-cell signatures, we deconvoluted large breast cancer cohorts to stratify them into nine clusters, termed 'ecotypes', with unique cellular compositions and clinical outcomes. This study provides a comprehensive transcriptional atlas of the cellular architecture of breast cancer.

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.