Elucidating the Transcriptional States of Spermatogenesis-Joint Analysis of Germline and Supporting Cell, Mice and Human, Normal and Perturbed, Bulk and Single-Cell RNA-Seq.

In studying the molecular underpinning of spermatogenesis, we expect to understand the fundamental biological processes better and potentially identify genes that may lead to novel diagnostic and therapeutic strategies toward precision medicine in male infertility. In this review, we emphasized our perspective that the path forward necessitates integrative studies that rely on complementary approaches and types of data. To comprehensively analyze spermatogenesis, this review proposes four axes of integration. First, spanning the analysis of spermatogenesis in the healthy state alongside pathologies. Second, the experimental analysis of model systems (in which we can deploy treatments and perturbations) alongside human data. Third, the phenotype is measured alongside its underlying molecular profiles using known markers augmented with unbiased profiles. Finally, the testicular cells are studied as ecosystems, analyzing the germ cells alongside the states observed in the supporting somatic cells. Recently, the study of spermatogenesis has been advancing using single-cell RNA sequencing, where scientists have uncovered the unique stages of germ cell development in mice, revealing new regulators of spermatogenesis and previously unknown cell subtypes in the testis. An in-depth analysis of meiotic and postmeiotic stages led to the discovery of marker genes for spermatogonia, Sertoli and Leydig cells and further elucidated all the other germline and somatic cells in the testis microenvironment in normal and pathogenic conditions. The outcome of an integrative analysis of spermatogenesis using advanced molecular profiling technologies such as scRNA-seq has already propelled our biological understanding, with additional studies expected to have clinical implications for the study of male fertility. By uncovering new genes and pathways involved in abnormal spermatogenesis, we may gain insights into subfertility or sterility.

PD1 ligand functionality a biomarker of response to anti PD1 treatment in patients with HNSCC.

Therapies targeting the PD-1/PD-L1 pathway have transformed head and neck squamous cell carcinoma (HNSCC) treatment. However, predicting the response to anti-PD-1 therapy remains a clinical challenge. This study evaluated the functional binding of PD-1 ligands in 29 HNSCC patients and compared it to the standard PD-L1 Combined Positive Score (CPS). The assessment of PD-1 ligands' functionality advances the current ability to predict the response of HNSCC patients to anti-PD-1 therapy.

High FOXA1 levels induce ER transcriptional reprogramming, a pro-metastatic secretome, and metastasis in endocrine-resistant breast cancer.

Aberrant activation of the forkhead protein FOXA1 is observed in advanced hormone-related cancers. However, the key mediators of high FOXA1 signaling remain elusive. We demonstrate that ectopic high FOXA1 (H-FOXA1) expression promotes estrogen receptor-positive (ER+) breast cancer (BC) metastasis in a xenograft mouse model. Mechanistically, H-FOXA1 reprograms ER-chromatin binding to elicit a core gene signature (CGS) enriched in ER+ endocrine-resistant (EndoR) cells. We identify Secretome14, a CGS subset encoding ER-dependent cancer secretory proteins, as a strong predictor for poor outcomes of ER+ BC. It is elevated in ER+ metastases vs. primary tumors, irrespective of ESR1 mutations. Genomic ER binding near Secretome14 genes is also increased in mutant ER-expressing or mitogen-treated ER+ BC cells and in ER+ metastatic vs. primary tumors, suggesting a convergent pathway including high growth factor receptor signaling in activating pro-metastatic secretome genes. Our findings uncover H-FOXA1-induced ER reprogramming that drives EndoR and metastasis partly via an H-FOXA1/ER-dependent secretome.

Functional binding of PD1 ligands predicts response to anti-PD1 treatment in patients with cancer.

Accurate predictive biomarkers of response to immune checkpoint inhibitors (ICIs) are required for better stratifying patients with cancer to ICI treatments. Here, we present a new concept for a bioassay to predict the response to anti-PD1 therapies, which is based on measuring the binding functionality of PDL1 and PDL2 to their receptor, PD1. In detail, we developed a cell-based reporting system, called the immuno-checkpoint artificial reporter with overexpression of PD1 (IcAR-PD1) and evaluated the functionality of PDL1 and PDL2 binding in tumor cell lines, patient-derived xenografts, and fixed-tissue tumor samples obtained from patients with cancer. In a retrospective clinical study, we found that the functionality of PDL1 and PDL2 predicts response to anti-PD1 and that the functionality of PDL1 binding is a more effective predictor than PDL1 protein expression alone. Our findings suggest that assessing the functionality of ligand binding is superior to staining of protein expression for predicting response to ICIs.

Mutated HRAS Activates YAP1-AXL Signaling to Drive Metastasis of Head and Neck Cancer.

The survival rate for patients with head and neck cancer (HNC) diagnosed with cervical lymph node (cLN) or distant metastasis is low. Genomic alterations in the HRAS oncogene are associated with advanced tumor stage and metastasis in HNC. Elucidation of the molecular mechanisms by which mutated HRAS (HRASmut) facilitates HNC metastasis could lead to improved treatment options for patients. Here, we examined metastasis driven by mutant HRAS in vitro and in vivo using HRASmut human HNC cell lines, patient-derived xenografts, and a novel HRASmut syngeneic model. Genetic and pharmacological manipulations indicated that HRASmut was sufficient to drive invasion in vitro and metastasis in vivo. Targeted proteomic analysis showed that HRASmut promoted AXL expression via suppressing the Hippo pathway and stabilizing YAP1 activity. Pharmacological blockade of HRAS signaling with the farnesyltransferase inhibitor tipifarnib activated the Hippo pathway and reduced the nuclear export of YAP1, thus suppressing YAP1-mediated AXL expression and metastasis. AXL was required for HRASmut cells to migrate and invade in vitro and to form regional cLN and lung metastases in vivo. In addition, AXL-depleted HRASmut tumors displayed reduced lymphatic and vascular angiogenesis in the primary tumor. Tipifarnib treatment also regulated AXL expression and attenuated VEGFA and VEGFC expression, thus regulating tumor-induced vascular formation and metastasis. Our results indicate that YAP1 and AXL are crucial factors for HRASmut-induced metastasis and that tipifarnib treatment can limit the metastasis of HNC tumors with HRAS mutations by enhancing YAP1 cytoplasmic sequestration and downregulating AXL expression. SIGNIFICANCE: Mutant HRAS drives metastasis of head and neck cancer by switching off the Hippo pathway to activate the YAP1-AXL axis and to stimulate lymphovascular angiogenesis.

Multidimensional Molecular Profiling of Metastatic Triple-Negative Breast Cancer and Immune Checkpoint Inhibitor Benefit.

PURPOSE: In metastatic triple-negative breast cancer (mTNBC), consistent biomarkers of immune checkpoint inhibitor (ICI) therapy benefit remain elusive. We evaluated the immune, genomic, and transcriptomic landscape of mTNBC in patients treated with ICIs. METHODS: We identified 29 patients with mTNBC treated with pembrolizumab or atezolizumab, either alone (n = 9) or in combination with chemotherapy (n = 14) or targeted therapy (n = 6), who had tumor tissue and/or blood available before ICI therapy for whole-exome sequencing. RNA sequencing and CIBERSORTx-inferred immune population analyses were performed (n = 20). Immune cell populations and programmed death-ligand 1 expression were assessed using multiplexed immunofluorescence (n = 18). Clonal trajectories were evaluated via serial tumor/circulating tumor DNA whole-exome sequencing (n = 4). Association of biomarkers with progression-free survival and overall survival (OS) was assessed. RESULTS: Progression-free survival and OS were longer in patients with high programmed death-ligand 1 expression and tumor mutational burden. Patients with longer survival also had a higher relative inferred fraction of CD8+ T cells, activated CD4+ memory T cells, M1 macrophages, and follicular helper T cells and enrichment of inflammatory gene expression pathways. A mutational signature of defective repair of DNA damage by homologous recombination was enriched in patients with both shorter OS and primary resistance. Exploratory analysis of clonal evolution among four patients treated with programmed cell death protein 1 blockade and a tyrosine kinase inhibitor suggested that clonal stability post-treatment was associated with short time to progression. CONCLUSION: This study identified potential biomarkers of response to ICIs among patients with mTNBC: high tumor mutational burden; presence of CD8+, CD4 memory T cells, follicular helper T cells, and M1 macrophages; and inflammatory gene expression pathways. Pretreatment deficiencies in the homologous recombination DNA damage repair pathway and the absence of or minimal clonal evolution post-treatment may be associated with worse outcomes.

ESR1 mutant breast cancers show elevated basal cytokeratins and immune activation.

Estrogen receptor alpha (ER/ESR1) is frequently mutated in endocrine resistant ER-positive (ER+) breast cancer and linked to ligand-independent growth and metastasis. Despite the distinct clinical features of ESR1 mutations, their role in intrinsic subtype switching remains largely unknown. Here we find that ESR1 mutant cells and clinical samples show a significant enrichment of basal subtype markers, and six basal cytokeratins (BCKs) are the most enriched genes. Induction of BCKs is independent of ER binding and instead associated with chromatin reprogramming centered around a progesterone receptor-orchestrated insulated neighborhood. BCK-high ER+ primary breast tumors exhibit a number of enriched immune pathways, shared with ESR1 mutant tumors. S100A8 and S100A9 are among the most induced immune mediators and involve in tumor-stroma paracrine crosstalk inferred by single-cell RNA-seq from metastatic tumors. Collectively, these observations demonstrate that ESR1 mutant tumors gain basal features associated with increased immune activation, encouraging additional studies of immune therapeutic vulnerabilities.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Opposing immune and genetic mechanisms shape oncogenic programs in synovial sarcoma.

Synovial sarcoma (SyS) is an aggressive neoplasm driven by the SS18-SSX fusion, and is characterized by low T cell infiltration. Here, we studied the cancer-immune interplay in SyS using an integrative approach that combines single-cell RNA sequencing (scRNA-seq), spatial profiling and genetic and pharmacological perturbations. scRNA-seq of 16,872 cells from 12 human SyS tumors uncovered a malignant subpopulation that marks immune-deprived niches in situ and is predictive of poor clinical outcomes in two independent cohorts. Functional analyses revealed that this malignant cell state is controlled by the SS18-SSX fusion, is repressed by cytokines secreted by macrophages and T cells, and can be synergistically targeted with a combination of HDAC and CDK4/CDK6 inhibitors. This drug combination enhanced malignant-cell immunogenicity in SyS models, leading to induced T cell reactivity and T cell-mediated killing. Our study provides a blueprint for investigating heterogeneity in fusion-driven malignancies and demonstrates an interplay between immune evasion and oncogenic processes that can be co-targeted in SyS and potentially in other malignancies.

Acquired FGFR and FGF Alterations Confer Resistance to Estrogen Receptor (ER) Targeted Therapy in ER+ Metastatic Breast Cancer.

PURPOSE: To identify clinically relevant mechanisms of resistance to ER-directed therapies in ER+ breast cancer. EXPERIMENTAL DESIGN: We conducted a genome-scale functional screen spanning 10,135 genes to investigate genes whose overexpression confer resistance to selective estrogen receptor degraders. In parallel, we performed whole-exome sequencing in paired pretreatment and postresistance biopsies from 60 patients with ER+ metastatic breast cancer who had developed resistance to ER-targeted therapy. Furthermore, we performed experiments to validate resistance genes/pathways and to identify drug combinations to overcome resistance. RESULTS: Pathway analysis of candidate resistance genes demonstrated that the FGFR, ERBB, insulin receptor, and MAPK pathways represented key modalities of resistance. The FGFR pathway was altered via FGFR1, FGFR2, or FGF3 amplifications or FGFR2 mutations in 24 (40%) of the postresistance biopsies. In 12 of the 24 postresistance tumors exhibiting FGFR/FGF alterations, these alterations were acquired or enriched under the selective pressure of ER-directed therapy. In vitro experiments in ER+ breast cancer cells confirmed that FGFR/FGF alterations led to fulvestrant resistance as well as cross-resistance to the CDK4/6 inhibitor palbociclib. RNA sequencing of resistant cell lines demonstrated that FGFR/FGF induced resistance through ER reprogramming and activation of the MAPK pathway. The resistance phenotypes were reversed by FGFR inhibitors, a MEK inhibitor, and/or a SHP2 inhibitor. CONCLUSIONS: Our results suggest that FGFR pathway is a distinct mechanism of acquired resistance to ER-directed therapy that can be overcome by FGFR and/or MAPK pathway inhibitors.

Author Correction: A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The Genomic Landscape of Intrinsic and Acquired Resistance to Cyclin-Dependent Kinase 4/6 Inhibitors in Patients with Hormone Receptor-Positive Metastatic Breast Cancer.

Mechanisms driving resistance to cyclin-dependent kinase 4/6 inhibitors (CDK4/6i) in hormone receptor-positive (HR+) breast cancer have not been clearly defined. Whole-exome sequencing of 59 tumors with CDK4/6i exposure revealed multiple candidate resistance mechanisms including RB1 loss, activating alterations in AKT1, RAS, AURKA, CCNE2, ERBB2, and FGFR2, and loss of estrogen receptor expression. In vitro experiments confirmed that these alterations conferred CDK4/6i resistance. Cancer cells cultured to resistance with CDK4/6i also acquired RB1, KRAS, AURKA, or CCNE2 alterations, which conferred sensitivity to AURKA, ERK, or CHEK1 inhibition. Three of these activating alterations-in AKT1, RAS, and AURKA-have not, to our knowledge, been previously demonstrated as mechanisms of resistance to CDK4/6i in breast cancer preclinically or in patient samples. Together, these eight mechanisms were present in 66% of resistant tumors profiled and may define therapeutic opportunities in patients. SIGNIFICANCE: We identified eight distinct mechanisms of resistance to CDK4/6i present in 66% of resistant tumors profiled. Most of these have a therapeutic strategy to overcome or prevent resistance in these tumors. Taken together, these findings have critical implications related to the potential utility of precision-based approaches to overcome resistance in many patients with HR+ metastatic breast cancer.This article is highlighted in the In This Issue feature, p. 1079.

A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors.

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

Sensitive Detection of Minimal Residual Disease in Patients Treated for Early-Stage Breast Cancer.

PURPOSE: Existing cell-free DNA (cfDNA) methods lack the sensitivity needed for detecting minimal residual disease (MRD) following therapy. We developed a test for tracking hundreds of patient-specific mutations to detect MRD with a 1,000-fold lower error rate than conventional sequencing. EXPERIMENTAL DESIGN: We compared the sensitivity of our approach to digital droplet PCR (ddPCR) in a dilution series, then retrospectively identified two cohorts of patients who had undergone prospective plasma sampling and clinical data collection: 16 patients with ER+/HER2- metastatic breast cancer (MBC) sampled within 6 months following metastatic diagnosis and 142 patients with stage 0 to III breast cancer who received curative-intent treatment with most sampled at surgery and 1 year postoperative. We performed whole-exome sequencing of tumors and designed individualized MRD tests, which we applied to serial cfDNA samples. RESULTS: Our approach was 100-fold more sensitive than ddPCR when tracking 488 mutations, but most patients had fewer identifiable tumor mutations to track in cfDNA (median = 57; range = 2-346). Clinical sensitivity was 81% (n = 13/16) in newly diagnosed MBC, 23% (n = 7/30) at postoperative and 19% (n = 6/32) at 1 year in early-stage disease, and highest in patients with the most tumor mutations available to track. MRD detection at 1 year was strongly associated with distant recurrence [HR = 20.8; 95% confidence interval, 7.3-58.9]. Median lead time from first positive sample to recurrence was 18.9 months (range = 3.4-39.2 months). CONCLUSIONS: Tracking large numbers of individualized tumor mutations in cfDNA can improve MRD detection, but its sensitivity is driven by the number of tumor mutations available to track.

FOXA1 upregulation promotes enhancer and transcriptional reprogramming in endocrine-resistant breast cancer.

Forkhead box A1 (FOXA1) is a pioneer factor that facilitates chromatin binding and function of lineage-specific and oncogenic transcription factors. Hyperactive FOXA1 signaling due to gene amplification or overexpression has been reported in estrogen receptor-positive (ER+) endocrine-resistant metastatic breast cancer. However, the molecular mechanisms by which FOXA1 up-regulation promotes these processes and the key downstream targets of the FOXA1 oncogenic network remain elusive. Here, we demonstrate that FOXA1 overexpression in ER+ breast cancer cells drives genome-wide enhancer reprogramming to activate prometastatic transcriptional programs. Up-regulated FOXA1 employs superenhancers (SEs) to synchronize transcriptional reprogramming in endocrine-resistant breast cancer cells, reflecting an early embryonic development process. We identify the hypoxia-inducible transcription factor hypoxia-inducible factor-2α (HIF-2α) as the top high FOXA1-induced SE target, mediating the impact of high FOXA1 in activating prometastatic gene sets and pathways associated with poor clinical outcome. Using clinical ER+/HER2- metastatic breast cancer datasets, we show that the aberrant FOXA1/HIF-2α transcriptional axis is largely nonconcurrent with the ESR1 mutations, suggesting different mechanisms of endocrine resistance and treatment strategies. We further demonstrate the selective efficacy of an HIF-2α antagonist, currently in clinical trials for advanced kidney cancer and recurrent glioblastoma, in reducing the clonogenicity, migration, and invasion of endocrine-resistant breast cancer cells expressing high FOXA1. Our study has uncovered high FOXA1-induced enhancer reprogramming and HIF-2α-dependent transcriptional programs as vulnerable targets for treating endocrine-resistant and metastatic breast cancer.

The Combination of Whole Cell Lipidomics Analysis and Single Cell Confocal Imaging of Fluidity and Micropolarity Provides Insight into Stress-Induced Lipid Turnover in Subcellular Organelles of Pancreatic Beta Cells.

Modern omics techniques reveal molecular structures and cellular networks of tissues and cells in unprecedented detail. Recent advances in single cell analysis have further revolutionized all disciplines in cellular and molecular biology. These methods have also been employed in current investigations on the structure and function of insulin secreting beta cells under normal and pathological conditions that lead to an impaired glucose tolerance and type 2 diabetes. Proteomic and transcriptomic analyses have pointed to significant alterations in protein expression and function in beta cells exposed to diabetes like conditions (e.g., high glucose and/or saturated fatty acids levels). These nutritional overload stressful conditions are often defined as glucolipotoxic due to the progressive damage they cause to the cells. Our recent studies on the rat insulinoma-derived INS-1E beta cell line point to differential effects of such conditions in the phospholipid bilayers in beta cells. This review focuses on confocal microscopy-based detection of these profound alterations in the plasma membrane and membranes of insulin granules and lipid droplets in single beta cells under such nutritional load conditions.

Pooled Genomic Screens Identify Anti-apoptotic Genes as Targetable Mediators of Chemotherapy Resistance in Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) is often sensitive to initial treatment with platinum and taxane combination chemotherapy, but most patients relapse with chemotherapy-resistant disease. To systematically identify genes modulating chemotherapy response, we performed pooled functional genomic screens in HGSOC cell lines treated with cisplatin, paclitaxel, or cisplatin plus paclitaxel. Genes in the intrinsic pathway of apoptosis were among the top candidate resistance genes in both gain-of-function and loss-of-function screens. In an open reading frame overexpression screen, followed by a mini-pool secondary screen, anti-apoptotic genes including BCL2L1 (BCL-XL) and BCL2L2 (BCL-W) were associated with chemotherapy resistance. In a CRISPR-Cas9 knockout screen, loss of BCL2L1 decreased cell survival whereas loss of proapoptotic genes promoted resistance. To dissect the role of individual anti-apoptotic proteins in HGSOC chemotherapy response, we evaluated overexpression or inhibition of BCL-2, BCL-XL, BCL-W, and MCL1 in HGSOC cell lines. Overexpression of anti-apoptotic proteins decreased apoptosis and modestly increased cell viability upon cisplatin or paclitaxel treatment. Conversely, specific inhibitors of BCL-XL, MCL1, or BCL-XL/BCL-2, but not BCL-2 alone, enhanced cell death when combined with cisplatin or paclitaxel. Anti-apoptotic protein inhibitors also sensitized HGSOC cells to the poly (ADP-ribose) polymerase inhibitor olaparib. These unbiased screens highlight anti-apoptotic proteins as mediators of chemotherapy resistance in HGSOC, and support inhibition of BCL-XL and MCL1, alone or combined with chemotherapy or targeted agents, in treatment of primary and recurrent HGSOC. IMPLICATIONS: Anti-apoptotic proteins modulate drug resistance in ovarian cancer, and inhibitors of BCL-XL or MCL1 promote cell death in combination with chemotherapy.

Subclonal cooperation drives metastasis by modulating local and systemic immune microenvironments.

Most human tumours are heterogeneous, composed of cellular clones with different properties present at variable frequencies. Highly heterogeneous tumours have poor clinical outcomes, yet the underlying mechanism remains poorly understood. Here, we show that minor subclones of breast cancer cells expressing IL11 and FIGF (VEGFD) cooperate to promote metastatic progression and generate polyclonal metastases composed of driver and neutral subclones. Expression profiling of the epithelial and stromal compartments of monoclonal and polyclonal primary and metastatic lesions revealed that this cooperation is indirect, mediated through the local and systemic microenvironments. We identified neutrophils as a leukocyte population stimulated by the IL11-expressing minor subclone and showed that the depletion of neutrophils prevents metastatic outgrowth. Single-cell RNA-seq of CD45+ cell populations from primary tumours, blood and lungs demonstrated that IL11 acts on bone-marrow-derived mesenchymal stromal cells, which induce pro-tumorigenic and pro-metastatic neutrophils. Our results indicate key roles for non-cell-autonomous drivers and minor subclones in metastasis.

A Functional Landscape of Resistance to MEK1/2 and CDK4/6 Inhibition in NRAS-Mutant Melanoma.

Combinatorial inhibition of MEK1/2 and CDK4/6 is currently undergoing clinical investigation in NRAS-mutant melanoma. To prospectively map the landscape of resistance to this investigational regimen, we utilized a series of gain- and loss-of-function forward genetic screens to identify modulators of resistance to clinical inhibitors of MEK1/2 and CDK4/6 alone and in combination. First, we identified NRAS-mutant melanoma cell lines that were dependent on NRAS for proliferation and sensitive to MEK1/2 and CDK4/6 combination treatment. We then used a genome-scale ORF overexpression screen and a CRISPR knockout screen to identify modulators of resistance to each inhibitor alone or in combination. These orthogonal screening approaches revealed concordant means of achieving resistance to this therapeutic modality, including tyrosine kinases, RAF, RAS, AKT, and PI3K signaling. Activated KRAS was sufficient to cause resistance to combined MEK/CDK inhibition and to replace genetic depletion of oncogenic NRAS. In summary, our comprehensive functional genetic screening approach revealed modulation of resistance to the inhibition of MEK1/2, CDK4/6, or their combination in NRAS-mutant melanoma. SIGNIFICANCE: These findings reveal that NRAS-mutant melanomas can acquire resistance to genetic ablation of NRAS or combination MEK1/2 and CDK4/6 inhibition by upregulating activity of the RTK-RAS-RAF and RTK-PI3K-AKT signaling cascade.