# Starmaya: The first arabica F1 coffee hybrid produced using genetic male sterility

Frederic Georget, Lison Marie, Edgardo Alpizar, Philippe Courtel, Mélanie Bordeaux, Martin Hidalgo, Pierre Marraccini, Jean-Christophe Breitler, Eveline Dechamp, Clément Poncon, et al.

► To cite this version:

Frederic Georget, Lison Marie, Edgardo Alpizar, Philippe Courtel, Mélanie Bordeaux, et al. Star-maya: The first arabica F1 coffee hybrid produced using genetic male sterility. Frontiers in Plant Science, 2019, 10, pp.1344. 10.3389/fpls.2019.01344. hal-05174321

HAL Id: hal-05174321

https://hal.science/hal-05174321v1

Submitted on 21 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Starmaya: The First Arabica F1 Coffee Hybrid Produced Using Genetic Male Sterility

Frédéric Georget<fnref n="1" />,²*, Lison Marie<fnref n="1" />,², Edgardo Alpizar³, Philippe Courtel³,
Mélanie Bordeaux⁴, Jose Martin Hidalgo⁴, Pierre Marraccini<fnref n="1" />,², Jean-christophe Breitler<fnref n="1" />,²,
Eveline Déchamp<fnref n="1" />,², Clément Poncon³, Hervé Etienne<fnref n="1" />,² and Benoît Bertrand<fnref n="1" />,²

In the present paper, we evaluated the implementation of a seed production system based on the exploitation of male sterility on coffee. We studied specifically the combination between CIR-SM01 and Marsellesa® (a Sarchimor line), which provides a hybrid population called Starmaya. We demonstrated that the establishment of seed garden under natural pollination is possible and produces a sufficient amount of hybrid seeds to be multiplied efficiently and economically. As expected for F1 hybrid, the performances of Starmaya are highly superior to conventional cultivars. However, we observed some heterogeneity on Starmaya cultivar in the field. We confirmed by genetic marker analysis that the off-types were partly related to the heterozygosity of the CIR-SM01 clone and could not be modified. Regarding the level of rust resistance of Starmaya cv., we saw that it could be improved if Marsellesa was more fully fixed genetically. If so, we should be able to decrease significantly the percentage of rust incidence of Starmaya from 15 to 5%, which would be quite acceptable at a commercial level. Starmaya represents the proof of concept for the mass propagation of Arabica F1 hybrid seeds using male sterility. Finally, we discuss the possibility to increase the number of hybrid varieties produced by seed, exploring some initiatives to identify male sterility markers to induce male sterility on any conventional cultivar. This would definitively open up the universe of known Arabica cultivars to be used in breeding new F1 hybrids.

**Keywords:** Arabica coffee, F1 hybrid seed, male sterility, micropropagation, commercial production

## INTRODUCTION

Due to self-fertilization of the allotetraploid *Coffea arabica* (2n = 4X = 44) species, Arabica cultivars are traditionally disseminated by seeds. This is not possible for F1 Arabica hybrids, as they are highly heterozygous and their characteristics segregate in F2 progeny produced by self-pollination. The advantages of using Arabica hybrids were first demonstrated in Kenya with the complex hybrid cultivar Ruiru 11, developed in the 1970s and released in 1986 to be cultivated in Kenya (Walyaro 1983; Van der Vossen, 1985), then in Tanzania (Teri et al., 2005; Van der Vossen et al., 2015). Breeding strategies aiming to develop new F1 hybrid cultivars (from two parents) of Arabica coffee were proposed simultaneously in Ethiopia (Bellachew, 1997) and in Central America (Bertrand et al., 1997; Bertrand et al., 2011). In Central America, the strategy was based on crossing American cultivars with wild Ethiopian accessions. Like many other hybrid plants, *C. arabica* F1 hybrids possess genetic and agronomic advantages (Gallais, 2009), such as higher and more stable yields,

## OPEN ACCESS

**Edited by:**

Marcelino Perez De La Vega,
Universidad de León, Spain

**Reviewed by:**

Aaron P. Davis,
Royal Botanic Gardens, Kew,
United Kingdom
Eveline Teixeira Caixeta,
Embrapa Café, Brazil

**\*Correspondence:**

Frédéric Georget
frederic.georget@cirad.fr

**Specialty section:**

This article was submitted to
Plant Breeding,
a section of the journal
Frontiers in Plant Science

**Received:** 15 May 2019

**Accepted:** 27 September 2019

**Published:** 22 October 2019

**Citation:**

Georget F, Marie L, Alpizar E, Courtel P, Bordeaux M, Hidalgo JM, Marraccini P, Breitler J-c, Déchamp E, Poncon C, Etienne H and Bertrand B (2019) Starmaya: The First Arabica F1 Coffee Hybrid Produced Using Genetic Male Sterility. *Front. Plant Sci.* 10:1344. doi: 10.3389/fpls.2019.01344

<footnote n="1">¹CIRAD, UMR IPME, Montpellier, France, ²IPME, Université de Montpellier, IRD, CIRAD, Montpellier, France, ³Plant material, ECOM, Exportadora Atlantic, Managua, Nicaragua, ⁴FONDATION NICAFRANCE, Managua, Nicaragua</footnote>
more vigor, disease resistance, better cup quality, and adaptability to agroforestry systems (Bertrand et al., 2005; Bertrand et al., 2011). These authors showed that F1 hybrids produced 11-47% higher yield than the best cultivars along with significantly higher or identical 100-bean weight and performed identically for fertility. To assess whether using hybrids represents substantial genetic progress on terms of productivity and full-sun cropping systems, the new F1 hybrids were grown in the same conditions as the best American cultivars. The results showed that yields of hybrids were earlier and superior to those of American cultivars. The hybrids were also more stable than the American cultivars in all environments tested. In the agroforestry system, the mean yield of hybrids was 58% higher than that of American cultivars, while the mean yield of hybrids in the full-sun system was 34% higher. So, by combining more productivity, better adaptability to agroforestry system, more resistance, and better cup quality, F1 hybrids are a promising alternative to conventional cultivars. Once created, several methods exist to multiply these new F1 hybrids on a large scale.

As previously reported in Kenya and Ethiopia, the first method consisted in producing F1 seeds by hand pollination (Opile and Agwanda, 1993; Bellachew, 1997). In that case, and because Arabica flowers are hermaphrodite and autogenic, it was necessary to emasculate them before pollinating by hand. For coffee, this approach has proven difficult to implement on a large scale. The second method was to produce clonal F1 hybrid plants by different micropropagation technics. In Tanzania for example, through the full participation of some 600 farmer groups, the F1 hybrids were multiplied by rooted cuttings and grafting (Magesa et al., 2013). On the other hand, the ECOM/Centre de Coopération Internationale en Recherche Agronomique Pour le Développement (CIRAD) consortium produced an average of one million F1 Arabica hybrid coffee plants per year by somatic embryogenesis (SE) in Nicaragua (Etienne et al., 2012; Etienne et al., 2018). Although this approach worked technically, it was not economically viable due to its high costs linked to manpower and the infrastructures required (laboratories and nurseries) (Van der Vossen et al., 2015). We recently showed that the cost of F1 coffee plants could be reduced by the half, if plants were produced by a combination of SE and rooted mini-cuttings technics (Georget et al., 2015b). A 6-month-old nursery hybrid plant propagated from a somatic seedling and ready for the field, which has a production cost of around 0.60 USD, whereas the same hybrid plant produced by rooted mini-cuttings has a production cost of around 0.27 USD. However, the price of such F1 hybrids produced first by SE, then re-multiplied by rooted minicuttings, is still 1.5 times the cost of traditional seedling, which still remains too high and constitutes a major deterrent for small coffee farmers (Georget et al., 2017).

Lastly, the third method consisted in producing F1 seeds in seed gardens using a sterile male line as the female plant (Georget et al., 2015a). Male sterility is common in flowering plants, but its application in hybrid breeding and seed production is limited because of the inability to propagate a pure sterile male line for commercial hybrid seed production (Chang et al., 2016). There is therefore a limited number of publications reporting on the use of sterile male genotypes to produce hybrid seeds at commercial level. Numerous hybrid crops that are significant components

of world production (>50%) such as tomato and maize are produced by hand emasculation, while spinach and the cucurbits hybrid seeds are produced by monoecy and the brassicas by self-incompatibility (Havey, 2004). For other crops, such as rice, wheat, and sorghum, commercial quantities of hybrid seeds can only be obtained by using sterile male female parents created by chemical or genetic manipulations (Wu et al., 2016). In this case, the successful application of such systems for large-scale hybrid seed production depends on whether the sterile male female parent line can be multiplied efficiently and economically (Perez-Prat and Van Lookeren Campagne, 2002). Compared with other methods, these hybrid seeds offer the advantages of being produced in quantity, easily handled, stored, and transported, and are much more affordable for farmers. With an annual renewal rate estimated at over one billion trees for *C. arabica* (B. Bertrand, personal communication), this type of propagation by seeds should allow rapid and large-scale dissemination of Arabica F1 hybrids.

In this paper, we evaluated the implementation of a seed production system using the *C. arabica* sterile male CIR-SM01, and we evaluated the performance of the first F1 hybrid, “Starmaya,” obtained by male sterility in comparison with conventional cultivars.

## MATERIALS AND METHODS

### C. arabica Sterile Male CIR-SM01

In 2001, researchers from CIRAD found a natural tall growth habit mutant (called CIR-SM01) of *C. arabica* without any pollen production, by screening Ethiopian plants at the Centro Agronómico Tropical de Investigación y Enseñanza (CATIE) germplasm bank for coffee in Costa Rica. Plants of CIR-SM01 were then multiplied *in vitro* in clonal form by SE (Atlantic S.A. laboratory, Sebacco, Nicaragua). Somatic embryos derived from plantlets were first acclimatized in plastic tunnels for 2 months in trays, then hardened in nurseries for 4 months before planting in the field and further used as mother receptors (Figures 1A, B). After this process, it was checked that the sterile male phenotype was not affected by temperature or photoperiod-sensitive, as described in rice for example (Zhang et al., 1994). This mutation displayed complete and stable sterility, whatever the environmental conditions and the age of the plant (juvenile or mature). The molecular basis of this mutation is still unknown, but we found in a previous study (Dufour et al., 1997) that phenotypic expression of this male sterility is functional, occurring during microsporogenesis, and is the product of a malfunction of the tapetum in the early microspore stages. At the same time, we planted four F2 populations (1,250 plants), and we checked that the mutation was recessive, giving a proportion of 18.3% sterile male plants. This proportion may have been the product of a Mendelian recessive epistatic interaction with a 9:4:3 ratio.

### Marsellesa® Cultivar of *C. arabica*

The Marsellesa® cultivar (https://www.upov.int/pluto/en/) is a Sarchimor cultivar selected (F8) with complete resistance to
leaf rust disease (coffee leaf rust caused by *Hemileia vastatrix*) and to coffee berry disease (caused by *Colletotrichum kahawae*), as described in the World Coffee Research variety catalog (https://varieties.worldcoffeeresearch.org/). Compared with a traditional dwarf cultivar, such as Caturra red, Marsellesa® produces up to 10% more, with good organoleptic qualities. Seeds from Marsellesa® were taken from a plantation located on the “La Cumplida” farm (Matagalpa, Nicaragua) and sprouted for 45 days in a germinator (**Figures 1C, D**). The resulting seedlings were transferred to nursery bags for 4 months to be fully acclimatized, then planted in the field to be further used as pollen donors.

## Creation of the Starmaya Hybrid Cultivar

In a preliminary study, we first investigated the possibility of producing F1 hybrid seeds by crossing (by hand pollination) the sterile male CIR-SM01 separately with several *C. arabica* cultivars with a dwarf growth habit (e.g., Caturra red, Catuai 44, IAPAR59, and Marsellesa®) used as pollen donors and evaluating the agronomic performances of the corresponding F1 progenies under field growing conditions (data not shown). In this trial, all plants were fertile with better vegetative vigor than the control (Caturra red). Of the four combinations tested, the best cross was the “CIR-SM01 × Marsellesa®” F1 population, which we called “Starmaya.” This candidate cultivar

yielded 30% more of green beans than Marsellesa®, with a standard to good beverage quality and rust resistance. The combinations with Caturra red and Catuai 44 were discarded due to high susceptibility to rust, and the combination with IAPAR59 presented a cup quality inferior to the Caturra red standard.

This first step was carried out on a few plants, so we decided to evaluate the agronomic performance of this new cultivar on a larger scale by comparing it with conventional ones and to validate the concept of commercially producing Starmaya hybrid seeds using male sterility and natural pollination.

## Set up of the Starmaya Seed Garden

To produce hybrid seeds, we planted CIR-SM01 together with Marsellesa®, assuming that the pollen of Marsellesa® would be able to pollinate the sterile male naturally. The proportion was approximately one pollen donor for four sterile male plants used as mothers. The first experimental seed garden, isolated from alien pollen contamination, was planted in 2009 on the “La Cumplida” farm (950 m above sea level, district of Matagalpa, Nicaragua) on half a hectare (2,000 trees), with 1,600 CIR-SM01 somatic seedling plants and 400 Marsellesa® plants as pollen donors. The planting density was 1.25 m along the row and 2 m between rows, i.e., 4,000 plants per hectare (**Figure 2**).

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Production of Starmaya F1 *C. arabica* hybrid. (A) Pre-germinated embryos of sterile-male CIR-SM01 produced by somatic embryogenesis. (B) Plantlets of sterile-male CIR-SM01 derived from somatic embryos and acclimatized in trays during 4 months. (C) Seeds of *C. arabica* Marsellesa® variety. (D) 2-month-old plants of Marsellesa® variety obtained after seed germination.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Experimental seed garden used to produce seeds of Starmaya F1 hybrid of *C. arabica*. (A) Scheme of planting in seed garden to produce seeds of Starmaya F1 hybrid. The fertile plants (in blue) producing pollen correspond to Marsellesa® pure line variety. The sterile male plants corresponding to CIR-SM01 produce Starmaya seeds. In seed garden, the ratio is approximately one plant of Marsellesa for 4 plants of sterile male (SM) CIR-SM01. (B) Seed garden established in the field of "La Cumplida" farm located in Matagalpa (Nicaragua). (C) First seeds of Starmaya F1 hybrid harvested in seed garden in 2011. (D) First flowering (in 2014) of Starmaya F1 hybrid plots established at "La Cumplida" farm. (E) Starmaya plants in 2015. (F) Aspect of producing branches of Starmaya producing the first commercial yield in 2015.</figcaption></figure>

## Multi-Site Trial to Assess the Performances of the Starmaya F1 Hybrid

A multi-site trial was conducted across seven environments in Nicaragua for 5 years (2013 to 2018). The objective was to test the Starmaya F1 hybrid in multiple environments, in order to compare it with two conventional cultivars, Caturra red and Marsellesa®. The agronomic management plan was identical at each site. The trials were established at seven sites ranging from 710 to 1,250 m above sea level. Implementation was identical for each site. Each site had five blocks, and each cultivar plot comprised 20 plants of each genotype. The Starmaya plot and the two control plots were randomly distributed in each block. The spacing between rows was 2.0 m, and the spacing between trees along a row was 1.5 m, corresponding to a planting density of around 3,300 trees per hectare. The sites are described in Table 1.

Productivity was recorded over 3 years for each variety plot. The growth parameters were measured for five plants per

genotype per block on 2-year-old plants. Rust resistance was evaluated for 10 plants per genotype per block in 2017 and 2018.

Rust was evaluated on a scale of 0 to 4, with 0 = no symptoms, 1 = some rust spots without sporulation, 2 = majority of leaves with

<table><caption>**TABLE 1** | List of parcel locations where the hybrid Starmaya cultivar was tested in multi-environment trials in Nicaragua.</caption><thead><tr><td>Location</td><td>Department</td><td>Name of plot</td><td>Elevation (m.a.s.l.)*</td></tr></thead><tbody><tr><td>Las Colinas</td><td>Boaco</td><td>San Jorge</td><td>710</td></tr><tr><td>Las Joyas</td><td>Jinotega</td><td>La Ceiba</td><td>715</td></tr><tr><td>La cueva del Tigre</td><td>Matagalpa</td><td>La Laguna</td><td>850</td></tr><tr><td>Zaragoza</td><td>Jinotega</td><td>El Guabal</td><td>1030</td></tr><tr><td>Las Marias</td><td>Nueva Segovia</td><td>Lote Siete</td><td>1190</td></tr><tr><td>La Aurora</td><td>Matagalpa</td><td>Los Cipes</td><td>1240</td></tr><tr><td>Albania</td><td>Matagalpa</td><td>Validacion</td><td>1250</td></tr></tbody></table>

<footnote>*m.a.s.l., meters above sea level.</footnote>
rust spots without sporulation, 3 = majority of leaves with rust sporulation, and 4 = majority of leaves with rust sporulation and with defoliation. For plant morphology, the volume was estimated by the formula LP × H, where LP = length of the two longest primary branches (centimeters) and H = the height of the tree (centimeters).

For bean characteristics, we estimated the weight of 100 healthy green beans (W100), and we screened the samples (5 kg for each genotype) through sieves (S14 to S20). Cup quality was assessed by six expert testers from Illy Caffè, with eight descriptors (aroma, flavor, aftertaste, acidity, body, balance, overall, and defects). A final score was assigned to each sample according to the Specialty Coffee Association (SCA) protocol (http://www.scaa.org/?page=resources&d=cupping-standards).

To estimate the stability of Starmaya across the range of environments, productivity data were analyzed following the additive main effects and multiplicative interaction (AMMI) model. This model takes into account genotype and environment interactions (Purchase et al., 2000; Bose et al., 2014; Darai et al., 2017), with the Agricolae package in R. AMMI analyses were used to partition genotype × environment deviations into different interaction principal component axes (IPCA). AMMI analysis first fits additive effects for host genotypes and environments by the usual additive analysis of variance procedure and then fits multiplicative effects for G × E (genotype × environment) by a principal component analysis (PCA).

AMMI Model:

$$Y_{ij} = \mu + gi + ej + \sum_{k=1}^{n} \lambda_k \alpha_{ik} \gamma_{jk} + e_{ij}$$

where μ is the general mean. Yij is the yield of the i<sub>th</sub> genotype in the j<sub>th</sub> environment. gi is the i<sub>th</sub> genotype mean deviation. e<sub>j</sub> is the j<sub>th</sub> environment mean deviation. λ<sub>k</sub> is the square root of the eigenvalue of the PCA axis k, α<sub>jk</sub> and γ<sub>jk</sub> are the principal component scores for PCA axis k of the i<sub>th</sub> genotype and the j<sub>th</sub> environment, respectively. ε<sub>ij</sub> is the residual (Zobel, 1998).

The AMMI stability value (ASV) described by Purchase et al. (2000) was calculated as follows:

$$ASV = \sqrt{\left[ \frac{\text{IPCA1 sum of square}}{\text{IPCA2 sum of square}} (\text{IPCA1 score}) \right]^2 + (\text{IPCA2 score})^2}$$

where SSIPCA1/SSIPCA2 is the weight given to the IPCA1 value by dividing the interaction principal component analysis axis I (IPCA1) sum of squares by the IPCA2 sum of squares. The higher the IPCA score, either negative or positive, the more specifically adapted a genotype is to certain environments. Lower ASV scores indicate a more stable genotype across environments.

## Setup of a Starmaya Commercial Field

F1 seeds of Starmaya and Marsellesa® were planted at La Cumplida (2012) in a field of 8 ha (1.25 m along rows × 2.25 m between rows). The 8 ha were divided into four blocks of 2 ha. Each block was subdivided in two equal plots of 1 ha for each

genotype. We obtained a little production in 2014, but the first commercial production was in 2015. Four years after planting (2016), plant height, stem diameter, internode number, primary branch mortality, and primary branch length were measured on 300 plants randomly chosen for each genotype in the four blocks. Yield was measured as the quantity of green beans produced for each plot and expressed in bags of 46 kg ha⁻¹.

Resistance to leaf rust was also evaluated by analyzing 300 plants of both varieties randomly chosen in the four previously described blocks. Rust resistance incidence was assessed during the rainy season, considered as the peak of fungus infestation observed on the farm on susceptible cultivars (i.e., Caturra red cv.). A plant was considered as susceptible if sporulation zones were observed on the abaxial surfaces of the leaves. The color of young leaves and shape of the plants (dwarf or tall) were observed three times and repeated if necessary when the information was ambiguous, until a consensus was reached. Marsellesa® is considered completely resistant to the main races prevailing in Nicaragua, while CIR-MS01 is susceptible to this fungus.

The coffee beans of Marsellesa® and Starmaya were classified according to their defects, size distribution (sieve S13 to S20), and sensorial analyses. After the 2015–2016 harvest, 100 kg of cherries of each variety was processed by the wet method and was dried in the sun to approximately 11% moisture. The beans were then passed through a series of sieves, thereby being divided into classes of small (sieves 13 and 14), medium (15 to 17), and large (sieves 19 and 20) beans. The percentages of defects and exportable parchment coffee were evaluated for each cultivar following procedures routinely used by Exportadora Atlantic S.A factories at Sebaco in Nicaragua (ECOM group).

Sensorial analyses of Marsellesa® and Starmaya were performed twice with two trained tasters, based on the SCA protocol and using scales from 0 to 10 for the evaluated attributes. The Starmaya and Marsellesa® cultivars where then compared with Caturra red cultivar harvested from a commercial plot located close to the experiment (200 m). The data corresponded to the mean of the two sensorial evaluations.

## Simple Sequence Repeat Analysis

In order to identify the homozygous/heterozygous allele patterns, as well as any existence of alien alleles, DNA was extracted from the leaves of 10 plants randomly selected for both parents (♂ Marsellesa® and ♀CIR-SM01) and Starmaya F1 hybrids for further simple sequence repeat (SSR) analysis. In this case, leaf samples of CIR-SM01 were taken from plants grown in the seed garden and those of Starmaya and Marsellesa® from plants grown in the agronomic trial. The SSR primer pairs used corresponded to those previously reported by Anthony et al. (2002) and Gichuru et al. (2008) (Table 2).

The polymerase chain reaction (PCR) reactions were performed in a 15-μl final volume comprising 30 ng of genomic DNA, 2 μmol ml⁻¹ of magnesium dichloride, 300 nmol ml⁻¹ of deoxyribonucleotide triphosphates, and 7.5 μl of 2× PCR buffer (Type it, Qiagen), 1.0 μM each of forward and reverse primer (10 μM). Amplifications were carried out in a thermal cycler (Eppendorf) programmed at 94°C for 5 min (initial denaturation), followed by 94°C for 30 s, annealing temperature (depending on
<table><caption>**TABLE 2** | List of primers used in simple sequence repeat analysis. The primer sequences, the range size (in base pairs) of amplified fragment, and annealing temperatures used in polymerase chain reaction reactions are indicated for each primer pair.</caption><thead><tr><th>Primer name</th><th>Oligo name</th><th>Sequence</th><th>Fragment size range (bp)</th><th>Annealing temperature</th></tr></thead><tbody><tr><td rowspan="2">Sat 11</td><td rowspan="2">CO220</td><td>F</td><td>ACCCGAAAGAAAGAACCAA</td><td rowspan="2">143-145</td><td rowspan="2">35</td></tr><tr><td>R</td><td>CCACACAACTCCTCCTATTC</td></tr><tr><td rowspan="2">Sat 225</td><td rowspan="2">CO221</td><td>F</td><td>CATGCCATCATCAATTCCAT</td><td rowspan="2">265-298</td><td rowspan="2">35</td></tr><tr><td>R</td><td>TTACTGCTCATCATTCCGCA</td></tr><tr><td rowspan="2">Sat 24</td><td rowspan="2">CO223</td><td>F</td><td>GGCTCGAGATATCTGTTTAG</td><td rowspan="2">148-172</td><td rowspan="2">35</td></tr><tr><td>R</td><td>TTTAATGGGCATAGGGTCC</td></tr><tr><td rowspan="2">Sat 32</td><td rowspan="2">CO227</td><td>F</td><td>AACTCTCCATTCCCGCATTC</td><td rowspan="2">105-133</td><td rowspan="2">35</td></tr><tr><td>R</td><td>CTGGGTTTTCTGTGTTCTCG</td></tr><tr><td rowspan="2">Sat 47</td><td rowspan="2">CO229</td><td>F</td><td>TGATGGACAGGAGTTGATGG</td><td rowspan="2">117-150</td><td rowspan="2">35</td></tr><tr><td>R</td><td>TGCCAATCTACCTACCCCTT</td></tr><tr><td rowspan="2">Sat 235</td><td rowspan="2">CO222</td><td>F</td><td>TCGTTCTGTCATTAAATCGTCAA</td><td rowspan="2">227-261</td><td rowspan="2">35</td></tr><tr><td>R</td><td>GCAAATCATGAAAATAGTTGGTG</td></tr><tr><td rowspan="2">Sat 254</td><td rowspan="2">CO225</td><td>F</td><td>ATGTTCTTCGCTTCGCTAAC</td><td rowspan="2">202-220</td><td rowspan="2">35</td></tr><tr><td>R</td><td>AAGTGTGGGAGTGCTGCAT</td></tr><tr><td rowspan="2">Sat 29</td><td rowspan="2">CO226</td><td>F</td><td>GACCATTACATTTCACACAC</td><td rowspan="2">119-135</td><td rowspan="2">35</td></tr><tr><td>R</td><td>GCATTTTGTTGCACACTGTA</td></tr><tr><td rowspan="2">Sat 244</td><td rowspan="2">CO303</td><td>F</td><td>GCATGTGCTTTTTGATGTCGT</td><td rowspan="2">278-308</td><td rowspan="2">35</td></tr><tr><td>R</td><td>GCATACTAAGGAATTATCTGACTGCT</td></tr><tr><td rowspan="2">CBD-Sat 207</td><td rowspan="2">CO308</td><td>F</td><td>GAAGCCGTTTCAAGCC</td><td rowspan="2">83-93</td><td rowspan="2">35</td></tr><tr><td>R</td><td>CAATCTCTTTCCGATGCTCT</td></tr></tbody></table>

each pair) for 30 s, and 72°C for 1 min (extension) for 35 cycles, followed by a final extension step at 72°C for 5 min. PCR samples were then run on a capillary electrophoresis 3130XL with an internal standard (homemade). The SSR analysis was provided by the ADNid company (http://www.adnid.fr, Montpellier, France).

## RESULTS AND DISCUSSION

### Potential Starmaya F1 Hybrid Seed Production of the Seed Garden

Based on the results obtained from the CIR-SM01 and Marsellesa® seed garden field, we found that natural pollination was efficient for producing Starmaya F1 hybrid seeds. In terms of seed production, an average of 100–150 kg of F1 hybrid seeds (data not shown) was harvested per year from this experiment (0.5 ha). These quantities amounted to potential production of approximately half a million hybrid seeds per hectare of seed garden, which should enable the planting of over 100–120 ha of Starmaya F1 hybrid coffee each year. The sterile male CIR-SM01 produced poorly compared with commercial cultivars, which usually produce around two million seeds under similar agronomic conditions.

The cost of the seed produced by this sterile male progenitor was estimated to be about 0.07 USD per unit, corresponding to half the price of a coffee hybrid plantlet propagated by SE (Georget et al., 2017). This difference is mainly due to the fact that seed gardens are less expensive to set up and maintain than *in vitro* clonal propagation laboratories. Moreover, only a few micropropagation laboratories in the world are able to commercially produce such Arabica F1 hybrids, and none of them produces more than 2 million somatic embryo-derived plantlets per year, which represents a key constraint for the democratization of using such hybrids.

No reliable statistic exists on the seed requirements for renovating coffee plantations worldwide. According to our calculations, and considering a renovation rate of 2 to 5% per year for Arabica crops estimated to cover 7 million hectares, it would be necessary to produce 250 to 500 million trees annually. Of course, this greatly depends on the coffee prices of international markets. Under a high coffee price, more than 1 billion plants could be required yearly for renovation worldwide. However, this requirement should not exceed 200–400 million plants per year when coffee prices are low. Consequently, 50 to 100 million plants of *C. arabica* F1 hybrids should be produced per year to have a major impact with the dissemination of such hybrids at the international level (Bertrand B., personal communication). In that case, this could be achieved by using 100 to 200 ha of seed garden with a sterile male coffee accession such as CIR-SM01. According to our current knowledge, it would be impossible to produce the same amount of F1 hybrids by SE.

### Starmaya Multi-Site Variety Trial

Figure 3 shows the average yields for each genotype in kilograms of green beans per plant. The graph shows that the Starmaya hybrid produced more than traditional Caturra red and Marsellesa® cultivars. A Tukey's honestly significant difference test was applied to determine significant differences between averages. Starmaya produced significantly more compared with Caturra red and Marsellesa® (30% more per plant).

By AMMI models, we evaluated and provided a stability index, the ASV. Genotypes with a low ASV are considered stable. The ASV allowed us to identify that Starmaya was the more stable genotype having the lowest ASV value (ASV = 0.308), while Marsellesa® and Caturra red, with an ASV of 0.807 and 0.819, respectively, appeared less stable.
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Average weight of green coffee (kilogram) per plant, for each variety. Treatment with the same letter are not significantly different at P ≤ 0.05 (Tukey's honestly significant difference test).</figcaption></figure>

The tree volume of Starmaya (Figure 4) (0.89 m³) was greater than that of Marsellesa® and Caturra red (0.59 m³). The differences are mainly due to lateral length and reflect strong vegetative vigor. Starmaya was found to be less susceptible to rust than Caturra red, but more than Marsellesa®, due to the presence of some susceptible plants. As rust resistance background in Marsellesa® is conferred by dominant genes, we concluded that some Starmaya plants do not inherit rust resistance.

The W100 appeared higher for the Starmaya hybrid compared with the conventional cultivars (Figure 5). Starmaya had an

average weight of 100 healthy green beans 25% higher than Caturra red. The higher W100 weight indicated better filling during bean development.

Bean grade also plays an important role in marketing. The purchase price of bean is related to the density and size of the bean, and its healthy physical appearance. Beans are classified in 64th of an inch increments, with sizes between 14 to >19. Beans with screen sizes up to 15 are considered the highest quality in terms of size. Figure 6 shows, for each cultivar, the percentage of beans of size 16 to 20 and the percentage of defective beans.

<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Relationship between productivity and two parameters: tree volume and susceptibility to rust. Treatment with the same letter are not significantly different at P ≤ 0.05 (Tukey's honestly significant difference test).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Average weight (gram) of 100 healthy green beans per variety (W100). Treatment with the same letter are not significantly different at P ≤ 0.05 (Tukey's honestly significant difference test).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | Percentage of green beans of size 16 to 20 and defective green beans by variety. Treatment with the same letter are not significantly different at P ≤ 0.05 (Tukey's honestly significant difference test).</figcaption></figure>

Starmaya gave a higher percentage of green beans of size 16 to 20 than the two controls and 1% more defective beans.

As concerns sensorial perception, according to the SCA method, coffees with a score of less than 80 are considered to be of inferior quality. Coffee beverages that have a score between

80 and 84.99 are very good, and beyond 85, they are excellent or exceptional. Cup quality varies depending on genotypes and altitudes. **Table 3** shows that Starmaya was better rated than Caturra red and obtained scores over or equal to 80, being very good at the highest altitude (>83).

<table><caption>**TABLE 3** | Score of sensorial perception according to the Specialty Coffee Association method and made by Illy Caffè (Italy).</caption><thead><tr><th>Genotypes/Elevation</th><th>710 m.a.s.l*</th><th>1,190 m.a.s.l</th><th>1,240 m.a.s.l</th><th>1,250 m.a.s.l</th></tr></thead><tbody><tr><td>Caturra red</td><td>72.5</td><td>77.7</td><td>78.9</td><td>74.4</td></tr><tr><td>Starmaya</td><td>79.2</td><td>80</td><td>80</td><td>83.1</td></tr></tbody></table>

<footnote>* m.a.s.l, meters above sea level.</footnote>
# Phenotypic and Yield Evaluation of Starmaya Under Agronomic Conditions

The phenotypic characteristics of Marsellesa® and Starmaya were evaluated in year 4 and are reported in **Table 4**. For all the phenotypic parameters studied (except the number of primary branch internodes at level 15 from the apex), the comparison showed that Starmaya plants were more vigorous than those of Marsellesa®. For example, Starmaya plants were always significantly taller and had a larger collar diameter than Marsellesa®. Moreover, Starmaya plants also had longer primary branches and larger internodes than those of Marsellesa®.

Regarding the yield assessment, the results presented in **Table 5** show that Starmaya produced 47% more than Marsellesa® in 2015–2016 (first commercial harvest) and 29% more for the second commercial harvest (2016–2017). The yield of Starmaya was still better than that of Marsellesa® for the third harvest (2017–2018). After three large yields, the experimental field was pruned. Considering the average harvests, Starmaya produced up to 35%

<table><caption>**TABLE 4** | Phenotypic comparison of *C. arabica* Marsellesa® Sarchimor cultivar and Starmaya F1 hybrid (CIR-SM01 × Marsellesa®).</caption><thead><tr><th>Parameters</th><th>Marsellesa®</th><th>Starmaya</th></tr></thead><tbody><tr><td>Stem diameter at collar (mm)</td><td>38.6 ± 4.0<sup>b</sup></td><td>44.0 ± 6.1<sup>a</sup></td></tr><tr><td>Plant height (cm)</td><td>177.5 ± 16.5<sup>b</sup></td><td>206.3 ± 15.9<sup>a</sup></td></tr><tr><td>Nb total of internodes</td><td>33.3 ± 3.2<sup>b</sup></td><td>37.9 ± 2.7<sup>a</sup></td></tr><tr><td>Primary branch length (cm)</td><td>40.5 ± 7.0<sup>b</sup></td><td>48.3 ± 9.6<sup>a</sup></td></tr><tr><td>level 10 from the apex (Branch 1)</td><td></td><td></td></tr><tr><td>Nb of primary internode (Branch 1)</td><td>10.2 ± 1.8<sup>b</sup></td><td>11.33 ± 2.1<sup>a</sup></td></tr><tr><td>Primary branch length (cm)</td><td>41.2 ± 7.7<sup>b</sup></td><td>51.9 ± 9.0<sup>a</sup></td></tr><tr><td>level 10 from the apex (Branch 2)</td><td></td><td></td></tr><tr><td>Nb of primary internode (Branch 2)</td><td>10.5 ± 1.9<sup>b</sup></td><td>11.7 ± 2.1<sup>a</sup></td></tr><tr><td>Primary branch length (cm)</td><td>59.1 ± 6.8<sup>b</sup></td><td>64.2 ± 9.3<sup>a</sup></td></tr><tr><td>level 15 from the apex (Branch 1)</td><td></td><td></td></tr><tr><td>Nb of primary internode (Branch 1)</td><td>14.9 ± 1.7<sup>a</sup></td><td>14.6 ± 1.7<sup>a</sup></td></tr><tr><td>Primary branch length (cm)</td><td>58.3 ± 7.4<sup>b</sup></td><td>64.0 ± 9.6<sup>a</sup></td></tr><tr><td>level 15 from the apex (Branch 2)</td><td></td><td></td></tr><tr><td>Nb of primary internode (Branch 2)</td><td>13.8 ± 1.7<sup>a</sup></td><td>14.6 ± 1.7<sup>a</sup></td></tr></tbody></table>

Measurements were performed on 300 plants from each cultivar evaluated in 2016, 4 years after their field planting (2012). The comparison was based on an analysis of variance. Values followed by different letters are significantly different at $P \le 0.05$ (Fischer test of least significance difference).

<table><caption>**TABLE 5** | Comparison of the productivity of Marsellesa® and Starmaya cultivars during three harvest campaigns.</caption><thead><tr><th rowspan="2">Harvest campaigns</th><th colspan="2">Productivity (no. of bags of 46-kg green coffee beans per hectare)</th></tr><tr><th>Marsellesa® (Sarchimor)</th><th>Starmaya F1 hybrid (CIR-SM01 × Marsellesa®)</th></tr></thead><tbody><tr><td>2015–2016</td><td>30.1<sup>b</sup></td><td>44.5<sup>a</sup></td></tr><tr><td>2016–2017</td><td>31.1<sup>b</sup></td><td>40.4<sup>a</sup></td></tr><tr><td>2017–2018</td><td>30.1<sup>b</sup></td><td>38.4<sup>a</sup></td></tr><tr><td>Average</td><td>30.4 ± 0.6<sup>b</sup></td><td>41.1 ± 3.1<sup>a</sup></td></tr></tbody></table>

The productivity (measured in the number of bags of 46-kg green coffee beans per hectare) was evaluated during 3 years. Values followed by different letters are significantly different at $P \le 0.05$ (Fischer test of least significance difference). The comparison was based on an analysis of variance based on data collected on five blocks of 1 ha each.

more than Marsellesa® under commercial growing conditions. These results confirmed the superiority of the Starmaya F1 hybrid compared with the Marsellesa® cultivar. The Starmaya F1 hybrids grown from seeds displayed similar vigor and yields compared with other F1 hybrid clones propagated by SE (Bertrand et al., 2005).

# Evaluation of Off-Types for the Tree Phenotype and Leaf Rust Resistance of Starmaya F1 Hybrids

Phenotypic observations led to the identification of mean values of 3 and 16% of off-type plants in the fields of Marsellesa® and Starmaya, respectively (Table 6). These off-type plants were characterized by tall plants, brown young leaves, and rust susceptibility. These unexpected phenotypes were probably a consequence of pollen contamination by wind and/or insects. The significant percentage of off-type plants observed in Starmaya fields could also be explained by the fact that one of the two parents (or both) was not sufficiently genetically fixed, therefore, still containing residual heterozygosity. Consequently, the cross between the two parents gave rise to a heterogeneous hybrid population instead of a homogeneous F1 progeny. Regarding leaf rust resistance, field observations identified 1 and 8% of plants susceptible to this disease in Marsellesa® and Starmaya, respectively. As Marsellesa® is a Sarchimor, coming from a cross between “Hibrido de Timor CIFC 832/2” and Villa Sarchi CIFC 971/10, it carried at least seven dominant resistance genes (SH1, SH2, SH4, SH5, SH6, SH7, and SH9) conferring a large spectrum of resistance to leaf rust (Eskes, 1989). The 1% of susceptible plants observed in the Marsellesa® field argues in favor of this cultivar not being fully fixed. On the other hand, the rust susceptibility still observed in the Starmaya F1 hybrid probably came from an illegitimate offspring resulting from alien pollination or from residual Marsellesa® heterozygosity.

Whatever the case, the fact that it was still possible to identify 16% of off-types and 8% of leaf rust-susceptible plants in the Starmaya F1 hybrid population is a disadvantage of hybrids propagated from seeds when compared with clonal F1 hybrids. Although off-types exist in F1 hybrid clones propagated by SE and are the result of somaclonal variation (aneuploidy), we have shown that a series of measures in the laboratory and nursery drastically limits their incidence of occurrence in the field (<1%) (Bobadilla Landey et al., 2013).

# Bean Quality and Cup Quality

The green bean size for Marsellesa® and Starmaya was also investigated (Table 7). This analysis revealed that bean size distribution increased significantly in Starmaya compared with Marsellesa®. For example, 25% of Starmaya green beans were retained by sieve 19, while the percentage was only 2.5 and 5.8 for Marsellesa® and Caturra red, respectively. In addition, 70% of the total green beans from Caturra red and Marsellesa® were retained by sieves 16 to 18, while the same percentage was retained by sieves 17 to 19 for Starmaya. With 98% of green exportable coffee, the gain for Starmaya in terms of exportable green beans was thus
<table><caption>**TABLE 6** | Evaluation of off-type frequencies based on plant phenotype and leaf rust sporulation for Marsellesa® and Starmaya.</caption><thead><tr><th rowspan="2">Plant phenotype</th><th colspan="2">Without rust sporulation</th><th colspan="2">With rust sporulation</th></tr><tr><th>Marsellesa®</th><th>Starmaya</th><th>Marsellesa®</th><th>Starmaya</th></tr></thead><tbody><tr><td>Dwarf and green young leaves (normal)</td><td>294</td><td>254</td><td>2</td><td>23</td></tr><tr><td>Dwarf and brown young leaves (off-type)</td><td>6</td><td>25</td><td>0</td><td>0</td></tr><tr><td>Tall and green young leaves (off-type)</td><td>4</td><td>17</td><td>1</td><td>1</td></tr><tr><td>Tall and brown young leaves (off-type)</td><td>0</td><td>4</td><td>0</td><td>0</td></tr><tr><td>Total number of plants</td><td>300</td><td>300</td><td>3</td><td>24</td></tr><tr><td>Off-type frequency (%)</td><td>3</td><td>16</td><td>1</td><td>8</td></tr></tbody></table>

Evaluations were performed in 2016, on 300 plants from each cultivar 4 years after field planting (2012). Rust resistance incidence was assessed during the rainy season. Expected phenotypes correspond to normal plants (dwarf with green young leaves). Other phenotypes are considered as off-types.

<table><caption>**TABLE 7** | Size distribution (expressed in percentage) of green coffee beans of Catura red, Marsellesa®, and Starmaya varieties of C. Arabica (2015–2016 harvest).</caption><thead><tr><th rowspan="2">Sieve number</th><th colspan="3">Bean size distribution by size (%)</th></tr><tr><th>Caturra red (control)</th><th>Marsellesa®</th><th>Starmaya F1 hybrid</th></tr></thead><tbody><tr><td>n° 20 (big)</td><td>1.05</td><td>0.48</td><td>10.62</td></tr><tr><td>n° 19</td><td>5.74</td><td>2.48</td><td>25.05</td></tr><tr><td>n° 18</td><td>21.89</td><td>15.80</td><td>30.65</td></tr><tr><td>n° 17</td><td>33.86</td><td>38.45</td><td>16.91</td></tr><tr><td>n° 16 (medium)</td><td>14.46</td><td>22.91</td><td>8.02</td></tr><tr><td>n° 15</td><td>9.83</td><td>13.37</td><td>5.22</td></tr><tr><td>n° 14</td><td>3.80</td><td>4.62</td><td>1.11</td></tr><tr><td>n° 13 (small)</td><td>1.43</td><td>1.02</td><td>0.91</td></tr><tr><td>Defects</td><td>7.94</td><td>0.87</td><td>1.51</td></tr><tr><td>Total</td><td>100.00</td><td>100.00</td><td>100.00</td></tr><tr><td>Exportable (%)*</td><td>90.63</td><td>98.11</td><td>97.58</td></tr></tbody></table>

<table><caption>**TABLE 8** | Cup sensorial attributes of coffee beans from Caturra red, Marsellesa®, and Starmaya F1 hybrid of C. arabica.</caption><thead><tr><th>Sensorial attributes</th><th>Caturra red (control)</th><th>Marsellesa®</th><th>Starmaya F1 hybrid</th></tr></thead><tbody><tr><td>Quality</td><td>SHB</td><td>SHB+</td><td>SHB+</td></tr><tr><td>Aroma</td><td>7.63</td><td>7.75</td><td>7.63</td></tr><tr><td>Flavor</td><td>7.63</td><td>7.88</td><td>7.63</td></tr><tr><td>After taste</td><td>7.5</td><td>7.75</td><td>7.63</td></tr><tr><td>Acidity</td><td>7.5</td><td>8</td><td>7.75</td></tr><tr><td>Body</td><td>7.63</td><td>7.63</td><td>7.75</td></tr><tr><td>Uniformity</td><td>10</td><td>10</td><td>10</td></tr><tr><td>Clean cup</td><td>10</td><td>10</td><td>10</td></tr><tr><td>Balance</td><td>7.63</td><td>7.75</td><td>7.63</td></tr><tr><td>Sweetness</td><td>8.88</td><td>8.75</td><td>8.88</td></tr><tr><td>General impression</td><td>7.63</td><td>7.63</td><td>7.63</td></tr><tr><td>Final merite note</td><td>82</td><td>83.13</td><td>82.5</td></tr></tbody></table>

Sensorial attributes were determined for the Caturra red cultivar corresponding to the control, for the parent Marsellesa® and for the Starmaya F1 hybrid of C. arabica, SHB, strictly hard bean.

seven points better than Caturra red but quite similar to that of Marsellesa®. In the world, there is considerable economic interest in producing the highest possible proportion of large size exportable beans (Vaast et al., 2006).

In terms of cup quality, the results presented in **Table 8** show that the highest merit for aroma was obtained by Marsellesa® beans. Beverages prepared with Starmaya and Marsellesa® displayed the same slight fragrance of caramel, whereas Caturra red was characterized by a chocolate note. Marsellesa® obtained better scores for acidity than the other Arabica samples. The same trend was observed for the overall flavor and aftertaste scores. The Starmaya beverage also had more body than Catura red and Marsellesa® and more sweetness than Marsellesa®. Regarding the merit score, Marsellesa® obtained the highest value (83.13), corresponding to very good scores for coffee beans usually produced at this altitude in Nicaragua (i.e., 950 m above sea level), while Starmaya and Caturra red obtained 82.5 and 82.0, respectively.

## SSR Analysis of the Hybrid and Its Progenitors

In order to understand the frequency of off-types observed in Marsellesa® and Starmaya, a genotypic analysis was performed with SSR markers developed in earlier studies. Ten plants of Marsellesa®, of the sterile male CIR-SM01, and of the resulting F1

hybrid Starmaya were analyzed by SSR (Table 9). For Marsellesa®, all the tested plants showed the same allelic patterns for marker numbers Sat 11, Sat 225, Sat 235, Sat 24, and Sat 254, which scored with similar alleles. Two alleles were observed for markers Sat 29, Sat 32, Sat 47, and Sat 207. Taking into account that C. arabica is

<table><caption>**TABLE 9** | Simple sequence repeat allelic patterns observed for the two parents (Marsellesa® and sterile male CIR-SM01) of the F1 hybrid Starmaya.</caption><thead><tr><th rowspan="2">Primer name</th><th colspan="2">Simple sequence repeat allelic patterns</th></tr><tr><th>Marsellesa®</th><th>CIR-SM01</th></tr></thead><tbody><tr><td>Sat 11</td><td>A1A1 (143)</td><td>A2A2 (145)</td></tr><tr><td>Sat 225</td><td>B1B1 (298)</td><td>B2B3 (269/296)</td></tr><tr><td>Sat 24</td><td>C1C1 (161)</td><td>C2C3 (153/172)</td></tr><tr><td>Sat 32</td><td>D1D2 (105/129)</td><td>D1D3D4 (105/121/133)</td></tr><tr><td>Sat 47</td><td>E1E2 (117/150)</td><td>E1E3 (117/137)</td></tr><tr><td>Sat 235</td><td>F1F1 (227)</td><td>F2F3 (259)</td></tr><tr><td>Sat 254</td><td>G1G1 (202)</td><td>G2G2 (220)</td></tr><tr><td>Sat 29</td><td>H1H2 (119/133)</td><td>H1H3 (119/135)</td></tr><tr><td>Sat 244</td><td>I1I2 (278/306)</td><td>I1I3 (278/302)</td></tr><tr><td>Sat 207</td><td>J1J2 (83/89)</td><td>J1J3 (83/93)</td></tr></tbody></table>

This analysis was performed using 10 plants for each genotype. Values indicated in italics correspond to the length (in base pairs) of fragments amplified by polymerase chain reaction.

<footnote>(*) Exportable = Total (100) – (defect + Sieve n°13).</footnote>
an allotetraploid (2n = 4x = 44) arising from natural hybridization between *C. canephora* and *C. eugenioioides* (Lashermes et al., 1999), the observation of two bands was explained by the presence of two different loci in each sub-genome and two identical alleles for each locus for these markers. Under our study conditions, all the Marsellesa® plants had the same allelic pattern, which tended to prove its high homozygosity. However, the presence of residual Marsellesa® heterozygosity could not be ruled out, due to the limited number of SSR markers used (only 10 SSR markers).

As regard the sterile male CIR-SM01, and for all SSR markers tested, all the plants presented the same patterns. This result was expected because all these plants were clonally propagated by SE. Markers Sat 11 and Sat 254 scored with similar alleles (homozygosity), while clearly residual heterozygosity was observed for marker loci Sat 32, Sat 24, and Sat 225, which identified the presence of three different alleles. We concluded that CIR-SM01 displayed constitutive heterozygosity that cannot be modified even when propagating the clone by SE.

The Starmaya F1 hybrid was assumed to be heterozygous because it received one set of chromosomes transmitted by the pollen of Marsellesa® and the ovule of CIR-SM01. In fact, different patterns were observed for Starmaya (Table 10). For example, two different patterns were detected with markers Sat 24 and Sat 32 for plants P8 and P9 (C1C3/C1C2C3 and D1D2D4/D1D2D3, respectively) confirming that these markers were heterozygous in the sterile male parent.

For the Starmaya P1 plant, two different alleles from both parents (alleles C4 and D5, respectively) were observed with marker loci Sat 24 and Sat 32. This result clearly indicates that this plant came from a cross-pollination involving alien pollen (i.e., not from Marsellesa®).

On the other hand, no Starmaya F1 hybrid plants had SSR patterns identical to the sterile male CIR-SM01, indicating that “sterile male × sterile male” crosses did not occur, thereby confirming the strict sterile male character of this accession. To reduce the residual heterozygosity carried by the pollen of Marsellesa®, we decided in 2017 to create new seed gardens using cultivars in generation F10. On the other hand, to avoid cross-pollination, we used windbreak hedges around the seed

garden plot and protected the plants with an insect-proof net. Our current objective is to produce around one ton of Starmaya seeds in Nicaragua and Costa Rica in 2019, with a maximum of 5% off-types, which will be commercially acceptable.

## CONCLUSION

We demonstrated that Starmaya F1 hybrids had better agronomic characteristics in terms of vigor, bean size, and yield than the Marsellesa® cultivar used as a parent or the Caturra red control. The Starmaya F1 hybrids also had a good cup quality compared with other traditional cultivars, such as Marsellesa® and Caturra red. The results presented here are a proof of concept that it is possible to propagate a *C. arabica* F1 hybrid on a large scale from seeds, using the CIR-SM01 sterile male parent, by setting up seed gardens. The Starmaya cultivar is already included in the global catalog of the World Coffee Research foundation (https://varieties.worldcoffeeresearch.org/). The first commercial seeds (one ton) will be available in 2019. We believe that F1 hybrids produced using male sterility can produce low-cost seeds that will be readily available wherever Arabica is grown and that this innovation will reshape the coffee seed sector. Starmaya was recently recognized by the whole coffee industry as a promising cultivar representing a major breakthrough in the widespread production of F1 hybrids for farmers around the world (https://dailycoffeenews.com/2017/03/13/new-starmaya-hybrid-could- reshape-the-industry-says-world-coffee-research/). Under the BREEDCAFS H2020 European project (http://www.breedcafs.eu/), Starmaya has been introduced in Vietnam and Cameroon, where it is being compared with local varieties. On the other hand, new hybrids with the CIR-SM01 sterile male parent have been created and are under study.

The number of F1 hybrid varieties that can be created by seed is limited by the number of sterile male plants that can be found in coffee collections of wild trees collected in the species' center of origin in Ethiopia. A consortium (World Coffee Research-CIRAD-Nestlé) is trying to identify male-sterility markers to

<table><caption>TABLE 10 | Allelic distribution patterns observed in 10 plants of Starmaya F1 hybrid. The 10 plants of Starmaya are numbered P1 to P10.</caption><thead><tr><th rowspan="2">Simple sequence repeat primer</th><th colspan="10">Allelic distribution patterns of the Starmaya F1hybrid</th></tr><tr><th>P1</th><th>P2</th><th>P3</th><th>P4</th><th>P5</th><th>P6</th><th>P7</th><th>P8</th><th>P9</th><th>P10</th></tr></thead><tbody><tr><td>Sat 11</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td><td>A1A2</td></tr><tr><td>Sat 225</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2b3</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2B3</td><td>B1B2B3</td></tr><tr><td>Sat 24</td><td>C2C3C4*</td><td>C1C3</td><td>C1C3</td><td>C1C3</td><td>C1C3</td><td>C1C3</td><td>C1C3</td><td>C1C3</td><td>C1C2C3</td><td>C1C3</td></tr><tr><td>Sat 32</td><td>D1D3D5**</td><td>D1D2D3</td><td>D1D2D3</td><td>D1D2D4</td><td>D1D2D3</td><td>D1D2D3</td><td>D1D2D3</td><td>D1D2D3</td><td>D1D2D4</td><td>D1D2D3</td></tr><tr><td>Sat 47</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td><td>E1E2E3</td></tr><tr><td>Sat 235</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td><td>F1F2</td></tr><tr><td>Sat 254</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td><td>G1G2</td></tr><tr><td>Sat 29</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td><td>H1H2H3</td></tr><tr><td>Sat 244</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td><td>I1I2I3</td></tr><tr><td>Sat 207</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td><td>J1J2J3</td></tr></tbody></table>

<footnote>*C4 allele with a fragment length of 163 bp. **C5 allele with a fragment length of 111 bp.</footnote>

<footnote>The simple sequence repeat markers correspond to those described in **Table 2**. The results obtained for the 10 plants of F1 hybrid Starmaya were compared with those previously obtained for the two parents (Marsellesa® and Sterile Male CIR-SM01) and described in **Table 9**. The black codes identified in P1 plant indicate a pollen contamination. The gray codes identified in plants P8 and P9 indicate the presence of heterozygous alleles in the CIR-SM01 parent.</footnote>
create conventional sterile male cultivars. This would definitively
open up the universe of known Arabica cultivars for use in
breeding new F1 hybrids.

AUTHOR CONTRIBUTIONS

FG: corresponding author and the writer of the paper. LM: PHD student who statistically analyzed the results of the multilocal Starmaya trials in Nicaragua. EA: director of ECOM company in Nicaragua and facilitated the setup in the field of the experiment. PC: engineer in charge of producing planting material for the experiment in ECOM company (Nicaragua). MB: director of R&D in Nicafrance Fondation. JH: technician in charge of the harvest of the parcels in Nicafrance Fondation. PM: helped to revise the document. J-CB: helped to revise the document and to

interpret the SSR results. ED: produced the MS plants by Somatic
embryogenesis in the lab for the experiment. CP: general director
of Nicafrance Fondation who provided the field in Cumplida
farm to establish the Starmaya garden seed. HE: director of
Coffee adapt team and cowriter of this paper. BB: breeder who
created the Starmaya Hybrid and cowriter of this paper.

ACKNOWLEDGMENTS

We thank Ilycaffé for cup tasting and Exportadora Atlantic S.A (ECOM Company) for its technical support in the field in Nicaragua. This research was partially funded by the BREEDCAFS project, which was supported by the European Commission under the Horizon 2020—Research and Innovation Programme, H2020-SFS-2016-2, grant agreement number: 727934.

REFERENCES

Anthony, F., Combes, M. C., Astorga, C., Bertrand, B., Graziosi, G., and
Lashermes, P. (2002). The origin of cultivated *Coffea arabica* L. varieties
revealed by AFLP and SSR markers. *Theor. Appl. Genet.* 104 (5), 894–900. doi:
10.1007/s00122-001-0798-8

Bellachew, B. (1997). "Arabica coffee breeding in Ethiopia" in *Proceedings of the 17th Colloquium of International Coffee Science Association*, ed. ASIC (Nairobi, Kenya), July 20-25, 1997, 406-414.

Bertrand, B., Aguilar, G., Santacreo, R., Anthony, F., Etienne, H., Eskes, A. B. et al. (1997). "Comportement d'hybrides F1 de Coffea arabica pour la vigueur, la production et la fertilité en Amérique Centrale" in *Proceedings of the 17th colloquium of international coffee science association*, ed. ASIC (Nairobi, Kenya), July 20-25, 1997, 415-423.

Bertrand, B., Etienne, H., Cilas, C., Charrier, A., and Baradat, P. (2005). *Coffea arabica* hybrid performance for yield, fertility and bean weight: a review. *Euphytica* 141, 255–262. doi: 10.1007/s10681-005-7681-7

Bertrand, B., Alpizar, E., Lara, L., Santa Creo, R., Hidalgo, M., Quijano, J. M., et al. (2011). Performance of *Coffea arabica* F1 hybrids in agroforestry and full-sun cropping systems in comparison with American pure line cultivars: a review. *Euphytica* 181, 1471–1458. doi: 10.1007/s10681-011-0372-7

Bobadilla Landey, R., Cenci, A., Georget, F., Bertrand, B., Camayo, G., Dechamp, E., et al. (2013). High genetic and epigenetic stability in *Coffea Arabica* plants derived from embryogenic suspensions and secondary embryogenesis as revealed by AFLP, MSAP and the phenotypic variation rate. *PLoS One* 8, e56372. doi: 10.1371/journal.pone.0056372

Bose, L. K., Jambhulkar, N. N., and Singh, O. N. (2014). Additive main effects and multiplicative interaction (AMMI) analysis of grain yield stability in early duration rice. *J. Anim. Plant Sci.* 24 (6), 1885–1897.

Chang, Z., Chen, Z., Wang, N., Xie, G., Lu, J., Yan, W., et al. (2016). Construction of a male sterility system for hybrid rice breeding and seed production using a nuclear male sterility gene. *Proc. Nat. Acad. Sci.* 113 (49), 4145–4150. doi: 10.1073/pnas.1613792113

Darai, R., Sarker, A., Sah, R. P., Pokhrel, K., and Chaudhary, R. (2017). AMMI biplot analysis for genotype × environment interaction on yield trait of high Fe content lentil genotypes in Terai and Mid-Hill environment of Nepal. *Ann. Agric. Crop Sci.* 2 (1026), 1028–1032. doi: 10.26420/annagriccropsci.2017.1026

Dufour, M., Anthony, F., Bertrand, B., and Eskes, A. (1997). Identification de caféiers mâle-stériles de "Coffea Arabica" au CATIE, Costa Rica. *Plantations, Recherche, Développement* 4 (6), 4014–4007.

Eskes, A. B. (1989). Resistance. *Resistance*, 171–291. doi: 10.1201/9781351070775-6

Etienne, H., Breton, D., Breitler, J.-C., Bertrand, B., Déchamp, E., Awada, R., et al. (2018). Coffee somatic embryogenesis: how did research, experience gained and innovations promote the commercial propagation of elite clones from the two cultivated species? *Front. Plant Sci.* 9 1630. doi: 10.3389/fpls.2018.01630

Etienne, H., Bertrand, B., Montagnon, C., Landey, R. B., Dechamp, E., Jourdan, I., et al. (2012). Un exemple de transfert de technologie réussi dans le domaine de la micropropagation: la multiplication de Coffea arabica par embryogenèse somatique. *Cah. Agric.* 21 (2-3), 115-124. doi: 10.1684/agr.2012.0553

Gallais, A. (2009). Hétérosité et variétés hybrides en amélioration des plantes. *Editions Quae.*, 156-162.

Georget, F., Courtel, P., Garcia, E., Hidalgo, M., Alpizar, E., Breitler, J. C. et al. (2017). Somatic embryogenesis-derived coffee plantlets can be efficiently propagated by horticultural rooted mini-cuttings: a boost for somatic embryogenesis. *Sci. Hortic.* 216, 177-185. doi: 10.1016/j.scienta.2016.12.017

Georget, F., Alpizar, E., Courtel, P., Hidalgo, J. M., Déchamp, E., Poncon, C. et al. (2015a). "Arabica F1 hybrid seed production based on genetic male sterility" in *Proceeding of the 25th colloquium of international coffee science association*, ed. ASIC (Armenia, Colombia), September 8-13, 2014, 8-12.

Georget, F., Courtel, P., Malo Garcia, E., Alpizar, E., Hidalgo, J. M., Ponçon, C., et al. (2015b). "A booster for commercial propagation of C. arabica F1 hybrids: somatic embryo-derived plantlets can be efficiently propagated in nursery via rooted cuttings" in *Proceeding of the 25th colloquium of international coffee science association*, ed. ASIC (Armenia, Colombia), September 8-13, 2014.

Gichuru, E. K., Agwanda, C. O., Combes, M. C., Mutitu, E. W., Ngugi, E. C. K., Bertrand, B., et al. (2008). Identification of molecular markers linked to a gene conferring resistance to coffee berry disease (Colletotrichum kahawae) in *Coffea arabica*. *Plant Pathol.* 57 (6), 1117-1124. doi: 10.1111/j.1365-3059.2008.01846.x

Havey, M. J. (2004). "The use of cytoplasmic male sterility for hybrid seed production" in *Molecular biology and biotechnology of plant organelles* (Dordrecht: Springer), 623–634. doi: 10.1007/978-1-4020-3166-3_23

Lashermes, P., Combes, M. C., Robert, J., Trouslot, P., D'Hont, A., Anthony, F., et al. (1999). Molecular characterization and origin of the *Coffea arabica* L. genome. *Mol. Gen. Genet.* 261 (2), 259–266. doi: 10.1007/s004380050965

Magesa, J. M., Mushi, I. K., Shayo, G., NG'Homa, N. M., Mdemu, S., Tarimo, E. et al. (2013). "Strengthening producer organizations to speed up the multiplication of improved hybrid coffee varieties in Tanzania" in *Proceeding of the 24th colloquium of international coffee science association*, ed. ASIC (San José, Costa Rica), November 11-16, 2012, 711-720.

Opile, W. R., and Agwanda, C. O. (1993). Propagation and distribution of cultivar Ruiru 11. *Kenya Coffee* 58, 1496-1508.

Perez-Prat, E., and Van Lookeren Campagne, M. M. (2002). Hybrid seed production and the challenge of propagating male-sterile plants. *Trends Plant Sci.* 7 (5), 199-203. doi: 10.1016/S1360-1385(02)02252-5

Purchase, J. L., Hatting, H., and Van Deventer, C. S. (2000). Genotypex environment interaction of winter wheat (Triticum aestivum L.) in South Africa: II. Stability analysis of yield performance. *S. Afr. J. Plant Soil* 17 (3), 101-107. doi: 10.1080/02571862.2000.10634878
Teri, J. M., Kilambo, D. L., Mtenga, D. J., Nyange, N. E., Nzallawahe, T. S., Chipungahelo, G. S. et al. (2005). "Improved Arabica varieties for the benefit of Tanzania coffee growers" in *Proceeding of the 20th colloquium of international coffee science association*, ed. ASIC (Bangalore, India), October 11-17, 2004, 1187-1191.

Van der Vossen, H., Bertrand, B., and Charrier, A. (2015). Next generation variety development for sustainable production of arabica coffee (*Coffea arabica* L.): a review. *Euphytica* 204, 243–256. doi: 10.1007/s10681-015-1398-z

Van der Vossen, H. (1985). "Coffee selection and breeding," in *Coffee* (Boston: Springer), 48–96. doi: 10.1007/978-1-4615-6657-1_3

Vaast, P., Bertrand, B., Perriot, J. J., Guyot, B., and Genard, M. (2006). Fruit thinning and shade improve bean characteristics and beverage quality of coffee (*Coffea arabica* L.) under optimal conditions. *J. Sci. Food Agric.* 86, 197–204. doi: 10.1002/jsfa.2338

Walyaro, D. J. A. (1983). *Considerations in breeding for improved yield and quality in arabica coffee (Coffea arabica L.). Doctoral dissertation*. Wageningen University & Research: Agricultural University Wageningen.

Wu, Y., Fox, T. W., Trimnell, M. R., Wang, L., Xu, R. J., Cigan, A. M., et al. (2016). Development of a novel recessive genetic male sterility system for hybrid seed production in maize and other cross-pollinating crops. *Plant Biotechnol. J.* 14, 1046–1054. doi: 10.1111/pbi.12477

Zhang, Q., Shen, B. Z., Dai, X. K., Mei, M. H., Maroof, M. S., and Li, Z. B. (1994). Using bulked extremes and recessive class to map genes for photoperiod-sensitive genic male sterility in rice. *Proc. Nat. Acad. Sci.* 91, 8675–8679. doi: 10.1073/pnas.91.18.8675

Zobel, J. (1998). "How reliable are the results of large-scale information retrieval experiments?" in *Proceedings of the 21st annual international ACM SIGIR conference on Research and development in information retrieval* (Melbourne, Australia), August 24-28, 1998, 307–314. doi: 10.1145/290941.291014

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Georget, Marie, Alpizar, Courtel, Bordeaux, Hidalgo, Marraccini, Breitler, Déchamp, Poncon, Etienne and Bertrand. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.