The Third Technical Consultation Meeting (3rd TCM) Research for Stock Enhancement of Sea Turtles (Japanese Trust Fund IV Program) 7 October 2008 REGIONAL ANALYSIS ON STOCK IDENTIFICATION OF GREEN AND HAWKSBILL TURTLES IN THE SOUTHEAST ASIAN REGION Wahidah Mohd Arshaad, Syed Abdullah Syed Abdul Kadir 2 Marine Fishery Resources Development and Management Department, Taman Perikanan Chendering,2080 Kuala Terengganu, Terengganu. 2 Turtle and Marine Ecosysten Centre, 200 Dungun, Terengganu. (wahidah@seafdec.org.my) INTRODUCTION Knowledge of population dynamics is largely obtained from long-term mark-recapture studies of females tagged while nesting on the beaches. These studies show that breeding female turtles display high fidelity to the same nesting beaches (Hendrickson 98; Carr & Ogren 960) and Carr (967) hypothesized that mature nesting female turtles were selecting their natal beach to deposit eggs. Studies in the southern Great Barrier Reef (sgbr) demonstrate that green turtles also display fidelity to resident feeding grounds throughout their adult lives (Limpus et al. 992) with females reaching sexual maturity at about 40 years (Limpus & Chaloupka 997). Detection of subpopulations of sea turtles is a very essential ecological aspect to promote their conservation and enhancement because the resources-protect schemes should be made to the each population unit. Mitochondrial DNA has proved particularly effective for detecting population structure in marine turtles. Analysis of mtdna structure in Atlantic green turtle population supported the natal homing hypothesis, as geographically distant were found to have heterogeneous mtdna frequencies (Bowen & Avise 996). ) To detect subpopulations of nesting green and hawksbill turtles in Southeast Asian region OBJECTIVE 2) To identify the genetic markers of different management units
Table : Number of samples green turtles collected for the Trust Fund IV. METHOD Country Location N Brunei Brunei Beach 4 Indonesia Derawan Sangalaki Panjang 22 9 6 Malaysia Myanmar Philippines Thailand Vietnam Redang Sarawak Turtle s Sabah Turtle s Tameahla Coco Panikian APO Reef NP Baguan Huyong Khram Con Dao Minh Chau 3 9 4 23 9 34 Total 332 Figure : Sampling sites of green turtle for Trust Fund IV study. Table 2: Number of samples hawksbill turtles collected for the Trust Fund IV. Thamee Hla Coco Huyong Minh Chau Khram Con Dao Redang Sarawak Turtle s APO Reef National Park Brunei Beach Philiphines Turtle s Panikian Sabah Turtle Derawan/Sangalaki/Panjang Country Location N Brunei Brunei Beach 4 Indonesia Kimar Belitung 38 Malaysia Melaka Sabah Turtle s 29 Myanmar Coco 4 Philipphine Dava0 Gulf APO Reef Natural Park Misamis Oriental Bataan Thailand Khram 9 Total 36 Figure 2: Sampling sites of hawksbill turtle for Trust Fund IV study. Coco Khram Melaka APO Reef National Park Brunei Beach Misamis Oriental Davao Gulf Sabah Turtle SAMPLES ANALYSIS Belitung 2
DNA was extracted from skin tissue PCR mtdna DNA control region - primer TCR and TCR6 (green), LTEi9 and H90 (hawksbill) - target an area of 389 bp and 740 bp of the control region Sequenced PCR product forward and reverse by automated sequencer Statistical analysis - Align sequences (Clustal X/W) - Blast with GenBank database for similar sequences - MEGA 4 (phlogenetic tree) - Arlequin 3. (estimate h and π, population differentiation & quantify genetic diversity) RESULT Figure 3: Neighbour-joining tree of Kimura 2-parameter distance for green turtle haplotypes based on mtdna control region sequences, using the hawksbill haplotype as an outgroup. GREEN SEA TURTLE Table 3: Frequencies of green turtles (Chelonia mydas) mtdna haplotypes. Clade I Clade II Clade III C CMP G G G G G Haplotypes M2 C Country Location N R T M C4 D2 C3 M3 C7 P2 B T2 T3 A2 A3 A6 P 8 4 G2 G3 G4 G 2 3 4 G7 No. 3 Brunei Brunei Beach 4 2 Indonesia Derawan 22 3 8 6 Sangalaki 38 8 8 3 3 Panjang 6 7 4 2 6 Pengumbahan 23 7 2 6 Enu 28 2 27 Redang 73 3 9 Malaysia 9 42 4 3 Paka 4 2 Pahang 2 8 4 2 Perak 3 2 2 Sarawak 62 6 40 4 Sipadan 40 24 8 4 3 Sabah Turtle s 8 47 0 3 Tameahla 8 4 Myanmar 2 9 Coco 4 7 3 2 2 7 Philippine Panikian 9 9 APO Reef NP 4 4 Phi Turtle s 62 44 7 3 Huyong 9 0 4 Thailand 8 Khram 8 6 3 Con Dao 42 40 3 Vietnam Minh Chau Total 643 293 3 20 9 6 48 24 2 2 4 26 3 2 2 4 3 * Haplotypes in black () was published by Moritz et al. (2002) and haplotypes in red (9) are new. Figure 4: Neighbour-joining tree of Kimura 2-parameter distance for all sample rookeries of green turtle haplotypes based on mtdna control region sequences, using the hawksbill haplotype as an outgroup. 3
Table 4: Polymorphic sites in the 38 bp sequence Polymorphic sites (47) Haplo 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 types 4 4 6 6 6 6 7 7 7 0 4 8 9 0 4 6 7 7 9 9 9 0 0 2 3 6 7 7 () 6 7 6 8 0 2 9 0 7 8 0 4 3 8 4 7 8 8 3 3 7 9 0 0 6 7 8 3 4 3 7 8 9 4 6 7 8 7 8 D2 G A G A C C A - T G T A T A C G G G G G G T G A A G C C C T T G G C A G T A C T - - - - G A A CMP8....... -................. A.............. - - - -... R....... -. A............................. G C A G... M....... -. A.............................. - - - -. G. T....... -. A C.................... A A....... - - - -... C4....... -. A.............................. - - - -... C3....... -. A.............................. - - - -... M3 A...... -. A.............................. - - - -... C7....... -. A.......................... C... - - - -... C.... T.. -. A.............................. - - - -... M2....... -. A.............................. - - - -... C4....... -. A.......... A................... - - - -... P2....... -. A.............................. - - - - A.. G2....... -......... A...................... - - - -... G3....... -................................ - - - - A.. G4. G..... -.... C........................... - - - -... G....... -............................... C - - - -... G....... -........ G.. A.................... - - - -... G2....... -................................ - - - -... G3....... -..... G.......................... - - - -... G....... -.......... A..................... - - - -... B A. A. T T G A. A. G.. T A.... A C. G G.. T T C........ T C - - - -... T2 A. A. T T G A. A. G.. T A.... A C. G G.. T T C... G... T C - - - -... T3 A. A. T T G A. A. G.. T A.... A C. G G.. T T C...... C. T C - - - -... G4 A. A. T T T -... G.. T A.... A A. G G.. T T C........ T C - - - -... A2.... T T G - C..... T A.... A C A G G. T. T C... T. T C G T. - - - -.. G A3... G T T G - C A.... T A.... A C A G G. T T T C... T. A C G T. - - - -.. G P... G T T G - C A.... T A.... A C A G G. T T T C... T. A C G T. - - - -... A6.... T T G - C A.... T A.... A C A G G. T T T C... T. A C G T. - - - -.. G G7... G T T G - C..... T A.... A C A G G. T T..... T. A C G T. - - - -... Table : Haplotype (h) and nucleotide (π) diversities by nesting beaches. No. Country Nesting beach N h hd π Brunei Brunei Beach 4 3 0.833 0.0026 2 Indonesia Derawan 22 7 0.836 0.026 3 Sangalaki 38 6 0.79 0.006 4 Panjang 6 6 0.767 0.00 Pengumbahan 23 2 0.403 0.00 6 Enu 28 2 0.07 0.0039 7 Malaysia Redang 73 9 0.647 0.0 8 Paka 2 0.33 0.0062 9 Pahang 2 4 0.6 0.002 0 Perak 2 0.248 0.0006 Sarawak 62 0.2 0.0072 2 Sipadan 40 4 0.99 0.0084 3 Sabah Turtle s 8 4 0.39 0.0009 4 Myanmar Tamaehla 4 0.63 0.006 Coco 7 0.662 0.0087 6 APO Reef Natural 4 0 0 Philiphine Park 7 Panikian 9 0 0 8 Philiphine Turtle s 62 3 0.427 0.002 9 Thailand Huyong 9 3 0.73 0.002 20 Khram 7 0.687 0.00247 2 Vietnam Minh Chau 0 0 22 Con Dao 42 3 0.094 0.0002 Overall 628 0.442 0.004 Table 6: Haplotype (h) and nucleotide (π) diversities by management units (Exact test of sample differentiation based on haplotype frequencies, Raymond M. and F. Rousset. 99). No. Country Nesting beach N h hd π Indonesia Derawan 22 7 0.836 0.026 2 Sangalaki 38 6 0.79 0.006 3 Pengumbahan 23 2 0.403 0.00 4 Enu 28 2 0.07 0.0039 Malaysia/ Redang+Paka+Perak 03 0 0.39 0.0092 6 Philippine Pahang 2 4 0.6 0.002 7 Sarawak 62 0.2 0.0072 8 Sipadan 40 4 0.99 0.0084 9 Sabah+Philippine TI 20 4 0.377 0.00 0 Myanmar Tamaehla+Coco 60 7 0.68 0.003 Thailand Huyong+Khram 49 8 0.640 0.0243 2 Vietnam Minh Chau+Con Dao 3 3 0.07 0.0020 Overall 60 0.028 0.0073 Figure : Location of 2 genetically distinct breeding stocks or management units (Exact test of sample differentiation based on haplotype frequencies Raymond M. and F. Rousset. 99). Thamee Hla Coco Huyong Khram Con Dao Redang Pahang Sarawak Turtle s Pengumbahan Minh Chau Sangalaki Philiphines Turtle s Sipadan Derawan Sabah Turtle Enu Table 7 : Frequencies of hawksbill turtles mtdna haplotypes. HAWKSBILL TURTLE Clade I Clade II Country Location N N 2 Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ Ei_ H 2 3 4 6 7 8 9 0 2 3 4 Brunei Brunei Beach 4 4 2 3 Indonesia Kimar Belitung 38 9 7 2 2 Malaysia Melaka 29 29 29 Sabah Turtle s 20 9 8 Myanmar Coco 4 4 2 3 Philipphine Dava0 Gulf 2 7 APO Reef Natural Park 4 2 Misamis Oriental Bataan Thailand Khram 9 4 9 3 2 3 Total 36 88 2 9 2 48 3 3 2 2 *Note: N = No. of samples received, N 2 = No. of samples analyzed, H = No. of haplotypes. 4
Figure 6: Neighbour-joining tree of Kimura 2-parameter distance for all green turtle haplotypes based on mtdna control region sequences. Figure 7: Neighbour-joining tree of Kimura 2-parameter distance for all sample countries of hawksbill turtle haplotypes based on 6 mtdna control region sequences. Haplotypes () 6 Table 8: Polymorphic sites in the 740 bp sequence Polymorphic sites (43) 2 2 2 2 2 2 2 2 2 3 3 3 3 3 4 4 4 6 4 6 7 8 2 3 3 3 4 8 8 8 9 2 2 4 4 7 7 9 2 6 6 6 2 8 0 2 7 7 7 3 0 4 4 2 6 0 4 2 4 6 0 4 6 4 8 4 4 9 6 3 4 6 8 8 0 7 8 9 9 Ei_ A - A G A G T T G A C T C A T T A A G A C A T A A A A T G A T A C G A G A A G T G A A G Ei_2. - G......................................... Ei_3. -..................................... C.... Ei_4. -.................................. G G.... G. Ei_. -.......... T....................... G G.... G. Ei_6. -.......... T... G.......... G..... C.......... Ei_7. -.......... T... G....... G................... Ei_8. -................. G..... G.................. Ei_9 G T. A C T C C A. T C T. C C. C A G T. C.. G C. G..... G A.. A C A. G A Ei_0 G T. A C T C C A. T C.. C C. C A G T. C... C. G..... G A.. A C A. G A Ei_ G T. A C T C C A. T C T G C C. C A G T G C... C....... G A.. A C A G G A Ei_2 G T. A C T C C A. T C T. C C. C A G T G C... C....... G A.. A C A. G A Ei_3 G T. A C T C C A. T C T. C C. C A G T. C G.. C G...... G A.. A C A. G A Ei_4 G T. A C T C C A. T C. G C C. C A G T. C... C....... G A.. A C A. G A Ei_ G T.. C T C. A G T C... C. C. G T. C... C.. G C A G..... A. A. G. 6 0 Table 9: Haplotype (h) and nucleotide (π) diversities by nesting beaches. No. Country Nesting Beach N H h π. Brunei Brunei Beach 4 3 0.83333 0.0968 2. Indonesia Kimar Belitung 9 2 0.38889 0.0002 3. Malaysia Melaka 29 0.0 0.0 Sabah Turtle s 20 0.6636 0.020392 4. Myanmar Coco 4 3 0.83333 0.027440. Philiphine Dava0 Gulf 8 7 0.96429 0.08942 APO Reef Natural Park Misamis Oriental Bataan 6. Thailand Khram 4 3 0.6044 0.00233 Total 88 *Note: H = No. of haplotypes, h = haplotype diversity and π = nucleotide diversity. GREEN TURTLE haplotypes were detected which was published by Moritz et al. 2002 and 9 are new. Dominant haplotypes are C3, found at all samples sites except Enu, Indonesia and Panikian, Philiphine. DISCUSSION C3 dominant haplotype for Pengumbahan (7/23), Redang (42/73), Paka (4/), Pahang (8/2), Perak (3/), Sipadan (24/40), Tameahla (8/), Coco (7/), Huyong (0/9), Khram (/), Vietnam (/3).
Dominant: C4, Sarawak (40/62), D2, Sabah (47/8) & Philippine Turtle (44/62), C4, Enu (27/28), A3, Panikian (9/9). haplotypes are share between nesting beach and 9 haplotypes are unique to individual nesting beach. HAWKSBILL TURTLE Due to the limited of samples number, the output of this study will not give a true figure of hawksbill population genetic. 3 haplotypes are share between nesting beach and 2 haplotypes are unique to individual rookeries. Melaka, no haplotype diversity. Table 0: Polymorphic sites of haplotypes found in Melaka (392 bp) Figure 6: Phylogenetic tree of haplotypes found in Melaka. Polymorphic sites (43) Haplo- 2 2 3 3 3 3 3 3 3 3 3 3 types 2 7 9 8 8 9 0 3 6 7 7 8 8 8 9 9 () 2 8 2 2 6 9 7 3 8 2 3 4 6 7 0 H2 C G G G G T C T A C A G T C G C A H3. A. A.. T C G........ H2... A A. T.......... H3. A. A A G T.......... Ei_9 T A A A A. T.. A G A C T A G T CONCLUSION This had study identified 2 genetically breeding aggregations of green turtle throughout Southeast Asia (Management Units). The population diversity overall for green turtle are still high except for Enu and Vietnam. More effort should be focused on these nesting populations to keep population size high and maintain this diversity. RECOMMENDATION 6
Research on feeding ground must done to know the feeding ground of haplotype C3 for green turtles. More work should do on identify the genetic of hawksbill turtles. AKNOWLEDGEMENT Japanese Trust Fund (Trust Fund IV, Research for stock enhancement of sea turtles) All Technical Members Countries (Brunei, Cambodia, Indonesia, Malaysia, Thailand,The Philiphines and Vietnam. Dr. Juanita Joseph (University Malaysia Terengganu) Nurul Azrin Yusof and Nik Zuraini Omar. Contract staffs, Ms. Zuhailah, Mr. Abdul Hadi and Ms. Tengku Noor Aminah. All practical students. 7