Matrix variate skew normal distributions