Characterizing the properties and the evolution of the first stars and galaxies is a challenging task for traditional galaxy surveys since they are sensitivity limited and can only detect the brightest light sources. Three-dimensional intensity mapping (IM) of transition lines can be a valuable alternative to study the high redshift Universe given that this technique avoids sensitivity limitation problems by measuring the overall emission of a line, with a low resolution, without resolving its sources. While 21cm line IM surveys probe neutral hydrogen gas and can, therefore, be used to probe the state of the IGM and the evolution of the ionization field during the Epoch of Reionization (EoR). IM surveys of other lines, such as CO, CII, Ly-alpha or H-alpha, can be used to probe the galaxies which emitted most of the ionizing radiation responsible for the EoR. These lines will trace the different ISM gas phases, the excitation state of this gas, its metallicity, etc. This study addresses IM of multiple transition lines and how it can be used to probe the EoR and to constrain the redshift evolution of galaxy properties.