A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

A cellular hierarchy in melanoma uncouples growth and metastasis.

Although melanoma is notorious for its high degree of heterogeneity and plasticity1,2, the origin and magnitude of cell-state diversity remains poorly understood. Equally, it is unclear whether growth and metastatic dissemination are supported by overlapping or distinct melanoma subpopulations. Here, by combining mouse genetics, single-cell and spatial transcriptomics, lineage tracing and quantitative modelling, we provide evidence of a hierarchical model of tumour growth that mirrors the cellular and molecular logic underlying the cell-fate specification and differentiation of the embryonic neural crest. We show that tumorigenic competence is associated with a spatially localized perivascular niche, a phenotype acquired through an intercellular communication pathway established by endothelial cells. Consistent with a model in which only a fraction of cells are fated to fuel growth, temporal single-cell tracing of a population of melanoma cells with a mesenchymal-like state revealed that these cells do not contribute to primary tumour growth but, instead, constitute a pool of metastatic initiating cells that switch cell identity while disseminating to secondary organs. Our data provide a spatially and temporally resolved map of the diversity and trajectories of melanoma cell states and suggest that the ability to support growth and metastasis are limited to distinct pools of cells. The observation that these phenotypic competencies can be dynamically acquired after exposure to specific niche signals warrant the development of therapeutic strategies that interfere with the cancer cell reprogramming activity of such microenvironmental cues.

Comparative single-cell RNA-sequencing profiling of BMP4-treated primary glioma cultures reveals therapeutic markers.

BACKGROUND: Glioblastoma (GBM) is the most aggressive primary brain tumor. Its cellular composition is very heterogeneous, with cells exhibiting stem-cell characteristics (GSCs) that co-determine therapy resistance and tumor recurrence. Bone Morphogenetic Protein (BMP)-4 promotes astroglial and suppresses oligodendrocyte differentiation in GSCs, processes associated with superior patient prognosis. We characterized variability in cell viability of patient-derived GBM cultures in response to BMP4 and, based on single-cell transcriptome profiling, propose predictive positive and early-response markers for sensitivity to BMP4. METHODS: Cell viability was assessed in 17 BMP4-treated patient-derived GBM cultures. In two cultures, one highly-sensitive to BMP4 (high therapeutic efficacy) and one with low-sensitivity, response to treatment with BMP4 was characterized. We applied single-cell RNA-sequencing, analyzed the relative abundance of cell clusters, searched for and identified the aforementioned two marker types, and validated these results in all 17 cultures. RESULTS: High variation in cell viability was observed after treatment with BMP4. In three cultures with highest sensitivity for BMP4, a substantial new cell subpopulation formed. These cells displayed decreased cell proliferation and increased apoptosis. Neuronal differentiation was reduced most in cultures with little sensitivity for BMP4. OLIG1/2 levels were found predictive for high sensitivity to BMP4. Activation of ribosomal translation (RPL27A, RPS27) was up-regulated within one day in cultures that were very sensitive to BMP4. CONCLUSION: The changes in composition of patient-derived GBM cultures obtained after treatment with BMP4 correlate with treatment efficacy. OLIG1/2 expression can predict this efficacy, and upregulation of RPL27A and RPS27 are useful early-response markers.

Cross-species analysis of enhancer logic using deep learning.

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Robust gene expression programs underlie recurrent cell states and phenotype switching in melanoma.

Melanoma cells can switch between a melanocytic and a mesenchymal-like state. Scattered evidence indicates that additional intermediate state(s) may exist. Here, to search for such states and decipher their underlying gene regulatory network (GRN), we studied 10 melanoma cultures using single-cell RNA sequencing (RNA-seq) as well as 26 additional cultures using bulk RNA-seq. Although each culture exhibited a unique transcriptome, we identified shared GRNs that underlie the extreme melanocytic and mesenchymal states and the intermediate state. This intermediate state is corroborated by a distinct chromatin landscape and is governed by the transcription factors SOX6, NFATC2, EGR3, ELF1 and ETV4. Single-cell migration assays confirmed the intermediate migratory phenotype of this state. Using time-series sampling of single cells after knockdown of SOX10, we unravelled the sequential and recurrent arrangement of GRNs during phenotype switching. Taken together, these analyses indicate that an intermediate state exists and is driven by a distinct and stable 'mixed' GRN rather than being a symbiotic heterogeneous mix of cells.

PlcRa, a new quorum-sensing regulator from Bacillus cereus, plays a role in oxidative stress responses and cysteine metabolism in stationary phase.

We characterized a new quorum-sensing regulator, PlcRa, which is present in various members of the B. cereus group and identified a signaling heptapeptide for PlcRa activity: PapRa(7). We demonstrated that PlcRa is a 3D structural paralog of PlcR using sequence analysis and homology modeling. A comparison of the transcriptomes at the onset of stationary phase of a ΔplcRa mutant and the wild-type B. cereus ATCC 14579 strain showed that 68 genes were upregulated and 49 genes were downregulated in the ΔplcRa mutant strain (>3-fold change). Genes involved in the cysteine metabolism (putative CymR regulon) were downregulated in the ΔplcRa mutant strain. We focused on the gene with the largest difference in expression level between the two conditions, which encoded -AbrB2- a new regulator of the AbrB family. We demonstrated that purified PlcRa bound specifically to the abrB2 promoter in the presence of synthetic PapRa(7), in an electrophoretic mobility shift assay. We further showed that the AbrB2 regulator controlled the expression of the yrrT operon involved in methionine to cysteine conversion. We found that the ΔplcRa mutant strain was more sensitive to hydrogen peroxide- and disulfide-induced stresses than the wild type. When cystine was added to the culture of the ΔplcRa mutant, challenged with hydrogen peroxide, growth inhibition was abolished. In conclusion, we identified a new RNPP transcriptional regulator in B. cereus that activated the oxidative stress response and cysteine metabolism in transition state cells.

Efficient gene expression profiling of laser-microdissected melanoma metastases.

Comparing the transcriptomes of primary and metastatic tumour tissues is a useful strategy for studying tumour progression. One factor limiting the interpretation of tissue-based transcriptomic data is the lack of cell-type purity. Laser capture microdissection (LCM) has been shown to be useful for overcoming this limitation. We established an efficient protocol for gene expression profiling of LCM and matched metastatic melanomas using a transgenic mouse model. This optimized workflow combines microsurgical recovery of mouse lungs, appropriate tissue freezing, laser microdissection of homogeneous tumour cell populations from cryosections, isolation of high-quality RNA and gene expression analysis. The RNA isolated from laser-microdissected material was not contaminated by stroma cells, was of excellent quality, and the synthesis of cDNAs was homogeneous and highly reproducible. Subsequent custom-based Taqman-low-density-array (TLDA)-based gene expression profiling identified stronger expression of five genes (M-MITF, TYR, STAT3, CCND1 and PAX3) in primary than metastatic melanoma. We detected only minor transcriptomic differences between primary and metastatic melanoma tissue. This optimized workflow could be very valuable for various studies requiring cell type-specific transcriptomic analysis.

Maintaining RNA integrity in a homogeneous population of mammary epithelial cells isolated by Laser Capture Microdissection.

BACKGROUND: Laser-capture microdissection (LCM) that enables the isolation of specific cell populations from complex tissues under morphological control is increasingly used for subsequent gene expression studies in cell biology by methods such as real-time quantitative PCR (qPCR), microarrays and most recently by RNA-sequencing. Challenges are i) to select precisely and efficiently cells of interest and ii) to maintain RNA integrity. The mammary gland which is a complex and heterogeneous tissue, consists of multiple cell types, changing in relative proportion during its development and thus hampering gene expression profiling comparison on whole tissue between physiological stages. During lactation, mammary epithelial cells (MEC) are predominant. However several other cell types, including myoepithelial (MMC) and immune cells are present, making it difficult to precisely determine the specificity of gene expression to the cell type of origin. In this work, an optimized reliable procedure for producing RNA from alveolar epithelial cells isolated from frozen histological sections of lactating goat, sheep and cow mammary glands using an infrared-laser based Arcturus Veritas LCM (Applied Biosystems®) system has been developed. The following steps of the microdissection workflow: cryosectioning, staining, dehydration and harvesting of microdissected cells have been carefully considered and designed to ensure cell capture efficiency without compromising RNA integrity. RESULTS: The best results were obtained when staining 8 µm-thick sections with Cresyl violet® (Ambion, Applied Biosystems®) and capturing microdissected cells during less than 2 hours before RNA extraction. In addition, particular attention was paid to animal preparation before biopsies or slaughtering (milking) and freezing of tissue blocks which were embedded in a cryoprotective compound before being immersed in isopentane. The amount of RNA thus obtained from ca.150 to 250 acini (300,000 to 600,000 µm2) ranges between 5 to 10 ng. RNA integrity number (RIN) was ca. 8.0 and selectivity of this LCM protocol was demonstrated through qPCR analyses for several alveolar cell specific genes, including LALBA (α-lactalbumin) and CSN1S2 (αs2-casein), as well as Krt14 (cytokeratin 14), CD3e and CD68 which are specific markers of MMC, lymphocytes and macrophages, respectively. CONCLUSIONS: RNAs isolated from MEC in this manner were of very good quality for subsequent linear amplification, thus making it possible to establish a referential gene expression profile of the healthy MEC, a useful platform for tumor biomarker discovery.

Enterococcal leucine-rich repeat-containing protein involved in virulence and host inflammatory response.

Enterococcus faecalis is an important nosocomial pathogen associated with high morbidity and mortality for patients who are immunocompromised or who have severe underlying diseases. The E. faecalis genome encodes numerous surface-exposed proteins that may be involved in virulence. This work describes the characterization of the first internalin-like protein in E. faecalis, ElrA, belonging to the recently identified WxL family of surface proteins. ElrA contains an N-terminal signal peptide for export, a leucine-rich repeat domain that may interact with host cells, and a C-terminal WxL domain that interacts with the peptidoglycan. Disruption of the elrA gene significantly attenuates bacterial virulence in a mouse peritonitis model. The elrA deletion mutant also displays a defect in infection of host macrophages and a decreased interleukin-6 response in vivo. Finally, elrA expression is induced in vivo. Altogether, these results demonstrate a role for ElrA in the E. faecalis infectious process in vivo and suggest that this surface protein may contribute to E. faecalis virulence by stimulating the host inflammatory response.