We consider the problem of estimating the graph structure associated with a GaussianMarkov random field (GMRF) from i.i.d. samples. We study the performance of study theperformance of the ℓ1-regularized maximum likelihood estimator in the high-dimensionalsetting, where the number of nodes in the graph p, the number of edges in the graph s andthe maximum node degree d, are allowed to grow as a function of the number of samplesn. Our main result provides sufficient conditions on (n, p, d) for the ℓ1-regularized MLEestimator to recover all the edges of the graph with high probability. Under some conditionson the model covariance, we show that model selection can be achieved for sample sizesn = (d2 log(p)), with the error decaying as O(exp(−c log(p))) for some constant c. Weillustrate our theoretical results via simulations and show good correspondences betweenthe theoretical predictions and behavior in simulations.