We investigate multivariate subordination of Lévy processes which was first introduced by Barndorff-Nielsen et al. [O.E. Barndorff-Nielsen, F.E. Benth, and A. Veraart, Modelling electricity forward markets by ambit fields, J. Adv. Appl. Probab. (2010)], in a Hilbert space valued setting which has been introduced in Pérez-Abreu and Rocha-Arteaga [V. Pérez-Abreu and A. Rocha-Arteaga, Covariance-parameter Lévy processes in the space of trace-class operators, Infin. Dimens. Anal. Quantum Probab. Related Top. 8(1) (2005), pp. 33–54]. The processes are explicitly characterized and conditions for integrability and martingale properties are derived under various assumptions of the Lévy process and subordinator. As an application of our theory we construct explicitly some Hilbert space valued versions of Lévy processes which are popular in the univariate and multivariate case. In particular, we define a normal inverse Gaussian Lévy process in Hilbert space. The resulting process has the property that at each time all its finite dimensional projections are multivariate normal inverse Gaussian distributed as introduced in Rydberg [T. Rydberg, The normal inverse Gaussian Lévy process: Simulation and approximation, Commun. Stat. Stochastic Models 13 (1997), pp. 887–910].