We present a sample of 38 intervening damped Lyman {alpha} (DLA) systems identified towards 100 z>3.5 quasars, observed during the XQ-100 survey. The XQ-100 DLA sample is combined with major DLA surveys in the literature. The final combined sample consists of 742 DLAs over a redshift range approximately 1.6<z_abs_<5.0. We develop a novel technique for computing {Omega}_HI_^DLA^ as a continuous function of redshift, and we thoroughly assess and quantify the sources of error therein, including fitting errors and incomplete sampling of the high column density end of the column density distribution function. There is a statistically significant redshift evolution in {Omega}_HI_^DLA^ (>=3{sigma}) from z~2 to z~5. In order to make a complete assessment of the redshift evolution of {Omega}_HI_, we combine our high-redshift DLA sample with absorption surveys at intermediate redshift and 21-cm emission line surveys of the local universe. Although {Omega}_HI_^DLA^, and hence its redshift evolution, remains uncertain in the intermediate-redshift regime (0.1<z_abs_<1.6), we find that the combination of high-redshift data with 21-cm surveys of the local universe all yield a statistically significant evolution in {Omega}_HI_ from z~0 to z~5 (>=3{sigma}). Despite its statistical significance, the magnitude of the evolution is small: a linear regression fit between {Omega}_HI_ and z yields a typical slope of ~0.17x10^-3^, corresponding to a factor of ~4 decrease in {Omega}_HI_ between z=5 and z=0.