Gauss-Bonnet holographic fluid is a useful theoretical laboratory to study the effects of curvature-squared terms in the dual gravity action on transport coefficients, quasinormal spectra and the analytic structure of thermal correlators at strong coupling. To understand the behavior and possible pathologies of the Gauss-Bonnet fluid in 3 + 1 dimensions, we compute (analytically and non-perturbatively in the Gauss-Bonnet coupling) its second-order transport coefficients, the retarded two-and three-point correlation functions of the energy-momentum tensor in the hydrodynamic regime as well as the relevant quasinormal spectrum. The Haack-Yarom universal relation among the second-order transport coefficients is violated at second order in the Gauss-Bonnet coupling. In the zero-viscosity limit, the holographic fluid still produces entropy, while the momentum diffusion and the sound attenuation are suppressed at all orders in the hydrodynamic expansion. By adding higher-derivative electromagnetic field terms to the action, we also compute corrections to charge diffusion and identify the non-perturbative parameter regime in which the charge diffusion constant vanishes. 4 We use notations and conventions of [2]. See Appendix B and footnote 91 on page 128 of ref.[8] for clarification of sign conventions appearing in the literature. 5 Note that all transport coefficients in H are "dynamical" in terminology of ref. [19]. 6 As advertised in ref.[30] and shown below (and, independently, in ref. [48] using fluid-gravity duality methods), the Haack-Yarom relation does not hold non-pertutbatively in the Gauss-Bonnet coupling. 7 It appears that at weak coupling, the relation H = 0 does not hold. We briefly review the results at weak coupling in Appendix A. 8 We would like to thank P. Kovtun and M. Rangamani for a discussion of these issues.