Supplementary Figure 1. Lkb1-deficient lung ADC progressively transdifferentiates into SCC. (a) A scheme showing the progression pattern of atypical adenomatous hyperplasia/epithelial hyperplasia (AAH/EH), adenoma, ADC and SCC in Kras/Lkb1mice post Ad-Cre treatment. (b) H&E staining and p63 immunohistochemical staining on lung tissues from Kras/Lkb1 1
mice at indicated time points post Ad-Cre treatment. Scale bar: 150 μm. (c) Statistic analysis of p63 positive lesions in the lungs from Kras/Lkb1 mice at 2, 4, 6, 8, 10 weeks post Ad-Cre treatment (n=8 each group). Data were shown as mean ± s.e.m. t test, **P < 0.01. (d) A scheme showing the strategy of serial transplantation of mouse lung ADC with Lkb1 deficiency in nude mice. (e) Representative histology, p63, Ttf1 immunohistochemical staining and collagen staining for tumor xenografts at different passages of serial transplantation in nude mice. Scale bar: 50 μm. (f) Statistic analysis of p63 positive cells on lung tumor sections at passage 1, 2 and 3 from the serial transplatation experiment in nude mice. 10 serial sections were statistically analyzed in each passage. Data were shown as mean ± s.e.m. t test, ***P < 0.001. (g) Representative histology and immunohistochemical staining of p63 and SP-C on mad-sccs from Kras/Lkb1 mice at 8 and 12 weeks post Ad-Cre treatment. Scale bar: 50 μm. (h) Statistic analysis of SP-C + /p63 + cell ratio in the mad-sccs from Kras/Lkb1 mice at 8 (n=10) and 12 weeks (n=10) post Ad-Cre treatment. 2
Supplementary Figure 2. Type II pneumocyte lineage-derived lung ADC with Lkb1 deficiency transdifferentiates into SCC. (a) Representative immunofluorescent photoes for LacZ (Green) and Sox2 (Red) staining in cyro section of SP-C-CreERT2/Rosa26R-LacZ mouse at 2 weeks post tamoxifen administration. Scale bar: 150 μm. (b) Statistic analysis of the percentage of LacZ positive cells in both type II pneumocytes and basal cells. (c) Schematic model (upper panel) and representative histologies (lower panel) for AAH, adenoma, ADC and SCC from SKL mice at a serial of time points post 3
tamoxifen administration. Scale bar: 150 μm. (d) Immunohistochemical staining of SP-C on lung tissues containing AAH, adenoma and adenocarcinoma from SKL mice at 4-8 weeks post tamoxifen administration. Scale bar: 150 μm. (e) Statistic analysis of the ADC or SCC incidence in SKL mice at 7 weeks post tamoxifen administration. (f) Representative immunostaining of SP-C, p63 and X-gal staining in ADC and SCC from SRKL mouse at 9 weeks post tamoxifen administration. Scale bar: 150 μm. (g) Representative immunohistochemical staining for p63 and SP-C in lung mad-scc from SKL mice at 9 weeks post tamoxifen administration. Scale bar: 150 μm. (h) Quantification of tumor percentage for mixed mad-scc and typical SCC in SKL mice at indicated time points post tamoxifen administration. n=8 for each time point. (i) Representative histology, p63, SP-C immunohistochemical staining on mad-scc lesions in Kras/Lkb1 mice with 8 and 12 weeks Ad-Cre treatment. Scale bar: 50 μm. (j) Statistic analysis of SP-C+/p63+ cell ratio in the mice with 9-12 (n=8) and 13-16 (n=10) weeks Ad-Cre treatment.. 4
Supplementary Figure 3. Validation of a serial of genes differentially expressed in ADC and SCC from Kras/Lkb1 mouse model. Real-time PCR validation of genes differentially expressed in ADC (n=5) and SCC (n=7) from Kras/Lkb1 mouse model according to microarray data analyses. Data were shown as mean ± s.e.m. t test, *P <0.05, **P < 0.01,***P<0.001. 5
Supplementary Figure 4. Detection of Lox expression and collagen deposition in human lung ADC and SCC samples. (a) Real-time analysis of Lox expression in human lung ADC (n=85) and SCC (n=85) samples. Data were shown as mean ± s.e.m. t test, *P<0.05. (b) Representative photos showing immunohistochemical staining for Lox in human lung ADC and SCC samples. (c) Statistic analysis of Lox expression in human lung ADC and SCC samples. Pearson χ2 test. (d) Representative photos showing Massons Trichrome staining of collagen in human lung ADC and SCC samples. (e) Statistic analysis of collagen deposition in human lung ADC and SCC samples. Pearson χ2 test. 6
Supplementary Figure 5. Lentivirus-mediated Gene Expression in Kras/Lkb1 Mice. (a) A scheme of lentivirus-mediated gene expression (Lenti-Cre and Lenti-Cre-MC) in Kras/Lkb1 mouse model. (b) Representative H&E staining, p63 and SP-C immunohistochemical staining of ADC and SCC from Kras/Lkb1 mice virally infected with Lenti-Cre. Scale bar: 150 μm. 7
Supplementary Figure 6. Lentivirus-mediated Lox expression in Kras/Lkb1 mice. Representative H&E staining as well as SP-C immunohistochemical staining in lung ADC from Kras/Lkb1 mice virally infected with Lenti-Cre or Lenti-Cre-Lox. Scale bar: 10 μm. 8
Supplementary Figure 7. Pharmacological inhibition of LOX enzymatic activity promotes ADC to SCC transdifferentiation. (a) Detection of Lox enzymatic activity in sera from Kras/Lkb1 mice treated with BAPN (n=12), DPA (n=10) or Saline (n=14). Data were shown as mean ± s.e.m. t test ***P < 0.001. (b) Detection of collagen deposition by Massons Trichrome staining in lung tumors from Kras/Lkb1 9
mice treated with BAPN, DPA or saline. Scale bar: 10 μm. (c) Representative histology and immunohistochemical staining of p63 and SP-C on mad-sccs from Kras/Lkb1 mice with or without 4 weeks BAPN treatment after 4 weeks Ad-Cre treatment. Scale bar: 50 μm. (d) Statistic analysis of SP-C + /p63 + cell ratio in the mad-sccs from Kras/Lkb1 mice with or without 4 weeks BAPN treatment after 4 weeks Ad-Cre treatment (n=12 each group). (e) A scheme showing the BAPN treatment strategy in SKL mouse model. (f-g) Quantification of average tumor number per mouse for ADC and SCC (F) and mad-scc (G) from SKL mice treated with BAPN (n=14) or saline (n=18). Data were shown as mean ±s.e.m. t test, *P <0.05, **P < 0.01. (h-i) Quantification of average tumor size per mouse for ADC and SCC (h) and individual tumor size for SCC (i) from SKL mice post BAPN (n=12) or saline (n=12) treatment. Scale bar: 150 μm. Data were shown as mean ± s.e.m. t test, *P <0.05. (j-k) Quantification of average tumor number per mouse for ADC and SCC (j) and mad-scc (k) from Kras/Lkb1 mice at 4 weeks post DPA (n=14) or saline (n=16) treatment. Data were shown as mean ± s.e.m. t test, *P <0.05, ***P < 0.001. (l) Quantification of average tumor size per mouse for ADC and SCC from Kras/Lkb1 mice at 4 weeks post DPA (n=14) or saline (n=16) treatment. Data were shown as mean ± s.e.m. (m) Quantification of individual tumor size of SCC in Kras/Lkb1 mice post DPA (n=10) or saline (n=14) treatment. Data were shown as mean ± s.e.m. t test. (n-o) Representative Ki67 immunohistochemical staining (n) and statistical analysis (o) in lung ADC and SCC from Kras/Lkb1 mice post BAPN (n=12), DPA (n=10) or saline (n=14) treatment. More than 200 high-power fields (HPF) were analyzed each mouse. Scale bar: 150 μm. t test, ***P < 0.001. (p-q) Representative Cleaved caspase 3 immunohistochemical staining (p) and statistical analysis (q) in lung ADC and SCC from Kras/Lkb1 mice post BAPN (n=12), DPA (n=10) or saline (n=14) treatment. More than 200 high-power fields (HPF) were analyzed each mouse. Scale bar: 150 μm. Data were shown as mean ± s.e.m, t test, ***P < 0.001. 10
Supplementary Figure 8. Activation of DNp63α partially promotes ADC to SCC transdifferentiation. (a) Lung ADC (n=5) and SCC (n=7) dissected from Kras/Lkb1mice at 10 weeks post Ad-Cre treatment were subjected to RNA preparation and real time PCR analyses. Data were shown as mean ± s.e.m.t test. *P <0.05.(b) Western blotting confirmed the ectopic expression of DNp63α in Kras/p53cells. (c) Immunofluorescence staining showed that transient DNp63αexpression partially activates Krt5 and Krt14 expression in Kras/p53cells. Scale bar: 150μm.(d) Distinct morphologies of Kras/p53cells with stablednp63αexpression in contrast to the parental cells. Scale bar: 75μm.
Supplementary Figure 9. Full gel scans for mouse Actin, SPC, Lox and p63(dnp63α) blots are shown. 12
Supplementary Table 1. The 20 most significantly deregulated signal pathways between ADC and SCC from Kras/Lkb1 mice. Fisher s Exact Test was used to generate the p value. Pathway SCC1 SCC2 SCC3 SCC4 SCC5 ADC1 ADC2 ADC3 ADC4 ADC5 Accuracy Validated transcriptional targets of deltanp63 1.03E-07 1.83E-05 4.95E-07 4.00E-03 2.27E-03-1.86E-05-1.63E-03-1.07E-06-3.18E-05-7.99E-06 100.00% isoforms Validated transcriptional targets of TAp63 3.47E-03 1.04E-02 3.15E-03 2.90E-02 3.39E-05-7.43E-06-3.97E-04-1.70E-04-5.05E-04-9.78E-03 100.00% isoforms Tight junction interactions 2.07E-01 3.91E-02 8.74E-03 2.62E-03 1.16E-01-2.81E-03-1.74E-02-1.22E-02-2.10E-02-8.22E-02 100.00% Cell junction organization 2.92E-02 5.85E-02 7.80E-03 7.58E-05 5.49E-03-4.08E-04-8.28E-04-1.29E-02-9.32E-03-4.32E-02 100.00% Direct p53 effectors 3.61E-03 1.64E-03 5.85E-05-1.48E-01-2.02E-01-4.08E-03-1.03E-02-1.95E-06-2.14E-04-6.88E-06 100.00% p53 signaling pathway 1.32E-04 8.11E-03 1.14E-04 1.91E-03 1.57E-05-2.12E-03-9.20E-03-1.08E-06-9.03E-07-1.13E-04 100.00% Apoptosis 1.86E-02 3.47E-02 4.98E-02 4.23E-03 5.69E-02-1.40E-01-6.99E-01-2.54E-03-2.15E-02-6.73E-02 100.00% Muscle cell TarBase 1.11E-02 7.38E-03 1.29E-01-1.05E-01 1.54E-02-2.81E-01-6.39E-02-1.01E-03-8.61E-04-1.23E-02 100.00% Interferon alpha/beta signaling 1.33E-02 4.69E-07 2.47E-03 1.18E-03 3.52E-03-2.38E-01-3.46E-05-7.52E-04-3.33E-02-8.31E-03 100.00% Apoptosis Modulation by HSP70 1.23E-01 1.83E-03 2.05E-02 3.23E-02 2.13E-01-5.05E-03-5.97E-02-3.00E-05-5.65E-06-4.57E-02 100.00% Nucleotide Metabolism 6.99E-02 5.77E-02 6.77E-04 1.31E-02 1.26E-01-4.14E-01 1.00E+00-1.06E-02-1.64E-02-1.72E-03 100.00% Staphylococcus aureus infection -3.46E-02-1.08E-01-4.04E-04-1.60E-07-1.86E-06 1.00E-07 1.08E-07 3.86E-05 8.66E-03 4.10E-06 100.00% NCAM1 interactions 3.95E-01-1.13E-01-1.38E-01-9.70E-02-8.19E-04 1.14E-05 5.69E-04 1.65E-03 2.21E-03 5.23E-06 100.00% Protein digestion and absorption -1.26E-01-2.76E-03-9.64E-02-5.40E-01-1.90E-01 2.38E-05 3.08E-03 2.78E-02 1.21E-02 4.44E-06 100.00% NCAM signaling for neurite out-growth 1.00E+00-2.25E-01-2.69E-01-3.50E-01-1.03E-02 1.95E-04 7.50E-03 1.89E-02 2.43E-02 1.05E-05 100.00% Adipogenesis -2.87E-01-9.52E-02-7.73E-04-4.26E-02-2.57E-03 1.58E-02 1.31E-02 3.65E-07 1.03E-05 1.90E-02 100.00% Signaling by PDGF 1.00E+00-6.74E-01-2.68E-01-3.41E-01-7.88E-02 1.66E-03 7.16E-04 1.55E-01 3.70E-02 6.85E-05 100.00% Integrin cell surface interactions 7.40E-01-9.10E-02-1.20E-01-4.93E-03-2.54E-02 8.93E-07 1.89E-02 1.60E-02 2.11E-02 6.62E-05 100.00% Initial triggering of complement -1.20E-03-1.60E-01-1.79E-01-1.38E-02-1.01E-02 3.81E-03 8.51E-03 4.18E-01 1.00E+00 8.79E-05 100.00% Beta3 integrin cell surface interactions 1.61E-01 3.56E-01-4.35E-03-4.26E-01-2.11E-02 4.61E-04 1.42E-05 3.33E-02 4.02E-02 8.98E-07 100.00% 13
Supplememtary Table 2. Genes significantly deregulated in SCC relative to ADC from Kras/Lkb1 mice. Fisher s Exact Test was used to generate the p value. Description Gene T score P value Lox -5.12922 0.002159 Lox family Loxl1-3.20176 0.018559 Loxl2-2.74239 0.033629 Loxl3-3.36864 0.015068 Col4a1-3.29365 0.016538 Col4a2-3.67688 0.010369 Col4a3-2.11822 0.078476 Collagen members Col4a5-4.32509 0.004955 Col4a6-2.15132 0.07497 Col14a1-2.44139 0.050376 Col3a1-2.62693 0.03922 Col6a1-4.32782 0.00494 Col6a2-5.42459 0.001626 Timp2-2.63373 0.038865 ECM assembly Timp3-4.90136 0.002707 Timp4-2.96583 0.025093 Vim -1.09753 0.31449 Snai1-3.09847 0.021158 EMT Fn1-2.19785 0.07031 Cdh1 1.781081 0.125189 Cdh3 3.41668 0.014202 Trp63 15.61573 4.37E-06 Sox2 6.19363 0.000816 Squamous markers Krt5 56.22009 2.13E-09 Krt6a 31.89217 6.32E-08 Krt14 8.250988 0.000171 ECM degradation Mmp9 2.23276 0.067011 Mmp14 2.451192 0.04971 Perp 9.258257 8.97E-05 Pkp1 10.77405 3.78E-05 Pkp3 7.439305 0.000304 Desmosome Jup 8.517886 0.000143 Dsc2 9.086583 9.97E-05 Dsc3 2.638301 0.038628 DSG2 3.705926 0.010018 Dsp 8.074286 0.000193 Cldn7 4.063547 0.006624 Tight junction Cldn8 3.761811 0.009379 Cldn4 3.762967 0.009366 Cldn23 2.459445 0.049157 14
Gap junction Gjb2 4.901054 0.002708 Gjd2 2.00267 0.092085 Pcdhb15 15.61573 4.37E-06 Cadm4 3.29845 0.01644 Pvrl4 3.37506 0.01494 Others Mpzl2 3.6082 0.01125 Efs 4.2837 0.00518 Pcdh7 4.69655 0.00333 Csf3r 9.5956 7.33E-05 Thbs1 2.61905 0.03963 15
Supplementary Table 3. The primers used for regular PCR and genomic sequencing Gene Primers for PCR (5'-3') Primers for sequencing (5'-3') Pi3k-ca (E542) Pi3k-ca (E545) Pi3k-ca(H1047) p53 (CDS) NRF2 (CDS) Keap1 (CDS) Cul3 (CDS) p16 (CDS) p19 (CDS) AAGGAGGAGCACTGTCCGTT AAGGAGGAGCACTGTCCGTT TGGGCCACTTCGTCTCTGGA TGGGCCACTTCGTCTCTGGA AAGGAGGAGCACTGTCCGTT AAGGAGGAGCACTGTCCGTT TGGGCCACTTCGTCTCTGGA TGGGCCACTTCGTCTCTGGA GATGTGTTACAAGGCTTACCT GATGTGTTACAAGGCTTACCT CTTGCTCAAGTCCTAATGTT CTTGCTCAAGTCCTAATGTT GGTAGCGACTACAGTTAGGGGGCA GGTAGCGACTACAGTTAGGGGGCA CAGCAGAAGGGACCGGGAGGA CAGCAGAAGGGACCGGGAGGA GCCTCACCTCTGCTGCAAGT GCCTCACCTCTGCTGCAAGT TTTCACATCACAGTAGGAAG TTTCACATCACAGTAGGAAG CTGTGCTTAGTCACCGTGA GCTAGCAGAGGAACTGTGTCT CTTTAGTACAGAGAAGCAGT TGTTCCACGCGTGCATCGAC CCTCCGAGTGCGAGCCGCA CCTCCGAGTGCGAGCCGCA CCTCCGAGTGCGAGCCGCA CCTCCGAGTGCGAGCCGCA GTAGACAGAGGAGCAATAAGA CACTGAATCTCCGCGAGGAA CACTGAATCTCCGCGAGGAA AGCCACATGCTAGACACGCT AGCCACATGCTAGACACGCT TGGGGGCGGCGCTTCTCACCT TGGGGGCGGCGCTTCTCACCT GGCTGAGGCCGGATTTAGCT GGCTGAGGCCGGATTTAGCT 16
Supplementary Table 4. Description and frequency of LKB1 genetic alterations Exonic deletion n(%) Deletion of whole exon 2(1.98) Deletion of exon 1 5(4.95) Deletion of exon10 2(1.98) Deletion of exon 1-4 1(0.99) Deletion of exon 1-6 1(0.99) Deletion of exon 5-8 1(0.99) Deletion of exon 7-8 1(0.99) 17
Supplementary Table 5. The primers used for real-time PCR Gene Forward primer Reverse primer β-actin TGAGCGCAAGTACTCTGTGTGGAT ACTCATCGTACTCCTGCTTGCTGA Brf2 CCTTCCTGGCATGGCAGTCTCT CACCACCGACCGCTTGTTAAGTT Cldn8 GCAACCTACGCTCTTCAAATGG TTCCCAGCGGTTCTCAAACAC cmyc CAACGACAGCAGCTCGCCCA AGCCCGACTCCGACCTCTTGG Col1a1 GCTCCTCTTAGGGGCCACT CCACGTCTCACCATTGGGG Col4a1 CTGGCACAAAAGGGACGAG ACGTGGCCGAGAATTTCACC Col4a5 GGAGAACGGGGGTTTCCAG CTCCCTTGGTTCCATTGCATC Col4a6 GAACTGGCAGAATCGGGACAG TCCAATGGGACCCTTATCTCC Col6a1 CTGCTGCTACAAGCCTGCT CCCCATAAGGTTTCAGCCTCA Col6a2 GCTCCTGATTGGGGGACTCT CCAACACGAAATACACGTTGAC Dcn GTCATCTTCGAGTGGTGCAGT CGGGTGGAAAATCCCAGGG Dlx5 TCTCTAGGACTGACGCAAACA GTTACACGCCATAGGGTCGC DLX5 AAGCGCCACCAACCAGCCAG CGTTCCGGCAAGGCGAGGTA Dsc2 ATGGCGGCTGTGGGATCTAT GCAAGGATCGCAAGGGTCAA Dsg2 CGTGGTTGAAGGCATTCATTTC TAGCTGCTTGACCAGTGTCTT Gjb2 ATCCTCGGGGGTGTCAACAA AGACAAAATCGGCTTGCTCATC Hif1a ACCTTCATCGGAAACTCCAAAG CTGTTAGGCTGGGAAAAGTTAGG Itga6 GCAGAAGCACTCCCGCTGCA TTGCCCCCTGGACCTTGGCT Jup TGGCAACAGACATACACCTACG GGTGGTAGTCTTCTTGAGTGTG Klf5 CCTCCGTCCTATGCCGCTACAA GCTTCTCGCCCGTATGAGTCCT KLF5 TCGGATGAGCTGACCCGCCA TCAGAGCGCGAGAAGCTGCG Krt14 CAGCGGCCCACTGAGATCAAAG CTTGGTCCGGAAGTCATCGGCA Krt15 TTTGGAGGCAGCTCTACCCGAGG CTGCCCCCGAGGAGACAAACC Krt17 CGGGGCCAGCCAGAGACTAC TCGGCTCTGGCCAGGGTCAG Krt5 GGCCCACAGAGACTGCTTCTTT AACATTTTGGGGTCTGGGTCAC Krt6a TCCCAGATGTCATGGCTGCAGAA AGATGAGCAGTGGGGACCCATG Lox TCTTCTGCTGCGTGACAACC GAGAAACCAGCTTGGAACCAG Mmp14 CAGTATGGCTACCTACCTCCAG GCCTTGCCTGTCACTTGTAAA Mmp9 TGTCTGGAGATTCGACTTGAAGTC TGGTGTGCCCTGGAACTCA Oct1 AGCTCTTGCTTCTAGTGGCTC CTGGCTGTAGGTGCAGAGTTC p63 GTGTTGGTGTGGCACAGGGG TCTTCCCCACAGCTCTGGCT Perp ATCGCCTTCGACATCATCGC CCCCATGCGTACTCCATGAG Sccro AGAGTGTAAAAGGATCGTTGGAC CTGGATCGAGGGCCAGATCA 18
Sox2 TGCTGCCTCTTTAAGACTAGGGCT CGGGCGAAGTGCAATTGGGA Spc ATGGACATGAGTAGCAAAGAGGT CACGATGAGAAGGCGTTTGAG Timp1 GCAACTCGGACCTGGTCATAA CGGCCCGTGATGAGAAACT Timp2 TCAGAGCCAAAGCAGTGAGC GCCGTGTAGATAAACTCGATGTC Timp3 CTTCTGCAACTCCGACATCGT GGGGCATCTTACTGAAGCCTC Trim29 AGAATGGCACTAAAGCAGACAG AAATAGGCCACTCTTCCCCTC TRIM29 GGGGTGAGTGGAGTGCACCG GAGCTGCCTTGGACGACGGG 19