A genome-wide association study (GWAS), which uses information on single nucleotide polymorphisms (SNPs) from many accessions, has become a powerful approach to gene identification. A metabolome GWAS (mGWAS), which relies on phenotypic information based on metabolite accumulation, can identify genes that contribute to primary and secondary metabolite contents. In this study, we carried out a mGWAS using seed metabolomic data from Arabidopsis thaliana accessions obtained by liquid chromatography-mass spectrometry to identify SNPs highly associated with the contents of metabolites such as glucosinolates. These SNPs were present in genes known to be involved in glucosinolate biosynthesis, thus confirming the effectiveness of our analysis. We subsequently focused on SNPs detected in an unknown methyltransferase gene associated with N-methylhistidine content. Knockout and overexpression of A. thaliana lines of this gene had significantly decreased and increased N-methylhistidine contents, respectively. We confirmed that the overexpressing line exclusively accumulated histidine methylated at the pi position, not at the tau position. Our findings suggest that the identified methyltransferase gene encodes a key enzyme for N-methylhistidine biosynthesis in A. thaliana.