In this paper, we consider N -dimensional real Wishart matrices Y in the class W R (Σ, M ) in which all but one eigenvalues of Σ is 1. Let the non-trivial eigenvalue of Σ be 1 + τ , then as N , M → ∞, with M/N = γ 2 finite and non-zero, the eigenvalue distribution of Y will converge into the Marchenko-Pastur distribution inside a bulk region. When τ increases from zero, one starts to see a stray eigenvalue of Y outside of the support of the Marchenko-Pastur density. As the this stray eigenvalue leaves the bulk region, a phase transition will occur in the largest eigenvalue distribution of the Wishart matrix. In this paper we will compute the asymptotics of the largest eigenvalue distribution when the phase transition occur. We will first establish the results that are valid for all N and M and will use them to carry out the asymptotic analysis. In particular, we have derived a contour integral formula for the Harish-Chandra Itzykson-Zuber integral O(N ) e tr(XgY g T ) g T dg when X, Y are real symmetric and Y is a rank 1 matrix. This allows us to write down a Fredholm determinant formula for the largest eigenvalue distribution and analyze it using orthogonal polynomial techniques. As a result, we obtain an integral formula for the largest eigenvalue distribution in the large N limit characterized by Painlevé transcendents. The approach used in this paper is very different from a recent paper [23], in which the largest eigenvalue distribution was obtained using stochastic operator method. In particular, the Painlevé formula for the largest eigenvalue distribution obtained in this paper is new.