In this paper, we consider the following Choquard equation \(\begin{aligned} -\varepsilon ^{N}\Delta _{N}u+V(x)|u|^{N-2}u=\varepsilon ^{\mu -N}\left( I_\mu *F(u)\right) f(u) \quad {\text{ in }\quad \mathbb {R}^N}, \end{aligned}\) where \(N\ge 3\) , \(I_\mu =|x|^{-\mu }\) with \(0<\mu <N\) , \(\Delta _{N}u=\textrm{div}(|\nabla u|^{N-2}\nabla u)\) denotes the N-Laplacian operator, V(x) is a continuous real function on \(\mathbb {R}^N\) , F(s) is the primitive of f(s) and \(\varepsilon \) is a positive parameter. Assuming that the nonlinearity f(s) has critical exponential growth in the sense of Trudinger–Moser inequality, we establish the existence, multiplicity and concentration of solutions by variational methods and Ljusternik–Schnirelmann theory, which extends the works of Alves and Figueiredo (J Differ Equ 246:1288–1311, 2009) to the problem with Choquard nonlinearity, Alves et al. (J Differ Equ 261:1933–1972, 2016) to higher dimension.