The three Minkowski functionals and the recently defined length of the skeleton are estimated for the co-added first-year Wilkinson Microwave Anisotropy Probe (WMAP) data and compared with 5000 Monte Carlo simulations, based on Gaussian fluctuations with the a priori best-fit running-index power spectrum and WMAP-like beam and noise properties. Several power spectrum-dependent quantities, such as the number of stationary points, the total length of the skeleton, and a spectral parameter, , are also estimated. While the area and length Minkowski functionals and the length of the skeleton show no evidence for departures from the Gaussian hypothesis, the northern hemisphere genus has a 2 that is large at the 95% level for all scales. For the particular smoothing scale of 3N40 FWHM it is larger than that found in 99.5% of the simulations. In addition, the WMAP genus for negative thresholds in the northern hemisphere has an amplitude that is larger than in the simulations with a significance of more than 3 . On the smallest angular scales considered, the number of extrema in the WMAP data is high at the 3 level. However, this can probably be attributed to the effect of point sources. Finally, the spectral parameter is high at the 99% level in the northern Galactic hemisphere, while perfectly acceptable in the southern hemisphere. The results provide strong evidence for the presence of both nonGaussian behavior and an unexpected power asymmetry between the northern and southern hemispheres in the WMAP data. Subject headingg s: cosmic microwave background -cosmology: observations -methods: statistical