Light variability is an essential attribute of young stellar objects (YSOs) at the pre-main-sequence (PMS) stage, which can be treated as an indicator to trace star formation. We present a statistical study on the infrared NEOWISE variability of new and large-sized PMS stars as well as Herbig Ae/Be (HAeBe) and classical Be (CBe) stars identified from Gaia Data Release 2. A total of 3792 (including 714 HAeBe), 417, and 473 NEOWISE variable sources were detected from the catalog of 8470 PMS (including 1361 HAeBe), 693 CBe, and 1309 uncertain (possibly belonging to either HAeBe or CBe) stars, respectively. We found that 80% of the PMS variables are dominated by irregular infrared variability, whereas variables in the CBe and uncertain categories exhibit comparable contributions from secular (linear, curved, and periodic) and stochastic (burst, drop, and irregular) variables. This, in combination with infrared W1-W2 colors, W4 luminosities, and the spectral energy distribution slope at near- and mid-infrared wavelengths, suggests different evolutionary stages of the PMS (including HAeBe), uncertain, and CBe stars, from less to more evolved. A statistical analysis of variable source properties (e.g., period and fractional amplitude) further reveals a potential sequence possibly associated with disk evolution for different types of light curves (from linear, to curved, periodic, stochastic, and nonvariable stars). Possible mechanisms for the variable sources are suggested to be related to either extinction via circumstellar matter, or contribution of multiple systems surrounding the YSOs, as opposed to variable accretion (e.g., episodic accretion or accretion burst).