A systematic study of the change in proton single-particle occupancies in two neutrinoless double beta decay candidates, 130Te to 130Xe and 136Xe to 136Ba, has been performed.Final states in 129Sb, 129I, 135I and 135Cs have been populated using the (d,3He) single proton removal reaction. The deuterons were accelerated to 101 MeV using the coupled azimuthally varying field and Ring cyclotrons at the Research Center for Nuclear Physics, University of Osaka. The outgoing ejectiles were momentum analysed using the Grand Raiden magnetic spectrometer. Absolute cross sections were measured for states up to 7 MeV in excitation. Transferred angular momenta were identified through a comparison of angular distributions with those calculated using the distorted-wave Born approximation. Spectroscopic factors were extracted from the experimental cross sections. In addition to the 0v2B decay candidates the reaction was also performed on 128Te, 132Xe, 134Xe and 138Ba as a consistency check. The occupancies of the nuclei were determined from the spectroscopic factors, the reaction model was normalised using a common normalisation factor across all targets. The change in occupancies between the 0v2B decay candidates and their daughters were then determined. The change in occupancies were then compared with those calculated by interacting shell model, interacting boson model and quasiparticle random phase approximation. This comparison showed that whilst the three theories were qualitatively able to reproduce the change in occupancies, quantitatively there are significant discrepancies. These are the same models that are used to determine the nuclear matrix elements governing the rate of 0v2B decay.