ABSTRACT In the manuscript, a 3.8 yr optical quasi-periodic oscillations (QPOs) is reported in blue quasar SDSS J132144+033055 (= SDSS J1321) at z = 0.269, based on 16.3yr-long light curve from both CSS and ZTF directly described by a sinusoidal function. The 3.8 yr QPOs can be confirmed through the Generalized Lomb–Scargle periodogram with confidence level higher than 5σ, through properties of the phase-folded light curve and the WWZ technique. Moreover, the collected Pan-STARRS light curves well follow the sinusoidal function described best fitting results to the Catalina Sky Survey (CSS) and Catalina Sky Survey (ZTF) light curves. The optical QPOs strongly indicate a central binary black hole (BBH) system in SDSS J1321, with expected space separation smaller than 0.018 pc, through the estimated upper limit of total BH mass 3.3 × 109 M⊙ through the correlation between BH mass and continuum luminosity. Meanwhile, we check disc precession applied to explain the optical QPOs. However, under the disc precession assumption, the determined optical emission regions from central BH have sizes about 37RG similar as the sizes 35RG of the expected NUV emission regions through the correlation between disc size and BH mass, indicating the disc precession is not preferred. And due to undetected radio emissions, jet precession can be ruled out. Furthermore, only 0.1 per cent probability can determined as the QPOs mis-detected through CAR process randomly created light curves related to intrinsic AGN activities, re-confirming the optical QPOs with significance level higher than 3σ. Therefore, combining long-term light curves from CSS and ZTF can lead to more QPOs candidates in the near future.