The problem of scheme and gauge dependence of the factorization property of the renormalization group β-function in the SU(Nc) QCD generalized Crewther relation (GCR), which connects the flavor non-singlet contributions to the Adler and Bjorken polarized sum rule functions, is investigated at the mathcal{O}left({a}_s^4right) level of perturbation theory. It is known that in the gauge-invariant renormalization overline{mathrm{MS}} -scheme this property holds in the QCD GCR at least at this order. To study whether this factorization property is true in all gauge-invariant schemes, we consider the MS-like schemes in QCD and the QED-limit of the GCR in the overline{mathrm{MS}} -scheme and in two other gauge-independent subtraction schemes, namely in the momentum MOM and the on-shell OS schemes. In these schemes we confirm the existence of the β-function factorization in the QCD and QED variants of the GCR. The problem of the possible β-factorization in the gauge-dependent renormalization schemes in QCD is studied. To investigate this problem we consider the gauge non-invariant mMOM and MOMgggg-schemes. We demonstrate that in the mMOM scheme at the mathcal{O}left({a}_s^3right) level the β-factorization is valid for three values of the gauge parameter ξ only, namely for ξ = −3, −1 and ξ = 0. In the mathcal{O}left({a}_s^4right) order of PT it remains valid only for case of the Landau gauge ξ = 0. The consideration of these two gauge-dependent schemes for the QCD GCR allows us to conclude that the factorization of RG β-function will always be implemented in any MOM-like renormalization schemes with linear covariant gauge at ξ = 0 and ξ = −3 at the mathcal{O}left({a}_s^3right) approximation. It is demonstrated that if factorization property for the MS-like schemes is true in all orders of PT, as theoretically indicated in the several works on the subject, then the factorization will also occur in the arbitrary MOM-like scheme in the Landau gauge in all orders of perturbation theory as well.