This study presents the results of the identifi-cation and quantification of 12 isoflavones (prunetin, irilone, pseudobaptigenin, glycitein, daidzin, genistin, daidzein, pratensein, puerarin, biochanin A, formononetin and genistein) in 23 species of Trifolium (an MSPD method and analyzed with HPLC coupled with a diode-array detector. The evaluation of molecular phylog-eny of the IFS gene and the relation with isoflavone content was also performed. Five species (T. subterraneum subsp. brachycalycinum, T. alexandrinum, T. pratense, T. sub-terraneum and T. lappaceum) were identified with high levels of biochanin A (431–83 mg/kg), formononetin (72–365 mg/kg) and genistein (9–509 mg/kg), which could be utilized as alternative sources for the nutraceutical industry. Genetic phylogeny for the IFS gene was found in the species studied, with 20 out of 23 species having been divided into two clades, while the remaining three were genetically distant. Based on our results, we confirm the direct correlation between IFS gene polymorphism and isoflavones content in species of Trifolium particularly noted for formononetin. Therefore, the IFS gene can be utilized for screening Trifolium genotypes for formonone-tin. The relation of the three isoflavones' contents and the molecular phylogeny of plants determined by the IFS sequences, as a screening marker for plants with high isoflavone contents in Trifolium species, are to the best of our knowledge described for the first time.