This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.Abstract This study aims to establish a system for the rapid discrimination of Zoysia species using metabolite fingerprinting of FT-IR spectroscopy combined with multivariate analysis. Whole cell extracts from leaves of 19 identified Zoysia japonica, 6 identified Zoysia sinica, and 38 different unidentified Zoysia species were subjected to Fourier transform infrared spectroscopy (FT-IR). PCA (principle component analysis) and PLS-DA (partial least square discriminant analysis) from FT-IR spectral data successfully divided the 25 identified turf grasses into two groups, representing good agreement with species identification using molecular markers. PC (principal component) loading values show that the 1,100~950 cm -1 region of the FT-IR spectra are important for the discrimination of Zoysia species. A dendrogram based on hierarchical clustering analysis (HCA) from the PCA and PLS-DA data of turf grasses showed that turf grass samples were divided into Zoysia japonica and Zoysia sinica in a species-dependent manner. PCA and PLS-DA from FT-IR spectral data of Zoysia species identified and unidentified by molecular markers successfully divided the 49 turf grasses into Z. japonica and Z. sinica. In particular, PLS-DA and the HCA dendrogram could mostly discriminate the 47 Z. japonica grasses into two groups depending on their origins (mountainous areas and island area). Considering these results, we suggest that FT-IR fingerprinting combined with multivariate analysis could be applied to discriminate between Zoysia species as well as their geographical origins of various Zoysia species.