Roots of kudzu (Pueraria lobata) are a rich source of isoflavone O- and C-glycosides. Although O-glycosylation of (iso)flavonoids has been well characterized at the molecular level, no plant isoflavonoid C-glycosyltransferase genes have yet been isolated. To address the biosynthesis of kudzu isoflavonoids, we generated 6,365 high-quality expressed sequence tags (ESTs) from a subtraction cDNA library constructed using RNA from roots that differentially accumulate puerarin. The ESTs were clustered into 722 TCs and 3,913 singletons, from which 15 family I glycosyltransferases (UGTs) were identified. Hierarchical clustering analysis of the expression patterns of these UGTs with isoflavone synthase (IFS) in a range of tissues identified UGTs with potential functions in isoflavone glycosylation. The open reading frames of these UGTs were expressed in E. coli for functional analysis, and one was shown to preferentially glycosylate isoflavones at the 7-O-position. In addition, ESTs corresponding to chalcone synthase, chalcone reductase, chalcone isomerase (CHI) and 2-hydroxyisoflavanone dehydratase were identified. Recombinant CHI proteins had high activities with both 6'-deoxy- and 6'-hydroxy chalcones, typical of Type II CHIs. Establishment of this EST database and identification of genes associated with kudzu isoflavone biosynthesis and glycosylation provide a new resource for metabolic engineering of bioactive kudzu isoflavones.