A systems biology approach was applied to investigate the mechanisms of chromosomal instability in melanoma cell lines. Chromosomal instability was quantified using array comparative genomic hybridization to identify somatic copy number alterations (deletions and duplications). Primary human melanocytes displayed an average of 8.5 alterations per cell primarily representing known polymorphisms. Melanoma cell lines displayed 25 to 131 alterations per cell, with an average of 68, indicative of chromosomal instability. Copy number alterations included approximately equal numbers of deletions and duplications with greater numbers of hemizygous (−1,+1) alterations than homozygous (−2,+2). Melanoma oncogenes, such as BRAF and MITF, and tumor suppressor genes, such as CDKN2A/B and PTEN, were included in these alterations. Duplications and deletions were functional as there were significant correlations between DNA copy number and mRNA expression for these genes. Spectral karyotype analysis of three lines confirmed extensive chromosomal instability with polyploidy, aneuploidy, deletions, duplications, and chromosome rearrangements. Bioinformatic analysis identified a signature of gene expression that was correlated with chromosomal instability but this signature provided no clues to the mechanisms of instability. The signature failed to generate a significant (P = 0.105) prediction of melanoma progression in a separate dataset. Chromosomal instability was not correlated with elements of DNA damage response (DDR) such as radiosensitivity, nucleotide excision repair, expression of the DDR biomarkers γH2AX and P‐CHEK2, nor G1 or G2 checkpoint function. Chromosomal instability in melanoma cell lines appears to influence gene function but it is not simply explained by alterations in the system of DDR. Environ. Mol. Mutagen. 55:457–471, 2014. © 2014 Wiley Periodicals, Inc.