We use the path-integral Monte Carlo (PIMC) methodology and state-of-the-art two-body and three-body potentials to calculate the fourth virial coefficients D(T) of 4He and 3He as features