We derive upper and lower bounds on the expectation of $f(\mathbf{S})$ under dependence uncertainty, i.e. when the marginal distributions of the random vector $\mathbf{S}=(S_1,\dots,S_d)$ are known but their dependence structure is partially unknown. We solve the problem by providing improved \FH bounds on the copula of $\mathbf{S}$ that account for additional information. In particular, we derive bounds when the values of the copula are given on a compact subset of $[0,1]^d$, the value of a functional of the copula is prescribed or different types of information are available on the lower dimensional marginals of the copula. We then show that, in contrast to the two-dimensional case, the bounds are quasi-copulas but fail to be copulas if $d>2$. Thus, in order to translate the improved \FH bounds into bounds on the expectation of $f(\mathbf{S})$, we develop an alternative representation of multivariate integrals with respect to copulas that admits also quasi-copulas as integrators. By means of this representation, we provide an integral characterization of orthant orders on the set of quasi-copulas which relates the improved \FH bounds to bounds on the expectation of $f(\mathbf{S})$. Finally, we apply these results to compute model-free bounds on the prices of multi-asset options that take partial information on the dependence structure into account, such as correlations or market prices of other traded derivatives. The numerical results show that the additional information leads to a significant improvement of the option price bounds compared to the situation where only the marginal distributions are known.