A theoretical description of femtosecond two-dimensional electronic spectroscopy of multichromophoric systems is presented. Applying the stationary phase approximation to the calculation of photon echo spectra and taking into account exciton relaxation processes, we obtain an analytic expression for numerical simulations of time- and frequency-resolved 2D photon echo signals. The delocalization of one-exciton states, spatial overlaps between the probability densities of different excitonic states, and their influences on both one- and two-dimensional electronic spectra are studied. The nature of the off-diagonal cross-peaks and the time evolution of both diagonal and off-diagonal peak amplitudes are discussed in detail by comparing experimentally measured and theoretically simulated 2D spectra of the natural Fenna-Matthews-Olson (FMO) photosynthetic light-harvesting complex. We find that there are two noncascading exciton energy relaxation pathways.