The variability data for the BL Lacertae object TXS 1902+556 in the optical and γ-ray wavebands were obtained from the 0.76 m Katzman Automatic Imaging Telescope and the Fermi Large Area Telescope (Fermi-LAT), covering periods of 14.4 years and 14.7 years, respectively. The variability properties were systematically analyzed, with particular emphasis on the first comprehensive investigation of radiation variation in the optical waveband. Four well-established techniques were employed for this purpose: the Lomb-Scargle periodogram, REDFIT program, Jurkevich method, and discrete correlation function (DCF) approach. The optical waveband exhibits quasi-periodic oscillations (QPO) with a timescale of PO = 276.8 ± 6.1 days at a significance level 3.87σ, while the γ-ray waveband does not exhibit any significant periodicity. However, it should be noted that the QPO timescale is consistent with the Sun-gaps in the optical light curve within 2σ uncertainties. The optical QPO behavior is most likely attributed to the helical motion of the jet driven by the orbital motion in a supermassive black hole binary system. Moreover, we have provided an explanation for the absence of QPO in the γ-ray light curves. Furthermore, utilizing the DCF method, a weak correlation between the variability in the optical and γ-ray wavebands was observed, suggesting that the emission of TXS 1902+556 may be generated through a combination of synchrotron self-Compton (SSC) and external Compton (EC) processes, or a leptonic-hadronic hybrid process.