Multi-Omic Analysis of CIC's Functional Networks Reveals Novel Interaction Partners and a Potential Role in Mitotic Fidelity.

CIC encodes a transcriptional repressor and MAPK signalling effector that is inactivated by loss-of-function mutations in several cancer types, consistent with a role as a tumour suppressor. Here, we used bioinformatic, genomic, and proteomic approaches to investigate CIC's interaction networks. We observed both previously identified and novel candidate interactions between CIC and SWI/SNF complex members, as well as novel interactions between CIC and cell cycle regulators and RNA processing factors. We found that CIC loss is associated with an increased frequency of mitotic defects in human cell lines and an in vivo mouse model and with dysregulated expression of mitotic regulators. We also observed aberrant splicing in CIC-deficient cell lines, predominantly at 3' and 5' untranslated regions of genes, including genes involved in MAPK signalling, DNA repair, and cell cycle regulation. Our study thus characterises the complexity of CIC's functional network and describes the effect of its loss on cell cycle regulation, mitotic integrity, and transcriptional splicing, thereby expanding our understanding of CIC's potential roles in cancer. In addition, our work exemplifies how multi-omic, network-based analyses can be used to uncover novel insights into the interconnected functions of pleiotropic genes/proteins across cellular contexts.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation.

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Discovery of 3-Benzyl-1-( trans-4-((5-cyanopyridin-2-yl)amino)cyclohexyl)-1-arylurea Derivatives as Novel and Selective Cyclin-Dependent Kinase 12 (CDK12) Inhibitors.

Cyclin-dependent kinase 12 (CDK12) plays a key role in the coordination of transcription with elongation and mRNA processing. CDK12 mutations found in tumors and CDK12 inhibition sensitize cancer cells to DNA-damaging reagents and DNA-repair inhibitors. This suggests that CDK12 inhibitors are potential therapeutics for cancer that may cause synthetic lethality. Here, we report the discovery of 3-benzyl-1-( trans-4-((5-cyanopyridin-2-yl)amino)cyclohexyl)-1-arylurea derivatives as novel and selective CDK12 inhibitors. Structure-activity relationship studies of a HTS hit, structure-based drug design, and conformation-oriented design using the Cambridge Structural Database afforded the optimized compound 2, which exhibited not only potent CDK12 (and CDK13) inhibitory activity and excellent selectivity but also good physicochemical properties. Furthermore, 2 inhibited the phosphorylation of Ser2 in the C-terminal domain of RNA polymerase II and induced growth inhibition in SK-BR-3 cells. Therefore, 2 represents an excellent chemical probe for functional studies of CDK12 and could be a promising lead compound for drug discovery.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer.

Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.

A transgenic mouse model demonstrating the oncogenic role of mutations in the polycomb-group gene EZH2 in lymphomagenesis.

The histone methyltransferase EZH2 is frequently mutated in germinal center-derived diffuse large B-cell lymphoma and follicular lymphoma. To further characterize these EZH2 mutations in lymphomagenesis, we generated a mouse line where EZH2(Y641F) is expressed from a lymphocyte-specific promoter. Spleen cells isolated from the transgenic mice displayed a global increase in trimethylated H3K27, but the mice did not show an increased tendency to develop lymphoma. As EZH2 mutations often coincide with other mutations in lymphoma, we combined the expression of EZH2(Y641F) by crossing these transgenic mice with Eµ-Myc transgenic mice. We observed a dramatic acceleration of lymphoma development in this combination model of Myc and EZH2(Y641F). The lymphomas show histologic features of high-grade disease with a shift toward a more mature B-cell phenotype, increased cycling and gene expression, and epigenetic changes involving important pathways in B-cell regulation and function. Furthermore, they initiate disease in secondary recipients. In summary, EZH2(Y641F) can collaborate with Myc to accelerate lymphomagenesis demonstrating a cooperative role of EZH2 mutations in oncogenesis. This murine lymphoma model provides a new tool to study global changes in the epigenome caused by this frequent mutation and a promising model system for testing novel treatments.

Signatures of mutational processes in human cancer.

All cancers are caused by somatic mutations; however, understanding of the biological processes generating these mutations is limited. The catalogue of somatic mutations from a cancer genome bears the signatures of the mutational processes that have been operative. Here we analysed 4,938,362 mutations from 7,042 cancers and extracted more than 20 distinct mutational signatures. Some are present in many cancer types, notably a signature attributed to the APOBEC family of cytidine deaminases, whereas others are confined to a single cancer class. Certain signatures are associated with age of the patient at cancer diagnosis, known mutagenic exposures or defects in DNA maintenance, but many are of cryptic origin. In addition to these genome-wide mutational signatures, hypermutation localized to small genomic regions, 'kataegis', is found in many cancer types. The results reveal the diversity of mutational processes underlying the development of cancer, with potential implications for understanding of cancer aetiology, prevention and therapy.

The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

Mutational processes molding the genomes of 21 breast cancers.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

DriverNet: uncovering the impact of somatic driver mutations on transcriptional networks in cancer.

Simultaneous interrogation of tumor genomes and transcriptomes is underway in unprecedented global efforts. Yet, despite the essential need to separate driver mutations modulating gene expression networks from transcriptionally inert passenger mutations, robust computational methods to ascertain the impact of individual mutations on transcriptional networks are underdeveloped. We introduce a novel computational framework, DriverNet, to identify likely driver mutations by virtue of their effect on mRNA expression networks. Application to four cancer datasets reveals the prevalence of rare candidate driver mutations associated with disrupted transcriptional networks and a simultaneous modulation of oncogenic and metabolic networks, induced by copy number co-modification of adjacent oncogenic and metabolic drivers. DriverNet is available on Bioconductor or at http://compbio.bccrc.ca/software/drivernet/.

Recurrent somatic DICER1 mutations in nonepithelial ovarian cancers.

BACKGROUND: Germline truncating mutations in DICER1, an endoribonuclease in the RNase III family that is essential for processing microRNAs, have been observed in families with the pleuropulmonary blastoma-family tumor and dysplasia syndrome. Mutation carriers are at risk for nonepithelial ovarian tumors, notably sex cord-stromal tumors. METHODS: We sequenced the whole transcriptomes or exomes of 14 nonepithelial ovarian tumors and noted closely clustered mutations in the region of DICER1 encoding the RNase IIIb domain of DICER1 in four samples. We then sequenced this region of DICER1 in additional ovarian tumors and in certain other tumors and queried the effect of the mutations on the enzymatic activity of DICER1 using in vitro RNA cleavage assays. RESULTS: DICER1 mutations in the RNase IIIb domain were found in 30 of 102 nonepithelial ovarian tumors (29%), predominantly in Sertoli-Leydig cell tumors (26 of 43, or 60%), including 4 tumors with additional germline DICER1 mutations. These mutations were restricted to codons encoding metal-binding sites within the RNase IIIb catalytic centers, which are critical for microRNA interaction and cleavage, and were somatic in all 16 samples in which germline DNA was available for testing. We also detected mutations in 1 of 14 nonseminomatous testicular germ-cell tumors, in 2 of 5 embryonal rhabdomyosarcomas, and in 1 of 266 epithelial ovarian and endometrial carcinomas. The mutant DICER1 proteins had reduced RNase IIIb activity but retained RNase IIIa activity. CONCLUSIONS: Somatic missense mutations affecting the RNase IIIb domain of DICER1 are common in nonepithelial ovarian tumors. These mutations do not obliterate DICER1 function but alter it in specific cell types, a novel mechanism through which perturbation of microRNA processing may be oncogenic. (Funded by the Terry Fox Research Institute and others.).

Raman microscopy-based cytochemical investigations of potential niche-forming inhomogeneities present in human embryonic stem cell colonies.

Measuring spatial and temporal patterns of cytochemical variation in human embryonic stem cell (hESC) colonies is necessary for understanding the role of cellular communication in spontaneous differentiation, the mechanisms of biological niche creation, and structure-generating developmental processes. Such insights will ultimately facilitate directed differentiation and therewith promote advances in tissue engineering and regenerative medicine. However, the patterns of cytochemical inhomogeneities of hESC colonies are not well studied and their causes are not fully understood. We used Raman spectroscopic mapping to contrast supracellular variations in cytochemical composition across pluripotent and partly differentiated hESC colonies to gain a better understanding of the early-stage (i.e., 5 days) effects of the differentiation process on the nature and evolution of these patterns. Higher protein-to-nucleic acid ratios, a differentiation status indicator observed previously using Raman spectroscopy, confirmed reported results that spontaneous differentiation is more pronounced on the edges of a colony than elsewhere. In addition, pluripotent and partly differentiated colonies also showed higher lipid concentrations relative to nucleic acids at colony edges in contrast to relative glycogen concentrations, which were up to 400% more pronounced in the colony centers compared to their edges. Pluripotent and partly differentiated colonies differed, with the latter having higher average protein-to-nucleic acid and lipid-to-nucleic acid ratios but a lower glycogen-to-nucleic acid ratio. In both cases, cell density, pluripotency, and high glycogen appeared to vary in tandem. Spatial variations in glycogen- and protein-to-nucleic acid ratios have features on the order of 100 µm and larger. These dimensions are consistent with those reported for stem cell niches and suggest that cytochemical inhomogeneities may provide colony-level information about niches and niche formation. These results demonstrate Raman mapping to be a potentially useful technique for revealing the complexities in the spatial organization of hESC cultures and thus for monitoring the evolution of engineered hESC niches.

Absolute quantification of intracellular glycogen content in human embryonic stem cells with Raman microspectroscopy.

We present a method to perform absolute quantification of glycogen in human embryonic stem cells (hESCs) in situ based on the use of Raman microspectroscopy. The proposed quantification method was validated by comparison to a commonly used commercial glycogen assay kit. With Raman microspectroscopy, we could obtain the glycogen content of hESCs faster and apparently more accurately than with the kit. In addition, glycogen distributions across a colony could be obtained. Raman spectroscopy can provide reliable estimates of the in situ glycogen content in hESCs, and this approach should also be extensible to their other biochemical constituents as well as to other cell types.

Comprehensive analysis of mammalian miRNA* species and their role in myeloid cells.

Processing of pre-miRNA through Dicer1 generates an miRNA duplex that consists of an miRNA and miRNA* strand. Despite the general view that miRNA*s have no functional role, we further investigated miRNA* species in 10 deep-sequencing libraries from mouse and human tissue. Comparisons of miRNA/miRNA* ratios across the miRNA sequence libraries revealed that 50% of the investigated miRNA duplexes exhibited a highly dominant strand. Conversely, 10% of miRNA duplexes showed a comparable expression of both strands, whereas the remaining 40% exhibited variable ratios across the examined libraries, as exemplified by miR-223/miR-223* in murine and human cell lines. Functional analyses revealed a regulatory role for miR-223* in myeloid progenitor cells, which implies an active role for both arms of the miR-223 duplex. This was further underscored by the demonstration that miR-223 and miR-223* targeted the insulin-like growth factor 1 receptor/phosphatidylinositol 3-kinase axis and that high miR-223* levels were associated with increased overall survival in patients with acute myeloid leukemia. Thus, we found a supporting role for miR-223* in differentiating myeloid cells in normal and leukemic cell states. The fact that the miR-223 duplex acts through both arms extends the complexity of miRNA-directed gene regulation of this myeloid key miRNA.

Retinoblastoma-binding proteins 4 and 9 are important for human pluripotent stem cell maintenance.

OBJECTIVE: The molecular mechanisms that maintain human pluripotent stem (PS) cells are not completely understood. Here we sought to identify new candidate PS cell regulators to facilitate future improvements in their generation, expansion, and differentiation. MATERIALS AND METHODS: We used bioinformatic analyses of multiple serial-analysis-of-gene-expression libraries (generated from human PS cells and their differentiated derivatives), together with small interfering RNA (siRNA) screening to identify candidate pluripotency regulators. Validation of candidate regulators involved promoter analyses, Affymetrix profiling, real-time PCR, and immunoprecipitation. RESULTS: Promoter analysis of genes differentially expressed across multiple serial-analysis-of-gene-expression libraries identified E2F motifs in the promoters of many PS cell-specific genes (e.g., POU5F1, NANOG, SOX2, FOXD3). siRNA analyses identified two retinoblastoma binding proteins (RBBP4, RBBP9) as required for maintenance of multiple human PS cell types. Both RBBPs were bound to RB in human PS cells, and E2F motifs were present in the promoters of genes whose expression was altered by decreasing RBBP4 and RBBP9 expression. Affymetrix and real-time PCR studies of siRNA-treated human PS cells showed that reduced RBBP4 or RBBP9 expression concomitantly decreased expression of POU5F1, NANOG, SOX2, and/or FOXD3 plus certain cell cycle genes (e.g., CCNA2, CCNB1), while increasing expression of genes involved in organogenesis (particularly neurogenesis). CONCLUSIONS: These results reveal new candidate positive regulators of human PS cells, providing evidence of their ability to regulate expression of pluripotency, cell cycle, and differentiation genes in human PS cells. These data provide valuable new leads for further elucidating mechanisms of human pluripotency.

The testosterone-dependent and independent transcriptional networks in the hypothalamus of Gpr54 and Kiss1 knockout male mice are not fully equivalent.

BACKGROUND: Humans and mice with loss of function mutations in GPR54 (KISS1R) or kisspeptin do not progress through puberty, caused by a failure to release GnRH. The transcriptional networks regulated by these proteins in the hypothalamus have yet to be explored by genome-wide methods. RESULTS: We show here, using 1 million exon mouse arrays (Exon 1.0 Affymetrix) and quantitative polymerase chain reaction (QPCR) validation to analyse microdissected hypothalamic tissue from Gpr54 and Kiss1 knockout mice, the extent of transcriptional regulation in the hypothalamus. The sensitivity to detect important transcript differences in microdissected RNA was confirmed by the observation of counter-regulation of Kiss1 expression in Gpr54 knockouts and confirmed by immunohistochemistry (IHC). Since Gpr54 and Kiss1 knockout animals are effectively pre-pubertal with low testosterone (T) levels, we also determined which of the validated transcripts were T-responsive and which varied according to genotype alone. We observed four types of transcriptional regulation (i) genotype only dependent regulation, (ii) T only dependent regulation, (iii) genotype and T-dependent regulation with interaction between these variables, (iv) genotype and T-dependent regulation with no interaction between these variables. The results implicate for the first time several transcription factors (e.g. Npas4, Esr2), proteases (Klk1b22), and the orphan 10-transmembrane transporter TMEM144 in the biology of GPR54/kisspeptin function in the hypothalamus. We show for the neuronal activity regulated transcription factor NPAS4, that distinct protein over-expression is seen in the hypothalamus and hippocampus in Gpr54 knockout mice. This links for the first time the hypothalamic-gonadal axis with this important regulator of inhibitory synapse formation. Similarly we confirm TMEM144 up-regulation in the hypothalamus by RNA in situ hybridization and western blot. CONCLUSIONS: Taken together, global transcriptional profiling shows that loss of GPR54 and kisspeptin are not fully equivalent in the mouse hypothalamus.

P-cadherin expression as a prognostic biomarker in a 3992 case tissue microarray series of breast cancer.

P-cadherin is a calcium-dependent cell-cell adhesion glycoprotein. P-cadherin expression is restricted to the myoepithelial cells in normal breast tissue, and aberrant staining has also been described in invasive tumors. Several small studies have reported P-cadherin as a marker of poor outcome in breast cancer patients but its prognostic significance in relation to other variables has not been established in a large series of breast cancers. A tissue microarray was constructed from 3992 cases of invasive breast carcinoma, and P-cadherin expression was evaluated using immunohistochemistry. Median follow-up was 12.5 years. The immunohistochemistry-based definitions of cancer subtypes were luminal (ER+ or PR+/HER2-), luminal/HER2+ (ER+ or PR+/HER2+), HER2+ (ER-/PR-/HER2+), and basal (ER-/PR-/HER2-/CK5/6+ or EGFR+). Clinical covariate and biomarker associations were assessed using contingency tables, and Pearson's χ(2) or Fisher's exact test. Survival associations were assessed using Kaplan-Meier plots, logrank and Breslow tests, and Cox proportional hazards regression analysis. P-cadherin was expressed in 34.8% (1290/3710, 50% cut point) of cases. P-cadherin staining was strongly associated with HER2+ and basal carcinoma subtypes (P<0.0005). P-cadherin-positive patients showed significantly poorer short-term (0-10 years) overall survival, disease-specific survival, distant relapse-free interval, and locoregional relapse-free interval in univariable models (P<0.05). In multivariable Cox models containing standard clinical covariates and cancer subtypes, P-cadherin did not show independent prognostic value. P-cadherin expression was positively associated with histological grade, chemotherapy, Ki-67, EGFR, CK5/6, p53, YB-1, and HER2 expression (P<0.002), and negatively associated with age at diagnosis, ER, PR, and Bcl-2 expression (P<0.0005). This study shows the value of P-cadherin as a marker of poor prognosis. The large sample size of this series clarifies contradictory findings of many smaller studies. P-cadherin positivity is associated with high-grade tumor subtypes and well-established markers of poor prognosis, and may represent a promising antibody therapeutic target.

PPM1H is a p27 phosphatase implicated in trastuzumab resistance.

The HER2 oncogene is overexpressed or amplified in 20% of breast cancers. HER2-positive cancer historically portends a poor prognosis, but the HER2-targeted therapy trastuzumab mitigates this otherwise ominous distinction. Nevertheless, some patients suffer disease recurrence despite trastuzumab, and metastatic disease remains largely incurable due to innate and acquired resistance. Thus, understanding trastuzumab resistance remains an unmet medical need. Through RNA interference screening, we discovered that knockdown of the serine/threonine phosphatase PPM1H confers trastuzumab resistance via reduction in protein levels of the tumor suppressor p27. PPM1H dephosphorylates p27 at threonine 187, thus removing a signal for proteasomal degradation. We further determined that patients whose tumors express low levels of PPM1H trend towards worse clinical outcome on trastuzumab. Identifying PPM1H as a novel p27 phosphatase reveals new insight into how cancer cells destabilize a well-recognized tumor suppressor. Furthermore, low PPM1H expression may identify a subset of HER2-positive tumors that are harder to treat.

Somatic mutations at EZH2 Y641 act dominantly through a mechanism of selectively altered PRC2 catalytic activity, to increase H3K27 trimethylation.

Next-generation sequencing of follicular lymphoma and diffuse-large B-cell lymphoma has revealed frequent somatic, heterozygous Y641 mutations in the histone methyltransferase EZH2. Heterozygosity and the presence of equal quantities of both mutant and wild-type mRNA and expressed protein suggest a dominant mode of action. Surprisingly, B-cell lymphoma cell lines and lymphoma samples harboring heterozygous EZH2(Y641) mutations have increased levels of histone H3 Lys-27-specific trimethylation (H3K27me3). Expression of EZH2(Y641F/N) mutants in cells with EZH2(WT) resulted in an increase of H3K27me3 levels in vivo. Structural modeling of EZH2(Y641) mutants suggests a "Tyr/Phe switch" model whereby structurally neutral, nontyrosine residues at position 641 would decrease affinity for unmethylated and monomethylated H3K27 substrates and potentially favor trimethylation. We demonstrate, using in vitro enzyme assays of reconstituted PRC2 complexes, that Y641 mutations result in a decrease in monomethylation and an increase in trimethylation activity of the enzyme relative to the wild-type enzyme. This represents the first example of a disease-associated gain-of-function mutation in a histone methyltransferase, whereby somatic EZH2 Y641 mutations in lymphoma act dominantly to increase, rather than decrease, histone methylation. The dominant mode of action suggests that allele-specific EZH2 inhibitors should be a future therapeutic strategy for this disease.

The orphan adhesion-GPCR GPR126 is required for embryonic development in the mouse.

Adhesion-GPCRs provide essential cell-cell and cell-matrix interactions in development, and have been implicated in inherited human diseases like Usher Syndrome and bilateral frontoparietal polymicrogyria. They are the second largest subfamily of seven-transmembrane spanning proteins in vertebrates, but the function of most of these receptors is still not understood. The orphan Adhesion-GPCR GPR126 has recently been shown to play an essential role in the myelination of peripheral nerves in zebrafish. In parallel, whole-genome association studies have implicated variation at the GPR126 locus as a determinant of body height in the human population. The physiological function of GPR126 in mammals is still unknown. We describe a targeted mutation of GPR126 in the mouse, and show that GPR126 is required for embryonic viability and cardiovascular development.

Synthetic lethality screens reveal RPS6 and MST1R as modifiers of insulin-like growth factor-1 receptor inhibitor activity in childhood sarcomas.

The insulin-like growth factor-1 receptor (IGF1R) is emerging as a promising therapeutic target in human cancers. In the high-risk childhood sarcomas Ewing family tumor and rhabdomyosarcoma, IGF1R-blocking antibodies show impressive antitumor activity in some but not all patients, and acquired resistance is observed. Because tumor IGF1R mutations are not described, the basis of IGF1R inhibitor resistance remains unknown. We hypothesized that compensatory signaling cascades bypassing targeted IGF1R inhibition might be involved. To test this systematically, we performed small interfering RNA (siRNA) screens in sarcoma cell lines to identify IGF1R pathway components or related protein tyrosine kinase (PTK) networks that modulate the antitumor efficacy of the BMS-536924 IGF1R kinase inhibitor. This strategy revealed (a) that sarcoma cells are exquisitely sensitive to loss of distal rather than proximal IGF1R signaling components, such as ribosomal protein S6 (RPS6); (b) that BMS-536924 fails to block RPS6 activation in resistant sarcoma cell lines; and (c) that siRNA knockdown of the macrophage-stimulating 1 receptor tyrosine kinase (MST1R; also known as RON) restores BMS-536924 efficacy, even in highly drug-resistant cell lines. We confirmed MST1R expression across a broad panel of childhood sarcomas, and found that loss of MST1R by RNA interference blocks downstream RPS6 activation when combined with BMS-536924 in vitro. These findings underscore the importance of fully understanding PTK networks for successful clinical implementation of kinase inhibitor strategies.