Special case of sunset: reduction and epsilon-expansion

We consider two loop sunset diagrams with two mass scales m and M at the threshold and pseudotreshold that cannot be treated by earlier published formula. The complete reduction to master integrals is given. The master integrals are evaluated as series in ratio m/M and in epsilon with the help of differential equation method. The rules of asymptotic expansion in the case when q^2 is at the (pseudo)threshold are given.