  We suggest simple ways of implementing Peccei-Quinn (PQ) symmetry to solve
the strong CP problem in renormalizable SUSY $SO(10)$ models with a minimal
Yukawa sector. Realistic fermion mass generation requires that a second pair of
Higgs doublets survive down to the PQ scale. We show how unification of gauge
couplings can be achieved in this context. Higgsino mediated proton decay rate
is strongly suppressed by a factor of $(M_{\rm PQ}/M_{\rm GUT})^2$, which
enables all SUSY particles to have masses of order TeV. With TeV scale SUSY
spectrum, $p \rightarrow \overline{\nu} K^+$ decay rate is expected to be in
the observable range. Lepton flavor violating processes $\mu \rightarrow
e\gamma$ decay and $\mu-e$ conversion in nuclei, induced by the Dirac neutrino
Yukawa couplings, are found to be within reach of forthcoming experiments.
