We report an ab initio study on the spectroscopy of the open-shell diatomic molecule yttrium oxide, YO. The study considers the six lowest doublet states, $X\,{}^{2}\Sigma^{+}$, $A\,{}^{2}\Pi$, $C\,{}^{2}\Pi$, $A'\,{}^{2}\Delta$, $B\,{}^{2}\Sigma^{+}$, $D\,{}^{2}\Sigma^{+}$ and a few higher-lying quartet states using high levels of electronic structure theory and accurate nuclear motion calculations. The coupled cluster singles, doubles, and perturbative triples, CCSD(T), and multireference configuration interaction (MRCI) methods are employed in conjunction with a relativistic pseudopotential on the yttrium atom and a series of correlation-consistent basis sets ranging in size from triple-$\zeta$ to quintuple-$\zeta$ quality. Core-valence correlation effects are taken into account and complete basis set limit extrapolation is performed for CCSD(T). Spin-orbit coupling is included through the use of both MRCI state-interaction with spin-orbit (SI-SO) approach and four-component relativistic equation-of-motion CCSD calculations. Using the ab initio data for bond lengths ranging from 1.0 to 2.5 A, we compute 6 potential energy, 12 spin-orbit, 8 electronic angular momentum, 6 electric dipole moment and 12 transition dipole moment (4 parallel and 8 perpendicular) curves which provide a complete description of the spectroscopy of the system of six lowest doublet states. The Duo nuclear motion program is used to solve the coupled nuclear motion Schr\"{o}dinger equation for these six electronic states. The spectra of $^{89}$Y$^{16}$O simulated for different temperatures are compared with several available high resolution experimental studies; good agreement is found once minor adjustments are made to the electronic excitation energies.