Large-scale Monte Carlo simulations are employed to study phase transitions in the three-dimensional compact Abelian Higgs model in adjoint representations of the matter field, labeled by an integer q, for $q=2,3,4,5.$ We also study various limiting cases of the model, such as the ${Z}_{q}$ lattice gauge theory, dual to the three-dimensional (3D) spin model, and the 3D $\mathrm{XY}$ spin model which is dual to the ${Z}_{q}$ lattice gauge theory in the limit $\stackrel{\ensuremath{\rightarrow}}{q}\ensuremath{\infty}.$ In addition, for benchmark purposes, we study the square lattice eight-vertex model, which is exactly solvable and features nonuniversal critical exponents. We have computed the first, second, and third moments of the action to locate the phase transition of the compact Abelian Higgs model in the parameter space $(\ensuremath{\beta},\ensuremath{\kappa}),$ where $\ensuremath{\beta}$ is the coupling constant of the matter term and $\ensuremath{\kappa}$ is the coupling constant of the gauge term. We have found that for $q=3,$ the three-dimensional compact Abelian Higgs model has a phase-transition line ${\ensuremath{\beta}}_{\mathrm{c}}(\ensuremath{\kappa})$ which is first order for $\ensuremath{\kappa}$ below a finite tricritical value ${\ensuremath{\kappa}}_{\mathrm{tri}}$ and second order above. The $\ensuremath{\beta}=\ensuremath{\infty}$ first order phase transition persists for finite $\ensuremath{\beta}$ and joins the second order phase transition at a tricritical point $({\ensuremath{\beta}}_{\mathrm{tri}},{\ensuremath{\kappa}}_{\mathrm{tri}})=(1.23\ifmmode\pm\else\textpm\fi{}0.03,1.73\ifmmode\pm\else\textpm\fi{}0.03).$ For all other integer $q>~2$ we have considered, the entire phase-transition line ${\ensuremath{\beta}}_{c}(\ensuremath{\kappa})$ is critical. We have used finite-size scaling of the second and third moments of the action to extract critical exponents $\ensuremath{\alpha}$ and $\ensuremath{\nu}$ without invoking hyperscaling, for the $\mathrm{XY}$ model, the ${Z}_{2}$ spin and lattice gauge models, as well as the compact Abelian Higgs model for $q=2$ and $q=3.$ In all cases, we have found that for practical system sizes, the third moment gives scaling of superior quality compared to the second moment. We have also computed the exponent ratio for the $q=2$ compact $U(1)$ Higgs model along the critical line, finding a continuously varying ratio $(1+\ensuremath{\alpha})/\ensuremath{\nu},$ as well as continuously varying $\ensuremath{\alpha}$ and $\ensuremath{\nu}$ as $\ensuremath{\kappa}$ is increased from $0.76$ to $\ensuremath{\infty},$ with the Ising universality class $(1+\ensuremath{\alpha})/\ensuremath{\nu}=1.763$ as a limiting case for $\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\beta}}\ensuremath{\infty},\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\kappa}}0.761,$ and the $\mathrm{XY}$ universality class $(1+\ensuremath{\alpha})/\ensuremath{\nu}=1.467$ as a limiting case for $\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\beta}}0.454,\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\kappa}}\ensuremath{\infty}.$ However, the critical line exhibits a remarkable resilience of ${Z}_{2}$ criticality as $\ensuremath{\beta}$ is reduced along the critical line. Thus, the three-dimensional compact Abelian Higgs model for $q=2$ appears to represent a fixed-line theory defining a new universality class. We relate these results to a recent microscopic description of zero-temperature quantum phase transitions within insulating phases of strongly correlated systems in two spatial dimensions, proposing the above to be the universality class of the zero-temperature quantum phase transition from a Mott-Hubbard insulator to a charge-fractionalized insulator in two spatial dimensions, which thus is that of the 3D Ising model for a considerable range of parameters.