Analysis of Genetic Diversity of Neem Using RAPD Markers

The present investigation was carried out during 2013-14. RAPD molecular markers were used to evaluate the genetic diversity in populations of Azadirachta indica A. Juss (neem) from different regions of north Karnataka. Out of the 05 random decamer primers used, 03 yielded polymorphic banding patterns. In total, 23 different DNA bands were reproducibly obtained, out of which 18 (78.26%) were polymorphic. The polymorphisms were scored and data was plotted as a binary matrix to be used for phylogenetic analysis using PyElph 1.4 and PAST software. Cluster analysis based on Jaccard’s similarity coefficient using UPGMA grouped all the 20 populations into two major groups. The principal component analysis observed on an average 37% of the total variation. Keyword: Azadirachta indica, genetic variation, DNA, RAPD, PCR. Research Article Cite this article: Bhandare, P. P., Waghmare, T. E., & Naik, G. R. (2021). Analysis of Genetic Diversity of Neem Using RAPD Markers, Ind. J. Pure App. Biosci. 9(2), 66-76. doi: http://dx.doi.org/10.18782/25822845.7658 This article is published under the terms of the Creative Commons Attribution License 4.0. Bhandare et al. Ind. J. Pure App. Biosci. (2021) 9(2), 66-76 ISSN: 2582 – 2845 Copyright © March-April, 2021; IJPAB 67 RAPD markers have successfully been used to estimate levels of relatedness among the individuals, studies of mating systems, and seed dispersal and seedling establishment in natural populations. The objective of this study is to evaluate the genetic variations and genetic relationships among the neem accessions collected from four ecological regions on the basis of RAPD. MATERIALS AND METHODS Plant material Twenty seed samples of Azadirachta indica were collected from four different regions, namely Bidar, Gulbarga, Raichur from Hyderabad-Karnataka (HK) region and Zaheerabad, Telangana. All the populations were studied for the genetic variation and genetic relationship by RAPD. Isolation of genomic DNA 200mg of Seeds were ground to powder with liquid nitrogen using a mortar and pestle. It was suspended in CTAB (Cethyltrimethyl ammonium bromide) extraction buffer modified method of Rahmatollah et al. 2014. The extraction buffer contained (2% CTAB, 5 M NaCl, 0.5 M EDTA, 1 M Tris-HCl, pH 8.0, and 0.1% 2-mercaptoethanol, 0.1% activated charcoal). The sample was transferred to the Eppendorf tube and vertexes at 5-10 times. Then tubes were centrifuged at 10,000 rpm for 15 min. Supernatant was transferred to the new Eppendorf tubes and 800 μl of chloroform: isoamylalcohol (24:1) was added. The content were mixed gently by vertexing and centrifuged at 10,000 rpm for 15 min. The aqueous layer was collected and precipitated by adding 0.8 volumes of cold isopropanol. Mixed by inverting subsequently, the tubes were incubated at -20 0 C for 1 hour. The DNA was pelleted by centrifugation at 10,000 rpm for 15 min. The supernatants were removed and the DNA pellet washed with 700 μl of 70% ethanol, repeating this wash for two times. The pellet was dried for 10 minutes and then dissolved in 200 μl of TE buffer. After DNA extraction, the existence and concentration of DNA were checked by agarose gel electrophoresis. Purification of Genomic DNA by RNase and Proteinase treatment The isolated genomic DNA was further purified by treating with RNase to remove RNA. To 100 μl of isolated DNA sample, 5 μl of RNase (5 mg/ml) was added and incubated at 37 0 C for 1 h. After incubation with RNase, Proteinase K treatment was done to make the DNA free from proteins. To the RNase treated 100 μl DNA, 0.25 μl of Proteinase K (20 mg/ml) was added and incubated at 37 0 C for 30 min. The reaction mixture was made up to 500 μl with TE buffer and extracted with saturated phenol, chloroform, isoamyl alcohol treatment. Analysis of RAPD RAPD-PCR analysis of isolated DNA was performed by five random decamer primes, which were procured from Eurofins Genomics India Pvt Ltd., Bengaluru. The primers and their sequences are listed in Table 1. The molecular biology chemicals required for RAPD-PCR analysis are Taq DNA polymerase, Taq DNA assay buffer (100 mM Tris-HCl (pH 9.0), 500 mM KCl, 15 mM MgCl2 and 0.1% gelatin), dNTPs (dATP, dCTP, dGTP and dTTP) were purchased from Bangalore Genei Pvt. Ltd., Bangalore. The DNA isolated from 20 neem ecotypes were diluted to a concentration of 20 ng/μl, RAPDPCR amplifications were performed in 20 μl reaction mixture, as described by Dhillon et al. (2007). The reaction mixture contained following components in 0.5 ml PCR tube (Table 2). The reaction mixture was thoroughly mixed and centrifuged the samples briefly to bring down the contents of the tube. To this 50 μl of mineral oil was overlaid to prevent the evaporation of the contents during the amplification process. The DNA thermal cycler Corbett was used with following standardized program (Table 3). The RAPD-PCR amplified products were visualized by ethidium bromide and banding pattern were stored for further analysis under gel documentation system and photographs were taken. Bhandare et al. Ind. J. Pure App. Biosci. (2021) 9(2), 66-76 ISSN: 2582 – 2845 Copyright © March-April, 2021; IJPAB 68 Data Analysis For data analysis of RAPD-PCR, only distinct visible and reproducible bands were manually scored, „1‟ for presence and „0‟ for the absence and the binary data were used for statistical analysis by using the software “PAST”. The sizes of the fragments (molecular weight in base pairs) were estimated by using 100 bp ladder markers, which was run along with the amplified products. A genetic dissimilarity matrix was calculated according to Jaccard„s similarity index (JSI) which estimated all pairwise differences in the amplification product and cluster analysis was done by the UPGMA method using the variance algorithm in PhyElp 1.4 software. Table 1: List of RAPD primers used for the amplification of genomic DNA of neem seed Primer code Nucleotide sequence (5’→3’)


INTRODUCTION
Various types of molecular markers are utilized to evaluate DNA polymorphism. RAPD a PCR based technique is commonly used by plant biologists to perform a number of tasks, including the genetic fingerprinting of plant varieties, determining similarities among inbred varieties, mapping of plant genomes, and establishing phylogeny among plant species. Although, RAPD is criticized for its low reproducibility (Hansen et al., 1998& Virk et al., 2000, it is overcome by optimization of the RAPD reaction and conditions. Large scale characterization of plant species in varying geoclimatic conditions can be performed using various parameters such as seed morphometric traits and isozymes (Kundu, 1999). However, environmental factors as well as the developmental stages of the plants influence such traits. DNA based markers provide new tool for ecological and genetic studies of evolutionary processes (Cruzan, 1998).

Purification of Genomic DNA by RNase and Proteinase treatment
The isolated genomic DNA was further purified by treating with RNase to remove RNA. To 100 μl of isolated DNA sample, 5 µl of RNase (5 mg/ml) was added and incubated at 37 0 C for 1 h. After incubation with RNase, Proteinase K treatment was done to make the DNA free from proteins. To the RNase treated 100 μl DNA, 0.25 μl of Proteinase K (20 mg/ml) was added and incubated at 37 0 C for 30 min. The reaction mixture was made up to 500 μl with TE buffer and extracted with saturated phenol, chloroform, isoamyl alcohol treatment. Analysis of RAPD RAPD-PCR analysis of isolated DNA was performed by five random decamer primes, which were procured from Eurofins Genomics India Pvt Ltd., Bengaluru. The primers and their sequences are listed in Table 1. The molecular biology chemicals required for RAPD-PCR analysis are Taq DNA polymerase, Taq DNA assay buffer (100 mM Tris-HCl (pH 9.0), 500 mM KCl, 15 mM MgCl 2 and 0.1% gelatin), dNTPs (dATP, dCTP, dGTP and dTTP) were purchased from Bangalore Genei Pvt. Ltd., Bangalore. The DNA isolated from 20 neem ecotypes were diluted to a concentration of 20 ng/μl, RAPD-PCR amplifications were performed in 20 μl reaction mixture, as described by Dhillon et al. (2007). The reaction mixture contained following components in 0.5 ml PCR tube ( Table 2). The reaction mixture was thoroughly mixed and centrifuged the samples briefly to bring down the contents of the tube. To this 50 μl of mineral oil was overlaid to prevent the evaporation of the contents during the amplification process. The DNA thermal cycler Corbett was used with following standardized program (Table 3).
The RAPD-PCR amplified products were visualized by ethidium bromide and banding pattern were stored for further analysis under gel documentation system and photographs were taken.

Data Analysis
For data analysis of RAPD-PCR, only distinct visible and reproducible bands were manually scored, "1" for presence and "0" for the absence and the binary data were used for statistical analysis by using the software "PAST". The sizes of the fragments (molecular weight in base pairs) were estimated by using 100 bp ladder markers, which was run along with the amplified products. A genetic dissimilarity matrix was calculated according to Jaccard"s similarity index (JSI) which estimated all pairwise differences in the amplification product and cluster analysis was done by the UPGMA method using the variance algorithm in PhyElp 1.4 software.

Isolation of genomic DNA from neem seeds
The isolated DNA was subjected to check their purity at A 260 /A 280 which were found 1.91 that indicated the isolated DNA was contaminated with RNA which were further purified by RNase treatment. The DNA yielded from the neem was 8-10 μg/g of seed material. The electrophoresis pattern of isolated genomic DNA after RNase treatment is shown in the Fig. 1.

Genetic diversity analysis by RAPD Marker
In the present study 20 neem ecotypes (Table  4) were subjected for RAPD analysis with 5 random decamer primers, out of which 3 primers were produced considerable polymorphism while 2 primers were either not able to produce the bands or had weaker amplification. The high intensity and reproducible bands were only selected for polymorphism analysis among the ecotypes. Differing band intensities were not taken into account to avoid errors. Five primers used OPA-5, OPC-8, OPA-19, OPA-03 and OPQ-17. No bands were generated in the primer OPC-8 and OPA-19. Primer OPA-5 gives 8 different bands, in which 6 are polymorphic band. Primer OPA-03 gives 7 different bands, in which 5 are polymorphic band. Primer OPQ-17 gives 8 different bands. Out of these 7 bands are polymorphic band. The 5 primers generated different bands with size ranging from 200-2000 bp. No PCR product was generated through the primer OPC-8 and OPA-19, indicating that no binding sites for this primer. Out of the 23 bands, 18 Scorable bands (78.26%) were found to be polymorphic. These results gave an average of polymorphic bands per primer was shown in Table 5. The gel images were analyzed for the presence and absence of bands and this data was plotted as a binary matrix to be used for phylogenetic analysis using PyElph 1.4 and PAST software. PyElph 1.4 and PAST (PAleontological STatistics) free software tool, runs on standard Windows computers.

DISCUSSION
Ecological and geographical differentiations are important factors, influencing strategies for breeding and sampling tree crops. With the objective of selecting and maintaining elite trees to bring about an overall genetic improvement in neem for enhanced output of oil, provenances of neem have been collected from four regions in which three regions are from Karnataka and one is from Telangana state and are at present being maintained at Gulbarga University campus, Kalaburagi. Geographically isolated population accumulates genetic differences as they adapt to different environmental conditions (Sarwat et al., 2008). Genetic diversity analysis based on morphological and biochemical parameters has plenty of limitations due to the influence of various environmental factors. Molecular markers based genetic diversity analysis is beneficial because they are unaffected by environmental conditions. Molecular studies using RAPD primers have been used to examine genetic variation within a species, to analyse particular genotypes, for cultivar identification and to study the clonal structure of several tree species (Kresovich et al., 1992;Vierling & Nguyen, 1992;& Hsaio & Rieseberg, 1994). High variations in the RAPD patterns have also been reported in the black and red spruce varieties, where the variation is attributed to the interaction between the plantlets and the microenvironment (Nkongolo et al., 1998). The pairwise comparison of the RAPD profiles based on both shared and unique amplification products was made to generate a similarity matrix. Similarity indices estimated on the basis of three primers ranged from 0.00 to 1.00. Such a wide range of similarity coefficient values suggest that the neem germplasm collection represents a genetically diverse population. The similar data have been reported by Singh et al. (1999) based on their study using amplified fragment length polymorphism (AFLP). One of the major contributory factors to the high degree of polymorphism observed in neem may be on account of its evolutionary status as an outcrossing angiosperm (Kundu, 1999;& Dhillon et al., 2005). The high diversity revealed by RAPD is that out-breeding woody plants retain considerable variability (Hamrick, 1990). This was further supported by molecular marker studies (Hamrick, 1989). The maintenance of a high genetic variance within populations was favored by the genetic system of the species like effective gene flow, out breeding, mutation, high genetic load, etc (Matyas, 1996).
The RAPD study revealed that there is a genetic diversity in neem. On average, 78.94% polymorphism was detected with each RAPD fingerprint, indicating a high marker index. RAPD analysis revealed that the genetic similarity varied from a high of 100% (1.00) to a low of 0% (0.00) within the 20 neem ecotypes. Such a wide range in similarity coefficient values suggests that the neem germplasm collection represents a genetically diverse population. One of the major contributory factors to the high degree of polymorphism observed in neem may be on account of its evolutionary status as an outcrossing angiosperm. Similar observations have been made with coconut by Perera et al. (1998) where the level of genetic diversity was shown to correlate with the breeding nature of the palms. The tall out-breeding variety exhibited a higher diversity than the inbreeding types.
Cluster analysis has clearly indicated that there are two major groups based on the similarity, values classified in neem population. Genetic relationships among neem accessions were also visualized by performing principal component analysis based on RAPD data. The principal component analysis observed on an average 37% of the total variation. This result is coherent with the dendrogram generated by employing UPGMA and is a further confirmation of the genetic similarities. Cluster analysis of the neem genotypes employing UPGMA as well as principal component analysis led to the segregation of the accessions into the distinct groups, which reflected their geographical distribution. Ecological and geographical differentiation are important factors, which influence breeding and sampling strategies of tree crops (Namkoong, 1986), which further help in understanding the population structure. The reason for the grouping of samples to one cluster collected from different sub-zones may be due to human intervention, which makes partitioning and distribution of variability complex. This assumption has been further supported by the Ranade et al. (2002), that neem has distributed throughout India from only a few groups of founder population.
Variation in genetic diversity within species is usually related with geographic range, mode of reproduction, mating systems, seed dispersal and fecundity (Loveless, 1992). Muller et al. (1992) in their review of the results of studies on genetic variation in various coniferous and angiosperm tree species observed that species with geographically distinct ranges tend to show a moderate too high-inter and intra-population genetic variation. High genetic diversity detected in the present study may be due to all these prevalent background factors as populations of CPTs studied are widely distributed in different eco-geographic regions. Similar conclusions were reached by Singh et al. (2002) while assessing genetic diversity in 37 neem accessions from different agroclimatic regions of India using AFLP and SAMPL molecular markers. They reported that neem germplasm within India constitutes a broad genetic base. The RAPD banding pattern in ecotypes of HK area also show wide genetic variation. This has resulted in fluctuation yield attributes of seeds and seed oil etc. It is necessary to identify the candidate plus trees with highly productivity genotypes and mass cultivation on wastelands for future economic and industrial development.