It is shown that the parametric spectral statistics in the critical random matrix ensemble with multifractal eigenvector statistics are identical to the statistics of correlated 1D fermions at finite temperatures. For weak multifractality the effective temperature of fictitious 1D fermions is proportional to T ef f ∝ (1 − dn)/n ≪ 1, where dn is the fractal dimension found from the n-th moment of inverse participation ratio. For large energy and parameter separations the fictitious fermions are described by the Luttinger liquid model which follows from the Calogero-Sutherland model. The low-temperature asymptotic form of the two-point equal-parameter spectral correlation function is found for all energy separations and its relevance for the low temperature equal-time density correlations in the Calogero-Sutherland model is conjectured. PACS number(s): 72.15. Rn, 72.70.+m, 72.20.Ht, The spectral statistics in complex quantum systems are signatures of the underlying dynamics of the corresponding classical counterpart. The spectral statistics in chaotic and disordered systems in the limit of infinite dimensionless conductance g is described by the classical random matrix theory of Wigner and Dyson [1] (WD statistics). The WD statistics possess a remarkable property of universality: it depends only on the symmetry class with respect to the time-reversal transformation T . The three symmetry classes correspond to the lack of T -invariance (the unitary ensemble, β = 2); the Tinvariant systems with T 2 = 1 (the orthogonal ensemble, β = 1), and the T -invariant systems with T 2 = −1 (the symplectic ensemble, β = 4), respectively. The physical ground behind this universality is the structureless eigenfunctions in the ergodic regime which implies the invariance of the eigenfunction statistics with respect to a unitary transformation of the basis.In real disordered metals the eigenfunctions are not basis-invariant. The basis-preference reaches its extreme form for the strongly impure metals where all eigenfunctions are localized in the coordinate space but delocalized in the momentum space. In this case the spectral statistics is Poissonian in the thermodynamic (TD) limit.For low-dimensional systems d = 1, 2, where all states are localized in the TD limit, one can observe the smooth crossover from the WD to the Poisson spectral statistics as a function of the parameter ξ/L, where ξ is the localization radius and L is the system size. The dependence of the spectral correlation functions on the energy variable s = E/∆ (∆ is the mean level separation) is non-universal for finite L/ξ but all of them tend to the Poisson limit as L/ξ → ∞.In systems of higher dimensionality d > 2 the situation is different because of the presence of the Anderson localization transition at a critical disorder W = W c . In the metal phase W < W c the dimensional conductance g(L) → ∞ as L → ∞ and one obtains the WD spectral statistics in the TD limit. In the insulator state g(L) → 0 at L → ∞ and the limiting statistics is Poissonian. However, there i...