We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

We present a 30h integration Very Large Telescope X-shooter spectrum of the Lyman series forest towards the $z = 7.084$ quasar ULAS J1120+0641. The only detected transmission at ${\rm S/N}>5$ is confined to seven narrow spikes in the Ly$\alpha$ forest, over the redshift range $5.858<z<6.122$, just longward of the wavelength of the onset of the Ly$\beta$ forest. There is also a possible detection of one further unresolved spike in the Ly$\beta$ forest at $z=6.854$, with ${\rm S/N}=4.5$. We also present revised Hubble Space Telescope F814W photometry of the source. The summed flux from the transmission spikes is in agreement with the F814W photometry, so all the transmission in the Lyman series forest may have been detected. There is a Gunn-Peterson (GP) trough in the Ly$\alpha$ forest from $z=6.122$ all the way to the quasar near zone at $z=7.04$. The trough, of comoving length $240\,h^{-1}$Mpc, is over twice as long as the next longest known GP trough. We combine the spectroscopic and photometric results to constrain the evolution of the Ly$\alpha$ effective optical depth with redshift, $\tau_{\rm GP}^{\rm eff}$ extending a similar analysis by Simpson et al. We find $\tau_{\rm GP}^{\rm eff} \propto (1+z)^{\xi}$ where $\xi = 11.2^{+0.4}_{-0.6}$, for $z > 5.5$. The data nevertheless provide only a weak limit on the volume-weighted hydrogen intergalactic (IGM) neutral fraction at $z\sim 6.5$, $x_{{\rm HI}} > 10^{-4}$, similar to limits at redshift $z\sim6$ from less distant quasars. The new observations cannot extend measurements of the neutral fraction of the IGM to higher values because absorption in the Ly$\alpha$ forest is already saturated near $z\sim6$. For higher neutral fractions, other methods such as measuring the red damping wing of the IGM will be required.

We present new observations of the highest-redshift quasar known, ULAS J1120+0641, redshift $z=7.084$, obtained in the optical, at near-, mid-, and far-infrared wavelengths, and in the sub-mm. We combine these results with published X-ray and radio observations to create the multiwavelength spectral energy distribution (SED), with the goals of measuring the bolometric luminosity $L_{\rm bol}$, and quantifying the respective contributions from the AGN and star formation. We find three components are needed to fit the data over the wavelength range $0.12-1000\,\mu$m: the unobscured quasar accretion disk and broad-line region, a dusty clumpy AGN torus, and a cool 47K modified black body to characterise star formation. Despite the low signal-to-noise ratio of the new long-wavelength data, the normalisation of any dusty torus model is constrained within $\pm40\%$. We measure a bolometric luminosity $L_{\rm bol}=2.6\pm0.6\times10^{47}\,$erg$\,$s$^{-1}=6.7 \pm 1.6\times10^{13}L_{\odot}$, to which the three components contribute $31\%,32\%,3\%$, respectively, with the remainder provided by the extreme UV $<0.12\,\mu$m. We tabulate the best-fit model SED. We use local scaling relations to estimate a star formation rate (SFR) in the range $60-270\,{\rm M}_\odot$/yr from the [C$\,{\scriptsize \rm II}$] line luminosity and the $158\,\mu$m continuum luminosity. An analysis of the equivalent widths of the [C$\,{\scriptsize \rm II}$] line in a sample of $z>5.7$ quasars suggests that these indicators are promising tools for estimating the SFR in high-redshift quasars in general. At the time observed the black hole was growing in mass more than 100 times faster than the stellar bulge, relative to the mass ratio measured in the local universe, i.e. compared to ${M_{\rm BH}}/{M_{\rm bulge}} \simeq 1.4\times10^{-3}$, for ULAS J1120+0641 we measure ${\dot{M}_{\rm BH}}/{\dot{M}_{\rm bulge}} \simeq 0.2$.