Cancers contain a wide diversity of cell types that are defined by differentiation states, genetic mutations and altered epigenetic programmes that impart functional diversity to individual cells. Elevated tumour cell heterogeneity is linked with progression, therapy resistance and relapse. Yet, imaging of tumour cell heterogeneity and the hallmarks of cancer has been a technical and biological challenge. Here we develop optically clear immune-compromised rag2(E450fs) (casper) zebrafish for optimized cell transplantation and direct visualization of fluorescently labelled cancer cells at single-cell resolution. Tumour engraftment permits dynamic imaging of neovascularization, niche partitioning of tumour-propagating cells in embryonal rhabdomyosarcoma, emergence of clonal dominance in T-cell acute lymphoblastic leukaemia and tumour evolution resulting in elevated growth and metastasis in BRAF(V600E)-driven melanoma. Cell transplantation approaches using optically clear immune-compromised zebrafish provide unique opportunities to uncover biology underlying cancer and to dynamically visualize cancer processes at single-cell resolution in vivo.

Multiplexing samples in sequencing experiments is a common approach to maximize information yield while minimizing cost. In most cases the number of samples that are multiplexed is determined by financial consideration or experimental convenience, with limited understanding on the effects on the experimental results. Here we set to examine the impact of multiplexing ChIP-seq experiments on the ability to identify a specific epigenetic modification. We performed peak detection analyses to determine the effects of multiplexing. These include false discovery rates, size, position and statistical significance of peak detection, and changes in gene annotation. We found that, for histone marker H3K4me3, one can multiplex up to 8 samples (7 IP + 1 input) at ~21 million single-end reads each and still detect over 90% of all peaks found when using a full lane for sample (~181 million reads). Furthermore, there are no variations introduced by indexing or lane batch effects and importantly there is no significant reduction in the number of genes with neighboring H3K4me3 peaks. We conclude that, for a well characterized antibody and, therefore, model IP condition, multiplexing 8 samples per lane is sufficient to capture most of the biological signal.

The efficacy of the highly selective RET inhibitor selpercatinib is now established in RET-driven cancers, and we sought to characterize the molecular determinants of response and resistance. We find that the pre-treatment genomic landscape does not shape the variability of treatment response except for rare instances of RAS-mediated primary resistance. By contrast, acquired selpercatinib resistance is driven by MAPK pathway reactivation by one of two distinct routes. In some patients, on- and off-target pathway reactivation via secondary RET solvent front mutations or MET amplifications are evident. In other patients, rare RET-wildtype tumor cell populations driven by an alternative mitogenic driver are selected for by treatment. Multiple distinct mechanisms are often observed in the same patient, suggesting polyclonal resistance may be common. Consequently, sequential RET-directed therapy may require combination treatment with inhibitors targeting alternative MAPK effectors, emphasizing the need for prospective characterization of selpercatinib-treated tumors at the time of monotherapy progression.

The two major melanoma histologic subtypes, superficial spreading and nodular melanomas, differ in their speed of dermal invasion but converge biologically once they invade and metastasize. Herein, we tested the hypothesis that distinct molecular alterations arising in primary melanoma cells might persist as these tumors progress to invasion and metastasis. Ribosomal protein S6 kinase, 90 kDa, polypeptide 1 (RSK1; official name RPS6KA1) was significantly hyperactivated in human melanoma lines and metastatic tissues derived from nodular compared with superficial spreading melanoma. RSK1 was constitutively phosphorylated at Ser-380 in nodular but not superficial spreading melanoma and did not directly correlate with BRAF or MEK activation. Nodular melanoma cells were more sensitive to RSK1 inhibition using siRNA and the pharmacological inhibitor BI-D1870 compared with superficial spreading cells. Gene expression microarray analyses revealed that RSK1 orchestrated a program of gene expression that promoted cell motility and invasion. Differential overexpression of the prometastatic matrix metalloproteinase 8 and tissue inhibitor of metalloproteinases 1 in metastatic nodular compared with metastatic superficial spreading melanoma was observed. Finally, using an in vivo zebrafish model, constitutive RSK1 activation increased melanoma invasion. Together, these data reveal a novel role for activated RSK1 in the progression of nodular melanoma and suggest that melanoma originating from different histologic subtypes may be biologically distinct and that these differences are maintained as the tumors invade and metastasize.