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