I have a short RNAseq time series dataset (two time pointsl T0 and TN) and I have two conditions Undrugged and Drugged. For each of the drug conditions the timepoints are paired (i.e one T0 corresponds to one T1). There are two T0/TN-paired replicates for each drug condition. So we have (Drugged_T0_A, Drugged_T1_A), (Drugged_T0_B, Drugged_T1_B), (Undrugged_T0_C, Undrugged_T1_C), (Undrugged_T0_D and Undrugged_T1_D). Is it possible to set up the test to identify genes that are altered more in the Drugged state (going from T0 to T1) than in the Undrugged state? I'm pretty sure this can be done with DESeqs glm functionality. Essentially I want to find genes where (keeping the paired structure intact) log(Drugged_TN/Drugged_T0) << log(Undrugged_TN/Undrugged_T0) or log(Drugged_TN/Drugged_T0) >> log(Undrugged_TN/Undrugged_T0). Thanks So Much Simon