About this Journal Submit a Manuscript Table of Contents
The Scientific World Journal
Volume 2012 (2012), Article ID 939820, 11 pages
http://dx.doi.org/10.1100/2012/939820
Research Article

Genetic Diversity of Arabica Coffee (Coffea arabica L.) in Nicaragua as Estimated by Simple Sequence Repeat Markers

1Department of Plant Breeding and Biotechnology, Swedish University of Agricultural Sciences, P.O. Box 101, 230 53 Alnarp, Sweden
2Department of Plant Protection, National Agrarian University, P.O. Box 453, Managua, Nicaragua

Received 12 February 2012; Accepted 28 February 2012

Academic Editors: A. Belyayev and G. Kunze

Copyright © 2012 Mulatu Geleta et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Coffea arabica L. (arabica coffee), the only tetraploid species in the genus Coffea, represents the majority of the world’s coffee production and has a significant contribution to Nicaragua’s economy. The present paper was conducted to determine the genetic diversity of arabica coffee in Nicaragua for its conservation and breeding values. Twenty-six populations that represent eight varieties in Nicaragua were investigated using simple sequence repeat (SSR) markers. A total of 24 alleles were obtained from the 12 loci investigated across 260 individual plants. The total Nei’s gene diversity ( 𝐻 𝑇 ) and the within-population gene diversity ( 𝐻 𝑆 ) were 0.35 and 0.29, respectively, which is comparable with that previously reported from other countries and regions. Among the varieties, the highest diversity was recorded in the variety Catimor. Analysis of variance (AMOVA) revealed that about 87% of the total genetic variation was found within populations and the remaining 13% differentiate the populations ( 𝐹 𝑆 𝑇 = 0 . 1 3 ; 𝑃 < 0 . 0 0 1 ). The variation among the varieties was also significant. The genetic variation in Nicaraguan coffee is significant enough to be used in the breeding programs, and most of this variation can be conserved through ex situ conservation of a low number of populations from each variety.

1. Introduction

Coffea arabica L. (arabica coffee) is a self-fertile allotetraploid species that belongs to the genus Coffea in the family Rubiaceae [1, 2]. Out of the 103 species in the genus, arabica coffee is the only tetraploid species ( 2 𝑛 = 4 𝑥 = 4 4 ), the remaining species being diploid with 2 𝑛 = 2 𝑥 = 2 2 chromosomes [3]. Arabica coffee originated from a relatively recent hybridization between Coffea canephora (robusta coffee) and C. eugenioides or their ecotypes in the plateaus of Central Ethiopia [2, 4]. Coffee is mainly grown in tropical and subtropical regions and is an important cash crop in more than 60 countries in South and Central America, Asia, and Africa with an acreage of over 11 million ha [5].

Coffee production is an important economic activity in Central America and accounts for about 10% of the world coffee production [6]. In Nicaragua, large-scale coffee production was started in the 1850s, and since 1870, coffee is the main export crop [7, 8]. Most of the coffee production in the country comes from arabica coffee and the most cultivated varieties are Caturra, Catuai, Bourbon, and Typica [9]. More than 70% of Nicaraguan coffee is produced at elevations between 600 and 1500 m above sea level (asl) in the north central part of the country, where it is considered optimal for coffee production; and the rest is produced below 600 m asl in the south pacific region [10]. During 2010, the total production in the country was about 78 kilo tonnes [11].

Several studies have shown that the genetic diversity of arabica coffee is low when compared to that of robusta coffee [2, 1219] due to its narrow genetic base associated with autogamy, evolutionary history, and domestication. This narrow genetic base has been reflected in different forms that include the lack of resistant genotypes to various pests and diseases (e.g., [2022]). The genetic base of arabica coffee in the American content is even narrower, as it represents only a small subset of the genetic variations present within the arabica coffee gene pool [14, 20] and are more prone to various pests and diseases [20, 23, 24]. Thus, enhancement of its resistance to pests and diseases is becoming a crucial priority for economic and sustainable coffee production. This is being done through crossing arabica coffee with other coffee species, particularly robusta coffee [12, 21, 23, 25, 26] and through selection of genotypes of interest from the arabica coffee gene pool [27, 28].

Detecting and quantifying genetic variation in crop species is important for successful conservation of genetic resources and plant breeding. Molecular marker techniques, such as random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), and simple sequence repeats (SSR, also known as microsatellites) have been used for genetic diversity analysis in wild and cultivated coffee [13, 14, 1619, 2931]. However, there is little information on the genetic diversity of arabica coffee varieties in Nicaragua. Hence, the present study was conducted to estimate the genetic diversity and population genetic structure of arabica coffee in Nicaragua using SSR markers.

2. Materials and Methods

2.1. Plant Material and DNA Extraction

Coffee seeds from twenty-six populations representing eight arabica coffee varieties were used in this study. Each population was represented by ten individual coffee trees. Fresh coffee berries were collected between December 2009 and February 2010 from the main coffee growing provinces of Nicaragua (Table 1). The berries were dried up at room temperature and processed to obtain seeds. The seeds were then grown in pots in a greenhouse at a mean temperature of 28°C. Individually sampled leaf tissue from the plants grown in the greenhouse was placed in 2 mL Eppendorf microcentrifuge tubes and immediately frozen in liquid nitrogen and stored at −80°C until DNA extraction. After the frozen samples were milled using a Retsch MM400 shaker (Haan, Germany), DNA was extracted using a modified CTAB procedure, as described in Bekele et al. [32]. DNA quality and concentration was measured using a Nanodrop ND-1000 spectrophotometer (Saveen Werner, Sweden).

tab1
Table 1: List of 26 arabica coffee populations used in this study representing eight coffee varieties in Nicaragua, and their geographic descriptions and SSR-based estimates of gene diversity ( 𝐻 L o c i ) and percent polymorphic loci (PPL).
2.2. SSR-PCR

Twenty-five SSR primer-pairs were initially screened for good amplification, polymorphism, specificity to their target loci, and suitability of the allele size for multiplexing. This led to the selection of twelve primer-pairs for final analysis (Table 2). The forward primers of selected primer-pairs were fluorescently 5′-labeled with either 6FAM, VIC, NED, or PET fluorescent dyes. The reverse primers were PIG-tailed with “GCTTCT” to avoid a nontemplated addition of a single nucleotide by Taq DNA polymerase to the PCR product, as described in Ballard et al. [34].

tab2
Table 2: List of primer-pairs used to amplify the SSR loci used in this study and their annealing temperature: the repeat motifs of the loci and observed fragment sizes of the alleles.

The PCR reactions were carried out in a volume of 25 μL containing 25 ng genomic DNA, 0.3 μM forward and reverse primers, 2 mM MgCl2, 0.3 mM dNTPs, 1 U Taq DNA polymerase (Sigma, Germany), and 1 × PCR buffer (10 mM Tris-HCl, pH 8.3 and 50 mM KCl). The reactions were performed using the GeneAMP PCR system 9700 thermocycler using the following temperature profiles: initial denaturation at 95°C for 3 min, followed by six touchdown cycles of denaturation at 94°C for 30 sec, annealing at 𝑋 - 𝑌 °C (−1°C/cycle) for 30 sec and extension at 72°C for 45 sec, and then 32 cycles of denaturation at 94°C for 30 sec, annealing at 𝑌 °C for 30 sec, and extension at 72°C for 45 sec, and a 20 min final extension step at 72°C. The annealing temperature ( 𝑇 𝑎 ) was changed based on the melting temperature ( 𝑇 𝑚 ) of each primer-pair (Table 1).

For each locus, amplification was confirmed by running 5 μL of the PCR products on 1.5% ethidium bromide containing agarose gels. The PCR products of the twelve primer-pairs were multiplexed into two panels, each of which containing six PCR products. In each panel, the size difference between the PCR products labeled with the same fluorescent dyes was at least 80 bp to avoid overlapping. The multiplex PCR products were then analyzed using an ABI Prism 3730 DNA Analyzer (Applied Biosystems) at Genomics Core Facility of the University of Gothenburg, Sweden.

2.3. Genotyping and Data Analysis

The allele peaks were visually inspected and then analyzed using PEAK SCANNER V1.0 software (Applied Biosystems) based on the internal Genescan-500 LIZ size standard. Each peak was considered as an allele at a codominant locus and the genotype of each individual at each locus was recorded. The Free Tree-Freeware program [35] was used to generate Nei’s standard genetic distance and for cluster analysis and bootstrapping. TreeView (Win32) 1.6.6 program [36] was used to view the trees. Analysis of molecular variance (AMOVA) was conducted using Arlequin ver. 3.5-2 [37].

3. Results

3.1. Total and within-Population Genetic Variation

Out of the 12 loci analyzed, eight loci were polymorphic whereas only one allele was detected across the 260 individuals analyzed in each of the remaining four loci. The four monomorphic loci were 838, DCM06, Sat235, and SSR06 (Table 2). The overall gene diversity for each polymorphic locus varied from 0.01 (Cam35) to 0.55 (SSR09). In addition to SSR09, the other loci with a relatively high level of gene diversity were CM5, Sat207, and Cam03 with 𝐻 𝑇 of 0.50, 0.50 and 0.54, respectively (Table 3).

tab3
Table 3: The partitioning of total gene diversity into within and among variety components for eight polymorphic SSR loci, and the number of alleles observed at each locus.

The total gene diversity ( 𝐻 𝑇 ) and the within-populations gene diversity ( 𝐻 𝑆 ), estimated based on Nei’s gene diversity [38], were 0.353 and 0.291, respectively (Table 3). The genetic diversity of each population ( 𝐻 L o c i ), which is the average gene diversity across the eight polymorphic loci, and the percent polymorphic loci (PPL) were also analyzed. 𝐻 L o c i ranged from 0.23 to 0.47, whereas %PL ranged from 0.33 to 0.58 (Table 1). At the variety level, the mean Nei’s gene diversity ranged from 0.24 (variety Maracaturra) to 0.37 (variety Catimor) with corresponding lowest and highest %PL of 0.33 and 0.52. Variety Catimor showed the highest gene diversity in five of the eight polymorphic loci (Table 3). The overall mean gene diversity and %PL per population were 0.29 and 0.42, respectively.

Population-specific rare alleles, with frequencies ranging from 0.025 to 0.1, were detected in five of the 26 populations. An 89 bp allele unique to population B3 was detected at locus Sat207 (Figure 1) at a frequency of 0.025. Similarly, a 97 bp allele was detected in populations CA3 and CM5 at this locus at the same frequency. The other populations bearing unique and rare alleles were B4 and CT7. These populations had unique alleles at loci SSR03 (frequency = 0.1) and Cam35 (frequency = 0.025), respectively.

fig1
Figure 1: The electrophoretogram showing the alleles at the SSR locus Sat207.
3.2. Genetic Variation among Populations and Groups

The genetic differentiation of populations was estimated based on gene diversity ( 𝐺 𝑆 𝑇 ; [39]) and AMOVA ( 𝐹 𝑆 𝑇 ; [37]). The overall mean 𝐺 𝑆 𝑇 and 𝐹 𝑆 𝑇 were 0.23 and 0.13, respectively (Tables 3 and 4). The estimates of 𝐺 𝑆 𝑇 varied from 0.00 (CM5) to 0.62 (Cam16) when calculated for each locus. The values of both 𝐺 𝑆 𝑇 and 𝐹 𝑆 𝑇 at locus CM5 were zero although each population has a gene diversity of 0.5. This is due to the fact that all individuals were heterozygous for the two alleles at this locus. The same is true for locus Sat207 except that few individuals in three of the 26 populations had additional rare alleles. At these two loci, both the total and within-population gene diversity was high but with no differentiation between the populations.

tab4
Table 4: SSR-based AMOVA for the 26 populations of arabica coffee: (A) without grouping the populations, (B) by grouping the populations according to varieties, (C) by grouping the populations according to regions of origin, and (D) by grouping the populations according to province of origin.

Overall, AMOVA revealed a highly significant genetic variation among populations ( 𝑃 < 0 . 0 0 0 1 ; Table 4) accounting for 13.5% of the total variation. The differentiation among varieties was also significant ( 𝐹 𝐶 𝑇 = 0 . 0 8 ; 𝑃 = 0 . 0 2 3 ) contributing 7.9% to the total genetic variation. The presence of rare alleles in four of the eight varieties contributed to the significant differentiation obtained. On the other hand, AMOVA revealed no significant variation among the two coffee growing regions and among the eight provinces ( 𝑃 > 0 . 4 ; Table 4). The pairwise AMOVA in the 26 populations revealed that each population was significantly differentiated from at least four populations. The most differentiated populations were CM2, CM3, and CM4, all of which belong to the variety Catimor. Pairwise 𝐹 𝑆 𝑇 showed that CM2 and CM3 were significantly differentiated from each other as well as from all other populations. Population CM4 was significantly differentiated from all populations except from B2 (Table 5). Population CA5 was significantly differentiated from only four populations (B2, CM2, CM3, and CM4). At variety level, Catimor, Catuai rojo, and Pacas were differentiated from each other and all other varieties (Pairwise 𝐹 𝑆 𝑇 ; Table 6).

tab5
Table 5: The Nei’s standard genetic distance (above the diagonal) and the pairwise 𝐹 𝑆 𝑇 (below the diagonal) between the 26 pairs of arabica coffee populations.
tab6
Table 6: The Nei’s standard genetic distance (above the diagonal) and the pairwise 𝐹 𝑆 𝑇 (below the diagonal) between the eight arabica coffee varieties.
3.3. Genetic Distance and Cluster Analysis

The Nei’s standard genetic distance between populations ranged from less than 0.001 (e.g., CA1 versus CA2) to 0.392 (B3 versus CM2) with the overall mean of 0.060 (Table 5). The genetic distance between the varieties ranged from 0.001 (Caturra versus Bourbon and Caturra versus Catuai amarillo) to 0.121 (Catimor versus Catuai amarillo) with the overall mean of 0.031 (Table 6). The Nei’s genetic distance-based cluster analysis revealed five clusters supported by moderate-to-high bootstrap values. Cluster I contained two populations from the variety Catimor (CM2 and CM3) with a bootstrap support of 100%. The 98% bootstrap supported Cluster II contained three populations (CM4, B2, and CT2), which belong to three different varieties. Similarly, cluster III comprised two populations (CR2 and P) from the Catuai-rojo and Pacas varieties with a 62% bootstrap support. Populations CM1, CT5, M, and PA, each of which belongs to different varieties, were placed under Cluster IV with a 65% bootstrap support. Cluster V is the largest cluster comprising 15 populations that were subclustered into three groups. However, the bootstrap support for the subclusters was low. In all clusters except cluster I, populations from more than one variety were clustered together showing a poor clustering of populations according to their varieties (Figure 2). At the variety level, Maracaturra and Pacamaras formed cluster I, whereas Catuai Amarillo, Bourbon, and Caturra formed cluster II with a 94% bootstrap support (Figure 3). The remaining three varieties remained solitary.

939820.fig.002
Figure 2: UPGMA phenogram for the 26 coffee populations based on Nei’s standard genetic distance. Numbers in front of the branches are bootstrap values.
939820.fig.003
Figure 3: UPGMA phenogram for the eight coffee varieties based on Nei’s standard genetic distance. Numbers in front of the branches are bootstrap values.

4. Discussion

4.1. The SSR Loci and Alleles

The difference in number of nucleotides between alleles obtained at the polymorphic loci in the present study indicates that the source of polymorphism was mainly the variation in number of repeat motifs of the SSRs. For example, CAM35 is a hexanucleotide repeat SSR [33], and the size of the alleles obtained in the present study was 207 bp and 213 bp. Similarly, the size of the two alleles of CM5, a trinucleotide repeat SSR locus [15], was 91 bp and 94 bp. However, the difference in size of the alleles observed at locus Sat207 appeared to be due to a combination of differences in the repeat motif and other types of variation, such as indels in the flanking sequences, as the difference in length between the alleles varied from three to four nucleotides (82 bp, 89 bp, 93 bp, 97 bp; Figure 1).

In addition to their application for analysis of genetic diversity, SSR markers have several other applications that include their use as markers for desirable traits. Among the SSRs used in the present study, Sat207 and Sat235 were reported to be tightly linked to locus Ck-1 that carries a major gene conferring resistance to the coffee berry disease (CBD) with Sat235 more closely linked to the gene than Sat207 [23]. CBD is a fungal disease caused by Colletotrichum kahawae that may cause severe damage in arabica coffee.

The SSR Sat235 was monomorphic across the 26 populations, and it is less likely that it can be a useful marker for genetic linkage analysis of Ck-1 in Nicaraguan coffee. On the other hand, Sat207 was polymorphic with two major alleles (ca 82 bp and 93 bp) and two rare alleles (ca 89 bp and 97 bp). Taking into consideration the amphidiploid nature of arabica coffee, it is most likely that the 82 bp allele on one hand and the other three alleles on the other hand originated from different progenitor genomes of arabica coffee (Figure 1).

If variation exists at the Ck-1 locus in Nicaraguan arabica coffee that gives resistance to CBD, the polymorphism detected at Sat207 is worth considering during the development of molecular markers linked to the resistance trait. Since arabica coffee is generally considered susceptible to CBD (e.g., [23, 40]), resistant genotypes should be rare and thus it would be interesting to evaluate the genotypes carrying the two rare alleles for resistance to this disease. Developing CBD-resistant arabica coffee varieties through identification of mutants is a method of choice, as it is simple and straight forward as compared to transferring resistance genes from other coffee species that requires crossing with donor genotypes followed by backcrossing to restore desirable traits. Considering that alleles of the same size at locus sat207 are identical by decent, it would also be interesting to compare the allele linked to resistance to CBD [23] with the allele introgressed to arabica coffee from robusta coffee [12], as this helps to assign the alleles to the two progenitor genomes.

Another interesting locus to discuss is CM5. Two alleles were detected at this locus and all the 260 individual plants studied were heterozygous for the two alleles. Baruah et al. [15] also identified only two alleles at this locus in arabica coffee. Several authors have reported a high cross-species transferability of SSR markers including the EST-SSRs within the genus Coffea (e.g., [15, 18, 41]). Given that arabica coffee is autogamous [3], the 100% heterozygosity obtained at this locus can only be explained by its amphidiploid nature. The two alleles should have been originated from different arabica coffee ancestral genomes. Baruah et al. [15] obtained 70% heterozygosity in arabica coffee at this locus, unlike the present study, suggesting that the two arabica coffee genomes of some genotypes carried the same alleles due to homoplasy. At this and other similar loci, fixed heterozygosity is the result when the two homoeologous loci are monomorphic and homozygous within the studied populations. The results clearly suggest the lack of recombination between the chromosomes of the two ancestral genomes due to the amphidiploid nature of arabica coffee.

Cubry et al. [19] obtained only two alleles per locus in the study that involved sixty SSR loci, and based on this they treated their data as diploid species data. However, a maximum of two alleles per locus is not always the case in arabica coffee, as shown in the present study. For example, three alleles were obtained at locus CAM03 in most individual plants analyzed. Three alleles per genotype were also observed at locus Sat207. These SSRs are reliable evidence that shows the presence of loci bearing nonrecombining alleles in arabica coffee representing the homoeologous loci from the two progenitor genomes. This supports the amphidiploid nature of the allotetraploid arabica coffee previously reported based on cytological evidence [2, 42].

4.2. The within- and among-Population Genetic Variation

The narrow genetic base of arabica coffee caused by rigorous selection during domestication and breeding has been reported by several authors (e.g., [1416, 18, 19]). For example, Cubry et al. [19] reported a mean of 2.1 alleles per locus for arabica coffee, which was the lowest among the Coffea species they studied. This is comparable with 2.3 alleles/locus obtained in the present study. Similarly, Moncada and McCouch [16] reported a mean of 1.9 alleles per locus.

In the present study, the mean Nei’s total ( 𝐻 𝑇 ) and within-population ( 𝐻 𝑆 ) gene diversity were estimated to be 0.35 and 0.29, respectively. The estimates for these parameters were 0.22 and 0.07, in that order, for the RAPD-based study of the Ethiopian arabica coffee by Anthony et al. [13]. Similarly, Cubry et al. [19] reported a mean gene diversity of 0.30 for the arabica coffee material they studied using SSR markers, whereas Aga et al. [43] reported an 𝐻 𝑇 of 0.37 using ISSR markers in Ethiopian forest coffee. Thus, the level of genetic variation in Nicaraguan arabica coffee is comparable to that previously reported from several countries and regions. The presence of the major SSR alleles across all the populations in very high frequencies in the present study suggests a narrow gene pool of arabica coffee in Nicaragua in line with previous reports. This suggests some difficulties in finding genotypes bearing desirable traits, such as resistance to diseases and pests within the domesticated arabica coffee gene pool.

4.3. The Arabica Coffee Varieties in Nicaragua

The cluster analysis of the SSR data for the 26 populations revealed that, in most cases, the clustering pattern of the populations was not in line with their varietal classification. The principal coordinate analysis (PCoA) of the 260 individual plants (data not shown) revealed the presence of divergent genotypes in populations B2, B3, CM3, CT2, CT5, and PA, which partly explains the poor clustering of populations according to their variety of origin. Given that arabica coffee is an autogamous species; such a poor clustering pattern of populations according to variety of origin is somewhat unexpected. However, the processes through which these varieties were developed may partly explain the lack of a clear differentiation between the populations of the different varieties. For example, the variety Caturra was developed from mutant genotypes of the variety Bourbon (http://www.coffeeresearch.org/coffee/varietals.htm). A significant differentiation between these two varieties may not be expected due to the relatively short time elapsed since the development of the variety Bourbon, especially at selectively neutral loci that include most of the SSRs used in the present study. Similarly, Catuai was the result of a cross between Mundo Novo and Caturra and thus there may not be a clear genetic differentiation between the Caturra and Catuai varieties at this stage. A relatively close relationship between the Bourbon, Caturra, and Catuai varieties can be observed from Figure 3. However, other possible factors, such as some degree of gene flow between varieties through cross pollination, might have also contributed to the population genetic structure obtained.

Among the eight coffee varieties we studied, the highest genetic diversity was recorded in the variety Catimor (GD = 0.37 and PPL = 0.52). This variety is interesting not only because of its high genetic diversity but also because two of its populations (CM2 and CM3) were significantly differentiated from all the other populations (Table 5). At locus 471, the alleles recorded in CM2 and CM3 were different from those in the other populations, excluding CM4. The relatively high diversity in Catimor can be partly explained by the fact that it was the result of a cross between the variety Caturra of arabica coffee and the Timor hybrid, which is a natural hybrid between arabica and robusta coffee [21, 23]. Robusta coffee has been reported to have a relatively high genetic diversity compared to arabica coffee in several studies (e.g., [16, 19]). Thus, a wise use of the genetic diversity in the locally adapted populations of the variety Catimor in coffee breeding programs in Nicaragua is very important. The presence of allelic variation at several loci in the other coffee varieties suggests the significance of using the existing genetic variation in these varieties in the hybrid breeding program to develop superior and improved varieties.

Root-knot Nematodes (RKN) of the genus Meloidogyne cause major damage in coffee worldwide, and coffee breeding for durable resistance to RKN is now a major goal in coffee producing countries [28]. Meloidogyne exigua and M. incognita are known RKN attacking arabica coffee in Nicaragua [44]. The best method to reduce the damage caused by RKN in coffee is through developing resistant varieties. Bertrand et al. [45] reported resistant arabica coffee genotypes to M. arabicida in Costa Rica, which suggests that identification of resistant genotypes to M. exigua and M. incognita from arabica coffee in Nicaragua may be possible. In addition, the resistance gene Mex-1 identified in robusta coffee and successfully introgressed into arabica coffee by Noir et al. [25] suggests the need to give attention to arabica coffee varieties developed through hybridization of the two cultivated Coffea species, such as the variety Catimor for pest and disease resistance.

Overall, the level of genetic diversity of arabica coffee in Nicaragua is generally low and is comparable to that previously reported for arabica coffee from other countries and regions. Therefore, it should be promoted through crossing with other closely related species such as robusta coffee. In addition, the presence of rare alleles in some populations suggests the need to explore such populations in order to identify mutants bearing desirable traits. The significant differentiation between most Nicaraguan arabica coffee varieties suggests that varieties grown in the country should be analyzed for resistance/tolerance to major biotic and abiotic stresses. On the other hand, the absence of a significant differentiation between the coffee populations based on regions of origin suggests that germplasm collecting missions should prioritize the representation of coffee varieties over coffee growing regions in the country.

Acknowledgments

The authors wish to thank the Swedish International Development Agency (Sida) for financing this work. They are grateful to Anna Zborowska for her assistance in the laboratory work. None of the authors has financial relations with any of the commercial identities mentioned in the paper.

References

  1. S. N. Raina, Y. Mukai, and M. Yamamoto, “In situ hybridization identifies the diploid progenitor species of Coffea arabica (Rubiaceae),” Theoretical and Applied Genetics, vol. 97, no. 8, pp. 1204–1209, 1998. View at Publisher · View at Google Scholar · View at Scopus
  2. P. Lashermes, M. C. Combes, J. Robert et al., “Molecular characterisation and origin of the Coffea arabica L. Genome,” Molecular and General Genetics, vol. 261, no. 2, pp. 259–266, 1999. View at Publisher · View at Google Scholar · View at Scopus
  3. A. P. Davis, R. Govaerts, D. M. Bridson, and P. Stoffelen, “An annotated taxonomic conspectus of the genus Coffea (Rubiaceae),” Botanical Journal of the Linnean Society, vol. 152, no. 4, pp. 465–512, 2006. View at Publisher · View at Google Scholar · View at Scopus
  4. K. Tesfaye, T. Borsch, K. Govers, and E. Bekele, “Characterization of Coffea chloroplast microsatellites and evidence for the recent divergence of C. arabica and C. eugenioides chloroplast genomes,” Genome, vol. 50, no. 12, pp. 1112–1129, 2007. View at Publisher · View at Google Scholar · View at Scopus
  5. J. M. Waller, M. Bigger, and R. A. Hillocks, “Coffee as a crop and commodity,” in Coffee Pests, Diseases and Their Management, pp. 1–17, CABI, Nosworthy Way Oxfordshire , UK, 2007.
  6. MIFIC (Ministerio de Fomento, Industria y Comercio), “El Café en Nicaragua. Dirección de Políticas Externas (DPCE),” Departamento de Análisis Económico, Managua, Nicaragua, pp. 47, 2005.
  7. A. E. Kuhl, “Nicaragua y su café,” Primera edición, Hispamer, Managua, Nicaragua, pp. 376, 2004.
  8. T. D. Beuchelt and M. Zeller, “Profits and poverty: certification's troubled link for Nicaragua's organic and fairtrade coffee producers,” Ecological Economics, vol. 70, no. 7, pp. 1316–1324, 2011. View at Publisher · View at Google Scholar · View at Scopus
  9. C. Rivas, “El café en Nicaragua: Análisis y descripción del comportamiento del rubro,” Banco Central de Nicaragua, Managua, Nicaragua, pp. 45, 2008.
  10. F. Guharay, J. Monterrey, D. Momnteroso, and C. Staver, “Manejo integrado de plagas en el cultivo del café,” CATIE, Managua, Nicaragua, pp. 272, 2000.
  11. ICO, 2011, http://www.ico.org/trade_statistics.asp.
  12. P. Lashermes, M. C. Combes, C. Ansaldi, E. Gichuru, and S. Noir, “Analysis of alien introgression in coffee tree (Coffea arabica L.),” Molecular Breeding, vol. 27, no. 2, pp. 223–232, 2011. View at Publisher · View at Google Scholar · View at Scopus
  13. F. Anthony, B. Bertrand, O. Quiros et al., “Genetic diversity of wild coffee (Coffea arabica L.) using molecular markers,” Euphytica, vol. 118, no. 1, pp. 53–65, 2001. View at Publisher · View at Google Scholar · View at Scopus
  14. F. Anthony, M. C. Combes, C. Astorga, B. Bertrand, G. Graziosi, and P. Lashermes, “The origin of cultivated Coffea arabica L. varieties revealed by AFLP and SSR markers,” Theoretical and Applied Genetics, vol. 104, no. 5, pp. 894–900, 2002. View at Publisher · View at Google Scholar · View at Scopus
  15. A. Baruah, V. Naik, P. S. Hendre, R. Rajkumar, P. Rajendrakumar, and R. K. Aggarwal, “Isolation and characterization of nine microsatellite markers from Coffea arabica L., showing wide cross-species amplifications,” Molecular Ecology Notes, vol. 3, no. 4, pp. 647–650, 2003. View at Publisher · View at Google Scholar · View at Scopus
  16. P. Moncada and S. McCouch, “Simple sequence repeat diversity in diploid and tetraploid Coffea species,” Genome, vol. 47, no. 3, pp. 501–509, 2004. View at Publisher · View at Google Scholar · View at Scopus
  17. M. Silvestrini, M. G. Junqueira, A. C. Favarin et al., “Genetic diversity and structure of Ethiopian, Yemen and Brazilian Coffea arabica L. accessions using microsatellites markers,” Genetic Resources and Crop Evolution, vol. 54, no. 6, pp. 1367–1379, 2007. View at Publisher · View at Google Scholar · View at Scopus
  18. R. K. Aggarwal, P. S. Hendre, R. K. Varshney, P. R. Bhat, V. Krishnakumar, and L. Singh, “Identification, characterization and utilization of EST-derived genic microsatellite markers for genome analyses of coffee and related species,” Theoretical and Applied Genetics, vol. 114, no. 2, pp. 359–372, 2007. View at Publisher · View at Google Scholar · View at Scopus
  19. P. Cubry, P. Musoli, H. Legnaté et al., “Diversity in coffee assessed with SSR markers: Structure of the genus Coffea and perspectives for breeding,” Genome, vol. 51, no. 1, pp. 50–63, 2008. View at Publisher · View at Google Scholar · View at Scopus
  20. B. Bertrand, M. X. Peña Durán, F. Anzueto et al., “Genetic study of Coffea canephora coffee tree resistance to Meloidogyne incognita nematodes in Guatemala and Meloidogyne sp. nematodes in El Salvador for selection of rootstock varieties in Central America,” Euphytica, vol. 113, no. 2, pp. 79–86, 2000. View at Publisher · View at Google Scholar · View at Scopus
  21. B. Bertrand, F. Anthony, and P. Lashermes, “Breeding for resistance to Meloidogyne exigua in Coffea arabica by introgression of resistance genes of Coffea canephora,” Plant Pathology, vol. 50, no. 5, pp. 637–643, 2001. View at Publisher · View at Google Scholar · View at Scopus
  22. L. Villain, J. Aribi, G. Réversat, and F. Anthony, “A high-throughput method for early screening of coffee (Coffea spp.) genotypes for resistance to root-knot nematodes (Meloidogyne spp.),” European Journal of Plant Pathology, vol. 128, no. 4, pp. 451–458, 2010. View at Publisher · View at Google Scholar · View at Scopus
  23. E. K. Gichuru, C. O. Agwanda, M. C. Combes et al., “Identification of molecular markers linked to a gene conferring resistance to coffee berry disease (Colletotrichum kahawae) in Coffea arabica,” Plant Pathology, vol. 57, no. 6, pp. 1117–1124, 2008. View at Publisher · View at Google Scholar · View at Scopus
  24. T. Tscharntke, Y. Clough, S. A. Bhagwat et al., “Multifunctional shade-tree management in tropical agroforestry landscapes—a review,” Journal of Applied Ecology, vol. 48, no. 3, pp. 619–629, 2011. View at Publisher · View at Google Scholar · View at Scopus
  25. S. Noir, F. Anthony, B. Bertrand, M. C. Combes, and P. Lashermes, “Identification of a major gene (Mex-1) from Coffea canephora conferring resistance to Meloidogyne exigua in Coffea arabica,” Plant Pathology, vol. 52, no. 1, pp. 97–103, 2003. View at Publisher · View at Google Scholar · View at Scopus
  26. E. Alpizar, H. Etienne, and B. Bertrand, “Intermediate resistance to Meloidogyne exigua root-knot nematode in Coffea arabica,” Crop Protection, vol. 26, no. 7, pp. 903–910, 2007. View at Publisher · View at Google Scholar · View at Scopus
  27. T. Leroy, F. Ribeyre, B. Bertrand et al., “Genetics of coffee quality,” Brazilian Journal of Plant Physiology, vol. 18, no. 1, pp. 229–242, 2006. View at Publisher · View at Google Scholar · View at Scopus
  28. E. V. S. Albuquerque, R. M. D. G. Carneiro, P. M. Costa et al., “Resistance to Meloidogyne incognita expresses a hypersensitive-like response in Coffea arabica,” European Journal of Plant Pathology, vol. 127, no. 3, pp. 365–373, 2010. View at Publisher · View at Google Scholar · View at Scopus
  29. E. Aga, Molecular genetic diversity study of forest coffee tree (Coffea arabica L.) populations in Ethiopia: implications for conservation and breeding, Ph.D. thesis, Swedish University of Agricultural Sciences, Acta Universitatis Agriculturae Sueciae, 2005, pp. 1652–6880, 2005:79.
  30. M. P. Maluf, M. Silvestrini, L. M. de C. Ruggiero, O. G. Filho, and C. A. Colombo, “Genetic diversity of cultivated Coffea arabica inbred lines assessed by RAPD, AFLP and SSR markers systems,” Scientia Agricola, vol. 62, no. 4, pp. 366–373, 2005.
  31. E. S. N. Vieira, E. V. de R. V. von Pinho, M. G. G. Carvalho, D. G. Esselink, and B. Vosman, “Development of microsatellite markers for identifying Brazilian Coffea arabica varieties,” Genetics and Molecular Biology, vol. 33, no. 3, pp. 507–514, 2010. View at Scopus
  32. E. Bekele, M. Geleta, K. Dagne et al., “Molecular phylogeny of genus Guizotia (Asteraceae) using DNA sequences derived from ITS,” Genetic Resources and Crop Evolution, vol. 54, no. 7, pp. 1419–1427, 2007. View at Publisher · View at Google Scholar · View at Scopus
  33. P. S. Hendre, R. Phanindranath, V. Annapurna, A. Lalremruata, and R. K. Aggarwal, “Development of new genomic microsatellite markers from robusta coffee (Coffea canephora Pierre ex A. Froehner) showing broad cross-species transferability and utility in genetic studies,” BMC Plant Biology, vol. 8, article 51, 2008. View at Publisher · View at Google Scholar · View at Scopus
  34. L. W. Ballard, P. S. Adams, Y. Bao, et al., “Strategies for genotyping: effectiveness of tailing primers to increase accuracy in short tandem repeat determinations,” Journal of Biomolecular Techniques, vol. 13, no. 1, pp. 20–29, 2002.
  35. A. Pavlícek, S. Hrdá, and J. Flegr, “FreeTree—freeware program for construction of phylogenetic trees on the basis of distance data and Bootstrap/Jackknife analysis of the tree robustness, application in the RAPD analysis of genus frenkelia,” Folia Biologica, vol. 45, no. 3, pp. 97–99, 1999. View at Scopus
  36. R. D. M. page, “TreeView: an application to display phylogenetic trees on personal computers,” Computer Applications in the Biosciences, vol. 12, no. 4, pp. 357–358, 1996. View at Scopus
  37. L. Excoffier, G. Laval, and S. Schneider, “Arlequin ver. 3.0: an integrated software package for population genetics data analysis,” Evolutionary Bioinformatics Online, vol. 1, pp. 47–50, 2005.
  38. M. Nei, “Estimation of average heterozygosity and genetic distance from a small number of individuals,” Genetics, vol. 89, no. 3, pp. 583–590, 1978. View at Scopus
  39. M. Nei, “Analysis of gene diversity in subdivided populations,” Proceedings of the National Academy of Sciences of the United States of America, vol. 70, no. 12, pp. 3321–3323, 1973. View at Scopus
  40. L. Hein and F. Gatzweiler, “The economic value of coffee (Coffea arabica) genetic resources,” Ecological Economics, vol. 60, no. 1, pp. 176–185, 2006. View at Publisher · View at Google Scholar · View at Scopus
  41. V. Poncet, M. Rondeau, C. Tranchant et al., “SSR mining in coffee tree EST databases: Potential use of EST-SSRs as markers for the Coffea genus,” Molecular Genetics and Genomics, vol. 276, no. 5, pp. 436–449, 2006. View at Publisher · View at Google Scholar · View at Scopus
  42. W. R. Clarindo and C. R. Carvalho, “First Coffea arabica karyogram showing that this species is a true allotetraploid,” Plant Systematics and Evolution, vol. 274, no. 3-4, pp. 237–241, 2008. View at Publisher · View at Google Scholar · View at Scopus
  43. E. Aga, E. Bekele, and T. Bryngelsson, “Inter-simple sequence repeat (ISSR) variation in forest coffee trees (Coffea arabica L.) populations from Ethiopia,” Genetica, vol. 124, no. 2, pp. 213–221, 2005. View at Publisher · View at Google Scholar · View at Scopus
  44. I. Herrera, T. Bryngelsson, A. Monzón, and M. Geleta, “Identification of coffee root-knot nematodes based on perineal pattern, SCAR markers and nuclear ribosomal DNA sequences,” Nematologia Mediterranea, vol. 39, no. 2, pp. 101–110, 2011.
  45. B. Bertrand, G. Ramirez, P. Topart, and F. Anthony, “Resistance of cultivated coffee (Coffea arabica and C. canephora) trees to corky-root caused by Meloidoqyne arabicida and Fusarium oxysporum, under controlled and field conditions,” Crop Protection, vol. 21, no. 9, pp. 713–719, 2002. View at Publisher · View at Google Scholar · View at Scopus