Abstract Isoflavonoids are commonly found in leguminous plants. Glycitein is one of the isoflavones produced by soybean. The genes encoding the enzymes in the isoflavone biosynthetic pathway have mostly been identified and characterized. However, the gene(s) for isoflavone O-methyltransferase (IOMT), which catalyzes the last step of glycitein biosynthesis, has not yet been identified. In this study, we conducted multi-omics analyses of fungal-inoculated soybean and indicated that glycitein biosynthesis was induced in response to biotic stress. Moreover, we identified a unique type of IOMT, which participates in glycitein biosynthesis. Soybean seedlings were inoculated with Aspergillus oryzae or Rhizopus oligosporus and sampled daily for 8 d. Multi-omics analyses were conducted using liquid chromatography–tandem mass spectrometry and RNA sequencing. Metabolome analysis revealed that glycitein derivatives increased following fungal inoculation. Transcriptome co-expression analysis identified two candidate IOMTs that were co-expressed with the gene encoding flavonoid 6-hydroxylase (F6H), the key enzyme in glycitein biosynthesis. The enzymatic assay of the two IOMTs using respective recombinant proteins showed that one IOMT, named as GmIOMT1, produced glycitein. Unlike other IOMTs, GmIOMT1 belongs to the cation-dependent OMT family and exhibited the highest activity with Zn2+ among cations tested. Moreover, we demonstrated that GmIOMT1 overexpression increased the levels of glycitein derivatives in soybean hairy roots when F6H was co-expressed. These results strongly suggest that GmIOMT1 participates in inducing glycitein biosynthesis in response to biotic stress.