An asymptotic reduction method is introduced to construct a rod theory for a linearized general anisotropic elastic material for space deformation. The starting point is Taylor expansions about the central line in rectangular coordinates, and the goal is to eliminate the two cross-section spatial variables in order to obtain a closed system for displacement coefficients. This is first achieved, in an ‘asymptotically inconsistent’ way, by deducing the relations between stress coefficients from a Fourier series for the lateral traction condition and the three-dimensional (3D) field equation in a pointwise manner. The closed system consists of 10 vector unknowns, and further refinements through elaborated calculations are performed to extract bending and torsion terms and to obtain recursive relations for the first- and second-order displacement coefficients. Eventually, a system of four asymptotically consistent rod equations for four unknowns (the three components of the central-line displacement and the twist angle) are obtained. Six physically meaningful boundary conditions at each edge are obtained from the edge term in the 3D virtual work principle, and a one-dimensional rod virtual work principle is also deduced from the weak forms of the rod equations.