Flavonoids present in food, botanicals, and body fluids occur as complex mixtures, and data on their combinatorial estrogenic effects are sparse. Human cell lines that permanently express estrogen receptor (ER) alpha and ERbeta proteins were developed for the measurement of the global estrogenicity of flavonoids in such complex mixtures. The presence of estrogenic ligands, known and unknown, in these mixtures can be detected by activation of an ER-driven luciferase reporter gene. We also examined the effect of hydroxylation on the estrogenic activities of four common flavonoids-apigenin, kaempferol, luteolin, and quercetin, alone and in combination. An inverse relationship was observed between the number of hydroxyl groups in flavonoids and ERalpha bioactivity. When submaximal doses of apigenin, luteolin, kaempferol, genistein, and estradiol were combined in binary and higher order mixtures, the experimental estrogenic effects matched those obtained by summing effects extrapolated from dose-response curves of individual compounds. The estrogenic activities of mixtures containing quercetin were observed to deviate from additivity, suggesting that it was a partial agonist/antagonist. Our assay reveals superagonistic, additive, and antagonistic ERalpha or ERbeta actions of flavonoids and adds to our understanding of the estrogenic effects of phytoestrogens in complex mixtures.