We characterize asymptotic collective behavior of rectangular random matrices, the sizes of which tend to infinity at different rates. It appears that one can compute the limits of all noncommutative moments (thus all spectral properties) of the random matrices we consider because, when embedded in a space of larger square matrices, independent rectangular random matrices are asymptotically free with amalgamation over a subalgebra. Therefore, we can define a "rectangular-free convolution", which allows to deduce the singular values of the sum of two large independent rectangular random matrices from the individual singular values. This convolution is linearized by cumulants and by an analytic integral transform, that we called the "rectangular R-transform".