We prove Li-Yau and Harnack inequalities for systems of linear reaction-diffusion equations. By introducing an additional discrete spatial variable, the system is rewritten as a scalar diffusion equation with an operator sum. For such operators in a mixed continuous and discrete setting, we introduce the hybrid curvature-dimension condition $CD_{hyb} (κ,d)$, which is a combination of the Bakry-Émery condition $CD(κ,d)$ and one of its discrete analogues, the condition $CD_Υ(κ,d)$. We establish a hybrid tensorisation principle and prove that under $CD_{hyb} (0,d)$ with $d<\infty$ a differential Harnack estimate of Li-Yau type holds, from which a Harnack inequality can be deduced by an integration argument.