More than 30 % of $gamma$-ray sources detected in the last source catalog of the Fermi satellite have no observational counterpart at other frequencies. A significant fraction of these sources is positionally in agreement with star-forming regions dominated by associations of T Tauri stars. Rho Ophiuchi, which is one of the closest star-forming regions, matches the unidentified Fermi source 4FGL J1625.3-2338. In this work we modeled the spectral energy distribution considering some dominant radiative processes in T Tauri stars. Accounting for a total of 22 Class III T Tauri stars in the region, integrated $gamma$-ray luminosity in the 100 MeV to 100 GeV energy range is consistent with the observed in the catalog for the 4FGL J1625.3-2338 source.