MOLECULAR IDENTIFICATION AND GENETIC RELATIONSHIPS AMONG COFFEE SPECIES (COFFEA L) INFERRED FROM ISSR AND SRAP MARKER ANALYSES

Transcription

1 Arch. Biol. Sci., Belgrade, 63 (3), , 2011 DOI: /ABS M MOLECULAR IDENTIFICATION AND GENETIC RELATIONSHIPS AMONG COFFEE SPECIES (COFFEA L) INFERRED FROM ISSR AND SRAP MARKER ANALYSES MANOJ KUMAR MISHRA, SANDHYARANI NISHANI and JAYARAMA Central Coffee Research Institute, Coffee Research Station P.O. Chikmagalur , India Abstract The identification and genetic relationships of 23 coffee species and one coffee-related species Canthium diccocum were studied using ISSR and SRAP markers. The average polymorphism information content of SRAP primers (0.81) was lower than ISSR primers (0.86), whereas the average resolving power of the SRAP primers (9.74) is higher than the ISSR primers (8.64). The genetic similarity among the species ranged from 0.30 to 0.89 using ISSR and 0.11 to 0.90 using SRAP marker systems. Based on marker analysis, all twenty three coffee species were clustered into two major groups. Both the markers amplified species-specific fragments and are useful in genetic diversity analysis of coffee. Key words: Coffee species, genetic relationships, molecular identification, inter simple sequence repeats (ISSR), sequencerelated amplified polymorphism (SRAP). UDC :577.2 INTRODUCTION The genus Coffea L. belongs to the Rubiaceae family and comprises 103 species (Davis et al. 2006). The majority of coffee species occurs naturally in Africa, Madagascar and the Mascarenes, predominantly restricted to humid evergreen forest, but some species are found in seasonally dry deciduous forest and/r bush land (Maurin et al. 2007). The importance of coffee as an agricultural commodity relies mainly upon two varieties, Coffea arabica and C. canephora, which contribute about 65% and 35% of total production, respectively. In spite of the commercial and social importance of the genus, the genetic relationship between the majorities of coffee species is not extensively studied and their taxonomic status is poorly understood. Understanding the genetic relationship between coffee species is not only important for resolving taxonomic ambiguity but also important for the genetic improvement program. Conventionally, morphological descriptors such as growth habit, leaf type, floral characteristics and fruit morphology are used to characterize the various species. However, developing morphological descriptors for any particular species/cultivar has severe limitations as these characteristics are influenced by environmental conditions. In contrast to the morphological markers, DNA-based marker techniques are more efficient, precise and reliable for discriminating between closely related species and cultivars. In previous studies, genetic diversity and phylogenetic relationships between different coffee species was carried out by using both Random Amplified Polymorphic DNA (RAPD) analysis as well as Restriction Fragment Length Polymorphism of chloroplast and mitochondrial DNA (Berthou, 1983; Cross et al. 1998; Lashermes et al. 1993; Lashermes et al. 1996; Orozco-Castillo, et al. 1996). Recently, Davis et al. (2006) published a detailed annotated taxonomic conspectus of the genus Coffea and included five indigenous coffee species from India, of which C. bengalensis, C. travancorensis, C. wightiana are placed under the genus Psilanthus and an- 667

2 668 MANOJ KUMAR MISHRA ET AL. other two species, C. khasiana and C. jenkinsii, are placed under the genus Nostolachma, both under Coffeeae. Based on the analysis of the internal transcribed spacers (ITS) of nuclear DNA, Lashermes et al. (1997) observed limited sequence divergence between Psilanthus and Coffea and concluded that Psilanthus should not be recognized as a separate genus from Coffea. Based on the comparative sequence analysis of plastid DNA, Cross et al. (1998) concurred with Lashermes et al. (1997) about the close relationship between Coffea and Psilanthus although their tree topology shows an unresolved relationship between two species of Psilanthus (P. mannii and P. ebracteolatus) studied by them and Coffea. However, both these studies did not include a representative of closely related genera for comparison. Further studies by Andreasen et al. (1999), Andreasen and Bremer (2000) and Davies et al. (2007) showed that Psilanthus is nested within Coffea and stressed the importance of species level studies to resolve the phylogenetic relationship. Maurin et al. (2007) observed that the separation of Coffea and Psilanthus based on morphological features is not convincing and suggested the inclusion of molecular data from other species of Psilanthus to establish the relationship between these two genera. Based on the morphological features, Nostolachma is placed under Coffeeae (Davis et al. 2007) although molecular data from additional species may prove to be useful. Although various researchers have studied the taxonomic and phylogenetic relationship between various coffee species, there are no reports yet available for the species-specific markers that are useful for distinguishing various species. In this paper, we have used two PCRbased molecular markers, Inter Simple Sequence Repeats (ISSR) and Sequence Related Amplified Polymorphism (SRAP), in developing species specific markers for 23 coffee species available in coffee germplasm in India. To our knowledge, both these markers are used in coffee for developing species diagnostic markers for the first time. The main objectives of the study are to generate a molecular database and identification tag for each species that may be useful for conservation and systematic studies as well as coffee breeding programs. MATERIALS AND METHODS Plant materials Leaf material was collected from 10 individuals of twenty three coffee species (Table 1) including five indigenous coffee species from India (C. bengalensis, C. travancorensis, C. wightiana, C. khasiana and C. jenkinsii), maintained at the germplasm plot of Central Coffee Research Institute, Chikmagalur, Karnataka and Regional Coffee Research Station Thandigudi, Tamil Nadu. Mature seeds were also collected and grown in earthen pots under net house conditions. The leaf materials were used for DNA isolation. DNA extraction Genomic DNA was extracted from fresh young leaves using a modified CTAB method. About 200 mg of leaf tissue was ground to a fine powder in liquid nitrogen, it was transferred to a 30 ml tube containing 5 ml preheated extraction buffer (2% CTAB (w/v), 100mM Tris-HCL (ph 8.0), 25mM EDTA, 2M Nacl, 0.1 % beta-mercaptoethanol). The tubes were incubated at 60 0 C for 1 h with occasional shaking. After incubation, the tubes were cooled to room temperature and centrifuged at 6000 rpm for 20 min. The supernatant was transferred to a new tube and extracted twice with equal volumes of chloroform-isoamyl alcohol (24:1). The supernatant was transferred to 2 ml tubes, precipitated with 0.7 vol of isopropanol at room temperature for 30 min., and then centrifuged at 8000 rpm for 20 min at 4 C. The pellet formed after centrifugation was washed with 75% (v/v) ethanol for 10 min and dissolved in 60 µl of Tris-EDTA (1-10mM). The concentration of DNA was measured using 0.8% agarose gel stained with ethidium bromide as well as by UV spectrophotometry at 260 nm and a 280 nm. The resuspended DNA was then diluted in sterile distilled water to 10ng/ µl concentration for use in amplification reactions. ISSR, SRAP analysis A total of 20 ISSR primers from the University of British Columbia and 60 SRAP primer combina-

3 Molecular based diversity analysis of coffee species 669 tions (8 forward primers and 14 reverse primers) synthesized by Sigma, India, were initially screened to determine the suitability of each primer for the study. Primers were selected for further analysis based on their ability to detect clear and distinct polymorphic amplification products within the species of Coffea. 12 ISSR primers and 21 SRAP primer combinations that had a high level of polymorphism and the best readability were used for PCR amplification (Table 2). The PCR reaction was carried out in a palm cycler (Corbett Research). ISSR PCR was conducted in a 20 µl reaction mixture containing 1x reaction buffer (10 mm Tris- HCL ph 8.8, 50 mm KCL 0.08% Nonidet P40), 30 ng DNA, 200 µm dntp mixture, 2.5 mm MgCl 2, 3 µm ISSR primer,0.2 µl of formamide and 1.0 U Taq DNA polymerase. The ISSR amplification conditions were: 3 min initial denaturation at 95ºC; 30 cycles consisting of 1 min denaturation at 94ºC; 1.30 min primer annealing at 55ºC and 2 min extension at 72ºC and a final 10 min extension at 72 º C. SRAP analysis was performed by adapting the procedure described by Li and Quiros (2001) with minor modifications: 20 µl reaction mixture containing 1 x reaction buffer (75mM Tris-HCL ph 8.8, 20 mm (NH 4 ) 2 SO 4, 0.01% Tween 20), 30 ng template DNA, 200 µm dntp mixture, 2.5 mm MgCl 2, 3 µm each of forward and reverse primer, and 1.0 U Taq DNA polymerase. The SRAP amplification program was a 4 min initial denaturation at 96º C; 5 cycles consisting of 1 min denaturation at 94ºC, 1.15 min primer annealing at 35ºC; and 2 min extension at 72ºC followed by 30 cycles consisting of 1 min denaturation at 94ºC, 1.15 min primer annealing at 50ºC; and 2 min extension at 72 ºC; and a final extension of 15 min at 72ºC. The PCR products of both ISSR and SRAP were run on 2% (w/w) agarose gels containing 0.5 µg ethidium bromide/ml in 1X TAE buffer and then visualized and photographed using the UV-transilluminator (SYNGENE) and documented using the Gene Snap software program. Both ISSR and SRAP PCR reactions were repeated at least twice to confirm the reproducibility of each PCR band. Data analysis The ISSR and SRAP amplified bands were scored for the presence (1) or absence (0). The total number of bands, the distribution of bands across all species, polymorphic bands, species-specific bands and average number bands per primer were calculated. The value of each primer was assessed using two indices; PIC, which is the same as a diversity index (Botstein et al. 1980; Milbourne et al. 1997) and resolving power (Rp) (Prevost and Wilkinson, 1999). PIC or DI was estimated as PIC= (1-p 2 i )/n, where n is the number of band positions analyzed in all the species, p i is the frequency of the banding pattern. The resolving power of a primer is Rp = I b where I b (band informativeness) takes the value of 1- [2x (0.5-p)] and p is the ratio of six species sharing the band. A pairwise similarity matrix was constructed using the Dice similarity coefficient (Sneath and Sokal, 1973). The relationship between the species was displayed as a dendrogram constructed using NTSYS PC 2.10e software (Rohlf, 1995) based on Unweighted Pair Group Method using Arithmetic averages (UPG- MA). Statistical support of the clusters was assessed by means of 1000-bootstrap replicates. RESULTS ISSR polymorphism and species identification Out of the 20 ISSR primers initially screened, 12 primers were found to be polymorphic and produced clear and reproducible amplification patterns. These 12 primers could produce 185 distinct reproducible bands across the 24 species with an average of per primer (Table 3). The size of the amplified products ranged from 100 to 2900 bp. Of the total 185 amplified bands, 174 (93.06%) were polymorphic, with an average of 14.5 polymorphic fragments per primer. All the ISSR primers except UBC-881 and UBC-855 showed 100% polymorphism. The primers UBC-881 and UBC-855 showed 8.3% and 90.90% of polymorphism respectively. The average poly-

4 670 MANOJ KUMAR MISHRA ET AL. Fig.1. The extent of polymorphism observed among various coffee species listed in table 1 by ISSR primer 842 Lanes M: Molecular ruler (Gene ruler 100 bp ladder plus). morphism percentage recorded with ISSR primers is The resolving power of the 12 ISSR primers tested ranged from 4.08 (UBC- 840) to (UBC-881) with an average of Similarly, the polymorphism information content (PIC) of the 12 ISSR primers ranged from (UBC- 881) to (UBC-840) with average of Fig. 2 Dendrogram generated using unweighted pair group method with arithmetic average (UPGMA) analysis showing relationships among different coffee species using ISSR data. Out of the 12 ISSR primers, only one ISSR polymorphic primer (UBC 842) could discriminate all the species independently (Fig. 1). Twelve ISSR primers generated 40 unique fragments of which 8, 7 and 5 fragments are obtained in C. khasiana, C. jenkinsii and C. arabica, respectively. Further, the ISSR primer UBC-836 has amplified species diagnostic fragments in five different species (Table 4). The five Indian coffee species generated 19 unique fragments of which maximum numbers of eight unique fragments were generated by C. khasiana, whereas only one unique fragment was generated by both by C. travancorensis and C. wightiana. The genetic similarity derived from the data of the ISSR marker analysis varied from 0.27 between C. bengalensis and C. liberica to 0.89 between C. canephora and C. congensis (Table 5). Among the Indian species, C. travancorensis showed maximum genetic similarity with C. wightiana, based on the ISSR marker analysis. Furthermore, Canthium diccocum, which is used as a reference species, showed maximum genetic similarity (0.65) with C. racemosa. The dendrogram based on ISSR data was constructed by UPGMA analysis, grouping all of the coffee species into two major clusters (Fig. 2). The

5 Molecular based diversity analysis of coffee species 671 Fig.3 The extent of polymorphism observed among various coffee species listed in Table 1 by SRAP primer combinations Me3 and Em12 Lanes M: Molecular ruler (Gene ruler 100 bp ladder plus). first major cluster divided into two minor clusters of which the first minor cluster again divided into two sub-minor clusters. The first sub-minor cluster consisted of C. congensis, four different varieties of C. canephora, C. arabica and C. liberica. The second sub-minor cluster was represented by C. eugenioides, C. salvatrix, C. kapakata, C. arnoldiana, and C. dewevrei var. excelsa. The second minor cluster included only the Indian species C. jenkinsii. The second major cluster divided in to two minor clusters of which the Fig.4 Dendrogram generated using the unweighted pair group method with arithmetic average (UPGMA) analysis showing relationships among different coffee species using SRAP data. Fig.5 Dendrogram generated using the unweighted pair group method with arithmetic average (UPGMA) analysis showing relationships among different coffee species using ISSR and SRAP data.

6 672 MANOJ KUMAR MISHRA ET AL. Table 1. Coffee species used in the study with their code, origin and conservation status Coffee Species Species code Place of origin/distribution Conservation status Coffea congensis A. Froehner C1 West Central Africa Congo Least concern C. canephora Pierre ex A. Froehner C2 West Tropical Africa Least concern C. canephora var. laurentii (De Wild.) A.Chev. C3 West Tropical Africa Least concern C. canephora var. ugandae (Cramer) A.Chev. C4 West Tropical Africa Least concern C. canephora var. quillon Philipp. C5 West Tropical Africa Least concern C. arabica cv. Kents L. C6 Ethiopia (Indian cultivar) Vulnerable C. eugenioides Moore C7 West Central Africa Least concern C. zanguebariae Lour. C8 East Tropical Africa Vulnerable C. racemosa Lour. C9 Southern Tropical Africa Near Threatened C. salvatrix Swynn.&Philipson C10 East Tropical Africa Near Threatened C. kapakata (A.Chev.) Bridson C11 West Angola Vulnerable C. stenophylla G.Don C12 West Tropical Africa Least concern C. abeokutae Cramer Ex De Wild. C13 West Tropical Africa Least concern C. liberica Bull.ex Hiern C14 West Tropical Africa Least concern C. dewevrei De Wild. & T.Durand C15 Democratic Republic of Congo Vulnerable C. arnoldiana De Wild. C16 West Central Africa Least concern C. dewevrei var. aruwimiensis (De Wild.) A.Chev. C17 West Central Africa Least concern C. dewevrei var. Excelsa (A.Chev.) C18 West Central Africa Least concern C. bengalensis Roxb.ex Schult.) C19 India endangered C. travancorensis Wight &Arn C20 India endangered C. khasiana (Korth.) Hook.f. C21 India endangered C. jenkinsii Hook.f. C22 India endangered C. wightiana Wall. Ex Wight &Arn C23 India endangered Canthium diccocum (Gaertn) C24 South Asia Threatened Table 2. Primer sequences used for SRAP analysis of coffee species Forward primer (5 3 ) Reverse primer (5 3 ) Me1TGAGTCCAAACCGGATA Em2 GACTGCGTACGAATTTGC Me2 TGAGTCCAAACCGGAGC Em3 GACTGCGTACGAATTGAC Me3 TGAGTCCAAACCGGAAT Me4 TGAGTCCAAACCGGACC Me6 TGAGTCCAAACCGGACA Me9 TGAGTCCAAACCGGAGG Me10 TGAGTCCAAACCGGAAA Me11 TGAGTCCAAACCGGAAC Em4 GACTGCGTACGAATTTGA Em5 GACTGCGTACGAATTAAC Em6 GACTGCGTACGAATTGCA Em7 GACTGCGTACGAATTCAA Em9 GACTGCGTACGAATTCAG Em10 GACTGCGTACGAATTCAT Em11 GACTGCGTACGAATTCTA Em12 GACTGCGTACGAATTCTC Em13 GACTGCGTACGAATTCTG Em14 GACTGCGTACGAATTCTT Em15 GACTGCGTACGAATTGAT Em16 GACTGCGTACGAATTGTC first minor cluster is again divided in to two sub-minor clusters. The first sub-minor cluster consisted of C. zanguebariae, C. racemosa, C. stenophylla, C. abeokutae, C. dewevrei and C. dewevrei var. aruwimiensis including C. diccocum and the second sub-minor cluster included C. bengalensis, C. travancorensis and C. wightiana whereas C. khasiana is represented in the second minor cluster independently.

7 Molecular based diversity analysis of coffee species 673 Table 3. Polymorphism obtained by RAPD, ISSR and SRAP analysis in coffee species Primer/ primer Total Number of bands in No. of Polymorphic bands morphism Percentage of poly- Size range (bp) combinations bands each species RP PIC ISSR UBC (5.25) UBC (4.0) UBC (3.20) UBC (3.58) UBC (3.5) UBC (3.21) UBC (2.04) UBC (3.63) UBC (4.08) UBC (4.96) UBC (3.33) UBC (11.08) Total (51.86) Average (4.32) SRAP Me1-Em (6.16) Me1-Em (5.958) Me1-Em (4.13) Me1-Em (3.16) Me1-Em (4.04) Me2-Em (5.12) Me2-Em (5.70) Me2-Em (6.08) Me2-Em (3.04) Me2-Em (3.45) Me3-Em (6.91) Me3-Em (5.70) Me3-Em (3.75) Me3-Em (3.33) Me3-Em (4.0) Me4-Em (1.91) Me4-Em (6.83) Me6-Em (5.79) Me6-Em (6.46) Me6-Em (6.16) Me10-Em (5.0) Total (102.67) Average (4.88) SRAP polymorphism and species identification Sixty SRAP primer combinations were initially screened against 23 coffee species along with the reference species, C. decorum, of which 21 primer combinations are found to be highly polymorphic and produce clear amplification patterns. These 21 primers could produce 336 distinct scorable bands with

8 674 MANOJ KUMAR MISHRA ET AL. Тable 4. Species- diagnostic ISSR and SRAP markers in coffee species Species ISSR marker SRAP marker C. congensis --- Me6-Em3-700, Me6-Em8-550, Me3-Em C. canephora --- Me6-Em5-1030, Me6-Em5-75 C. canephora var. laurentii UBC 855 Me1-Em4-250, Me3-Em3-375, Me2-Em3-500, 250 C. canephora var. ugandae Me6-Em C. canephora var. quillon UBC-855 Me4-Em1 1350, Me6-Em3-1500, Me6-Em5-800 C. arabica cv. Kents UBC-841 UBC-826 UBC-811 UBC-880 Me1-Em7-700, Me4-Em1 600, Me6-Em5-1500, Me10- Em , Me3-Em11-550, Me3-Em C. eugenioides UBC-841 Me1-Em4-150, Me1-Em6-800, Me2-Em3-350, Me2-Em4-375, Me2-Em6-375, Me3-Em7 500, 300, Me4-Em1 850, C. zanguebariae UBC-836 C. racemosa Me3-Em7 600, C. salvatrix Me1-Em4-200, Me2-Em6-550, Me3-Em4-1400, 520, Me4-Em2 UBC , Me6-Em3-500, Me6-Em5-150, Me1-Em4 240, Me2-Em12 UBC , 550, 400, Me2-Em10 550, Me3-Em ,1100,900 C. kapakata Me4-Em2 1030, Me6-Em C. stenophylla Me6-Em5-300, 220, Me10-Em13-800, Me3-Em C. abeokutae Me3-Em C. liberica UBC-811 UBC-880 Me1-Em6-900, Me1-Em7-650, Me6-Em5-200, Me2-Em10 220, C. dewevrei Me2-Em C. arnoldiana Me2-Em6-900, Me6-Em5-1300, Me6-Em8-700, Me2-Em , C. dewevrei var. aruwimiensis C. dewevrei var. excelsa UBC-836 Me3-Em4-1600, Me1-Em12 680, Me10-Em13-330, Me2-Em C. bengalensis UBC-811 Me1-Em6-300, Me2-Em12 400, Me10-Em , 750, 350, UBC-880 Me2-Em10 750, Me3-Em C. travancorensis Me4-Em2 300, Me6-Em5-600 C. khasiana C. jenkinsii UBC-836 UBC-842 UBC-855 UBC-880 UBC-836 UBC-840 UBC-834 UBC-855 UBC-810 Me2-Em3-100, Me3-Em7 350, Me6-Em8-500, Me2-Em12 175, Me10-Em13-600, Me3-Em Me1-Em2-100, Me1-Em7-450, Me3-Em4-75, Me4-Em2 400, Me6-Em3-820, Me6-Em5-250, Me1-Em12 720, 400,, Me2- Em12 120, Me10-Em13-450, 170, Me2-Em10 600, Me3-Em11-600, 500, Me3-Em , 650 C. wightiana UBC-836 Me1-Em4-550, Me2-Em12 480, Me10-Em13-270, Me3-Em Canthium diccocum UBC-855 Me3-Em3-1450, Me4-Em1 700, Me3-Em the number of amplified fragments ranging from 10 (Me1-Em7) to 23 (Me1-Em12) with an average of 16 per primer combination (Table 3). The size of the amplified products ranged from 50 to 3500 bp in different primer combinations. Of the total 336 amplified bands, 323 were polymorphic, with an average polymorphic fragments per primer combination. The percent of polymorphism ranged from 83.33% to a maximum of 100% with an average of 96.12% polymorphism. Of the total 21 primer combinations,

9 Molecular based diversity analysis of coffee species 675 Table. 5 Similarity co-efficient of coffee species using ISSR marker (DICE) C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C C C C C C C C C C C C C C C C C C C C C C C C primer combinations showed 100% polymorphism. The resolving power (RP) of 21 SRAP primer combinations ranged from 3.83 (Me4-Em1) to (Me3-Em3) with an average of Similarly, the polymorphism information content (PIC) or the genetic diversity of 21 SRAP primer combinations ranged from 0.61 (Me1-Em2) to 0.96 (Me1-Em12, Me3-Em11) with an average of Only three SRAP primer pairs could discriminate all the species independently, although none of these primers could amplify species-specific fragments for all the 23 species (Fig. 3). 20 SRAP primer combinations generated 75 unique fragments in 23 coffee species out of which 5 Indian species generated 30 unique fragments. Among the Indian species, the maximum number of unique fragments was generated by C. jenkinsii (16) followed by C. khasiana (8) and the least number of unique fragments was in C. travancorensis (1) which could be used as a unique fingerprinting tools. The genetic similarity derived from the data of the SRAP marker analysis varied from 0.11 between C. wightiana and C. congensis to 0.90 between C. canephora and C. canephora var. laurentii (Table 6). Among the Indian species, maximum similarity was observed between C. travancorensis and C. wightiana and the least genetic similarity was obtained between C.

10 676 MANOJ KUMAR MISHRA ET AL. Table.6 Similarity co-efficient of Indian coffee species using SRAP marker (DICE) C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C C C C C C C C C C C C C C C C C C C C C C C C khasiana and C. jenkinsii based on the SRAP marker analysis. The UPGMA clustering algorithm from the SRAP analysis grouped the 23 coffee species and C. diccocum into two major clusters (Fig. 4). The first major cluster was represented by C. wightiana and the second major was divided into two minor clusters of which the first minor cluster consisted of C. jenkinsii alone. The second minor cluster divided in to two sub-minor clusters. The first sub-minor cluster included C. congensis, C. canephora varieties, C. arabica and C. liberica, whereas the second sub-minor cluster further dived into two groups. The first group included C. zanguebariae, C. racemosa, C. abeokutae, C. dewevrei, C. dewevrei var. aruwimiensis, C. stenophylla, C. arnoldiana and C. dewevrei var. excelsa, whereas the second group included three Indian species viz. C. bengalensis, C. travancorensis, C. khasiana along with Canthium diccocum. A strict consensus tree based on both ISSR and SRAP data were constructed (Fig. 5). This dendrogram shows maximum similarity to ISSRbased dendrogram except that the position of C. arnoldiana and C. dewevrei var. excelsa was taken over by C. zanguebariae and C. racemosa, respectively, in the first major sub-cluster. The genetic similarities obtained using both ISSR and SRAP data has revealed the highest similarities (0.90) between C. canephora and C. canephora var. laurentii and the lowest similarities (0.36) between C. liberica and C. bengalensis (Table 7).

11 Molecular based diversity analysis of coffee species 677 Table. 7 Similarity co-efficient of coffee species using combined ISSR and SRAP markers (DICE) C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 C22 C23 C24 C C C C C C C C C C C C C C C C C C C C C C C C DISCUSSION Analysis of crop genetic diversity is very important for breeding and conservation programs, and molecular markers offer an approach to unveil the genetic diversity among different species and cultivars based on nucleic acid polymorphisms. In this study, both ISSR and SRAP marker systems were simultaneously used to investigate the genetic diversity among 23 coffee species available in the germplasm collection of India. The results showed that both ISSR and SRAP markers were suitable for genetic diversity analysis in coffee by amplifying several species specific diagnostic markers. Species specific ISSR primers were generated in eucalyptus (Balasaravanan et al. 2006) and oak (Carvalho et al. 2009). In the present study, species-diagnostic ISSR markers were identified in 12 coffee species including all the indigenous coffee species from India except C. travancorensis (Table 4). Similar problems were experienced while developing species-specific ISSR markers in eucalyptus (Balasaravanan et al. 2006). This could be due to the occurrence of a high divergence among the particular species population and/or the low number of ISSR primers used in this study. In contrast to the ISSR markers, SRAP is more effective and amplified species specific markers in all the coffee species except C. zanguebariae, and C. dewevrei var. aruwimiensis (Table. 4). Both genetic factors and selection pressure influence genetic diversity (Sun and Lin 2003). Since various coffee species have their origin in different agro-climatic zones in their native habitat, with varying selection pressure during the course of evolution, it is therefore not surprising to find wide polymorphism among different coffee species. The percentage of polymorphic bands detected by SRAP primer

12 678 MANOJ KUMAR MISHRA ET AL. combinations (96.12%) was higher compared to ISSR (93.06%). Similarly, the number of polymorphic bands detected by SRAP primer combinations (15.38) was higher than that obtained by the ISSR primer (14.5). This indicated that SRAP markers are more efficient than ISSR markers in genetic diversity analysis of coffee. The average PIC of the SRAP primers (0.82) was less than that of the ISSR primers (0.84), whereas the average RP of the SRAP primers (9.74) was higher than that of the ISSR primers (8.64). The difference in PIC and RP values among ISSR and SRAP markers is expected because both markers differ in their operational principle and each marker targets a different region of the genome. This is evident from the UPGMA clustering of the coffee species using ISSR and SRAP markers. While both the C. travancorensis and C. wightiana from the southern peninsular India showed maximum similarity between themselves (Table 5) and placed close to each other in the dendrogram (Fig. 2) using ISSR markers, using a SRAP marker, C. travancorensis shared maximum similarities with C. bengalensis (Table 6) and placed closely in the dendrogram (Fig. 4). Interestingly, the genetic similarities observed between various indigenous coffee species with other coffee species using an ISSR marker varied from 0.3 to 0.6 and 0.1 to 0.6 using SRAP markers (Tables 5 and 6). This difference could be explained by the fact that ISSR markers target the region within the microsatellite repeats whereas SRAP markers preferentially detect polymorphism in coding sequences which are usually conserved among closely related cultivars and species with low mutation rate. Therefore, simultaneous use of different types of molecular markers may be useful in generating all- sided information. In conclusion, this study provides information on genetic relatedness not only among different coffee species but also on their relatedness with the indigenous coffee species which are either placed under the genus Psilanthus or Nostolachma. The present study has clearly demonstrated a close relationship between Indian indigenous coffee species and other wild coffee species and supports the contention of their inclusion in the genus Coffea rather than placing them under the genus Psilanthus or Nostolachma, as suggested by previous workers (Cross et al., 1998; Lashermes et al. 1997; Andreasen et al. 1999, Andreasen and Bremer, 2003 and Davies et al. 2007). Further understanding of the level and partitioning of genetic variation within the species would provide an important input in designing appropriate breeding exercises and conservation strategies. REFERENCES Andreasen, K., and B. Bremer (2000). Combined phylogenetic analysis in the Rubiaceae-Ixoroidae: morphology, nuclear and chloroplast DNA data. American Journal of Botany 87, Andreasen, K., Baldwin, B.G., and B. Bremer (1999). Phylogenetic utility of the nuclear rdna ITS region in the subfamily Ixoroideae (Rubiaceae): comparisons with cpdna and rbcl sequence data. Plant Systematics and Evolution 217, Balasaravanan, T., Chezhian, P., Kamalakannan, R., Yasodha, R., Varghese, M., Gurumurthi, K., and M. Ghosh (2006). Identification of species-diagnostic ISSR markers for six Eucalyptus species. Silvae Genetica 55(3), Berthou, F., Mathieu, C., and F. Vedel (1983). Chloroplast and mitochondrial DNA variation as an indicator of phylogenetic relationships in the genus Coffea L. Theor Appl Genet 65, Botstein, D., White, R.L., Skolnick, M.H., and R.M. Davies (1980). Construction of a genetic map in man using restricted length polymorphism. American Journal of Human Genetics 32, Carvalho, A., Paula, A., Guedes-Pinto, H., Martins, L., Carvalho, J., and J. Lima-Britto (2009). Preliminary genetic approach based on both cytogenetic and molecular characterisations of nine oak species. Plant biosystem. 143, S Cross, J., Combes M.C., Trouslet, P., Anthony, F., Hamon, S., Charrier, A., and P, Lashermes. (1998). Phylogenetic analysis of chloroplast DNA variation in Coffea L. Molecular Phylogenetics and Evolution 9, Cros, J., Combes, M.C., Trouslot, P., Anthony, F., Hamon, S., Charrier, A., and P. Lashermes (1998). Phylogenetic analysis of chloroplast DNA variation in Coffea L. Molecular Phylogenetics and Evolution 9, Davis, A.P., Chester, M., Maurin, O., and M.F. Fay (2007) Searching for the relatives of coffea (Rubiaceae, Ixoroideae): the circumscription and phylogeny of Coffeeae based on plas-

13 Molecular based diversity analysis of coffee species 679 tid sequence data and morphology. American Journal of Botany 94 (3), Davis, A.P., Governs, R., Birdsong, D.M., and P. Stoffelen (2006). An annotated taxonomic conspectus of the genus Coffee L. (Rubiaceae). Botanical Journal of Linnaean Society 152, Lashermes, P., Cross J., Marmey, P., and A. Charrier (1993). Use of random amplified DNA markers to analyze genetic variability and relationships of Coffea species. Genet Res Crop Evol 40, Lashermes, P., Cross J., Combes, M.C., Trouslot, P., Anthony, F., Hamon, S., and A. Charrier (1996). Inheritance and restriction fragment length polymorphism of chloroplast DNA in the genus Coffea L. Theor Appl Genet 93, Lashermes, P., Combes, M.C., Trouslot, P., and A. Charrier (1997). Phylogenetic relationship of coffee tree-species (Coffea L.) as inferred from ITS sequences of nuclear ribosomal DNA. Theoretical and Applied Genetics 94, Li, G., and C. F. Quiros (2001). Sequence related amplified polymorphism (SRAP), a new marker system on a simple PCR reaction: its application to mapping and gene tagging in Brassica. Theoretical and Applied Genetics 107(1), Maurin, O., Davis, A. P., Chester, M., Mvungi, E.F., Jaufeerally- Fakim, Y., and M. F. Fay (2007). Towards a phylogeny for Coffea (Rubiaceae): Identifying Well supported Lineages based on Nuclear and Plastid DNA sequences. Annals of Botany 100, Milbourne, D., Meyer, R., Bradshaw, J.E., Baird, E., Bonar, N., Provan, J, Powell, W., and R. Waugh (1997). Comparison of PCR based marker systems for the analysis of genetic relationships in cultivated potato. Molecular breeding 3 (2), Orozco-Castillo, C., Chalmers, K. J., Powell, W., and R. Waugh (1996). RAPD and organelle specific PCR re-affirms taxonomic relationships within the genus Coffea. Plant Cell Reports 15, Prevost, A., and Wilkinson, M. J. (1999). A new system of comparing PCR primers applied to ISSR finger printing of potato cultivars. Theoretical and Applied Genetics 98, Rohlf, F. J. (1995). NTSYS-pc Numerical taxonomy and multivariate analysis system version 1.80, Exterior software, Setauket, New York. Sneath, P. H. A., and R. R. Sokal (1973). Numerical taxonomy, The principles and practice of numerical classification; W.H. Freeman & Co, San Francisco, California, P.573. Sun, Y., and F. C. Lin (2003). Analysis of genetic diversity in natural germplasm of Lentinula edode in China using RAPD technique (Chinese), Mycosystema 22 (3),