Neuroblastoma (NB) is the most commonly occurring solid tumor in children. The disease usually arises in the adrenal medulla, and it is characterized by a remarkable heterogeneity in its progression. Most NB patients with an advanced disease have massive bone marrow infiltration at diagnosis. Lung metastasis represents a widely disseminated stage and is typically considered to be a terminal event. Much like other malignancies, NB progression is a complex, multistep process. The expression, function, and significance of the various factors involved in NB progression must be studied in relevant in vivo and in vitro models. Currently, models consisting of metastatic and nonmetastatic cell variants of the same genetic background exist for several types of cancer; however, none exists for NB. In the present study, we describe the generation of a NB metastasis model. SH-SY5Y and MHH-NB-11 NB cells were inoculated orthotopically into the adrenal glands of athymic nude mice. Neuroblastoma cells metastasizing to the lungs were isolated from mice bearing adrenal tumors. Lung metastatic variants were generated by repeated cycles of in vivo passage. Characterization of these variants included cellular morphology and immunophenotyping in vitro, aggressiveness in vivo, and various biologic parameters in vitro. The NB metastatic variant in each model displayed unique properties, and both metastatic variants demonstrated a metastatic phenotype in vivo. These reproducible models of human NB metastasis will serve as an unlimited source of transcriptomic and proteomic material. Such models can facilitate future studies on NB metastasis and the identification of novel NB biomarkers and targets for therapy.