  We briefly summarize our recent study on the number of constituent quark
(NCQ) scaling of hadron elliptic flows in high multiplicity p-Pb collisions at
$\sqrt{s_{NN}}=$ 5.02 TeV. With the inclusion of hadron production via the
quark coalescence model at intermediate $p_T$, the viscous hydrodynamics at low
$p_T$, and jet fragmentation at high $p_T$, our $Hydro-Coal-Frag$ model
provides a nice description of the $p_T$-spectra and differential elliptic flow
$v_2(p_T)$ of pions, kaons and protons over the $p_T$ range from 0 to 6 GeV.
Our results demonstrate that including the quark coalescence is essential for
reproducing the observed approximate NCQ scaling of hadron $v_2$ at
intermediate $p_T$ in experiments, indicating strongly the existence of
partonic degrees of freedom and the formation of quark-gluon plasma in high
multiplicity p--Pb collisions at the LHC.
