Sesame (Sesamum indicum L.) family Pedaliaceae, is one of the popular high-value, multipurpose oil seed crop grown widely in tropical and subtropical regions. Production of Sesame in India is encouraging, as more than 30% of world production is contributed by India alone and identified as the first rank in total production. The average yield of sesame on global scale is 5.1 quintals/ha while, current world production is estimated at about 4.04 million tons annually. India placed the second in the world after Myanmar with 18.20 lakh ha and 6.10 lakh tons production respectively. The average yield of sesame on global scale is 5.11 q/ha, while in India, it is 3.30 q/ha which is very low. It is widely used against various diseases including cancer, cold, colic etc.  . Sesame oil contains a unique compound known as lignans. Lignans comprise sesamin, sesamolin, and a small amount of sesamol  . Lignans are also phytoestorgens and their conversion to enterolactone is very important in preventing hormone-dependent cancers like breast and prostate and cardiovascular diseases. Sesamum has also a unique anti-oxidant compound “sesamol” and poly-unsaturated fatty acid    . Even, Sesamum has great economic importance but the farmers do not attracted much for cultivation as the low productivity and low economic return. Even India is self-suffi- cient in production, the productivity is decline due to unavailability of high yielding varieties    . Sesamum species are represented by a large number of genotypes/ accessions which are facing genetic erosion. Genetic diversity analysis is utmost important for breeding program. Morphological, physiological, agronomic traits and genetic analysis through molecular marker are the prerequisite to know the genome organization of wild land races and related cultivated ones. Recently, molecular marker technology helps the identification of plant genetic resources. Moreover, the precise number of cultivars/accessions is still unverified since problems of mislabelling are very often detected. Therefore, it is essential to establish the strategies for the preservation of local genotypes. Till now, very scanty systematic work has been made to identify the indigenous as well as exotic collections along with the local germplasms to study the genetic diversity  . Different PCR and non-PCR based molecular markers have been utilized to identify the varieties/accessions  -  . The information of genetic variation between indigenous and exotic collection of sesame will provide better management strategy for crop improvement programme under adverse climatic conditions   . The combination of both agro-morphological and molecular investigation would be more appropriate to analyse the germplasms of sesame. The present investigation was to identify the phylogenetic analysis of 33 cultivars of Sesamum representing through agro-morphological and molecular markers.
2. Materials and Methods
2.1. Plant Materials
The genetically pure seed material of 33 Sesamum varieties representing different geographical locations in India was collected from the germplasm centre of AICRP in Sesamum at Orissa University of Agriculture & Technology, Bhubaneswar, India (Table 1). The seeds were sown in the pots and grown in green house with 30˚C ± 2˚C with 75% relative humidity to study the morphological, biochemical and molecular analysis. The experimental layout was a randomized block design with three replications. The Soil was mixed with cow-dung at a ratio of 3:1. The plants were irrigated every three days interval. The morphological data were taken every 7 days interval.
2.2. DNA Preparation
DNA was isolated using CTAB method following the protocol of Doyle and Doyle  , with minor modifications. 1.0 - 1.5 g of young non-senescent leaves were ground in liquid nitrogen. Then they were incubated in CTAB buffer (3% w/v CTAB, 100 mM Tris-HCl, 20 mM EDTA, 1.4M NaCl, 2% v/v β-mercaptoethanol, 2% w/v polyvinyl pyrrolidine, pH 8.0) for 2 h at 65˚C. The homogenate was then extracted with an equal volume mixture of chloroform:isoamylalcohol (24:1) and centrifuged at 9000 rpm for 10 min. The upper aqueous layer was recovered and precipitated with pre-chilled isopropanol. The pallet was suspended with Tris-EDTA buffer (pH 8.0). The crude DNA was treated with RNase and incubated for 30 min at 37˚C and again extracted with 1 volume phenol and subsequently with 1 volume of chloroform: isoamylalcohol (24:1). The supernatant were collected and precipitated with 3 M sodium acetate and pre- chilled ethanol. The DNA pellet was washed with 70% ethanol, dried, and re-suspended in TE buffer. The high molecular weight DNA was checked for quality and quantity electrophoretically using 0.8% agarose gel against a known amount of DNA taken as standard.
2.3. Primer Screening
Thirty synthesized ISSR primers (M/S Merck Bioscience, India) were initially screened using six genotypes to determine the suitability of each primer for the study. Primers were selected for further analysis based on their ability to detect distinct, clearly resolved, and polymorphic amplified products within the genotypes. To ensure reproducibility, the primers generating no, weak, or complex patterns were discarded. A few well-amplified fragments that were not reproducible across 2 replicates of DNA extraction were also discarded from analysis.
2.4. PCR Amplifications and Electrophoresis of PCR Products
For ISSR experiment, the initial optimization of PCR was conducted including concentration of template DNA, primer, MgCl2, number of PCR cycles, and annealing temperature. The PCR reaction had a total volume of 25 μl containing 20 ng templates DNA, 100 mM each dNTPs, 20 ng of oligonucleotides synthesized primer (M/S Emerck Bioscience, India), 2.5 mM MgCl2, 1x Taq buffer (10 mM Tris-HCl) pH 9.0, 50 mM KCl, 0.01% gelatin], and 0.5 U Taq DNA polymerase (M/S Emerck Bioscience, India). DNA amplification was performed in a thermal cycler (Peqlab, Germany) programmed for a preliminary 5 min denaturation step at 94˚C, followed by 40 cycles of denaturation at 94˚C for 20 s, annealing temperature depending on the primer (50˚C - 56˚C) for 30 s and extension at 72˚C for 45 s, and finally at 72˚C for 5 min. Amplification products were separated alongside a low range molecular weight marker (M/S Emerck Bioscience, India) on a 2% (w/v) agarose gel electrophoresis in 1x TAE (Tris Acetate- EDTA) buffer stained with ethidium bromide and visualized under UV light. Gel photographs were scanned through Gel Doc System (UVITECH, UK) and the amplification product sizes were evaluated using the software Quantity one (BioRad, California, USA).
2.5. Analysis of Amplification Profile
Clearly defined ISSR bands that behave as dominant markers were scored for the presence (1) or absence (0) for all the species/cultivars and entered into a data matrix. The genetic relationships among the species/cultivar were determined by calculating the Dice coefficient, estimated as S = 2NAB/NA + NB. Where NAB is the number of amplified product common to both A and B. NA and NB corresponds to number of amplified product in A and B, respectively. Diversity patterns were represented in the form of a dendogram that was generated by subjecting the genetic similarity matrix to Unweighted Pair-group Method Arithmetic average (UPGMA) cluster analysis with software NTSYS-pc, Version 2.0  .
3. Results and Discussion
All the varieties grew and developed well. The day of emergence after sowing was ranged from 3 - 5 days and percentage of germination varied from 50% - 100%. The data relating to plant height, stem length to the first capsule, days to maturity, number of fruiting branches and number of capsules/plant were presented in Table 1. The present study also illustrates to provide additional and rapid molecular marker suitable for identification of genetic resources and also to assess phylogenetic relationship among the Sesamum varieties occurring in India. Thirty ISSR primers were used to assess the 33 genotypes of Sesamum is presented in Table 2. Out of 30 primers, 18 primers generated clear multiplex banding profiles (good/good with less bands), among which ten primers produced the best ISSR profiles. On modification of annealing temperature, MgCl2, primer concentration and template concentration did not improve their patterns. In addition, the results also showed that the most of the primers based on GA/AG and CA/GT dinucleotides core repeats generated good banding profiles. The amplification of ISSR profiles was consistent across three replicate DNA extractions from three samples, over with 98% of scorable fragments reproducible. This was in concordance with previous studies an oilseed where the reliability of the protocol was notified between PCRs and DNA extraction  . A higher concentration of MgCl2 (2.5 mM) showed good results. This may be due to non-specific amplification because of reduced enzyme activity  . The concentration of magnesium chloride affects the specificity and yield of reaction by increasing the stringency of primer annealing or has a direct effect on Taq DNA polymerase  . The selected primers showed highly polymorphic banding profile. ISSR amplification for all samples resulted in multiple band fingerprint profile for the selected ISSR primers (Figure 1). The average number of identified fragments per primer was 6.33, with a range from 4 to 8, while the average
Table 1. Agromorphological characteristics of 33 genotypes of Sesamum indicum L.
Table 2. Polymorphism detected by using the eighteen selected ISSR primers.
number of polymorphic fragments per primer was 5.38 with a range from 2 to 8. Out of the 114 scorable fragments, 97 were polymorphic revealing 85.08% polymorphism across the 33 varieties studied (Table 2). In this study primers based on GA/AG and GT/TG dinucleotides core repeats generated good profiles, which seem to indicate that the more frequent microsatellite in Sesamae contains the repeated dinucleotides (AG/GA) n and (GT/TG) n. Our results are in accordance with the SSR analysis reported earlier  . Baydar  reported that the sesame varieties with non-branching character having higher yield advantages for mechanized cultivation, because of their uniform maturity in a short growing period. The varieties were ranked with low to high seed yield as reported  . Sharma et al.  reported that RAPD and ISSR markers may be converted to co-dominant sequence characterized amplified region/sequence tagged site (SCAR/STS) markers to develop variety specific markers. The present study indicates that there was a distant variation in DNA amplification of 33 varieties of Sesamum collected from Eastern India. The matrix calculated for all possible pair wise comparisons between varieties showed that the index value varied from a minimum of
0.45 between “SI-2138-2” and “S-0140” to a maximum of 0.87 between “Usha and KIS-398” and “Kanaka” and “KIS-398” (Table 3). All the others ones have displayed
Figure 1. ISSR banding pattern of 33 varieties of Sesamum using synthesis primer: USB 843 (A), USB-835 (B), UBC-818 (C) and UBC-825 (D). M-low range DNA marker. No. 1 - 33 is the varieties of Sesamum as in Table 1.
Table 3.Similarity coefficient among the 33 varieties of Sesamum indicum obtained from ISSR markers. (V1-V33 is the varieties as indicated in Table 1).
various intermediate levels of similarity. Cluster analysis showing a dendrogram using Jaccard’s coefficient under UPGMA method divided into two major clusters at 0.45 similarity level: A1 and A2 (Figure 2). Cluster A1 had 2 varieties and cluster A2 had 31 varieties. Cluster A1 had two varieties “RJS-44” and “SI-339”. Similarly, A2 major cluster was divided into two minor groups i.e. B1 and B2 with 56% similarity with each other. The minor cluster B1 had only two varieties i.e. “S-0140” and “SI-1025”. Another minor cluster B2 was further divided into two sub-minor clusters i.e. C1 and C2 having 65% similarity with each other. The minor sub cluster C1 having four varieties (“SI- 205”, “IC-1634-3”, “S-0337” and “SI-1926”) with 70% similarity with each other... However, “SI-205” and “SI-1926” have 76% similarity among themselves. The minor sub cluster C2 was further divided into two clusters i.e. D1 and D2 with 67% similarity. In this study, it was observed that the var. Usha and Kanaka have highest similarity (87%) among themselves and similar days of maturity and day to 50% flowering. But, morphologically they were very close similarity among themselves. The varieties “Nirmala” and “Kalika” have 82% similarity with each other and 80% similarity with “Amit”. All these three varieties were morphologically alike with high germination ability and 1000 seed weight ranged from 2.46 to 2.82 g. Our results indicate that the genetic relationship among the Sesamae varieties, inferred by ISSR markers, were in accordance with their morphological characters. The principal component analysis (PCA) and 2D analysis were reflected in Figure 3 & Figure 4. The molecular markers consistently show that pollinator wasp species are deeply and consistently well resolved phylogenetically, both within and among the varieties that are associated with different sesamae germplasms as reported  . The Sesamae is a short day plant and is normally
Figure 2. UPGMA dendrogram based on the similarity coefficient, showing the clustering pattern among the 33 varieties of Sesamum. V1 - V33 is the varieties of Sesamum as in Table 1.
Figure 3. Principal Component analysis based on the similarity coefficient, showing the distribution pattern among the 33 varieties of Sesamum. V1 - V33 is the 33 varieties of Ses- amum as in Table 1.
Figure 4. 2D-analysis based on the similarity coefficient, showing the distribution pattern among the 33 varieties of Sesamum. V1 - V33 is the 33 varieties of Sesamum as in Table 1.
self pollinated, although cross pollinated ranging from 5 to over 50% occurs. Most of the varieties show an indeterminate growth habit.
This study indicates that the information for variety identification and the presence of genotypes in the distribution of parental species for Sesamae have been confirmed with both agro-mophological and ISSR markers. This analysis is quick and reproducible, and can generate sufficient polymorphism to identify the varieties, although most ISSR alleles are dominant rather than co-dominant. This investigation could be refined to include more markers and individual variety analysis for detailed characterization of Sesamae taxa that would be essential for future breeding and crop improvement program.
The authors also acknowledge to Officer-in-Charge, AICRP Sesamum, OUAT, Bhubaneswar for providing seed materials for the experiment.
 FAO (2000) The State of Food and Agriculture, 2000.
 Morris, J.B. (2002) Food, Industrial, Nutraceutical, and Pharmaceutical Uses of Sesame Genetic Resources. In: Janick, J. and Whipkey, A., Eds., Trends in New Crops and New Uses, ASHS Press, Alexandria, 153-156.
 Fukuda, Y., Nagata, M., Osawa, T. and Namiki, M. (1986) Contribution of Lignin Analogues to Antioxidative Activity of Refined Unroasted Sesame Seed Oil. Journal of the American Oil Chemists’ Society, 63, 1027-1031.
 Ram, S.G., Sundaravelpandian, K., Kumar, M., Vinod, K.K., Bapu, J.K. and Raveendran, T.S. (2006) Pollen-Pistil Interaction in the Inter-Specific Crosses of Sesamum sp. Euphytica, 152, 379-385.
 Anthony, O., Ogunshakin, R., Vaghela, S. and Patel, B. (2015) Towards Sustainable Intensification of Sesame-Based Cropping Systems Diversification in North Western India. Journal of Food Security, 3, 1-5.
 Kim, D.H., Zur, G., Danin-Poleg, Y., Lee, S., Shim, K., Kang, C. and Kashi, Y. (2002) Genetic Relationships of Sesame Germplasm Collection as Revealed by Inter-Simple Sequence Repeats. Plant Breeding, 121, 259-262.
 Bhat, K.V., Babrekar, P.P. and Lakhanpaul, S. (1999) Study of Genetic Diversity in Indian and Exotic Sesame (Sesamum indicum L.) Germplasm Using Random Amplified Polymorphic DNA (RAPD) Markers. Euphytica, 110, 21-34.
 Ercan, A.G., Taskin, M. and Turgut, K. (2004) Analysis of Genetic Diversity in Turkish Sesame (Sesamum indicum L.) Populations Using RAPD Markers. Genetic Resources and Crop Evolution, 51, 599-607.
 Dixit, A., Jin, M.H., Chung, J.W., Yu, J.W., Chung, H.K., Ma, K.H., Park, Y.J. and Cho, E.G. (2005) Development of Polymorphic Microsatellite Markers in Sesame (Sesamum indicum L.). Molecular Ecology Notes, 5, 736-738.
 Laurentin, H.E. and Karlovsky, P. (2006) Genetic Relationship and Diversity in a Sesame (Sesamum indicum L.) Germplasm Collection Using Amplified Fragment Length Polymorphism (AFLP). BMC Genetics, 7, 10-18.
 Laurentin, H.E. and Karlovsky, P. (2007) AFLP Fingerprinting of Sesame (Sesamum indicum L.) Cultivars Identification, Genetic Relationship and Comparison of AFLP Informativeness Parameters. Genetic Resources Crop Evolution, 54, 1437-1446.
 Ali, G.M., Yasumoto, S. and Seki-Katsuta, M. (2007) Assessment of Genetic Diversity in Sesame (Sesamum indicum L.) Detected by Amplified Fragment Length Polymorphism Markers. Electronic Journal of Biotechnology, 10, 12-23.
 Wu, K., Yang, M., Liu, H., Tao, Y., Mei, J. and Zhao, Y. (2014) Genetic Analysis and Molecular Characterization of Chinese Sesame (Sesamum indicum L.) Cultivars Using Insertion-Deletion and Simple Sequence Repeat (SSR) Markers. BMC Genetics, 15, 35-45.
 Zhang, P., Zhang, H.Y., Guo, W.Z., Zheng, Y.Z., Wei, L.B. and Zhang, T.Z. (2007) Genetic Diversity Analysis of Sesamum indicum L. Germplasms Using SRAP and EST-SSR Markers. Genes and Genomics, 32, 207-215.
 Pandey, S.K., Das, A., Rai1, P. and Dasgupta, T. (2015) Morphological and Genetic Diversity Assessment of Sesame (Sesamum indicum L.) Accessions Differing in Origin. Physiology and Molecular Biology of Plants, 21, 519-529.
 Banerjee, P.P. and Kole, P.C. (2009) Analysis of Genetic Architecture for Some Physiological Characters in Sesame (Sesamum indicum L.). Euphytica, 168, 11-22.
 Tabatabaei, I., Pazouki, L., Bihamta, M.R., Mansoori, S., Javaran, M.J. and Niinemets, U. (2011) Genetic Variation among Iranian Sesame (Sesamum indicum L.) Accessions vis-à-vis Exotic Genotypes on the Basis of Morpho-Physiological Traits and RAPD Markers. Australian Journal of Crop Science, 5, 1396-1407.
 Charters, Y.M., Robertson, A. and Wilkinson, M.J. (1996) PCR Analysis of Oil Seed Rape Cultivars (Brassica napus L. ssp. Ollifera) Using 5’-Anchored Simple Sequence Repeat (SSR) Primers. Theoretical and Applied Genetics, 92, 442-447.
 Hopkins, K.L. and Hilton, A.C. (2001) Optimization of Random Amplification of Polymorphic DNA Analysis for Molecular Sub Typing of Escherichia coli. Letters in Applied Microbiology, 32, 126-130.
 Saiki, R.K. (1989) The Design and Optimization of the PCR, In: Erlich, H.A., Ed., PCR Technology: Principles and Applications for DNA Amplification, Palgrave Macmillan, New York, 7-16.
 DuePham, T., ThiNguyen, T.D., Carlsson, A.S. and Buui, T.M. (2010) Morphological Evaluation of Sesame (Sesamum indicum L.) Varieties from Different Region. Australian Journal of Crop Science, 4, 498-504.
 Sharma, S.N., Kumar, V. and Mathur, S. (2009) Comparative Analysis of RAPD and ISSR Markers for Characterization of Sesame (S. indicum L.) Genotypes. Journal of Plant Biochemistry and Biotechnology, 18, 37-43.