The two whole chitons were surface dried, then morphological features were examined using a Zeiss Discovery V12 stereomicroscope and photographed using an AxioCam MRc digital camera (Carl Zeiss MicroImaging Inc., Germany). The preserved lengths were ~9.5 mm and 6 mm respectively. Morphologic description was according to published matrix of
Sigwart (2009) and key in
Sigwart and Sirenko (2012). Total genomic DNA was isolated from foot tissue using a CTAB (hexadecyltrimethylammonium bromide) protocol (
Marko, 2002). Polymerase chain reaction (PCR) amplification of the mitochondrial COI gene was carried out using primers LCO1490 (5′-GGTCAACAAATCATAAAGATATTGG-3′) and HCO2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA-3′) (
Folmer et al., 1994). The 28S rDNA was amplified with primers 28SF1 (5′-ACCCGCTGAATTTAAGCATAT-3′) modified from the primer F63mod in
Medina et al. (2001) and 28SR1 (5′-CCATTTAAAGTTTGAGAATAGGT-3′) which is the reverse complement of primer 28S rd4.8a in
Giribet et al. (2006). The 25 µL volume reaction mixture contained 100 ng genomic DNA, 0.15 μmol/L of each primer, 12.5 µL 2×EasyTaq PCR SuperMix (-dye) (TransGen Biotech, Beijing). The reaction cycling profile involved an initial denaturation at 95°C for 5 min, followed by 30 cycles of denaturation at 94°C for 45 s, annealing at 50°C (for COI)/52°C (for 28S rDNA) for 45 s and extension at 72°C for 2 min, and a final extension at 72°C for 10 min. PCR products were purified using the E.Z.N.A.
® Gel Extraction Kit (Omega Bio-tek, Inc., GA, USA), and then were bidirectional sequenced using the PCR primers and the BigDye
® Terminator Cycle Sequencing Kit (v3.1) on an ABI PRISM
® 3 730 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). The forward and reverse sequences were assembled using SeqMan (Lasergene Version 7; DNASTAR, Inc., Madison, WI, USA), and the assembled files were checked by eye. Then homologous sequences were tested using NCBI BLAST (National Center for Biotechnology Information basic local alignment search tool) to confirm that they corresponded with known polyplacophoran sequences deposited in GenBank. To reconstruct the two chitons phylogenetic positions in the order Lepidopleurida, COI and 28S rDNA sequences used for molecular phylogenetic study (
Sigwart et al., 2011) were included. Multiple alignments were performed using ClustalW (
Thompson et al., 1994), then the shortest common length was used for phylogenetic reconstruction. Neighbour joining (NJ) trees of the two genes were reconstructed with the Kimura 2-parameter (K2P) model of nucleotide substitution. The nodal reliabilities were assessed using 1 000 non-parametric bootstrap replicates. Sequence alignments, data sets preparations and phylogenetic analyses were carried out using MEGA 6 (
Tamura et al., 2013).