We apply the notion of discrete supersymmetry based on matrix factorization to quantum systems consisting of coupled bosonic oscillators to construct isospectral bosonic quantum networks. By using the algebra that arises due to the indistinguishabili