This paper proposes a novel Pearson-type quasi maximum likelihood estimator (QMLE) of GARCH(p, q) models. Unlike the existing Gaussian QMLE, Laplacian QMLE, generalized non-Gaussian QMLE, or LAD estimator, our Pearsonian QMLE (PQMLE) captures not just the heavy-tailed but also the skewed innovations. Under strict stationarity and some weak moment conditions, the strong consistency and asymptotic normality of the PQMLE are obtained. With no further efforts, the PQMLE can be applied to other conditionally heteroskedastic models. A simulation study is carried out to assess the performance of the PQMLE. Two applications to four major stock indexes and two exchange rates further highlight the importance of our new method. Heavy-tailed and skewed innovations are often observed together in practice, and the PQMLE now gives us a systematical way to capture these two co-existing features.