  We compute the beta functions of the three Standard Model gauge couplings to
four-loop order in the modified minimal subtraction scheme. At this order a
proper definition of $\gamma_5$ in $D=4-2\epsilon$ space-time dimensions is
required; however, in our calculation we determine the $\gamma_5$-dependent
terms by exploiting relations with beta function coefficients at lower loop
orders.
