We developed lentiviral vectors that express Cre-recombinase (Lenti-Cre)
10 and determined the dose that results in 5 to 20 lung tumours per
KrasLSL-G12D/+;p53flox/flox mouse after intratracheal administration. These mice lived 8–14 months after tumour initiation and developed macroscopic metastases to the draining lymph nodes, pleura, kidneys, heart, adrenal glands, and liver (
Supplementary Fig. 1). Because lentiviruses integrate stably into the genome, the integration site was a unique molecular identifier that unambiguously linked primary tumours to their related metastases (). We used linker-mediated PCR (LM-PCR) to determine the genomic sequence directly 3′ of the integrated lentiviral genome followed by a specific PCR for the lentiviral integration site (). To have samples of sufficient quantity and purity for our analyses, we derived cell lines from primary tumours and metastases. Cell lines were pure tumour cells as determined by recombination of the
p53floxed alleles (data not shown). The clonal relationship of these cell lines was established using LM-PCR or Southern blot analysis for the lentiviral genome ( and data not shown). We termed cell lines derived from verified metastatic primary lung tumours T
Met.
Gene expression profiling was performed on cell lines from twenty-three lung tumours and metastases (nine metastases, seven T
Met primary tumours, and seven potentially non-metastatic primary tumours). Using unsupervised consensus clustering
11, we identified four cell lines from likely non-metastatic tumour samples that had highly concordant gene expression and were separate from all T
Met and metastasis (Met) samples (
Supplementary Fig. 2). Therefore, we surmised that these could represent non-metastatic primary tumours and classified them as T
nonMet. These T
nonMet cell lines consistently formed fewer tumour nodules in the liver after intrasplenic injection despite equivalent proliferation rates in cell culture ( and
Supplementary Fig. 2).
Significant gene expression alterations distinguished T
nonMet from T
Met and Met-derived cell lines ( and
Supplementary Table 1), many of which were validated by qRT-PCR, flow cytometry, and western blotting (data not shown). A gene expression signature generated by comparing T
nonMet to T
Met/Met samples predicted patient outcome in human lung adenocarcinoma gene expression datasets
12,13, suggesting the possibility of evolutionarily-conserved molecular mechanisms of tumour progression (
Supplementary Fig. 2). Thus, we integrated mouse and human data by comparing the differences in expression between T
nonMet and T
Met/Met samples with the association of human gene expression and patient survival (). Two genes were particularly notable from this analysis: the NK-related homeobox transcription factor
Nkx2-1 and the Nkx2-1 target gene surfactant protein B (
Sftpb; ). Nkx2-1 regulates lung development and is expressed in Type II pneumocytes and bronchiolar cells in the adult
14–16.
Nkx2-1 expression was >10-fold higher in T
nonMet samples, and higher
NKX2-1 expression in human tumours correlated with longer survival. Of note,
NKX2-1 is focally amplified in ~10% of human lung adenocarcinoma, with functional data supporting oncogenic activity
6–9,17. Conversely, most immunohistochemical analyses of NKX2-1 in this disease suggest an association between NKX2-1-negative tumours and poor patient outcome
17,18. Thus, we focused on validating and characterizing the function of this transcription factor in suppressing tumour progression and metastasis.
We confirmed reduced Nkx2-1 mRNA and protein in T
Met and Met cell lines without evidence of focal genomic loss of this region (,
Supplemental Fig. 4, and data not shown). Nkx2-1 was consistently downregulated in high-grade poorly differentiated tumours from our mouse model ( and
Supplementary Fig. 3). Nkx2-1 expression was also reduced in advanced Kras
G12D-driven lung adenocarcinomas with p53
R270H or p53
R172H point mutations
4,19. Using our LM-PCR assay, we identified three primary lung tumours as metastatic based on the presence of metastases with the same lentiviral integration site ( and data not shown). These tumours each contained poorly-differentiated areas that were Nkx2-1
neg (
Supplementary Fig. 6). Interestingly, Nkx2-1 expression was low/absent in almost all lymph node and distant macrometastases, though some micrometastases were Nkx2-1
pos or Nkx2-1
mixed (
Supplementary Fig. 3). Whether certain micrometastases were seeded by Nkx2-1
pos cells or reverted to an Nkx2-1
pos phenotype due to cues from their new environment is unknown.
In human lung adenocarcinoma
12,13 the expression of
NKX2-1 correlated with a mouse T
nonMet gene expression signature (
Supplementary Fig. 3). Additionally, the T
nonMet signature was anti-correlated with an embryonic stem cell signature, supporting the notion that T
Met/Met cells have transitioned to a less differentiated and more stem-like state
20 (
Supplementary Fig. 3).
The correlative mouse and human data were consistent with Nkx2-1 being either a marker or a functional regulator of tumour progression. Nkx2-1 expression in a T
Met cell line greatly suppressed tumour formation after intravenous transplantation (, and
Supplementary Fig. 5). Moreover, of the tumours that formed after injection of T
Met-Nkx2-1 cells, many were either Nkx2-1
neg or Nkx2-1
mixed (). In general, tumours that continued to express Nkx2-1 were well-differentiated, while Nkx2-1
neg tumours often displayed solid architecture or areas of poorly-differentiated cells ( and
Supplementary Fig. 5). Intrasplenic transplantation unveiled a similar diminution of tumour formation by T
Met–
Nkx2-1 cells (
Supplementary Fig. 5). Nkx2-1 expression did not alter proliferation or cell death in cell culture, or affect established tumour proliferation
in vivo (
Supplementary Fig. 5 and data not shown), but dramatically reduced the ability of these cells to grow in anchorage-independent conditions and initiate tumours after subcutaneous transplantation ( and
Supplementary Fig. 5)
To further elucidate the function of Nkx2-1, we knocked-down Nkx2-1 in T
nonMet cell lines (
Supplementary Fig. 7).
Nkx2-1 knockdown allowed the formation of more liver nodules after intrasplenic injection and more lung nodules after intravenous transplantation (). Nkx2-1 knockdown did not alter proliferation or cell death in cell culture (
Supplementary Fig. 7) but enhanced the cells’ ability to form colonies under anchorage-independent conditions and tumours after subcutaneous transplantation ( and
Supplementary Fig. 7). Reexpression of an shRNA-insensitive
Nkx2-1 cDNA (
Nkx2-1*) reverted the phenotypic alterations elicited by
shNkx2-1 confirming that the effects of the
shNkx2-1 were specifically due to Nkx2-1 knockdown (
Supplementary Fig. 7). Finally, we induced tumours in
KrasLSL-G12D/+;p53flox/flox mice with either Lenti-Cre or a lentiviral vector expressing both Nkx2-1 and Cre (Lenti-Nkx2-1/Cre). Expression of exogenous Nkx2-1 limited tumour progression resulting in fewer tumours of advanced histopathological grades ().
To discover Nkx2-1 regulated genes, we compared gene expression in T
nonMet and T
nonMet–
shNkx2-1 cells. Overlapping this gene list with the genes expressed at different levels in T
nonMet versus T
Met/Met cells uncovered high priority candidate genes (
Supplementary Fig. 8). We elected to focus on Hmga2 given its role in altering global gene expression through the regulation of chromatin structure and its association with embryonic and adult stem cell states
21–24 as well as malignant tumours of diverse origins
25–29. Hmga2 was derepressed by Nkx2-1 knockdown in T
nonMet cells and regions of
KrasG12D/+;p53Δ/Δ tumours that lacked Nkx2-1 expression were almost universally Hmga2
pos (). Importantly, Nkx2-1
neg areas of known metastatic primary tumours and metastases were also Hmga2
pos (
Supplementary Fig. 9 and data not shown). Additionally, Hmga2 was downregulated in T
Met cells after expression of
Nkx2-1 cDNA and in T
nonMet-
shNkx2-1 cells after reexpression of Nkx2-1* (data not shown).
Although Hmga2 can be regulated by the Let7 family of miRNAs
21,25,30, Let7 levels,
Lin28 expression, and Let7 activity were equivalent in T
nonMet, T
Met, and Met cell lines and were unaltered in T
nonMet-
shNkx2-1 cells (
Supplementary Fig. 10 and data not shown).
Hmga2 promoter activity was derepressed in T
nonMet-
shNkx2-1 cells and repressed in T
Met-
Nkx2-1 cells, indicating that expression of Hmga2 in lung adenocarcinoma cells is regulated at least in part through differential promoter activity (
Supplementary Fig. 10).
We hypothesized that lung adenocarcinomas progress from an Nkx2-1
posHmga-2
neg to an Nkx2-1
negHmga-2
pos state. However, metastatic and non-metastatic tumours could be fundamentally distinct at the time of initiation. Hmga2 is highly expressed in embryonic lung but not in any normal adult lung cells, and early after initiation,
KrasG12D/+;p53Δ/Δ tumours were uniformly Nkx2-1
posHmga-2
neg ( and
Supplementary Fig. 11).
KrasG12D/+;p53-proficient tumours, which maintain their differentiated phenotype and never metastasize even late after tumour initiation
5, were almost universally Nkx2-1
posHmga-2
neg (
Supplementary Fig. 11). Poorly differentiated areas of
KrasG12D/+;p53Δ/Δ tumours with reduced Nkx2-1 expression were never found as
in situ lesions and were almost always associated with lower grade Nkx2-1-expressing areas (
Supplementary Fig. 6). Finally, we induced
KrasG12D/+;p53Δ/Δ tumours with a pool of lentiviral vectors that contain nucleotide barcodes. Amplification and sequencing of the lentivirus-encoded barcodes from adjacent low-grade Nkx2-1
posHmga2
neg and high-grade Nkx2-1
negHmga2
pos areas showed that these areas were clonally related (
Supplemental Fig. 12). While alternate mechanisms leading to the generation of clonally-related but phenotypically-distinct tumour cell populations are possible, including the expansion of rare Nkx2-1
neg cells that pre-exist within the tumour, we believe that our data strongly suggest that lung adenocarcinomas undergo a transition from an Nkx2-1
posHmga-2
neg state to a more aggressive Nkx2-1
negHmga-2
pos state. Our data additionally indicate that an Nkx2-1-dependent gene expression program is a key regulator of this transition.
We next analyzed the expression of NKX2-1 and HMGA2 in human adenocarcinoma. Although the expression patterns were diverse, two important conclusions could be made. First, tumours of NKX2-1
posHMGA2
neg and NKX2-1
negHMGA2
pos phenotypes exist within the spectrum of human lung adenocarcinoma ( and
Supplementary Fig. 10). Second, there was a trend towards well-differentiated tumours being NKX2-1
pos/HMGA2
neg whereas moderately and poorly differentiated tumours were more often represented by other combinations of these proteins. Most notably, the moderately and poorly differentiated groups contained NKX2-1
neg/HMGA-2
pos tumours (). These results underscore the diversity within this single human tumour type and suggest that our genetically defined model likely represents, at the molecular level, a subset of these tumours.
Next we knocked-down
Hmga2 in T
nonMet-
shNkx2-1 cells and found that their metastasis seeding potential was greatly reduced after transplantation ( and
Supplementary Fig. 13). Additionally, Hmga2 knockdown in a metastasis-derived cell line reduced its anchorage-independent growth and tumour seeding ability after transplantation ( and
Supplementary Fig. 13). A future challenge will be to understand the molecular mechanism by which Hmga2 controls lung adenocarcinoma metastatic potential. The expansion of Nkx2-1
negHmga-2
pos regions within primary lung tumours suggests the acquisition of phenotypes that are advantageous to the primary tumours and also increase the probability of metastatic spread.
That NKX2-1 can have both oncogenic and tumour suppressive functions in lung cancer presumably illustrating context-dependent functions within individual tumours of the same type. Lung adenocarcinomas may differ in their cell of origin, mutation spectrum, or gene expression leading to distinct requirements for continued NKX2-1 expression and different capacity to tolerate or benefit from NKX2-1 downregulation. Our studies uncovered the molecular and cellular basis for the association of NKX2-1 expression with good patient outcome
17,18 and HMGA2 expression with poor patient outcome
26,27. Our results emphasize the power of genetically-engineered mouse models of advanced disease, used in conjunction with human studies, to elucidate mechanisms that control cancer progression and metastatic spread. Through this approach we identified one molecular mechanism by which a highly prevalent tumour type can progresses to its malignant state.