In this paper, we consider a nonlinear fractional diffusion equations with a Riemann-Liouville derivative. First, we establish the global existence and uniqueness of mild solutions under some assumptions on the input data. Some regularity results for the mild solution and its derivatives of fractional orders are also derived. Our key idea is to combine the theories of Mittag-Leffler functions, Banach fixed point theorem and some Sobolev embeddings.