Akemann and Strahov recently introduced a product of two coupled Gaussian random matrices with the identity coupling matrix and showed that the singular values are distributed as a determinantal point process(cf. arXiv:1504.02047). We turn to the product with a general coupling matrix A and obtain a similar determinantal point process. Interestingly, a phase transition phenomenon for the least singular values is observed as the coupling matrix varies in four different regimes.