Effect of STAT3 activity on KRAS-mediated transformation

A. Western blot analysis of STAT3 expression in KP MEFs (Contr) transduced with CRISPR/CAS9 lentivirus expressing CAS9 or gRNAs targeting STAT3 (ST3KO). Two independent knockout clones are shown. ERK1/2 is a loading control.

B. Western blot analysis of total and Y705-phosphorylated STAT3 in KP MEFs transduced with lentivirus expressing control (Vector), wild-type (ST3WT) or mutant STAT3 alleles as noted. ERK1/2 is a loading control.

C. Transduced control (Vector) and KP MEFs were evaluated for foci formation. Cells expressing the indicated STAT3 genotypes were co-cultured with 103 p53KO feeder MEFs and macroscopic colonies were counted after two weeks (n=3 for each cell type, *p<0.05). Values correspond to mean ± s.d.

D. Representative images of tissue culture plates stained with Giemsa to detect foci formed by KP MEFs expressing vector, STAT3 WT or hyperactive STAT3 Y640F.

E. Representative images of transformed foci visualized by fluorescence microscopy. Foci formed by KP MEFs co-expressing GFP with vector alone, STAT3 WT, STAT3 Y705F, or STAT3 Y640F are shown.

STAT3 GOF mutations reduce tumor development in mice

A. Subcutaneous tumor formation in nude mice by KP MEFs expressing vector alone, STAT3 KO, STAT3 K658Y or STAT3 Y640F (104 cells per injection site, n=4 for each cell type). Error bars represent s.d.

B. Western blot analysis of Y705-phosphorylated STAT3 in KP MEFs expressing vector alone, and low or high levels of STAT3 Y640F protein. Cells were fractionated by FACS for GFP. ERK1/2 is a loading control.

C. Subcutaneous tumor formation in nude mice by cells from (B) (104 cells per implant site, n=4 for each cell type). Error bars represent s.d.

D. Quantification of tumor-initiating cell (TIC) frequency in control (Vector) and STAT3 Y640F-expressing KP MEFs by extreme limiting dilution assays (ELDA) in nude mice.

E. Focus formation by KP MEFs expressing vector alone, STAT3 Y640F or double mutants of STAT3 Y640F with DNA binding domain (DBD) mutations, V463Δ and R382W, or phospho-inactivating Y705F. Results represent replicates from two independent experiments (n=4 for each cell type). Values correspond to mean ± s.d..

STAT3 is a genetic modifier of EMT

A. Representative H&E staining of pancreatic tumor sections derived from KPC cells of the indicated KRAS, STAT3, and SMAD4 genotypes.

B. Scatter plot representing differentially expressed genes in parental (STAT3 intact) and STAT3 KO KPC cells.

C. Pie charts representing gene expression profiles of the top 100 upregulated genes in parental control (STAT3 intact) and STAT3 KO KPC cells with percentages linked to specific cellular functions.

D. GO enrichment analysis of the top 100 upregulated genes in STAT3 KO cells compared to STAT3 intact (control) cells.

E. Analysis of overlapping differentially expressed genes (n=243) in STAT3 KO and KRAS KO KPC cell lines relative to parental control cells.

F. Mouse RAS dependency index (RDI) of KPC cells of the indicated KRAS, STAT3, and SMAD4 genotypes.

STAT3 and SMAD4 play opposing roles in pancreatic tumorigenesis

A. Representative H&E staining of pancreatic tumor sections derived from KPC cells of the indicated genotypes, SMAD4 KO (n=12), SMAD4/STAT3 DKO (n=3) and TGFBR2/STAT3 DKO (n=3). At least two independent clones were used. Scale bar = 100μm

B. Western blot analysis of control (Intact), STAT3 KO (ST3KO), STAT3/SMAD4 double KO (ST3/SM4DKO), or STAT3/TGFBR2 double KO (ST3/TRB2DKO) KPC cells. Expression of STAT3, SMAD4, TGFBR2, CDH1 (E-cadherin), and VIM (Vimentin) is shown. ERK1/2 is a loading control.

C. Summary of pancreatic tumor development in nude mice by KPC cells of the indicated STAT3, SMAD4, and TGFBR2 genotypes presented.

D. Pearson correlation heatmap comparing gene expression of SMAD4 KO (SMAD4KO_up, down) and STAT3 KO (STAT3KO_up, down) KPC cells with human PDACs from TCGA (n=168). Signature scores were calculated using the top up- and down-regulated genes in STAT3 and SMAD4 KO KPC cells. The number of genes for each comparison is shown.

E. Pearson correlation heatmap comparing gene expression signature scores of genes regulated in STAT3 KO KPC cells with pancreatic tumors from the COMPASS database (n=92) classified by RAS dependency index (RDI), KRAS-dependent signature (KRAS-sig), KRAS-independent signature (RSK-sig), epithelial (EPI) or mesenchymal (MES) gene signature score.

A. Schematic representation of STAT3 domain structure and the mutants used in this study.

B. Western blot analysis of KP MEFs harboring STAT3 mutations reveals unperturbed PI3K/AKT and MAPK/ERK signaling.

C. Loss of STAT3 does not affect growth and viability of KRASG12D p53KO MEFs. The cells were cultured for two weeks in 2D adherent cultures. Cumulative cell numbers are shown. Values correspond to mean ± s.d.

D. Focus formation assays of p53null MEFs overexpressing STAT3 WT, STAT3 Y640F or STAT3 D661Y.

E. Western blot analyses of total STAT3 and Y705-phosphorylated STAT3 in KP MEFs expressing vector alone, and single and double mutations of STAT3 as noted. ERK1/2 is a loading control.

F. Transactivation by STAT3 mutants of a STAT3-responsive promoter driving a luciferase reporter. Hep3B cells were co-transfected with the p3XGAS-Hsp70-Luc reporter plasmid and mutant STAT3 alleles as noted and stimulated with 20 ng/ml human IL-6 for 24 hours (n=3 for each cell type). Relative luciferase values are presented as mean ± s.d.

G. GO enrichment of differentially expressed genes in KP MEF cells expressing STAT3 Y640F compared to control KP cells provided by Novogene Corp. from RNASeq data.

H. Western blot analysis of SMAD4KO and TGFBR2KO KP MEFs.

I. Results of focus formation assays by control KP MEFs, and MEFs ablated for SMAD4 and TGFBR2 expression with or without expression of STAT3Y640F as indicated (n>6 for each). Values correspond to average and s.d.

A. Western blots of parental KPC pancreatic ductal adenocarcinoma cells, and two isogenic knockout clones for STAT3, SMAD4 or TGFBR2 generated by CRISPR/Cas9 gene editing. ERK1/2 is a loading control.

B. Top up- and down-regulated subset of genes in STAT3 KO KPC cells compared to parental control cells.

C. Heatmaps of differentially expressed genes in KPC parental cells and the indicated KRAS, STAT3, and SMAD4 knockout genotypes. TGF-β super family genes and transcription factors (TFs) are shown.

D. EMT scores of KPC cells of the indicated KRAS, STAT3, and SMAD4 genotypes based on epidermal (EPI) and mesenchymal (MES) genes derived from Serresi et al. Sci. Adv. 2021;7(9) shown to right.

E. GO enrichment analysis of top up-regulated genes in KRAS KO cells compared to parental KPC cells. The normalized FDR (Padj) was provided by Novogene Corp. from RNASeq data.

F. RAS dependency indices of KPC pancreatic whole tumors based on bulk RNA-seq of tumors formed from parental KPC (Contr), KRAS KO or STAT3 KO cells. Murine RAS dependency index (mRDI) was derived from Ischenko et al. Nature Communications, 2021;12(1):1482.

A. Top up- or down-regulated genes in murine KPC STAT3 KO or SMAD4 KO cells derived from RNA-Seq.

B. Signature genes used to designate TCGA and COMPASS PDAC individual tumors as epithelial (EPI) or mesenchymal (MES) derived from T.Z. Tan et al. EMBO Mol Med. 2014;6(10):1279-93.

C. Pearson correlation coefficients of STAT3 KO gene signature scores (top) or SMAD4 KO gene signature scores (bottom) with curated human pancreatic tumor samples of expressing wild-type SMAD4 from TCGA (n=112) classified by Ras dependency index (RDI), KRAS-dependent signature (KRAS-sig), KRAS-independent signature (RSK-sig), epithelial (EPI) or mesenchymal (MES). The number of genes for each comparison is shown.

D. Genes used to define a RAS Dependency Index (RDI) derived from A.Singh et al. Cancer Cell, 2009;15(6):489-500, and genes used to define a KRAS dependent signature (KRAS_sig) or an RSK dependent signature (RSK_sig) derived from T.L.Yuan et al. Cell Rep. 2018;22(7):1889-902.

E. STAT3 Y705 phosphorylation levels in cell lines from the Cancer Cell Line Encyclopedia (n=898).