Dirac fermion in topological materials exhibits intriguing nonlinear optical responses. However, their direct correlation with the linearly dispersed band remains elusive experimentally. Here, we take topological semimetal ZrSiS as a paradigm, unveiling three unique nonlinear optical signatures of Dirac fermion. These signatures include strong quadrupolar response, quantum interference effect, and exponential divergent four-wave mixing (FWM), all of which are described by the prominent third-order nonlinear optical susceptibility. Resonantly enhanced by linear bands, quadrupolar second harmonic generation in centrosymmetric bulk overwhelms the electric-dipole contribution at the surface with inherent inversion symmetry breaking. Furthermore, owing to the interference between multiple resonant transition pathways within linear bands, difference-frequency FWM is several orders of magnitude stronger than sum-frequency FWM and third harmonic generation. The difference-frequency FWM further displays an inverse-square divergence toward degenerate excitation, whose scaling law perfectly matches with the long-sought behavior of Dirac fermion. These signatures lay the solid foundation toward the practical applications of topological materials in nonlinear optoelectronics and photonics.