We present a new release of the turboTDDFT code featuring an implementation of hybrid functionals, a recently introduced pseudo-Hermitian variant of the Liouville-Lanczos approach to time-dependent density-functional perturbation theory, and a newly developed Davidson-like algorithm to compute selected interior eigenvalues/vectors of the Liouvillian super-operator. Our implementation is thoroughly validated against benchmark calculations performed on the cyanin (C$_{21}$O$_{11}$H$_{21}$) molecule using the Gaussian09 and turboTDDFT 1.0 codes.