Assessing multivariate gene-metabolome associations with rare variants using Bayesian reduced rank regression