BRAFV600E induces reversible mitotic arrest in human melanocytes via microRNA-mediated suppression of AURKB
Abstract
Benign melanocytic nevi frequently emerge when an acquired BRAFV600E mutation triggers unchecked proliferation and subsequent arrest in melanocytes. Recent observations have challenged the role of oncogene-induced senescence in melanocytic nevus formation, necessitating investigations into alternative mechanisms for the establishment and maintenance of proliferation arrest in nevi. We compared the transcriptomes of melanocytes from healthy human skin, nevi, and melanomas arising from nevi and identified a set of microRNAs as highly expressed nevus-enriched transcripts. Two of these microRNAs—MIR211-5p and MIR328-3p—induced mitotic failure, genome duplication, and proliferation arrest in human melanocytes through convergent targeting of AURKB. We demonstrate that BRAFV600E induces a similar proliferation arrest in primary human melanocytes that is both reversible and conditional. Specifically, BRAFV600E expression stimulates either arrest or proliferation depending on the differentiation state of the melanocyte. We report genome duplication in human melanocytic nevi, reciprocal expression of AURKB and microRNAs in nevi and melanomas, and rescue of arrested human nevus cells with AURKB expression. Taken together, our data describe an alternative molecular mechanism for melanocytic nevus formation that is congruent with both experimental and clinical observations.
eLife digest
Lots of people have small dark patches on their skin known as moles. Most moles form when individual cells known as melanocytes in the skin acquire a specific genetic mutation in a gene called BRAF. This mutation causes the cells to divide rapidly to form the mole. After a while, most moles stop growing and remain harmless for the rest of a person’s life.
Melanoma is a type of skin cancer that develops from damaged melanocytes. The same mutation in BRAF that is found in moles is also present in half of all cases of melanoma. Unlike in moles, the melanoma-causing mutation makes the melanocytes divide rapidly to form a tumor that keeps on growing indefinitely. It remains unclear why the same genetic mutation in the BRAF gene has such different consequences in moles and melanomas.
To address this question, McNeal et al. used genetic approaches to study melanocytes from moles and melanomas. The experiments identified some molecules known as microRNAs that are present at higher levels in moles than in melanomas. Increasing the levels of two of these microRNAs in melanocytes from human skin stopped the cells from growing and dividing by inhibiting a gene called AURKB. This suggested that these microRNAs are responsible for halting the growth of moles.
Introducing the mutated form of BRAF into melanocytes also stopped cells from growing and dividing by inhibiting AURKB. However, changing the environment surrounding the cells reversed this effect and allowed the melanocytes to resume dividing. In this way the mutated form of BRAF acts like a switch that allows melanocytes in skin cancers to start growing again under certain conditions.
Further experiments found that a drug called barasertib is able to inhibit the growth of melanoma cells with the mutant form of BRAF. Future work will investigate whether it is possible to use this drug and other tools to stop skin cancer tumors from growing, and possibly even prevent skin tumors from forming in the first place.
Introduction
Cutaneous melanoma is a potentially deadly skin cancer arising from the pigment-producing melanocytes of the human epidermis. An activating mutation in the BRAF proto-oncogene (BRAFV600E) drives over half of all cutaneous melanomas (Garnett and Marais, 2004; Davies et al., 2002; Pollock et al., 2003). Yet, when a melanocyte acquires a BRAFV600E mutation, the cell does not immediately transform to cancer. Instead, it usually undergoes clonal proliferation followed by stable arrest resulting in a benign skin lesion known as a melanocytic nevus or ‘mole’ (Pollock et al., 2003; Shain et al., 2015; Bastian, 2014; Dankort et al., 2009). Despite the continued expression of BRAFV600E, the majority of melanocytic nevi remain innocuous for the lifespan of the individual, suggesting that nevus cells have robust intrinsic defenses against hyperproliferation. By inducing hyperproliferation and subsequent arrest, the BRAFV600E mutation elicits divergent, biphasic phenotypes within a single cell—a poorly understood phenomenon. Characterization of the mechanisms and environmental factors that distinguish the two phenotypes could illuminate candidate strategies for chemoprevention of melanocytic nevus formation or interception of melanoma initiation.
The BRAFV600E mutation is clonally present in ~80% of human melanocytic nevi (Pollock et al., 2003) and melanocyte-specific expression of the oncogene is sufficient to induce nevus formation in mice (Dankort et al., 2009; Damsky et al., 2015; Ruiz-Vega et al., 2020), suggesting that BRAFV600E drives proliferation arrest in melanocytes. The prevailing theory to explain this seemingly paradoxical role for BRAFV600E is oncogene-induced senescence (OIS) (Michaloglou et al., 2005). Cellular senescence is defined as the permanent cell-cycle arrest of previously replication-competent cells (He and Sharpless, 2017; Roy et al., 2020). Senescence is associated with a variety of molecular hallmarks including elevated expression of p16INK4A and other cyclin-dependent-kinase inhibitors, TP53, H2AX, lysosomal beta-galactosidase and senescence-associated secretory proteins (SASPs) as well as a DNA content profile indicative of arrest in the G0/G1 phases of the cell cycle (Roy et al., 2020; Serrano et al., 1997; Campisi and d’Adda di Fagagna, 2007; Collado and Serrano, 2006; Afshari et al., 1993; Chatsirisupachai et al., 2019). The concept of OIS as a barrier to tumorigenesis derived from observations of a durable proliferation arrest induced by overexpression of oncogenic RAS or RAF in immortalized human fibroblasts (Serrano et al., 1997; Zhu et al., 1998).
More recent studies have questioned whether the term ‘senescence’ aptly describes the proliferation arrest of melanocytic nevi. Cell cycle re-entry that accompanies nevus recurrence (King et al., 2009), eruption (Burian and Jemec, 2019), or transformation to primary melanoma suggests that the nevus arrest phenotype is reversible rather than permanent. Similarly, the expression of senescence markers does not readily distinguish human melanocytic nevi from primary or transformed melanocytes in humans (16INK4A, TP53, H2AX, and beta-galactosidase) or mice (any known hallmark of senescence) (Ruiz-Vega et al., 2020; Cotter et al., 2007; Tran et al., 2012). Setting aside the permanence of OIS, melanocytic nevus arrest was largely thought to be driven by induction of p16INK4A expression from the CDKN2A gene. Indeed, early studies demonstrated that BRAFV600E induces p16INK4A expression and that melanocytic nevi frequently express high levels of p16INK4A (Michaloglou et al., 2005). Despite this, p16INK4A expression is neither ubiquitous across melanocytic nevi, nor uniform within the cells of a single nevus (Oaxaca et al., 2020; Zeng et al., 2018), and neither knockdown of the 16INK4A transcript nor ablation of the CDKN2A locus prevents the onset or rescues BRAFV600E-induced proliferation arrest (Michaloglou et al., 2005; Zeng et al., 2018; Haferkamp et al., 2009; McNeal et al., 2015). BRAFV600E also stimulates expression of p15INK4B—the translational product of the CDKN2B gene (McNeal et al., 2015) which is situated adjacent to the CDKN2A locus—suggesting BRAFV600E might orchestrate arrest via multiple cyclin-dependent kinase inhibitors that induce arrest in the G1 phase of the cell cycle, prior to genome synthesis. While the loss of the CDKN2A/B locus is a formative event in melanoma progression, genetic evolution studies of melanoma progression suggest that selection against the CDKN2A/B locus occurs after BRAFV600E melanocytes have already escaped nevus-associated arrest (Zeng et al., 2018; Shain et al., 2018). Collectively, these data complicate the straightforward model wherein BRAFV600E induces the expression of cell cycle regulators leading to proliferation arrest, which is later circumvented by the genetic loss of those regulators.
Given these observations, we reasoned that the arrest response of human melanocytes to BRAFV600E might be distinct from previous reports of BRAFV600E-induced OIS in other cell types. Specifically, we hypothesized that nevus-associated arrest is conditional and coordinated by reversible changes in gene expression. Here, we identify a transcriptional program that is consistently elevated in melanocytic nevi and includes microRNAs (miRNAs) previously shown to have diagnostic value (Babapoor et al., 2016; Torres et al., 2020). In parallel, we observe that, unlike RAF-induced OIS in fibroblasts, BRAFV600E-induced arrest in human melanocytes is reversible and conditional. Merging these two observations, we discover that in melanocytes, regulation of miRNAs coupled to differentiation status permits toggling between BRAFV600E-driven hyperproliferative to arrested phenotypes.
Results
Melanocytic nevus-enriched miRNAs induce mitotic failure and proliferation arrest in primary human melanocytes
Since melanocytic nevi do not consistently express known senescence markers, we sought to characterize the transcripts that are specifically elevated within nevus melanocytes (Figure 1A). We analyzed previously published data sets comparing RNA expression between benign and malignant melanocytic lesions (Shain et al., 2018; Torres et al., 2020). Each clinical specimen contained matched RNA samples derived from melanocytic nevi and melanomas that arose from those nevi. We reasoned that downregulated transcripts within these newly formed melanomas could represent the most immediate barriers to melanocyte growth. Of the differentially expressed genes (adjusted p value<0.05), 6 of the top 10 with elevated expression in nevi were non-coding miRNAs (Figure 1B and Supplementary file 1). Among these were miRNAs previously shown as enriched in nevi compared to melanomas, including MIR211-5p, MIR125A-5p, MIR125B-5p, and MIR328-3p (Torres et al., 2020; Latchana et al., 2016). To determine whether elevated expression of these miRNAs distinguished nevus melanocytes from normal epidermal melanocytes, we profiled melanocytes from fresh healthy skin and nevi (Supplementary file 2). Three of the miRNAs (MIR211-5p, MIR328-3p, and MIR125B-5p) were also more highly expressed in melanocytic nevi as compared to healthy melanocytes or melanomas—an expression pattern consistent with enrichment in arrested melanocytes—whereas one miRNA (MIR125A-5p) presented a pattern consistent with lineage-specific expression that is lost upon transformation (Figure 1C). We conclude that MIR211-5p, MIR328-3p, and MIR125B-5p are transcripts specifically enriched in arrested melanocytic nevi.
To investigate whether elevated expression of MIR211-5p, MIR328-3p, or MIR125B-5p induce melanocytic arrest, we introduced synthesized mimetics of the mature miRNAs into primary human melanocytes and tracked proliferation over 7 days in culture. The introduction of either MIR211-5p or MIR328-3p resulted in significantly fewer melanocytes compared to a non-targeting control (Figure 1D). Increased concentrations of these mimics resulted in a dose-dependent decrease in EdU uptake, consistent with a cell cycle defect (Figure 1E). To better characterize the effect of MIR211-5p or MIR328-3p expression on cell cycle and cell viability, we performed time-lapsed quantitative phase imaging (QPI). We and others have previously established QPI as a method for adherent cell cytometry that can simultaneously measure melanocyte proliferation, cell death, and cell cycle stage with single-cell resolution (Hejna et al., 2017; Barker et al., 2020; Falck Miniotis et al., 2014). Images were taken every hour for 5 days post-nucleofection (Figure 1F). Cells expressing MIR211-5p or MIR328-3p exhibited significantly reduced growth curves (Figure 1G). We observed no evidence of increased cell death (Figure 1H), and confirmed this result with Annexin V staining (Figure 1I). Experiments performed using three additional melanocyte preparations counted over 7 days post nucleofection yielded similar results (Figure 1J). Thus, MIR211-5p and MIR328-3p induce proliferation arrest in human melanocytes.
Further characterization of the QPI images revealed a distinct increase in intra-cellular dry mass in proliferation-arrested conditions (Figure 1K). The dry mass of a cell is a measurement of total biomass (DNA, RNA, proteins, lipids, etc.) and is quantitatively measured by QPI. Relative dry mass fluctuates with cell cycle stage in a predictable manner (Falck Miniotis et al., 2014). The 20–40% increase we observed here is characteristic of cells that have doubled their DNA content but failed to complete mitosis. We confirmed this finding by measuring DNA content of the arrested populations (Figure 1L). We observed significant increases of cells with both 4 N DNA content and greater than 4 N DNA content, indicative of cytokinesis failure following DNA replication (Figure 1M). Consistent with this interpretation, time-lapsed imaging revealed incidences of tripolar mitosis and anaphase bridging in miRNA-expression melanocytes—both processes associated with cytokinesis failure (Normand and King, 2010)—which were not observed in control conditions or standard culture (Video 1 and Figure 1—figure supplement 1). We conclude that MIR211-5p and MIR328-3p induce mitotic failure and proliferation arrest when introduced into primary human melanocytes.
Inhibition of AURKB by MIR211-3p and MIR328-5p restricts proliferation in nevi
To investigate the mechanism by which MIR211-5p and MIR328-3p induce arrest in primary human melanocytes, we applied our previously established pipeline for identification of relevant miRNA targets (Figure 2A; Judson et al., 2013). First, we nucleofected MIR211-5p, MIR328-3p, or control RNA mimics into primary human melanocytes freshly isolated from four independent donors and conducted RNA sequencing (Supplementary file 3). We next selected all genes that were significantly downregulated by either miRNA as compared to the control RNA, and further refined the set to the computationally predicted targets of MIR211-5p and/or MIR328-3p. Finally, we nucleofected siRNAs targeting the 115 resulting genes and assayed for reduced EdU incorporation that would phenocopy MIR211-5p or MIR328-3p expression (Supplementary file 4). As expected, computationally predicted targets of MIR211-5p and MIR328-3p were enriched in the set of downregulated genes following MIR211-5p or MIR328-3p nucleofection, respectively (Figure 2B). AURKB, GPR3, and MAPKAPK2 showed the greatest magnitude of repression due to miRNA expression with significant EdU reduction upon siRNA-mediated knockdown (Figure 2C). Of these, AURKB (predicted MIR211-5p target) and GPR3 (predicted MIR328-3p target) were significantly inhibited upon nucleofection of either miRNA (Figure 2D). We reasoned these genes might represent convergent nodes driving proliferation arrest. Further exploration of the AURKB transcript revealed a potential MIR328-3p binding site just before the 3′ untranslated region (Figure 2E; Dweep et al., 2014). In contrast, we were unable to identify a candidate MIR211-5p binding site in the GPR3 transcript, suggesting indirect inhibition of this gene by MIR211-5p.
To determine whether decreased AURKB and/or GPR3 expression was associated with melanocytic nevi, we acquired 23 clinical specimens and measured both the transcript and protein expression of these genes (Figure 2F–H). The expression of both AURKB and GPR3 protein and mRNA was enriched in melanomas as compared to nevi (Figure 2G–H). Taken together, these data demonstrate an inverse expression between two miRNAs (MIR211-5p and MIR328-3p) and two targeted mRNAs (AURKB and GPR3) that toggles concurrently with the transition of arrested melanocytic nevi to melanomas. To further investigate whether either mRNA was sufficient to rescue miRNA-mediated arrest, primary melanocytes were transduced with constructs expressing either zsGreen (control), AURKB, or GPR3 (Figure 2I-K), then subsequently nucleofected with either MIR211-5p or MIR328-3p. Expression of either mRNA partially rescued MIR211-5p-induced arrest, while AURKB, but not GPR3, partially rescued MIR328-3p-induced arrest (Figure 2L). Thus, while other MIR211-5p and MIR328-3p targets, including GPR3, may play a role in the induction of melanocyte proliferation arrest, these data suggest that AURKB inhibition is a critical convergent node of both miRNAs in facilitating proliferation arrest in nevi. This interpretation is further supported by the well-established role of AURKB as an essential orchestrator of successful mitosis (D’Avino and Capalbo, 2015).
BRAFV600E induces a reversible and conditional proliferation arrest in human melanocytes
We hypothesized that AURKB inhibition by MIR211-5p and MIR328-3p may be the underlying mechanism by which BRAFV600E induces proliferation arrest in melanocytes. BRAFV600E expression has been demonstrated to result in an accumulation of growth-arrested fibroblasts with 2 N DNA content, presumably in a permanent G0/G1 senescent state (Michaloglou et al., 2005; Serrano et al., 1997; Zhu et al., 1998). In contrast, upon overexpression of MIR211-5p or MIR328-3p, we observed an accumulation of arrested melanocytes with ≥4 N DNA content, and a significant decrease in cells with 2 N DNA content (Figure 1M). To determine the DNA content profile of primary human melanocytes under BRAFV600E-induced arrest, we utilized an established lentiviral construct for dose-responsive doxycycline-inducible BRAFV600E expression (Figure 3A; McNeal et al., 2015). We monitored melanocyte growth under a titration of doxycycline and identified the lowest concentration (15.6 ng/ml) that induced sustained arrest after 40 hr exposure (Figure 3B). Similar to our observations for miRNA-induced arrest, relative per cell dry mass increased by 50% concurrent with BRAFV600E induced-arrest (Figure 3C). Cell cycle profiling confirmed a significant and dose-dependent accumulation of cells with ≥4 N DNA content and a decrease in cells with 2 N DNA content (Figure 3D–E). The cell cycle profiles were notably dissimilar to that of melanocytes treated with a pharmacologic CDK4/6 inhibitor (PD0332991, Palbociclib) (Clark et al., 2016) known to induce G0/G1 arrest. Treatment with a pharmacologic AURKB inhibitor (AZD2811, Barasertib) (Helfrich et al., 2016) that induces G2/M arrest and mitotic failure also resulted in an accumulation of melanocytes with ≥4 N DNA content and a decrease in cells with 2 N DNA content.
In addition to the accumulation of 2 N cells, previous reports of RAF expression in fibroblasts show that transient expression can induce a durable arrest phenotype (Zhu et al., 1998). To test the durability of BRAFV600E-induced arrest in human melanocytes, we induced doxycycline-dependent BRAFV600E expression for 48 hr and monitored cell number for an additional 2 days to confirm arrest. Media was then replaced to remove doxycycline. Within 12 hr of doxycycline removal, the previously arrested melanocytes began to proliferate (Figure 3F), demonstrating that BRAFV600E-induced proliferation arrest in melanocytes is reversible and dependent on continued oncogene expression.
Nevi recurrence after excision (King et al., 2009) and eruption upon external stimuli (Burian and Jemec, 2019), together with our observation that BRAFV600E induced arrest is reversible in melanocytes, suggested that BRAFV600E-induced arrest may depend on secondary signals from the microenvironment. Tetradecanoylphorbol acetate (TPA) is a protein kinase C (PKC) activator commonly used as a mitogen in primary melanocyte media to substitute for endothelin receptor type B activation (Hsu et al., 2005). In contrast, TPA also inhibits the cell cycle of melanoma cells (Coppock et al., 1995; Arita et al., 1998). BRAFV600E expression permits growth of the Melan-A cell line in the absence of TPA (Wellbrock et al., 2004). We therefore wanted to determine the relative contributions of TPA and BRAFV600E to the BRAFV600E-induced primary melanocyte arrest in vitro. We induced BRAFV600E expression in primary human melanocytes in two medias—one containing TPA (+TPA) or one containing the endothelin receptor type B ligand, endothelin-1 (TPA-free) and observed divergent dose-responsive phenotypes (Figure 3G). In the presence of TPA, BRAFV600E expression inhibited cell proliferation. By contrast, in TPA-free media, BRAFV600E expression induced proliferation, a phenotype that persisted for the remaining duration of the primary cultures (4–6 weeks). We next subjected melanocytes to the arrest-inducing conditions of doxycycline treatment in +TPA media for 4 days, then replaced with doxycycline-containing, TPA-free media. We observed a rapid shift to the proliferative phenotype even in the continued presence of doxycycline (Figure 3H). To determine whether the observed reversibility of BRAFV600E-induced growth arrest was due to insufficient duration or degree of expression of the oncogene, we exposed melanocytes to a high concentration of doxycycline (60 ng/ml) for 30 days prior to switching medias (Figure 3I). Melanocytes in +TPA media remained arrested for the duration of doxycycline treatment then grew again upon removal of either doxycycline or TPA (Figure 3J–K). From these results, we conclude that BRAFV600E expression induces divergent and interconvertible phenotypes—either hyperproliferation or a reversible arrest—in human melanocytes conditional on extrinsic stimuli.
BRAFV600E-induced arrest is dependent on melanocyte growth conditions and differentiation state
To investigate the relative contributions of BRAFV600E and TPA to the expression of MIR211-5p and MIR328-3p, we conducted small RNA sequencing of melanocytes grown in +TPA or TPA-free media and with or without BRAFV600E expression. The expression of MIR328-3p was relatively stable, demonstrating less than twofold changes across conditions (Figure 4A–B, Supplementary file 5). In contrast, MIR211-5p was one of the most differentially expressed genes in +TPA conditions compared to TPA-free conditions (adj. p=6.8E−23) (Figure 4A–B, Supplementary file 5). Induction of BRAFV600E in +TPA conditions did not further alter expression of either miRNA. Direct comparison of the arrested melanocytes in +TPA +Dox conditions to the proliferative melanocytes in −TPA +Dox conditions revealed MIR211-5p as the most significantly downregulated miRNA in the proliferative condition (Figure 4B). Thus, in comparison to TPA, BRAFV600E induction had minimal effect on the expression of the miRNAs. The most notable effect was a decrease of MIR328-3p expression coinciding with the induction of BRAFV600E in TPA-free conditions (Figure 4A). While MIR328-3p downregulation likely contributes to increased proliferation in this condition, the more pronounced effect of TPA on MIR211-5p expression suggests that the mechanisms governing whether BRAFV600E induces arrest or proliferation are a function of the cell’s transcriptional state rather than directly downstream of BRAFV600E itself.
To better characterize the two transcriptional states influencing the BRAFV600E phenotype, we conducted mRNA sequencing of the same specimens. Principal component (PC) analysis revealed that 97.4 % of the variance between samples was explained by two PCs—PC1 (81.75 variance) separated transcriptomes based upon exposure to TPA and PC2 (15.7% variance) separated transcriptomes based upon BRAFV600E expression (Figure 4C, Supplementary file 5). We conducted differential expression (DE) analyses to compare all samples exposed to +TPA versus TPA-free media, regardless of BRAFV600E expression (Figure 4D, top) and to compare all samples expressing BRAFV600E, regardless of TPA exposure (Figure 4D, bottom). As expected from the PC analysis and similar to the miRNA profiles, most gene expression variance was due to media conditions with a smaller, but consistent, set of genes activated by BRAFV600E.
One initially surprising result of this analysis was the static expression of both MIR211-5p and MIR328-3p in +TPA conditions, regardless of BRAFV600E expression. However, as inhibitors of translation, miRNA function depends on the expression level of its targets (Mukherji et al., 2011; Ovando-Vázquez et al., 2016), effectively functioning as buffers of gene expression that prevent aberrantly high levels of a transcript (Ebert and Sharp, 2012). BRAFV600E has been previously demonstrated as an upstream activator of AURKB transcription in melanoma cells (Sharma et al., 2013). We therefore reasoned that MIR211-5p expression in +TPA conditions might dampen AURKB activation by BRAFV600E. Consistent with this hypothesis, the addition of doxycycline induced AURKB expression in TPA-free/MIR211-5p low conditions, but this effect was not observed in +TPA/MIR211-5p high conditions (Figure 4E).
We next conducted gene set enrichment analyses (GSEA) comparing TPA-regulated genes and BRAFV600E- regulated genes to Molecular Signature Database Hallmark gene sets (Liberzon et al., 2015), the KEGG PATHWAY database (Kanehisa et al., 2016), and signatures important in BRAFV600E signaling (Ryu et al., 2011), human melanocyte differentiation (Belote et al., 2021), or melanoma cell state (Hoek and Goding, 2010; Tirosh et al., 2016; Tsoi et al., 2018). Among the pathways uniquely enriched with BRAFV600E expression were the expected increases in signatures downstream of BRAFV600E, MAPK (KRAS) signaling, and glycolysis (Haq et al., 2014 ; Figure 4F and Supplementary file 6). Gene sets uniquely enriched in +TPA media included signatures associated with cell cycle regulation, such as the G2-M checkpoint, oxidative phosphorylation, the MITF program, and differentiated melanocytes. In contrast, TPA-free media-induced gene sets are associated with cell adhesion and melanocyte progenitor and stem cells. Corroborating this observation, melanocytes grown in TPA-free media expressed higher levels of AXL (Figure 4G–H), lower levels of MLANA, a melanocyte differentiation marker, and presented greater activity of aldehyde dehydrogenase, enzymatic activity associated with stemness (Figure 4I). These results indicate that exposure to TPA induces a more differentiated melanocyte phenotype, consistent with previous reports (Prince et al., 2003; Vidács et al., 2021; Kormos et al., 2011).
To look for synergist effects of both BRAFV600E expression and media conditions, we performed GSEA analyses on each set of pair-wise conditions (Supplementary file 6 and summarized in Figure 4—figure supplement 1). Generally, most gene sets identified in Figure 4F were regulated by only one of the two stimuli. For example, genes downstream of BRAFV600E signaling were only influenced by doxycycline exposure. Interesting exceptions to this trend were the inflammatory response, which required both stimuli; cell cycle regulation gene sets, which were activated by TPA alone but further elevated when supplemented with doxycycline; and epithelial to mesenchymal transition genes, which were activated by either stimulus. A subset of melanocyte dedifferentiation gene sets, specifically neural crest genes, the melanoma invasion program, and the AXL program, but not markers of melanocyte stem cells, were activated by either TPA-free media or BRAFV600E expression, albeit the activation was more pronounced in the former. This observation suggests that BRAFV600E expression alone does induce some programs associated with melanocyte dedifferentiation, but with the retention of MIR211-5p expression, the induction is insufficient to overcome growth arrest.
Collectively, our transcriptomic analyses support a model whereby growth in +TPA conditions induces a more differentiated state in melanocytes coinciding with an increase in MIR211-5p expression (Figure 4J). This increased MIR211-5p combined with relatively stable MIR328-3p expression, dampens entry into mitosis and/or successful completion of mitosis at least in part through AURKB inhibition. The attenuation of cell division results in an accumulation of cells in G2/M—an effect that is more pronounced in conditions stimulating accelerated growth, such as with BRAFV600E expression. To test this model, we transduced doxycycline-inducible BRAFV600E melanocytes in +TPA conditions with lentiviral constructs constitutively expressing either AURKB (Figure 2I–J) or miRNA sponges (Figure 4K) that inhibit either MIR211-5p or MIR328-3p. AURKB expression had no baseline effect on melanocyte growth, but rescued BRAFV600E-induced arrest (Figure 4L), suggesting the inability of the oncogene to induce AURKB in this context serves as a bottleneck in cell division. Also consistent with the model, sponges against either MIR211-5p or MIR328-3p increased baseline division rate and rescued BRAFV600E-induced arrest, demonstrating that both miRNAs contribute to inhibition of melanocyte division, with effects that are more pronounced in the presence of growth signals.
AURKB inhibition drives melanocytic nevus-associated proliferation arrest
To further investigate whether this mechanism of arrest occurs in clinical nevus specimens, we examined an independent cohort of 15 nevi for evidence of mitotic failure. Histopathologic analysis revealed that all nevi within the cohort exhibited binucleation (Figure 5A). The percent observed binucleated melanocytes ranged from 3% to 12% (mean=6.3%). It is important to note that since binucleation can only be observed in the fraction of cells sectioned perpendicular to the plane of the adjacent nuclei, 3–12% is likely an underestimate. We also observed giant multinucleated cells in 10 of the 15 nevi, indicative of multiple rounds of DNA synthesis and cytokinetic failure and consistent with previous reports (Rogers et al., 2016).
We sought to determine whether AURKB inhibition was necessary and sufficient for the multi-nucleated arrest we observed in melanocytic nevi. We obtained three fresh skin biopsies containing nevi, microdissected the nevus portion, and isolated melanocytes (Figure 5B). Nevus-derived melanocytes appeared healthy in culture. Consistent with the histopathologic cohort, co-staining with melanocyte marker, MLANA, and Hoechst revealed frequent 4 N and 8 N melanocytes (Figure 5C). Transduction of nevus melanocytes with an AURKB and mCherry expressing lentiviral vector at low multiplicity of infection showed that the mCherry positive cells divided significantly more than their mCherry negative counterparts (Figure 5D–E) in cells from all three nevus donors. This result suggests that expression of AURKB is also rate limiting for the proliferation of human nevus melanocytes, as we observed in primary melanocytes expressing BRAFV600E.
To test the sufficiency of AURKB inhibition to limit BRAFV600E hyperproliferation in melanoma, we assessed normalized growth rate inhibition of three cultures—primary human melanocytes, an established BRAFV600E human melanoma cell line, and a short-term BRAFV600E human melanoma culture—when exposed to a pharmacologic AURKB inhibitor (Barasertib). Normalized growth rate measurements (such as GR50) incorporate baseline cell division rate and are demonstrably more reproducible than the commonly calculated half-maximal inhibitor concentration (IC50) when comparing different cell cultures (Hafner et al., 2016). An additional advantage of GR50 is the ability to distinguish growth rate inhibition (GR50 approaches 0) versus cell death (GR50 is negative). Since GR50 calculations require proliferation curves, we conducted live QPI over 48 hr with a range of Barasertib concentrations (0–300 nM). Primary melanocytes were largely resistant to the compound with reduced growth apparent only at supraphysiological levels (Figure 5F). In contrast, both melanoma cultures responded to tenfold lower dosage of the compound with GR50 approaching zero, indicating primarily arrest, not death. In these conditions, the per cell dry mass also increased significantly, consistent with a G2/M arrest (Figure 5G). Taken together, these data demonstrate that AURKB inhibition limits the rate of BRAFV600E-induced proliferation in human melanocytic nevi and melanomas.
Discussion
Cellular senescence is conventionally defined as a permanent cell-cycle arrest associated with a variety of molecular markers and an irreversible arrest in the G0/G1 phase of the cell cycle (He and Sharpless, 2017; Roy et al., 2020; Serrano et al., 1997). The vast majority of melanocytic nevi are stably arrested and harbor BRAFV600E, an oncogenic mutation that induces cellular senescence when expressed in a variety of mammalian cells (Michaloglou et al., 2005; Serrano et al., 1997; Zhu et al., 1998). While these observations seemingly implicate OIS as the mechanism driving nevus-associated proliferation arrest, a substantial body of evidence has challenged this paradigm in recent years (Ruiz-Vega et al., 2020; King et al., 2009; Burian and Jemec, 2019; Cotter et al., 2007; Tran et al., 2012; Oaxaca et al., 2020; Zeng et al., 2018; Haferkamp et al., 2009; McNeal et al., 2015). Here, we provide further evidence to challenge the model that acquisition of BRAFV600E induces premature cellular senescence in melanocytes and uncover an alternative mechanism driving nevus formation and stability. We propose that acquisition of the BRAFV600E mutation does not necessitate growth arrest, but instead permits melanocytes to toggle between hyperproliferation and mitotic arrest. Our transcriptomic analyses demonstrate that external stimuli can modulate the differentiation state of BRAFV600E melanocytes, orchestrate the expression level of MIR211-5p, and inform which phenotype is triggered by the oncogene—before or after the oncogene is introduced.
The senescence model is largely rooted in seminal work demonstrating that sustained ectopic expression of BRAFV600E induces proliferative arrest in primary human melanocytes (Michaloglou et al., 2005)—an experiment reproduced by several independent labs (Haferkamp et al., 2009; McNeal et al., 2015; Leikam et al., 2008; Lu et al., 2016), including our own (Zeng et al., 2018). In each of these studies, melanocytes were cultured in TPA-containing media. Thus, while observations presented here—that BRAFV600E -induced proliferative arrest is reversible and conditional on the presence of TPA— provide new insights that challenge the senescence model, they do not contradict the data presented in these previous reports.
One inherent limitation of using any in vitro approach to investigate the mechanisms underlaying clinical phenomena is an inherent incongruity in the magnitude of duration. Melanocytic nevi frequently persist for several decades, whereas in vitro experiments are generally conducted on the order of days to weeks. In previous in vitro studies, the periods of time melanocytes were monitored subsequent to BRAFV600E introduction ranged from 1 to 3 weeks. In studies utilizing an inducible RAF in human fibroblasts, 24 hr of expression was sufficient to drive a durable proliferative arrest, even upon removal of the oncogene (Zhu et al., 1998). In this study, we induced growth arrest with high levels of BRAFV600E for 30 days and observed the arrest was still reversible. It remains possible that an even longer exposure to BRAFV600E would elicit a permanent proliferative arrest consistent with senescence. It is noteworthy, however, that the identified mechanisms driving proliferative arrest in vitro were also supported by analyses of clinical nevi specimens.
While supported by analyses of both primary human melanocytes and clinical specimens, one limitation to this report is the lack of in vivo confirmation using animal models. One well-established genetically engineered mouse model of melanoma expresses BRAFV600E from the endogenous locus in tyrosinase-expressing cells with spatiotemporal control. The predominant phenotype of BRAFV600E activation in this model is melanocytic hyperplasia, which increases in prominence with BRAFV600E dose (Dankort et al., 2009). While growth-arrested melanocytic ‘nevus-like’ lesions also form in this model, they are not associated with markers of senescence (Ruiz-Vega et al., 2020). More recently, an elegant study utilizing a zebrafish model of melanoma demonstrated that driving BRAFV600E expression from promoters associated with melanocyte progenitors, but not differentiated melanocytes, resulted in melanoma initiation (Arianna et al., 2021). These observations are all consistent with our data and support the conclusion that BRAFV600E-induced proliferation in melanocytes is dependent on differentiation state. Further studies will be required to confirm whether the microRNA mechanism reported here is conserved between human and model systems.
The presented model offers an explanation for several clinical phenomena. First, melanocytic nevi form when BRAFV600E drives hyperproliferation followed by proliferation arrest. Second, nevus eruption involves the expansion of previously stable nevi over a short period of time. Third, incomplete excision of melanocytic nevi can result in subsequent regrowth of the benign lesion (King et al., 2009). Each event requires interconversion of the BRAFV600E melanocyte between a proliferative and arrested phase. Our data demonstrate that BRAFV600E melanocytes can indeed oscillate between proliferative and arrested states and are consistent with each of these phenomena. Our interpretation assumes that once acquired, BRAFV600E expression is constant and that the environmental context of the melanocytes is variable, serving as a ‘gatekeeper’ to BRAFV600E-induced proliferation. However, if BRAFV600E expression also fluctuates, then an equally plausible interpretation of our data is that BRAFV600E serves as a ‘gatekeeper’ to environmentally induced proliferation. Interestingly, nevus eruption has been observed in metastatic melanoma patients treated with single agent selective BRAF inhibitors (Cohen et al., 2013; Chen et al., 2014; Chu et al., 2012; Dalle et al., 2013; Zimmer et al., 2012), consistent with our observation that proliferation arrest in melanocytes requires sustained expression of the oncogene. A more thorough understanding of the responsible environmental stimuli in physiologic conditions, discussed in more detail below, would help to clarify which of the mechanisms—variable BRAFV600E, variable environment, or both—underlie nevus formation, eruption, and recurrence.
It is important to note that while our data support a model of reversible proliferation arrest in nevi, we do not suggest this reversal is commonplace. Unquestionably, the vast majority of melanocytic nevi present with sustained arrest. However, we do propose that the duration of arrest is dynamic and modulated by exposure to external signals. Since our data suggest that loss of MIR211-5p and MIR328-3p expression and restoration of AURKB expression is concurrent with progression to melanoma, identification of the source of cell-extrinsic signals might suggest new strategies for chemoprevention or therapy. We have identified TPA, a known activator of PKC signaling, as one extrinsic factor that induces expression of MIR211-5p and a more dedifferentiated state. Previous studies have identified TPA as either a mitogen, inhibitor of the G1/S transition, or inhibitor of the G2/M transition, dependent on the context (Coppock et al., 1995; Arita et al., 1998; Arita et al., 1992; Stavroulaki et al., 2008; Chao-Hsing and Hsin-Su, 1991). We demonstrate that in vitro, TPA activation of MIR211-5p is required for BRAFV600E-induced mitotic failure. Yet, the inclusion of TPA in primary melanocyte culture is artificial and the environmental stimulus that regulates MIR211-5p expression in skin remains unknown. One hypothesis is that PKC activation from keratinocyte interactions serves as the extrinsic stimulus. MIR211-5p is reported to be transcriptionally regulated by both MITF (Mazar et al., 2010) and UV exposure (Su et al., 2020). However, TPA-free media containing an endothelin receptor type B ligand, endothelin-1, did not elicit MIR211-5p expression or BRAFV600E arrest, suggestive that the downstream signaling pathways regulated by the two molecules are not identical. PKC activation can also result in increased MAPK signaling (Herrera et al., 1998; Ueda et al., 1996) raising the possibility that internal MAPK fluctuations may play a role in toggling between BRAFV600E-induced proliferation versus arrest. The term ‘mitogenic window’ describes the specific range of oncogene doses that elicit a hyperproliferative phenotype. The mitogenic window for MAPK signaling activation can be remarkably fine. Here, we conducted a twofold dilution series of doxycycline but observed no evidence of a mitogenic window. However, it remains possible that we simply missed the ‘sweet spot’ in TPA conditions, and that a finer dilution series would reveal a mitogenic window even in the presence of TPA. Regardless, our data demonstrate that in TPA conditions the mitogenic window for BRAFV600E expression is, at least, exceedingly narrow, whereas in non-TPA conditions the window is wide open. An important future direction will be to identify the nature and physiological source of this signal for melanocytic nevi in human skin.
Diminished MIR211-5p/MIR328-3p expression concurrent with increased AURKB expression occurs in a substantial portion of melanomas. AURKB expression was previously reported as significantly elevated in melanoma metastases (Sharma et al., 2013) and here we report increased expression in primary melanomas. Similarly, MIR211-5p expression is consistently observed as lost in melanoma as compared to melanocytic nevi across independent studies (Babapoor et al., 2016; Torres et al., 2020; Latchana et al., 2016). Indeed, MIR211-5p expression levels classify nevi from melanomas with a high degree of accuracy (Babapoor et al., 2016; Torres et al., 2020). We have reported on the diagnostic potential of MIR328-3p expression as well, although multi-study validation was impeded by the absence of hybridization probes in earlier microarray-based profiling data sets (Torres et al., 2020). In vitro, we have now demonstrated that inhibition of either miRNA is sufficient to block BRAFV600E-induced arrest in melanocytes. These data suggest that mitotic failure serves as a barrier to melanoma progression and is consistent with the genome duplication and copy number alterations associated with invasive melanoma (Shain et al., 2015; Vergara et al., 2021). In melanoma, the specific cause of this genomic instability is unknown, but reentry of cells into cell cycle after mitotic failure has been shown to cause genome duplication and accumulation of copy number alterations in other cancers (Telentschak et al., 2015). Further experiments will be necessary to determine whether dysregulation of the spindle checkpoint in nevi contributes directly to the copy number variations present in melanomas.
In summary, we have identified that BRAFV600E induces a reversible arrest in human melanocytes orchestrated by MIR211-5p/MIR328-3p regulation of AURKB and conditional on the melanocyte differentiation state. We present a mechanistic model that allows for melanocytic nevus formation, eruption, and recurrence. Moreover, our observations provide an explanation for the genetic duplication and instability inherent to invasive melanomas and open new directions for novel chemopreventative and therapeutic strategies.
Materials and methods
Transcriptomic profiling
Request a detailed protocolFor FFPE samples, histopathologic review, microdissection, targeted exon sequencing, phylogenetic analysis, and RNA and small RNA sequencing of each tumor area were previously described (phs001550.v2.p1) (Zeng et al., 2018; Shain et al., 2018; Torres et al., 2020). For this study, small RNA reads were aligned to human reference (hg37) with Bowtie (Langmead et al., 2009) and small RNA reference groups (miRBase21) were counted. For mRNA sequencing, reads were aligned to human reference (hg37) with Bowtie2 and counted with HTSeq. DE analysis was performed using DESeq2 (Love et al., 2014) with p values adjusted by the Benjamini-Hochberg method (p-adj). Previously published RNA sequencing data sets from Shain et al., 2018 and Torres et al., 2020 were re-analyzed here with DESeq2 and combined for visualization. For RNA profiling from primary melanocytes, total RNA was extracted using TRIzol Reagent (Thermo Fisher Scientific, 15596-026) and further purified with the RNeasy Power Clean Pro Cleanup Kit (Qiagen, 13997-50) to remove melanin. For small RNAseq of nevus melanocytes sequencing libraries were constructed with the TailorMix Small RNA Library Preparation Kit (SeqMatic, CA) and sequencing was performed on the Illumina HiSeq2500 platform at single-end 50 bp. For small RNAseq of cultured BRAFV600E transduced melanocytes, sequencing libraries were constructed with the Qiagen QIAseq miRNA Library Prep Kit, and sequencing was performed on the NovaSeq SP platform at paired-end 50 bp. For mRNAseq of microRNA mimic nucleofected melanocytes, 150 bp paired-end sequencing was conducted by GeneWiz. For mRNAseq of cultured BRAFV600E transduced melanocytes, sequencing libraries were constructed with the Illumina TruSeq Stranded mRNA Library Prep Kit, and sequencing was performed on the NovaSeq S4 platform at paired-end 150 bp.
Cell derivation and culture
Request a detailed protocolBenign human melanocytic nevi were excised with informed consent from patient donors at the UCSF Dermatology clinic (San Francisco, CA) or HCI Dermatology clinic (Salt Lake City, UT) according to Institutional Review Board approved protocols. Nevus tissue was minced and enzymatically dissociated in 1 mg/ml collagenase (Sigma-Aldrich, 11213857001) and 3.3 mg/ml dispase (Sigma-Aldrich, 4942078001) in DMEM (Thermo Fisher Scientific, 10569044) for 1 hr at 37°C. Melanocytes were further isolated by 5 days exposure to 10 µg/ml G418 (InvivoGen, ant-gn-1). BRAF status was confirmed via sanger sequencing (Quintarabio) using the primer set: BRAF forward: 5′-GCA CGA CAG ACT GCA CAG GG -3′; BRAF reverse: 5′-AGC GGG CCA GCA GCT CAA TAG-3′. BRAF wild-type (normal) human melanocytes were isolated from de-identified and IRB consented neonatal foreskins or adult skin. Foreskin tissue was incubated overnight at 4°C in dispase and epithelia were mechanically separated from the dermis the following morning. Epithelial tissue was minced and incubated in 0.25% trypsin (Gibco 25200056) for 4 min at 37°C. Trypsin was quenched and tissue centrifuged at 500×g for 5 min at room temperature. The cell/tissue pellet was resuspended in Melanocyte medium (Thermo Fisher Scientific, M254500) with HMGS (Thermo Fisher Scientific, S0025) and plated in low volume. Melanocytes were expanded in Melanocyte medium with HMGS and assayed in either Melanocyte medium with HMGS (+TPA) or Melanocyte medium with HMGS-2 (TPA-free, Thermo Fisher Scientific, S0165) as indicated.
The 501Mel human melanoma line (Gift from Dr. Boris Bastian, CVCL_4633, authenticated via STR profiling) were cultured in RPMI with 10% fetal bovine serum (FBS; VWR, S107G), 1% penicillin-streptomycin (Gibco, 15140122), 1% L-glutamine (Gibco, 25030149), and 1% non-essential amino acids (Gibco, 11140050). HCIMel019 was derived from patient-derived xenograft tumors propagated in mice. An HCIMel019 (P2) subcutaneous tumor was resected from mouse, minced in digestion buffer (100 μM HEPES (Gibco, 15630-080), 5% FBS (DENVILLE, FB5001-H), 20 μg/ml gentamicin (Gibco, 15710-064), 1× insulin (Gibco, 51500-056), and 1 mg/ml collagenase IV (Gibco, 17104-019) in DMEM (Gibco, 11965-092)), and digested overnight at 37°C. Cells were filtered through a 100 µm filter (Falcon, 352360) and red blood cells were removed with RBC lysis buffer (0.5 M EDTA, 0.5 M KHCO3 [Sigma-Aldrich, 237205], and 5 M NH4CL [Sigma-Aldrich, A9434]). Remaining cells were washed with phosphate-buffered saline (PBS; Gibco, 10010023), and cultured at 37°C and 5% CO2 in Mel2 media consisting of 80% MCDB153 (Sigma-Aldrich, M7403), 20% L15 (Gibco, 11415-064), 2% FBS (DENVILLE, FB5001-H), 1.68 mM CaCL (Sigma-Aldrich, C4901), 1× insulin (Gibco, 51500-056) 5 ng/ml EGF (Sigma-Aldrich, E9644), 15 μg/ml Bovine Pituitary Extract (Gibco 13028-014), and 1× Pen/strep (Gibco, 15070-063). Cell cultures were tested monthly for mycoplasma contamination (ATCC).
Small RNA nucleofection
Request a detailed protocolHuman melanocytes were trypsinized (Gibco, 25300062), quenched, and centrifuged at 300×g. Cells were resuspended in R Buffer (Thermo Fisher Scientific, MPK1025) at 10,000 cells/µL. About 10 µl of cell slurry was mixed with miRIDIAN microRNA mimics (Dharmacon: hsa-MIR211-5p C-300566-03-0005, hsa-MIR328-3p C-300695-03-0005, MIRControl 1 CN-001000-01-05, or MIRControl 2 CN-002000-01-05); or On-TargetPlus siRNA (Dharmacon, custom library); at 4 µM final concentration and nucleofected using the NEON Transfection System and protocol (Thermo Fisher Scientific, MPK5000).
Generation of Lentiviral vectors
Request a detailed protocolpTRIPZ-diBRAFV600E was a gift from Todd Ridky (McNeal et al., 2015). pLVX-AURKB (Addgene #153316) and pLVX-GPR3 (Addgene #153317) were generated by subcloning the respective human cDNA (from Addgene #100142 and #66350) into the MluI and BamHI sites of the pLVX-Che-hi3 vector (a gift of Sanford Simon) (Takacs et al., 2017). pLVX-anti-MIR211-5p (Addgene #153318), pLVX-anti-MIR328-3p (Addgene #153319), and pLVX-Che-zsGreen (Addgene #153320) were generated by inserting zsGreen with or without a 3′UTR into the MluI and XbaI sites of the pLVX-Che-hi3 vector. The 3′UTRs contained seven tandem microRNA binding sites to report and inhibit microRNA function, as previously described (Judson et al., 2013). ZsGreen was subcloned from pHIV-zsGreen gift from Bryan Welm & Zena Werb (Addgene #18121) (Welm et al., 2008).
Lentiviral transduction
Request a detailed protocol2.75×106 HEK293T cells were plated on 10 cm tissue culture dishes and grown for ~24 hr in DMEM with 10% FBS. For each 10 cm plate, 5 µg lentiviral vector, 3.3 µg of pMLV-GagPol, and 1.7 µg of pVSV-G packaging plasmids were added to 500 µl of jetPRIME buffer and 20 µl of jetPRIME transfection reagent (Polyplus, 712-60), and transfected according to the manufacturer’s instructions. 48 hr post-transfection, viral supernatant was collected and filtered using 0.45 µm syringe filters (Argos 4395-91). Human melanocytes seeded at 1.0×105–2.0×105 cells/well density in six-well plates were incubated in viral supernatant with 10 µg/ml polybrene (Sigma-Aldrich, TR-1003) and centrifuged at 300×g for 60 min at room temperature. The viral supernatant was removed and replaced with growth media. Transduced cells were either selected with puromycin (1 µg/mL for 5 days, Sigma-Aldrich P8833) or sorted for mCherry expression using a BD FACSAria II. Cells expressing pTRIPZ-diBRAFV600E were treated with doxycycline (Sigma-Aldrich, D9891) at indicated concentrations.
Live quantitative phase imaging
Request a detailed protocolLive quantitative imaging was performed using either the HoloMonitor M4 imaging cytometer (Phase Holographic Imaging, Lund, Sweden) or the Livecyte platform (Phasefocus, Sheffield, UK). Analyses of cell proliferation, dry cell mass, and death were acquired with the M4 platform and analyzed using HStudio (v2.6.3) as previously described (Hejna et al., 2017). For each experiment, human melanocytes were seeded into six-well plates (Sarstedt, 83.3920) at 100,000 cells/well and either live-imaged or serially imaged as indicated. Analysis of growth rate was acquired with the M4 platform and analyzed using App Suite (v3.2.0.60) as previously described (Hafner et al., 2016). For each experiment, 100,000 melanocytes, 60,000 501Mel cells, or 150,000 HCIMel019 cells were plated per well and media containing indicated concentrations of barasertib (AURKB inhibitor; AZD1152-HQPA | AZD2811, Selleckchem, A1147) was added. Cells were imaged for 48–72 hr. Analyses of fluorescent reporters coupled with QPI were conducted using the Livecyte platform and were analyzed using Analyze (v3.1) (Phasefocus, Sheffield, UK).
EdU assays
Request a detailed protocolNormal human melanocytes were nucleofected with microRNA mimicsor infected with lentiviral vectors (as above) and seeded in 48-well plates at a density of 50,000 cells/well. Four days post-seeding, EdU was added to culture media at a 10 µM final concentration. 24 hr after EdU addition, media was removed and cells were stained for EdU incorporation and nuclei using the Click-iT EdU Imaging Kit (Thermo Fisher Scientific, C10337) and protocol. Images of EdU and nuclei staining were acquired using the Evos FL microscope and quantified with Fiji.
Western blotting
Request a detailed protocolProtein was collected using RIPA Buffer (Thermo Fisher Scientific, 89901) supplemented with HALT Protease and Phosphatase Inhibitor (Thermo Fisher Scientific, 78446). Immunoblotting was carried out as previously described. Membranes were incubated overnight at 4°C with primary antibodies at the following dilutions: anti-HSP90 (CST, 4874) 1:1000, anti-BRAFV600E (Spring Bioscience Corp, E19292) 1:1000, anti-phosphoERK1/2 (CST, 4970) 1:1000, and anti-AURKB (Abcam, ab2254) 1:1000. Membranes were washed 4× with TBS and 0.5% Tween20, incubated with HRP-conjugated secondary antibody (1:2000) for 30 min at room temperature, and visualized with Lumina Forte Western HRP substrate (Millipore, WBLUF0500).
Flow analyses
Request a detailed protocolFor assessing DNA content, trypsinized melanocytes (as above) were resuspended in 400 PBS. 1 ml –20°C 200 proof ethanol was added dropwise while gently vortexing to achieve 70% final concentration for fixation and incubated overnight at –30°C. Cells were centrifuged at 800×g for 5 min, washed with fresh PBS, and incubated in 500 µL FxCycle PI/RNase Staining Solution (Invitrogen, F10797) for 20 min at 37°C. Data were collected on the Fortessa (BD) at low speed and analyzed with FlowJo v10.7.1. Annexin V/PI staining was performed by trypsinizing melanocytes as above. Cells were washed with PBS and resuspended with Annexin V binding buffer (BioLegend, 422201) at a concentration of 1×106 cell/ml. 5 µl APC Annexin V antibody (BioLegend, 640919) and 10 µl PI were added. Cells were gently vortexed and incubated in the dark at room temperature for 15 min. Cells were resuspended in 400 µl of Annexin V Binding Buffer. Data were collected on the FACs Verse (BD) and were analyzed with FlowJo v10.7.1. Aldehyde dehydrogenase activity was measured first by trypsinizing melanocytes as above and then using the ALDEFLUOR Kit (Stemcell Technologies, 01700) as described in the manufacture’s protocol.
Human nevus and melanoma tissue immunohistochemistry
Request a detailed protocolThe archived clinical specimens used in the study were procured with IRB approval from the UCSF Dermatopathology Pathology Service archives. For assessment of AURKB and GPR3 expression, 11 melanoma and 11 melanocytic nevi with previous diagnosis were de-identified and re-verified histologically (UEL). For assessment of DNA content, 15 additional melanocytic nevi with previous diagnosis were re-verified histologically (UEL). Tissue was fixed in 10% neutral-buffered formalin, processed, embedded in paraffin, and stained with hematoxylin and eosin. 4 µm FFPE sections were stained with AURKB (1:400 dilution, Abcam ab2254) or GPR3 (1:125 dilution, Abnova H00002827-M01) antibodies. UEL reviewed immunohistochemical stains with semiquantitative grading for GPR3 (0=none; 1=patchy positive; 2=strong positive) and AURKB (0=none; 1=rare positive; 2=scattered positive; 3 = frequent positive; 4=many positive). For in vitro immunofluorescence, cells were fixed and stained as previously described with AXL (1:250 dilution, Cell Signaling Technology #8661) or MLANA (1:100, Abcam ab731).
Statistical analyses
Request a detailed protocolFor differential gene expression, p values were calculated with the DESeq2 (v1.30.1) default Wald test adjusted by the Benjamini-Hochberg method using a 5% false discovery rate (FDR) (Benjamini et al., 2001). Pathway analyses were analyzed using the fast gene set enrichment package (Korotkevich et al., 2016) in R with a 10% FDR. Gene set enrichment analyses were conducted using GSEA (v4.0.3, Broad Institute) Preranked tool with 1000 permutations. For in vitro experiments, pilot studies were initially conducted in triplicate. The required minimal sample size to assess the difference between independent means with independent standard deviations with an alpha error probability of 0.05 and a power of 0.95 was calculated using G*Power (v. 3.1.9.4). P values were calculated using either paired or unpaired two-tailed t-tests or Wilcoxon tests via Prism 8 (GraphPad) as indicated in figure legends.
Experimental set-up
Request a detailed protocolFor cohorts of clinical specimens, sample number (n) refers to the number of specimens, each from a different patient. For in vitro experiments, n refers to independent experiments conducted on different days. In cases where fresh primary melanocytes were used, the cells in each experiment are derived from an independent donor. In cases where cell lines are used, the cells in each experiment represent a different passage and/or day of that line. For different conditions within experiments, replicate wells were plated and arbitrarily chosen for each condition (e.g., not treated versus treated; different concentrations of a compounds, etc.). All data were included.
Data availability
The sequencing data produced in this manuscript is available from GEO under the accession number GSE150849. Source data for other large datasets are provided as supplementary files. Other data are included in the manuscript and figures.
-
NCBI Gene Expression OmnibusID GSE150849. microRNAs Restrain Proliferation in BRAFV600E Melanocytic Nevi.
-
NCBI Gene Expression OmnibusID phs001550.v2.p1. The genetic and transcriptomic evolution of melanoma.
References
-
vestigation of the role of G1/S cell cycle mediators in cellular senescenceExperimental Cell Research 209:231–237.https://doi.org/10.1006/excr.1993.1306
-
Growth of human melanocyte cultures supported by 12-O-tetradecanoylphorbol-13-acetate is mediated through protein kinase C activationCancer Research 52:4514–4521.
-
The molecular pathology of melanoma: An integrated taxonomy of melanocytic neoplasiaAnnual Review of Pathology 9:239–271.https://doi.org/10.1146/annurev-pathol-012513-104658
-
Human melanocyte development and melanoma dedifferentiation at single-cell resolutionNature Cell Biology 23:1035–1047.https://doi.org/10.1038/s41556-021-00740-8
-
Controlling the false discovery rate in behavior genetics researchBehavioural Brain Research 125:279–284.https://doi.org/10.1016/s0166-4328(01)00297-2
-
Eruptive Melanocytic Nevi: A ReviewAmerican Journal of Clinical Dermatology 20:669–682.https://doi.org/10.1007/s40257-019-00444-8
-
Cellular senescence: When bad things happen to good cellsNature Reviews. Molecular Cell Biology 8:729–740.https://doi.org/10.1038/nrm2233
-
A study of the effects of phorbol 12-myristate-13-acetate on cell differentiation of pure human melanocytes in vitroArchives of Dermatological Research 283:119–124.https://doi.org/10.1007/BF00371620
-
Diverse cutaneous side effects associated with BRAF inhibitor therapy: A clinicopathologic studyJournal of the American Academy of Dermatology 67:1265–1272.https://doi.org/10.1016/j.jaad.2012.04.008
-
Appearance of New Vemurafenib-associated Melanocytic Nevi on Normal-appearing Skin: Case Series and a Review of Changing or New Pigmented Lesions in Patients with Metastatic Malignant Melanoma After Initiating Treatment with VemurafenibThe Journal of Clinical and Aesthetic Dermatology 6:27–37.
-
The power and the promise of oncogene-induced senescence markersNature Reviews. Cancer 6:472–476.https://doi.org/10.1038/nrc1884
-
Absence of senescence-associated beta-galactosidase activity in human melanocytic nevi in vivoThe Journal of Investigative Dermatology 127:2469–2471.https://doi.org/10.1038/sj.jid.5700903
-
Tracking of second primary melanomas in vemurafenib-treated patientsJAMA Dermatology 149:488–490.https://doi.org/10.1001/jamadermatol.2013.21
-
Braf(V600E) cooperates with Pten loss to induce metastatic melanomaNature Genetics 41:544–552.https://doi.org/10.1038/ng.356
-
MiRWalk database for miRNA-target interactionsMethods in Molecular Biology 1182:289–305.https://doi.org/10.1007/978-1-4939-1062-5_25
-
Oncogene-induced senescence does not require the p16(INK4a) or p14ARF melanoma tumor suppressorsThe Journal of Investigative Dermatology 129:1983–1991.https://doi.org/10.1038/jid.2009.5
-
Molecular pathways: BRAF induces bioenergetic adaptation by attenuating oxidative phosphorylationClinical Cancer Research 20:2257–2263.https://doi.org/10.1158/1078-0432.CCR-13-0898
-
Barasertib (AZD1152), a small molecule Aurora B inhibitor, inhibits the growth of SCLC cell lines in vitro and in vivoMolecular Cancer Therapeutics 15:2314–2322.https://doi.org/10.1158/1535-7163.MCT-16-0298
-
Cancer stem cells versus phenotype-switching in melanomaPigment Cell & Melanoma Research 23:746–759.https://doi.org/10.1111/j.1755-148X.2010.00757.x
-
Cultivation of normal human epidermal melanocytes in the absence of phorbol estersMethods in Molecular Medicine 107:13–28.https://doi.org/10.1385/1-59259-861-7:013
-
MicroRNA-based discovery of barriers to dedifferentiation of fibroblasts to pluripotent stem cellsNature Structural & Molecular Biology 20:1227–1235.https://doi.org/10.1038/nsmb.2665
-
KEGG as a reference resource for gene and protein annotationNucleic Acids Research 44:457–462.https://doi.org/10.1093/nar/gkv1070
-
MicroRNA dysregulation in melanomaSurgical Oncology 25:184–189.https://doi.org/10.1016/j.suronc.2016.05.017
-
Oncogenic BRAF-Mediated Melanoma Cell InvasionCell Reports 15:2012–2024.https://doi.org/10.1016/j.celrep.2016.04.073
-
MicroRNAs can generate thresholds in target gene expressionNature Genetics 43:854–859.https://doi.org/10.1038/ng.905
-
Understanding cytokinesis failureAdvances in Experimental Medicine and Biology 675:27–55.https://doi.org/10.1007/978-1-4419-6199-0_3
-
p16 Range of expression in dermal predominant benign epithelioid and spindled nevi and melanomaJournal of Cutaneous Pathology 47:815–823.https://doi.org/10.1111/cup.13721
-
Biologically distinct subsets of neviGiornale Italiano Di Dermatologia e Venereologia 151:365–384.
-
Global Analysis of BRAFV600E Target Genes in Human Melanocytes Identifies Matrix Metalloproteinase-1 as a Critical Mediator of Melanoma GrowthThe Journal of Investigative Dermatology 131:1579–1583.https://doi.org/10.1038/jid.2011.65
-
The Genetic Evolution of Melanoma from Precursor LesionsThe New England Journal of Medicine 373:1926–1936.https://doi.org/10.1056/NEJMoa1502583
-
Identification of aurora kinase B and Wee1-like protein kinase as downstream targets of (V600E)B-RAF in melanomaThe American Journal of Pathology 182:1151–1162.https://doi.org/10.1016/j.ajpath.2012.12.019
-
Exposure of normal human melanocytes to a tumor promoting phorbol ester reverses growth suppression by transforming growth factor betaJournal of Cellular Physiology 214:363–370.https://doi.org/10.1002/jcp.21207
-
Role of the p53‑TRPM1/miR‑211‑MMP9 axis in UVB‑induced human melanocyte migration and its potential in repigmentationTernational Journal of Molecular Medicine 45:1017–1026.https://doi.org/10.3892/ijmm.2020.4478
-
MicroRNA Ratios Distinguish Melanomas from NeviThe Journal of Investigative Dermatology 140:164–173.https://doi.org/10.1016/j.jid.2019.06.126
-
Absence of distinguishing senescence traits in human melanocytic neviThe Journal of Investigative Dermatology 132:2226–2234.https://doi.org/10.1038/jid.2012.126
-
Protein kinase C activates the MEK-ERK pathway in a manner independent of Ras and dependent on RafThe Journal of Biological Chemistry 271:23512–23519.https://doi.org/10.1074/jbc.271.38.23512
-
Phenotypic plasticity of melanocytes derived from human adult skinPigment Cell & Melanoma Research 31:1–14.https://doi.org/10.1111/pcmr.13012
-
V599EB-RAF is an Oncogene in MelanocytesCancer Research 64:2338–2342.https://doi.org/10.1158/0008-5472.CAN-03-3433
-
Senescence of human fibroblasts induced by oncogenic RafGenes & Development 12:2997–3007.https://doi.org/10.1101/gad.12.19.2997
Article and author information
Author details
Funding
National Institutes of Health (DP5OD019787)
- Robert L Judson-Torres
National Cancer Institute (R01CA229896)
- Robert L Judson-Torres
National Cancer Institute (5F31CA236377)
- Andrew S McNeal
National Cancer Institute (P30CA042014)
- Robert L Judson-Torres
Program for Breakthrough Biomedical Research (Sandler Fellowship)
- Robert L Judson-Torres
5 for the Fight (5 for the Fight Fellowship)
- Robert L Judson-Torres
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors acknowledge the use of the HCI Shared Resources for Research Informatics (RI), Cancer Biostatistics (CB), High-Throughput Genomics and Bioinformatics Analysis (GBA), and the Biorepository and Molecular Pathology (BMP) Research Histology Section supported by P30CA042014 awarded to HCI from the National Cancer Institute.
Ethics
Human subjects: The archived clinical specimens used in the study were procured with IRB approval from the UCSF Dermatopathology Pathology Service archives (11-07951). Fresh human tissue samples used in this study were obtained with informed consent and consent to publish from each patient. Samples were first coded before being provided to the researchers. The study protocol conformed to the ethical guidelines of the Declaration of Helsinki (1975). Samples were processed following standard operating procedures approved by IRB at University of California, San Francisco (16-20299) and at the Huntsman Cancer Institute (00124195).
Copyright
© 2021, McNeal et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,540
- views
-
- 431
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
Clonal hematopoiesis of indeterminate potential (CHIP) allows estimation of clonal dynamics and documentation of somatic mutations in the hematopoietic system. Recent studies utilizing large cohorts of the general population and patients have revealed significant associations of CHIP burden with age and disease status, including in cancer and chronic diseases. An increasing number of cancer patients are treated with immune checkpoint inhibitors (ICIs), but the association of ICI response in non-small cell lung cancer (NSCLC) patients with CHIP burden remains to be determined. We collected blood samples from 100 metastatic NSCLC patients before and after ICI for high-depth sequencing of the CHIP panel and 63 samples for blood single-cell RNA sequencing. Whole exome sequencing was performed in an independent replication cohort of 180 patients. The impact of CHIP status on the immunotherapy response was not significant. However, metastatic lung cancer patients showed higher CHIP prevalence (44/100 for patients vs. 5/42 for controls; p = 0.01). In addition, lung squamous cell carcinoma (LUSC) patients showed increased burden of larger clones compared to lung adenocarcinoma (LUAD) patients (8/43 for LUSC vs. 2/50 for LUAD; p = 0.04). Furthermore, single-cell RNA-seq analysis of the matched patients showed significant enrichment of inflammatory pathways mediated by NF-κB in myeloid clusters of the severe CHIP group. Our findings suggest minimal involvement of CHIP mutation and clonal dynamics during immunotherapy but a possible role of CHIP as an indicator of immunologic response in NSCLC patients.
-
- Cancer Biology
The prevalence and mortality rates of colorectal cancer (CRC) are increasing worldwide. Radiation resistance hinders radiotherapy, a standard treatment for advanced CRC, leading to local recurrence and metastasis. Elucidating the molecular mechanisms underlying radioresistance in CRC is critical to enhance therapeutic efficacy and patient outcomes. Bioinformatic analysis and tumour tissue examination were conducted to investigate the CPT1A mRNA and protein levels in CRC and their correlation with radiotherapy efficacy. Furthermore, lentiviral overexpression and CRISPR/Cas9 lentiviral vectors, along with in vitro and in vivo radiation experiments, were used to explore the effect of CPT1A on radiosensitivity. Additionally, transcriptomic sequencing, molecular biology experiments, and bioinformatic analyses were employed to elucidate the molecular mechanisms by which CPT1A regulates radiosensitivity. CPT1A was significantly downregulated in CRC and negatively correlated with responsiveness to neoadjuvant radiotherapy. Functional studies suggested that CPT1A mediates radiosensitivity, influencing reactive oxygen species (ROS) scavenging and DNA damage response. Transcriptomic and molecular analyses highlighted the involvement of the peroxisomal pathway. Mechanistic exploration revealed that CPT1A downregulates the FOXM1-SOD1/SOD2/CAT axis, moderating cellular ROS levels after irradiation and enhancing radiosensitivity. CPT1A downregulation contributes to radioresistance in CRC by augmenting the FOXM1-mediated antioxidant response. Thus, CPT1A is a potential biomarker of radiosensitivity and a novel target for overcoming radioresistance, offering a future direction to enhance CRC radiotherapy.