Background: Non-small cell lung cancer (NSCLC) is the major type of lung cancer with high morbidity and poor prognosis. Erlotinib, an inhibitor of epidermal growth factor receptor (EGFR), has been clinically applied for NSCLC treatment. Nevertheless, the erlotinib acquired resistance of NSCLC occurs inevitably in recent years. Methods: Through analyzing two microarray datasets, erlotinib resistant NSCLC cells microarray (GSE80344) and NSCLC tissue microarray (GSE19188), the differentially expressed genes (DEGs) were screened via R language. DEGs were then functionally annotated by Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, which up-regulated more than 2-folds in both datasets were further functionally analyzed by Oncomine, GeneMANIA, R2, Coremine, and FunRich. Results: We found that matrix metalloproteinase 1 (MMP1) may confer the erlotinib therapeutic resistance in NSCL C. MMP1 highly expressed in erlotinib-resistant cells and NSCLC tissues, and it associated with poor overall survival. In addition, MMP1 may be associated with COPS5 and be involve in an increasing transcription factors HOXA9 and PBX1 in erlotinib resistance. Conclusions: Generally, these results demonstrated that MMP1 may play a crucial role in erlotinib resistance in NSCLC, and MMP1 could be a prognostic biomarker for erlotinib treatment.