We compute the β 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 γ_{5} in D=4-2ε space-time dimensions is required; however, in our calculation we determine the γ_{5}-dependent terms by exploiting relations with β function coefficients at lower loop orders.