Q-systems provide an efficient way of solving Bethe equations. We formulate here Q-systems for both the isotropic and anisotropic open Heisenberg quantum spin-1/2 chains with diagonal boundary magnetic fields. We check these Q-systems using novel Wronskian-type formulas (relating the fundamental Q-function and its dual) that involve the boundary parameters.