Kinematics, kinetics and biomechanics of human gait are widely investigated fields of research. The biomechanics of locomotion has been described characterizing muscle activations and synergistic control, i.e, spatial and temporal patterns of coordinated muscle groups and joints. Both kinematic synergies and muscle synergies have been extracted from locomotion data, showing that in healthy people 4-5 synergies underlie human locomotion; such synergies are in general robust across subjects and might be altered in pathological gait, depending on the severity of the impairment. In this work, for the first time, we apply the mixed matrix factorization algorithm to locomotion data of 15 healthy participants to extract hybrid kinematic-muscle synergies, and show that they allow to directly link task space variables (i.e., kinematics) to the neural structure of muscle synergies. We show that kinematic-muscle synergies can describe the biomechanics of motion at a better extent than muscle synergies or kinematic synergies alone. Morevoer, this study shows that at a functional level, modular control of the lower limb during locomotion is underlied by an increased number of functional synergies with respect to standard muscle synergies and account for different biomechanical roles that each synergy may have within the movement. Kinematic-muscular synergies may have impact in future work for a deeper understanding of modular control and neuro-motor recovery in the medical and rehabilitation fields, as they associate neural and task space variables in the same factorization, including the evaluation of post-stroke, Parkinson and cerebral palsy patients, and in other fields such as sports.