This paper is dedicated to the optimisation of Non-Binary LDPC codes when associated to high-order modulations. To be specific, we propose to specify the values of the non-zero NB-LDPC parity matrix coefficients depending on the corresponding check node equation and the Euclidean distance of the coded modulation. In other words, we explore the joint optimisation of the modulation mapping and the non-binary matrix. The performance gains announced by a theoretical analysis based on the Union Bound are confirmed by simulations results. We obtain an 0.2-dB gain in the high SNR regime compared to other stateof-the-art matrices.