6 5 2 01 0 1 0 South China Fisheries Science Vol. 6, No. 5 Oct., 2010 doi: 10. 3969/ j. issn. 1673-2227. 2010. 05. 005 COI 16S rrna 1 2 3 1,,, ( 1., 266003; 2., 325400; 3., 116015) : 5 DNA COI 16S rrna, ( Perna viridis) ( Septifer virgatus), ( NJ) ( MP) ( ML) ( BI), 5 [ Mytilus edulis, ( M. galloprovincialis), ( M. trossulus), ( M. coruscus), ( M. californianus) ] 2,, M. edulis, : ; I; 16S rrna; ; ; : Q 244; Q 522 +. 6 : A : 1673-2227 - ( 2010) 05-0027 - 10 Molecular phylogenetics of genus Mytilus based on COI and 16S rrna sequences MAO Yangli 1, CAI Houcai 2, LI Chengjiu 3, GAO Tianxiang 1 ( 1. College of Fisheries, Ocean University of China, Qingdao 266003, China; 2. Nanji Islands National Marine Nature Reserve Administration, Pingyang 325400, China; 3. Fisheries Seed Bureau of Liaoning, Dalian 116015, China) Abstract: In the present study, we compare the partial sequences of mitochondrial COI and 16S rrna of five genus Mytilus to study their phylogenetic relationship. Taking Perna viridis and Septifer virgatus as outgroups, phylogenetic trees are constructed based on neighbor-joining, maximum parsimony, maximum likelihood, Bayesian inference and combined Bayesian analyses methods. The result demonstrates that Mytilus edulis, M. galloprovincialis, M. trossulus, M. coruscus and M. californianus form into two clades. Among the five species of genus Mytilus, M. californianus is the most primitive, then comes M. coruscus. M. galloprovincialis and M. edulis show very close genetic relationship, and M. coruscus is closer to M. californianus than to M. edulis, M. galloprovincialis and M. trossulus. The molecular phylogeny analysis provides theoretical reference for future comparative study on species evolution, immigration and breeding strategies of Mytilus. Key words: genus Mytilus; cytochrome oxidase I; 16S rrna; phylogeny; combined sequences analyses; divergent time : 2010-05-12; : 2010-05 -27 : ( 200905019-2) : ( 1986 - ),,, E-mail: maoyangli123 @ 126. com :, E-mail: gaozhang@ ouc. edu. cn
28 6 ( Mollusca) COI 16S rrna ( Bivalvia) ( Anisomyaria), ( Mytilidae),, 1, [ 1 ] 1. 1, 2 2008 5, Mytilus edulis ( M. gallprovincialis) ( M. trossulus) ( M. californianus), 3 ( M. coruscus) M. edulis 2009 7, 95%,, [ ( Perna viridis), Bank ( 1) [ 2] TO [ 3] [ 12 ], - SHEN [ 4] DNA DNA, COI 2 2, 4 1. 2. 2 PCR, COI COI-F ( 5 -GGTCAACAAAT,, GGTGACCAAAAAATCA-3 ) [ 13 ],, COI 16S rrna, SCARLAT- 1. 2. 1 DNA 100 mg ( Septifer virgatus) ] Gen- CATAAAGATATTGG) COI-R ( 5 -TAAACTTCAG- 16S rrna-ra: 5 -CGCCTGTTTAT- CAAAAACAT-3 ; 16S rrna-rb: 5 -CCGGTCT- DNA ( mtdna) GAACTCAGATCACGT-3 [ 14 ] PCR DNA, 50 L, 10 PCR [ 200 mmol L - 1 Tris-HCl, ph 8. 4; 200 mmol L - 1, KCl; 100 mmol L - 1 ( NH 4 ) 2 SO 4 ; 15 mmol L - 1 [ 5] MgCl 2 ] 5 L, dntp 200 mol L - 1, I ( cytochrome oxidase subunit I, COI) 0. 2 mol L - 1, Taq 1. 25 U (, COI ), DNA 20 ng, Milli-Q H 2 O 50 L [ 6], Eppendorf PCR, COI [ 7] [ 8 ] 94 40 s, 52 30 s, 72 [ 9] 16S rrna 60 s, 35 94 4,, min, 72 7 min 16S rrna [ 10 ], 16S 95 45 s, 50 45 s, 72 45 rrna, s, 35 95 3 min, [ 11 ] 72 10 min 10, 3 2009 2 ; 2 2009 4 1. 2 DNA 100 L TE
5 : COI 16S rrna 29 1 GenBank Tab. 1 Sample species, sampling sites and dates, and their molecular markers and GenBank accession No. species sampling site sampling date 16S rrna abbreviation No. source COI abbreviation No. source M. galloprovincialis 2008. 10 Qingdao M. galloprovincialis1 M. galloprovincialis2 M. galloprovincialis1 M. galloprovincialis2 2009. 02 Dalian M. galloprovincialis3 M. galloprovincialis5 M. galloprovincialis3 M. galloprovincialis5 - - M. galloprovincialis6 AY497292 M. galloprovincialis6 AY497292 M. coruscus Wenzhou 2009. 04 M. coruscus1 M. coruscus2 M. coruscus1 M. coruscus2 2009. 07 Fuding M. coruscus3 M. coruscus5 M. coruscus3 M. coruscus5 - - M. coruscus6 M. coruscus10 GQ472146 GQ472149, AF317545 M. coruscus6 M. coruscus10 FJ495279, FJ495282, FJ495284 FJ495286 M. edulis - - M. edulis1 M. edulis2 AF317054, AY484747 M. edulis1 M. edulis2 AY484747, NC_ 006161 M. trossulus - - M. trossulus1 M. trossulus3 MTU22879 MTU22881 M. trossulus1 M. trossulus2 AF242033, AF242034 M. californianus - - M. californianus AF317057, U68770 M. californianus MCU68776 P. viridis S. virgatus - - P. viridis AB265680 P. viridis DQ343574 - - S. virgatus GQ472165 S. virgatus AB076941 DNA, MEGA 2. 0 [ 15 ] 2 L PCR 1% ( Ts) ( Tv), ( U = 5 V cm - 1 ) UNIQ- 10 DNA ( ) PCR 1. 2. 3 ABI 3700 Ts/Tv ( parsim-informative site) PAUP * 4. 0b10 [ 16 ] COI ( chi- square test, 2 ) PAUP * 4. 0b10 Modeltest DNA, PCR 3. 7 [ 1 7] COI 16S rrna PAUP * 4. 0b10 MrBayes 3. 0 [ 18 ] 1. 3 DNAStar ( DNASTAR, Inc ), ( NJ) ( MP) SeqMan, COI ( ML) ( Bayesian Inference, BI) 3 16S rrna EditSeq, PAUP * 4. 0b10, K2P 1 COI ;, 4 ; MegAlign,, 3 1 ; 10
30 6 10 4 20 10 4, 100 1 16S rrna MrModeltest 2. 0 [ 2 0], ( 2 500 ), PAUP * 4. 0b10 2 4, ( 2), P, ( partition homogeneity test) [ 19 ] COI 1 2 3,, 2 Tab. 2 Partitioning strategies used in this study partitioning strategies partitioning method P 1 P 2 P 3 1 COI; 16S rrna COI 1 2 ; COI 3 ; 16S rrna P 4 COI 3 ; 16S rrna 2 2. 1 DNA 5 COI 432 bp COI 16S rrna 0. 014 432 bp 121, 107 0. 006, M. edulis ; A T C G 26. 4% 35. 0% 16. 8% 21. 7%, A + T 61. 4%, AT / 2. 79, TC,, TA CG 3, COI 2 : 1 = 75. 218, df = 63, P 2 = 0. 139; 2 = 50. 226, df = 63, P = 2 0. 878; 3 = 42. 428, df =63, P =0. 978 COI 3 1 5 16S rrna 460 bp 460 bp 68, 56 ; A T C G 29. 6% 30. 6% 15. 8% 24. 0%, A + T 60. 2%, AT / 2. 53,, K TC AG,, 2 COI TA, CG 2. 2 3 COI 16S rrna 0. 014 0. 219 0. 006 0. 078 M. edulis,, 0. 219 ( COI ) 0. 078 ( 16S rrna) 0. 127 ( COI) 0. 023 ( 16S rrna),, ( synon- ymous substitution), 3,, COI 143, M. edulis 1 ( amino acid substitution) R = K/2T, K, T 2, > 16S rrna, COI ( 0. 203) > 16S rrna ( 0. 070) COI COI 16S rrna 16S rrna ( 3)
5 : COI 16S rrna 31 3 K2-P ( : 16S rrna, : COI) Tab. 3 Average K2-P genetic distance based on 16S rrna region ( below diagonal) and COI region ( above diagonal) among different species of genus Mytilus species M. galloprovincialis M. coruscus M. edulis M. trossulus M. californianus M. galloprovincialis - 0. 203 0. 014 0. 139 0. 178 M. coruscus 0. 070-0. 219 0. 187 0. 127 M. edulis 0. 006 0. 078-0. 154 0. 188 M. trossulus 0. 028 0. 072 0. 033-0. 146 M. californianus 0. 068 0. 023 0. 071 0. 071-1 COI MP ML BI (, ) Fig. 1 Phylogenetic tree of Mytilus species constructed with neighbor-joining, maximum parsimony, maximum likelihood and Bayes analyses based on COI ( numbers on the tree represent posterior probability values and supporting values; the same case in the following figure) ( 350 ) ( Late Miocene) ( Early Pliocene) M. edulis ( Late Pleistocene) [ 21], COI 16S rrna 2. 3, M. edulis 498 827 2. 3. 1 COI 16S rrna, 427 COI 16S rrna, NJ MP ML 753, 31. 8 BI ( 1 2), 63. 6,
32 6 2 16S rrna MP ML BI Fig. 2 Phylogenetic tree of Mytilus species constructed with neighbor-joining, maximum parsimony, maximum likelihood and Bayes analyses based on 16S rrna 5 2 ( A B), A M. edulis, ; B 1, M. edulis 1, M. edulis 1 M. edulis 3 2. 3. 2 COI 16S rrna,, 12S rrna [ 2 2] COI 16S rrna 2 16S rrna,,, COI, PAUP COI 16S rrna 2 ( partition 460 bp 16S rrna, 5 homogeneity test) [ 19 ], 2 [ 2 3],, 68, 432 bp COI P 0. 07, 2, 121 901,, 5 COI 4 16S rrna COI 16S rrna, P 3, COI ( 3) 16S rrna
5 : COI 16S rrna 33 4 Tab. 4 Partitions, optimal model of sequence evolution and total number of characters in each partition partition evolutionary model number of characters in partition COI GTR + G 432 COI 1 HKY + G 144 COI 2 HKY 144 COI 3 HKY + G 144 COI 1 2 HKY + G 288 16S rrna GTR + G 469 GTR + I + G 901 Fig. 3 3 COI 16S rrna ( P 3 ) 50% ( 50% ) 50% consensus tree of Mytilus species constructed with Bayesian analysis ( strategy P 3 ) based on COI and 16S rrna sequences. Bayesian posterior probabilities are indicated above the branch ( bootstrap scores larger than 50% were shown) 2 AT DNA A + T, 61. 4% 60. 2%, G + C, 16S rrna [ 24-26] 12S rrna COI, [ 2 7] [ 2 8],
34 6 [ 29-32 ], M. edulis mtdna 2 ( F, M ) [ 40], M. edulis KOEHN [ 3 3] 5, VARVIO [ 41] M. edulis Mc- ( geological-palaeontological) DONALD [ 28], M. edulis,, 3 2, 1,, McDONALD [ 28] 3 3 ( lineage), NJ MP ML BI 1, 350 [ 2 1] 3 M. edulis, M. edulis,, 498 827, 1, ( > 80) ( 1 3) M. edulis 0. 014 ORT POGSON [ 42 ] M. edulis ( COI) 0. 006 ( 16S rrna) ( 7. 64 7. 61 MYA), 0. 021 ( COI ) 0. 007 ( 16S rrna), 1 ( 1 3 ) M. edulis ( Early Pliocene) ( Late Pleistocene) 15. 4% ( COI), RIGINOS [ 34 ], COI 16S rrna 17. 1% ( 4. 89% / MY) [ 35 ], rrna 16,, M. edulis, HILBISH [ 36 ], 3 :,, 1, M. edulis GERARD [ 37 ] 16S rrna COI :,, 16S rrna WENNE SKIBINSKI [ 39] DNA 3, 2 M. edulis [ 38 ], [ 38 ] 16S rrna 282 bp, 469 bp, COI 2, WENNE SKIBINSKI [ 3 9] ( relict), 427 753, 31. 8 63. 6 ( Late Miocene ),! [ 1] EYNON C M, SKIBINSKI D O F. The evolutionary relationships between three species of mussel ( Mytilus) based on anonymous DNA polymorphisms [ J]. J Exp Mar Biol Ecol, 1996, 203 ( 1) : 1-10. [ 2] JOHANNESSON K, KAUTSKY N, TEDENGREN M. Genotypic and phenotypic differences between Baltic and North Sea populations of Mytilus edulis evaluated through reciprocal transplantations. Genetic variation [ J]. Mar Ecol Prog Ser, 1990, 59 ( 3 ) : 211-219. [ 3] SCARLATTO O A. Bivalve mollusks of temperate latitudes of the western portion of the Pacific Ocean [ M]. Moscow: Opred Faune SSSR, 1981, 126: 1-461. [ 4] SHEN Yubang, LI Jiale, FENG Bingbing. Genetic analysis of cul- II
5 : COI 16S rrna 35 tured and wild populations of Mytilus coruscus based on mitochondrial DNA [ J]. Zool Res, 2009, 30 ( 3 ) : 240-246. [ 5 ] AVISE J C. Phylogeography [ M]. London: Harvard University Press, 2000: 1-36. [ 6] HOEH W R, STEWART D T, SUTHERLAND B W, et al. Cytochrome c oxidase sequence comparisons suggest an unusually high rate of mitochondria DNA evolution in Mytilus ( Mollusca: Bivalvia) [ J]. Mol Biol Evol, 1996, 13 ( 2) : 418-421. [ 7],,,. 3 16 S rrna COI [ J]. 2009, 5 ( 3) : 46-521., 2005, 29 ( 8) : 27-32. [ 8] KING T L, EACKLES M S, GJETVAJ B, et al. Intraspecific phylogeography of Lasmigona subviridis ( Bivalvia: Unionidae) : conservation implications of range discontinuity [ J]. Mol Ecol, 1999, 8 ( S1 ) : 65-78. [ 9] AN H S, JEE Y J, MIN K S, et al. Phylogenetic analysis of six species of Pacific abalone ( Haliotidae) based on DNA sequences of 16 S rrna and cytochrome c oxidase subunit I mitochondrial genes [ J]. Mar Biotech, 2005, 7 ( 4) : 373-380. [ 10 ],,,. RNA 3 [ J]., 1997, 24 ( 4 ) : 296-304. [ 11],,,. 8. :, 2001 : 59-63. 16S rrna [ J]. [ 27],,,. 16S rrna, 2004, 14 ( 5) : 514-521. [ 12] SAMBROOK J, FRITSCH E F, MANIATIS T, et al. Molecular cloning: a laboratory manual [ M]. New York: Cold Spring Harbor Laboratory Press, 1996 : 464-468. [ 13] FOLMER O, BLACK M, HOEH W, et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates [ J]. Mol Mar Biol Biotech, 1994, 3 ( 5) : 294-299. [ 14] PALUMBI S R, MARTIN A P, ROMANO S, et al. The simple fool s guide to PCR. version 2. 0 [ M]. Honolulu: Zoology, University of Hawaii, 1991. Department of [ 15] KUMAR S, TAMURA K, JAKOBSEN I B, et al. MEGA: molecular evolutionary genetics analysis software [ J]. 2001, 17 ( 12) : 1244-1245. Bioinformatics, [ 16 ] SWOFFORD D L. PAUP *. Phylogenetic analysis using parsimony ( * and other methods), version 4 [ M]. Sunderland, MA: Sinauer Assoc, 2002. [ 17 ] POSADA D, CRANDALL K A. Modeltest: testing the model of DNA substitution [ J]. Bioinformatics, 1998, 14 ( 9 ) : 817-818. [ 18] RONQUIST F, HUELSENBECK J P. Mrbayes 3: Bayesian phylogenetic inference under mixed models [ J]. Bioinformatics, 2003, 19 ( 12) : 1572-1574. [ 19] FARRIS J, KALLERSJO M, KLUGE A, et al. Testing significance of incongruence [ J]. Cladistics, 1995, 10: 315-319. [ 20] NYLANDER J. MrModeltest 2. 0. Program distributed by the author [ M]. Sweden: Evolutionary Biology Centre, Uppsala University, 2004. [ 21] RIGINOS C, HICKERSON M J, HENZLER C M, et al. Differential patterns of male and female mtdna exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis [ J]. Evol, 2004, 58 ( 11) : 2438-2451. [ 22],,,. 4 16S rrna COI Cyt b [ J]., [ 23 ] HOWLAND D E, HEWITT G M. Phylogeny of the coleoptera based on mitochondrial cytochrome oxidase I sequence data [ J]. Insect Mol Biol, 1995, 4 ( 3) : 203-205. [ 24 ] ANDERSON F E. Phylogeny and historical biogeography of the loliginid squids ( Mollusca: Cephalopoda) based on mitochondrial DNA sequence data [ J]. Mol Phylogenet Evol, 2000, 15 ( 2 ) : 191-214. [ 25],,,. 12 S rrna [ J]., 2000, 30 ( 1) : 43-47. [ 26],,,. DNA 16S rrna [ G]. [ J]., 2009, 5 ( 6) : 47-53. [ 28] McDONALD J H, SEED R, AND KOEHN R K. Allozymes and morphometric characters of three species of Mytilus in the Northern and Southern Hemispheres [ J]. Mar Biol, 1991, 111 ( 3 ) : 323-333. [ 29] KOEHM R K. The genetics and taxonomy of species in the genus Mytilus [ J]. Aquac, 1991, 94 ( 2 /3) : 125-145. [ 30] MCDONALD J H, KOEHM R K. The mussels Mytilus galloprovincialis and M. trossulus on the Pacific of North America [ J]. Mar Biol, 1988, 99 ( 1) : 111-118. [ 31] EDWARDS C A, SKIBINSKI D O F. Genetic variation of mitochondrial DNA in mussel ( Mytilus edulis and M. galloprovincialis) populations from South West England and South Wales [ J]. Mar Biol, 1987, 94 ( 4) : 547-556. [ 32] VARVIO S L, KOEHN R, VAINOLA R. Evolutionary genetics of the Mytilus edulis complex in the North Atlantic region [ J]. Mar Biol, 1988, 98 ( 1) : 51-60. [ 33] KOEHN R K, HALL J G, INNES D J, et al. Genetic differentiation in Mytilus edulis in eastern North America [ J]. 1984, 79 ( 2) : 117-126. Mar Biol, [ 34] RIGINOS C, HICKERSON M J, HENZLER C M, et al. Differential patterns of male and female mtdna exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis [ J]. Evol, 2004, 58 ( 11) : 2438-2451. [ 35],,,.
36 6 [ J]., 2009, 31 ( 11 ) : 1-11. [ 36] HILBISH T J, MULLINAX A, DOLVEN S I, et al. Origin of the antitropical distribution pattern in marine mussels ( Mytilus spp. ) : routes and timing of transequatorial migration [ J]. 2000, 136 ( 1) : 69-77. Mar Biol, [ 37] GERARD K, BIERNE N, BORSA P, et al. Pleistocene separation of mitochondrial lineages of Mytilus spp. mussels from northern and southern hemispheres and strong genetic differentiation among southern populations [ J]. Mol Phylogenet Evol, 2008, 49 ( 1) : 84-91. [ 38],,. 16S rrna mussel, Mytilus californianus [ J]. Genet, 2007, 177 ( 2 ) : [ J]., 2009, 33 ( 12 ) : 50-55. in European populations of the mussel, Mytilus trossulus [ J]. Mar Biol, 1995, 122 ( 4) : 619-624. [ 40] ZOUROS E, BALL A O, SAAVEDRA C, et al. A unusual type of mitochondrial DNA inheritance in the blue mussel Mytilus [ J]. Proc Natl Acad Sci USA, 1994, 91 ( 16) : 7463-7467. [ 41] VARVIO S L, KOEHN R K, VAINOLA R. Evolutionary genetics of the Mytilus edulis complex in the North Atlantic region [ J]. Mar Biol, 1988, 98 ( 1) : 51-60. [ 42] ORT B S, POGSON G. H. Molecular population genetics of the male and female mitochondrial DNA molecules of the California Sea 1087-1099. [ 39] WENNE R, SKIBINSKI D O F. Mitochondrial DNA heteroplasmy