# An innovative automated active compound screening system allows high-throughput optimization of somatic embryogenesis in Coffea arabica

Rayan Awada, Dorothée Verdier, Solène Froger, Eric Brulard, Simone de Faria

Maraschin, Hervé Etienne, David Breton

► To cite this version:

Rayan Awada, Dorothée Verdier, Solène Froger, Eric Brulard, Simone de Faria Maraschin, et al. An innovative automated active compound screening system allows high-throughput optimization of somatic embryogenesis in Coffea arabica. Scientific Reports, 2020, 10 (1), pp.810. 10.1038/s41598-020-57800-6. http://dx.doi.org/10.1038/s41598-020-57800-6

HAL Id: hal-03454869

https://hal.umontpellier.fr/hal-03454869

Submitted on 29 Nov 2021

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
OPEN

# An innovative automated active compound screening system allows high-throughput optimization of somatic embryogenesis in *Coffea arabica*

Rayan Awada<fnref n="1" />,²,³*, Dorothée Verdier<fnref n="1" />, Solène Froger<fnref n="1" />, Eric Brulard<fnref n="1" />, Simone de Faria Maraschin⁴, Hervé Etienne²,³ & David Breton<fnref n="1" />

Somatic embryogenesis (SE) faces many challenges in fulfilling the growing demand for elite materials. A high-throughput approach is required to accelerate the optimization of SE protocols by multiplying experimental conditions within a limited time period. For the first time in plant micropropagation, we have developed a miniaturized and automated screening system to meet high-throughput standards. *Coffea arabica* embryo regeneration, classically achieved in 250-ml Erlenmeyer flasks, was successfully miniaturized in 24-well plates, allowing a volume downscaling factor of 100 and a space saving of 53 cm²/well. Cell clusters were ground and filtered to fit the automated pipetting platform, leading to fast, reproducible and uniform cluster distribution (23.0 ± 5.5 cell clusters/well) and successful regeneration (6.5 ± 2.2 embryos/well). Pilot screening of active compounds on SE was carried out. Compounds belonging to the histone deacetylase inhibitor family were tested for embryo regeneration efficiency. Cells treated with 1 µM Trichostatin A showed a marked 3-fold increase in the number of regenerated embryos. When re-tested in 250-ml flasks, the same enhancement was obtained, thereby validating the miniaturized and automated screening method. These results showed that our screening system is reliable and well suited to screening hundreds of compounds, offering unprecedented perspectives in plant micropropagation.

The rapid increase in environmental stresses, especially those associated with climate change, and growing demand from the food and timber industries, have put enormous pressure on trees<fnref n="1" />. Developing new technologies for the clonal propagation, improvement and breeding of trees can help solve these problems². This has been achieved in part using biotechnology methodologies, such as *in vitro* propagation, genetic transformation, and marker-assisted breeding for the gradual genetic improvement of woody plants³,⁴.

For *in vitro* propagation, shoot proliferation methods using adventitious shoots and axillary buds have been widely applied, though research is still largely focused on using somatic embryogenesis (SE) in woody plants⁵. Today, SE is considered as a very promising tool for the rapid and large-scale propagation of elite varieties in many species⁵,⁶, and is particularly important for those that have a long life cycle and are difficult to propagate by conventional methods⁷.

However, to date, very few examples of commercial SE applications are available. The industrial application of SE in hardwoods remains limited to a few species, such as the hybrid yellow-poplar, cocoa, and coffee⁸. SE technology is still work-in-progress, mainly due to low biological efficiency and high genotypic effects resulting in high production costs⁶,⁹.

To guarantee the security of the supply chain in response to market demand for elite materials, it is necessary to maximize SE-derived plant production and quality, i.e. to improve SE protocols. Whatever the plant species, a number of common bottlenecks have been noted throughout the main three steps of the current SE process –(1)

<footnote n="1">¹Nestlé R&D Center Tours, 101 avenue Gustave Eiffel, 37097 Cedex 2, Tours, France. ²CIRAD, UMR IPME, F-34398, Montpellier, France. ³UMR IPME, Université de Montpellier, CIRAD, IRD, F-34398, Montpellier, France. ⁴Bioalternatives SAS, 1 bis Rue des Plantes, 86160, Gençay, France. *email: rayan.awada@cirad.fr</footnote>
the induction of embryogenic callus, (2) its proliferation, and (3) its differentiation into embryos. For instance, a strong genotypic effect has often been reported in the first step¹⁰. Explants of some interesting clones either do not respond to auxin stimulus, or respond by directly producing embryos without passing through a proliferating callus step. For the second step, calli of some clones lose their embryogenic potential during proliferation¹¹. For the third step, maximizing embryo yields is a necessity to secure the supply chain and reduce production costs¹².

Until now, one of the conventional methods for improving SE protocols has consisted in applying active compounds. Adding small (<500 Da) bioactive compounds (nutrients, growth regulators, epigenetic regulators, etc.) to the culture medium that have a potential effect on SE targets, and are capable of modulating the physiological state of cells, is a possible way of maximizing embryo production and quality, and of removing genetic or stress-induced recalcitrance¹³. The ability to use small compounds to improve tissue culture responses eliminates the need to create and market transgenic plants¹⁴. This is important in the food sector, where consumers are reluctant to consume transgenic products.

However, testing the effect of new compounds on SE efficiency under standard conditions is generally a laborous, intensive and tricky operation, because it is time-consuming, space-demanding and costly. To date, the literature on SE in several species has shown that only a limited number of compounds have been chosen to be tested, on the basis of prior hypotheses (targeted selection)¹⁵⁻¹⁷ and when tested, only a very limited range of concentrations was chosen. The above inconveniences also hampered the possibility of testing these compounds in combinations. The resulting low-throughput research strategy has led to overall slow technical progress on SE protocols in recent years.

At the same time, high-throughput strategies are gaining increasing attention, particularly in the fields of biofuel production and pharmaceuticals¹⁸,¹⁹. These strategies consist in testing a large number of compounds at different concentrations (untargeted selection) in a miniaturized and automated system, to systematically screen for novel modifiers of a biological phenomenon. They combine large-scale chemistry and biology data, along with bioinformatics, which is required for data mining²⁰. The effectiveness of these strategies is enhanced by the fact that most active compounds are small molecules that modulate target proteins and/or pathways of a determine biological process²¹. In the last decade, several collections of bioactive compounds became available for the research community, enabling the screening of hundreds of libraries for novel activities²². Nevertheless, to our knowledge, this kind of approach has never been applied to optimizing culture media for micropropagation, particularly for SE.

In plants and mammals, epigenetic regulation and histone modifications occur widely during cellular differentiation and development²³,²⁴. Nic-Can et al.²⁵ showed that genes related to SE, such as *LEAFY COTYLEDON1*, *BABY BOOM1* and *WUSCHEL-RELATED HOMEOBOX 4*, are epigenetically regulated in *Coffea canephora*. Boutilier et al.¹⁴ showed that the switch to haploid embryogenesis is controlled by the activity of histone deacetylases (HDAC). Blocking HDAC activity with HDAC inhibitors (HDACi) in *Brassica napus*, *Brassica rapa*, *Arabidopsis thaliana* and *Capsicum annuum* male gametophytes resulted in a large increase in the proportion of cells undergoing embryogenic growth¹⁴.

Coffee (*C. arabica* & *C. canephora* cv Robusta) has a great economic impact in many producing countries, especially in South America²⁶. Today, coffee SE is one of the most advanced technologies. Research on SE by two leading groups - the CIRAD/ECOM alliance and Nestlé - has led over the last 13 to 15 years to the successful industrialization and commercial application of SE for the two cultivated species, through the large-scale dissemination of Arabica F1 hybrids and Robusta clones⁸. Even though the processes established in the 2 species display good biological efficiency, with low genotypic effects²⁷,²⁸ and controlled somaclonal variation¹¹,²⁹, SE production costs remain high and still cannot compete with the production costs of conventional seedlings or cuttings³⁰.

While the micropropagation of *coffee* *in vitro* plantlets can today reach 2–3 M SE-derived plants per year, SE faces many challenges in fulfilling growing market demand for elite coffee varieties, estimated at around 50–100 M plants per year⁸. This market demand for more resilient varieties comes from the recent loss of leaf rust resistance in all the conventional cultivated line varieties for Arabica, and the increasing impact of the harmful effects of climate change. For Robusta, there is also an urgent need to replace unselected clones with improved ones, particularly in Vietnam – the world's main coffee producer – where aging coffee orchards need to be totally renewed⁸.

As for other woody species, optimizing coffee SE protocols is a long and laborious process. Yet, increasing production efficiency *in vitro* would make SE technology more cost-effective, hence the interest in finding new active compounds able to increase embryogenic callus formation and embryo yield. A high-throughput approach might be a solution, given the possibility of multiplying experimental conditions in a limited time period.

We focus in this paper on one of the main three SE steps in *C. arabica*: embryo regeneration from embryogenic cell clusters. This SE step is classically achieved in 250-ml Erlenmeyer flasks. Here, we describe how we miniaturized this step in 24-well plates and made it compatible with an automated pipetting platform. Cell cluster miniaturization was also required. Once miniaturization and automation were successful, we carried out pilot screening using compounds belonging to the HDACi family, for their potential role as active compounds to increase embryo yield. Lastly, treatments showing a significant positive effect on embryo yield were tested again under standard *in vitro* conditions to confirm treatment efficiency and validate the automated and miniaturized screening method.

## Results

### Grinding and filtration of embryogenic tissues to overcome the pipetting constraint.

Preliminary profiling of long-term cultured Arabica embryogenic cell clusters showed that their diameter fell within the range of 100–1000 µm and followed a statistical normal distribution, with a maximum peak around 500–600 µm (Fig. 1a). For automation purposes, uniformity and reproducibility are required. Hence, in order for the cell clusters to fit the pipetting tips, filtration at 700 µm was needed to exclude a diameter greater than 700 µm,
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Successful miniaturization of Arabica embryogenic cell clusters. (a) Profiling of cell clusters according to their diameter before grinding and filtration. Data of 210 cell clusters were plotted. The black curves show that the plotted values are normally distributed. (b) Profiling of cell clusters according to their diameter after grinding and filtration at 700 µm. Data of 1210 cell clusters were plotted. Viability of cell clusters before (c) and after grinding (d) following double staining with FDA and PI. FDA fluoresces in green in living cells and PI fluoresces in red in dead cells.</figcaption></figure>

which amounted to approximately 40% of the whole cell suspension. To benefit from all the available embryogenic material at our disposal, we decided to add a grinding step to reduce and homogenize cell cluster diameter. Figure 1b shows the effect of grinding on the particle diameter distribution. Cell cluster diameters maintained their statistical normal distribution, but with a shift of the maximum peak towards the range of 100–150 µm. Viability tests (Fig. 1c,d) showed similar staining profiles between ground and intact cell clusters, with a prevailing green light emission indicating that cells had high viability. Some peripheral cells emitted red light, indicating their dead status in both ground and intact clusters, though the red emission intensity was more noticeable in ground cells, due to the effect of the grinding blades.

The regeneration capacity of ground material had first to be tested under standard *in vitro* culture conditions (Erlenmeyer flasks) (Fig. 2). After 9 weeks of culturing, an average of 3,039 ± 485 normal torpedo-shaped embryos were obtained per gram of inoculated callus (embryos. g⁻¹ callus) representing a biomass of 0.43 ± 0.12 g (n = 6 flasks/condition/replicate). These values were significantly lower compared to those obtained with intact cell clusters, i.e. 4,611 ± 685 embryos. g⁻¹ callus and a biomass of 0.83 ± 0.18 g (***P < 0.001) (Fig. 2a,b). However, it should be noted that better synchronization of embryo development was observed. Embryos generated from ground cell clusters were morphologically more uniform, since very few proliferating clusters remained, limiting spatial competition (Fig. 2c,d).

**Achieving uniform cell cluster distribution after automated pipetting in miniaturized culture containers.** In order to obtain uniform distribution of cell clusters at the start of culturing, two parameters were quantified to verify uniformity of distribution: the number of cell clusters in each well, and the total area they occupied. The values of 4 independent replicates (n = 24 wells/replicate) were uniform (P > 0.05) and similarly distributed between plates (1, 2, 3, 4) and between rows (A, B, C, D), having an average of 23.0 ± 5.6 clusters/well and occupying an area of 1.39 ± 0.12 mm²/well (Table 1, Fig. 3a).

The preliminary results showed medium evaporation in peripheral wells (loss of 40% of total volume after 4 weeks compared to central wells) causing non-reproducible embryo regeneration (data not shown). Parafilm-sealed plates were then placed in a culture chamber with 90% RH, which drastically limited the evaporation rate to less than 2% in all wells, and led to uniformity between tested replicates (P > 0.05). Similar distribution of biomass parameters was obtained between peripheral wells (22.5 ± 5.39 cell clusters/well; area occupied 1.39 ± 0.13 mm²) and central wells (23.9 ± 5.63 clusters/well; area occupied 1.39 ± 0.11 mm²) (Table 1).
<figure><img src="image_3.png" /></figure>

**Figure 2.** Assessing the effect of cell cluster grinding on embryo regeneration in Erlenmeyer flasks. (a) Embryo yield (embryos. g⁻¹ callus) obtained after 9 weeks of embryo regeneration. (b) Biomass (g) obtained after 9 weeks of embryo regeneration. Bars show the mean of triplicate samples (n = 6 flasks/condition/replicate) and error bars represent the SD. Differences between both groups were analysed with t-tests (***P < 0.001). (c) Appearance of torpedo-shaped embryos (white structures) obtained after 9 weeks of embryo differentiation from intact cell clusters. Remaining embryogenic clusters can still be observed (brown structures). (d) Appearance of torpedo-shaped embryos regenerated after 9 weeks of differentiation from ground cell clusters. Better synchronization of embryo development is observed.

**Maintenance of regeneration capacity and embryo yield reproducibility in a miniaturized environment.** The efficiency of embryogenic cell cluster conversion to embryos was validated in the 24-well plate culture system. The regenerated embryos were normally shaped, with a conversion rate of 6.5 ± 2.1 torpedo-shaped embryos/well after 9 weeks (Table 1, Fig. 3b). The values of 4 independent replicates (n = 24 wells/replicate) followed a statistical normal distribution. High uniformity in embryo yield rates was found between plates, rows and peripheral and central wells (P > 0.05), confirming the very limited impact of automated pipetting on cluster distribution, whatever the well position. However, it should be noted that expressing this embryo yield in wells by gram of inoculated callus resulted in a mean of 13,000 ± 4,200 embryos, which was higher than the embryo yield obtained in a flask culture system from ground or intact cell clusters, i.e. 3,039 ± 485 embryos and 4,611 ± 685 embryos, respectively. However, for the same inoculated cell cluster density, a delay in embryo development was found in the miniaturized system, where the first globular-shaped embryos emerged after 6 weeks (Fig. 4d–f), as opposed to 4 weeks in Erlenmeyer flasks.

**Successful pilot screening for the effect of active compounds on somatic embryogenesis.** After satisfactory implementation of a miniaturized and automated embryo regeneration assay, pilot screening was carried out by testing histone deacetylase inhibitor (HDACi) treatments. Four compounds belonging to the HDACi family – Oxamflatin, SAHA, Scriptaid and Trichostatin A – were tested at 3 different concentrations (0.1, 1 and 10 µM) as potential active compounds affecting *C. arabica* embryo differentiation. The regenerated embryo numbers of 3 independent replicates (n = 18 wells/HDACi/concentration/replicate) were uniform between plates and followed a statistical normal distribution between culture wells. Only clusters treated with 1 µM TSA showed a significant increase in the number of torpedo-shaped embryos obtained after 9 weeks (18.3 ± 3.9 embryos/well vs 6.0 ± 2.5 for untreated control, ***P < 0.001) (Fig. 4a,b) and the appearance of the regenerated torpedo-shaped embryos was comparable between clusters treated or not with TSA (Fig. 4b). The area occupied by cell clusters in a well was recorded over a 6-week culture period up to the emergence of globular-shaped embryos (Fig. 4c,d). Comparing the resulting curves established with cell clusters treated with 1 µM TSA and untreated clusters showed that cluster growth started after 3 weeks and exponentially increased up to the sixth week of culturing (R² = 0.93), when the first globular-shaped embryos emerged (Fig. 4d). From 4 weeks under embryo regeneration conditions, the growth of cell clusters treated with 1 µM TSA was significantly better compared to the untreated control (total area of 35.0 ± 3.9 mm² vs 21.9 ± 3.4 mm² after 6 weeks, ***P < 0.001), explaining the larger number of regenerated embryos obtained after 9 weeks (Fig. 4a).
<table><thead><tr><th rowspan="3">Parameters</th><th colspan="10">Assessment of distribution uniformity</th></tr><tr><th colspan="4">Between plates</th><th colspan="4">Between rows</th><th colspan="2">Between central and peripheral wells</th></tr><tr><th>Plate 1</th><th>Plate 2</th><th>Plate 3</th><th>Plate 4</th><th>Row A</th><th>Row B</th><th>Row C</th><th>Row D</th><th>Central wells</th><th>Peripheral wells</th></tr></thead><tbody><tr><td>Number of cell clusters per well</td><td>22.1 ± 5.7a</td><td>24.4 ± 6.0a</td><td>22.0 ± 5.5a</td><td>23.3 ± 5.0a</td><td>21.1 ± 4.7a</td><td>24.5 ± 5.4a</td><td>22.5 ± 5.0a</td><td>23.8 ± 6.6a</td><td>23.9 ± 5.6a</td><td>22.5 ± 5.4a</td></tr><tr><td>Total area occupied (mm²)</td><td>1.40 ± 0.10a</td><td>1.41 ± 0.12a</td><td>1.37 ± 0.13a</td><td>1.38 ± 0.12a</td><td>1.39 ± 0.11a</td><td>1.38 ± 0.10a</td><td>1.41 ± 0.12a</td><td>1.38 ± 0.15a</td><td>1.39 ± 0.11a</td><td>1.39 ± 0.13a</td></tr><tr><td>Embryo yield per well</td><td>6.9 ± 2.5a</td><td>6.7 ± 2.4a</td><td>6.2 ± 1.7a</td><td>6.3 ± 2.2a</td><td>7.0 ± 2.9a</td><td>6.5 ± 1.8a</td><td>6.2 ± 2.2a</td><td>6.4 ± 1.8a</td><td>6.3 ± 2.0a</td><td>6.6 ± 2.3a</td></tr></tbody></table>

**Table 1.** Assessing the uniformity and reproducibility of automated distribution of Arabica cell clusters in 24-well plates. Uniformity of distribution was revealed between plates, between rows, and between central and peripheral wells. Three parameters were quantified: the number of clusters per well and the total area they occupied (mm²) at the beginning of the regeneration step, and the number of regenerated embryos per well after 9 weeks. Values represent the mean ± SD of 4 independent replicates (n = 24 wells/replicate). For each parameter, differences between plates and between rows were analysed with a one-way ANOVA test, while a t-test was used to analyse differences between central and peripheral wells. The letter “a” shows no significant difference between groups (P > 0.05).

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

**Figure 3.** Successful uniform distribution of Arabica cell clusters resulting in a uniform number of regenerated embryos. Programming the automated pipetting platform resulted in uniform distribution of cell clusters between wells at the beginning of the regeneration step as well as a uniform number of embryos regenerated after 9 weeks. (a) Appearance of inoculated cell clusters in a 24-well plate after distribution using the automated pipetting platform after cell cluster grinding and filtration. (b) Appearance of embryos regenerated from ground clusters in a “24-well plate” miniaturized culture system.

To validate the screening carried out in the automated and miniaturized system, the best treatment obtained with 1 µM TSA, with an exposure time of 24 hours, was re-tested under standard *in vitro* conditions, i.e. 250-ml Erlenmeyer flasks. The TSA treatment enhanced embryo regeneration in the same way compared to the untreated embryogenic clusters (9,325 ± 2,240 embryos. g⁻¹ callus vs 4,673 ± 1,614 embryos. g⁻¹ callus, respectively; P < 0.001; Fig. 5a). The beneficial effects of TSA were also visible through the increased overall biomass (1.12 ± 0.26 g vs 0.69 ± 0.15 g, P < 0.001; Fig. 5b), and the increased length of the treated torpedo-shaped embryos (896 ± 463 µm vs 760 ± 429 µm, respectively, P < 0.001; Fig. 5c). Extended exposure to 1 µM TSA, which was also tested (7, 14 d), did not significantly improve the embryo regeneration yield compared to an exposure time of 24 h, but biomass and torpedo mean size significantly increased with an exposure of 14 d to 1 µM TSA.

**Treating cell clusters with TSA does not affect plant regeneration.** The phenotypic data recorded on 160 young plants ready for greenhouse transfer did not reveal any significant differences between plants regenerated from cell clusters treated, or not, with TSA. An ANOVA did not indicate any treatment effect on either main root or stem lengths (data not shown). All plants displayed orthotropic development, with uniform pigmentation and leaf shape (Fig. 6a). A Multiple Component Analysis (MCA) performed on the scores obtained for the main phenotypic categorical variables (number of pairs of leaves, number of roots, number of stems, leaf colour index, Table 2) provided 19 components reaching 3.8 inertia (data not shown). However, the first
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Successful pilot screening of active compounds affecting SE in the automated and miniaturized system. Four compounds belonging to the HDACi family (Oxamflatin, Scriptaid, SAHA and TSA) were tested in 24-well plates at three different concentrations (0.1, 1 and 10 µM), as potential active compounds affecting *C. arabica* embryo differentiation. Cell clusters were exposed for 24 h to treatment before washing. (a) Embryo regeneration was scored after 9 weeks. The “CTRL” group was untreated (control). Bars show the mean of triplicate samples (n = 18 wells/HDACi/concentration/replicate) and error bars represent the SD. Differences between treated clusters and the control were analysed with a one-way ANOVA test (***P < 0.001). (b) Picture showing the increased embryo yield in the wells treated with 1 µM TSA (bottom) vs control wells (top). (c) The total area occupied by cell clusters was monitored during the embryo regeneration period. Plotted values are the mean of triplicate samples (n = 18 wells/treatment/replicate) and error bars represent the SD. Differences between clusters treated with 1 µM TSA and the control were analysed with a t-test (***P < 0.001). Cell cluster growth followed an exponential trendline (R² = 0.93) up to the sixth week. (d–f) Pictures of an emerged globular-shaped, heart-shaped and torpedo-shaped embryo after 6, 7 and 8 weeks under embryo regeneration conditions, respectively.</figcaption></figure>

two axes accounted for only 18.9% of total inertia. Plotting categorical data according to axes 1 and 2 revealed a strong contribution to global variability for 2 plants scored “376” and one scored “382” for their colour index (Supplementary Data). As a result, no distinction could be detected between TSA treatments so we decided to
<figure><img src="image_6.png" /></figure>

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

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

**Figure 5.** Validation in Erlenmeyer flasks of active compounds highlighted in the pilot screening. Treatment with 1 µM TSA was re-tested in standard 250-ml Erlenmeyer flasks with an exposure time varying between 1, 7 and 14 days. (a) Embryo yield (embryos. g⁻¹ callus) obtained after 9 weeks. (b) Biomass (g) obtained after 9 weeks. (c) Embryo length. Bars show the mean of triplicate samples (n = 6 flasks/exposure time/replicate for (a,b) and n = 100 embryos/exposure time/replicate for (c)) and error bars represent SD. The “CTRL” group was untreated (control). Differences between treated clusters and the control were analysed with a one-way ANOVA test for (a,b), while a Kruskal-Wallis test was used for (c) followed by a Dunn test as a post hoc. Means with different letters are significantly different (*P* < 0.001).

discard these data and run another MCA. This time, 17 components were detected with a total level of 3.4 inertia (data not shown). The contribution of the first two axes to this inertia reached 19.6%. This time, a graphical analysis revealed 2 main groups of categorical data and a smaller third one including atypical, double stem and unrooted plants (Fig. 6e). This group did not include TSA treatment or control conditions, indicating a lack of correlation between variables. Considering the upper left of the graph, plants from untreated tissues and tissues treated for 7 days with TSA correlated with plants bearing 2 to 4 pairs of leaves, 1 to 2 roots and various leaf colour indexes (357, 364, 371). Conversely, on the right of the graph, plants arising from cell clusters treated for 1 and 14 days correlated with well-developed plants: 5 to 7 pairs of leaves, 3 to 4 roots, and a stable colour index (377). Thus, treating cell clusters with 1 µM TSA did not show any significant deleterious effect on regenerated plants.

## Discussion

We proved in this study that a major step of SE processes, i.e. embryo regeneration from embryogenic cell clusters, can be successfully miniaturized and automated allowing easy, quick, cheap and reliable testing of numerous compounds at varied concentrations. Switching from conventional 250-ml flask culture containers filled with 100 ml of culture medium to 24-well plates with 1 ml of available nutritive medium allowed a miniaturization factor of 100. With this strong downscaling, consideration could be given to testing the effectiveness of a much larger number of molecules at different concentrations, while reducing the space required (2.5 cm² required per well compared to 55.4 cm² required per flask and the possibility of stacking multi-well plates). In addition, automation led to a considerable gain in speed and reduced manpower, where only 1 working full-time employee (FTE) was needed to supervise the pipetting platform. For instance, the pilot screening carried out using HDACi treatments required 36 plates, 864 wells, i.e. 864 independent culture environments. The same experiment in flasks was simply unmanageable, as 1 FTE can only inoculate 100 flasks a day, not to mention the space required, whereas inoculation of several hundred wells can be done within hours. Many authors have also shown the importance of miniaturization and automation in animal or plant cell screening systems³¹⁻³³. For instance, Rodriguez-Furlan et al.³³ also showed that the 24-well microplate format used to screen for compounds modulating *A. thaliana* seedling growth required a minimum of sample handling and allowed automatic acquisition of quantitative data and rapid monitoring of morphological traits (leaf and root growth). This container format, along with automation, allowed the screening of 10,000 structurally diverse compounds in a limited time period.

While animal cell screening systems are mainly derived from drug discovery programmes and cancer research¹⁸,³⁴, plant cell screening systems were only recently used to study compound effects on growth and development. To date, in plants, only *A. thaliana* seedlings have been used in such studies, because of their small size,
<figure><img src="image_9.png" /><figcaption>**Figure 6.** *In vitro* plant phenotypic characterization following 1 µM TSA treatment during the coffee SE process. (a) Plants regenerated from untreated embryogenic callus. (b–d) Plants regenerated from TSA-treated embryogenic clusters for 1, 7 and 14 days, respectively. (e) Graphic representation of a Multiple Component Analysis (MCA) performed on plant phenotypic variables. Data from 157 plants were plotted according to the first 2 axes of the MCA. The number of pairs of leaves are represented by white squares, the number of roots by crosses, the number of stems by black spots, the leaf colour index by triangles and the *in vitro* treatment prior to regeneration by black squares.</figcaption></figure>

allowing development in multi-well containers³,³⁵. However, automated screening systems for active compounds have never been used in plant propagation procedures. Our system provided the first evidence in plants that miniaturized containers do not hamper embryogenic cell cluster growth and differentiation, and thus enable screening for the effectiveness of compounds in micropropagation, under real conditions and on economically important crops.

Our results showed that even after miniaturizing the embryo regeneration step, mechanisms occurring during conventional-scale SE were preserved, as reliable production of well-shaped embryos was achieved. SE miniaturization did not engender any developmental heterogeneity. However, the type of grinding chosen had significant effects by lowering embryo yields in flasks. Since this biological efficiency was reproducible during this step, we did not consider this difference to be a drawback, because the main objective of this miniaturized system was not large-scale embryo production, or to study molecular mechanisms underlying embryogenesis, but to use it on SE as a tool to screen for active compounds, with a view to speeding up protocol optimization.
<table><caption>**Table 2.** Phenotypical distribution of the 160 acclimated plants as a function of the TSA treatments *in vitro*. Each plant was individually scored for their leaf pairs (lp) number, roots (r) number, stems (a) number and colour index. Treatments consisted in growing *in vitro* embryogenic clusters in contact with 1 µM TSA for 0 (control), 1, 7 and 14 days.</caption><thead><tr><th colspan="10">Phenotypic scores</th></tr><tr><th></th><th>2lp</th><th>3lp</th><th>4lp</th><th>5lp</th><th>6lp</th><th>7lp</th><th></th><th></th><th></th></tr></thead><tbody><tr><td>Leaves</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Control</td><td></td><td>12</td><td>18</td><td>10</td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA1</td><td></td><td>13</td><td>14</td><td>11</td><td>2</td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA7</td><td></td><td>6</td><td>15</td><td>16</td><td>2</td><td>1</td><td></td><td>40</td><td></td></tr><tr><td>TSA14</td><td>1</td><td>17</td><td>16</td><td>6</td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>Total</td><td>1</td><td>48</td><td>63</td><td>43</td><td>4</td><td>1</td><td></td><td>160</td><td></td></tr><tr><td>Roots</td><td>0r</td><td>1r</td><td>2r</td><td>3r</td><td>4r</td><td></td><td></td><td>Total</td><td></td></tr><tr><td>Control</td><td>5</td><td>20</td><td>13</td><td>2</td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA1</td><td>1</td><td>20</td><td>15</td><td>3</td><td>1</td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA7</td><td>3</td><td>17</td><td>14</td><td>5</td><td>1</td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA14</td><td>3</td><td>19</td><td>14</td><td>4</td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>Total</td><td>12</td><td>76</td><td>56</td><td>14</td><td>2</td><td></td><td></td><td>160</td><td></td></tr><tr><td>Stems</td><td>1a</td><td>2a</td><td></td><td></td><td></td><td></td><td></td><td>Total</td><td></td></tr><tr><td>Control</td><td>36</td><td>4</td><td></td><td></td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA1</td><td>38</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA7</td><td>37</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>TSA14</td><td>36</td><td>4</td><td></td><td></td><td></td><td></td><td></td><td>40</td><td></td></tr><tr><td>Total</td><td>147</td><td>13</td><td></td><td></td><td></td><td></td><td></td><td>160</td><td></td></tr><tr><td>Colour index</td><td>357</td><td>364</td><td>370</td><td>371</td><td>376</td><td>377</td><td>382</td><td>Total</td><td></td></tr><tr><td>Control</td><td>1</td><td>11</td><td>13</td><td>6</td><td>1</td><td>8</td><td></td><td>40</td><td></td></tr><tr><td>TSA1</td><td></td><td>2</td><td></td><td>5</td><td></td><td>33</td><td></td><td>40</td><td></td></tr><tr><td>TSA7</td><td></td><td>3</td><td>4</td><td>4</td><td>1</td><td>27</td><td>1</td><td>40</td><td></td></tr><tr><td>TSA14</td><td></td><td>9</td><td>20</td><td>4</td><td></td><td>7</td><td></td><td>40</td><td></td></tr><tr><td>Total</td><td>1</td><td>25</td><td>37</td><td>19</td><td>2</td><td>75</td><td>1</td><td>160</td><td></td></tr></tbody></table>

A difference in embryo development timing was also noted between the Erlenmeyer flask and 24-well plate culture systems, which could be explained by obviously different physico-chemical environments. While rotary shaking was set to the same speed for both systems (120 rpm), it showed different effects, mainly with fast spatial movement of the nutritive medium in flasks, compared to slow movement in wells, due to their smaller diameter. Such lower shaking is known to affect medium O₂/CO₂ uptake³⁶,³⁷ and probably explained the differences between the 2 culture systems.

Pilot screening was carried out as proof-of-concept to show that the miniaturized automated culture system we developed was efficient for testing the effect of putatively active compounds on SE, and was reliable and reproducible. Having uniformity of cell cluster distribution is essential before initiating regeneration, to avoid false positives and negatives when screening for active compounds. The use of an automated scanning microscope enabled a considerable gain in speed and quality in terms of image acquisition and analysis. Proof-of-concept was achieved by proving the stimulation of coffee SE by TSA. A 1 µM TSA supply promoted the proliferation of embryogenic cells and led to a 3-fold increase in embryo yield. When re-tested under standard conditions, treatment with 1 µM TSA confirmed its efficiency, hence validating the whole screening system pipeline. Other authors also showed positive effects of TSA on embryogenic structure formation. Li *et al.*³⁸ showed that 0.5 µM TSA promoted totipotency in male gametophytes of *B. napus* and *A. thaliana*. Male gametophytes, which are rather recalcitrant in haploid embryo development in cultures, formed embryogenic cell clusters after TSA treatment. In agreement with these authors, we hypothesize that the 0.1 µM TSA treatment used in our study was at suboptimal concentration to have any observable effect on SE expression while a 10 µM TSA treatment was toxic to cell clusters preventing embryogenesis to occur. Li *et al.*³⁸ associated the TSA mode of action with hyper-acetylation of histones H3 and H4 to block histone deacetylase (HDAC) activity. Wójcikowska *et al.*³⁹ recently showed that treatment of *A. thaliana* explants with TSA induces SE without the exogenous application of auxin and, thus, the TSA mode of action is related to the activation of auxin-responsive transcription factor genes. In conifers, TSA caused the de-repression of *LEAFY COTYLEDON* genes and stimulated the initiation of SE from cotyledonary embryos in *Picea abies*⁴⁰ and in *Pinus Sylvestris*⁴¹, but inhibited embryo germination. For *C. arabica*, 1 µM TSA did not have any significant deleterious effect on plants regenerated from treated cell clusters, and thus can now be considered as an active compound for optimizing SE protocols. As increasing the time of exposure to 1 µM TSA did not differ in embryo yield, we put forward the hypothesis that the mechanism by which TSA triggers its response, i.e. blocking histone deacetylase (HDAC) activity, had already occurred on day 1.

This screening system proof-of-concept is a key step towards high-throughput active compound screening platforms. The next step will be to scale-up this technology by screening hundreds of libraries on robotic platforms. Such platforms are widely used in drug discovery programmes¹⁸,³⁶. They combine robotics, data processing
and control software, liquid handling devices, and sensitive detectors to quickly conduct thousands of tests⁴², allowing high-throughput screening (HTS) for active compounds able to modulate a particular molecular pathway. Current generations of platforms are so robust and user-friendly⁴³ that HTS can be used to investigate new areas in biology, such as micropropagation. For coffee SE, achieving the intended change of scale, i.e. HTS of hundreds of libraries to optimize SE protocols, still requires three important upgrades to our screening system: (1) new equipment, e.g. high-throughput pipetting platforms; the platforms described above are coupled with a robotic arm and enable full automation of the screening process; (2) new laboratory strategies, e.g. preparation and maintenance of large amounts of embryogenic tissue; (3) new data acquisition strategies; as embryo development occurs in all 3 dimensions, embryo quantification is still difficult to analyse precisely with an automated 2-dimensional scanning microscope. To that end, early SE-markers would be very helpful.

Further high-throughput experiments will require a statistical methodology called “design of experiments”, which helps experimenters to predict how input variables (e.g. active compounds and nutritive compounds) interact to create output responses (e.g. embryo yield) in a process or system (e.g. regeneration process)⁴⁴. Full factorial and fractional factorial experiments are among the most useful multi-factor experiments for scientific investigations¹⁹,⁴⁵, in contradiction with a one factor-at-a-time approach used until now to optimize culture media in micropropagation, which requires many test runs and fails to consider any possible interaction between factors.

Thanks to masses of information obtained from “-omics” technologies making it possible to describe, in a global way, the pathways involved in the success of the key SE steps⁶, high-throughput screening can already be piloted in a rational way by selecting libraries able to alter target-molecules playing key roles within these pathways. In many cultivated species, including coffee, sequences of genomes, transcriptomes, proteomes and metabolomes are now available⁴⁶⁻⁴⁹. This strategy combining these 2 approaches – “-omics” technologies and robotic screening of active compounds – is currently being implemented by our teams on both cultivated species, *C. arabica* and *C. canephora*, in order to speed up the optimization of SE cloning processes and the dissemination of improved varieties.

Our focus was turned to the SE differentiation phase, where research tends to maximize embryo yields, but other key phases of the coffee SE process could also be easily optimized as a result of the “automated screening of active compound” method, such as the explant-to-callus phase, or the proliferation phase.

These results clearly demonstrate the feasibility, in plants, of automated screening for active compounds in miniaturized SE processes. This proof-of-concept represents a key step towards high-throughput optimization of cell culture media in plant micropropagation. This new method is particularly adapted to SE processes carried out in liquid media and, more widely, to all cell culture-based propagation processes.

## Material and Methods

**Plant material and experimental design to set up a miniaturized and automated screening platform.** Embryogenic cell clusters of *C. arabica* cv. GPFA116 were initiated from leaf explants obtained from 1-year-old plants obtained by SE and grown in the greenhouse of the Nestlé Plant Science Research Unit (Tours, France). The embryogenic cultures were cultivated following the protocols described previously for *C. arabica*⁵⁰. Briefly, the embryogenic cultures proliferated in 250-ml Erlenmeyer flasks containing M'proliferation' liquid nutritive medium, i.e. a medium supplemented with 0.3 mg/l 2,4-D (2,4-dichlorophenoxyacetic acid) and 1 mg/l BA (6-benzylaminopurine) (Fig. 7a). To stimulate differentiation of early somatic embryos, cultures are “classically” transferred to 250-ml Erlenmeyer flasks containing DIF 'embryo regeneration' liquid medium lacking the auxin 2,4-D at a rate of 0.1 g of embryogenic cell clusters per 100 ml of medium for 9 weeks. For miniaturization purposes, 24-well plates (ref. 734-2325, VWR International, Radnor - USA) were chosen as containers for optimum early embryo development. One ml of DIF medium was distributed per well. For automation purposes, the Viaflo Assist automated pipetting platform (Integra Bioscience Corp., Hudson - USA) equipped with a 6-channel pipette (Voyager, Integra Bioscience) was used. Disposable and sterile tips (Cat#4445, 1250 µl), as well as reservoirs (Cat#4322, 100 ml), were also purchased from Integra Bioscience. One g of Arabica cell clusters was first ground with a commercial blender (model: HGB2WT from Waring, Durrington - USA) with a sterilized MC-3 bowl (Waring Products, Durrington - USA) in 100 ml of differentiation medium for 30 sec at 18,000 rpm before filtration on a 700 µm diameter filter, in order to have particle diameters fitting the pipetting tips (Fig. 7b,c). The cell clusters were then inoculated in 24-well plates at an estimated rate of 1 mg per ml per well (Fig. 7d). The 6-channel pipette was programmed in order to have uniform distribution of cell clusters between plates, between rows, and between peripheral and central wells. This programme included a mixing step before pipetting, in order to avoid particle sedimentation. The plates were sealed with Parafilm (M from Bemis, Neenah - USA) and placed in the dark in an incubation shaker (CERTOMAT CT plus, Sartorius Stedim Biotech, Goettingen - Germany) set to 90% RH, 25 °C and 120 rpm for 9 weeks up to the regeneration of torpedo-shaped embryos (Fig. 7e). Uniformity of cluster distribution was assessed at the beginning of culturing, and embryo yield after 9 weeks. This experiment was carried out in 4 replicates with n = 24 wells/replicate.

At the same time, ground and intact cell clusters were cultured in DIF medium under standard *in vitro* conditions (n = 6 flasks/condition/replicate). Erlenmeyer flasks were placed on shakers (120 rpm) at 25 °C in the dark for 9 weeks. The effect of grinding on regeneration capacity was evaluated according to the number of torpedo-shaped embryos obtained per gram of inoculated callus. Total fresh weight was also reported. The experiment was carried out in triplicate.

**Viability test on cell clusters before and after grinding.** A viability test was carried out before and after grinding cell clusters. This test consisted of double staining with fluorescein diacetate (FDA) and propidium iodide (PI). In living cells, FDA is degraded by cytosolic esterases allowing it to fluoresce in green (emission at 517 nm) in living cells and PI fluoresces in red (emission at 636 nm) in dead cells. A staining solution was
<figure><img src="image_10.png" /><figcaption>**Figure 7.** The pipeline established to screen for active compounds affecting coffee SE by miniaturizing and automating the embryo regeneration step. (a) Embryogenic cell clusters in M 'proliferation' medium. Cell clusters had a heterogeneous diameter range (50–1000 µm). (b) Grinding of cell clusters for 30 sec at 18,000 rpm in DIF 'regeneration' medium. (c) Cell cluster filtration on a 700-µm diameter filter to obtain a uniform diameter range (50–250 µm). (d) Automated distribution of ground cell clusters in miniaturized culture systems (24-well plates). The same automated platform was used to screen for potential active compounds affecting SE. (e) Plates were placed in the dark in incubation shakers (90% RH 25°C and 100 rpm) for 9 weeks up to embryo regeneration. (f) An automated inverted plate-reading microscope coupled with a digital camera and computer-assisted image analysis were used to count cell clusters, measure the total area they occupied in a well and monitor embryo development.</figcaption></figure>

prepared with 5 ml of DIF medium, 8 µl of FDA at 5 g/l (F7378, Sigma Aldrich) and 50 µl of IP at 2 g/l (P4170, Sigma Aldrich). Cell clusters were incubated with 1 ml of staining solution for 10 minutes in the dark before washing 3 times with PBS pH 7.2. Clusters were then suspended in DIF medium and placed between 2 microscope slides for observation.

**Image acquisition and analysis.** Images from 24-well plates were automatically acquired using an inverted microscope equipped with a motorized scanning xy stage (Axio Observer D1, Zeiss, Jena – Germany) and a digital camera (Axiocam 503, Zeiss) (Fig. 7f). The software (Zen, Zeiss) was programmed to find the position of each well using auto-focusing (4x objective) and 8 horizontal images by 8 consecutive rows (8 × 8) were taken in order to register each entire well. The resulting collection of 64 reconstructed images per well was later processed using ImageJ freeware (NIH, Maryland – USA) for cell cluster counting and measurement of their diameter, or the total area they occupied in a well.

For viability tests, images of cell clusters were acquired using the same microscope equipped with dichroic filter cubes for the detection of fluorochromes and connected to a UV light source (HXP 120 v, Zeiss).

Images of torpedo-shaped embryos were acquired using a digital camera (Moticam 5+, Motic, Hong-Kong) mounted on a trinocular stereomicroscope (SZH10, Olympus Optical, Hamburg – Germany) and controlled by Motic Images Plus (v.3.0, Motic) software for length measurements.

**Screening of active compounds affecting SE by carrying out HDACi treatments in plates.** A full 4 × 4 factorial design of 18 wells replicated 3 times was used to implement the miniaturized bioassay. Four HDACi, Oxamflatin, SAHA, Scriptaid and Trichostastin A (O3139, SML0061, S7817 and T8552, respectively, from Sigma-Aldrich) were tested. Stock solutions of 10 mM were prepared by dissolving the HDACi powders in DMSO, then immediately sterilized by filtration (0.2 µm, Sartorius AG) and diluted in DIF differentiation medium to obtain final concentrations of 0, 0.1, 1 and 10 µM for each HDACi. Ground and filtered Arabica embryogenic clusters were prepared to inoculate 12 24-well plates at an inoculation density of 1 mg/ml of DIF differentiation medium using the Integra Bioscience Corp., Hudson – USA). For each HDACi, 3 24-well plates were filled with differentiation medium supplemented with 0 (6 wells of row A), 0.1
(wells of row B), 1 (wells of row C) and 10µM HDACi (wells of row D) respectively (n = 18 wells/HDACi/concentration/replicate). The plates were sealed with Parafilm and placed in an incubation shaker (CERTOMAT CT plus, Sartorius AG) at 90% RH, 120 rpm, 25°C, in the dark for 24 h. The medium was then gradually replaced with fresh DIF medium by replacing half of the well volume for 6 consecutive times using the pipetting platform (Integra Bioscience Corp.) in order to prevent loss of cell clusters. The plates were placed again in the CERTOMAT for 9 weeks under the same culture conditions. The response to treatment was evaluated according to the number of torpedo-shaped embryos obtained. This experiment was carried out in triplicate.

**TSA re-testing assay under standard in vitro culture conditions.** TSA at an optimum concentration of 1 µM was re-tested under standard conditions. To that end, Arabica embryogenic cell clusters were inoculated in 24 250-ml Erlenmeyer flasks at a rate of 0.1 g per 100 ml of DIF differentiation medium. Four exposure times to 1 µM TSA were tested (0, 1, 7 and 14 days), at the end of which the medium was replaced with fresh DIF differentiation medium (n = 6 flasks/exposure time/replicate). Flasks were placed on shakers (120 rpm) at 25 °C in the dark for 9 weeks. The response to treatment was evaluated according to the number of torpedo-shaped embryos obtained per gram of inoculated callus. Total fresh weight and mean length of embryos were also recorded. This experiment was carried out in triplicate.

**Plant phenotyping assay post-TSA treatment.** Torpedo-shaped embryos from cell clusters treated, or not, with TSA were placed on germination medium in 140 × 21 mm, vented Petri dishes (ref. 391-1500, VWR International) for further development. Dishes were maintained at 25 °C, under a 12/12 photoperiod of artificial light (24 µmol/s, Green Power TLED DR/W/MB 18 W, Philips Lighting, Zurich - Switzerland) delivering 65 µmol.m⁻².s⁻¹ for 2 to 3 months. After cotyledon development, green embryos were transplanted in 145 × 101 × 60 mm plastic boxes (ref: 017002, D. Dutscher S.A., Issy-Les-Moulineaux - France) containing growth medium for shoot and root development under the same conditions as for the Petri dishes. After 4 months, and for each TSA treatment, 40 SE-derived plantlets ready for greenhouse transfer from 3 boxes were assessed as follows: main root and stem lengths, number of pairs of leaves (lp, from 2 to 7) and roots (r, 0-4), number of stems (a, 1 or 2), growth type (orthotropic or plagiotropic), leaf pigmentation (uniform or variegated), leaf colour index (using a green colour chart, Pantone Carlstadt USA) and leaf shape (normal or abnormal).

**Statistical analysis.** All analyses were performed in R⁵¹ using the default functions. First, homogeneity between replicates was verified using a Kruskal-Wallis test (P > 0.05), then normality of data distribution and homogeneity of variances were assessed. If these conditions were verified, the ANOVA test was used to calculate any significant differences (at *P* = 0.05) between the groups (if groups >2), or a t-test was used (if groups = 2). If data did not follow a normal distribution, a Kruskal-Wallis test was used to calculate any significant differences (at *P* = 0.05) between the groups, with a Dunn test as post hoc (dunn.test package).

Based on acclimation data, ANOVA were performed on numerical variables, such as main root and stem lengths, to detect any after-effect of the TSA treatment on plant development. For categorical variables, a multiple component analysis (MCA), performed with Minitab v. 18.1, was used to identify correlations between TSA treatments and phenotypes.

Received: 12 May 2019; Accepted: 7 January 2020;

Published online: 21 January 2020

**References**

1. Guan, Y., Li, S. -G., Fan, X. -F. & Su, Z. -H. Application of Somatic Embryogenesis in Woody Plants. *Front. Plant Sci.* 7 (2016).

2. Harfouche, A. et al. Accelerating the domestication of forest trees in a changing world. *Trends Plant Sci.* 17, 64–72 (2012).

3. Merkle, S. A. & Dean, J. F. Forest tree biotechnology. *Curr. Opin. Biotechnol.* 11, 298–302 (2000).

4. Lelu-Walter, M.-A. et al. Somatic Embryogenesis for More Effective Breeding and Deployment of Improved Varieties in Pinus spp.: Bottlenecks and Recent Advances. in *Somatic Embryogenesis: Fundamental Aspects and Applications* (eds. Loyola-Vargas, V. M. & Ochoa-Alejo, N.) 319–365 (Springer International Publishing, 2016).

5. Lelu-Walter, M.-A. et al. Somatic embryogenesis in forestry with a focus on Europe: state of-the-art, benefits, challenges and future direction. *Tree Genet. Genomes* 9, 883–899 (2013).

6. Campos, N. A., Panis, B. & Carpentier, S. C. Somatic Embryogenesis in Coffee: The Evolution of Biotechnology and the Integration of Omics Technologies Offer Great Opportunities. *Front. Plant Sci.* 8 (2017).

7. Isah, T. Induction of somatic embryogenesis in woody plants. *Acta Physiol. Plant.* 38, 118 (2016).

8. Etienne, H. et al. 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 (2018).

9. Guillou, C. et al. Indirect somatic embryogenesis of Theobroma cacao L. in liquid medium and improvement of embryo-to-plantlet conversion rate. *Vitro Cell. Dev. Biol. - Plant* 54, 377–391 (2018).

10. Park, Y. S. et al. Initiation of somatic embryogenesis in Pinus banksiana, P. strobus, P. pinaster, and P. sylvestris at three laboratories in Canada and France. *Plant Cell Tissue Organ Cult.* 86, 87–101 (2006).

11. Bobadilla Landey, R. et al. Assessment of genetic and epigenetic changes during cell culture ageing and relations with somaclonal variation in Coffea arabica. *Plant Cell Tissue Organ Cult.* PCTOC 122, 517–531 (2015).

12. Maximova, S. N. et al. Efficiency, genotypic variability, and cellular origin of primary and secondary somatic embryogenesis of Theobroma cacao L. *Vitro Cell. Dev. Biol. - Plant* 38, 252–259 (2002).

13. Lamberth, C., Jeanmart, S., Luksch, T. & Plant, A. Current Challenges and Trends in the Discovery of Agrochemicals. *Science* 341, 742–746 (2013).

14. Boutilier, K., Angenent, G. C., Castan, M. S. & Hui, L. Inventors; Stichting Wageningen Research, assignee. Haploid embryogenesis. United States patent US 20160212956A1. 2016 Jul 28".

15. Corredoira, E., Ballester, A., Ibarra, M. & Vieitez, A. M. Induction of somatic embryogenesis in explants of shoot cultures established from adult Eucalyptus globulus and E. saligna × E. maidenii trees. *Tree Physiol.* 35, 678–690 (2015).

16. Moon, H. -K., Kim, Y. -W., Hong, Y. -P. & Park, S. -Y. Improvement of somatic embryogenesis and plantlet conversion in Oplopanax elatus, an endangered medicinal woody plant. *SpringerPlus* 2 (2013).
17. Maciel, A. L. R. et al. Large-scale, high-efficiency production of coffee somatic embryos. Crop Breed. Appl. Biotechnol. 16, 102–107 (2016).

18. Szymański, P., Markowicz, M. & Mikiciuk-Olasik, E. Adaptation of High-Throughput Screening in Drug Discovery—Toxicological Screening Tests. Int. J. Mol. Sci. 13, 427–452 (2011).

19. Radzun, K. A. et al. Automated nutrient screening system enables high-throughput optimisation of microalgae production conditions. Biotechnol. Biofuels 8, 65 (2015).

20. Robert, S., Raikhel, N. V. & Hicks, G. R. Powerful Partners: Arabidopsis and Chemical Genomics. Arab. Book 7, e0109 (2009).

21. Hicks, G. R. & Raikhel, N. V. Small Molecules Present Large Opportunities in Plant Biology. Annu. Rev. Plant Biol. 63, 261–282 (2012).

22. Connor, C. J. O., Beckmann, H. S. G. & Spring, D. R. Diversity-oriented synthesis: producing chemical tools for dissecting biology. Chem. Soc. Rev. 41, 4444–4456 (2012).

23. Feng, S., Jacobsen, S. E. & Reik, W. Epigenetic Reprogramming in Plant and Animal Development. Science 330, 622–627 (2010).

24. He, X.-J., Chen, T. & Zhu, J.-K. Regulation and function of DNA methylation in plants and animals. Cell Res. 21, 442–465 (2011).

25. Nic-Can, G. I. et al. New Insights into Somatic Embryogenesis: Leafy Cotyledon1, Baby Boom1 And Wuschel-Related Homeobox 4 Are Epigenetically Regulated in Coffea canephora. PLOS ONE 8, e72160 (2013).

26. Bertrand, B., Montagnon, C., Georget, F., Charmetant, P. & Etienne, H. Création et diffusion de variétés de caféiers Arabica: quelles innovations variétales? Cah. Agric. 21, 77–88 (2012).

27. Etienne, H. et al. 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, 115–124 (2012).

28. Ducos, J. P., Labbe, G., Lambot, C. & Pétiard, V. Pilot scale process for the production of pre-germinated somatic embryos of selected robusta (Coffea canephora) clones. *Vitro Cell. Dev. Biol. - Plant* **43**, 652–659 (2007).

29. Ducos, J. P. et al. Agronomic performance of Coffea canephora P. trees derived from large-scale somatic embryo production in liquid medium. *Euphytica* **131**, 215–223 (2003).

30. Georget, F. et al. Somatic embryogenesis-derived coffee plantlets can be efficiently propagated by horticultural rooted mini-cuttings: A boost for somatic embryogenesis. Sci. Hortic. 216, 177–185 (2017).

31. Hughes, J., Rees, S., Kalindjian, S. & Philpott, K. Principles of early drug discovery. Br. J. Pharmacol. 162, 1239–1249 (2011).

32. Mayr, L. M. & Fuerst, P. The future of high-throughput screening. J. Biomol. Screen. 13, 443–448 (2008).

33. Rodriguez-Furlán, C., Miranda, G., Reggiardo, M., Hicks, G. R. & Norambuena, L. High throughput selection of novel plant growth regulators: Assessing the translatability of small bioactive molecules from Arabidopsis to crops. Plant Sci. 245, 50–60 (2016).

34. Gao, G., Chen, L. & Huang, C. Anti-cancer Drug Discovery: Update and Comparisons in Yeast, Drosophila, and Zebrafish. Curr. Mol. Pharmacol. 7, 44–51 (2014).

35. Halder, V. & Kombrink, E. Facile high-throughput forward chemical genetic screening by in situ monitoring of glucuronidase-based reporter gene expression in Arabidopsis thaliana. Front. Plant Sci. 6, 13 (2015).

36. Zhang, H., Lamping, S. R., Pickering, S. C. R., Lye, G. J. & Shamlou, P. A. Engineering characterisation of a single well from 24-well and 96-well microtitre plates. Biochem. Eng. J. 40, 138–149 (2008).

37. Etienne, H., Dechamp, E., Barry-Etienne, D. & Bertrand, B. Bioreactors in coffee micropropagation. Braz. J. Plant Physiol. 18, 45–54 (2006).

38. Li, H. et al. The Histone Deacetylase Inhibitor Trichostatin A Promotes Totipotency in the Male Gametophyte[W]. Plant Cell 26, 195–209 (2014).

39. Wójcikowska, B. et al. Trichostatin A Triggers an Embryogenic Transition in Arabidopsis Explants via an Auxin-Related Pathway. Front. Plant Sci. 9 (2018).

40. Uddenberg, D. et al. Embryogenic potential and expression of embryogenesis-related genes in conifers are affected by treatment with a histone deacetylase inhibitor. BMC Proc. 5, P151 (2011).

41. Abrahamsson, M., Valladares, S., Merino, I., Larsson, E. & von Arnold, S. Degeneration pattern in somatic embryos of Pinus sylvestris L. *Vitro Cell. Dev. Biol. - Plant* **53**, 86–96 (2017).

42. Kong, F., Yuan, L., Zheng, Y. F. & Chen, W. Automatic Liquid Handling for Life Science: A Critical Review of the Current State of the Art. J. Lab. Autom. 17, 169–185 (2012).

43. Michael, S. et al. A Robotic Platform for Quantitative High-Throughput Screening. *Assay Drug Dev. Technol.* **6**, 637–657 (2008).

44. Ilzarbe, L., Álvarez, M. J., Viles, E. & Tanco, M. Practical applications of design of experiments in the field of engineering: a bibliographical review. Qual. Reliab. Eng. Int. 24, 417–428 (2008).

45. Altekar, M. et al. Assay Optimization: A Statistical Design of Experiments Approach. Clin. Lab. Med. 27, 139–154 (2007).

46. Coffea arabica genome. Available at, https://worldcoffeeresearch.org/work/coffea-arabica-genome/ (Accessed: 6th March 2019).

47. Ivamoto, S. T. et al. Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis. PLoS ONE 12 (2017).

48. Campos, N. A., Paiva, L. V., Panis, B. & Carpentier, S. C. The proteome profile of embryogenic cell suspensions of Coffea arabica L. *Proteomics* 16, 1001–1005 (2016).

49. Businge, E., Brackmann, K., Moritz, T. & Egertsdotter, U. Metabolite profiling reveals clear metabolic changes during somatic embryo development of Norway spruce (Picea abies). Tree Physiol. 32, 232–244 (2012).

50. Etienne, H. Somatic Embryogenesis Protocol: Coffee (Coffea arabica L. and C. canephora P.). in Protocols for Somatic Embryogenesis in Woody Plants (eds. Jain, S. M. & Gupta, P. K.) Springer International Publishing, 2005).

51. R-Core-Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2018).

## Acknowledgements

Rayan Awada benefits from a PhD grant and funding from Nestlé and CIRAD.

## Author contributions

D.B. and S.F.M. designed the research. R.A., D.V., S.F. and E.B. performed the experiments. R.A. and D.B. analysed the data. R.A., D.B. and H.E. wrote the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-57800-6.

Correspondence and requests for materials should be addressed to R.A.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020