  We initiate the study of Regge theory in a bottom-up holographic model for
QCD in the Veneziano limit, where the backreaction of the quarks to the gluon
dynamics is included. We determine the parameters of the model by carrying out
a precise fit to the meson spectrum in QCD. The spectrum for spin-one and
pseudoscalar mesons is well reproduced. We then generalise the model to incluce
higher spin fields in the bulk trajectories dual to the Pomeron and meson Regge
trajectories at the boundary. With this setting, we fit the masses of the
mesons with spins $J=2$, $3$, and $4$,as well as the experimental data of the
total cross-sections $\sigma(\gamma \gamma \to X)$, $\sigma(\gamma p \to X)$
and $\sigma(p p \to X)$. For the cross sections we obtain a
$\chi^2_{\mathrm{d.o.f.}}$ of $0.74$ for a total of 199 experimental points.
