The $J/\Psi$ photo-production reactions on the nucleon can provide information on the roles of gluons in determining the $J/\Psi$-nucleon ($J/\Psi$-N) interactions and the structure of the nucleon. The information on the $J/\Psi$-N interactions is needed to test lattice QCD (LQCD) calculations and to understand the nucleon resonances such as $N^*(P_c)$ recently reported by the LHCb Collaboration. In addition, it is also needed to investigate the production of nuclei with hidden charms and to extract the gluon distributions in nuclei. The main purpose of this article is to review six reaction models of $\gamma + p \rightarrow J/\Psi +p$ reactions which have been and can be applied to analyze the data from Thomas Jefferson National Accelerator Facility (JLab). The formulae for each model are given and used to obtain the results to show the extent to which the available data can be described. The models presented include the Pomeron-exchange model of Donnachie and Landshoff (Pom-DL) and its extensions to include $J/\Psi$-N potentials extracted from LQCD (Pom-pot) and to also use the constituent quark model (CQM) to account for the quark substructure of $J/\Psi$ (Pom-CQM). The other three models are developed from applying the perturbative QCD approach to calculate the two-gluon exchange using the generalized parton distribution (GPD) of the nucleon (GPD-based), two- and three-gluon exchanges using the parton distribution of the nucleon ($2g+3g$), and the exchanges of scalar ($0^{++}$) and tensor ($2^{++}$) glueballs within the holographic formulation (holog). The results of investigating the excitation of the nucleon resonances $N^*(P_c)$ in the $\gamma + p \rightarrow J/\Psi +p$ reactions are also given.