genetic diversity of kelampayan using dominant dna

2 downloads 0 Views 1MB Size Report
Three ISSR primers, namely (GTG)6, (AG)10 and (AC)10, that yielded reproducible ... tree that gives economic returns within 8 to 10 y (Joker 2000). The tree can ...
GENETIC DIVERSITY OF KELAMPAYAN USING DOMINANT DNA MARKERS BASED ON INTER-SIMPLE SEQUENCE REPEATS IN SARAWAK W.S. Ho1, K.S. Liew1, A.B. Elias1, F.B. Fedrick1, M.A. Mohamad Khairil1, S.L. Phui2 & A. Julaihi2 1

2

Forest Genomics and Informatics Laboratory (fGiL) Department of Molecular Biology Faculty of Resource Science and Technology Universiti Malaysia Sarawak 94300 Kota Samarahan, Sarawak

Applied Forest Science and Industry Development (AFSID) Sarawak Forestry Corporation 93250 Kuching, Sarawak Email: [email protected]

Abstract Neolamarckia cadamba (Roxb.) Bosser, or locally known as kelampayan, is a fast-growing timber species which produces one of the best sources of raw material for the plywood industry and also for the pulp and paper industry. It has been selected as one of the promising plantation tree species for large-scale planted forest development in Sarawak. Therefore, the molecular characterization of this indigenous tropical tree species is needed to maintain its high quality. Intersimple sequence repeats (ISSR) markers were used in this study to determine the genetic diversity of kelampayan in three progeny trial blocks at the Landeh Nature Reserve, Semengok, Sarawak. The seeds were collected from the selected mother trees located at the Pasai Bon, Niah and Lawas seed production areas (SPAs) in Sarawak. Three ISSR primers, namely (GTG)6, (AG)10 and (AC)10, that yielded reproducible, informative and scorable fragments were chosen for ISSR analysis. A total of 64 loci were generated of which 45.3–74.6% of the loci were identified as polymorphic bands with the size ranging from 500 bp to 2 kb among 247 kelampayan progenies selected in the present study. Molecular diversity based on Shannon’s diversity indices (I) among 247 trees ranged from 0.268 to 0.350. In general, the kelampayan trees in the three progeny trial blocks exhibited a high level of molecular diversity and DNA polymorphism compared with its natural populations. This preliminary information will form the base for kelampayan tree improvement and conservation programmes.

332

Introduction The Sarawak state government has set a target to establish one million hectares of planted forest by the year 2020. The state government has issued 43 licenses for planted forests (LPFs) and about 250 000 ha of planted forests have been established to date. Planted forests are a viable and sustainable method to produce high-value commercial timber to meet the current increase in global demand for wood. Datuk Len Talif Salleh (State Forestry Department Director) stated that planted forests can reduce dependency on natural forest and generate about 10 to 20 times more wood volume than the natural forests (Borneo Post Online 2010). Besides, planted forests also can provide a number of social and environmental services such as rehabilitation of degraded lands, soil and water protection, sequestering and storing carbon, conservation of biological diversity, creating rural employment, helping communities raise their standard of living and contributing to sustainable development (FAO 2010). In order to achieve this goal, approximately 30 million high quality planting materials are required annually for planted forest development in Sarawak. In this regard, forest genetics and tree improvement research will help respond to the need to develop adequate tools for producing good quality seedlings that are of faster growth, high yield and high wood quality in the shortest of time at a reasonable cost. With the advent in DNA sequencing, data analysis and PCR technology, a diverse array of DNA marker systems have been described in the literature such as SSR, RAPD, ISSR, RFLP and AFLP (Semagn et al. 2006). These DNA markers will greatly facilitate the selection of quality planting materials for planted forest development in Sarawak. Besides, such DNA markers also have proven their utility in the analysis of phylogenetic relationships, population structure, mating system, gene flow, parental assignment, introgressive hybridization, marker-assisted selection and genetic linkage (Kumar et al. 2009). ISSR technique is a PCR-based method to amplify DNA fragments between two closely spaced and oppositely oriented SSRs (Moreno et al. 1998). Compared to other DNA fingerprinting methods, ISSR markers offer several advantages such as: a) fast, inexpensive and no radioactive handling facilities required, b) does not require the knowledge of flanking sequence, c) highly dispersed throughout the genome, and d) highly polymorphic (Galvan et al. 2003). A study carried out by Moreno et al. (1998) demonstrated that the reproducibility of ISSR (91.8%) was higher than RAPD (85.8%) due to the use of longer ISSR primers (16–25 mers) which permits much more stringent annealing conditions. ISSR also generates a large number of markers by amplifying multiple microsatellite loci, thus allowing screening of a large number of samples in a single gel (Nagaraju et al. 2001). Therefore, these markers are useful in studies on genetic diversity (Balasaravanan et al. 2005, Okun et al. 2007, Chezhian et al. 2010), gene tagging (Ammiraju et al. 2001), cultivar identification (Wong et al. 2009), genome mapping (Zietkiewicz et al. 1994, Godwin et al. 1997) and phylogenetic analysis (Dogan et al. 2007). Neolarmarckia cadamba (Roxb.) Bosser, or locally known as kelampayan, is one of the selected indigenous species for planted forest development in Sarawak (Figure 1). Kelampayan is characterized as a large, deciduous and fast-growing tree that gives economic returns within 8 to 10 y (Joker 2000). The tree can grow 333

up to 40–45 m tall and 100–160 cm in diameter. Kelampayan is one of the best sources of raw material for the plywood industry, besides pulp and paper production (Joker 2000). It also can be used for other purposes such as picture frame, disposable chopstick, wooden sandal, general utility furniture and plywood. Moreover, the leaves and bark of kelampayan are reported to possess various medicinal uses (Patel & Kumar 2008, Mondal et al. 2009).

Figure 1

Morphological characteristics of N. cadamba. (a) Mature tree; (b) Twig with inflorescence; (c) Fruit

Despite kelampayan being an economically and pharmaceutically important timber species, genetic information of this species is still scanty until now. Thus, we applied ISSR markers to determine the genetic diversity of kelampayan in three progeny trial blocks at the Landeh Nature Reserve, Semengok, Sarawak.

Materials and Methods A total of 247 kelampayan samples (83, 74 and 90 samples from Blocks 1, 2 and 3; respectively) were collected from the Kelampayan Provenance Trial Plot at Landeh Natural Reserve, Sarawak (Figure 1, Table 1). Total genomic DNA was extracted from the fresh leaf samples using a modified CTAB method (Doyle & Doyle 1990). The quality and quantity of the extracted DNA were estimated spectrophotometrically and verified using a 0.8% agarose gel. The DNA was then subjected to ISSR-PCR amplifications.

334

Table 1 Origins of kelampayan tree according to lines in three different blocks Line

Tree origin Block 1

Block 2

Block 3

1

Pasai Bon

Niah

Lawas

2

Niah

Pasai Bon

Pasai Bon

3

Pasai Bon

Pasai Bon

Lawas

4

Niah

Niah

Niah

5

Lawas

Lawas

Pasai Bon

6

Niah

Pasai Bon

Lawas

7

Lawas

Lawas

Niah

8

Lawas

Niah

Niah

9

Pasai Bon

Lawas

Pasai Bon

Block 1 Block 2 Block 3

Figure 1 Kelampayan progeny trial plot at Landeh Nature Reserve, Sarawak

Three microsatellite primers, namely (AC)10, (AG)10 and (GTG)6, were used in this study to amplify the ISSR region. PCR was carried out using a Mastercycler Gradient PCR (Eppendorf, Germany). Amplification was carried out in 25-µl reaction volume containing 1X PCR buffer (10 mM Tris-HCl at pH 8.8 and 50 mM KCl), 2.5mM MgCl2, 0.2 mM of each dNTPs (dATP, dCTP, dTTP and dGTP), 0.5 unit of Taq DNA polymerase (Promega, USA), 10.0 pmol µl-1 of primer and 2 ng µl-1 of genomic DNA. The thermal cycling profile was programmed at 94 °C for 2 min as the initial denaturation step, 39 cycles of 30 s at 94 °C, 30 s at 59.1 °C for (GTG)6, at 57.8 °C for 335

(AC)10 and at 56.6 °C for (AG)10, 1 min at 72 °C and final extension step at 72 °C for 10 min. The PCR products were then examined on 1.5% agarose gel and 1 kb DNA ladder (Promega, USA) was ran simultaneously. The gel was documented using the Geliance 200 Imaging System (Perkin Elmer, USA). The DNA bands produced at different loci were determined and named for each DNA sample. Banding profiles generated were converted into binary data matrices on the basis of present (1) or absent (0) of bands. Data scoring was based on several criteria: (1) locus was assumed as independent or non-allelic, (2) there was no bias in scoring monomorphic fragments versus polymorphic fragments, (3) amplified loci were expected to be in the range of 500 bp to 2000 bp, and (4) the similarity of fragment size was assumed to be the indicator of homology. POPGENE Version 1.32 software was used to estimate the genetic diversity of N. cadamba using Shannon’s diversity Index, H’=-∑piln pi-[(S-1/(2N)] where pi stands for frequency of the ISSR fragment and the second fraction of the formula is the correction factor (Lewinton 1972). Clustering of all kelampayan samples was also performed based on shared allele distance, DSA (Chakraborty & Jin 1993) using PowerMarker Version 3.25 and neighbour-joining tree was constructed using MEGA Version 4 (Tamura et al. 2007).

Results and Discussion The genetic diversity of kelampayan was successfully determined using dominant DNA markers based on inter-simple sequence repeats. Three ISSR primers produced 64 loci across all the 247 individuals ranging from 500 bp to 2.0 kb in size. Out of 64 loci, 29 loci were generated by (GTG)6 primer, 16 loci by (AC)10 and 19 loci by (AG)10. The genetic diversity of kelampayan is summarized in Table 2. The percentages of polymorphic loci were in the range of 45.3–74.6% with an average of 58.2%. Shannon’s diversity indices (I) ranged 0.268–0.350. It was found that the genetic diversity of kelampayan in the present study was higher compared with its natural populations, with Shannon’s diversity indices ranging 0.154–0.235 and polymorphic loci ranging 41.3–59.4% (Tiong et al. 2010). Kelampayan progenies from Block 1 (0.350) were found to be more diverse compared with Blocks 2 and 3. Low percentages of polymorphic loci indicate that the genetic diversity of kelampayan is low when compared with other plant species such as Tectona grandis plus tree with 95.5% polymorphic loci (Narayanan et al. 2007), 100% in Asparagus acutifolius (Sica et al. 2005), 85.7% in Swertia chirayita (Joshi & Dhawan 2007) and Glycyrrhiza uralensis with 92.2% polymorphic loci (Yao et al. 2008) amplified by ISSR primers. Table 2 Summary of mean Shannon’s diversity indeices and percentage of polymorphic loci in three different blocks Block

Shannon’s diversity index (I)

Percentage of polymorphic loci (%)

Block 1

0.350

74.6

Block 2

0.270

45.3

Block 3

0.268

54.7

Mean

0.296

58.2 336

Kelampayan progenies from the three different blocks were grouped into 6, 4 and 5 clusters for Blocks 1, 2 and 3 respectively (Figure 2). Each cluster consists of kelampayan progenies that originated from the selected mother trees in Pasai Bon, Niah and Lawas seed production areas (SPAs). In other words, all kelampayan progenies were genetically closely related. A study carried out by Tiong et al. (2010) demonstrated that the coefficient of population differentiation of six natural forests of kelampayan was low (Gst = 0.2013) compared with other species such as Ceriops tagal (Gst = 0.529) (Ge & Sun 2001), C. decandra (0.882) (Tan et al. 2005), Hagenia abyssinica (Gst = 0.25), (Feyissa et al. 2007) and Taxus fauna (Gst = 0.5842) (Shah et al. 2008). The low genetic differentiation value was also observed in other species such as Shorea leprosula, Gst = 0.085 (Lee et al. 2000), Larix potaninii, Gst = 0.116 (Yu et al. 2006) and Calocedrus macrolepis, Gst = 0.042 (Wang et al. 2003). Furthermore, relatively low genetic diversity or high level of genetic differentiation was also recorded in several conifers (Wang et al. 2003).

L3 L2T16 T9

L2 T6

L6T4 L2T5 8

3

T1

L2

T7

L2

T3 9 L5 L7T

L3 T3

L5 T1 6

L2 T4

L2 L4T 8 T9

T6 L1

T5

L9T11

L8

L5T

L6T2

L2T10

T7 L1 L1T2

T3 L2

L8T7

L1T 9

(a)

0

T1

L9

T6

2 T1 L6

L6

L7T

14

L2T1

L7T6 L4T9

L6T14

L8T6 T 13 L3T4 L7T10

L3T13

L7T3

L4T3

T2 L8

L7

L9T12

L6T3 L2T11

L9T15 L1T11

4 T6 T1 L4 L3

0 L4T1 5 L3T1

T8 L2

L4T12

L1 T6

L7

T4

2 12 T1 7T L8 L

T6

1

L4T1

L8T

5

L9

7 L7T 5 16 T1 T L6 L8

8 6T11 LL6T L6T9 L6T7 L 6 L T T1 6T 14 3 10

L2 L 4 T1 T1 5 2

L9T2

T9

T3 L3 1 L3T

L8T7

L3

T7 L3

13

10 L4T4T1 L

L4T12 L7T2 L7T1

L3T8

L3T L4T2

L7T3

L8T4 L8T L8T15 8

L8T14 5 L9T 1 L7T1 4 T L9

12 L6T T5 L6 2 3 L6T L5TL5T10T2 L6 L5T7 6 L4T13 L4T1 L4T3 L5T4 L 5T 6 L4T1 4 L3 L3TL3T11 T1 10 2

L 2T 10 L2T4 L2T5 L2T3 L2T1 L2T2 L1T12 L3T14 16 L3T15 L3T

L6

L6T 11 L4T1 L6T 0 4 L6 L4 T7 T1 4

T5 L6 T7 10 L9 T 6 L 9 L4T 2 T 4 L L6T9 L6T8 L4T1 L7T13 L3T1 L1T7 L1T L8T2 L3T 2 9

L1T7

L2 T1 L2T T7 1 8 L2 T1 L2T31 6

L8T L8 3 T 7 L8 T1 1

L5 L8 T10 T6 L1 T L8T L7T L9T144 3 2 L3T1 L2T 4 10 L7 L7T14 T5 L9T8 L2T11

L1T2 L1T3

L2

T6

L1T1 L7T3 0 L1T8 L2T3 L2T13 L3T7 L5T14 L5T7 L8T4 6 L5T1 4 13 L 5 T L 9T 9 5 T 6 1 7 T1 L L7T L3

L2

T2 L3 9 5 T L95T13 L5T L 12 L3T T4 L3 L2T1 L9T10 L4T8 L8T1 L1T12 L7T1 L8T13 L6T1 6

L1T 9

T6

L1T5 L1T8

L1

L1 T1 L9 1 L3TT1 L1T56 L2T14

5 T1 L5

(c)

T4 L1

L3 L L4T T11 L9T 4T4 12 6

L3T5

T1

L1T161T9 L

(b)

T5

L3

T9 L8 15 T2 L7T L8 13 L 7T 6 L7T1 T6 L7 L7T11

L1

L5T4 L5T1 2

L3T5

13

L4 T1 5

L1 T8

L9 T3

L2 T1 5

L8 T1 2

T1

L1T12

T5

L8

L5T 1 L5T9 5 L4T1 6 L5T6

L7

L4T

1 L6T

L8 T1 3 T4 L5 T1 1

L6T8 L8T8 11 L7T 6

L1

L9T

L3T12

0.02

0.02

Figure 2 Neighbour-joining trees showing relationship among kelampayan trees by using 64 ISSR loci (a) Block 1, (b) Block 2 and (c) Block 3 337

The low level of genetic differentiation among populations of kelampayan may be due to long-distance gene flow that has occurred over fairly large geographical areas (224–519 km) (Table 3). As explained by Hamrick et al. (1992), long-lived woody species with large geographic ranges, outcrossing breeding systems, and wind- or animal-dispersed seeds typically display less variation among populations. Further studies are needed because there is lack of evidence regarding the reproductive biology as well as pollen and seed dispersal mechanisms of kelampayan. Table 3 Geographical distances between three different kelampayan seed production areas in Sarawak Pasai Bon Pasai Bon Niah

295 km

Lawas

519 km

Niah

Lawas

295 km

519 km 224 km

224 km

Conclusion In this study, ISSR analysis was proven as a powerful tool for assessing genetic diversity of kelampayan progenies collected from the three progeny trial blocks at the Landeh Nature Reserve, Semengok, Sarawak. This preliminary information will form the base for kelampayan tree improvement and conservation programmes.

Acknowledgments The authors would like to thank all the laboratory assistants and foresters involved in this research programme for their excellent field assistance in sample collection. This work was part of the joint Industry-University Partnership Programme, a research programme funded by the Sarawak Forestry Corporation (SFC) and Universiti Malaysia Sarawak (UNIMAS).

References Ammiraju, J.S.S., Dholakia, B.B., Santra, D.K., Singh, H., Lagu, M.D., Tamhankar, S.A., Dhaliwal, H.S., Rao, V.S., Gupta, V.S. & Ranjekar, P.K. 2001. Identification of inter-simple sequence repeat (ISSR) markers associated with seed size in wheat. Theoretical and Applied Genetics 102: 726–732. Balasaravanan, T., Chezhian, P., Kamalakannan, R., Ghosh, M., Yasodha, R., Varghese, M. & Gurumurthi, K. 2005. Determination of inter- and intra-species genetic relationships among six Eucalytus species based on inter-simple sequence repeats (ISSR). Tree Physiology 25: 1295–1302. Borneo Post Online. 2010. 25 percent of planted forest target achieved. Retrieved on March 22, 2010 from http://www.theborneopost.com/?p=7718 Chakraborty, R. & Jin, L. 1993. Determination of relatedness between individuals by DNA fingerprinting. Human Biology 65: 875–895. 338

Chezhian, P., Yasodha, R. & Ghosh, M. 2010. Genetic diversity analysis in a seed orchard of Eucalyptus tereticornis. New Forests 40: 85–99. Dogan, B., Duran, A. & Hakki, E.E. 2007. Phylogenetic analysis of Jurinea (Asteraceae) species from Turkey based on ISSR amplification. Annales Botanici Fennici 44: 353–358. Doyle, J.J. & Doyle, J.L. 1990. A rapid total DNA preparation procedure for fresh plant Tissue. Focus 12: 13–15. FAO 2010. Planted forests in sustainable forest management—a statement of principles. Retrieved July 31, 2010 from Forestry Department, Food and Agriculture Organization of the United Nations website: http://www.fao.org/ forestry/plantedforests Feyissa, T., Nybom, H., Bartish, I.V. & Welander, M. 2007. Analysis of genetic diversity in the endangered tropical tree species Hagenia abyssinica using ISSR markers. Genetic Resource Crop Evolution 54: 947–958. Galvan, M. Bornet, B., Balatti, P.A. & Branchard, M. 2003. Inter-simple sequence repeat (ISSR) markers as a tool for the assessment of both genetic diversity and gene pool origin in common bean (Phaseolus vulgaris L.). Euphytica 132: 297–301. Ge, J.X. & Sun, M. 2001. Population genetic structure of Ceriops tagal (Rhizophoraceae) in Thailand and China. Wetlands Ecology and Management 9: 203–209. Godwin, I.D., Aitken, E.A. & Smith, L.W. 1997. Application of inter-simple sequence repeat (ISSR) markers to plant genetics. Electrophoresis 18(9): 1524–1528. Hamrick, J.L., Murawski, D.A. & Nason, J.D. 1992. Factors influencing levels of genetic diversity in woody plant species. New Forests 6: 95–124. Joker, D. 2000. Neolamarckia cadamba (Roxb.) Bosser. Seed Leaflet 6: 11–45. Joshi, P. & Dhawan, V. 2007. Assessment of genetic fidelity of micropropagated Swertia chirayita plantlets by ISSR marker assay. Biologia Plantarum 51(1): 22–26. Kumar, P., Gupta, V.K., Misra, A.K., Modi, D.R. & Pandey, B.K. 2009. Potential of molecular markers in plant biotechnology. Plant Omics Journal 2(4): 141– 162. Lee, S.L., Wickneswari, R., Mahani, M.C. & Zakri, A.H. 2000. Genetic diversity of a tropical tree species, Shorea leprosula Miq. (Dipterocarpaceae), in Malaysia: implications for conservation of genetic resources and tree improvement. Biotropica 32(2): 213–224. Lewinton, R.C. 1972. The apportionment of human diversity. Evol. Biol. 6: 381–398. Mondal, S., Dash, G.K. & Acharyya, S. 2009. Analgesics, anti-inflammatory and antipyretic studies of Neolamarckia cadamba barks. Journal of Pharmacy Research 2(6): 1133–1136. Moreno, S., Pedro, J. & Ortiz, J.M. 1998. Inter-simple sequence repeats PCR for characterization of closely related grapevine germplasm. Euphytica 101: 117–125. Nagaraju, J., Reddy, K.D., Nagaraja, G.M. & Sethuraman, B.N. 2001. Comparison of multilocus RFLPs and PCR-based marker systems for genetic analysis of the silkworm, Bombyx mori. Heredity 86: 588–597. Narayanan, C., Wali, S.A., Shukla, N., Kumar, R., Mandal, A.K. & Ansari, S.A. 2007. RAPD and ISSR markers for molecular characterization of teak (Tectona grandis) plus trees. Journal of Tropical Forest Science 19(4): 218–225. 339

Okun, D.O., Kenya, E.U., Oballa, P.O., Odee, D.W. & Muluvi, G.M. 2007. Analysis of genetic diversity in Eucalyptus grandis (Hill ex Maiden) seed sources using inter-simple sequence repeats (ISSR) molecular markers. Proceedings of African Crop Science Conference 8: 661–665. Patel, D. & Kumar, V. 2008. Pharmacognostical studies of Neolamarckia cadamba (roxb.) Bosser leaf. International Journal of Green Pharmacy 2(1): 26–27. Tamura, K., Dudley, J., Nei, M. & Kumar, S. 2007. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24: 1596–159. Tan, F.X., Huang, Y.L., Ge, X.J., Su, G.H., Ni, X.W. & Shi, S.H. 2005. Population genetic structure and conservation implications of Ceriops decandra in the Malay Peninsula and North Australia. Aquatic Botany 81: 175–188. Tiong, S.Y., Chew, S.F., Ho, W.S. & Julaihi, A. 2010. Genetic diversity of kelampayan (Neolamarckia cadamba) in Sarawak using ISSR markers. In Proceedings of the 3rd Biotechnology Colloquium 2010: How Far Have We Gone? Department of Molecular Biology, Faculty of Resource Science and Technology, UNIMAS. 24–25 April 2010. Semagn, K., Bjornstad, A. & Ndjiondjop, M.N. 2006. An overview of molecular marker methods for plants. African Journal of Biotechnology 5(25): 2540–2568. Shah, A., Li, D.Z., Gao, L.M., Li, H.T. & Moller, M. 2008. Genetic diversity within and among populations of the endangered species Taxus fauna (Taxaceae) from Pakistan and implication for its conservation. Biochemical Systematics and Ecology 36: 183–193. Sica, M., Gamba, G., Montieri, S., Gaudio, L. & Aceto, S. 2005. ISSR markers show differentiation among Italian population of Asparagus acutifolius L. BioMed Central 6: 17. Wang, D.L., Li, Z.C., Hao, G., Chiang, T.Y. & Ge, X.J. 2003. Genetic diversity of Calocedrus macrolepis (Cupressaceae) in southwestern China. Biochemical Systematics and Ecology 32: 797–807. Wong, S.C., Yiu, P.H., Bong, S.T. W., Lee, H.H., Neoh, P.N.P. & Rajan, A. 2009. Analysis of Sarawak bario rice diversity using microsatellite markers. American Journal of Agricultural and Biological Sciences 4(4): 298–304. Yao, H., Zhao, Y., Chen, D.F., Chen, J.K. & Zhou, T.S. 2008. ISSR primer screening and preliminary evaluation of genetic diversity in wild populations of Glycyrrhiza uralensis. Biologia Planatarium 52(1): 117–120. Yu, X.M., Zhou,Q., Qian, Z.Q., Li, S. & Zhao, G.F. 2006. Analysis of genetic diversity and population differentiation of Larix potaninii var. chinensis using microsatellite DNA. Biochemical Genetics 44(11–12): 491–501. Zietkiewicz, E., Rafalski, A. & Labuda, D. 1994. Genome fingerprinting by simple sequence Repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 20(2): 176–183.

340