We consider the approximate recovery of multivariate periodic functions from a discrete set of function values taken on a rank-s integration lattice. The main result is the fact that any (non-)linear reconstruction algorithm taking function values on a rank-s lattice of size M has a dimension-independent lower bound of 2 −(α+1)/2 M −α/2 when considering the optimal worst-case error with respect to function spaces of (hybrid) mixed smoothness α > 0 on the d-torus. We complement this lower bound with upper bounds that coincide up to logarithmic terms. These upper bounds are obtained by a detailed analysis of a rank-1 lattice sampling strategy, where the rank-1 lattices are constructed by a componentby-component (CBC) method. This improves on earlier results obtained in [25] and [27]. The lattice (group) structure allows for an efficient approximation of the underlying function from its sampled values using a single one-dimensional fast Fourier transform. This is one reason why these algorithms keep attracting significant interest. We compare our results to recent (almost) optimal methods based upon samples on sparse grids. This paper deals with the reconstruction of multivariate periodic functions from a discrete set of M function values along rank-1 lattices. Such lattices are widely used for the efficient numerical integration of multivariate periodic functions since the 1950ies [1,21,29,35,6] and represent a well-distributed set of points in [0, 1) d . A rank-1 lattice with M ∈ N points and generating vector z ∈ Z d is given byIn this paper we will show that restricting the set of available discrete information to samples from a rank-s lattice, cf.[35], seriously affects the rate of convergence of a corresponding worst-case error with respect to classes of functions with (hybrid) mixed smoothness α > 0.To be more precise, for any (possibly nonlinear) reconstruction procedure from sampled values along rank-s lattices we can find a function in the periodic Sobolev spaces H α mix such that the L 2 (T d ) mean square error is at least 2 −(α+1)/2 M −α/2 . In contrast to that, it has been proved recently, that the sampling recovery from (energy) sparse grids leads to much better convergence rates, namely M −α in the main term, see [33,41,4].Subsequently, we study particular reconstructing algorithms, which are based on the naive approach of approximating the potentially "largest" Fourier coefficients (integrals) with the same rank-1 lattice rule. Despite the lacking asymptotical optimality, recovery from so-called reconstructing rank-1 lattices, cf. [15,18], has some striking advantages.First, the matrix of the underlying linear system of equations has orthogonal columns due to the group structure [2] and the reconstructing property of the used rank-1 lattices. Consequently, the computation is stable, cf. [17,15].Second, the CBC strategy [14, Tab. 3.1] provides a search method for a reconstructing rank-1 lattice which allows for the computation of the approximate Fourier coefficients belonging to frequencies lying on ...