Estrogen receptor (ER) biology reflects the actions of estrogens through the two receptors, ERα and ERβ, although little is known regarding the preference for formation of ER homo- vs. heterodimers, and how this is affected by the level of ligand occupancy and preferential ligand affinity for one of the ER subtypes. In this report, we use a split optical reporter-protein complementation system to demonstrate the physical interaction between ERα and ERβ in response to different ER ligands in cells and, for the first time, by in vivo imaging in living animals. The genetically encoded reporter vectors constructed with the ligand-binding domains of ERα and ERβ, fused to split firefly or Renilla luciferase (Fluc or hRluc) fragments, were used for this study. This molecular proteomic technique was used to detect ERα/ERα or ERβ/ERβ homodimerization, or ERα/ERβ heterodimerization induced by ER subtype-selective and nonselective ligands, and selective ER modulators (SERM), as well as in dimers in which one mutant monomer was unable to bind estradiol. The SERM-bound ERα and ERβ form the strongest dimers, and subtype-preferential homodimerization was seen with ERα-selective ligands (methyl piperidino pyrazole/propyl pyrazole triol) and the ERβ-selective ligands (diarylpropionitrile/tetrahydrochrysene/genistein). We also demonstrated that a single ligand-bound monomer can form homo- or heterodimers with an apo-monomer. Xenografts of human embryonic kidney 293T cells imaged in living mice by bioluminescence showed real-time ligand induction of ERα/ERβ heterodimerization and reversal of dimerization upon ligand withdrawal. The results from this study demonstrate the value of the split luciferase-based complementation system for studying ER-subtype interactions in cells and for evaluating them in living animals by noninvasive imaging. They also probe what combinations of ERα and ERβ dimers might be the mediators of the effects of different types of ER ligands given at different doses.