We consider the chemotaxis-fluid system n t + u · ∇n = ∇ · (D(n)∇n) − ∇ · (nS(x, n, c) · ∇c), c t + u · ∇c = ∆c − nf (c),in a bounded convex domain Ω ⊂ R 3 with smooth boundary, where φ ∈ W 1,∞ (Ω) and D, f and S are given functions with values in [0, ∞), [0, ∞) and R 3×3 , respectively.In the existing literature, the derivation of results on global existence and qualitative behavior essentially relies on the use of energy-type functionals which seem to be available only in special situations, necessarily requiring the matrix-valued S to actually reduce to a scalar function of c which, along with f , in addition should satisfy certain quite restrictive structural conditions.The present work presents a novel a priori estimation method which allows for removing any such additional hypothesis: Besides appropriate smoothness assumptions, in this paper it is only required that f is locally bounded in [0,for all n ≥ 0 with some k D > 0 and someIt is shown that then for all reasonably regular initial data, a corresponding initial-boundary value problem for (0.1) possesses a globally defined weak solution.The method introduced here is efficient enough to moreover provide global boundedness of all solutions thereby obtained in that, inter alia, n ∈ L ∞ (Ω × (0, ∞)). Building on this boundedness property, it can finally even be proved that in the large time limit, any such solution approaches the spatially homogeneous equilibrium (n 0 , 0, 0) in an appropriate sense, where n 0 := 1 |Ω| Ω n 0 , provided that merely n 0 ≡ 0 and f > 0 on (0, ∞). To the best of our knowledge, these are the first results on boundedness and asymptotics of large-data solutions in a three-dimensional chemotaxisfluid system of type (0.1).