Following our previous spectroscopic observations of $z>7$ galaxies with Gemini/GNIRS and VLT/XSHOOTER, which targeted a total of 8 objects, we present here our results from a deeper and larger VLT/FORS2 spectroscopic sample of Wide Field Camera 3 selected $z>7$ candidate galaxies. With our FORS2 setup we cover the 737-1070nm wavelength range, enabling a search for Lyman-$\alpha$ in the redshift range spanning 5.06 - 7.80. We target 22 $z$-band dropouts and find no evidence of Lyman-$\alpha$ emission, with the exception of a tentative detection ($<5\sigma$, which is our adopted criterion for a secure detection) for one object. The upper limits on Lyman-$\alpha$ flux and the broad-band magnitudes are used to constrain the rest-frame Equivalent Widths for this line emission. We analyse our FORS2 observations in combination with our previous GNIRS and XSHOOTER observations, and suggest that a simple model where the fraction of high rest-frame Equivalent Width emitters follows the trend seen at $z=3-6.5$ is inconsistent with our non-detections at $z\sim7.8$ at the 96\% confidence level. This may indicate that a significant neutral HI fraction in the intergalactic medium suppresses Lyman-$\alpha$, with an estimated neutral fraction $\chi_{HI}\sim0.5$, in agreement with other estimates.

By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts.