The physics of recombination lines (RLs) in the HeI singlet system is expected to be rela- tively simple, supported by accurate atomic models. We examine the intensities of HeI singlets λλ3614, 3965, 5016, 6678, 7281 and the triplet He I λ5876 in various types of ionized nebulae and compare them with theoretical predictions to test the validity of the “Case B” recombination scenario and the assumption of thermal homogeneity. Our analysis includes 85 spectra from Galactic and extragalactic HII regions, 90 from star-forming galaxies, and 218 planetary nebulae, all compiled by the DEep Spectra of Ionized REgions Database Extended (DESIRED-E) project. By evaluating the ratios He I λ7281/λ6678 and He I λ7281/λ5876, we determine Te (He I) and compare it with direct measurements of Te([O III] λ4363/λ5007). We find that Te(He I) is systematically lower than Te([O III]) across most objects and nebula types. Additionally, we identify a correlation between the abundance discrepancy factor (ADF(O2+)) and the difference Te([OIII]) - Te(HeI) for planetary nebulae. We explore two potential explanations: photon loss from n^1P → 11S transitions and temperature inhomogeneities. Deviations from “Case B” may indicate photon absorption by H I rather than He I and/or generalized ionizing photon escape, highlighting the need for detailed consideration of radiative transfer effects. If temperature inhomogeneities are widespread, identifying a common physical phenomenon affecting all ionized nebulae is crucial. Our results suggest that both scenarios can contribute to the observed discrepancies.