Studies of young stellar objects (YSOs) in the Galaxy have found that a significant fraction exhibit photometric variability. However, no systematic investigation has been conducted on the variability of extragalactic YSOs. Here we present the first variability study of massive YSOs in a $\sim 1.5\,\mathrm{deg^2}$ region of the Large Magellanic Cloud (LMC). The aim is to investigate whether the different environmental conditions in the metal-poor LMC ($\sim$ 0.4-0.5 Z_sun) have an impact on the variability characteristics. Multi-epoch near-infrared (NIR) photometry was obtained from the VISTA Survey of the Magellanic Clouds (VMC) and our own monitoring campaign using the VISTA telescope. By applying a reduced $\chi^2$-analysis, stellar variability was identified. We found 3062 candidate variable stars from a population of 362 425 stars detected. Based on several Spitzer studies, we compiled a sample of high-reliability massive YSOs: a total of 173 massive YSOs have NIR counterparts ($K_{\mathrm{s}}\sim 18.5\,$mag) in the VMC catalogue, of which 39 display significant ($>3\sigma$) variability. They have been classified as eruptive, fader, dipper, short-term variable and long-period variable YSOs based mostly on the appearance of their $K_{\mathrm{s}}$ band light curves. The majority of YSOs are aperiodic, only five YSOs exhibit periodic lightcurves. The observed amplitudes are comparable or smaller than those for Galactic YSOs (only two Magellanic YSOs exhibit $\Delta K_{\mathrm{s}}>1\,$mag), not what would have been expected from the typically larger mass accretion rates observed in the Magellanic Clouds.

Detailed studies of intermediate/low mass pre-main sequence (PMS) stars outside the Galaxy have so far been conducted only for small targeted regions harbouring known star formation complexes. The VISTA Survey of the Magellanic Clouds (VMC) provides an opportunity to study PMS populations down to solar masses on a galaxy-wide scale. Our goal is to use near-infrared data from the VMC survey to identify and characterise PMS populations down to ~1 M_sun across the Magellanic Clouds. We present our colour-magnitude diagram method, and apply it to a ~1.5 deg^2 pilot field located in the Large Magellanic Cloud. The pilot field is divided into equally-sized grid elements. We compare the stellar population in every element with the population in nearby control fields by creating K_s/(Y-K_s) Hess diagrams; the observed density excesses over the local field population are used to classify the stellar populations. Our analysis recovers all known star formation complexes in this pilot field (N44, N51, N148 and N138) and for the first time reveals their true spatial extent. In total, around 2260 PMS candidates with ages $\lesssim$ 10 Myr are found in the pilot field. PMS structures, identified as areas with a significant density excess of PMS candidates, display a power-law distribution of the number of members with a slope of -0.86+-0.12. We find a clustering of the young stellar populations along ridges and filaments where dust emission in the far-infrared (FIR) (70 micron - 500 micron) is bright. Regions with young populations lacking massive stars show a lesser degree of clustering and are usually located in the outskirts of the star formation complexes. At short FIR wavelengths (70 micron, 100 micron) we report a strong dust emission increase in regions hosting young massive stars, which is less pronounced in regions populated only by less massive ($\lesssim$ 4 M_sun) PMS stars.