For centrosymmetric materials such as monolayer graphene, no optical second harmonic generation (SHG) is generally expected because it is forbidden under the electric-dipole approximation. Yet we observed a strong, doping induced SHG from graphene, with its highest strength comparable to the electric-dipole allowed SHG in non-centrosymmetric 2D materials. This novel SHG has the nature of an electric-quadrupole response, arising from the effective breaking of inversion symmetry by optical dressing with an in-plane photon wave vector. More remarkably, the SHG is widely tuned by carrier doping or chemical potential, being sharply enhanced at Fermi edge resonances, but vanishing at the charge neutral point that manifests the electron-hole symmetry of massless Dirac Fermions. The striking behavior in graphene, which should also arise in graphene-like Dirac materials, expands the scope of nonlinear optics, and holds the promise of novel optoelectronic and photonic applications. 3 Second harmonic generation (SHG) is the most fundamental second-order nonlinear optical process, described by ( ) ( ) ( ) ( ) ( ) [1]. In this process, the output signal is frequency doubled from the incident photon field of ( ).Here, ( ) ( ) is the rank-three nonlinear susceptibility tensor and depends on the incident frequency and photon wave vector q. Since q is typically small, a Taylorthe leading electric-dipole (ED) term, and ( ) ( ) is the often neglected electric-quadrupole/ magnetic-dipole term [1,2] or the EQ response for simplicity. For the electric-dipole allowed SHG to exist, the breaking of inversion symmetry is essential. Hence, SHG is a sensitive probe to symmetry-governed phenomena such as ferroelectricity [3], valley pseudospin [4], and phase transitions [5]. For 2D materials such as hexagonal boron nitride, transition metal dichalcogenide and monochalcogenide, their atomic lattices in monolayer form are non-centrosymmetric, giving rise to the electric-dipole allowed SHG. In fact, SHG has become an indispensable tool to characterize their crystal orientation, stacking symmetry and electronic features [6-10]. Yet for an isolated monolayer graphene, it is centrosymmetric and no electric-dipole SHG is allowed. The third-order optical nonlinearity such as third-harmonic generation [11-15] and four-wave mixing [14,16] was often regarded as the dominant nonlinear process in graphene. Only weak SHG was observed on supported monolayers, which was attributed to the inversion symmetry breaking by the underlying substrate [17] or an in-plane electric current [18,19]. Therefore, graphene provides a 4 unique platform to study unusual SHG responses such as valley polarization induced SHG [20,21] and the EQ response beyond the electric-dipole approximation [22,23], as theoretically proposed. In this work, we exclusively investigate the EQ response of SHG in graphene by introducing an in-plane photon wave vector at oblique incidence, which effectively break the overall inversion symmetry of the system [22-25]. By comparing with the SHG respo...