The high level phylogenetic relationships within Harmonia. axyridis Pallas (Coleoptera: Coccinellidae) was investigated based on nucleotide data from two molecular markers (12S rRNA, 16S rRNA), along with those of four species obtained from GenBank, including Adalia bipunctata, Calvia quatuordecimguttata, Coccinella septempunctata, Propylea quattuordecimpunctata as outgroup taxa to construct NJ (Neighbor Joining), MP (Maximum Parsimony) and Bayesian trees using Mega 4.0, PAUP (Version 4.0 b10) and Mrbayes (version 3.1.2) software packages. Nucleotides composition, pairwise distances, ratio of transition / transversion, considerable variable sites of those segments were analyzed. The results from phylogenetic trees indicated that H .axyridis var. spectabilis Fald. and H. axyridis ab. lunata Hem. had very close relationships. H. axyrisdis ab.½123456789-undevigintisignata Fald. and H. axyridis ab. succinea Hope. were closely related. H. axyridis var. axyridis Pallas. was an original group, and it had derived to be a independent forma earlier than the other 11 formas within H. axyridis according to the molecular phylogeny analyzing results. Previous hypotheses on the evolution of different color forms of H. axyridis Pallas based on morphological data was supported by our molecular systematic analyzing results.