E-ISSN 2549-8703 I P-ISSN 2302-7282 BIOTROPIKA Journal of Tropical Biology https://biotropika.ub.ac.id/ Vol. 9 | No. 2 | 2021 | DOI: 10.21776/ub.biotropika.2021.009.02.09 GENETIC DIVERSITY OF imbricatus AT BUKIT TAPAK, TABANAN, BALI BASED ON RAPD MARKER

KERAGAMAN GENETIK Dacrycarpus imbricatus DI BUKIT TAPAK, TABANAN, BALI BERDASARKAN PENANDA RAPD

Made Pharmawati 1)*, I Made Anom Sutrisna Wijaya2)

Received : May, 29 2021 ABSTRACT Dacrycarpus imbricatus Blume is a member of the family. In Bali, D. Accepted : July, 13 2021 imbricatus was found in Bukit Tapak, Tabanan Regency. This species is one of the dominant species in Bukit Tapak. This study aimed to determine the genetic variation of D. imbricatus in Bukit Tapak using molecular markers RAPD (Random Amplified Polymorphic DNA). The genetic diversity of D. imbricatus needs to be studied to obtain Authors affiliation: the information used for the conservation of this species. Leaf samples were taken from 1) Biology Department, Faculty of Bukit Tapak, Candikuning, Baturiti, Tabanan Regency, Bali. DNA was extracted using Mathematics and Natural the CTAB method, followed by extraction using chloroform: isoamyl alcohol. DNA Sciences, Universitas Udayana, precipitation was carried out using ethanol. RAPD analysis was performed using Kampus Bukit Jimbaran, Bali polymerase chain reaction (PCR) using four primers. PCR products were visualized 2) Department of Agricultural and using agarose gel electrophoresis and ethidium bromide staining. The results showed that Biosystems Engineering, the amplified DNA bands ranged from 1 to 5 bands with DNA band sizes ranging from Faculty of Agricultural 230 bp - 422 bp. Only OPA4 and UBC106 primers can be used to detect D. imbricatus Technology, Universitas diversity based on the H, D, and R values. The detected genetic variation is low, as Udayana, Kampus Bukit indicated by an average polymorphism of 32.5% and similarities between samples Jimbaran, Bali ranging from 0.51 to 1.

Keywords: Bukit Tapak, Dacrycarpus imbricatus, polymorphism, RAPD Correspondence email: ABSTRAK *[email protected] Dacrycarpus imbricatus Blume adalah anggota dari famili Podocarpaceae. Di Bali, D. imbricatus ditemukan di Bukit Tapak, Kabupaten Tabanan. Spesies ini merupakan salah satu spesies yang dominan di Bukit Tapak. Penelitian ini bertujuan untuk menentukan variasi genetik D. imbricatus di Bukit Tapak menggunakan penanda molekuler RAPD (Random Amplified Polymorphic DNA). Keragaman genetik D. imbricatus perlu dipelajari untuk mendapatkan informasi yang digunakan untuk konservasi spesies ini. Sampel daun diambil dari Bukit Tapak, Candikuning, Baturiti, Kabupaten Tabanan, Bali.

Ekstraksi DNA menggunakan metode CTAB dilanjutkan dengan menggunakan

kloroform: isoamylalkohol. Pengendapan DNA dilakukan menggunakan etanol. Analisis

RAPD dilakukan dengan menggunakan polymerase chain reaction menggunakan empat

primer. Produk PCR diamati menggunakan elektroforesis gel agarosa dan pewarnaan How to cite: etidium bromida. Hasil penelitian menunjukkan bahwa pita DNA yang teramplifikasi Pharmawati, M, IMAS Wijaya. berkisar 1 – 5 pita dengan ukuran pita DNA berkisar 230 bp – 422 bp. Keragaman D. 2021. Genetic diversity of imbricatus hanya terdeteksi bila menggunakan primer OPA4 dan UBC106 berdasarkan nilai H, D dan H. Variasi genetik yang terdeteksi tergolong rendah yang ditunjukkan Dacrycarpus imbricatus at Bukit dengan rata-rata polimorfisme 32,5% dan similarity antar sampel berkisar 0.51 sampai Tapak, Tabanan, Bali based on 1. RAPD marker. Journal of Tropical Biology 9 (2): 157-161. Kata kunci: Bukit Tapak, Dacrycarpus imbricatus, polimorfisme, RAPD

INTRODUCTION In Bali, the habitat of D. imbricatus is in Bukit Tapak, Tabanan Regency. Bukit Tapak Forest is Dacrycarpus imbricatus Blume or also known located at an altitude of 1520-1600 m asl [2]. This as Podocarpus imbricatus is a member of in Bali is known as the local name cemara Gymnosperms and belongs to the Podocarpaceae pandak. D. imbricatus is one of the dominant family. This plant is a native species of species in Bukit Tapak. Bukit Tapak forest is and is distributed in Sumatra, Java, Sulawesi, classified as a secondary forest [3]. D. imbricatus Lesser Sunda, and Papua. This species grows in the is a long-lived pioneer species that grew due to past highlands from 500 to 2,000 m asl (above sea level) disturbances in its territory [4]. Conservation of D. [1]. imbricatus is needed to prevent the loss of this species in the area due to species change.

Pharmawati & Wijaya 157 https://biotropika.ub.ac.id/ Information on the genetic diversity of D. imbricatus is important for developing conservation strategies. Preservation of genetic diversity is necessary for the survival of a species over a long period. Loss of genetic diversity will decrease the ability to adapt to environmental changes [5]. Genetic diversity studies can identify alleles that can contribute to an organism's ability to survive in its habitat or survive in a varied habitat. Loss of genetic diversity can lead to extinction in small populations due to decreased environmental adaptability [6, 7]. The detection of genetic diversity can be carried out by several methods. The genetic diversity of has traditionally been observed based on Figure 2. D. imbricatus tree in Bukit Tapak morphological characters or growth responses. In this method, individuals are differentiated based on DNA extraction. The leaves were washed with their phenotypic character, which can be soap and running water until clean, then 0.1 g of influenced by the environment [8]. The leaves were crushed using a mortar and pestle. As development of molecular markers provides a much as 1 ml of CTAB buffer (2% w / v CTAB, sensitive method for detecting diversity at the DNA 1.4 M NaCl, 50 mM EDTA, 100 mM Tris-HCl (pH level. One of the techniques used in DNA 8), 2% v / v 2- mercaptoethanol) was added [10, fingerprinting in evaluating genetic diversity is 11]. Incubation of the mixture was done at 65ºC for RAPD (Random Amplified Polymorphic DNA). 1.5 hours. Following incubation, centrifugation The RAPD technique detects DNA polymorphisms was done at 14,000 rpm for 5 minutes. using single primers with random sequences [9]. The supernatant was collected to a new tube and RAPD analysis is useful for detecting nucleotide added 500 µl of chloroform: isoamyl alcohol (24: polymorphisms that show genetic diversity. This 1). The mixture was vortexed and centrifuged for 5 study is aimed to analyze the genetic variation of minutes at 14,000 rpm. The top layer was removed D. imbricatus in Bukit Tapak, Tabanan, Bali by and transferred to a new tube, and purification using RAPD. The results of this study can be used using chloroform: isoamyl alcohol (24:1) was as a basis for developing a conservation strategy. repeated. The top layer was transferred to a new

tube, and 500 µl of chilled ethanol was added and METHODS incubated overnight. After that, the sample was Sampling of plant materials. D. imbricatus centrifuged at 8,000 rpm for five minutes. The leaves were taken from Bukit Tapak, Tabanan pellets were washed using 70% ethanol and Regency, Bali (Figure 1). The leaves used were centrifuged for 5 minutes at 8,000 rpm. The pellets fully developed leaves. The samples used were 16 were air-dried and resuspended in 100 µl of sterile individuals’ trees of D. imbricatus (Figure 2). H2O [10, 11]. DNA was visualized with agarose gel (1%) electrophoresis in TAE buffer for 45 minutes at 100 volts. DNA staining was done with ethidium bromide. Visualization was performed using UV- transilluminator [12]. The DNA concentration was estimated by comparing with 50 ng and 100 ng of lambda () DNA. RAPD. PCR-RAPD was conducted in 20 µl volume with 25 ng DNA, 1 x PCR buffer, 200 µM dNTP, 3 mM MgCl2, 1 mM primer, and 1 unit taq polymerase. Eight primers from operon primers (OP) and the University of British Columbia (UBC) were screened (Table 1).

Figure 1. Map of the location of Bukit Tapak (arrow) based on Google Earth

158 Biotropika: Journal of Tropical Biology | Vol. 9 No. 2 | 2021 https://biotropika.ub.ac.id/ Table 1. Name of primers screened and their cells are exposed to lysis buffer (CTAB) longer so sequences that a better DNA lysis process occurred [17]. The EDTA concentration in the CTAB lysis Primer Name Sequence (5’-3’) buffer used in this study was also higher than that OPA2 TGCCGAGCTG of Doyle & Doyle [10]. Modifications have been OPA4 AATCGGGCTG made previously using 50 mM EDTA [11]. In the OPB8 GTCCACACGG lysis buffer, EDTA functions as a chelator of Mg2+, OPD11 AGCGCCATTG where Mg2+ is a cofactor of the enzyme nuclease OPD14 CTTCCCCAAG [18]. By increasing the EDTA concentration, it is OPH1 GGTCGGAGAA expected that the amount of Mg2+ was chelated UBC106 CGTCTGCCCG more, thereby inhibiting nuclease activity. UBC250 CGACAGTCCC Diversity analysis with RAPD. DNA amplification by PCR-RAPD was successfully The PCR cycle was 1x initial denaturation at performed using four primers (OPA4, OPD14, 95ºC for 5 minutes, 40 x denaturation at 95ºC for 1 UBC106, and UBC250) of the eight primers tested. minute, primer annealing at 37oC for 1 minute, Only the four primers mentioned above produced extension at 72ºC for 1 minute. The final extension clear PCR products, while the other four primers was carried out for 10 minutes at 72ºC. PCR (OPA2, OPB8, OPD11, and OPH1) produced products were visualized using 1.8% agarose gel products with a very strong background smear so electrophoresis with ethidium bromide staining. that the DNA bands were not clearly observed. Data analysis. Existing DNA bands were score Figure 4. shows the PCR-RAPD results using UBC 1, whereas the absence of DNA band was scored 0. 106 primers with a 1kb ladder as a size marker. Primary efficiency analysis was performed using Primer selection in diversity analysis plays an iMEC [13]. The kinship between individuals was important role. The analysis of the RAPD primer analyzed using NTSYS 2.1 [14], with the UPGMA efficiency is shown in Table 2, while the number method with a simple matrix similarity coefficient. of amplified DNA bands, the size range of DNA bands, and the degree of polymorphism was shown RESULTS AND DISCUSSION in Table 3. DNA extraction. DNA with concentrations between 80 ng/µl to 175 ng/µl was obtained. The resulting DNA has a large size indicated by the position of the DNA band at the top of the gel. In addition to these DNA bands, there is DNA that appears smeared (Figure 3). The smeared DNA indicated that the DNA had been degraded into fragments of various sizes. DNA smearing is caused by nuclease activity during extraction [15, Figure 4. PCR-RAPD products using UBC 106 16]. primer. Numbers 1 – 16 are the individual sample numbers. 1 kb ladder was used as a size marker

Of the four primers used in PCR-RAPD, only two primers (OPA4 and UBC106) produced polymorphic DNA bands, while the other two primers (OPD14 and UBC250, produced monomorphic DNA bands. Figure 3. Results of D. imbricatus DNA Based on the primary efficiency analysis, electrophoresis. Numbers 1 – 16 are the individual OPD14 primer and UBC250 primer did not have sample numbers. The DNA was compared with  the ability to differentiate between individuals. DNA at 50 ng and 100 ng This could be seen from the value of H, D, and R, which are 0 (Table 2). The value of H was the In the early stages of DNA extraction, the probability of an individual being heterozygous for crushed leaves were incubated in a CTAB buffer at a particular locus in the population, the value of D o a temperature of 65 C for 1.5 hours. This is the probability that two individuals randomly incubation period was much longer than the plant have different patterns, while the value of R is the DNA extraction method by Doyle & Doyle [10]. primary ability to distinguish a number of large The long incubation time in the CTAB buffer was genotypes [13]. Primer UBC106 had the highest done because the leaves of D. imbricatus were H, D, and R values indicating that UPB106 is the quite hard. By prolonging the incubation time, the best primer in differentiating individuals.

Pharmawati & Wijaya 159 https://biotropika.ub.ac.id/ Table 2. The efficiency of each primer in the analysis of the diversity of D. imbricatus in Bukit Tapak Primer H D R Name OPA4 0.219 0.238 0.5 OPD14 0 0 0 UBC106 0.447 0.564 1.875 UBC250 0 0 0 Average 32.5 H=heterozygosity index; D=discriminating power; R=resolving power Figure 5. Genetic relationship between individuals D. imbricatus. Number 1 - 16 are the individual Table 3. The number and size range of DNA bands sample number. The coefficient showed in the and the percentage of polymorphisms from PCR figure was the coefficient of similarity products CONCLUSION Primer No. of Band Polymorphism Name DNA Size (%) The population of D. imbricatus in Bukit Tapak, Bands (bp) Tabanan, Bali showed a low level of diversity, OPA4 2 422-250 50 which has an average polymorphism of 32.5% and OPD14 1 415 0 a similarity index between samples between 0.51 UBC106 5 408-230 80 to 1. Further research is needed by increasing the UBC250 2 410-230 0 number of samples, the number of primers, and testing for the occurrence of inbreeding pressure. Average 32.5

ACKNOWLEDGMENT RAPD analysis of D. imbricatus showed an average polymorphism of 32.5%. The low The authors thank Sutomo and Wawan Sujarwo percentage of polymorphic DNA bands indicates from the Bali Botanic Garden for their assistance low genetic diversity and vice versa [19]. Table 3 in sample collection. This study was funded by the supported the analysis shown in Table2. Primer Directorate of Higher Education, Republic UBC106 showed the highest percentage of Indonesia, through a Higher Education Excellence polymorphism. Based on Table 2 and Table 3, Research Grant. primer OPA4 and UBC106 were suitable for detecting the genetic variation of D. imbricatus. REFERENCES The kinship between individuals of the 16 [1] Thomas P (2013) Dacrycarpus imbricatus. samples tested is shown in Figure 5. It can be seen The IUCN Red List of Threatened Species. that the similarity coefficient ranges from 0.51 to e.T42445A2980614. 1. This showed a fairly high level of similarity http://dx.doi.org/10.2305/IUCN.UK.2013- between individuals. Some individuals have a 1.RLTS.T42445A2980614.en. similarity coefficient of 1, namely between individuals 1, 2, 4, and 5, between individuals 3, [2] Fardila D, Sutomo (2011) Efek tepi koridor 7, and 13, between individuals 10 and 11, also jalan di Hutan Bukit Pohen, Cagar Alam between individuals 14, 15, and 16. A high value Batukahu, Bali. Berkala Penelitian.Hayati 17: of the similarity coefficient means the population 9–13. did not have enough diversity [20]. Based on the [3] Sutomo, Erosi NK, Undaharta, Bangun TM, average polymorphism and the tree of kinship, the Lugrayana IN (2012) Studi awal komposisi results of this study indicate that the genetic dan dinamika vegetasi pohon,hutan Gunung diversity of D. imbricatus in Bukit Tapak was low Pohen Cagar Alam Batukahu Bali. Bumi and needs to be improved. Lestari 12(2): 366-381. [4] van Steenis CGGJ (1972) The mountain flora of Java. E. J. Brill, Leiden. [5] van Zonneveld M, Scheldeman X, Escribano P, Viruel MA, van-Damme P, Garcia W, Tapia C, Romero J, Siguen M, Hormaza IJ (2012) Mapping Genetic Diversity of Cherimoya (Annona cherimola Mill.): Application of spatial analysis for

160 Biotropika: Journal of Tropical Biology | Vol. 9 No. 2 | 2021 https://biotropika.ub.ac.id/ conservation and use of plant genetic Technology 33(4): 753-756. resources. PLoS ONE 7(1):e29845. doi.org/10.1590/S0101-20612013000400022. doi:10.1371/journal.pone.0029845. [18] El-Ashram S, Al Nasr I, Suo, X. (2016). [6] Martin MA, Herrera MA, Martín LM (2012) Nucleic acid protocols: Extraction and In situ conservation and landscape genetics in optimization. Biotechnology reports forest species. Journal of Natural Resources (Amsterdam, Netherlands) 12: 33–39. and Development 2: 1-5. doi.org/10.1016/j.btre.2016.10.001. [7] Nonić M., Šijačić-Nikolić M. (2019) Genetic [19] Probojati RT, Wahyudi D, Hapsari L (2019). diversity: Sources, threats, and conservation. Clustering analysis and genome inference of In: Leal Filho W., Azul A., Brandli L., Özuyar pisang raja local cultivars (Musa spp.) from P., Wall T. (Eds) Life on land. Encyclopedia Java Island by Random Amplified of the UN Sustainable Development Goals. Polymorphic DNA (RAPD) Marker. Journal Springer Cham, Switzerland. of Tropical Biodiversity and Biotechnology doi.org/10.1007/978-3-319-71065-5_53-1. 4(2): 42- 53. [8] Choudhary N, Ahuja U, Chawla V, Jain RK, [20] atyawan D, Tasma IM (2011) Genetic Kumari P, Batan KR. (2011) Morphological diversity analysis of Jatropha curcas and molecular variability in weedy rices of provenances using Randomly Amplified Haryana. Asian Journal of Agricultural Polymorphic DNA Markers Jurnal Research. 5: 250-259. AgroBiogen 7(1):47-55. [9] Welsh J, M. McClelland M (1990) Fingerprinting genomes using PCR with arbitrary primers. Nucleic Acids Research 18(24):7213-7218. [10] Doyle JJ, Doyle JL (1990) A rapid total DNA preparation procedure for fresh plant tissue. Focus 12: 13-15. [11] Pharmawati M, Yan G. McFarlane I. J. (2004) Application of RAPD and ISSR markers to analyse molecular relationships in Grevillea (Proteaceae). Australian Systematic Botany 17: 49-61. doi.org/10.1071/SB03016. [12] Sambrook J. Fritsch EF, Maniatis T (1989). Molecular Cloning. A Laboratory Manual.: Cold Spring Harbor, Lab.Press,USA. [13] Amiryousefi A, Hyvönen J, Poczai P. (2018) iMEC: Online marker efficiency calculator. Application in Plant Science 6(6):e01159. doi: 10.1002/aps3.1159. [14] Rohlf F.J (2000). NTSYSpc. Numerical taxonomy and multivariate analysis system. Version 2.1. Exeter Publishing Setauket, New York. [15] Sahu SK, Thangaraj M, Kathiresan K (2012) DNA extraction protocol for plants with high levels of secondary metabolites and polysaccharides without using liquid nitrogen and phenol. ISRN Molecular Biology 2012: 205049,. doi:10.5402/2012/205049. [16] Aboul-Maaty NAF, Oraby HAS (2019). Extraction of high-quality genomic DNA from different plant orders applying a modified CTAB-based method. Bulletin of the National Research Center 43: 25 doi.org/10.1186/s42269-019-0066-1. [17] Jain SA, de Jesus FT, Marchioro GM, de Araújo, ED (2013). Extraction of DNA from honey and its amplification by PCR for botanical identification. Food Science and Pharmawati & Wijaya 161