Does log-transforming a slightly-skewed y variable make any sense? I'm separately analyzing 5 subsets of my data and running multiple linear regressions with the same outcome variable in each. (I am also running a single regression with all the subsets and using the subset variable as a predictor, but it is useful to run separate regressions for each subset because the results are quite different for each subset, in ways you would not be able to tell from the overall regression.)
The outcome variable is continuous and, prior to the log-transformation, I was modeling it in its original form. It is usually positive but in some cases zero (for that reason, I add a constant when log-transforming it.)
The variable is highly skewed in 4 subsets so I'm log-transforming it with a constant, chosen separately for each subset to lower the skewness of the variable in that subset to acceptable levels.
The thing is, the variable in the 5th subset is only slightly skewed (.28), well within the acceptable range. However, if I have a single table illustrating all 5 regressions, it's going to be hard for any reader to compare the magnitude of coefficients across the different subsets if 4 of the subsets have a logged y variable and one of them doesn't.
So my question is, does it make any sense to log-transform the y variable for the slightly skewed subset too? Would this make the results more accurate?
I tried doing so, adding a constant to reduce the skewness all the way to .006. When I run the regression with the log-transformed y variable, the results are mainly the same, but now two additional x variables have statistically significant effects. So it really does change the overall results.
On the one hand, presenting my results this way would be easier to interpret for the reader, and now the outcome variable is utterly devoid of skewness. But on the other hand, it wasn't necessary to log-transform that y variable for that subset. Is there a risk I actually made my results less accurate? Is there any other reason why this is a bad idea?

A separate model for each subset, each with a separate transformation of the same outcome variable, isn't a good approach. Your worry

...if I have a single table illustrating all 5 regressions, it's going to be hard for any reader to compare the magnitude of coefficients across the different subsets if 4 of the subsets have a logged y variable and one of them doesn't

holds even "with a constant, chosen separately for each subset" added prior to log transformation. A regression coefficient associated with an outcome of $\log(y+1)$ isn't on the same scale as one associated with an outcome of $\log(y+10)$.
The "results ... quite different for each subset" can be handled by going beyond a simple coefficient for subset to include interactions between other predictors and subset. In a single combined model you could examine the overall distribution of residuals and then decide how to proceed. As you note, it's the distribution of residuals that matters, not that of the original y values.
If you want to go the y-transformation route, choose a common transformation that provides an acceptable overall distribution in the combined model. That could be a single choice of constant prior to log transformation. Or it could be a different type of transformation entirely, like a square root transformation. You don't have to end up with an exactly normal distribution of residuals; they just have to be well enough behaved. You can also get robust estimators of coefficient (co)variances that can handle deviations from the simplest linear regression assumptions.
Perhaps better, avoid y-transformation entirely. Continuous outcomes can be modeled with ordinal regression. A proportional-odds model is an extension of standard non-parametric tests to regression models. Frank Harrell explains this in Chapter 13 of his course notes and book, and provides an orm() function in his rms package that's designed to handle continuous outcomes efficiently.
