Recent measurements of various charm-hadron ratios in $pp$, $p$-Pb and Pb-Pb collisions at the LHC have posed challenges to the theoretical understanding of heavy-quark hadronization. The $Lambda_c/D^0$ ratio in $pp$ and $p$-Pb collisions shows larger values than those found in $e^+e^-$ and $ep$ collisions and predicted by Monte-Carlo event generators based on string fragmentation, at both low and intermediate transverse momenta ($p_T$). In AA collisions, the $D_s/D$ ratio is significantly enhanced over its values in $pp$, while the $Lambda_c/D^0$ data indicates a further enhancement at intermediate $p_T$. Here, we report on our recent developments for a comprehensive description of the charm hadrochemistry and transport in $pp$ and $AA$ collisions. For $pp$ collisions we find that the discrepancy between the $Lambda_c/D^0$ data and model predictions is much reduced by using a statistical hadronization model augmented by a large set of missing states in the charm-baryon spectrum, contributing to the $Lambda_c$ via decay feeddown. For $AA$ collisions, we develop a 4-momentum conserving resonance recombination model for charm-baryon formation implemented via event-by-event simulations that account for space-momentum correlations (SMCs) in transported charm- and thermal light-quark distributions. The SMCs, together with the augmented charm-baryon states, are found to play an important role in describing the baryon-to-meson enhancement at intermediate momenta. We emphasize the importance of satisfying the correct (relative) chemical equilibrium limit when computing the charm hadrochemistry and its momentum dependence with coalescence models.