In this article, a characterization of the class of Herglotz-Nevanlinna functions in several variables is given in terms of an integral representation. The conditions on the representing measure are discussed in detail, and, furthermore, the properties of the symmetric extension of a Herglotz-Nevanlinna function are also investigated.