We study the three-dimensional Anderson model of localization with anisotropic hopping, i.e., weakly coupled chains and weakly coupled planes. In our extensive numerical study we identify and characterize the metal-insulator transition using energy-level statistics. The values of the critical disorder W c are consistent with results of previous studies, including the transfermatrix method and multifractal analysis of the wave functions. W c decreases from its isotropic value with a power law as a function of anisotropy. Using high accuracy data for large system sizes we estimate the critical exponent ν = 1.45 ± 0.2. This is in agreement with its value in the isotropic case and in other models of the orthogonal universality class. The critical level statistics which is independent of the system size at the transition changes from its isotropic form towards the Poisson statistics with increasing anisotropy. 71.30.+h, 72.15.Rn, 73.20.Dx