DEseq2 was also used to calculate fold changes and p-values and perform optional covariate correction.