Properties of molecules solvated in ionic liquids (ILs) are strongly affected by solvent environment. For this reason, to give reliable results, ab initio calculations on solutes in ILs, including ions constituting ionic liquid itself, have to self-consistently account for the change of both electronic and classical solvation structure in ILs. Here, we study the electronic structure of the methyl-methylimidazolium ion in the bulk liquid of [mmim][CI] by using the self-consistent field coupling of Kohn-Sham density functional theory and three-dimensional molecular theory of solvation (aka 3D-RISM) with the closure approximation of Kovalenko and Hirata. The KS-DFT/3D-RISM-KH method yields the 3D distribution of the IL solvent species around the [mmim] solute, underlying the most important peculiarities of this kind of systems such as the stacking interaction between neighboring cations, and reproduces the enhancement of the dipole moment resulting from the polarization of the cation by the solvent in a very good agreement with the results of an ab initio MD calculation. The KS-DFT/3D-RISM-KH method offers an accurate and computationally efficient procedure to perform ab initio calculations on species solvated in ionic liquids.