A metastasis map of human cancer cell lines

1Human cancer cell lines have been a driving force in cancer research, leading to the discovery of oncogenic mechanisms and therapeutic targets–. However, large-scale characterization of cell lines has been limited to rudimentary readouts such as viability in cell culture, because more complex phenotypes—such as behaviours in vivo—have not been tractable at scale. By contrast, most studies of metastasis rely on only a small number of experimental models–, thereby making it difficult to extrapolate findings to genetically diverse human tumours.

2Ideally, it would be possible to construct a map of organ-specific metastatic potential of hundreds of human cancer cell lines using xenograft models, so that the molecular features of the cell lines could be related to their ability to survive and proliferate in organ-specific microenvironments. However, the prospect of in vivo testing of each cell line individually is unattractive, because it is labour-intensive and expensive, as well as because of the difficulty in sufficiently controlling for variability between animal experiments. We proposed that if cell lines were labelled with molecular barcodes and injected into recipient mice as a pool, internally controlled, metastatic potential could be assessed in a highly scalable manner.

3To test the feasibility and reliability of in vivo barcoding to monitor growth in different tissues in mice, we performed a pilot study using four breast cancer cell lines (Fig. 1a, Extended Data Fig. 1, Supplementary Note 1). Each cell line was engineered to express a unique 26-nucleotide barcode, together with luciferase for in vivo imaging and either GFP or mCherry to facilitate subsequent cell sorting and measurement of reproducibility within a single mouse (Extended Data Fig. 1a, Supplementary Table 1). The 8 barcoded lines were injected as a pool into the left ventricle of 5–6-week-old NOD-SCID-gamma (NSG) mice so as to focus our analysis on the ability of tumour cells to exit circulation and undergo expansion in distant organs. Bioluminescence imaging (BLI) revealed metastatic lesions throughout the body (Extended Data Fig. 1b). Five weeks after injection, brain, lung, liver, kidney and bone were collected, human tumour cells were isolated by fluorescence-activated cell sorting (FACS) using GFP or mCherry, and barcodes were quantified using RNA sequencing (RNA-seq) (Extended Data Fig. 1c–g). Whereas barcode abundances were similar pre-injection, some barcodes were enriched in specific organs (Extended Data Fig. 1h). Different cell lines exhibited distinct patterns of metastatic spread, but each cell line showed highly similar pattern of spread across multiple mice independent of whether GFP or mCherry versions were used, demonstrating the reproducibility of this pooled approach (Extended Data Fig. 1d). For example, HCC1954 was most strongly detected in brain, whereas extracranial metastases were dominated by MDAMB231. Barcodes quantified by bulk RNA-seq were independently validated by quantitative PCR with reverse transcription (RT–qPCR) and single-cell RNA-seq (Extended Data Fig. 1i–m, Supplementary Note 1).

4Having validated the method, we next characterized the metastatic behaviours of all 21 basal-like breast cancer cell lines in the Cancer Cell Line Encyclopedia (CCLE) (Extended Data Fig. 1a–d). Basal-like breast cancers are known to have diverse metastatic abilities in patients. Reflecting this diversity, the cell lines showed disparate metastatic patterns: pan-metastatic, metastatic preferentially to particular organs or not metastatic (Fig. 1b, Supplementary Table 2). Notably, one cell line (BT20) was detected in multiple organs, but at very low abundance in all of them, reflecting its ability to colonize but not expand. To validate the patterns of metastasis observed in the pooled in vivo system, we selected eight cell lines for individual characterization, and observed similar results from the pooled and individual screens (Extended Data Fig. 1n, o).

5Having demonstrated its feasibility in breast cancer, we attempted to expand the mapping of metastatic potential to human cancer cell lines from diverse lineages. To facilitate higher-throughput profiling, we used cell lines barcoded for use with the PRISM method, which was developed for in vitro drug-sensitivity screening. A simplified workflow enabled the quantitative detection of barcodes from crude tissue lysates without the need for FACS-based tumour cell purification (Extended Data Fig. 2, Supplementary Note 2). We applied this method to 503 cell lines spanning 21 lineages to develop a first-generation Metastasis Map (MetMap) (Fig. 2a). The data and interactive visualization are publicly accessible at https://pubs.broadinstitute.org/metmap.

6To test the robustness of the MetMap dataset, we tested cell lines in two formats: in one, we injected all 498 cell lines as a single pool; in the other, we injected 5 pools of 25 lines, with each pool being injected into different mice (referred to as MetMap500 and MetMap125, respectively) (Fig. 2b). We similarly varied cell numbers, mouse age and cohort size to determine whether results varied substantially with these parameters. We observed strong correlation of the metastatic potential despite differences in experimental conditions (Fig. 2c), suggesting that the approach is extremely robust. We also note that intracardiac injection enabled the evaluation of many more cell lines in vivo compared with subcutaneous injection. Specifically, we recovered an average of 197 cell lines per mouse following intracardiac injection, whereas an average of 42 cell lines were recovered following subcutaneous injection (Extended Data Fig. 3a–c). We suspect that this difference is explained by the local competition for nutrients and other microenvironmental factors in the subcutaneous setting, whereas the spatial separation of tumour cells delivered through the intracardiac route minimizes such competition. A similarly reduced diversity was observed in the orthotopic setting, where injection of a pool of nine breast cancer cell lines into the mammary fat pad resulted in a single cell line dominating the resulting tumour (Extended Data Fig. 3d).

7To determine whether the MetMap reflects the metastatic behaviour of human cancers, we analysed available clinical annotations of the cell lines (Fig. 3a–e, Extended Data Fig. 4). We found statistically significant associations with tumour lineage, the site from which the cell line was derived (primary tumour versus metastatic lesion) and patient age. There was no association between metastatic potential and gender or ethnicity. As expected, metastatic potential was higher in certain tumour types, such as melanoma and pancreatic cancer, which also tend to develop metastasis in the human disease setting. By contrast, cell lines derived from brain tumours were generally non-metastatic, reflective of their tendency to not undergo haematogenous spread,. Similarly, the DU145 prostate cancer cell line, derived from a brain metastasis lesion, exhibited brain metastasis in our experimental system. Cell lines derived from metastases showed higher metastatic potential than lines derived from primary tumours, although lines derived from primary tumours known to later give rise to metastases in patients were metastatic in the MetMap (Fig. 3b), consistent with previously reported suggestions that metastatic potential is already encoded in primary tumours–. The association between decreased metastatic potential and increased patient age was unexpected (Fig. 3c), and its basis remains to be determined.

8Perhaps most importantly, extensive variation in metastatic potential was observed within individual lineages, making it possible to search for associations between metastasis propensity and genomic features of the tumours. Of note, metastatic potential was not simply explained by proliferation rate or mutational burden (Fig. 3f–h, Extended Data Fig. 4f, g), suggesting that more subtle molecular determinants of metastasis were involved.

9To develop mechanistic insights, we focused on breast cancer and its potential for brain metastasis (Fig. 1b), because brain metastasis is a feature of some—but not all—breast cancers, and little is known about the underlying factors that could inform therapeutic approaches,. We therefore undertook a systematic and unbiased comparison of the molecular features that distinguished brain metastatic versus non-metastatic lines, using genomic data available for each of the cell lines.

10At the level of somatic mutations, PIK3CA was the top associated correlate: 4 out of 7 brain metastatic lines contained a PIK3CA mutation, compared with 0 out of 14 non-metastatic or weakly metastatic lines (false discovery rate (FDR) = 0.0034) (Fig. 4a, Extended Data Fig. 5a). A fifth line, HCC70, has a loss-of-function mutation in PTEN. PI3K is a principal downstream mediator of ERBB2 (also known as HER2), which itself has been reported to be associated with brain metastasis in humans,. Indeed, two of the brain metastatic cell lines (JIMT1 and HCC1954) also contain typical ERBB2 gene amplifications (Extended Data Fig. 5a).

11At the level of DNA copy number, we observed an association between metastatic potential and deletions of chromosome 8p12–8p21.2 (referred to as 8p) (FDR = 0.0017) (Fig. 4b). Five out of seven brain metastatic breast cancer cell lines contained deletions in this region, compared with 0 out of 14 non-metastatic lines (Extended Data Fig. 5a). A sixth metastatic line, JIMT1, has a small deletion within this commonly deleted region.

12To ascertain the clinical relevance of these associations, we analysed clinical breast cancer datasets for which metastasis information was available. We observed a strong correlation between 8p copy number and gene expression in the METABRIC and TCGA datasets, (Extended Data Fig. 6a), thereby validating 8p expression as a surrogate for copy number in datasets for which copy number data were not available. Coordinated expression of 8p genes stratified tumours into two clusters, with the low-expressing cluster showing enrichment in brain metastasis and lower brain metastasis-free survival (Extended Data Fig. 6b). Whereas 8p loss was more frequent in basal-subtype breast cancer (known to have poor prognosis), 8p loss remained significantly associated with brain metastases within basal tumours. A similar trend was seen in other subtypes, but the sample size was too small to reach statistical significance. Concordant with these findings, the 8p-low signature was strongly enriched in brain metastasis lesions compared with extracranial metastases or primary tumours (Extended Data Fig. 6c, d). Similarly, we observed that response signatures, indicating PI3K activation are associated with brain metastases (Extended Data Fig. 6e–g). The PI3K-high signature tended to co-occur with the 8p-low signature, and the overlapping events captured the majority of patients with brain metastases (Extended Data Fig. 6h, i). These results established the validity of the MetMap experimental system for discovery.

13Confirming these genetic findings, expression analysis revealed enrichment of PI3K and ERBB2 signatures in the brain metastatic cell lines (Fig. 4c). Furthermore, we observed a strong association between brain metastatic potential and a lipid-synthesis signature (Fig. 4c), which has been reported in association with both PI3K activation and 8p-deletion,. To investigate a potential role of lipid metabolism in breast-cancer brain metastatic potential, we analysed the abundance of lipid metabolites across the cell lines. We observed increased levels of cholesterol species in highly brain metastatic cells (Fig. 4d). In addition to cholesterols, membrane lipids including phosphatidylcholine and sphingomyelin were similarly more abundant, as were metabolites associated with the pentose phosphate pathway, which can support cholesterol and lipid synthesis. By contrast, we observed global decreases in levels of triacylglycerols (TAGs) in brain metastatic cells (Fig. 4d). Non-brain metastatic cells had higher levels of TAGs and contained a fatty acid oxidation signature (Fig. 4c). Metabolite profiling of normal mouse tissues showed that the brain has markedly lower levels of TAGs compared with other tissues (Fig. 4e). This reflects brain physiology, whereby instead of storing fatty acids as TAGs, the brain accumulates specialized lipids to support neural activity and brain function. One possibility is that for breast cancer cells to survive in the brain microenvironment, where TAGs and other storage lipids present in other tissues are not abundant, they must access lipids via de novo synthesis or another route, in line with the seed-and-soil hypothesis.

14To further investigate the characteristics of breast cancer cell lines capable of brain metastasis, we analysed genome-wide CRISPR–Cas9 viability-screening data to identify gene vulnerabilities associated with the brain metastatic state. We identified SREBF1 as the top-correlated dependency with brain metastasis (FDR = 0.001) (Fig. 4f). SREBF1 is a pivotal transcription factor that mediates lipid synthesis downstream of the PI3K pathway,. SREBF1 was selectively required for growth of brain metastatic lines in culture compared with breast cancer lines with low or no brain metastatic potential. The association was specific to brain, as no association was observed between SREBF1 essentiality and metastasis to other organs (Fig. 4f). This SREBF1–breast-cancer brain metastasis association was also recovered in the MetMap500 dataset, indicating strong reproducibility of the finding (Extended Data Fig. 5b, c). Of note, the SREBF1 paralogue SREBF2 showed no association between its essentiality in culture and metastatic potential (Fig. 4g).

15To investigate the role of SREBF1 in affecting the lipid phenotype observed in brain metastatic cells, we performed lipidomics after knocking out SREBF1 in JIMT1 and HCC1806 cells using CRISPR–Cas9. SREBF1 knockout resulted in a marked shift in intracellular lipid content, including a decrease in levels of cholesterol, membrane lipids and diacylglycerols (Fig. 4h). SREBF1 knockout also resulted in an increase in intracellular TAG levels, presumably by scavenging TAGs from the lipid-rich serum added to the culture medium. To test this hypothesis, we repeated the experiment in culture medium prepared with delipidated serum, which prevented the increase in TAGs observed in SREBF1-knockout cells (Extended Data Fig. 7).

16To further explore the role of SREBF1, we performed RNA-seq following SREBF1 knockout and found SCD to be the most consistently downregulated gene (Fig. 4i). Consistent with this, SCD was the top co-dependency of SREBF1 across 688 cell lines in the genome-wide CRISPR–Cas9 viability screens (Fig. 4j). The next highest scoring SREBF1 co-dependency was SCAP, which encodes the upstream activator of SREBF1. Comparison of gene expression in breast cancer cells grown in vitro or in the brain similarly showed that in the brain, cells adopted gene-expression signatures of adipogenesis, fatty acid metabolism and xenobiotic metabolism (Extended Data Fig. 8, Supplementary Note 3). The enrichment of lipid-metabolism signatures (including upregulation of SREBF1 and SCD) was unique to brain compared with other sites of metastasis. Similar upregulation was also observed in brain metastases from patients compared with extracranial metastases or their matched primary tumours (Extended Data Fig. 9). Furthermore, the requirement for SREBF1, SCD, SCAP and other members of the lipid-metabolism pathway for brain metastasis formation was confirmed in both mini-pool and individual gene-knockout experiments (Fig. 5a–c, Supplementary Note 4). Together, these genetic, metabolic, transcriptomic and functional genomic evidence all point to an association between SREBF1-mediated lipid metabolism and brain metastasis.

17Given the observation that SREBF1 knockout resulted in a viability defect in vitro (Extended Data Fig. 10a), we compared the relative effect of knockout on metastasis to different organs, to determine whether the viability defect was preferentially observed in brain (Fig. 5d). Five weeks following intracardiac injection of SREBF1-knockout cells, we observed a marked defect in brain metastasis (196-fold reduction), compared with a modest defect in other organs (9–21 fold) (Fig. 5d). Histologic examination of brains from xenografted mice revealed large metastatic lesions in mice receiving wild-type cells, whereas those receiving SREBF1-knockout cells contained micrometastases (Extended Data Fig. 10b), suggesting that SREBF1 is not required for seeding the brain, but rather for proliferation in the brain microenvironment. Consistent with this hypothesis, injection of tumour cells into the carotid artery increased the probability of seeding the brain, but nevertheless a marked growth defect was still observed in SREBF1-knockout cells (Fig. 5e).

18To determine the generality of the SREBF1 requirement for breast cancer growth in the brain, we knocked out SREBF1 in additional brain metastatic lines including HCC1954, MDAMB231 and HCC1806 using CRISPR–Cas9. As with JIMT1, a significant inhibition in brain metastatic growth was also observed in these lines, although the magnitude and duration of growth inhibition varied (Extended Data Fig. 10c, d). The least responsive cell line was HCC1806, in which SREBF1-knockout cells displayed a brain growth defect for the first week, but then assumed a growth trajectory that paralleled wild-type cells. This restoration of growth was not explained by reversion of the genome editing, as brain metastases at the end of the experiment showed evidence of editing at the SREBF1 locus and minimal SREBF1 protein expression (Extended Data Fig. 10e, f). Instead, we found that the SREBF1-independent growth was associated with upregulation of the fatty acid transporter CD36 and the fatty acid-binding protein FABP6 (Extended Data Fig. 10g). Of note, culture of HCC1806 in mouse brain-slice-conditioned medium similarly resulted in upregulation of SCD and CD36 expression (Extended Data Fig. 10h, i). JIMT1 cells did not upregulate CD36 or FABP6 expression following SREBF1 knockout (Extended Data Fig. 10g), perhaps explaining their inability to survive in the brain. Together, these results further demonstrate the relationship between lipid metabolism and brain metastasis, as cells under the selective pressure of SREBF1 loss must acquire lipids by other means to survive in the brain microenvironment.

19This work describes MetMap as an approach for large-scale in vivo characterization of human cancer cell lines. The MetMap resource (available at https://pubs.broadinstitute.org/metmap) currently includes metastasis profiles of 500 cell lines spanning 21 tumour types, providing a large repertoire of models for exploration of metastasis mechanisms. A limitation of the use of human cell lines for such experiments is that they require the use of immunodeficient mice. The extent to which the immune system has a role in mediating patterns of metastasis remains to be determined.

20We followed up only a small proportion of the MetMap findings—specifically, breast cancer metastasis to brain. Multiple lines of experimental and clinical evidence pointed to a role of lipid metabolism in governing the ability of cells to survive in the brain microenvironment. The importance of lipid metabolism in cancer has been highlighted by a number of studies, but its role in brain metastasis has, to our knowledge, not been fully appreciated–. The possibility that interfering with lipid or cholesterol metabolism might abrogate metastatic growth in the brain is particularly intriguing. More generally, this work illustrates the complex interplay between cancer cell growth and the tissue microenvironment.

21No statistical methods were used to predetermine sample size. The experiments were not randomized. The investigators were not blinded to allocation during experiments and outcome assessment.

22Breast cell lines were cultured under the recommended conditions from CCLE (https://portals.broadinstitute.org/ccle). Cell line identities were confirmed by SNP fingerprinting as well as RNA-seq, and compared to the CCLE results. All cell lines were tested negative for mycoplasma. The fluorescence-luciferase-barcode (FLB) construct was engineered using the FUW lentiviral vector backbone (a gift from D. Baltimore; Addgene plasmid no. 14882). Barcodes 26 nucleotides in length were designed using barcode_generator.py (v.2.8; http://comailab.genomecenter.ucdavis.edu/index.php/), and cloned into the landing pad C-terminal to the TGA stop codon of fluorescence luciferase using Gibson assembly (New England Biolabs). Lentivirus preparation and cell infection were performed according to published protocols available at http://www.broadinstitute.org/rnai. Infected cells were analysed by FACS with a fixed gate for GFP or mCherry, using a Sony SH4800 sorter.

23Animal work was performed in accordance with a protocol approved by the Broad Institute Institutional Animal Care and Use Committee (IACUC). NSG female mice (The Jackson Laboratory) at 5–6 weeks of age were used. Cancer cells were suspended in PBS, 0.4% BSA and 100 μl of cell suspensions were injected into the left ventricle of anaesthetized mice (ketamine 100 mg kg−1; xylazine 10 mg kg−1). In vivo metastasis progression was monitored via real-time BLI using the IVIS SpectrumCT Imaging System (PerkinElmer) on a weekly basis. Mice were anaesthetized with inhaling isoflurane, injected intraperitoneally with D-Luciferin (150 mg kg−1), and imaged with the auto exposure setting in prone and supine positions. At the end point, ex vivo BLI was performed by submerging the excised organs in DMEM/F12 medium (Thermo Fisher Scientific) containing D-Luciferin for 10 min and imaged with the auto exposure setting. BLI analysis was performed using Living Image software (v.4.5, PerkinElmer). In the case of breast cancer cohort study (pilot, group 1 and group 2 in Fig. 1, Extended Data Fig. 1), cell lines were mixed at an equal ratio immediately before animal injection, and cell line pools containing 2 × 104 cells per barcoded line were injected. In the case of single breast cell line validation (Extended Data Fig. 1n), cell lines were injected individually at a density of 2 × 104 cells, to be comparable with the pooled experiments. In the case of MetMap125 (Fig. 2, Extended Data Fig. 2), PRISM pools of 25 cell lines were used, and 2.5 × 105 total cells were injected per mouse, corresponding to 1 × 104 cells per barcoded line. Five PRISM pools were injected separately into cohorts of 5–6-week-old NSG mice. In the case of MetMap500, 20 PRISM pools of 25 cell lines were combined to form a large pool of 498 cell lines. The large pool was injected into a cohort of 8–10-week-old NSG mice, with 2.5 × 105 cells per mouse, equivalent to a density of 500 cells per line. Mammary fat pad and subcutaneous injections were performed with Matrigel (Corning) support, at a matching density to their intracardiac assays, respectively (Extended Data Fig. 3). For all pooled cell line experiments, mice were euthanized 5 weeks after injection, in a time-matched manner, unless they displayed severe paralysis or poor body condition, in which case they were euthanized earlier. Intracartoid injection of JIMT1 was performed following a published protocol, at a density of 1 × 105 cells per mouse, similar to the intracardiac injection (Fig. 5e). Intracranial injection was performed as previously described, at a density of 1 × 103 cells per perturbation per animal (Fig. 5a–c, Extended Data Fig. 10c, d).

24Organs including brain, lung, liver, kidney were dissociated using gentleMACS Octo Dissociator with Heaters (Miltenyi Biotec). The optimized dissociation solutions and programs (Miltenyi Biotec) are listed in Supplementary Table 9. Bones (from both hind limbs) were chopped into fine pieces and incubated in the dissociation buffer with vigorous shaking. The dissociated cell suspensions were filtered using 100-μm filters, and washed with DMEM/F12 twice. Cell suspensions were then washed with staining buffer (PBS, 2 mM EDTA, 0.5% BSA), and incubated with mouse cell depletion beads according to the instructions (Miltenyi Biotec). Cell suspensions were subjected to negative selection using autoMACS Pro Separator (Miltenyi Biotec) to deplete mouse stroma. Brains were subjected to an additional myelin-debris-depletion step using myelin removal beads II (Miltenyi Biotec). The resultant cell suspensions were then analysed by FACS using a Sony SH4800 sorter, with the fixed gate for GFP or mCherry. DAPI staining was used to exclude dead cells. For bulk RNA-seq, cells were sorted to a single tube in PBS, 0.4% BSA and RNasin Plus RNase Inhibitor (Promega), centrifuged at 1,500 rpm for 10 min, and cell pellets were frozen at −80 °C for downstream use. For single-cell RNA-seq, single cells were sorted into 96-well plates containing cold TCL buffer (Qiagen) containing 1% β-mercaptoethanol, snap frozen on dry ice, and then stored at −80 °C. Ninety single cells were sorted per plate, the rest wells on the plate were used for negative and positive controls.

25Individual cell lines, cell line pools before injection, and cells isolated from metastases were analysed by RNA-seq. RNA extraction was performed using Quick-RNA MicroPrep according to the manufacturer’s instructions (Zymo Research). RNA was quantified using an RNA 6000 Pico Kit on a 2100 Bioanalyzer (Agilent). RNA samples from cell numbers lower than 500 were not measured but all were used as input for library preparation. cDNA was synthesized using Clontech SmartSeq v.4 reagents from up to 2 ng RNA input according to the manufacturer’s instructions (Clontech). Full-length cDNA was fragmented to a mean size of 150 bp with a Covaris M220 ultrasonicator and Illumina libraries were prepared from 2 ng of sheared cDNA using Rubicon Genomics Thruplex DNaseq reagents according to the manufacturer’s protocol. The finished double-stranded DNA (dsDNA) libraries were quantified by Qubit fluorometer, Agilent TapeStation 2200, and RT–qPCR using the Kapa Biosystems library-quantification kit. Uniquely indexed libraries were pooled in equimolar ratios and sequenced on Illumina NextSeq500 runs with paired-end 75-bp reads at the Dana-Farber Cancer Institute Molecular Biology Core Facilities. RT–qPCR quantification of barcodes was performed using Maxima First Strand cDNA Synthesis Kit, Taqman Fast Advanced Master Mix, custom synthesized Taqman probes, and QuantStudio 6 PCR System (ThermoFisher Scientific). Single-cell RNA-seq was performed as previously described.

26Because the RNA-seq library preparation sheared the cDNA randomly into small pieces, demultiplexed RNA-seq reads were mapped to the barcode references using Bowtie 2 local mode for barcode detection and quantification. Mapped reads were filtered with the criteria that reads (either 5′ or 3′) must cover over 50% of the barcodes from either end, and counted using samtools. Barcode percentage corresponding to cell composition was calculated for single cell lines, pre-injected cell mixtures and in vivo metastasis samples.

27For breast cohort study, metastatic potential of cell line j targeting organ i, Mi,j was calculated as:, in which ci is the total cancer cell number isolated from organ i, pj is the fractional proportion of cell line j estimated by barcode quantification, and n is the number of replicates of mice. To identify features that associate with brain metastatic potential, a two-class comparison method was used. The analysis was performed on mutation, copy number, expression, metabolite, and CRISPR-gene dependency (available at https://depmap.org/portal/). Copy number data were binarized using a cutoff of ≤−1 (loss) and ≥1 (gain).

28Potential mouse contaminating reads were removed by competitive mapping to the human/mouse hybrid genome using BBSplit (https://sourceforge.net/projects/bbmap/). Reads that uniquely mapped to the human genome were then used as input for mapping and gene-level counting with the RSEM package. Gene count estimates were normalized using the TMM method in edgeR. For differential analysis, to properly account for the cancer cell composition differences in each in vivo sample, an in silico modelled in vitro mixture was generated first. For each in silico metastasis model, the estimated expression ĝ of gene i is computed as a weighted average of the cell lines present in the corresponding in vivo sample: in which gi,j is the baseline in vitro expression of gene i in cell line j and pj is the fractional proportion of cell line j in the in vivo sample, as estimated by barcode quantification, and M is the number of cell lines present in the in vivo sample. The in vivo and in silico counterpart were then compared using a paired design for each organ in voom-limma. GSEA was performed using camera or GSEA-preranked method implemented in fgsea,. Single sample GSEA signature projection was performed using gsva package. Gene-signature datasets were from MSigDB (https://www.gsea-msigdb.org/).

29PRISM cell lines (source of each available at https://depmap.org/portal/) were adapted to the same culture conditions in phenol-red-free RPMI1640 medium (ThermoFisher Scientific) and barcoded as previously described. SNP fingerprinting authentication was performed before and after barcoding. Mycoplasma contamination was examined (MycoScope, Genlantis) and only negative lines were used for experiments. These included eight oestrogen-receptor-positive breast cancer cell lines. Despite the lack of phenol red (a weak oestrogen) these breast cell lines maintained ESR1 positivity and expression of a downstream marker of its activity, FOXA1. This is probably explained by the remaining oestrogens in the fetal bovine serum (FBS). PRISM cell lines were pooled on the basis of their in vitro doubling bins, at equal number, in the format of 25 lines per pool, and cryopreserved until use. Cells were thawed and recovered for 48 h before in vivo injection. To form the large pool of 498 cell lines, 20 PRISM pools were mixed at equal total number immediately before injection.

30After in vivo experiments, organs were subjected to tissue dissociation, mouse stroma depletion, and the dissociated cell pellets were frozen at −80 °C as described above. The pellets (≤50 mg dry weight) were lysed in 200 μl freshly prepared lysis buffer with proteinase K, heat digested at 60 °C, and denatured at 95 °C for 10 min. Twenty microlitres of the lysates was used for barcode amplification per 100 μl PCR volume (multiple technical replicates per sample). PCR was performed using the following conditions: 95 °C for 3 min; 98 °C for 20 s, 57 °C for 15 s, 72 °C for 10 s (30 cycles); 72 °C for 5 min; 4 °C stop. PCR libraries were pooled, purified using Select-a-Size DNA Clean & Concentrator Kit (Zymo Research), and quantified using Qubit dsDNA HS Assay Kit (ThermoFisher Scientific) and a 2100 Bioanalyzer (Agilent). The purified 2 nM of libraries with 20% spike-in PhiX DNA were sequenced on Illumina MiSeq or HiSeq at 800 K mm−2 cluster density.

31Demultiplexed sequencing reads were mapped to the barcode reference to generate a table of cell line barcode counts for each sample/condition. Sequencing-depth normalized read counts were used for calculation of relative metastatic potential. Relative metastatic potential of cell line j targeting organ i, rMi,j, was defined as:, in which ci,j is the read counts of cell line j from organ i, pj is the read counts of cell line j from pre-injected population, n is the number of replicate samples of mice, m is the number of replicates of pre-injected population. Confidence intervals were calculated using bootstrap resampling.

32CRISPR–Cas9 versions of cell lines were generated by infecting luciferized cells with Cas9-Blast lentivirus and selecting in 5 μg ml−1 blasticidin for 10 days with continuous passaging until non-infected controls were killed. For pooled in vivo screen, JIMT1–Cas9 cells were infected with a CRISPR guide library (Supplementary Table 10) in an arrayed-fashion in 6-well plates, and selected in 2 μg ml−1 puromycin for 4 days. At this time, non-infected controls were killed, and no growth defect was observed in the perturbed cell lines. Post antibiotic selection, cells were pooled and subjected to intracranial injection at 6 × 104 cells per mouse in 1 μl PBS. This was equivalent to 1 × 103 cells per guide on average per mouse. Intracranial growth was allowed to progress for 4 weeks, and brain tissues were processed adopting the workflow of PRISM in vivo assay, except that guides were amplified using primers targeting the guide vector. Demultiplexed sequencing reads were mapped to the guide reference to generate a table of barcode counts for each guide for each sample. Sequencing-depth was normalized using the upper-quartile method and relative depletion was quantified using a linear model in limma. For individual gene validation (Fig. 5c, Extended Data Fig. 10c, d), Cas9-expressing cells of different cell lines were infected with corresponding guides, selected in 2 μg ml−1 puromycin for 4 days, and subjected to intracranial injection at 1 × 103 cells per mouse in 1 μl PBS. Two independent guides per gene were tested, with one mouse per guide. Intracranial growth was monitored by BLI following injection.

33Positive ion mode analyses of polar and nonpolar lipids (C8-pos) were conducted using a liquid chromatography–mass spectrometry (LC–MS) system composed of a Shimadzu Nexera X2 U-HPLC (Shimadzu) coupled to an Exactive Plus orbitrap mass spectrometer (ThermoFisher Scientific). Cellular extracts were collected from 6-well plate culture, in LC-MS-grade isopropanol (Sigma-Aldrich) containing an internal standard 1,2-didodecanoyl-sn-glycero-3-phosphocholine (Avanti Polar Lipids). Extracts were centrifuged for 10 min at 10,000g to remove residual cellular debris. After centrifugation, supernatants were injected directly onto a 100 × 2.1 mm, 1.7-μm ACQUITY BEH C8 column (Waters). The column was eluted isocratically with 80% mobile phase A (95:5:0.1 v/v/v 10 mM ammonium acetate/methanol/formic acid) for 1 min followed by a linear gradient to 80% mobile phase B (99.9:0.1 v/v methanol/formic acid) over 2 min, a linear gradient to 100% mobile phase B over 7 min, then 3 min at 100% mobile phase B. Mass spectrometry analyses were performed using electrospray ionization in the positive ion mode using full scan analysis over 200 to 1,000 m/z at 70,000 resolution and 3 Hz data acquisition rate. Other mass spectrometry settings were as follows: sheath gas 50, in source collision-induced dissociation 5 eV, sweep gas 5, spray voltage 3 kV, capillary temperature 300 °C, S-lens RF 60, heater temperature 300 °C, microscans 1, automatic gain control target 106, and maximum ion time 100 ms. Lipid identities were determined on the basis of comparison to reference standards and reference plasma extracts and were denoted by the total number of carbons in the lipid acyl chain(s) and total number of double bonds in the lipid acyl chain(s).

34Protein lysates were prepared in RIPA lysis buffer (ThermoFisher Scientific) with cOmplete Mini EDTA-free Protease Inhibitor Cocktail (Roche). Western blot was performed using NuPAGE gel (ThermoFisher Scientific) with wet tank blotting (Bio-Rad) and Odyssey detection system (LI-COR). SREBF1 primary antibody (14088-1-AP, Proteintech), SCD (CD.E10) antibody (ab19862, Abcam), GAPDH (D16H11) XP rabbit monoclonal antibody (5174S, Cell Signaling), β-actin (8H10D10) mouse monoclonal antibody (3700S, Cell Signaling), and IRDye 800CW goat anti-mouse IgG (926-32210, LI-COR), IRDye 680RD goat anti-rabbit IgG (926-68071, LI-COR) secondary antibodies were used. Western blot was performed for cells cultured in different medium conditions. These include RPMI1640 with 10% FBS, with 10% delipidated FBS, with 10% human cerebrospinal fluid (991-19-P-5, Lee BioSolutions), or with 1% SM1 supplement (05711, STEMCELL Tech), or brain-slice-conditioned medium. Brain-slice-conditioned medium was prepared by submerging brain slices (150 μm) in RPMI1640 (no serum) for 48 h. Delipidated FBS was prepared as described.

35METABRIC, TCGA and MSK targeted-sequencing breast cancer datasets were downloaded from cBioPortal. The EMC-MSK dataset including 615 primary tumours (GSE2035, GSE2603, GSE5327 and GSE12276) and the 65-metastasis-sample dataset (GSE14020) were collected and processed as previously described. Paired primary breast tumour and brain metastasis RNA-seq was obtained from ref.. To exclude the confounding effect of brain stroma contamination in this dataset, a contamination indicator generated from GSE52604 was applied, and the contaminating effect was regressed out, generating a corrected gene matrix. PI3K-response signatures were from refs.,. Signature analysis was conducted as described. Hierarchical clustering and heatmaps were generated using gplots package. Other plots were generated using ggplot2. log-rank tests of survival curve difference were calculated using survival package. A multivariate Cox proportional hazards model was built using the coxph function (Extended Data Fig. 6h). Significance of overlap was calculated using chisq.test or fisher.test function.

36Further information on research design is available in the Nature Research Reporting Summary linked to this paper.

37Any methods, additional references, Nature Research reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-020-2969-2.