Identification of Hydroxypyrazine O-Methyltransferase Genes in Coffea arabica: A Potential Source of Methoxypyrazines That Cause Potato Taste Defect

The goal of this study is to identify Coffea arabica O-methyltransferase (OMT) genes involved in the biosynthesis of methoxypyrazines. High levels of 2-isopropyl-3-methoxypyrazine (IPMP) and 2-isobutyl-3-methoxypyrazine (IBMP) in coffee beans are associated with the potato taste defect (PTD). Among the 34 putative O-methyltransferase genes identified in the published genome of C. canephora, three genes are highly homologous to known hydroxypyrazine OMT genes. Genes of interest were amplified and sequenced from genomic DNA of single C. arabica beans grown in eight different locations, including regions with endemic PTD. Although C. arabica OMT target sequences were almost identical regardless of source location, individual beans shared numerous polymorphisms in each of the target genes. Two of the predicted C. arabica OMT enzymes were successfully expressed in Escherichia coli and purified, and one enzyme shows slow yet measurable turnover of both 3-isobutyl-2-hydroxypyrazine (IBHP) and 3-isopropyl-2- hydroxypyrazine (IPHP), supporting a possible role of the coffee plant in PTD.

Mechanistic insights from plant heteromannan synthesis in yeast

Heteromannan (HM) is one of the most ancient cell wall polymers in the plant kingdom, consisting of β-(1-4)-linked backbones of glucose (Glc) and mannose (Man) units. Despite the widespread distribution of HM polysaccharides, their biosynthesis remains mechanistically unclear. HM is elongated by glycosyltransferases (GTs) from the cellulose synthase-like A (CSLA) family. MANNAN-SYNTHESIS RELATED (MSR) putative GTs have also been implicated in (gluco)mannan synthesis, but their roles have been difficult to decipher in planta and in vitro. To further characterize the products of the HM synthases and accessory proteins, we chose a synthetic biology approach to synthesize plant HM in yeast. The expression of a CSLA protein in Pichia pastoris led to the abundant production of plant HM: up to 30% of glycans in the yeast cell wall. Based on sequential chemical and enzymatic extractions, followed by detailed structural analyses, the newly produced HM polymers were unbranched and could be larger than 270 kDa. Using CSLAs from different species, we programmed yeast cells to produce an HM backbone composed exclusively of Man or also incorporating Glc. We demonstrate that specific MSR cofactors were indispensable for mannan synthase activity of a coffee CSLA or modulated a functional CSLA enzyme to produce glucomannan instead of mannan. Therefore, this powerful platform yields functional insight into the molecular machinery required for HM biosynthesis in plants.

Use of a draft genome of coffee (Coffea arabica) to identify SNPs associated with caffeine content

Arabica coffee (Coffea arabica) has a small gene pool limiting genetic improvement. Selection for caffeine content within this gene pool would be assisted by identification of the genes controlling this important trait. Sequencing of DNA bulks from 18 genotypes with extreme high- or low-caffeine content from a population of 232 genotypes was used to identify linked polymorphisms. To obtain a reference genome, a whole genome assembly of arabica coffee (variety K7) was achieved by sequencing using short read (Illumina) and long-read (PacBio) technology. Assembly was performed using a range of assembly tools resulting in 76 409 scaffolds with a scaffold N50 of 54 544 bp and a total scaffold length of 1448 Mb. Validation of the genome assembly using different tools showed high completeness of the genome. More than 99% of transcriptome sequences mapped to the C. arabica draft genome, and 89% of BUSCOs were present. The assembled genome annotated using AUGUSTUS yielded 99 829 gene models. Using the draft arabica genome as reference in mapping and variant calling allowed the detection of 1444 nonsynonymous single nucleotide polymorphisms (SNPs) associated with caffeine content. Based on Kyoto Encyclopaedia of Genes and Genomes pathway-based analysis, 65 caffeine-associated SNPs were discovered, among which 11 SNPs were associated with genes encoding enzymes involved in the conversion of substrates, which participate in the caffeine biosynthesis pathways. This analysis demonstrated the complex genetic control of this key trait in coffee.

Identification of a seed maturation protein gene from Coffea arabica (CaSMP) and analysis of its promoter activity in tomato

A seed maturation protein gene (CaSMP) from Coffea arabica is expressed in the endosperm of yellow/green fruits. The CaSMP promoter drives reporter expression in the seeds of immature tomato fruits. In this report, an expressed sequence tag-based approach was used to identify a seed-specific candidate gene for promoter isolation in Coffea arabica. The tissue-specific expression of the cognate gene (CaSMP), which encodes a yet uncharacterized coffee seed maturation protein, was validated by RT-qPCR. Additional expression analysis during coffee fruit development revealed higher levels of CaSMP transcript accumulation in the yellow/green phenological stage. Moreover, CaSMP was preferentially expressed in the endosperm and was down-regulated during water imbibition of the seeds. The presence of regulatory cis-elements known to be involved in seed- and endosperm-specific expression was observed in the CaSMP 5'-upstream region amplified by genome walking (GW). Additional histochemical analysis of transgenic tomato (cv. Micro-Tom) lines harboring the GW-amplified fragment (~ 1.4 kb) fused to uidA reporter gene confirmed promoter activity in the ovule of immature tomato fruits, while no activity was observed in the seeds of ripening fruits and in the other organs/tissues examined. These results indicate that the CaSMP promoter can be used to drive transgene expression in coffee beans and tomato seeds, thus representing a promising biotechnological tool.

Ectopic expression of the Coffea canephora SERK1 homolog-induced differential transcription of genes involved in auxin metabolism and in the developmental control of embryogenesis

Somatic embryogenesis receptor-like kinase 1 (SERK1) is a membrane receptor that might serve as common co-regulator of plant cell differentiation processes by forming heterodimers with specific receptor-like kinases. The Coffea canephora SERK1 homolog (CcSERK1) was cloned in this work, and its early function in the transcription of embryogenesis master genes and of genes encoding proteins involved in auxin metabolism was investigated by externally manipulating its expression in embryogenic leaf explants, before the appearance of embryogenic structures. Overexpression of CcSERK1 early during embryogenesis caused an increase in the number of somatic embryos when the 55-day process was completed. Suppression of CcSERK1 expression by RNA interference almost abolished somatic embryogenesis. Real time-PCR experiments revealed that the transcription of the CcAGL15, CcWUS, CcBBM, CcPKL, CcYUC1, CcPIN1 and CcPIN4 homologs was modified in direct proportion to the expression of CcSERK1 and that only CcLEC1 was inversely affected by the expression levels of CcSERK1. The expression of the CcYUC4 homolog was induced to more than 80-fold under CcSERK1 overexpression conditions, but it was also induced when CcSERK1 expression was silenced. The level of CcTIR1 was not affected by CcSERK1 overexpression but was almost abolished during CcSERK1 silencing. These results suggest that CcSERK1 co-regulates the induction of somatic embryogenesis in Coffea canephora by early activation of YUC-dependent auxin biosynthesis, auxin transport mediated by PIN1 and PIN4, and probably auxin perception by the TIR1 receptor, leading to the induction of early-stage homeotic genes (CcAGL15, CcWUS, CcPKL and CcBBM) and repression of late-stage homeotic genes (CcLec1).

The coffee bean transcriptome explains the accumulation of the major bean components through ripening

The composition of the maturing coffee bean determines the processing performance and ultimate quality of the coffee produced from the bean. Analysis of differences in gene expression during bean maturation may explain the basis of genetic and environmental variation in coffee quality. The transcriptome of the coffee bean was analyzed at three stages of development, immature (green), intermediate (yellow) and mature (red). A total of more than 120 million 150 bp paired-end reads were collected by sequencing of transcripts of triplicate samples at each developmental stage. A greater number of transcripts were expressed at the yellow stage. As the beans matured the types of highly expressed transcripts changed from transcripts predominantly associated with galactomannan, triacylglycerol (TAG), TAG lipase, 11 S and 7S-like storage protein and Fasciclin-like arabinogalactan protein 17 (FLA17) in green beans to transcripts related to FLA1 at the yellow stage and TAG storage lipase SDP1, and SDP1-like in red beans. This study provides a genomic resource that can be used to investigate the impact of environment and genotype on the bean transcriptome and develop coffee varieties and production systems that are better adapted to deliver quality coffee despite climate variations.

Rust and Thinning Management Effect on Cup Quality and Plant Performance for Two Cultivars of Coffea arabica L

Beverage quality is a complex attribute of coffee ( Coffea arabica L.). Genotype (G), environment (E), management (M), postharvest processing, and roasting are all involved. However, little is known about how G × M interactions influence beverage quality. We investigated how yield and coffee leaf rust (CLR) disease (caused by Hemileia vastatrix Berk. et Br.) management affect cup quality and plant performance, in two coffee cultivars. Sensory and chemical analyses revealed that 10 of 70 attributes and 18 of 154 chemical volatile compounds were significantly affected by G and M. Remarkably, acetaminophen was found for the first time in roasted coffee and in higher concentrations under more stressful conditions. A principal component analysis described 87% of the variation in quality and plant overall performance. This study is a first step in understanding the complexity of the physiological, metabolic, and molecular changes in coffee production, which will be useful for the improvement of coffee cultivars.

Integrative analysis of the late maturation programme and desiccation tolerance mechanisms in intermediate coffee seeds

The 'intermediate seed' category was defined in the early 1990s using coffee (Coffea arabica) as a model. In contrast to orthodox seeds, intermediate seeds cannot survive complete drying, which is a major constraint for seed storage and has implications for both biodiversity conservation and agricultural purposes. However, intermediate seeds are considerably more tolerant to drying than recalcitrant seeds, which are highly sensitive to desiccation. To gain insight into the mechanisms governing such differences, changes in desiccation tolerance (DT), hormone contents, and the transcriptome were analysed in developing coffee seeds. Acquisition of DT coincided with a dramatic transcriptional switch characterised by the repression of primary metabolism, photosynthesis, and respiration, and the up-regulation of genes coding for late-embryogenesis abundant (LEA) proteins, heat-shock proteins (HSPs), and antioxidant enzymes. Analysis of the heat-stable proteome in mature coffee seeds confirmed the accumulation of LEA proteins identified at the transcript level. Transcriptome analysis also suggested a major role for ABA and for the transcription factors CaHSFA9, CaDREB2G, CaANAC029, CaPLATZ, and CaDOG-like in DT acquisition. The ability of CaHSFA9 and CaDREB2G to trigger HSP gene transcription was validated by Agrobacterium-mediated transformation of coffee somatic embryos.

Auxin Immunolocalization in Coffea canephora Tissues

Auxins are plant growth regulators that participate in a variety of biological mechanisms during the growth and development of plants. The most abundant natural auxin is indole-3-acetic acid (IAA). The physiological processes regulated by IAA depend on their temporal space accumulation in different tissues of a plant. This accumulation is regulated by its biosynthesis, conjugation, degradation, and transport. Therefore tools that allow us a qualitative and quantitative detection of IAA in plant tissues are very useful to understand the homeostasis of IAA during the life cycle of plants. In this protocol, the complete procedure for localization of IAA in different tissues of Coffea canephora is described using specific anti-IAA monoclonal antibodies.

An endoplasmic reticulum-localized Coffea arabica BURP domain-containing protein affects the response of transgenic Arabidopsis plants to diverse abiotic stresses

The Coffea arabica BURP domain-containing gene plays an important role in the response of transgenic Arabidopsis plants to abiotic stresses via regulating the level of diverse proteins. Although the functions of plant-specific BURP domain-containing proteins (BDP) have been determined for a few plants, their roles in the growth, development, and stress responses of most plant species, including coffee plant (Coffea arabica), are largely unknown. In this study, the function of a C. arabica BDP, designated CaBDP1, was investigated in transgenic Arabidopsis plants. The expression of CaBDP1 was highly modulated in coffee plants subjected to drought, cold, salt, or ABA. Confocal analysis of CaBDP1-GFP fusion proteins revealed that CaBDP1 is localized in the endoplasmic reticulum. The ectopic expression of CaBDP1 in Arabidopsis resulted in delayed germination of the transgenic plants under abiotic stress and in the presence of ABA. Cotyledon greening and seedling growth of the transgenic plants were inhibited in the presence of ABA due to the upregulation of ABA signaling-related genes like ABI3, ABI4, and ABI5. Proteome analysis revealed that the levels of several proteins are modulated in CaBDP1-expressing transgenic plants. The results of this study underscore the importance of BURP domain proteins in plant responses to diverse abiotic stresses.

Photosynthetic and metabolic acclimation to repeated drought events play key roles in drought tolerance in coffee

Over the last decades, most information on the mechanisms underlying tolerance to drought has been gained by considering this stress as a single event that happens just once in the life of a plant, in contrast to what occurs under natural conditions where recurrent drought episodes are the rule. Here we explored mechanisms of drought tolerance in coffee (Coffea canephora) plants from a broader perspective, integrating key aspects of plant physiology and biochemistry. We show that plants exposed to multiple drought events displayed higher photosynthetic rates, which were largely accounted for by biochemical rather than diffusive or hydraulic factors, than those submitted to drought for the first time. Indeed, these plants displayed higher activities of RuBisCO and other enzymes associated with carbon and antioxidant metabolism. Acclimation to multiple drought events involved the expression of trainable genes related to drought tolerance and was also associated with a deep metabolite reprogramming with concordant alterations in central metabolic processes such as respiration and photorespiration. Our results demonstrate that plants exposed to multiple drought cycles can develop a differential acclimation that potentiates their defence mechanisms, allowing them to be kept in an 'alert state' to successfully cope with further drought events.

Differential fine-tuning of gene expression regulation in coffee leaves by CcDREB1D promoter haplotypes under water deficit

Despite the importance of the DREB1D gene (also known as CBF4) in plant responses to water deficit and cold stress, studies analysing its regulation by transgenic approaches are lacking. In the current work, a functional study of three CcDREB1D promoter haplotypes (named HP15, HP16 and HP17) isolated from drought-tolerant and drought-sensitive clones of Coffea canephora was carried out in plants of C. arabica stably transformed by Agrobacterium tumefaciens by analysing their ability to regulate the expression of the uidA reporter gene in response to water deficit mimicked by polyethylene glycol (-2.0 MPa) and low relative humidity treatments. A deletion analysis of their corresponding 5'-upstream regions revealed increased specificity of β-glucuronidase activity in the polyethylene glycol and low relative humidity treatments, with high expression in leaf mesophyll and guard cells in full-length constructs. RT-qPCR assays also revealed that the HP16 haplotype (specific to clone tolerant to water deficit) had stronger and earlier activity compared with the HP15 and HP17 haplotypes. As most of the cis-regulatory elements involved in ABA-dependent and -independent networks, tissue specificity and light regulation are common to these haplotypes, we propose that their organization, as well as the nucleic acid polymorphisms present outside these boxes, may play a role in modulating activities of DREB1D promoters in guard cells.

A humidity shock leads to rapid, temperature dependent changes in coffee leaf physiology and gene expression

Climate change is expected to increase the frequency of above-normal atmospheric water deficits contemporaneous with periods of high temperatures. Here we explore alterations in physiology and gene expression in leaves of Coffea canephora Pierre ex A. Froehner caused by a sharp drop in relative humidity (RH) at three different temperatures. Both stomatal conductance (gs) and CO2 assimilation (A) measurements showed that gs and A values fell quickly at all temperatures after the transfer to low RH.  However, leaf relative water content measurements indicated that leaves nonetheless experienced substantial water losses, implying that stomatal closure and/or resupply of water was not fast enough to stop excessive evaporative losses.  At 27 and 35 °C, upper leaves showed significant decreases in Fv/Fm compared with lower leaves, suggesting a stronger impact on photosystem II for upper leaves, while at 42 °C, both upper and lower leaves were equally affected. Quantitative gene expression analysis of transcription factors associated with conventional dehydration stress, and genes involved with abscisic acid signalling, such as CcNCED3, indicated temperature-dependent, transcriptional changes during the Humidity Shock ('HuS') treatments.  No expression was seen at 27 °C for the heat-shock gene CcHSP90-7, but it was strongly induced during the 42 °C 'HuS' treatment. Consistent with a proposal that important cellular damage occurred during the 42 °C 'HuS' treatment, two genes implicated in senescence were induced by this treatment. Overall, the data show that C. canephora plants subjected to a sharp drop in RH exhibit major, temperature-dependent alterations in leaf physiology and important changes in the expression of genes associated with abiotic stress and senescence. The results presented suggest that more detailed studies on the combined effects of low RH and high temperature are warranted.

Detection of plant microRNAs in honey

For the first time in the literature, our group has managed to demonstrate the existence of plant RNAs in honey samples. In particular, in our work, different RNA extraction procedures were performed in order to identify a purification method for nucleic acids from honey. Purity, stability and integrity of the RNA samples were evaluated by spectrophotometric, PCR and electrophoretic analyses. Among all honey RNAs, we specifically revealed the presence of both plastidial and nuclear plant transcripts: RuBisCO large subunit mRNA, maturase K messenger and 18S ribosomal RNA. Surprisingly, nine plant microRNAs (miR482b, miR156a, miR396c, miR171a, miR858, miR162a, miR159c, miR395a and miR2118a) were also detected and quantified by qPCR. In this context, a comparison between microRNA content in plant samples (i.e. flowers, nectars) and their derivative honeys was carried out. In addition, peculiar microRNA profiles were also identified in six different monofloral honeys. Finally, the same plant microRNAs were investigated in other plant food products: tea, cocoa and coffee. Since plant microRNAs introduced by diet have been recently recognized as being able to modulate the consumer’s gene expression, our research suggests that honey’s benefits for human health may be strongly correlated to the bioactivity of plant microRNAs contained in this matrix.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

The Hemileia vastatrix effector HvEC-016 suppresses bacterial blight symptoms in coffee genotypes with the SH 1 rust resistance gene

A number of genes that confer resistance to coffee leaf rust (S_H 1-S_H 9) have been identified within the genus Coffea, but despite many years of research on this pathosystem, the complementary avirulence genes of Hemileia vastatrix have not been reported. After identification of H. vastatrix effector candidate genes (HvECs) expressed at different stages of its lifecycle, we established an assay to characterize HvEC proteins by delivering them into coffee cells via the type-three secretion system (T3SS) of Pseudomonas syringae pv. garcae (Psgc). Employing a calmodulin-dependent adenylate cyclase assay, we demonstrate that Psgc recognizes a heterologous P. syringae T3SS secretion signal which enables us to translocate HvECs into the cytoplasm of coffee cells. Using this Psgc-adapted effector detector vector (EDV) system, we found that HvEC-016 suppresses the growth of Psgc on coffee genotypes with the S_H 1 resistance gene. Suppression of bacterial blight symptoms in S_H 1 plants was associated with reduced bacterial multiplication. By contrast, HvEC-016 enhanced bacterial multiplication in S_H 1-lacking plants. Our findings suggest that HvEC-016 may be recognized by the plant immune system in a S_H 1-dependent manner. Thus, our experimental approach is an effective tool for the characterization of effector/avirulence proteins of this important pathogen.

Diterpenes biochemical profile and transcriptional analysis of cytochrome P450s genes in leaves, roots, flowers, and during Coffea arabica L. fruit development

Lipids are among the major chemical compounds present in coffee beans, and they affect the flavor and aroma of the coffee beverage. Coffee oil is rich in kaurene diterpene compounds, mainly cafestol (CAF) and kahweol (KAH), which are related to plant defense mechanisms and to nutraceutical and sensorial beverage characteristics. Despite their importance, the final steps of coffee diterpenes biosynthesis remain unknown. To understand the molecular basis of coffee diterpenes biosynthesis, we report the content dynamics of CAF and KAH in several Coffea arabica tissues and the transcriptional analysis of cytochrome P450 genes (P450). We measured CAF and KAH concentrations in leaves, roots, flower buds, flowers and fruit tissues at seven developmental stages (30-240 days after flowering - DAF) using HPLC. Higher CAF levels were detected in flower buds and flowers when compared to fruits. In contrast, KAH concentration increased along fruit development, peaking at 120 DAF. We did not detect CAF or KAH in leaves, and higher amounts of KAH than CAF were detected in roots. Using P450 candidate genes from a coffee EST database, we performed RT-qPCR transcriptional analysis of leaves, flowers and fruits at three developmental stages (90, 120 and 150 DAF). Three P450 genes (CaCYP76C4, CaCYP82C2 and CaCYP74A1) had transcriptional patterns similar to CAF concentration and two P450 genes (CaCYP71A25 and CaCYP701A3) have transcript accumulation similar to KAH concentration. These data warrant further investigation of these P450s as potential candidate genes involved in the final stages of the CAF and KAH biosynthetic pathways.

Identification by the DArTseq method of the genetic origin of the Coffea canephora cultivated in Vietnam and Mexico

BACKGROUND

The coffee species Coffea canephora is commercially identified as "Conilon" when produced in Brazil, or "Robusta" when produced elsewhere in the world. It represents approximately 40 % of coffee production worldwide. While the genetic diversity of wild C. canephora has been well studied in the past, only few studies have addressed the genetic diversity of currently cultivated varieties around the globe. Vietnam is the largest Robusta producer in the world, while Mexico is the only Latin American country, besides Brazil, that has a significant Robusta production. Knowledge of the genetic origin of Robusta cultivated varieties in countries as important as Vietnam and Mexico is therefore of high interest.

RESULTS

Through the use of Sequencing-based diversity array technology-DArTseq method-on a collection of C. canephora composed of known accessions and accessions cultivated in Vietnam and Mexico, 4,021 polymorphic SNPs were identified. We used a multivariate analysis using SNP data from reference accessions in order to confirm and further fine-tune the genetic diversity of C. canephora. Also, by interpolating the data obtained for the varieties from Vietnam and Mexico, we determined that they are closely related to each other, and identified that their genetic origin is the Robusta Congo - Uganda group.

CONCLUSIONS

The genetic characterization based on SNP markers of the varieties grown throughout the world, increased our knowledge on the genetic diversity of C. canephora, and contributed to the understanding of the genetic background of varieties from very important coffee producers. Given the common genetic origin of the Robusta varieties cultivated in Vietnam, Mexico and Uganda, and the similar characteristics of climatic areas and relatively high altitude where they are grown, we can state that the Vietnamese and the Mexican Robusta have the same genetic potential to produce good cup quality.

Abiotic stresses affect differently the intron splicing and expression of chloroplast genes in coffee plants (Coffea arabica) and rice (Oryza sativa)

Despite the increasing understanding of the regulation of chloroplast gene expression in plants, the importance of intron splicing and processing of chloroplast RNA transcripts under stress conditions is largely unknown. Here, to understand how abiotic stresses affect the intron splicing and expression patterns of chloroplast genes in dicots and monocots, we carried out a comprehensive analysis of the intron splicing and expression patterns of chloroplast genes in the coffee plant (Coffea arabica) as a dicot and rice (Oryza sativa) as a monocot under abiotic stresses, including drought, cold, or combined drought and heat stresses. The photosynthetic activity of both coffee plants and rice seedlings was significantly reduced under all stress conditions tested. Analysis of the transcript levels of chloroplast genes revealed that the splicing of tRNAs and mRNAs in coffee plants and rice seedlings were significantly affected by abiotic stresses. Notably, abiotic stresses affected differently the splicing of chloroplast tRNAs and mRNAs in coffee plants and rice seedlings. The transcript levels of most chloroplast genes were markedly downregulated in both coffee plants and rice seedlings upon stress treatment. Taken together, these results suggest that coffee and rice plants respond to abiotic stresses via regulating the intron splicing and expression of different sets of chloroplast genes.

Advances in genomics for the improvement of quality in coffee

Coffee is an important crop that provides a livelihood to millions of people living in developing countries. Production of genotypes with improved coffee quality attributes is a primary target of coffee genetic improvement programmes. Advances in genomics are providing new tools for analysis of coffee quality at the molecular level. The recent report of a genomic sequence for robusta coffee, Coffea canephora, is a major development. However, a reference genome sequence for the genetically more complex arabica coffee (C. arabica) will also be required to fully define the molecular determinants controlling quality in coffee produced from this high quality coffee species. Genes responsible for control of the levels of the major biochemical components in the coffee bean that are known to be important in determining coffee quality can now be identified by association analysis. However, the narrow genetic base of arabica coffee suggests that genomics analysis of the wild relatives of coffee (Coffea spp.) may be required to find the phenotypic diversity required for effective association genetic analysis. The genomic resources available for the study of coffee quality are described and the potential for the application of next generation sequencing and association genetic analysis to advance coffee quality research are explored. © 2016 Society of Chemical Industry.

Identification of candidate genes for drought tolerance in coffee by high-throughput sequencing in the shoot apex of different Coffea arabica cultivars

BACKGROUND

Drought is a widespread limiting factor in coffee plants. It affects plant development, fruit production, bean development and consequently beverage quality. Genetic diversity for drought tolerance exists within the coffee genus. However, the molecular mechanisms underlying the adaptation of coffee plants to drought are largely unknown. In this study, we compared the molecular responses to drought in two commercial cultivars (IAPAR59, drought-tolerant and Rubi, drought-susceptible) of Coffea arabica grown in the field under control (irrigation) and drought conditions using the pyrosequencing of RNA extracted from shoot apices and analysing the expression of 38 candidate genes.

RESULTS

Pyrosequencing from shoot apices generated a total of 34.7 Mbp and 535,544 reads enabling the identification of 43,087 clusters (41,512 contigs and 1,575 singletons). These data included 17,719 clusters (16,238 contigs and 1,575 singletons) exclusively from 454 sequencing reads, along with 25,368 hybrid clusters assembled with 454 sequences. The comparison of DNA libraries identified new candidate genes (n = 20) presenting differential expression between IAPAR59 and Rubi and/or drought conditions. Their expression was monitored in plagiotropic buds, together with those of other (n = 18) candidates genes. Under drought conditions, up-regulated expression was observed in IAPAR59 but not in Rubi for CaSTK1 (protein kinase), CaSAMT1 (SAM-dependent methyltransferase), CaSLP1 (plant development) and CaMAS1 (ABA biosynthesis). Interestingly, the expression of lipid-transfer protein (nsLTP) genes was also highly up-regulated under drought conditions in IAPAR59. This may have been related to the thicker cuticle observed on the abaxial leaf surface in IAPAR59 compared to Rubi.

CONCLUSIONS

The full transcriptome assembly of C. arabica, followed by functional annotation, enabled us to identify differentially expressed genes related to drought conditions. Using these data, candidate genes were selected and their differential expression profiles were confirmed by qPCR experiments in plagiotropic buds of IAPAR59 and Rubi under drought conditions. As regards the genes up-regulated under drought conditions, specifically in the drought-tolerant IAPAR59, several corresponded to orphan genes but also to genes coding proteins involved in signal transduction pathways, as well as ABA and lipid metabolism, for example. The identification of these genes should help advance our understanding of the genetic determinism of drought tolerance in coffee.

Sustained enhancement of photosynthesis in coffee trees grown under free-air CO2 enrichment conditions: disentangling the contributions of stomatal, mesophyll, and biochemical limitations

Coffee (Coffea spp.), a globally traded commodity, is a slow-growing tropical tree species that displays an improved photosynthetic performance when grown under elevated atmospheric CO2 concentrations ([CO2]). To investigate the mechanisms underlying this response, two commercial coffee cultivars (Catuaí and Obatã) were grown using the first free-air CO2 enrichment (FACE) facility in Latin America. Measurements were conducted in two contrasting growth seasons, which were characterized by the high (February) and low (August) sink demand. Elevated [CO2] led to increases in net photosynthetic rates (A) in parallel with decreased photorespiration rates, with no photochemical limitations to A. The stimulation of A by elevated CO2 supply was more prominent in August (56% on average) than in February (40% on average). Overall, the stomatal and mesophyll conductances, as well as the leaf nitrogen and phosphorus concentrations, were unresponsive to the treatments. Photosynthesis was strongly limited by diffusional constraints, particularly at the stomata level, and this pattern was little, if at all, affected by elevated [CO2]. Relative to February, starch pools (but not soluble sugars) increased remarkably (>500%) in August, with no detectable alteration in the maximum carboxylation capacity estimated on a chloroplast [CO2] basis. Upregulation of A by elevated [CO2] took place with no signs of photosynthetic downregulation, even during the period of low sink demand, when acclimation would be expected to be greatest.

Multiclass Classification of Agro-Ecological Zones for Arabica Coffee: An Improved Understanding of the Impacts of Climate Change

Cultivation of Coffea arabica is highly sensitive to and has been shown to be negatively impacted by progressive climatic changes. Previous research contributed little to support forward-looking adaptation. Agro-ecological zoning is a common tool to identify homologous environments and prioritize research. We demonstrate here a pragmatic approach to describe spatial changes in agro-climatic zones suitable for coffee under current and future climates. We defined agro-ecological zones suitable to produce arabica coffee by clustering geo-referenced coffee occurrence locations based on bio-climatic variables. We used random forest classification of climate data layers to model the spatial distribution of these agro-ecological zones. We used these zones to identify spatially explicit impact scenarios and to choose locations for the long-term evaluation of adaptation measures as climate changes. We found that in zones currently classified as hot and dry, climate change will impact arabica more than those that are better suited to it. Research in these zones should therefore focus on expanding arabica's environmental limits. Zones that currently have climates better suited for arabica will migrate upwards by about 500m in elevation. In these zones the up-slope migration will be gradual, but will likely have negative ecosystem impacts. Additionally, we identified locations that with high probability will not change their climatic characteristics and are suitable to evaluate C. arabica germplasm in the face of climate change. These locations should be used to investigate long term adaptation strategies to production systems.

The Greater Phenotypic Homeostasis of the Allopolyploid Coffea arabica Improved the Transcriptional Homeostasis Over that of Both Diploid Parents

Polyploidy impacts the diversity of plant species, giving rise to novel phenotypes and leading to ecological diversification. In order to observe adaptive and evolutionary capacities of polyploids, we compared the growth, primary metabolism and transcriptomic expression level in the leaves of the newly formed allotetraploid Coffea arabica species compared with its two diploid parental species (Coffea eugenioides and Coffea canephora), exposed to four thermal regimes (TRs; 18-14, 23-19, 28-24 and 33-29°C). The growth rate of the allopolyploid C. arabica was similar to that of C. canephora under the hottest TR and that of C. eugenioides under the coldest TR. For metabolite contents measured at the hottest TR, the allopolyploid showed similar behavior to C. canephora, the parent which tolerates higher growth temperatures in the natural environment. However, at the coldest TR, the allopolyploid displayed higher sucrose, raffinose and ABA contents than those of its two parents and similar linolenic acid leaf composition and Chl content to those of C. eugenioides. At the gene expression level, few differences between the allopolyploid and its parents were observed for studied genes linked to photosynthesis, respiration and the circadian clock, whereas genes linked to redox activity showed a greater capacity of the allopolyploid for homeostasis. Finally, we found that the overall transcriptional response to TRs of the allopolyploid was more homeostatic compared with its parents. This better transcriptional homeostasis of the allopolyploid C. arabica afforded a greater phenotypic homeostasis when faced with environments that are unsuited to the diploid parental species.

Aquaporins in Coffea arabica L.: Identification, expression, and impacts on plant water relations and hydraulics

Plant aquaporins (AQPs) are involved in the transport of water and other small solutes across cell membranes, and thus play major roles in the regulation of plant water balance, as well as in growth regulation and response to abiotic stress factors. Limited information is currently available about the presence and role of AQPs in Coffea arabica L., despite the economic importance of the species and its vulnerability to drought stress. We identified candidate AQP genes by screening a proprietary C. arabica transcriptome database, resulting in the identification of nine putative aquaporins. A phylogenetic analysis based on previously characterized AQPs from Arabidopsis thaliana and Solanum tuberosum allowed to assign the putative coffee AQP sequences to the Tonoplast (TIP) and Plasma membrane (PIP) subfamilies. The possible functional role of coffee AQPs was explored by measuring hydraulic conductance and aquaporin gene expression on leaf and root tissues of two-year-old plants (C. arabica cv. Pacamara) subjected to different experimental conditions. In a first experiment, we tested plants for root and leaf hydraulic conductance both before dawn and at mid-day, to check the eventual impact of light on AQP activity and plant hydraulics. In a second experiment, we measured plant hydraulic responses to different water stress levels as eventually affected by changes in AQPs expression levels. Our results shed light on the possible roles of AQPs in the regulation of C. arabica hydraulics and water balance, opening promising research lines to improve the sustainability of coffee cultivation under global climate change scenarios.

Large distribution and high sequence identity of a Copia-type retrotransposon in angiosperm families

Retrotransposons are the main component of plant genomes. Recent studies have revealed the complexity of their evolutionary dynamics. Here, we have identified Copia25 in Coffea canephora, a new plant retrotransposon belonging to the Ty1-Copia superfamily. In the Coffea genomes analyzed, Copia25 is present in relatively low copy numbers and transcribed. Similarity sequence searches and PCR analyses show that this retrotransposon with LTRs (Long Terminal Repeats) is widely distributed among the Rubiaceae family and that it is also present in other distantly related species belonging to Asterids, Rosids and monocots. A particular situation is the high sequence identity found between the Copia25 sequences of Musa, a monocot, and Ixora, a dicot species (Rubiaceae). Our results reveal the complexity of the evolutionary dynamics of the ancient element Copia25 in angiosperm, involving several processes including sequence conservation, rapid turnover, stochastic losses and horizontal transfer.

New Insights on Coffea miRNAs: Features and Evolutionary Conservation

Small RNAs influence the gene expression at the post-transcriptional level by guiding messenger RNA (mRNA) cleavage, translational repression, and chromatin modifications. In addition to model plants, the microRNAs (miRNAs) have been identified in different crop species. In this work, we developed a specific pipeline to search for coffee miRNA homologs on expressed sequence tags (ESTs) and genome survey sequences (GSS) databases. As a result, 36 microRNAs were identified and a total of 616 and 362 potential targets for Coffea arabica and Coffea canephora, respectively. The evolutionary analyses of these molecules were performed by comparing the primary and secondary structures of precursors and mature miRNAs with their orthologs. Moreover, using a stem-loop RT-PCR assay, we evaluated the accumulation of mature miRNAs in genomes with different ploidy levels, detecting an increase in the miRNAs accumulation according to the ploidy raising. Finally, a 5' RACE (Rapid Amplification of cDNA Ends) assay was performed to verify the regulation of auxin responsive factor 8 (ARF8) by MIR167 in coffee plants. The great variety of target genes indicates the functional plasticity of these molecules and reinforces the importance of understanding the RNAi-dependent regulatory mechanisms. Our results expand the study of miRNAs and their target genes in this crop, providing new challenges to understand the biology of these species.

The Multi-Resistant Reaction of Drought-Tolerant Coffee 'Conilon Clone 14' to Meloidogyne spp. and Late Hypersensitive-Like Response in Coffea canephora

Root-knot nematodes (RKN), Meloidogyne spp., have major economic impact on coffee production in Central and South America. Genetic control of RKN constitutes an essential part for integrated pest management strategy. The objective of this study was to evaluate the resistance of Coffea canephora genotypes (clones) to Meloidogyne spp. Sensitive and drought-tolerant coffee genotypes were used to infer their resistance using nematode reproduction factor and histopathology. Eight clonal genotypes were highly resistant to M. paranaensis. 'Clone 14' (drought-tolerant) and 'ESN2010-04' were the only genotypes highly resistant and moderately resistant, respectively, to both M. incognita races 3 and 1. Several clones were highly resistant to both avirulent and virulent M. exigua. Clone 14 and ESN2010-04 showed multiple resistance to major RKNs tested. Roots of 'clone 14' (resistant) and 'clone 22' (susceptible) were histologically studied against infection by M. incognita race 3 and M. paranaensis. Reduction of juvenile (J2) penetration in clone 14 was first seen at 2 to 6 days after inoculation (DAI). Apparent early hypersensitive reaction (HR) was seen in root cortex between 4 and 6 DAI, which led to cell death and prevention of some nematode development. At 12 to 20 DAI, giant cells formed in the vascular cylinder, besides normal development into J3/J4. From 32 to 45 DAI, giant cells were completely degenerated. Late, intense HR and cell death were frequently observed around young females and giant cells reported for the first time in coffee pathosystem. These results provide rational bases for future studies, including prospection, characterization, and expression profiling of genomic loci involved in both drought tolerance and resistance to multiple RKN species.

A revised time tree of the asterids: establishing a temporal framework for evolutionary studies of the coffee family (rubiaceae)

Divergence time analyses in the coffee family (Rubiaceae) have all relied on the same Gentianales crown group age estimate, reported by an earlier analysis of the asterids, for defining the upper age bound of the root node in their analyses. However, not only did the asterid analysis suffer from several analytical shortcomings, but the estimate itself has been used in highly inconsistent ways in these Rubiaceae analyses. Based on the original data, we here reanalyze the divergence times of the asterids using relaxed-clock models and 14 fossil-based minimum age constraints. We also expand the data set to include an additional 67 taxa from Rubiaceae sampled across all three subfamilies recognized in the family. Three analyses are conducted: a separate analysis of the asterids, which completely mirrors the original asterid analysis in terms of taxon sample and data; a separate analysis of the Gentianales, where the result from the first analysis is used for defining a secondary root calibration point; and a combined analysis where all taxa are analyzed simultaneously. Results are presented in the form of a time-calibrated phylogeny, and age estimates for asterid groups, Gentianales, and major groups of Rubiaceae are compared and discussed in relation to previously published estimates. Our updated age estimates for major groups of Rubiaceae provide a significant step forward towards the long term goal of establishing a robust temporal framework for the divergence of this biologically diverse and fascinating group of plants.

Influence of coffee genotype on bioactive compounds and the in vitro capacity to scavenge reactive oxygen and nitrogen species

The influence of green coffee genotype on the bioactive compounds and the in vitro antioxidant capacity against the principal reactive oxygen (ROO(•), H2O2, HO(•), and HOCl) and nitrogen (NO(•) and ONOO(-)) species of biological relevance was investigated. This is the first report on the capacity of green coffee to scavenge H2O2, HOCl, and NO(•). Variations in the contents of total chlorogenic acids (22.9-37.9 g/100 g), cinnamoyl-amino acid conjugates (0.03-1.12 g/100 g), trigonelline (3.1-6.7 g/100 g), and caffeine (3.9-11.8 g/100 g) were found. Hydrophilic extracts of Coffea canephora and Coffea kapakata were the most potent scavengers of ROO(•), H2O2, HO(•), NO(•), and ONOO(-) due to their chlorogenic acid contents, which were, on average, 30% higher than those found in Coffea arabica and Coffea racemosa. The results showed that genotype is a determinant characteristic in the bioactive compound contents and consequently in the antioxidant capacity of green coffee.

Identification of two metallothioneins as novel inhalative coffee allergens cof a 2 and cof a 3

Background

Dust of green coffee beans is known to be a relevant cause for occupational allergic disorders in coffee industry workers. Recently, we described the first coffee allergen (Cof a 1) establishing an allergenic potential of green coffee dust.

Objective

Our aim was to identify allergenic components of green coffee in order to enhance inhalative coffee allergy diagnosis.

Methods

A Coffea arabica pJuFo cDNA phage display library was created and screened for IgE binding with sera from allergic coffee workers. Two further coffee allergens were identified by sequence analysis, expressed in E. coli, and evaluated by Western blots. The prevalence of sensitization to recombinant Cof a 1, Cof a 2, and Cof a 3 and to commercially available extract was investigated by ELISA (enzyme-linked immunosorbent assay) respectively CAP (capacity test) screening in 18 sera of symptomatic coffee workers.

Results

In addition to the previously described chitinase Cof a 1, two Coffea arabica cysteine-rich metallothioneins of 9 and 7 kDa were identified and included in the IUIS Allergen Nomenclature as Cof a 2 and Cof a 3. Serum IgE antibodies to at least one of the recombinant allergens were found in 8 out of 18 symptomatic coffee workers (44%). Only 2 of the analysed sera (11%) had reacted previously to the commercial allergy test.

Conclusions

In addition to the previously described Cof a 1 we have identified two further coffee proteins to be type I coffee allergens (Cof a 2 and Cof a 3) which may have a relevant potential for the specific diagnosis and/or therapy of coffee allergy.

Differential regulation of caffeine metabolism in Coffea arabica (Arabica) and Coffea canephora (Robusta)

Caffeine is a metabolite of great economic importance, especially in coffee, where it influences the sensorial and physiological impacts of the beverage. Caffeine metabolism in the Coffea species begins with the degradation of purine nucleotides through three specific N-methyltransferases: XMT, MXMT and DXMT. A comparative analysis was performed to clarify the molecular reasons behind differences in caffeine accumulation in two Coffea species, namely Coffea arabica and Coffea canephora var. robusta. Three different genes encoding N-methyltransferase were amplified in the doubled haploid Coffea canephora: CcXMT1, CcMXMT1 and CcDXMT. Six genes were amplified in the haploid Coffea arabica: CaXMT1, CaXMT2, CaMXMT1, CaMXMT2, CaDXMT1, and CaDXMT2. A complete phylogenic analysis was performed to identify specific key amino acids defining enzymatic function for each protein identified. Furthermore, a quantitative gene-expression analysis was conducted on leaves and on maturing coffee beans, simultaneously analyzing caffeine content. In the different varieties analyzed, caffeine accumulation is higher in leaves than in the coffee bean maturation period, higher in Robusta than in Arabica. In Robusta, CcXMT1 and CcDXMT gene expressions are predominant and transcriptional activity is higher in leaves than in maturing beans, and is highly correlated to caffeine accumulation. In Arabica, the CaXMT1 expression level is high in leaves and CaDXMT2 as well to a lesser extent, while global transcriptional activity is weak during bean maturation, suggesting that the transcriptional control of caffeine-related genes differs within different organs and between Arabica and Robusta. These findings indicate that caffeine accumulation in Coffea species has been modulated by a combination of differential transcriptional regulation and genome evolution.

The coffee genome hub: a resource for coffee genomes

The whole genome sequence of Coffea canephora, the perennial diploid species known as Robusta, has been recently released. In the context of the C. canephora genome sequencing project and to support post-genomics efforts, we developed the Coffee Genome Hub (http://coffee-genome.org/), an integrative genome information system that allows centralized access to genomics and genetics data and analysis tools to facilitate translational and applied research in coffee. We provide the complete genome sequence of C. canephora along with gene structure, gene product information, metabolism, gene families, transcriptomics, syntenic blocks, genetic markers and genetic maps. The hub relies on generic software (e.g. GMOD tools) for easy querying, visualizing and downloading research data. It includes a Genome Browser enhanced by a Community Annotation System, enabling the improvement of automatic gene annotation through an annotation editor. In addition, the hub aims at developing interoperability among other existing South Green tools managing coffee data (phylogenomics resources, SNPs) and/or supporting data analyses with the Galaxy workflow manager.

Flower development in Coffea arabica L.: new insights into MADS-box genes

Coffea arabica L. shows peculiar characteristics during reproductive development, such as flowering asynchrony, periods of floral bud dormancy, mucilage secretion and epipetalous stamens. The MADS-box transcription factors are known to control several developmental processes in plants, including flower and fruit development. Significant differences are found among plant species regarding reproductive development and little is known about the role of MADS-box genes in Coffea reproductive development. Thus, we used anatomical and comparative molecular analyses to explore the flowering process in coffee. The main morphological changes during flower development in coffee were observed by optical and scanning electron microscopy. Flowering asynchrony seems to be related to two independent processes: the asynchronous development of distinct buds before the reproductive induction and the asynchronous development of floral meristems within each bud after the reproductive induction. A total of 23 C. arabica MADS-box genes were characterized by sequence comparison with putative Arabidopsis orthologs and their expression profiles were analyzed by RT-PCR in different tissues. The expression of the ABC model orthologs in Coffea during floral development was determined by in situ hybridization. The APETALA1 (AP1) ortholog is expressed only late in the perianth, which is also observed for the APETALA3 and TM6 orthologs. Conversely, the PISTILLATA ortholog is widely expressed in early stages, but restrict to stamens and carpels in later stages of flower development, while the expression of the AGAMOUS ortholog is always restricted to fertile organs. The AP1 and PISTILLATA orthologs are also expressed at specific floral organs, such as bracts and colleters, respectively, suggesting a potential role in the development of such structures. Altogether, the results from our comprehensive expression analyses showed significant differences between the spatiotemporal expression profiles of C. arabica MADS-box genes and their orthologs, which suggests differential functionalization in coffee. Moreover, these differences might also partially explain the particular characteristics of floral development in coffee, such as mucilage secretion and formation of epipetalous stamens.

Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next-generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single-nucleotide polymorphisms detected using RNA-seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.

Regulation of galactomannan biosynthesis in coffee seeds

The seed of Coffea arabica accumulates large amounts of cell wall storage polysaccharides (CWSPs) of the mannan family in the cell walls of the endosperm. The variability induced by the growing environment and extensive pairwise correlation analysis with stringent significance thresholds was used to investigate transcript-transcript and transcript-metabolite relationships among 26 sugar-related genes, and the amount of CWSPs and seven soluble low molecular weight carbohydrates in the developing coffee endosperm. A dense module of nine quantitatively co-expressed genes was detected at the mid-developmental stage when CWSPs accumulate. This module included the five genes of the core galactomannan synthetic machinery, namely genes coding for the enzymes needed to assemble the mannan backbone (mannan synthase, ManS), and genes that introduce the galactosyl side chains (galactosyltransferase, GMGT), modulate the post-depositional degree of galactose substitution (α-galactosidase), and produce the nucleotide sugar building blocks GDP-mannose and UDP-galactose (mannose-1P guanyltransferase and UDP-glucose 4'-epimerase, respectively). The amount of CWSPs stored in the endosperm at the onset of their accumulation was primarily and quantitatively modulated at the transcriptional level (i.e. positively correlated with the expression level of these key galactomannan biosynthetic genes). This analysis also suggests a role for sorbitol and raffinose family oligosaccharides as transient auxiliary sources of building blocks for galactomannan synthesis. Finally, a microarray-based analysis of the developing seed transcriptome revealed that all genes of the core galactomannan synthesis machinery grouped in a single cluster of 209 co-expressed genes. Analysis of the gene composition of this cluster revealed remarkable functional coherence and identified transcription factors that putatively control galactomannan biosynthesis in coffee.

Dynamics of the concentration of IAA and some of its conjugates during the induction of somatic embryogenesis in Coffea canephora

Most of the somatic embryogenesis (SE) process requires the presence, either before or during the embryogenic process, of at least one exogenous auxin. This exogenous auxin induces the presence of endogenous auxins, which appears to be essential for SE induction. We found that during the preincubation period of SE in Coffea canephora, there is an important increase in both free and conjugated indole-3-acetic acid (IAA), as well as indole-3-butyric acid. This increase is accompanied by an increase in the expression of YUCCA (CcYUC), TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS 1 (CcTAA1), and GRETCHEN HAGEN 3 (GH3) genes. On the other hand, most of the IAA compounds decreased during the induction of SE. The results presented in this research suggest that a balance between free IAA and its amide conjugates is necessary to allow the expression of SE-related genes.

Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures

Polyploidy has occurred throughout the evolutionary history of plants and led to diversification and plant ecological adaptation. Functional plasticity of duplicate genes is believed to play a major role in the environmental adaptation of polyploids. In this context, we characterized genome-wide homoeologous gene expression in Coffea arabica, a recent allopolyploid combining two subgenomes that derive from two closely related diploid species, and investigated its variation in response to changing environment. The transcriptome of leaves of C. arabica cultivated at different growing temperatures suitable for one or the other parental species was examined using RNA-sequencing. The relative contribution of homoeologs to gene expression was estimated for 9959 and 10,628 genes in warm and cold conditions, respectively. Whatever the growing conditions, 65% of the genes showed equivalent levels of homoeologous gene expression. In 92% of the genes, relative homoeologous gene expression varied < 10% between growing temperatures. The subgenome contributions to the transcriptome appeared to be only marginally altered by the different conditions (involving intertwined regulations of homeologs) suggesting that C. arabica's ability to tolerate a broader range of growing temperatures than its diploid parents does not result from differential use of homoeologs.

New insights into somatic embryogenesis: leafy cotyledon1, baby boom1 and WUSCHEL-related homeobox4 are epigenetically regulated in Coffea canephora

Plant cells have the capacity to generate a new plant without egg fertilization by a process known as somatic embryogenesis (SE), in which differentiated somatic cells can form somatic embryos able to generate a functional plant. Although there have been advances in understanding the genetic basis of SE, the epigenetic mechanism that regulates this process is still unknown. Here, we show that the embryogenic development of Coffea canephora proceeds through a crosstalk between DNA methylation and histone modifications during the earliest embryogenic stages of SE. We found that low levels of DNA methylation, histone H3 lysine 9 dimethylation (H3K9me2) and H3K27me3 change according to embryo development. Moreover, the expression of LEAFY cotyledon1 (LEC1) and BABY BOOM1 (BBM1) are only observed after SE induction, whereas WUSCHEL-related homeobox4 (WOX4) decreases its expression during embryo maturation. Using a pharmacological approach, it was found that 5-Azacytidine strongly inhibits the embryogenic response by decreasing both DNA methylation and gene expression of LEC1 and BBM1. Therefore, in order to know whether these genes were epigenetically regulated, we used Chromatin Immunoprecipitation (ChIP) assays. It was found that WOX4 is regulated by the repressive mark H3K9me2, while LEC1 and BBM1 are epigenetically regulated by H3K27me3. We conclude that epigenetic regulation plays an important role during somatic embryogenic development, and a molecular mechanism for SE is proposed.

Promoter analysis of the WRKY transcription factors CaWRKY1a and CaWRKY1b homoeologous genes in coffee (Coffea arabica)

KEY MESSAGE

The regulation of the CaWRKY1 homoeologous genes were analyzed through the characterization of their promoters. The pW1a promoter is proposed as a new tool for coffee plant biotechnologies. WRKY transcription factors are important elements of the plant immune response. The CaWRKY1 gene from Coffea arabica is induced by several biotic and abiotic stresses, including challenge by the rust fungus Hemileia vastatrix. Two homoeologous CaWRKY1 genes, named CaWRKY1a and CaWRKY1b, were previously identified in the C. arabica allotetraploid genome. To gain insight into the transcriptional regulation of these genes, their promoter sequences, named pW1a and pW1b, respectively, were cloned and characterized in this study. In silico analysis revealed some important defense-associated regulatory elements, including W-boxes and as-1 elements. Promoter activities were analyzed in transient assays conducted by agroinfiltration of tobacco leaves. Exogenous salicylic acid (SA) treatments increased promoter activities corroborating the presence of as-1 regulatory elements. Transactivation assays with the CaWRKY1 protein showed the reduction of both pW1a and pW1b promoter activities, indicating that the CaWRKY1 protein may negatively regulate its own promoters. Stable transgenic C. arabica lines expressing a pW1a::GUS construct were obtained by Agrobacterium-mediated transformation and high GUS activity was observed in leaves subjected to mechanical wounding. Hence, the ability of pW1a to drive transgene expression in coffee plants as well as to enhance expression in response to stresses opens possibilities for using this promoter as a new tool for biotechnological approaches in coffee plants.

Functional characterization of three Coffea arabica L. monoterpene synthases: insights into the enzymatic machinery of coffee aroma

The chemical composition of the coffee beverage is extremely complex, being made up of hundreds of volatile and non-volatile compounds, many of which are generated in the thermal reactions that occur during the roasting process. However, in the raw coffee bean there are also compounds that survive roasting and are therefore extracted into the beverage. Monoterpenes are an example of this category, as their presence has been reported in the coffee flower, fruit, seed, roasted bean and in the beverage aroma. The present work describes the isolation, heterologous expression and functional characterization of three Coffea arabica cDNAs coding for monoterpene synthases. RNA was purified from C. arabica (cv. Catuai Red) flowers, seeds and fruits at 4 successive ripening stages. Degenerate primers were designed on the most conserved regions of the monoterpene synthase gene family, and then used to isolate monoterpene synthase-like sequences from the cDNA libraries. After 5'- and 3'-RACE, the complete transcripts of 4 putative C. arabica monoterpene synthases (CofarTPS) were obtained. Gene expression in different tissues and developmental stages was analysed. After heterologous expression in Escherichia coli, enzyme activity and substrate specificity were evaluated in vitro by incubation of the recombinant proteins with geranyl pyrophosphate (GPP), geranylgeranyl pyrophosphate (GGPP) and farnesyl pyrophosphate (FPP), precursors respectively of mono-, di- and sesquiterpenes. The reaction products were characterized by HS-SPME GC-MS. CofarTPS1 was classified as a limonene synthase gene, while CofarTPS2 and 3 showed lower activity with the production of linalool and β-myrcene.

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

Embryogenic suspensions that involve extensive cell division are risky in respect to genome and epigenome instability. Elevated frequencies of somaclonal variation in embryogenic suspension-derived plants were reported in many species, including coffee. This problem could be overcome by using culture conditions that allow moderate cell proliferation. In view of true-to-type large-scale propagation of C. arabica hybrids, suspension protocols based on low 2,4-D concentrations and short proliferation periods were developed. As mechanisms leading to somaclonal variation are often complex, the phenotypic, genetic and epigenetic changes were jointly assessed so as to accurately evaluate the conformity of suspension-derived plants. The effects of embryogenic suspensions and secondary embryogenesis, used as proliferation systems, on the genetic conformity of somatic embryogenesis-derived plants (emblings) were assessed in two hybrids. When applied over a 6 month period, both systems ensured very low somaclonal variation rates, as observed through massive phenotypic observations in field plots (0.74% from 200,000 plant). Molecular AFLP and MSAP analyses performed on 145 three year-old emblings showed that polymorphism between mother plants and emblings was extremely low, i.e. ranges of 0-0.003% and 0.07-0.18% respectively, with no significant difference between the proliferation systems for the two hybrids. No embling was found to cumulate more than three methylation polymorphisms. No relation was established between the variant phenotype (27 variants studied) and a particular MSAP pattern. Chromosome counting showed that 7 of the 11 variant emblings analyzed were characterized by the loss of 1-3 chromosomes. This work showed that both embryogenic suspensions and secondary embryogenesis are reliable for true-to-type propagation of elite material. Molecular analyses revealed that genetic and epigenetic alterations are particularly limited during coffee somatic embryogenesis. The main change in most of the rare phenotypic variants was aneuploidy, indicating that mitotic aberrations play a major role in somaclonal variation in coffee.

Genetic structure and diversity of coffee (Coffea) across Africa and the Indian Ocean islands revealed using microsatellites

BACKGROUND AND AIMS

The coffee genus (Coffea) comprises 124 species, and is indigenous to the Old World Tropics. Due to its immense economic importance, Coffea has been the focus of numerous genetic diversity studies, but despite this effort it remains insufficiently studied. In this study the genetic diversity and genetic structure of Coffea across Africa and the Indian Ocean islands is investigated.

METHODS

Genetic data were produced using 13 polymorphic nuclear microsatellite markers (simple sequence repeats, SSRs), including seven expressed sequence tag-SSRs, and the data were analysed using model- and non-model-based methods. The study includes a total of 728 individuals from 60 species.

KEY RESULTS

Across Africa and the Indian Ocean islands Coffea comprises a closely related group of species with an overall pattern of genotypes running from west to east. Genetic structure was identified in accordance with pre-determined geographical regions and phylogenetic groups. There is a good relationship between morpho-taxonomic species delimitations and genetic units. Genetic diversity in African and Indian Ocean Coffea is high in terms of number of alleles detected, and Madagascar appears to represent a place of significant diversification in terms of allelic richness and species diversity.

CONCLUSIONS

Cross-species SSR transferability in African and Indian Ocean islands Coffea was very efficient. On the basis of the number of private alleles, diversification in East Africa and the Indian Ocean islands appears to be more recent than in West and West-Central Africa, although this general trend is complicated in Africa by the position of species belonging to lineages connecting the main geographical regions. The general pattern of phylogeography is not in agreement with an overall east to west (Mascarene, Madagascar, East Africa, West Africa) increase in genome size, the high proportion of shared alleles between the four regions or the high numbers of exclusive shared alleles between pairs or triplets of regions.

A structural basis for the biosynthesis of the major chlorogenic acids found in coffee

Chlorogenic acids (CGAs) are a group of phenolic secondary metabolites produced by certain plant species and an important component of coffee (Coffea spp.). The CGAs have been implicated in biotic and abiotic stress responses, while the related shikimate esters are key intermediates for lignin biosynthesis. Here, two hydroxycinnamoyl-coenzyme A shikimate/quinate hydroxycinnamoyl transferases (HCT/HQT) from coffee were biochemically characterized. We show, to our knowledge for the first time, that in vitro, HCT is capable of synthesizing the 3,5-O-dicaffeoylquinic acid diester, a major constituent of the immature coffee grain. In order to further understand the substrate specificity and catalytic mechanism of the HCT/HQT, we performed structural and mutagenesis studies of HCT. The three-dimensional structure of a native HCT and a proteolytically stable lysine mutant enabled the identification of important residues involved in substrate specificity and catalysis. Site-directed mutagenesis confirmed the role of residues leucine-400 and phenylalanine-402 in substrate specificity and of histidine-153 and the valine-31 to proline-37 loop in catalysis. In addition, the histidine-154-asparagine mutant was observed to produce 4-fold more dichlorogenic acids compared with the native protein. These data provide, to our knowledge, the first structural characterization of a HCT and, in conjunction with the biochemical and mutagenesis studies presented here, delineate the underlying molecular-level determinants for substrate specificity and catalysis. This work has potential applications in fine-tuning the levels of shikimate and quinate esters (CGAs including dichlorogenic acids) in different plant species in order to generate reduced or elevated levels of the desired target compounds.

CaPrx, a Coffea arabica gene encoding a putative class III peroxidase induced by root-knot nematode infection

Class III peroxidases (Prxs) are enzymes involved in a multitude of physiological and stress-related processes in plants. Here, we report on the characterization of a putative peroxidase-encoding gene from Coffea arabica (CaPrx) that is expressed in early stages of root-knot nematode (RKN) infection. CaPrx showed enhanced expression in coffee roots inoculated with RKN (at 12 h post-inoculation), but no significant difference in expression was observed between susceptible and resistant plants. Assays using transgenic tobacco plants harboring a promoter-β-glucuronidase (GUS) fusion revealed that the CaPrx promoter was exclusively active in the galls induced by RKN. In cross sections of galls, GUS staining was predominantly localized in giant cells. Up-regulation of GUS expression in roots of transgenic plants following RKN inoculation was observed within 16 h. Moreover, no increase in GUS expression after treatment with jasmonic acid was detected. Altogether, these results point to a putative role of this peroxidase in the general coffee response to RKN infection.

Characterization, high-resolution mapping and differential expression of three homologous PAL genes in Coffea canephora Pierre (Rubiaceae)

Phenylalanine ammonia lyase (PAL) is the first entry enzyme of the phenylpropanoid pathway producing phenolics, widespread constituents of plant foods and beverages, including chlorogenic acids, polyphenols found at remarkably high levels in the coffee bean and long recognized as powerful antioxidants. To date, whereas PAL is generally encoded by a small gene family, only one gene has been characterized in Coffea canephora (CcPAL1), an economically important species of cultivated coffee. In this study, a molecular- and bioinformatic-based search for CcPAL1 paralogues resulted successfully in identifying two additional genes, CcPAL2 and CcPAL3, presenting similar genomic structures and encoding proteins with close sequences. Genetic mapping helped position each gene in three different coffee linkage groups, CcPAL2 in particular, located in a coffee genome linkage group (F) which is syntenic to a region of Tomato Chromosome 9 containing a PAL gene. These results, combined with a phylogenetic study, strongly suggest that CcPAL2 may be the ancestral gene of C. canephora. A quantitative gene expression analysis was also conducted in coffee tissues, showing that all genes are transcriptionally active, but they present distinct expression levels and patterns. We discovered that CcPAL2 transcripts appeared predominantly in flower, fruit pericarp and vegetative/lignifying tissues like roots and branches, whereas CcPAL1 and CcPAL3 were highly expressed in immature fruit. This is the first comprehensive study dedicated to PAL gene family characterization in coffee, allowing us to advance functional studies which are indispensable to learning to decipher what role this family plays in channeling the metabolism of coffee phenylpropanoids.

Differentially expressed genes and proteins upon drought acclimation in tolerant and sensitive genotypes of Coffea canephora

The aim of this study was to investigate the molecular mechanisms underlying drought acclimation in coffee plants by the identification of candidate genes (CGs) using different approaches. The first approach used the data generated during the Brazilian Coffee expressed sequence tag (EST) project to select 13 CGs by an in silico analysis (electronic northern). The second approach was based on screening macroarrays spotted with plasmid DNA (coffee ESTs) with separate hybridizations using leaf cDNA probes from drought-tolerant and susceptible clones of Coffea canephora var. Conilon, grown under different water regimes. This allowed the isolation of seven additional CGs. The third approach used two-dimensional gel electrophoresis to identify proteins displaying differential accumulation in leaves of drought-tolerant and susceptible clones of C. canephora. Six of them were characterized by MALDI-TOF-MS/MS (matrix-assisted laser desorption-time of flight-tandem mass spectrometry) and the corresponding proteins were identified. Finally, additional CGs were selected from the literature, and quantitative real-time polymerase chain reaction (qPCR) was performed to analyse the expression of all identified CGs. Altogether, >40 genes presenting differential gene expression during drought acclimation were identified, some of them showing different expression profiles between drought-tolerant and susceptible clones. Based on the obtained results, it can be concluded that factors involved a complex network of responses probably involving the abscisic signalling pathway and nitric oxide are major molecular determinants that might explain the better efficiency in controlling stomata closure and transpiration displayed by drought-tolerant clones of C. canephora.

Quantitative bioluminometric method for DNA-based species/varietal identification in food authenticity assessment

A method is reported for species quantification by exploiting single-nucleotide polymorphisms (SNPs). These single-base changes in DNA are particularly useful because they enable discrimination of closely related species and/or varieties. As a model, quantitative authentication studies were performed on coffee. These involved the determination of the percentage of Arabica and Robusta species based on a SNP in the chloroplastic trnL(UAA)-trnF(GAA) intraspacer region. Following polymerase chain reaction (PCR), the Robusta-specific and Arabica-specific fragments were subjected to 15 min extension reactions by DNA polymerase using species-specific primers carrying oligo(dA) tags. Biotin was incorporated into the extended strands. The products were captured in streptavidin-coated microtiter wells and quantified by using oligo(dT)-conjugated photoprotein aequorin. Aequorin was measured within 3 s via its characteristic flash-type bioluminescent reaction that was triggered by the addition of Ca(2+). Because of the close resemblance between the two DNA fragments, during PCR one species serves as an internal standard for the other. The percentage of the total luminescence signal obtained from a certain species was linearly related to the percent content of the sample with respect to this species. The method is accurate and reproducible. The microtiter well-based assay configuration allows high sample throughput and facilitates greatly the automation.

Dipstick test for DNA-based food authentication. Application to coffee authenticity assessment

This paper reports DNA-based food authenticity assays, in which species identification is accomplished by the naked eye without the need of specialized instruments. Strongly colored nanoparticles (gold nanoparticles) are employed as reporters that enable visual detection. Furthermore, detection is performed in a low-cost, disposable, dipstick-type device that incorporates the required reagents in dry form, thereby avoiding multiple pipetting and incubation steps. Due to its simplicity, the method does not require highly qualified personnel. The procedure comprises the following steps: (i) PCR amplification of the DNA segment that flanks the unique SNP (species marker); (ii) a 15 min extension reaction in which DNA polymerase extends an allele-specific primer only if it is perfectly complementary with the target sequence; (iii) detection of the products of the extension reaction within a few minutes by the naked eye employing the dipstick. No purification is required prior to application of the extension products to the dipstick. The method is general and requires only a unique DNA sequence for species discrimination. The only instrument needed is a conventional thermocycler for PCR, which is common equipment in every DNA laboratory. As a model, the method was applied to the discrimination of Coffea robusta and arabica species in coffee authenticity assessment. As low as 5% of Robusta coffee can be detected in the presence of Arabica coffee.

454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction

Coffee (Coffea arabica L.), one of the key export and cash crops in tropical and subtropical countries, suffers severe losses from the rust fungus Hemileia vastatrix. The transcriptome of H. vastatrix was analysed during a compatible interaction with coffee to obtain an exhaustive repertoire of the genes expressed during infection and to identify potential effector genes. Large-scale sequencing (454-GS-FLEX Titanium) of mixed coffee and rust cDNAs obtained from 21-day rust-infected leaves generated 352 146 sequences which assembled into 22 774 contigs. In the absence of any reference genomic sequences for Coffea or Hemileia, specific trinucleotide frequencies within expressed sequence tags (ESTs) and blast homology against a set of dicots and basidiomycete genomes were used to distinguish pathogen from plant sequences. About 30% (6763) of the contigs were assigned to H. vastatrix and 61% (13 951) to C. arabica. The majority (60%) of the rust sequences did not show homology to any genomic database, indicating that they were potential novel fungal genes. In silico analyses of the 6763 H. vastatrix contigs predicted 382 secreted proteins and identified homologues of the flax rust haustorially expressed secreted proteins (HESPs) and bean rust transferred protein 1 (RTP1). These rust candidate effectors showed conserved amino-acid domains and conserved patterns of cysteine positions suggestive of conserved functions during infection of host plants. Quantitative reverse transcription-polymerase chain reaction profiling of selected rust genes revealed dynamic expression patterns during the time course of infection of coffee leaves. This study provides the first valuable genomic resource for the agriculturally important plant pathogen H. vastatrix and the first comprehensive C. arabica EST dataset.