We study gravitational waves propagating through an anisotropic Bianchi I dust-filled universe (containing the Einstein-de-Sitter universe as a special case). The waves are modeled as small perturbations of this background cosmological model and we choose a family of null hypersurfaces in this space-time to act as the histories of the wavefronts of the radiation. We find that the perturbations we generate can describe pure gravitational radiation if and only if the null hypersurfaces are shear-free. We calculate the gauge-invariant small perturbations explicitly in this case. How these differ from the corresponding perturbations when the background space-time is isotropic is clearly exhibited.