Metastasis is the biggest obstacle to esophageal squamous cell carcinoma (ESCC) treatment. Single-cell RNA sequencing analyses are applied to investigate lung metastatic ESCC cells isolated from pulmonary metastasis mouse model at multiple timepoints to characterize early metastatic microenvironment. A small population of parental KYSE30 cell line (Cluster S) resembling metastasis-initiating cells (MICs) is identified because they survive and colonize at lung metastatic sites. Differential expression profile comparisons between Cluster S and other subpopulations identified a panel of 7 metastasis-initiating signature genes (MIS), including CD44 and TACSTD2, to represent MICs in ESCC. Functional studies demonstrated MICs (CD44high) exhibited significantly enhanced cell survival (resistances to oxidative stress and apoptosis), migration, invasion, stemness, and in vivo lung metastasis capabilities, while bioinformatics analyses revealed enhanced organ development, stress responses, and neuron development, potentially remodel early metastasis microenvironment. Meanwhile, early metastasizing cells demonstrate quasi-epithelial-mesenchymal phenotype to support both invasion and anchorage. Multiplex immunohistochemistry (mIHC) staining of 4 MISs (CD44, S100A14, RHOD, and TACSTD2) in ESCC clinical samples demonstrated differential MIS expression scores (dMISs) predict lymph node metastasis, overall survival, and risk of carcinothrombosis.