Simultaneous Heterogeneity and Reduced-rank Learning for Multivariate Response Regression