The aim of this paper is to derive some Calderón-Zygmund estimates for local solutions u of nonlinear elliptic systems of the type \(\begin{aligned} \textrm{div}\textbf{A}(x, Du) = \textrm{div}|\textbf{G}|^{p - 2}\textbf{G} \qquad \textrm{in} \quad \Omega \subset \mathbb {R}^n, \end{aligned}\) where \(\Omega \subset \mathbb {R}^n\) is a bounded domain. We assume that \(u: \Omega \rightarrow \mathbb {R}^N\) belongs to a weighted Sobolev space \(W_{loc}^{1, p}\) , with \(p \in \left( \frac{2n}{n + 2}, 2\right) \) , \(\textbf{G}\) belongs to a weighted \(L_{loc}^{p}\) space and \(x \rightarrow \textbf{A}(x, \xi )\) has growth coefficients in the John and Nirenberg space BMO. As an application of similar techniques, we also consider weak solutions u of linear elliptic systems of the type \(\begin{aligned} \textrm{div}(A(x)Du) = \textrm{div}\, G, \end{aligned}\) where the matrix A(x) lies in the space BMO. In this case, an improvement of the admissible exponent range in the Calderón-Zygmund estimate is obtained.