We prove weighted estimates for rough bilinear singular integral operators with kernel K(y₁, y₂) = (Ω((y₁,y₂)/|(y₁,y₂)|)/|(y₁, y₂)|^2d), where y_i ∈ R^d and Ω ∈ L^∞(S^2d-1) with ∈t_S^2d-1Ω dσ = 0. The argument is by sparse domination of rough bilinear operators, via an abstract theorem that is a multilinear generalization of recent work by Conde-Alonso, Culiuc, Di Plinio and Ou, 2016. We also use recent results due to Grafakos, He, and Honzík, 2015, for the application to rough bilinear operators. In particular, since the weighted estimates are proved via sparse domination, we obtain some quantitative estimates in terms of the A_p characteristics of the weights. The abstract theorem is also shown to apply to multilinear Calderón-Zygmund operators with a standard smoothness assumption. Due to the generality of the sparse domination theorem, future applications not considered in this paper are expected.