The present paper is dedicated to the study of the discrepancies encountered in electron temperatures (Te) measurements carried out with Electron Cyclotron Emission (ECE) and Thomson Scattering (TS) diagnostics in the core of the JET tokamak. A large database of discharges has been collected, including high-performance scenarios performed with Deuterium only and Deuterium-Tritium mixtures. Discrepancies have been found between core Te measurements taken with an X-mode ECE interferometer (TECE) and a LIDAR TS system (TLID) for Te > 5 keV. Depending on the plasma scenario, TECE has been found to be systematically higher or lower than TLID. Discrepancies have also been observed between the peaks of the ECE spectrum in the second (X2) and third (X3) harmonic domains, even in high optical thickness conditions. These discrepancies can be interpreted as evidence of the presence of non-Maxwellian features in the electron energy distribution function (EEDF). In order to investigate the relation between the shape of the EEDF and the measured discrepancies, a model for bipolar perturbations of Maxwellian EEDF has been developed. The model allows analytical calculations of ECE absorption and emission coefficients, hence the comparison of modelled ECE spectra with experimental data. The different experimental results observed for the various JET scenarios have been found to be qualitatively reproducible by adapting the model parameters, suggesting that bipolar distortions of the bulk EEDF could play a role in giving rise to the reported discrepancies between ECE and TS measurements.