We have developed an efficient method for performing density matrix renormalization group (DMRG) simulations of the SU() Fermi-Hubbard chain, fully leveraging the SU() symmetry of the problem. This method extends a previously developed approach for the SU() Heisenberg model and relies on the systematic use of the semistandard Young tableaux (SSYT) basis in a DMRG algorithm “à la White.” Specifically, the method aligns the site-by-site growth process of the infinite-size part of the DMRG, in its original formulation, with the site-by-site construction of the SSYT (or Gelfand-like) basis, based on the chain of unitary subgroups . We give special emphasis to the calculation of the symmetry-resolved reduced matrix elements of the hopping terms between the left and the right block, which makes direct use of the basis of SSYT and of the Gelfand-Tsetlin coefficients, offering a computational advantage in scaling with compared to alternative methods that rely on summing over Clebsch-Gordan coefficients. Focusing on the model with homogeneous hopping between nearest neighbors, we have calculated the ground-state energy as a function of , i.e., the atom-atom interaction amplitude, up to for filling (one particle per site in average), and for one atom (respectively, hole) away from filling . It allows us to compute the charge gaps and to estimate in the thermodynamical limit, the critical value (separating the Mott insulator from the metallic phase), which is proved to be finite (i.e., ) and to increase with . Central charges are also extracted from the entanglement entropy using the Calabrese-Cardy formula and are consistent with the theoretical predictions: , expected from the Wess-Zumino-Witten conformal field theory in the spin sector for the Mott phase, and in the metallic phase, reflecting the presence of one additional (charge) gapless critical mode.