Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.

Cryptosexuality and the genetic diversity paradox in coffee rust, Hemileia vastatrix

BACKGROUND

Despite the fact that coffee rust was first investigated scientifically more than a century ago, and that the disease is one of the major constraints to coffee production--constantly changing the socio-economic and historical landscape of the crop--critical aspects of the life cycle of the pathogen, Hemileia vastatrix, remain unclear. The asexual urediniospores are regarded as the only functional propagule: theoretically, making H. vastatrix a clonal species. However, the well-documented emergence of new rust pathotypes and the breakdown in genetic resistance of coffee cultivars, present a paradox.

METHODS AND RESULTS

Here, using computer-assisted DNA image cytometry, following a modified nuclear stoichiometric staining technique with Feulgen, we show that meiosis occurs within the urediniospores. Stages of spore development were categorised based on morphology, from the spore-mother cell through to the germinating spore, and the relative nuclear DNA content was quantified statistically at each stage.

CONCLUSIONS

Hidden sexual reproduction disguised within the asexual spore (cryptosexuality) could explain why new physiological races have arisen so often and so quickly in Hemileia vastatrix. This could have considerable implications for coffee breeding strategies and may be a common event in rust fungi, especially in related genera occupying the same basal phylogenetic lineages.

Genomic expression dominance in the natural allopolyploid Coffea arabica is massively affected by growth temperature

• Polyploidy occurs throughout the evolutionary history of many plants and considerably impacts species diversity, giving rise to novel phenotypes and leading to ecological diversification and colonization of new niches. Recent studies have documented dynamic changes in plant polyploid gene expression, which reflect the genomic and functional plasticity of duplicate genes and genomes. • The aim of the present study was to describe genomic expression dominance between a relatively recently formed natural allopolyploid (Coffea arabica) and its ancestral parents (Coffea canephora and Coffea eugenioides) and to determine if the divergence was environment-dependent. Employing a microarray platform designed against 15,522 unigenes, we assayed unigene expression levels in the allopolyploid and its two parental diploids. For each unigene, we measured expression variations among the three species grown under two temperature conditions (26-22°C (day-night temperatures) and 30-26°C (day-night temperatures)). • More than 35% of unigenes were differentially expressed in each comparison at both temperatures, except for C. arabica vs C. canephora in the 30-26°C range, where an unexpectedly low unigene expression divergence (< 9%) was observed. • Our data revealed evidence of transcription profile divergence between the allopolyploid and its parental species, greatly affected by environmental conditions, and provide clues to the plasticity phenomenon in allopolyploids.

Identification and characterization of the Non-race specific Disease Resistance 1 (NDR1) orthologous protein in coffee

BACKGROUND

Leaf rust, which is caused by the fungus Hemileia vastatrix (Pucciniales), is a devastating disease that affects coffee plants (Coffea arabica L.). Disadvantages that are associated with currently developed phytoprotection approaches have recently led to the search for alternative strategies. These include genetic manipulations that constitutively activate disease resistance signaling pathways. However, molecular actors of such pathways still remain unknown in C. arabica. In this study, we have isolated and characterized the coffee NDR1 gene, whose Arabidopsis ortholog is a well-known master regulator of the hypersensitive response that is dependent on coiled-coil type R-proteins.

RESULTS

Two highly homologous cDNAs coding for putative NDR1 proteins were identified and cloned from leaves of coffee plants. One of the candidate coding sequences was then expressed in the Arabidopsis knock-out null mutant ndr1-1. Upon a challenge with a specific strain of the bacterium Pseudomonas syringae (DC3000::AvrRpt2), analysis of both macroscopic symptoms and in planta microbial growth showed that the coffee cDNA was able to restore the resistance phenotype in the mutant genetic background. Thus, the cDNA was dubbed CaNDR1a (standing for Coffea arabica Non-race specific Disease Resistance 1a). Finally, biochemical and microscopy data were obtained that strongly suggest the mechanistic conservation of the NDR1-driven function within coffee and Arabidopsis plants. Using a transient expression system, it was indeed shown that the CaNDR1a protein, like its Arabidopsis counterpart, is localized to the plasma membrane, where it is possibly tethered by means of a GPI anchor.

CONCLUSIONS

Our data provide molecular and genetic evidence for the identification of a novel functional NDR1 homolog in plants. As a key regulator initiating hypersensitive signalling pathways, CaNDR1 gene(s) might be target(s) of choice for manipulating the coffee innate immune system and achieving broad spectrum resistance to pathogens. Given the potential conservation of NDR1-dependent defense mechanisms between Arabidopsis and coffee plants, our work also suggests new ways to isolate the as-yet-unidentified R-gene(s) responsible for resistance to H. vastatrix.

Micro-collinearity and genome evolution in the vicinity of an ethylene receptor gene of cultivated diploid and allotetraploid coffee species (Coffea)

Arabica coffee (Coffea arabica L.) is a self-compatible perennial allotetraploid species (2n=4x=44), whereas Robusta coffee (C. canephora L.) is a self-incompatible perennial diploid species (2n=2x=22). C. arabica (C(a) C(a) E(a) E(a) ) is derived from a spontaneous hybridization between two closely related diploid coffee species, C. canephora (CC) and C. eugenioides (EE). To investigate the patterns and degree of DNA sequence divergence between the Arabica and Robusta coffee genomes, we identified orthologous bacterial artificial chromosomes (BACs) from C. arabica and C. canephora, and compared their sequences to trace their evolutionary history. Although a high level of sequence similarity was found between BACs from C. arabica and C. canephora, numerous chromosomal rearrangements were detected, including inversions, deletions and insertions. DNA sequence identity between C. arabica and C. canephora orthologous BACs ranged from 93.4% (between E(a) and C(a) ) to 94.6% (between C(a) and C). Analysis of eight orthologous gene pairs resulted in estimated ages of divergence between 0.046 and 0.665 million years, indicating a recent origin of the allotetraploid species C. arabica. Analysis of transposable elements revealed differential insertion events that contributed to the size increase in the C(a) sub-genome compared to its diploid relative. In particular, we showed that insertion of a Ty1-copia LTR retrotransposon occurred specifically in C. arabica, probably shortly after allopolyploid formation. The two sub-genomes of C. arabica, C(a) and E(a) , showed sufficient sequence differences, and a whole-genome shotgun approach could be suitable for sequencing the allotetraploid genome of C. arabica.

Expression and trans-specific polymorphism of self-incompatibility RNases in coffea (Rubiaceae)

Self-incompatibility (SI) is widespread in the angiosperms, but identifying the biochemical components of SI mechanisms has proven to be difficult in most lineages. Coffea (coffee; Rubiaceae) is a genus of old-world tropical understory trees in which the vast majority of diploid species utilize a mechanism of gametophytic self-incompatibility (GSI). The S-RNase GSI system was one of the first SI mechanisms to be biochemically characterized, and likely represents the ancestral Eudicot condition as evidenced by its functional characterization in both asterid (Solanaceae, Plantaginaceae) and rosid (Rosaceae) lineages. The S-RNase GSI mechanism employs the activity of class III RNase T2 proteins to terminate the growth of "self" pollen tubes. Here, we investigate the mechanism of Coffea GSI and specifically examine the potential for homology to S-RNase GSI by sequencing class III RNase T2 genes in populations of 14 African and Madagascan Coffea species and the closely related self-compatible species Psilanthus ebracteolatus. Phylogenetic analyses of these sequences aligned to a diverse sample of plant RNase T2 genes show that the Coffea genome contains at least three class III RNase T2 genes. Patterns of tissue-specific gene expression identify one of these RNase T2 genes as the putative Coffea S-RNase gene. We show that populations of SI Coffea are remarkably polymorphic for putative S-RNase alleles, and exhibit a persistent pattern of trans-specific polymorphism characteristic of all S-RNase genes previously isolated from GSI Eudicot lineages. We thus conclude that Coffea GSI is most likely homologous to the classic Eudicot S-RNase system, which was retained since the divergence of the Rubiaceae lineage from an ancient SI Eudicot ancestor, nearly 90 million years ago.

Genetic molecular analysis of Coffea arabica (Rubiaceae) hybrids using SRAP markers

In Coffea arabica (arabica coffee), the phenotypic as well as genetic variability has been found low because of the narrow genetic basis and self fertile nature of the species. Because of high similarity in phenotypic appearance among the majority of arabica collections, selection of parental lines for inter-varietals hybridization and identification of resultant hybrids at an early stage of plant growth is difficult. DNA markers are known to be reliable in identifying closely related cultivars and hybrids. Sequence Related Amplified Polymorphism (SRAP) is a new molecular marker technology developed based on PCR. In this paper, sixty arabica-hybrid progenies belonging to six crosses were analyzed using 31 highly polymorphic SRAP markers. The analysis revealed seven types of SRAP marker profiles which are useful in discriminating the parents and hybrids. The number of bands amplified per primer pair ranges from 6.13 to 8.58 with average number of seven bands. Among six hybrid combinations, percentage of bands shared between hybrids and their parents ranged from 66.29% to 85.71% with polymorphic bands varied from 27.64% to 60.0%. Percentage of hybrid specific fragments obtained in various hybrid combinations ranged from 0.71% to 10.86% and ascribed to the consequence of meiotic recombination. Based on the similarity index calculation, it was observed that F1 hybrids share maximum number of bands with the female parent compared to male parent. The results obtained in the present study revealed the effectiveness of SRAP technique in cultivar identification and hybrid analysis in this coffee species.

Agrobacterium-mediated genetic transformation of Coffea arabica (L.) is greatly enhanced by using established embryogenic callus cultures

BACKGROUND

Following genome sequencing of crop plants, one of the main challenges today is determining the function of all the predicted genes. When gene validation approaches are used for woody species, the main obstacle is the low recovery rate of transgenic plants from elite or commercial cultivars. Embryogenic calli have frequently been the target tissue for transformation, but the difficulty in producing or maintaining embryogenic tissues is one of the main problems encountered in genetic transformation of many woody plants, including Coffea arabica.

RESULTS

We identified the conditions required for successful long-term proliferation of embryogenic cultures in C. arabica and designed a highly efficient and reliable Agrobacterium tumefaciens-mediated transformation method based on these conditions. The transformation protocol with LBA1119 harboring pBin 35S GFP was established by evaluating the effect of different parameters on transformation efficiency by GFP detection. Using embryogenic callus cultures, co-cultivation with LBA1119 OD600 = 0.6 for five days at 20 °C enabled reproducible transformation. The maintenance conditions for the embryogenic callus cultures, particularly a high auxin to cytokinin ratio, the age of the culture (optimum for 7-10 months of proliferation) and the use of a yellow callus phenotype, were the most important factors for achieving highly efficient transformation (> 90%). At the histological level, successful transformation was related to the number of proembryogenic masses present. All the selected plants were proved to be transformed by PCR and Southern blot hybridization.

CONCLUSION

Most progress in increasing transformation efficiency in coffee has been achieved by optimizing the production conditions of embryogenic cultures used as target tissues for transformation. This is the first time that a strong positive effect of the age of the culture on transformation efficiency was demonstrated. Our results make Agrobacterium-mediated transformation of embryogenic cultures a viable and useful tool both for coffee breeding and for the functional analysis of agronomically important genes.

RBCS1 expression in coffee: Coffea orthologs, Coffea arabica homeologs, and expression variability between genotypes and under drought stress

BACKGROUND

In higher plants, the inhibition of photosynthetic capacity under drought is attributable to stomatal and non-stomatal (i.e., photochemical and biochemical) effects. In particular, a disruption of photosynthetic metabolism and Rubisco regulation can be observed. Several studies reported reduced expression of the RBCS genes, which encode the Rubisco small subunit, under water stress.

RESULTS

Expression of the RBCS1 gene was analysed in the allopolyploid context of C. arabica, which originates from a natural cross between the C. canephora and C. eugenioides species. Our study revealed the existence of two homeologous RBCS1 genes in C. arabica: one carried by the C. canephora sub-genome (called CaCc) and the other carried by the C. eugenioides sub-genome (called CaCe). Using specific primer pairs for each homeolog, expression studies revealed that CaCe was expressed in C. eugenioides and C. arabica but was undetectable in C. canephora. On the other hand, CaCc was expressed in C. canephora but almost completely silenced in non-introgressed ("pure") genotypes of C. arabica. However, enhanced CaCc expression was observed in most C. arabica cultivars with introgressed C. canephora genome. In addition, total RBCS1 expression was higher for C. arabica cultivars that had recently introgressed C. canephora genome than for "pure" cultivars. For both species, water stress led to an important decrease in the abundance of RBCS1 transcripts. This was observed for plants grown in either greenhouse or field conditions under severe or moderate drought. However, this reduction of RBCS1 gene expression was not accompanied by a decrease in the corresponding protein in the leaves of C. canephora subjected to water withdrawal. In that case, the amount of RBCS1 was even higher under drought than under unstressed (irrigated) conditions, which suggests great stability of RBCS1 under adverse water conditions. On the other hand, for C. arabica, high nocturnal expression of RBCS1 could also explain the accumulation of the RBCS1 protein under water stress. Altogether, the results presented here suggest that the content of RBCS was not responsible for the loss of photosynthetic capacity that is commonly observed in water-stressed coffee plants.

CONCLUSION

We showed that the CaCe homeolog was expressed in C. eugenioides and non-introgressed ("pure") genotypes of C. arabica but that it was undetectable in C. canephora. On the other hand, the CaCc homeolog was expressed in C. canephora but highly repressed in C. arabica. Expression of the CaCc homeolog was enhanced in C. arabica cultivars that experienced recent introgression with C. canephora. For both C. canephora and C. arabica species, total RBCS1 gene expression was highly reduced with WS. Unexpectedly, the accumulation of RBCS1 protein was observed in the leaves of C. canephora under WS, possibly coming from nocturnal RBCS1 expression. These results suggest that the increase in the amount of RBCS1 protein could contribute to the antioxidative function of photorespiration in water-stressed coffee plants.

The impact of cold on photosynthesis in genotypes of Coffea spp.--photosystem sensitivity, photoprotective mechanisms and gene expression

Environmental constraints disturb plant metabolism and are often associated with photosynthetic impairments and yield reductions. Among them, low positive temperatures are of up most importance in tropical plant species, namely in Coffea spp. in which some acclimation ability has been reported. To further explain cold tolerance, the impacts on photosynthetic functioning and the expression of photosynthetic-related genes were analyzed. The experiments were carried out along a period of slow cold imposition (to allow acclimation), after chilling (4°C) exposure and in the following rewarming period, using 1.5-year-old coffee seedlings of 5 genotypes with different cold sensitivity: Coffea canephora cv. Apoatã, Coffea arabica cv. Catuaí, Coffea dewevrei and 2 hybrids, Icatu (C. arabica×C. canephora) and Piatã (C. dewevrei×C. arabica). All genotypes suffered a significant leaf area loss only after chilling exposure, with Icatu showing the lowest impact, a first indication of a higher cold tolerance, contrasting with Apoatã and C. dewevrei. During cold exposure, net photosynthesis and Chl a fluorescence parameters were strongly affected in all genotypes, but stomatal limitations were not detected. However, the extent of mesophyll limitation, reflecting regulatory mechanisms and/or damage, was genotype dependent. Overnight retention of zeaxanthin was common to Coffea genotypes, but the accumulation of photoprotective pigments was highest in Icatu. That down-regulated photochemical events but efficiently protected the photosynthetic structures, as shown, e.g., by the lowest impacts on A(max) and PSI activity and the strongest reinforcement of PSII activity, the latter possibly reflecting the presence of a photoprotective cycle around PSII in Icatu (and Catuaí). Concomitant to these protection mechanisms, Icatu was the sole genotype to present simultaneous upregulation of caCP22, caPI and caCytf, related to, respectively, PSII, PSI and to the complex Cytb(6)/f, which could promote better repair ability, contributing to the maintenance of efficient thylakoid functioning. We conclude that Icatu showed the best acclimation ability among the studied genotypes, mostly due to a better upregulation of photoprotection and repair mechanisms. We confirmed the presence of important variability in Coffea spp. that could be exploited in breeding programs, which should be assisted by useful markers of cold tolerance, namely the upregulation of antioxidative molecules, the expression of selected genes and PSI sensitivity.

Molecular characterization and functional analysis of the β-galactosidase gene during Coffea arabica (L.) fruit development

β-D-Galactosidase (EC 3.2.1.23) has been detected in several plant species, and is characterized in different organs and tissues by its ability to hydrolyse terminal non-reducing β-D-galactosyl residues from β-D-galactoside polymers. In the present paper the cloning and the biochemical and molecular characterization of Coffea arabica β-galactosidase expressed in the pericarp and the endosperm of coffee fruits in all phases of ripeness are described. It was found that coffee β-galactosidase is not evenly transcribed throughout fruit ripening, oscillating with two distinct peaks: the first peak when immature fruits are at the active growing stage and the second when fully developed coffee fruits are completely ripe. Both in vitro enzymatic activity of coffee fruit protein extracts and in vivo histochemical assay of freshly harvested coffee fruits confirmed the uneven transcription of β-galactosidase as fruit maturation advanced. Partial genomic DNA sequencing indicated a complex arrangement of nine putative exons. In silico translation of the cloned coding sequences clearly revealed the cloned gene as β-galactosidase, with the presence of a signal peptide directing the enzyme to the apoplast. Two isoforms were distinguished by sequencing reverse transcription-PCR transcripts, one expressed in young and adult leaves and another in stems, petals, and coffee fruit endosperm and pericarp. Southern blot analysis indicates that there are at least two copies of this gene in the C. arabica genome that could explain the presence of two β-galactosidase isoforms.

Expression patterns of three α-expansin isoforms in Coffea arabica during fruit development

As a first step towards understanding the physiological role and regulation of the expansin gene (EXP) family in Coffea arabica fruits during growth and maturation, we identified 11 expansin genes, nine belonging to the α-expansin family (EXPA), one EXLA and one EXLB, through in silico analysis of expressed sequence tags (ESTs). Within the α-expansin family, three isoforms were selected for detailed examination based on their high expression in coffee fruits or because they were specifically induced during different fruit developmental stages, according to the EST information. The expression patterns were analysed in different fruit tissues (perisperm, endosperm and pericarp) of C. arabica cv. IAPAR-59 and C. arabica cv. IAPAR-59 Graúdo, the latter being a closely related cultivar with a larger fruit size. Accumulation of CaEXPA1 and CaEXPA3 transcripts was high in the perisperm (tissue responsible for coffee bean size) and in the early stages of pericarp development. Transcripts of CaEXPA2 were detected only in the pericarp during the later stages of fruit maturation and ripening. There was no detectable transcription of the three EXPs analysed in the endosperm. The observed differences in mRNA expression levels of CaEXPA1 and CaEXP3 in the perisperm of IAPAR-59 and IAPAR-59 Graúdo suggest the participation of these two isoforms in the regulation of grain size.

Expression of three galactinol synthase isoforms in Coffea arabica L. and accumulation of raffinose and stachyose in response to abiotic stresses

Galactinol synthase (EC 2.4.1.123; GolS) catalyzes the first step in the synthesis of raffinose family oligosaccharides (RFOs). Their accumulation in response to abiotic stresses implies a role for RFOs in stress adaptation. In this study, the expression patterns of three isoforms of galactinol synthase (CaGolS1-2-3) from Coffea arabica were evaluated in response to water deficit, salinity and heat stress. All CaGolS isoforms were highly expressed in leaves while little to no expression were detected in flower buds, flowers, plagiotropic shoots, roots, endosperm and pericarp of mature fruits. Transcriptional analysis indicated that the genes were differentially regulated under water deficit, high salt and heat stress. CaGolS1 isoform is constitutively expressed in plants under normal growth conditions and was the most responsive during all stress treatments. CaGolS2 is unique among the three isoforms in that it was detected only under severe water deficit and salt stresses. CaGolS3 was primarily expressed under moderate and severe drought. This isoform was induced only at the third day of heat and under high salt stress. The increase in GolS transcription was not reflected into the amount of galactinol in coffee leaves, as specific glycosyltransferases most likely used galactinol to transfer galactose units to higher homologous oligosaccharides, as suggested by the increase of raffinose and stachyose during the stresses.

An EST-based analysis identifies new genes and reveals distinctive gene expression features of Coffea arabica and Coffea canephora

BACKGROUND

Coffee is one of the world's most important crops; it is consumed worldwide and plays a significant role in the economy of producing countries. Coffea arabica and C. canephora are responsible for 70 and 30% of commercial production, respectively. C. arabica is an allotetraploid from a recent hybridization of the diploid species, C. canephora and C. eugenioides. C. arabica has lower genetic diversity and results in a higher quality beverage than C. canephora. Research initiatives have been launched to produce genomic and transcriptomic data about Coffea spp. as a strategy to improve breeding efficiency.

RESULTS

Assembling the expressed sequence tags (ESTs) of C. arabica and C. canephora produced by the Brazilian Coffee Genome Project and the Nestlé-Cornell Consortium revealed 32,007 clusters of C. arabica and 16,665 clusters of C. canephora. We detected different GC3 profiles between these species that are related to their genome structure and mating system. BLAST analysis revealed similarities between coffee and grape (Vitis vinifera) genes. Using KA/KS analysis, we identified coffee genes under purifying and positive selection. Protein domain and gene ontology analyses suggested differences between Coffea spp. data, mainly in relation to complex sugar synthases and nucleotide binding proteins. OrthoMCL was used to identify specific and prevalent coffee protein families when compared to five other plant species. Among the interesting families annotated are new cystatins, glycine-rich proteins and RALF-like peptides. Hierarchical clustering was used to independently group C. arabica and C. canephora expression clusters according to expression data extracted from EST libraries, resulting in the identification of differentially expressed genes. Based on these results, we emphasize gene annotation and discuss plant defenses, abiotic stress and cup quality-related functional categories.

CONCLUSION

We present the first comprehensive genome-wide transcript profile study of C. arabica and C. canephora, which can be freely assessed by the scientific community at http://www.lge.ibi.unicamp.br/coffea. Our data reveal the presence of species-specific/prevalent genes in coffee that may help to explain particular characteristics of these two crops. The identification of differentially expressed transcripts offers a starting point for the correlation between gene expression profiles and Coffea spp. developmental traits, providing valuable insights for coffee breeding and biotechnology, especially concerning sugar metabolism and stress tolerance.

A high-throughput data mining of single nucleotide polymorphisms in Coffea species expressed sequence tags suggests differential homeologous gene expression in the allotetraploid Coffea arabica

Polyploidization constitutes a common mode of evolution in flowering plants. This event provides the raw material for the divergence of function in homeologous genes, leading to phenotypic novelty that can contribute to the success of polyploids in nature or their selection for use in agriculture. Mounting evidence underlined the existence of homeologous expression biases in polyploid genomes; however, strategies to analyze such transcriptome regulation remained scarce. Important factors regarding homeologous expression biases remain to be explored, such as whether this phenomenon influences specific genes, how paralogs are affected by genome doubling, and what is the importance of the variability of homeologous expression bias to genotype differences. This study reports the expressed sequence tag assembly of the allopolyploid Coffea arabica and one of its direct ancestors, Coffea canephora. The assembly was used for the discovery of single nucleotide polymorphisms through the identification of high-quality discrepancies in overlapped expressed sequence tags and for gene expression information indirectly estimated by the transcript redundancy. Sequence diversity profiles were evaluated within C. arabica (Ca) and C. canephora (Cc) and used to deduce the transcript contribution of the Coffea eugenioides (Ce) ancestor. The assignment of the C. arabica haplotypes to the C. canephora (CaCc) or C. eugenioides (CaCe) ancestral genomes allowed us to analyze gene expression contributions of each subgenome in C. arabica. In silico data were validated by the quantitative polymerase chain reaction and allele-specific combination TaqMAMA-based method. The presence of differential expression of C. arabica homeologous genes and its implications in coffee gene expression, ontology, and physiology are discussed.

Use of the growing environment as a source of variation to identify the quantitative trait transcripts and modules of co-expressed genes that determine chlorogenic acid accumulation

Developing Coffea arabica seeds accumulate large amounts of chlorogenic acids (CGAs) as a storage form of phenylpropanoid derivatives, making coffee a valuable model to investigate the metabolism of these widespread plant phenolics. However, developmental and environmental regulations of CGA metabolism are poorly understood. In the present work, the expression of selected phenylpropanoid genes, together with CGA isomer profiles, was monitored throughout seed development across a wide set of contrasted natural environments. Although CGA metabolism was controlled by major developmental factors, the mean temperature during seed development had a direct impact on the time-window of CGA biosynthesis, as well as on final CGA isomer composition through subtle transcriptional regulations. We provide evidence that the variability induced by the environment is a useful tool to test whether CGA accumulation is quantitatively modulated at the transcriptional level, hence enabling detection of rate-limiting transcriptional steps [quantitative trait transcripts (QTTs)] for CGA biosynthesis. Variations induced by the environment also enabled a better description of the phenylpropanoid gene transcriptional network throughout seed development, as well as the detection of three temporally distinct modules of quantitatively co-expressed genes. Finally, analysis of metabolite-to-metabolite relationships revealed new biochemical characteristics of the isomerization steps that remain uncharacterized at the gene level.

Comparative sequence analyses indicate that Coffea (Asterids) and Vitis (Rosids) derive from the same paleo-hexaploid ancestral genome

The complete sequence of Vitis vinifera revealed that the rosid clade derives from a hexaploid ancestor. At present, no analysis of complete genome sequence is available for an asterid, the other large eudicot clade, which includes the economically important species potato, tomato and coffee. To elucidate the genomic history of asterids, we compared the sequence of an 800 kb region of diploid Coffea genome to the orthologous regions of V. vinifera, Populus trichocarpa and Arabidopsis thaliana. We found a very high level of collinearity between around 80 genes of the three rosid species and Coffea. Collinearity comparisons between orthologous and paralogous regions indicates that (1) the Coffea (and consequently all asterids) and rosids share the same hexaploid ancestor; (2) the diploidization process (loss of duplicated and redundant copies from the whole genome duplication) was very advanced in the most recent common ancestor of rosids and asterids. Finally, no additional polyploidization events were detected in the Coffea lineage. Differences in gene loss rates were detected among the three rosid species and linked to the divergence in protein sequences.

Somatic embryogenesis and vegetative cutting capacity are under distinct genetic control in Coffea canephora Pierre

The purpose of the study was to evaluate the possible genetic effect on vegetative propagation of Coffea canephora. Diversity for somatic embryogenesis (SE) ability was observed not only among two groups of C. canephora Pierre (Congolese and Guinean), but also within these different genetic groups. The results therefore showed that, under given experimental conditions, SE ability is depending on genotype. Furthermore the detection of quantitative trait loci (QTLs) controlling the SE and cutting abilities of C. canephora was performed on a large number of clones including accessions from a core collection, three parental clones and their segregating progenies. On the one hand we detected eight QTLs determining SE. Six positive QTLs for SE ability, whatever the criteria used to quantify this ability, were localized on one single chromosome region of the consensus genetic map. Two negative QTLs for SE ability (frequency of micro calli without somatic embryo) were detected on another linkage group. Deep analysis of the six QTLs detected for SE ability came to the conclusion that they can be assimilated to one single QTL explaining 8.6-12.2% of the observed variation. On the other hand, two QTLs for average length of roots and length of the longest sprouts of cuttings were detected in two linkage groups. These QTLs detected for cutting ability are explaining 12-27% of the observed variation. These observations led to conclude that SE and cutting abilities of C. canephora Pierre appeared to be genetic dependent but through independent mechanisms.

Stress metabolism in green coffee beans (Coffea arabica L.): expression of dehydrins and accumulation of GABA during drying

In order to produce tradeable standard green coffee, processed beans must be dried. The drying procedure affects the abundance of relevant aroma substances, e.g. carbohydrates. Using molecular tools, the corresponding metabolic basis is analyzed. A decrease in water potential of the still living coffee seeds induces massive drought stress responses. As a marker for these stress reactions, accumulation of a general stress metabolite, GABA (gamma-aminobutyric acid), and associated gene expression of drought stress-associated dehydrins were monitored. The results of this study indicate that metabolism in drying coffee beans is quite complex since several events trigger accumulation of GABA. The first peak of GABA accumulation during drying is correlated with expression of isocitrate lyase and thus with ongoing germination processes in coffee seeds. Two subsequent peaks of GABA accumulation correspond to maxima of dehydrin gene expression and are thought to be induced directly by drought stress in the embryo and endosperm tissue, respectively. Apart from the significance for understanding basic seed physiology, metabolic changes in coffee seeds during processing provide valuable information for understanding the role and effect of the steps of green coffee processing on the quality of the resulting coffee.

Transcriptional profile of genes involved in the biosynthesis of phytate and ferritin in Coffea

The present work aimed to study the control of the biosynthesis of the antinutritional factor phytate and its associated Fe-rich protein family, ferritin, in coffee. Phytate has the ability to chelate Fe, making it unavailable to human absorption. The Coffea genome databases were queried for genes associated with phytate metabolism and ferritin genes. The genetic framework for phytate biosynthesis and its reverse pathway was identified in silico analyses and indicate that Coffea phosphatidyl inositol kinase and monophosphatase families play nonredundant roles in phytate metabolism. The transcriptional profiles of phytate biosynthesis key-genes MYO-INOSITOL(3)P1 SYNTHASE, two genes coding for PHOSPHATIDYL INOSITOL KINASE, and three FERRITIN genes were temporally evaluated by qPCR in coffee seeds from two crop locations, Adamantina-SP and Ouro-Fino-MG, the last one traditionally associated with high-quality coffee beverage grain. A targeted metabolome profile of phytic acid contents throughout three fruit maturation stages in association with the transcriptional analysis was also obtained. Taken together, our data indicate that the investigated local conditions did not cause significant alterations in phytate biosynthesis. Futhermore, the temporal transcriptional profiling revealed that candidate gene expression is regulated independently of phytate accumulation. In contrast, the expression profile of ferritin-unit genes is affected by environmental conditions and genetic background. The roles of the investigated genes are discussed concerning the quality of coffee beverage.

Biochemical and molecular characterization of the antioxidative system of Coffea sp. under cold conditions in genotypes with contrasting tolerance

Low positive temperature (chilling) is frequently linked to the promotion of oxidative stress conditions, and is of particular importance in the coffee plant due to its severe impact on growth, development, photosynthesis and production. Nevertheless, some acclimation ability has been reported within the Coffea genus, and is possibly related to oxidative stress control. Using an integrated biochemical and molecular approach, the characterization of the antioxidative system of genotypes with different cold acclimation abilities was performed. Experiments were carried out using 1.5-year-old coffee seedlings of Coffea canephora cv. Apoatã, C. arabica cv. Catuaí, C. dewevrei and 2 hybrids, Icatu (C. arabicaxC. canephora) and Piatã (C. dewevreixC. arabica) subjected to a gradual cold treatment and a recovery period. Icatu showed the greatest ability to control oxidative stress, as reflected by the enhancement of several antioxidative components (Cu,Zn-SOD and APX activities; ascorbate, alpha-tocopherol and chlorogenic acids (CGAs) contents) and lower reactive oxygen species contents (H(2)O(2) and OH). Gene expression studies show that GRed, DHAR and class III and IV chitinases might also be involved in the cold acclimation ability of Icatu. Catuaí showed intermediate acclimation ability through the reinforcement of some antioxidative molecules, usually to a lesser extent than that observed in Icatu. On the other hand, C. dewevrei showed the poorest response in terms of antioxidant accumulation, and also showed the greatest increase in OH values. The difference in the triggering of antioxidative traits supports the hypothesis of its importance to cold (and photoinhibition) tolerance in Coffea sp. and could provide a useful probe to identify tolerant genotypes.

Producing a recombinant flavin-containing monooxygenase from Coffea arabica in Escherichia coli for screening of potential natural substrates

Only few biological functions have been related with flavin-containing monooxygenases (FMOs) in plants, such as specific roles in auxin biosynthesis, pathogen defense, and metabolism of glucosinolates. Biochemical characterization using recombinant proteins is a promising approach to determine the precise specificity of plant FMOs for potential natural substrates. FMOs may be very difficult to express in a soluble form due to their highly hydrophobic nature and this can be improved by fusing them to solubility-enhancing proteins, such as maltose-binding protein (MBP) and N-utilization substance A (NusA). Here we describe the expression of a recombinant FMO from Coffea arabica as a maltose-binding protein fusion in Escherichia coli and its purification by affinity chromatography, producing a ready-to-use protein for enzymatic activity assays.

Genotypic and environmental effects on coffee (Coffea arabica L.) bean fatty acid profile: impact on variety and origin chemometric determination

In a previous study, the effectiveness of chlorogenic acids, fatty acids (FA), and elements was compared for the discrimination of Arabica varieties and growing terroirs. Since FA provided the best results, the aim of the present work was to validate their discrimination ability using an extended experimental design, including twice the number of location x variety combinations and 2 years of study. It also aimed at understanding how the environment influences FA composition through correlation analysis using different climatic parameters. Percentages of correct classification of known samples remained very high, independent of the classification criterion. However, cross-validation tests across years indicated that prediction of unknown locations was less efficient than that of unknown genotypes. Environmental temperature during the development of coffee beans had a dramatic influence on their FA composition. Analysis of climate patterns over years enabled us to understand the efficient location discrimination within a single year but only moderate efficiency across years.

Proteomic evaluation of coffee zygotic embryos in two different stages of seed development

Coffee seed development is accompanied by severe modifications in water soluble proteins, several of these being associated to a specific developmental stage. For this reason, a proteomic approach has been used to describe spatial-temporal proteome modifications in zygotic embryos at different stages of seed development. Embryos from Coffea arabica seeds were harvested in two different developmental stages: stage 1 at 210 days after anthesis and stage 2 at 255 days. Total proteins were extracted and submitted to 2-DE. From these gels, several spots were identified by mass spectrometry including kinases, MYB transcription factor and enzymes involved in metabolic pathways. All proteins identified seem to affect coffee development in different ways, being directly involved in plant growth or used as an intermediate in some metabolic pathway that, indirectly, will influence coffee development. This is the first work using two-dimensional electrophoresis followed by mass spectrometry analyses that evaluates the expression of proteins during coffee zygotic embryos development. Data here reported supply some light over coffee development and could be used in a near future to improve coffee plants' growth and development by molecular strategies.

MoccaDB - an integrative database for functional, comparative and diversity studies in the Rubiaceae family

BACKGROUND

In the past few years, functional genomics information has been rapidly accumulating on Rubiaceae species and especially on those belonging to the Coffea genus (coffee trees). An increasing number of expressed sequence tag (EST) data and EST- or genomic-derived microsatellite markers have been generated, together with Conserved Ortholog Set (COS) markers. This considerably facilitates comparative genomics or map-based genetic studies through the common use of orthologous loci across different species. Similar genomic information is available for e.g. tomato or potato, members of the Solanaceae family. Since both Rubiaceae and Solanaceae belong to the Euasterids I (lamiids) integration of information on genetic markers would be possible and lead to more efficient analyses and discovery of key loci involved in important traits such as fruit development, quality, and maturation, or adaptation. Our goal was to develop a comprehensive web data source for integrated information on validated orthologous markers in Rubiaceae.

DESCRIPTION

MoccaDB is an online MySQL-PHP driven relational database that houses annotated and/or mapped microsatellite markers in Rubiaceae. In its current release, the database stores 638 markers that have been defined on 259 ESTs and 379 genomic sequences. Marker information was retrieved from 11 published works, and completed with original data on 132 microsatellite markers validated in our laboratory. DNA sequences were derived from three Coffea species/hybrids. Microsatellite markers were checked for similarity, in vitro tested for cross-amplification and diversity/polymorphism status in up to 38 Rubiaceae species belonging to the Cinchonoideae and Rubioideae subfamilies. Functional annotation was provided and some markers associated with described metabolic pathways were also integrated. Users can search the database for marker, sequence, map or diversity information through multi-option query forms. The retrieved data can be browsed and downloaded, along with protocols used, using a standard web browser. MoccaDB also integrates bioinformatics tools (CMap viewer and local BLAST) and hyperlinks to related external data sources (NCBI GenBank and PubMed, SOL Genomic Network database).

CONCLUSION

We believe that MoccaDB will be extremely useful for all researchers working in the areas of comparative and functional genomics and molecular evolution, in general, and population analysis and association mapping of Rubiaceae and Solanaceae species, in particular.

Physical mapping of rDNA and heterochromatin in chromosomes of 16 Coffea species: a revised view of species differentiation

The chromosome organization among 15 wild diploid Coffea species and cultivated tetraploid C. arabica was determined by fluorochrome banding (CMA, DAPI) and double fluorescence in-situ hybridization (FISH) of 5S and 18S rDNA achieved on the same chromosome plates. Two to five chromosome pairs (plus one putative chromosome B) are marked. Overall, there are two SAT-chromosome pairs for East African species and one for the Malagasy and the West and Central African species. 18S rDNA loci are telomeric and strongly marked the SAT-chromosome pairs. Generally, only one pericentromeric 5S rDNA locus characterized East African species, while an additional minor locus co-localized with the 18S rDNA-SAT locus for the Malagasy species and West and Central African species. A combination of rDNA FISH plus CMA and DAPI banding patterns enables identification of almost all the species, even those for which the genetic or botanical status is still being discussed. C. arabica clearly appears to be an allotetraploid species, including one genome from East Africa and one from West and Central Africa. However, since the minor 5S rDNA-SAT locus present in West/Central African genomes is not detected, two evolutionary hypotheses could be put forward for C. arabica. Considering only the diploid species, global trends are obvious in rDNA signal patterns, genome size variations, and geographic distribution of the species, but there are no clear evolutionary trends. However, complex interactions between these factors and environmental growing conditions exist, which have resulted in loss and gain of rDNA loci and probably also in copy repeat number variations in each rDNA family.

Microcollinearity in an ethylene receptor coding gene region of the Coffea canephora genome is extensively conserved with Vitis vinifera and other distant dicotyledonous sequenced genomes

BACKGROUND

Coffea canephora, also called Robusta, belongs to the Rubiaceae, the fourth largest angiosperm family. This diploid species (2x = 2n = 22) has a fairly small genome size of approximately 690 Mb and despite its extreme economic importance, particularly for developing countries, knowledge on the genome composition, structure and evolution remain very limited. Here, we report the 160 kb of the first C. canephora Bacterial Artificial Chromosome (BAC) clone ever sequenced and its fine analysis.

RESULTS

This clone contains the CcEIN4 gene, encoding an ethylene receptor, and twenty other predicted genes showing a high gene density of one gene per 7.8 kb. Most of them display perfect matches with C. canephora expressed sequence tags or show transcriptional activities through PCR amplifications on cDNA libraries. Twenty-three transposable elements, mainly Class II transposon derivatives, were identified at this locus. Most of these Class II elements are Miniature Inverted-repeat Transposable Elements (MITE) known to be closely associated with plant genes. This BAC composition gives a pattern similar to those found in gene rich regions of Solanum lycopersicum and Medicago truncatula genomes indicating that the CcEIN4 regions may belong to a gene rich region in the C. canephora genome. Comparative sequence analysis indicated an extensive conservation between C. canephora and most of the reference dicotyledonous genomes studied in this work, such as tomato (S. lycopersicum), grapevine (V. vinifera), barrel medic M. truncatula, black cottonwood (Populus trichocarpa) and Arabidopsis thaliana. The higher degree of microcollinearity was found between C. canephora and V. vinifera, which belong respectively to the Asterids and Rosids, two clades that diverged more than 114 million years ago.

CONCLUSION

This study provides a first glimpse of C. canephora genome composition and evolution. Our data revealed a remarkable conservation of the microcollinearity between C. canephora and V. vinifera and a high conservation with other distant dicotyledonous reference genomes. Altogether, these results provide valuable information to identify candidate genes in C. canephora genome and serve as a foundation to establish strategies for whole genome sequencing. Future large-scale sequence comparison between C. canephora and reference sequenced genomes will help in understanding the evolutionary history of dicotyledonous plants.

Identification of suitable internal control genes for expression studies in Coffea arabica under different experimental conditions

BACKGROUND

Quantitative data from gene expression experiments are often normalized by transcription levels of reference or housekeeping genes. An inherent assumption for their use is that the expression of these genes is highly uniform in living organisms during various phases of development, in different cell types and under diverse environmental conditions. To date, the validation of reference genes in plants has received very little attention and suitable reference genes have not been defined for a great number of crop species including Coffea arabica. The aim of the research reported herein was to compare the relative expression of a set of potential reference genes across different types of tissue/organ samples of coffee. We also validated the expression profiles of the selected reference genes at various stages of development and under a specific biotic stress.

RESULTS

The expression levels of five frequently used housekeeping genes (reference genes), namely alcohol dehydrogenase (adh), 14-3-3, polyubiquitin (poly), beta-actin (actin) and glyceraldehyde-3-phosphate dehydrogenase (gapdh) was assessed by quantitative real-time RT-PCR over a set of five tissue/organ samples (root, stem, leaf, flower, and fruits) of Coffea arabica plants. In addition to these commonly used internal controls, three other genes encoding a cysteine proteinase (cys), a caffeine synthase (ccs) and the 60S ribosomal protein L7 (rpl7) were also tested. Their stability and suitability as reference genes were validated by geNorm, NormFinder and BestKeeper programs. The obtained results revealed significantly variable expression levels of all reference genes analyzed, with the exception of gapdh, which showed no significant changes in expression among the investigated experimental conditions.

CONCLUSION

Our data suggests that the expression of housekeeping genes is not completely stable in coffee. Based on our results, gapdh, followed by 14-3-3 and rpl7 were found to be homogeneously expressed and are therefore adequate for normalization purposes, showing equivalent transcript levels in different tissue/organ samples. Gapdh is therefore the recommended reference gene for measuring gene expression in Coffea arabica. Its use will enable more accurate and reliable normalization of tissue/organ-specific gene expression studies in this important cherry crop plant.

Comparison of the Coffea canephora and C. arabica karyotype based on chromosomal DNA content

Nuclear genome size has been measured in various plants, seeing that knowledge of the DNA content is useful for taxonomic and evolutive studies, plant breeding programs and genome sequencing projects. Besides the nuclear DNA content, tools and protocols to quantify the chromosomal DNA content have been also applied, expanding the data about genomic structure. This study was conducted in order to calculate the Coffea canephora and Coffea arabica chromosomal DNA content, associating cytogenetic methodologies with flow cytometry (FCM) and image cytometry (ICM) tools. FCM analysis showed that the mean nuclear DNA content of C. canephora and C. arabica is 2C = 1.41 and 2.62 pg, respectively. The cytogenetic methodology provided prometaphase and metaphase cells exhibiting adequate chromosomes for the ICM measurements and karyogram assembly. Based on cytogenetic, FCM and ICM results; it was possible to calculate the chromosomal DNA content of the two species. The 1C chromosomal DNA content of C. canephora ranged from 0.09 (chromosome 1) to 0.05 pg (chromosome 11) and C. arabica from 0.09 (chromosome 1) to 0.03 pg (chromosome 22). The methodology presented in this study was suitable for DNA content measuring of each chromosome of C. canephora and C. arabica. The cytogenetic characterization and chromosomal DNA content analyses evidenced that C. arabica is a true allotetraploid originated from a cross between Coffea diploid species. Besides, the same analyses also reinforce that C. canephora is a possible progenitor of C. arabica.

Metabolic pathways in tropical dicotyledonous albuminous seeds: Coffea arabica as a case study

* The genomic era facilitates the understanding of how transcriptional networks are interconnected to program seed development and filling. However, to date, little information is available regarding dicot seeds with a transient perisperm and a persistent, copious endosperm. Coffea arabica is the subject of increasing genomic research and is a model for nonorthodox albuminous dicot seeds of tropical origin. * The aim of this study was to reconstruct the metabolic pathways involved in the biosynthesis of the main coffee seed storage compounds, namely cell wall polysaccharides, triacylglycerols, sucrose, and chlorogenic acids. For this purpose, we integrated transcriptomic and metabolite analyses, combining real-time RT-PCR performed on 137 selected genes (of which 79 were uncharacterized in Coffea) and metabolite profiling. * Our map-drawing approach derived from model plants enabled us to propose a rationale for the peculiar traits of the coffee endosperm, such as its unusual fatty acid composition, remarkable accumulation of chlorogenic acid and cell wall polysaccharides. * Comparison with the developmental features of exalbuminous seeds described in the literature revealed that the two seed types share important regulatory mechanisms for reserve biosynthesis, independent of the origin and ploidy level of the storage tissue.

Effectiveness of the fatty acid and sterol composition of seeds for the chemotaxonomy of Coffea subgenus Coffea

The chemotaxonomic relationships between Coffea (subgenus Coffea) species have been poorly studied to date and the compounds tested so far - chlorogenic acids, diterpenoids and purine alkaloids - did not enable the establishment of phylogenetic relationships analogous to those revealed by chloroplast and nuclear DNA studies. In the present study, the relationships between African Coffea species were assessed on the basis of their seed lipid composition. Fatty acids and sterols were determined in 59 genotypes belonging to 17 distinct Coffea species/origins. Principal Component Analysis of fatty acid and sterol data enabled easy identification of the few species for which one or several compounds could serve as a quantitative signature. Hierarchical Clustering classified the Coffea species in seven groups with both fatty acids and sterols. However, while groupings based on seed fatty acid composition showed remarkable ecological and geographical coherence, no phylogeographic explanation was found for the clusters retrieved from sterol data. When compared with previous phylogenetic studies, the groups deduced from seed fatty acid composition were remarkably congruent with the clades inferred from nuclear and plastid DNA sequences.

Effect of coffee alkaloids and phenolics on egg-laying by the coffee leaf miner Leucoptera coffeella

The recognized importance of coffee alkaloids and phenolics mediating insect-plant interactions led to the present investigation aiming to test the hypothesis that the phenolics chlorogenic and caffeic acids and the alkaloid caffeine and some of its derivatives present in coffee leaves affect egg-laying by the coffee leaf miner Leucoptera (=Perileucoptera) coffeella (Guérin-Méneville & Perrottet) (Lepidoptera: Lyonetiidae), one of the main coffee pests in the Neotropical region. These phytochemicals were, therefore, quantified in leaves from 12 coffee genotypes and their effect on the egg-laying preference by the coffee leaf miner was assessed. Canonical variate analysis and partial canonical correlation provided evidence that increased leaf levels of caffeine favour egg-laying by the coffee leaf miner. An egg-laying preference bioassay was, therefore, carried out to specifically test this hypothesis using increasing caffeine concentrations sprayed on leaves of one of the coffee genotypes with the lowest level of this compound (i.e. Hybrid UFV 557-04 generated from a cross between Coffea racemosa Lour. and C. arabica L.). The results obtained allowed the recognition of a significant concentration-response relationship, providing support for the hypothesis that caffeine stimulates egg-laying by the coffee leaf miner in coffee leaves.

An investigation of carotenoid biosynthesis in Coffea canephora and Coffea arabica

Carotenoids are essential components of the photosynthetic apparatus in a wide range of organisms. They participate in the adaptation of plastids to changing environmental light conditions and prevent photo-oxidative damage of the photosynthetic apparatus by detoxifying reactive oxygen species. We identified eight cDNAs from the carotenoid biosynthetic pathway (PSY, PDS, ZDS, PTOX, LCY-E, CRTR-B, ZEP and VDE) and two cDNA encoding carotenoid cleavage dioxygenase family members (NCED3 and CCD1) in Coffea canephora. We also obtained cDNA encoding several different fibrillin proteins involved in carotenoid sequestration (FIB). Expression of the coffee carotenoid genes was determined in leaf, branch and flower tissues using quantitative RT-PCR. Expression analysis of these genes in leaf tissue from osmotically stressed plants was also carried out. These experiments showed that the transcript levels of PTOX, CRTR-B, NCED3, CCD1 and FIB1 increased under these stress conditions, while LCY-E decreased, indicating that the metabolic flux towards the xanthophyll cycle branch of the carotenoid biosynthetic pathway may be favoured in leaves under drought conditions. Functional analysis of CcCRTR-B using an in vivo method employing Escherichia coli strains engineered to make carotenoids confirmed that the beta-carotene hydroxylase activity of CcCRTR-B generates beta-cryptoxanthin and zeaxanthin from beta-carotene. A similar approach was also used to show that CcCCD1 encoded a functional 9,10(9'10') carotenoid cleavage dioxygenase, and thus that this enzyme is capable of forming one or more apocarotenoids in vivo. Finally, high-performance liquid chromatography analysis of coffee leaves revealed the presence of alpha-carotene and suggests that Coffea arabica may have higher levels of alpha-carotene than C. canephora.

Agrobacterium rhizogenes-transformed roots of coffee (Coffea arabica): conditions for long-term proliferation, and morphological and molecular characterization

BACKGROUND AND AIMS

The aims of this study were to set up proliferation conditions for hairy roots of Coffea arabica regenerated after transformation by Agrobacterium rhizogenes strain A4-RS, and to carry out the morphological and molecular characterization of hairy root clones maintained over the long term.

METHODS

Auxin supply, light conditions and sucrose concentration were modified with the aim of establishing efficient root proliferation conditions. The morphological variability among 62 established hairy root clones was phenotyped by scanning the roots and analysing the images using 'whinRHIZO' software procedures. PCR analysis of integration in transformed root cells of rol and aux oncogenes from the T-DNA of the Ri plasmid was used to study the molecular variability among clones.

KEY RESULTS

Auxin supply was necessary to obtain and stimulate growth and branching, and IBA applied at 0.5 microm was the most efficient auxin. Significant differences were shown among the 62 clones for total root length and for the percentage of fine roots. These variables were stable across subcultures and could hence be used for efficient characterization of hairy root clones. The majority of hairy root clones (86 %) exhibited non-significant phenotype differences with non-transformed roots. Eight clones were significantly different from the non-transformed controls in that they possessed a low proportion of fine roots. Two other hairy root clones grew significantly faster than the other clones. The PCR analysis revealed a low variability in the integration of rol and aux oncogenes in transformed root cells. The T(R)-DNA was never integrated as aux1 and aux2 genes were not found, although rolB and rolC genes from the T(L)-DNA were always present.

CONCLUSIONS

The discovery of low morphological variability among coffee hairy roots together with the identification of morphological variables allowing easy identification of phenotypically altered clones represent two important results. They make hairy roots a possible, and efficient, tool for functional-genomic studies of coffee root genes.

Effects of shade on the development and sugar metabolism of coffee (Coffea arabica L.) fruits

Coffee fruits grown in shade are characterized by larger bean size than those grown under full-sun conditions. The present study assessed the effects of shade on bean characteristics and sugar metabolism by analyzing tissue development, sugar contents, activities of sucrose metabolizing enzymes and expression of sucrose synthase-encoding genes in fruits of coffee (Coffea arabica L.) plants submitted to full-sun (FS) and shade (SH) conditions. Evolution of tissue fresh weights measured in fruits collected regularly from flowering to maturation indicated that this increase is due to greater development of the perisperm tissue in the shade. The effects of light regime on sucrose and reducing sugar (glucose and fructose) contents were studied in fresh and dry coffee beans. Shade led to a significant reduction in sucrose content and to an increase in reducing sugars. In pericarp and perisperm tissues, higher activities of sucrose synthase (EC 2.4.1.13) and sucrose-phosphate synthase (SPS: EC 2.4.1.14) were detected at maturation in the shade compared with full sun. These two enzymes also had higher peaks of activities in developing endosperm under shade than in full sun. It was also noted that shade modified the expression of SUS-encoding genes in coffee beans; CaSUS2 gene transcripts levels were higher in SH than in FS. As no sucrose increase accompanied these changes, this suggests that sucrose metabolism was redirected to other metabolic pathways that need to be identified.

Development of new genomic microsatellite markers from robusta coffee (Coffea canephora Pierre ex A. Froehner) showing broad cross-species transferability and utility in genetic studies

BACKGROUND

Species-specific microsatellite markers are desirable for genetic studies and to harness the potential of MAS-based breeding for genetic improvement. Limited availability of such markers for coffee, one of the most important beverage tree crops, warrants newer efforts to develop additional microsatellite markers that can be effectively deployed in genetic analysis and coffee improvement programs. The present study aimed to develop new coffee-specific SSR markers and validate their utility in analysis of genetic diversity, individualization, linkage mapping, and transferability for use in other related taxa.

RESULTS

A small-insert partial genomic library of Coffea canephora, was probed for various SSR motifs following conventional approach of Southern hybridisation. Characterization of repeat positive clones revealed a very high abundance of DNRs (1/15 Kb) over TNRs (1/406 kb). The relative frequencies of different DNRs were found as AT >> AG > AC, whereas among TNRs, AGC was the most abundant repeat. The SSR positive sequences were used to design 58 primer pairs of which 44 pairs could be validated as single locus markers using a panel of arabica and robusta genotypes. The analysis revealed an average of 3.3 and 3.78 alleles and 0.49 and 0.62 PIC per marker for the tested arabicas and robustas, respectively. It also revealed a high cumulative PI over all the markers using both sib-based (10-6 and 10-12 for arabicas and robustas respectively) and unbiased corrected estimates (10-20 and 10-43 for arabicas and robustas respectively). The markers were tested for Hardy-Weinberg equilibrium, linkage dis-equilibrium, and were successfully used to ascertain generic diversity/affinities in the tested germplasm (cultivated as well as species). Nine markers could be mapped on robusta linkage map. Importantly, the markers showed ~92% transferability across related species/genera of coffee.

CONCLUSION

The conventional approach of genomic library was successfully employed although with low efficiency to develop a set of 44 new genomic microsatellite markers of coffee. The characterization/validation of new markers demonstrated them to be highly informative, and useful for genetic studies namely, genetic diversity in coffee germplasm, individualization/bar-coding for germplasm protection, linkage mapping, taxonomic studies, and use as conserved orthologous sets across secondary genepool of coffee. Further, the relative frequency and distribution of different SSR motifs in coffee genome indicated coffee genome to be relatively poor in microsatellites compared to other plant species.

Differential regulation of grain sucrose accumulation and metabolism in Coffea arabica (Arabica) and Coffea canephora (Robusta) revealed through gene expression and enzyme activity analysis

* Coffea arabica (Arabica) and Coffea canephora (Robusta) are the two main cultivated species used for coffee bean production. Arabica genotypes generally produce a higher coffee quality than Robusta genotypes. Understanding the genetic basis for sucrose accumulation during coffee grain maturation is an important goal because sucrose is an important coffee flavor precursor. * Nine new Coffea genes encoding sucrose metabolism enzymes have been identified: sucrose phosphate synthase (CcSPS1, CcSPS2), sucrose phosphate phosphatase (CcSP1), cytoplasmic (CaInv3) and cell wall (CcInv4) invertases and four invertase inhibitors (CcInvI1, 2, 3, 4). * Activities and mRNA abundance of the sucrose metabolism enzymes were compared at different developmental stages in Arabica and Robusta grains, characterized by different sucrose contents in mature grain. * It is concluded that Robusta accumulates less sucrose than Arabica for two reasons: Robusta has higher sucrose synthase and acid invertase activities early in grain development - the expression of CcSS1 and CcInv2 appears to be crucial at this stage and Robusta has a lower SPS activity and low CcSPS1 expression at the final stages of grain development and hence has less capacity for sucrose re-synthesis. Regulation of vacuolar invertase CcInv2 activity by invertase inhibitors CcInvI2 and/or CcInvI3 during Arabica grain development is considered.

Caffeine and related purine alkaloids: biosynthesis, catabolism, function and genetic engineering

Details of the recently elucidated biosynthetic pathways of caffeine and related purine alkaloids are reviewed. The main caffeine biosynthetic pathway is a sequence consisting of xanthosine-->7-methylxanthosine-->7-methylxanthine-->theobromine-->caffeine. Genes encoding N-methyltransferases involved in three of these four reactions have been isolated and the molecular structure of N-methyltransferases investigated. Pathways for the catabolism of caffeine have also been studied, although there are currently no reports of enzymatic and genetic studies having been successfully carried out. Metabolism of purine alkaloids in species including Camellia, Coffea, Theobroma and Ilex plants is summarised, and evidence for the involvement of caffeine in chemical defense and allelopathy is discussed. Finally, information is presented on metabolic engineering that has produced coffee seedlings with reduced caffeine content, and transgenic caffeine-producing tobacco plants with enhanced disease resistance.

Deciphering transcriptional networks that govern Coffea arabica seed development using combined cDNA array and real-time RT-PCR approaches

Due to its economic importance, Coffea arabica is becoming the subject of increasing genomic research and, in particular, the genes involved in the final chemical composition of the bean and the sensorial quality of the coffee beverage. The aim of the present study was to decipher the transcriptional networks that govern the development of the C. arabica seed, a model for non-orthodox albuminous seeds of tropical origin. For this purpose, we developed a transcriptomic approach combining two techniques: targeted cDNA arrays, containing 266 selected candidate gene sequences, and real-time RT-PCR on a large subset of 111 genes. The combination of the two techniques allowed us to limit detection of false positives and to reveal the advantages of using large real-time RT-PCR screening. Multivariate analysis was conducted on both datasets and results were broadly convergent. First, principle component analysis (PCA) revealed a dramatic re-programming of the transcriptional machinery between early cell division and elongation, storage and maturation phases. Second, hierarchical clustering analysis (HCA) led to the identification of 11 distinct patterns of gene expression during seed development as well as to the detection of genes expressed at specific developmental stages that can be used as functional markers of phenological changes. In addition, this study led to the description of gene expression profiles for quality-related genes, most of them formerly uncharacterised in Coffea. Their involvement in storage compound synthesis and accumulation during endosperm development and final metabolic re-adjustments during maturation is discussed.

Towards a Phylogeny for Coffea (Rubiaceae): identifying well-supported lineages based on nuclear and plastid DNA sequences

BACKGROUND AND AIMS

The phylogenetic relationships between species of Coffea and Psilanthus remain poorly understood, owing to low levels of sequence variation recovered in previous studies, coupled with relatively limited species sampling. In this study, the relationships between Coffea and Psilanthus species are assessed based on substantially increased molecular sequence data and greatly improved species sampling.

METHODS

Phylogenetic relationships are assessed using parsimony, with sequence data from four plastid regions [trnL-F intron, trnL-F intergenic spacer (IGS), rpl16 intron and accD-psa1 IGS], and the internal transcribed spacer (ITS) region of nuclear rDNA (ITS 1/5.8S/ITS 2). Supported lineages in Coffea are discussed within the context of geographical correspondence, biogeography, morphology and systematics.

KEY RESULTS

Several major lineages with geographical coherence, as identified in previous studies based on smaller data sets, are supported. Other lineages with either geographical or ecological correspondence are recognized for the first time. Coffea subgenus Baracoffea is shown to be monophyletic, but Coffea subgenus Coffea is paraphyletic. Sequence data do not substantiate the monophyly of either Coffea or Psilanthus. Low levels of sequence divergence do not allow detailed resolution of relationships within Coffea, most notably for species of Coffea subgenus Coffea occurring in Madagascar. The origin of C. arabica by recent hybridization between C. canephora and C. eugenioides is supported. Phylogenetic separation resulting from the presence of the Dahomey Gap is inferred based on sequence data from Coffea.

Use of fluorescence in situ hybridization as a tool for introgression analysis and chromosome identification in coffee (Coffea arabica L.)

Fluorescence in situ hybridization (FISH) was used to study the presence of alien chromatin in interspecific hybrids and one introgressed line (S.288) derived from crosses between the cultivated species Coffea arabica and the diploid relatives C. canephora and C. liberica. In situ hybridization using genomic DNA from C. canephora and C. arabica as probes showed elevated cross hybridization along the hybrid genome, confirming the weak differentiation between parental genomes. According to our genomic in situ hybridization (GISH) data, the observed genomic resemblance between the modern C. canephora genome (C) and the C. canephora-derived subgenome of C. arabica (Ca) appears rather considerable. Poor discrimination between C and Ca chromosomes supports the idea of low structural modifications of both genomes since the C. arabica speciation, at least in the frequency and distribution of repetitive sequences. GISH was also used to identify alien chromatin segments on chromosome spreads of a C. liberica-introgressed line of C. arabica. Further, use of GISH together with BAC-FISH analysis gave us additional valuable information about the physical localization of the C. liberica fragments carrying the SH3 factor involved in resistance to the coffee leaf rust. Overall, our results illustrate that FISH analysis is a complementary tool for molecular cytogenetic studies in coffee, providing rapid localization of either specific chromosomes or alien chromatin in introgressed genotypes derived from diploid species displaying substantial genomic differentiation from C. arabica.

Differential responses of Coffea arabica L. leaves and roots to chemically induced systemic acquired resistance

Coffea arabica is susceptible to several pests and diseases, some of which affect the leaves and roots. Systemic acquired resistance (SAR) is the main defence mechanism activated in plants in response to pathogen attack. Here, we report the effects of benzo(1,2,3)thiadiazole-7-carbothioic acid-S-methyl ester (BTH), a SAR chemical inducer, on the expression profile of C. arabica. Two cDNA libraries were constructed from the mRNA isolated from leaves and embryonic roots to create 1587 nonredundant expressed sequence tags (ESTs). We developed a cDNA microarray containing 1506 ESTs from the leaves and embryonic roots, and 48 NBS-LRR (nucleotide-binding site leucine-rich repeat) gene fragments derived from 2 specific genomic libraries. Competitive hybridization between untreated and BTH-treated leaves resulted in 55 genes that were significantly overexpressed and 16 genes that were significantly underexpressed. In the roots, 37 and 42 genes were over and underexpressed, respectively. A general shift in metabolism from housekeeping to defence occurred in the leaves and roots after BTH treatment. We observed a systemic increase in pathogenesis-related protein synthesis, in the oxidative burst, and in the cell wall strengthening processes. Moreover, responses in the roots and leaves varied significantly.

Functional characterization of two p-coumaroyl ester 3'-hydroxylase genes from coffee tree: evidence of a candidate for chlorogenic acid biosynthesis

Chlorogenic acid (5-CQA) is one of the major soluble phenolic compounds that is accumulated in coffee green beans. With other hydroxycinnamoyl quinic acids (HQAs), this compound is accumulated in particular in green beans of the cultivated species Coffea canephora. Recent work has indicated that the biosynthesis of 5-CQA can be catalyzed by a cytochrome P450 enzyme, CYP98A3 from Arabidopsis. Two full-length cDNA clones (CYP98A35 and CYP98A36) that encode putative p-coumaroylester 3'-hydroxylases (C3'H) were isolated from C. canephora cDNA libraries. Recombinant protein expression in yeast showed that both metabolized p-coumaroyl shikimate at similar rates, but that only one hydroxylates the chlorogenic acid precursor p-coumaroyl quinate. CYP98A35 appears to be the first C3'H capable of metabolising p-coumaroyl quinate and p-coumaroyl shikimate with the same efficiency. We studied the expression patterns of both genes on 4-month old C. canephora plants and found higher transcript levels in young and in highly vascularized organs for both genes. Gene expression and HQA content seemed to be correlated in these organs. Histolocalization and immunolocalization studies revealed similar tissue localization for caffeoyl quinic acids and p-coumaroylester 3'-hydroxylases. The results indicated that HQA biosynthesis and accumulation occurred mainly in the shoot tip and in the phloem of the vascular bundles. The lack of correlation between gene expression and HQA content observed in some organs is discussed in terms of transport and accumulation mechanisms.

The complete nucleotide sequence of the coffee (Coffea arabica L.) chloroplast genome: organization and implications for biotechnology and phylogenetic relationships amongst angiosperms

The chloroplast genome sequence of Coffea arabica L., the first sequenced member of the fourth largest family of angiosperms, Rubiaceae, is reported. The genome is 155 189 bp in length, including a pair of inverted repeats of 25,943 bp. Of the 130 genes present, 112 are distinct and 18 are duplicated in the inverted repeat. The coding region comprises 79 protein genes, 29 transfer RNA genes, four ribosomal RNA genes and 18 genes containing introns (three with three exons). Repeat analysis revealed five direct and three inverted repeats of 30 bp or longer with a sequence identity of 90% or more. Comparisons of the coffee chloroplast genome with sequenced genomes of the closely related family Solanaceae indicated that coffee has a portion of rps19 duplicated in the inverted repeat and an intact copy of infA. Furthermore, whole-genome comparisons identified large indels (> 500 bp) in several intergenic spacer regions and introns in the Solanaceae, including trnE (UUC)-trnT (GGU) spacer, ycf4-cemA spacer, trnI (GAU) intron and rrn5-trnR (ACG) spacer. Phylogenetic analyses based on the DNA sequences of 61 protein-coding genes for 35 taxa, performed using both maximum parsimony and maximum likelihood methods, strongly supported the monophyly of several major clades of angiosperms, including monocots, eudicots, rosids, asterids, eurosids II, and euasterids I and II. Coffea (Rubiaceae, Gentianales) is only the second order sampled from the euasterid I clade. The availability of the complete chloroplast genome of coffee provides regulatory and intergenic spacer sequences for utilization in chloroplast genetic engineering to improve this important crop.

Molecular characterization of an ethylene receptor gene (CcETR1) in coffee trees, its relationship with fruit development and caffeine content

To understand the importance of ethylene receptor genes in the quality of coffee berries three full-length cDNAs corresponding to a putative ethylene receptor gene (ETR1) were isolated from Coffea canephora cDNA libraries. They differed by their 3'UTR and contained a main ORF and a 5'UTR short ORF putatively encoding a small polypeptide. The CcETR1 gene, present as a single copy in the C. canephora genome, contained five introns in the coding region and one in its 5'UTR. Alternative splicing can occur in C. canephora and C. pseudozanguebariae, leading to a truncated polypeptide. C. pseudozanguebariae ETR1 transcripts showed various forms of splicing alterations. This gene was equally expressed at all stages of fruit development. A segregation study on an inter-specific progeny showed that ETR1 is related to the fructification time, the caffeine content of the green beans, and seed weight. Arabidopsis transformed etiolated seedlings, which over-expressed CcETR1, displayed highly reduced gravitropism, but the triple response was observed in an ethylene enriched environment. These plants behaved like a low-concentration ethylene-insensitive mutant thus confirming the receptor function of the encoded protein. This gene showed no induction during the climacteric crisis but some linkage with traits related to quality.

Introgressive hybridization between the allotetraploid Coffea arabica and one of its diploid ancestors, Coffea canephora, in an exceptional sympatric zone in New Caledonia

The importance of introgressive hybridization in plant evolution has long been recognized. Nevertheless, information on gene flow between allopolyploids and their diploid relatives is very limited, even though gene flow could play a major role in polyploid establishment and evolution. Here, we investigated the processes governing hybrid formation and introgression between the allotetraploid Coffea arabica and one of its ancestral diploid progenitors, C. canephora, in a sympatric zone of New Caledonia. The occurrence of a large assortment of hybridization events between the 2 coffee species is clearly established. First-generation hybrids (F1) and post-F1 hybrids were characterized. The involvement of unreduced gametes of C. canephora is suggested, because tetraploid F1 hybrid plants were detected. Moreover, although bidirectional mating was observed, only unidirectional gene flow from C. canephora to C. arabica was noted in post-F1 hybrids. Most of the collected post-F1 hybrid plants exhibited a high level of introgression, and the frequency of introgression observed among the different analyzed loci was homogeneous, suggesting no significant counterselection against introgressions from C. canephora. Overall, the New Caledonian central mountains appear to be a highly favourable environment for introgressive hybridization and a genetic diversity center for C. arabica.

2D-RNA-coupling numbers: A new computational chemistry approach to link secondary structure topology with biological function

Methods for prediction of proteins, DNA, or RNA function and mapping it onto sequence often rely on bioinformatics alignment approach instead of chemical structure. Consequently, it is interesting to develop computational chemistry approaches based on molecular descriptors. In this sense, many researchers used sequence-coupling numbers and our group extended them to 2D proteins representations. However, no coupling numbers have been reported for 2D-RNA topology graphs, which are highly branched and contain useful information. Here, we use a computational chemistry scheme: (a) transforming sequences into RNA secondary structures, (b) defining and calculating new 2D-RNA-coupling numbers, (c) seek a structure-function model, and (d) map biological function onto the folded RNA. We studied as example 1-aminocyclopropane-1-carboxylic acid (ACC) oxidases known as ACO, which control fruit ripening having importance for biotechnology industry. First, we calculated tau(k)(2D-RNA) values to a set of 90-folded RNAs, including 28 transcripts of ACO and control sequences. Afterwards, we compared the classification performance of 10 different classifiers implemented in the software WEKA. In particular, the logistic equation ACO = 23.8 . tau(1)(2D-RNA) + 41.4 predicts ACOs with 98.9%, 98.0%, and 97.8% of accuracy in training, leave-one-out and 10-fold cross-validation, respectively. Afterwards, with this equation we predict ACO function to a sequence isolated in this work from Coffea arabica (GenBank accession DQ218452). The tau(1)(2D-RNA) also favorably compare with other descriptors. This equation allows us to map the codification of ACO activity on different mRNA topology features. The present computational-chemistry approach is general and could be extended to connect RNA secondary structure topology to other functions.

Identification, characterization and utilization of EST-derived genic microsatellite markers for genome analyses of coffee and related species

Genic microsatellites or EST-SSRs derived from expressed sequence tags (ESTs) are desired because these are inexpensive to develop, represent transcribed genes, and often a putative function can be assigned to them. In this study we investigated 2,553 coffee ESTs (461 from the public domain and 2,092 in-house generated ESTs) for identification and development of genic microsatellite markers. Of these, 2,458 ESTs (all >100 bp in size) were searched for SSRs using MISA--search module followed by stackPACK clustering that revealed a total of 425 microsatellites in 331 (13.5%) non-redundant ESTs/consensus sequences suggesting an approximate frequency of 1 SSR/2.16 kb of the analysed coffee transcriptome. Identified microsatellites mainly comprised of di-/tri-nucleotide repeats, of which repeat motifs AG and AAG were the most abundant. A total of 224 primer pairs could be designed from the non-redundant SSR-positive ESTs (excluding those with only mononucleotide repeats) for possible use as potential genic markers. Of this set, a total of 24 (10%) primer pairs were tested and 18 could be validated as usable markers. Sixteen of these markers revealed moderate to high polymorphism information content (PIC) across 23 genotypes of C. arabica and C. canephora, while 2 markers were found to be monomorphic. All the markers also showed robust cross-species amplifications across 14 Coffea and 4 Psilanthus species. The apparent broad cross-species/genera transferability was further confirmed by cloning and sequencing of the amplified alleles. Thus, the study provides an insight about the frequency and distribution of SSRs in coffee transcriptome, and also demonstrates the successful development of genic-SSRs. It is expected that the potential markers described here would add to the repertoire of DNA markers needed for genetic studies in cultivated coffee and also related taxa that constitute the important secondary genepool for coffee improvement.

Coffee (Coffea sp.)

Coffee (Coffea sp.) is a perennial plant widely cultivated in many tropical countries. It is a cash crop for millions of small farmers in these areas. As compared with other tree species, coffee has long breeding cycles that make conventional breeding programs time consuming. For that matter, genetic transformation can be an effective technique to introduce a desired trait in an already "elite" variety, or to study a gene function and expression. In this chapter, we describe two Agrobacterium-mediated transformation techniques; the first with A. tumefaciens to introduce an insect resistance gene and the second with A. rhizogenes to study candidate gene expression for nematode resistance in transformed roots.