Characterization and genetic diversity of Coffea canephora accessions in a germplasm bank in Espírito Santo, Brazil

Abstract:

The state of Espírito Santo is the major producer of Coffea canephora in Brazil. Knowledge of genetic reserves is fundamental to plant breeding. Therefore, the present study characterized and analyzed the genetic diversity of 600 C. canephora accessions from the germplasm bank of Incaper based on 38 traits evaluated in 24-30-month-old plants. Further, the predominant descriptors or traits were identified, and high phenotypic variability was determined. Genetic distances for the grouped (Gower), quantitative, and qualitative datasets were 0.48, 0.61, and 0.92, respectively, with accessions 76 (Conilon) and 407 (Robusta) being the most divergent ones at Incaper. In clustering using the Tocher optimization method, 30 groups were formed, with three accessions introduced from Epamig’s Robusta collection being the most dissimilar ones. Graphical dispersion analysis using the principal coordinate method revealed the predominance of three groups formed by the Robusta, Conilon, and hybrid Robusta × Conilon genotypes.

Keywords:
Breeding; selection; Conilon coffee; genetic divergence; genetic resources.

INTRODUCTION

Coffee cultivation, as an important activity of the agricultural sector, is highly relevant in the context of the socioeconomic development of Brazil. In particular, in the state of Espírito Santo, coffee production has become a major component of the agribusiness sector, occupying an area of approximately 460.000 hectares (CONAB 2020). Coffea canephora is the second most cultivated species in the world, representing 40% of overall coffee production, and Espírito Santo stands out as the largest Brazilian producer of this species, designated in the State as Conilon coffee (Ferrão et al. 2019).

Coffee (Eudicotyledoneae) belongs to the genus Coffea of the angiosperm family Rubiaceae, with 124 species cataloged in the literature to date (Davis et al. 2011). Among these, only C. arabica and C. canephora are economically importance, while the others are of fundamental relevance in genetic breeding programs for hybridization and gene transfer (Ferrão et al. 2019, Mistro et al. 2019).

C. canephora is a diploid, perennial, allogamous species, with easy vegetative propagation and self-incompatibility of the gametophytic type (Conagin and Mendes 1961, Ferrão et al. 2017c, Moraes et al. 2018), and it is known to have denominations, such as Kouillou (Conilon) and Robusta (Charrier and Berthaud 1988). At the beginning of its expansion in the world, the species was known as Robusta. Based on records, it was introduced to Espírito Santo in Brazil in 1912, with the receipt of being the Kouillou variety, which was called Conilon. It became significantly widespread after the 70s, with the mother plants selected by the farmers themselves or introduced over the years, enabling the establishment of populations with wide genetic variability (Ferrão et al. 2019).

In 1985, the Institute of Research, Technical Assistance and Rural Extension of the State of Espírito Santo (Incaper), initiated a research program to explore the great variability of C. canephora through diverse selection and breeding strategies. In the 35 years of research, nine clonal and two seminal cultivars have been developed and made available to producers in Espírito Santo. These cultivars are the basis of the existing crops in the state. Incaper also maintains a germplasm bank of C. canephora in the field, with a large number of accessions, including cultivars, progenies, hybrids, and clones. Analysis of genetic variability among accessions, together with the maintenance and adequate characterization of germplasm, is vital to exploit the genetic diversity and trait reserve of the present collection for developing new genotypes through selection and recombination strategies (Ferrão et al. 2019). Characterization involves the evaluation of phenotypic descriptors defined a priori for the species, together with the assessments of biotic and abiotic factors, quality, and productive potential, among other variables (Ferrão et al. 2017a).

To this end, the present study characterized 600 C. canephora accessions from the germplasm bank of Incaper, maintained at three locations in Espírito Santo, for different agronomic traits based on specific descriptors, combined with the analysis of genetic variability and divergence.

MATERIAL AND METHODS

Genetic materials and evaluation sites

Incaper has maintained a C. canephora germplasm bank at the experimental farm in Marilândia, Espírito Santo, Brazil, since 1988. In 2016, actions to renew, duplicate, and expand this collection were initiated. Between January and May 2017, the collection was planted in the field in three coffee growing regions of Espírito Santo, including the experimental farms in Marilândia (FEM), Sooretama (FES), and Bananal do Norte (FEBN), located in the northwest, northeast, and south of the state, respectively. In the field, 3-5 plants of each of the 600 accessions were planted at a spacing of 3.0 m × 1.20 m in a single plot. Management and fertilization followed the crop recommendations by Incaper, and included various activities such as cycle pruning and irrigation for the maintenance of the initial establishment stand; no pest and disease control was performed (Ferrão et al. 2017b). The accessions are selected progenies from crops in different municipalities of Espírito Santo and southern Bahia, progenies from controlled crosses at Incaper and genotypes introduced from the Robusta collections.

Assessments, characterization, and data collection

From May to December 2019, 31 qualitative and 7 quantitative traits were evaluated at each location using 24-30-month-old plants. The qualitative traits were assessed based on coffee descriptors (Fazuoli et al. 1994, IPGRI 1996, Brasil 2000), evaluating variables related to the plant, branches, leaves, cycle, fruits, and seeds as well as those related to response to pests, diseases, and drought, using categorical scales (described in Tables 1 and 2). The following quantitative traits were derived from evaluations at harvest and post-harvest: grain yield at the first harvest (60 kg bag processed per hectare); percentage of floats; ratio of cherry coffee weight to processed coffee weight and that of natural coffee weight (dry at 11% moisture) to processed coffee weight (2 kg sample of cherry coffee obtained at the harvest); and percentage of average flat grains with sieves above 13, mocha grains with sieves above 10, and small grains (residue) (flat below 13 and mocha below 10).

Data analysis

The present data represent the mean of values from the three locations explored. For traits 1 to 26, which originated from multi-categorical evaluations based on categorical scales, the consistency of data from the three sites was analyzed. When divergence between data was observed, new data were collected and the results were revalidated.

Descriptive analysis was performed based on the data matrix of 600 accessions, which showed uniform development and 38 traits (qualitative and quantitative). The frequency of genotypes classified in each descriptor and their predominance in the germplasm were analyzed for each trait. Genetic divergence for 562 accessions was analyzed based on the complete set of qualitative and quantitative traits evaluated in the same phase using multivariate statistical procedures (Regazzi and Cruz 2020). The qualitative data (multi-categorical) were analyzed using the simple coincidence index method, and the quantitative data were analyzed using the standardized Euclidean mean distance. The matrices of the two datasets were grouped using the mean Gower method. In each dataset, the cophenetic correlations and genetic distances between the more and less similar genotypes were estimated. The genotypes were grouped with the joint Gower matrix using the modified Tocher optimization method, unweighted pair group method using arithmetic averages (UPGMA), and graphical dispersion analysis using the principal coordinate method. Statistical analyses were performed using GENES (Cruz 2016) and R (R Core Team 2019).

RESULTS AND DISCUSSION

Characterization of genotypes based on descriptors and complementary traits

Based on the data obtained through the evaluation of genotypes at the three sites, the average value of each trait was defined for each genotype, and descriptive analysis was performed to determine the predominance of each descriptor in Incaper’s germplasm bank. Table 1 presents the average frequencies (in percentage) of 26 phenotypic descriptors in the studied accessions: plants (n = 3), branches (n = 2), leaves (n = 8), fruits (n = 5), seeds (n = 6), and cycle (n = 2). Table 2 presents the average frequencies of 12 complementary traits in the studied accessions: response to diseases (n = 2), pests (n = 2), and drought (n = 1) and quantitative traits related to production (n = 1) and after harvest (n = 6). A summary of description of each trait, together with the discussion on the most relevant aspects, is presented below.

Table 1
Average frequencies of 26 phenotypic descriptors from the characterization of 600 Coffea canephora accessions in Incaper’s germplasm bank obtained through evaluation at three sites

Descriptor/Trait		Category	Code	Average (%)	Descriptor/Trait		Category	Code	Average (%)
Plant	Shape	Cylindrical (Cy)	1	3.01	Fruit	Size	Very small	1	0.00
		Conic (Co)	2	0.00			Small	3	4.72
		Cy-Co	[3]	96.99			Medium	[5]	72.85
		Inverted conic	4	0.00			Large	7	20.57
							Very large	9	1.85
	Height	Very low	1	3.17		Shape	Round	1	23.78
		Low	3	13.52			Oblong	[2]	62.39
		Medium	[5]	52.25			Elliptical	3	13.83
		High	7	21.70
		Very high	9	9.36
	Canopy diameter	Very small	1	0.67		Color (ripe bean phase)	Yellow	1	0.00
		Small	3	1.84			Orange-red	2	11.13
		Medium	5	15.03			Medium red	[3]	72.34
		Large	[7]	60.77			Dark red	4	16.53
		Very large	9	21.70
Branch	Primary and lateral: internode length	Short	3	4.34		Adherence to the branch	Low	[3]	58.18
		Medium	[5]	66.44			Medium	5	34.06
		Long	7	29.22			High	7	7.76
	Degree of plagiotropic branching	Low	3	2.34		Juiciness of the mesocarp (ripe fruits)	Dry	3	4.89
		Medium	[5]	53.09			Medium	[5]	65.94
		High	7	44.57			Succulent	7	29.17
Leaf	Length	Short	3	1.84	Seed	Length	Short	3	20.34
		Medium	[5]	66.61			Medium	[5]	52.37
		Long	7	31.55			Long	7	27.29
	Width	Narrow	3	5.01		Width	Narrow	3	12.54
		Medium	[5]	75.79			Medium	[5]	60.17
		Wide	7	19.20			Wide	7	27.29
	Shape	Elliptical	1	35.06		Thickness	Thin	3	5.59
		Oval	2	4.01			Medium	[5]	63.73
		Lanceolate	[3]	60.93			Thick	7	30.68
	Color in young phase	Green	1	7.68		Endosperm color	Yellow	1	39.12
		Bronze	2	34.89			Green	[2]	60.88
		Bronze-Green	[3]	57.43
		Purple	4	0.00
	Color in adult phase	Light green	1	1.34		Silver skin shade	Light	1	41.65
		Dark green	[2]	98.66			Dark	[2]	58.35
		Purple	3	0.00
	Curling of edges	Absent	1	0.00		Degree of silver skin adhesion	Weak	3	39.12
		Present	[9]	100			Medium	[5]	52.28
							Strong	7	8.60
	Leaf edge intensity	Low	3	8.85	Cycle	Ripening (over 50% ripe fruits)	Very early	1	4.41
		Medium	[5]	72.62			Early	[3]	46.95
		High	7	18.53			Medium	5	32.03
							Late	7	15.93
							Very late	9	0.68
	Depth of secondary vein	Shallow	3	10.85		First harvest	Early	[3]	53.36
		Medium	[5]	70.28			Medium	5	32.03
		Deep	7	18.86	Late	7	16.61

Table 2
Average frequencies of 12 phenotypic descriptors from the characterization of 600 Coffea canephora genotypes in Incaper’s germplasm bank obtained through evaluation at three sites

Descriptor/Trait	Category	Code	Average (%)
Response to diseases: Hemileia vastatrix	Resistant	1	31.55
	Moderately resistant	3	6.34
	Moderately susceptible	5	14.53
	Susceptible	7	41.07
	Highly susceptible	9	6.51
Response to diseases: Cercospora coffeicola	Resistant	1	5.34
	Moderately resistant	3	11.36
	Moderately susceptible	5	71.29
	Susceptible	7	11.85
	Highly susceptible	9	0.17
Response to pests: Leaf miners	Resistant	1	0.33
	Moderately resistant	3	16.19
	Moderately susceptible	5	64.44
	Susceptible	7	18.70
	Highly susceptible	9	0.33
Response to pests: Rosette cochineal	Resistant	1	24.28
	Moderately resistant	3	57.34
	Moderately susceptible	5	14.67
	Susceptible	7	3.54
	Highly susceptible	9	0.17
Response to drought	Low tolerance	1	35.68
	Moderate tolerance	2	56.11
	High tolerance	3	8.21
Average productivity (60 kg bags per hectare)	Low (<35)	1	68.56
	Medium (35.1-50)	3	25.93
	High (50.1-70)	5	4.97
	Very high (>70.1)	7	0.54
Percentage of floaters (%)	Low (0 -20)	1	3.74
	Medium (21.1-40)	3	71.30
	High (41.1-60)	5	23.89
	Very high (>60)	7	1.07
Improvement yield (ratio of cherry coffee weight to processed coffee weight)	High (3.1-4.0)	1	18.44
	Medium (4.1-5.0	3	53.37
	Low (5.1-6.0)	5	18.09
	Very low (>6.0)	7	10.00
Yield of natural coffee (ratio of dry natural coffee weight to processed coffee weigh)	High (1.5-2.0)	1	54.71
	Medium (2.01-2.5)	3	34.28
	Low (2.51-3.0)	5	8.17
	Very low (>3.0)	7	2.84
Percentage of flat beans, with sieves 13, 15, and 17	Very low (0-20)	1	0.18
	Low (2.1-40)	3	8.35
	Medium (40.1-60)	5	44.58
	high (60.1-.80)	7	46.89
Percentage of mocha beans, with sieves 10, 11, and 12	Very low (0-15)	1	8.20
	Low (15-30)	3	54.19
	Medium (30-45	5	31.73
	High (> 45)	7	5.88
Percentage of residue (bottom of the sieve)	0-15	1	65.84
	15.1-30	3	27.76
	30.1-45	5	4.80
	45.1-72	7	1.60

Plants: Regarding the shape of the plants, almost all accessions (96.99%) were classified as cylindrical-conic, with no variation in this characteristic among the three sites. Among the Robusta accessions, only two were classified as cylindrical. Regarding plant height, majority of the genotypes showed medium height (52.25%), followed by high (21.70%), low (13.52%), very high (9.36%), and very low (3.17%). Of the accessions classified as low at all evaluated sites, three clones were identified as Robusta, which generally has larger plants. Regarding the diameter of the plant, many accessions had a large diameter (60.77%), followed by very large (21.70%) and medium (15.03%). Over 82% of the accessions had a large or very large canopy diameter, and over 40% of such accessions were identified as Robusta. Regarding the internode length of the stem, 66.44% accessions were classified as medium, 29.22% as long, and 4.34% as short.

Branches: Most of the accessions had a medium internode length (66.44%), followed by long (29.22%) and short (4.34%). The degree of plagiotropic branching was high (53.09%) or medium (44.57%).

Leaves: Over 90% of the accessions had leaves classified as medium or long in terms of length and medium or wide in terms of width. In general, C. canephora plants bear three types of leaves: elliptical, oval, and lanceolate. Across the three sites studies, 60.93% leaves were lanceolate, 35.06% were elliptical, and 4.01% were oval, and there were variations in leaf shape among the three sites. Leaf color in the young phase was bronze-green (57.43%), bronze (34.89%), or green (7.68%), whereas that in the adult phase was predominantly dark green (98.66%). There were no purple leaves. Leaf edge intensity and depth of the secondary vein were classified as medium in over 70% accessions.

Fruits: Regarding fruit size, respectively 4.72%, 72.85%, 20.57%, and 1.85% were classified as small, medium, large, and very large. Most fruits were considered oblong (62.39%), followed by round (23.78%) and elliptical (13.83%). Regarding color, 72.34% accessions bore medium red fruits, 16.53% bore dark red fruit, 11.13% bore orange red fruit, and none bore yellow fruits. In ripe fruit, 58.18% accessions showed a low degree of adherence to the branch, 34.06% showed average adherence, and 7.76% showed strong adherence. Regarding juiciness of the mesocarp of ripe fruits, only 4.89% accessions were classified as being dry, while the rest were classified as being moderately juicy (65.94%) or succulent (29.14%).

Seeds: The seed length was medium, long, and short for 52.37%, 27.29%, and 20.34% accessions, respectively. The seed width was classified as medium, wide, and narrow for 60.17%, 27.29%, and 12.54% accessions, respectively. The seed thickness was classified as medium for 63.73% accessions. The endosperm color was green in 60.88% accessions and yellow in the rest (39.12%). The silver skin shade was dark and light in 58.35% and 41.65% accessions, respectively. The degree of adhesion of the silver skin was weak, medium, and strong in 39.12%, 52.28%, and 8.60% accessions, respectively.

Cycle: The ripening cycle, evaluated at the first harvest when over 80% fruits were ripe, showed early (46.95%) and medium (32.03%) maturation in most accessions, followed by late (15.93%), very early (4.41%), and very late (0.68%). Overall, fruit maturation in all accessions was earlier in FES, located in the northeastern region of Espírito Santo. The genotypes that matured before April 15, from April 15 to May 31, from June 1 to June 30, and from July 1 to July 31 were classified as very early, early, medium, late, and very late cycle, respectively.

Responses to diseases, pests, and drought: The responses of genotypes to diseases and pests were assessed on a scale from 1 to 9, where 1 indicates resistant genotypes and 9 indicates highly susceptible genotypes, and the average of the highest score at the three sites was calculated. Drought response was assessed on a scale from 1 to 3, where 1 indicates susceptible genotypes and 3 indicates tolerant genotypes. Disease incidence was higher at FES located in the northeastern region. Specifically, 37.89% genotypes were resistant (scores 1 to 3), 14.53% were moderately susceptible, and 47.58% were susceptible (scores 7 and 9) to Hemileia vastatrix infection. Majority of the accessions were moderately susceptibility to Cercospora coffeicola infection, with 71.29% accessions achieving a score of 5. Regarding pests, 64.44% accessions were moderately susceptible to coffee leaf miners (score 5), but majority were moderately resistant (57.34%) and some were resistant (24.28%) to rosette cochineals. Regarding drought tolerance, only 8.21% accessions were highly tolerant (score = 3), while the rest showed moderate (56.11%) to low (35.68%) tolerance. Disease and pest prevalence was high (at least at some locations) and was helpful in characterizing the germplasm in terms of resistance and susceptibility.

Quantitative traits (complementary): Quantitative traits were presented as the average of values at the three sites, representing the data of the first harvest. On average, 68.56% accessions showed productivity below 35 bags ha^-1 and only 0.54% showed productivity exceeding 70 bags ha^-1. The proportion of floaters varied from 0% to 100% in FEM and FES and from 0% to 60% in FEBN. The average proportion of floaters across the three sites was between 21% and 40%. Regarding post-harvest coffee yield, in majority of the accessions, the ratio of cherry coffee weight to processed coffee weight ranged from 4.1 to 5.0 (53.37%) and that of dry natural to processed coffee weight ranged from 1.5 to 2.0 (54.71%). For the classification of grains, the sum of flat grains retained in sieves 13, 15, and 17; the sum of mocha grains retained in sieves 10, 11, and 12; and the residue (grains retained at the bottom of these sieves) were considered. The average proportions of flat grains, mocha grains, and residue were 57%, 27%, and 17.19%, respectively. Given that a high proportion of flat grains is a desirable trait to select genetic material for plant breeding, most genotypes studied here appear suitable as genetic material, with the proportion of flat grains exceeding 50% and that of mocha grains being below 30%.

Overall, the frequency distributions exhibited a wide variability. Akpertey et al. (2019) characterized Robusta species based on morphology and found that few genotypes differed in terms of the most prevalent traits.

Genetic divergence of accessions in Incaper’s germplasm bank

Genetic divergence was analyzed using multivariate statistics with the data of 562 accessions evaluated for 38 traits (qualitative and quantitative). The pairs of genotypes with greater similarity and dissimilarity differed among the three data matrices (qualitative, quantitative, and grouped). Nonetheless, the magnitude of cophenetic correlation coefficients was comparable (0.614-0.652) (Table 3). The greatest genetic distance was observed in the qualitative data matrix (0.920), followed by the quantitative data (0.613) and grouped data (0.472) matrices. The shortest genetic distance was observed in the quantitative data matrix (0.027), between accessions 191 and 368, which were selected in different years and from different locations. Accessions 76 and 407 were identified as the most divergent in this study, with a distance of 0.920. These accessions were also phenotypically distinct at the field level, exhibiting typical traits of Conilon (76) and Robusta (407). This result also demonstrates that this pair of accessions presents the greatest genetic distance ever observed in the studies based on phenotypic and/or molecular data from Incaper's reasearch program (Fonseca et al. 2006, Ferrão et al. 2009, Ferrão et al. 2017d, Ferrão et al. 2018, Giles et al. 2018, Senra et al. 2020). Genotype 76 exhibits superior agronomic characteristics and has been used as material in clonal cultivars and seeds developed and released by Incaper.

Table 3
Cophenetic correlations and minimum and maximum genetic distances of 562 Coffea canephora accessions in Incaper’s germplasm bank collection evaluated based on the qualitative, quantitative, and grouped data matrices using the simple coincidence index, standardized Euclidean mean distance, and Gower method, respectively

Data matrix	Cophenetic correlation	Maximum distance	Accessions with the maximum distance	Minimum distance	Accessions with the minimum distance
Qualitative	0.614	0.920	76 (V05) and 407 (IAC 18)	0.034	293 (416/97) and 347 (475/97)
Quantitative	0.652	0.613	33 (32/86) and 547 (16/2006)	0.027	191 (71/89) and 368 (01/2008)
Grouped	0.618	0.472	167 (03Bahia) and 206 (198/89)	0.045	117 (87/89) and 182 (102/89)

Clustering of materials using the modified Tocher optimization method formed 30 groups (Table 4). The first group included the highest number of accessions (191) with the greatest distance (0.187). The most similar accessions were 117 (87/89) and 182 (102/89), representing genotypes selected in the same year (1989) from farms of producers in the northern region of Espírito Santo through phenotypic mass selection. In contrast, the most dissimilar accessions were 507, 470, and 463, representing Robusta genotypes introduced from of the Agricultural Research Company of Minas Gerais (Epamig), Minas Gerais, Brazil.

Table 4
Grouping of 562 Coffea canephora accessions in Incaper’s germplasm bank using the modified Tocher optimization method based on the matrix obtained using the Gower method

Group	Accessions1
1	117, 182, 177, 157, 271, 176, 45, 173, 221, 169, 197, 149, 57, 325, 244, 93, 97, 293, 350, 187, 36, 315, 356, 118, 92, 68, 279, 535, 19, 27, 369, 95, 103, 17, 52, 113, 13, 343, 346, 109, 34, 349, 110, 41, 26, 213, 288, 192, 543, 277, 38, 40, 347, 190, 331, 58, 85, 49, 12, 212, 46, 332, 21, 231, 186, 86, 379, 30, 371, 5, 262, 151, 251, 202, 382, 274, 31, 236, 183, 43, 116, 378, 232, 122, 220, 156, 32, 209, 374, 353, 53, 94, 82, 227, 126, 127, 341, 25, 158, 235, 132, 294, 60, 168, 108, 336, 73, 140, 299, 364, 544, 536, 24, 542, 65, 178, 259, 386, 91, 298, 217, 141, 146, 261, 48, 383, 7, 138, 345, 249, 354, 128, 296, 248, 253, 198, 264, 11, 120, 310, 191, 51, 348, 286, 537, 516, 78, 316, 42, 297, 225, 358, 165, 89, 189, 263, 250, 303, 272, 164, 560, 360, 323, 292, 154, 144, 276, 289, 340, 214, 387, 247, 321, 270, 327, 506, 56, 324, 174, 112, 538, 184, 76, 180, 105, 204, 188, 66, 278, 160, 258
2	63, 269, 355, 6, 59, 228, 390, 308, 295, 268, 319, 314, 333, 275, 339, 334, 359, 185, 254, 71, 218, 8, 267, 306, 362, 205, 257, 290, 203, 193, 181, 266, 260, 139, 98, 376
3	287, 380, 389, 155, 459, 283, 549, 367, 285, 394, 533, 426, 478, 365, 440, 511, 447, 512, 476, 398, 444, 531, 460, 518, 417, 434, 396, 395, 448, 513, 462, 445, 485, 472, 409, 412, 427, 457, 397, 408, 466, 473, 418, 491, 525, 552, 519, 393, 414, 483, 415, 532, 475, 520, 493, 504, 555, 527, 503, 411, 502, 424, 432, 150, 443, 438, 496
4	401, 402, 400, 404, 522, 410, 419, 523, 497, 425, 464, 498, 468, 488, 437, 406, 436, 416, 480, 553, 558, 431, 106, 428
5	50, 80, 4, 233, 28, 70, 67, 143, 55, 377, 172, 29, 81, 222, 88, 370, 307, 195, 255, 373, 361, 125, 230, 368, 162, 337, 133, 14, 3, 130, 562, 163, 74, 201, 148, 153, 304, 111, 72, 330, 326, 312, 381, 102, 252, 335, 280, 210, 352, 338, 385, 101, 161, 234, 170, 96
6	495, 556, 505, 557, 550, 554, 224
7	429, 433, 456, 405, 430, 450, 458, 455, 442, 422, 451, 517, 420, 515
8	500, 501, 509, 490, 528, 551, 152, 211, 454, 384, 391, 208, 366
9	124, 238, 237, 16, 546, 363, 77, 241, 83, 559, 134, 90
10	99, 179, 313, 240, 351, 100, 136, 44
11	75, 167, 246, 245, 39, 302
12	20, 484, 1, 216
13	492, 521, 487, 529, 471, 486, 467, 441, 469
14	145, 166, 196, 291, 273, 215, 79, 159, 357, 226, 9
15	423, 540, 541, 489, 372, 481, 514, 449, 526, 317, 561, 421, 69, 342, 530, 524, 23
16	446, 494, 407, 452, 439, 461, 301
17	47, 219, 37, 18, 147, 119, 375, 107, 35, 121, 129, 87
18	206, 539, 547, 242, 548, 104
19	477, 545, 388, 435, 243, 10
20	305, 322, 311, 62, 309, 61, 281, 64, 207, 22
21	123, 142, 329, 318, 137, 171, 256, 84, 194, 320
22	175, 199, 15, 508, 54, 499, 300, 200
23	344, 453, 284, 403, 399, 392, 135, 482, 474, 479
24	115, 465, 229, 114
25	2, 328
26	33, 510, 265
27	223, 239
28	131, 282, 534, 413
29	463, 470
30	507

• ¹
  Robusta accessions 392-461 were obtained from IAC; accessions 462-532 were introduced from Epamig; and accession 554-562 were selected by Incaper. Other accessions were identified as Conilon.

The results of grouping using UPGMA were consistent with those of grouping using the modified Tocher optimization method, but illustrative graphical representation was difficult due to the large number of genotypes. However, in dispersion analysis, the genotypes were distributed in the four quadrants of the graph (Figure 1) based on the principal coordinate method. The accessions classified as the Robusta-type were concentrated in the quadrants on the left and were represented by numbers 392-461 (introduced from IAC), 462-532 (introduced from Epamig), and 554-562 (selected by Incaper). The accession identified as the Conilon-type or Robusta × Conilon hybrids were concentrated on the right and in the middle of the graph, respectively. The component clones of Incaper’s clonal cultivars were spatially well distributed in the middle and on the left, representing the variability among them.

Figure 1
Graphical dispersion of 562 Coffea canephora accessions in Incaper’s germplasm bank assessed using principal coordinates analysis. Robusta and Conilon accessions are shown in blue and pink, respectively. Component clones of cultivars are identified with numbers.

Overall, our results corroborate the findings reported in literature. As such, many Conilon genotypes are classified in the SG1 subgroup, originating from interspecific hybridizations over the years between genetic materials from the Guinean and Congolese groups. Initially, Berthaud (1986), Montagnon et al. (1992), and Dusset et al. (1999) divided C. canephora into two variable groups, namely Guinean and Congolese. Subsequently, within the Congolese group, subgroups SG1, SG2, B, C, and Ug were identified. Previous studies with 519 different genetic materials, including Brazilian Conilon accessions, identified five distinct groups (Montagnon et al. 2012), with the Brazilian Conilon placed in the SG1 subgroup. Solórzano et al. (2017) analyzed 1,491 C. canephora plants grown in different living collections in Ecuador and compared them with 29 genotypes representing the major diversity groups identified within the species. The authors found that the materials introduced from Brazil mainly belonged to SG1. This reinforces the previous reports that within the SG1 group, there are genetic materials of the Kouillou (Conilon), Robusta, and Conilon× Robusta hybrid genotype.

The accessions in Incaper’s germplasm collection are numbered according to the order of their inclusion from 1986 to 2016. Many accessions are obtained via phenotypic selection from farms in Espírito Santo and southern Bahia performed by Incaper researchers, with the purpose of sourcing the natural genetic variability present in crops resulting from natural crossing over the years. The results showed random grouping of many accessions, regardless of the year and place of selection, characterizing the well distributed variability of the species in the region.

These data and knowledge related to the genetic structure indicate that the Brazilian germplasm of C. canephora represents only a small portion of the total diversity of the species (Ferrão et al. 2019). Thus, additional efforts to collect and introduce accessions are warranted to expand the genetic basis as well as to maintain and characterize germplasm collections for exploring the genetic potential of the species in the light of constant demand for improved cultivars adapted to environmental and technological changes. Importantly, Incaper’s clonal cultivars of Conilon coffee comprise a group of diverse genotypes with superior agronomic characteristics. Thus, to maintain the genetic basis of this species in the state of Espírito Santo, the producers must use the set of clones of each clonal cultivar and, whenever possible, multiple cultivars (clonal, seminal, or both) for plantation within their rural properties.

ACKNOWLEDGMENTS

The authors thank Consórcio Pesquisa Café, the Espírito Santo Research and Innovation Support Foundation (Fapes), and the National Council for Scientific and Technological Development (CNPq) for their financial support in the form of research and scientific research grants.

REFERENCES

• Akpertey A, Anim-Kwapong E and Ofori A (2019) Assessment of genetic diversity in Robusta coffee using morphological characters. International Journal of Fruit Science 19: 276-299.
• Berthaud J (1986) Les ressources génétiques pour L’amériolations des caféiers africains diploïds. Evaluation de larichesse génétique des populations sylvestres et de ses mécanismes organisateures. Conséquences pour l’application. ORSTOM, Paris, 379p. (Collection Trauvaux at Documents, 188).
• Brasil (2000) Secretaria de Apoio Rural e Cooperativismo. Serviço Nacional de Proteção de Cultivares. Portaria, nº 2, 17 de novembro de 2000. Diário Oficial da República Federativa do Brasil, Brasília, n. 223, p. 6-7, 21 nov. 2000, Seção 1.
• Charrier A and Berthaud J (1988) Principles and methods in coffee plant breeding: Coffea canephora Pierre. In Clarke RJ and Macrae R (eds) Coffee: agronomy. Elsevier Applied Science, London, p. 167-195.
• CONAB - Companhia Nacional de Abastecimento (2020) Follow-up of the Brazilian coffee crop in Brazil: second survey. Available at <Available at https://www.conab.gov.br/info-agro/safras/cafe >. Accessed on September 30, 2020.
  » https://www.conab.gov.br/info-agro/safras/cafe
• Conagin CHTM and Mendes AJT (1961) Pesquisas citológicas e genéticas em três espécies de Coffea: autoincompatibilidade em C. Canephora Pierre ex Froehner Bragantia 20: 787-804.
• Cruz CD (2016) Genes software-extended and integrated with the R, Matlab and Selegen. Acta Scientiarum Agronomy 38: 547-552.
• Davis AP, Tosh J, Ruch N and Fay MF (2011) Growing coffee: Psilanthus (Rubiaceae) subsumed on the basis of molecular and morphological data, implications for the size, morphology, distribution and evolutionary history of Coffea Botanical Journal of the Linnean Society 167: 1-21.
• Dusset S, Lasharmes P, Anthony F, Montagnon C, Touslot P, Cambes MC, Berthaud J, Noirot M and Haman S (1999) Le cafeir, Coffea canephora In Hamon P, Seguin M, Perrier X and Glaszmann JC (eds) Diversité génetique des plantes tropicales cultivées. CIRAD, France, p. 175-194
• Fazuoli LC, Silvarolla MB, Camargo CEO, Pommer CV, Chiavegato EJ, Dall’Ort’ FAC, NagaI H, Godoy IJ and Veiga RFA (1994) Descritores mínimos para o registro institucional de cultivares: Café. Instituto Agronômico, Campinas, 8p. (Documentos IAC, 46).
• Ferrão LFV, Ferrão RG, Ferrão MAG, Fonseca AFA and Garcia AAF (2017a) A mixed model to multiple harvest-location trials applied to genomic prediction in Coffea canephora Tree Genetics & Genomes 13: 1-15.
• Ferrão LFV, Ferrão RG, Ferrao MAG, Fonseca AFA, Carbonetto P, Stephens M and Garcia AAF (2018) Accurate genomic prediction of Coffea canephora in multiple environments using whole-genome statistical models. Heredity 1: 1-16.
• Ferrão MAG, Ferrão RG, Fonseca AFA, Verdin Filho AC and Volpi PS (2019) Origin, geographical dispersion, taxonomy and genetic diversity of Coffea canephora In Ferrão RG, Fonseca AFA, Ferrão MAG and DeMuner LH (eds) Conilon coffee. Incaper, Vitória, p. 85-110.
• Ferrão MAG, Fonseca AFA, Ferrão RG, Barbosa WM and Souza EMR (2009) Genetic divergence in conilon coffee revealed by RAPD markers. Crop Breeding and Applied Biotechnology 9: 67-74.
• Ferrão MAG, Fonseca AFA, Ferrão RG, Volpi PS, Verdin Filho AC, Riva-Souza EM, Comério M and Kaulz M (2017d) Variabilidade genética de progênies híbridas de Coffea Canephora In Congresso brasileiro melhoramento de plantas. SBMP, Foz do Iguaçu, 9p.
• Ferrão MAG, Souza EMR, Fonseca AFA and Ferrão RG (2017b) Autoincompatilidade e produção sustentável do café conilon. In Ferrão RG, Fonseca AFA, Ferrão MAG and DeMuner LH (eds) Café conilon. Incaper, Vitória , p.177-188.
• Ferrão RG, Fonseca AFA, Bragança SM, Ferrão MAG and DeMunner LH (2017c) Café conilon . Incaper, Vitória , 702p.
• Fonseca AFA, Sediyama T, Cruz CD, Sakiyama NS, Ferrão MAG and Ferrão RG (2006) Divergência genética em café conilon (Coffea canephora Pierre). Pesquisa Agropecuária Brasileira 41: 599-605.
• Giles JAD, Partelli FL, Ferreira A, Rodrigues JP, Oliosi GE and Silva FHL (2018) Genetic diversity of promising ‘conilon’ coffee clones based on morpho-agronomic variables. Anais da Academia Brasileira de Ciências 90: 2437-2446.
• IPGRI - International Plant Genetic Resources Institute (1996) Descriptores del café (Coffea spp. E Psilanthus spp.). IPGRI, Roma, 36p. Available at Available at https://www.bioversityinternational.org/e-library/publications/detail/descriptors-for-coffee-coffea-spp-and-psilanthus-spp/ Accessed on September 1, 2016.
  » https://www.bioversityinternational.org/e-library/publications/detail/descriptors-for-coffee-coffea-spp-and-psilanthus-spp/
• Mistro JC, Resende MDV, Fazuoli LC and Vencovsky R (2019) Effective population size and genetic gain expected in a population of Coffea canephora Crop Breeding and Applied Biotechnology 19: 1-17.
• Montagnon C, Cubry P and Leroy T (2012) Amélioration génétique du caféier Coffea canephora Pierre: connaissances acquises, stratégies et perspectives. Cahiers Agricultures 21: 143-153.
• Montagnon C, Leroy T and Yapo AB (1992) Diversité génotypique et phénotypeque de quelques groupes de caféiers (Coffea canephora Pierre) en collection. Conséquences sur leur utilisation en sélection. Café Cacao Thé 36: 187-198.
• Moraes MS, Teixeira AL, Ramalho AR, Espi ndula MC, Ferrão MAG and Rocha RB (2018) Characterization of gametophytic self-incompatibility of superior clones of Coffea canephora Genetics and Molecular Research 17: 1-8.
• R Core Team (2019) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available at < Available at https://www.R-project.org/ >. Accessed on June 10, 2019.
  » https://www.R-project.org/
• Regazzi AJ and Cruz CD (2020) Análise multivariada aplicada. Editora UFV, Viçosa, 401p.
• Senra JFB, Ferrão MAG, Mendonça RF, Fonseca AFA, Ferrao RG, Volpi PS, Comerio M and Silva M (2020) genetic variability of access of the active germplasm bank of Coffea canephora of incaper in southern Espírito Santo. Journal of Genetic Resources 6: 172-184.
• Solórzano RGL, Bellis, FD, Leroy T, Plaza L, Guerrrero H, Subia C, Calderón D, Fernández F, Guarzón I, Lopez D and Vera D (2017) Revealing the diversity of introduced Coffea canephora germplasm in Ecuador: towards a national strategy to improve robusta. The Scientific World Journal 2017: 12

Publication Dates

• Publication in this collection
  13 Sept 2021
• Date of issue
  2021

History

• Received
  17 Dec 2020
• Accepted
  14 May 2021
• Published
  20 June 2021