In view that most Global Navigation Satellite System (GNSS) users are still using single-frequency receivers due to the low costs, single-frequency Precise Point Positioning (PPP) has been attracting increasing attention in the GNSS community. For a long period, single-frequency PPP technology has mainly relied on the Global Positioning System (GPS). With the recent revitalisation of the Russian GLONASS constellation and two newly emerging constellations, BeiDou and Galileo, it is now feasible to investigate the performance of Four-Constellation integrated Single-Frequency PPP (FCSF-PPP) with GPS, GLONASS, BeiDou and Galileo measurements. In this study, a FCSF-PPP model is presented to simultaneously process observations from all four GNSS constellations. Datasets collected at 47 globally distributed four-system Multi-GNSS Experiment (MGEX) stations on seven consecutive days and a kinematic experimental dataset are employed to fully assess the performance of FCSF-PPP. The FCSF-PPP solutions are compared to GPS-only and combined GPS/GLONASS single-frequency PPP solutions. The results indicate that the positioning performance is significantly improved by integrating multi-constellation signals.