# The coffee genome hub: a resource for coffee genomes

Alexis Dereeper, Stéphanie Bocs, Mathieu Rouard, Valentin Guignon, Sébastien Ravel, Christine Tranchant-Dubreuil, Valérie Poncet, Olivier Garsmeur, Philippe Lashermes, Gaëtan Droc

## ► To cite this version:

Alexis Dereeper, Stéphanie Bocs, Mathieu Rouard, Valentin Guignon, Sébastien Ravel, et al. The coffee genome hub: a resource for coffee genomes. Nucleic Acids Research, 2015, 43 (D1), pp.D1028-D1035. 10.1093/nar/gku1108. hal-03219674

HAL Id: hal-03219674

https://hal.science/hal-03219674

Submitted on 6 May 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# The coffee genome hub: a resource for coffee genomes

Alexis Dereeper¹,*, Stéphanie Bocs², Mathieu Rouard³, Valentin Guignon³,
Sébastien Ravel¹, Christine Tranchant-Dubreuil⁴, Valérie Poncet⁴, Olivier Garsmeur²,
Philippe Lashermes¹ and Gaëtan Droc²,*

¹UMR Résistance des Plantes aux Bioagresseurs (RPB), Institut de Recherche pour le Développement (IRD), BP 64501, 34394 Montpellier Cedex 5, France, ²UMR Amélioration Génétique et Adaptation des Plantes Méditerranéennes et Tropicales (AGAP), CIRAD, F-34398 Montpellier, France, ³Bioversity International, Parc Scientifique Agropolis II, 34397 Montpellier Cedex 5, France and ⁴UMR Diversité Adaptation et Développement des plantes (DIADE), Institut de Recherche pour le Développement (IRD), BP 64501, 34394 Montpellier Cedex 5, France

Received September 11, 2014; Revised October 20, 2014; Accepted October 23, 2014

## ABSTRACT

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.

## INTRODUCTION

Coffee is the world's most widely traded tropical agricultural commodity. The ability to capture and efficiently use the abundant genetic resources in coffee breeding programs is considered as essential for sustainable coffee production.

Significant advances in our understanding of the coffee genome and its biology must be achieved in the next decades to increase quality, yield and protect the crop from major losses caused by insect pests, diseases and abiotic stress related to climatic change. Unravelling the genetic basis of the traits of interest is therefore a worthwhile goal where genomics can play a prominent role by developing links with breeding programs.

Advances in Next Generation Sequencing (NGS) technologies and the establishment of an international research consortium allowed to complete the project of the first fully sequenced coffee species, *Coffea canephora* (1). It is one of the diploid robusta varieties, which accounts for about 30% of the world's coffee production. *C. canephora* is also one of the parents of *C. arabica*, an allotetraploid derived from hybridization between *C. eugenioides* and *C. canephora*. Anticipating that additional *Coffea* genomes would be sequenced within the next few years, we built a dynamic and scalable crop-specific hub, the Coffee Genome Hub (CGH) which includes the first complete genome of *C. canephora*. These resources will be critical for re-sequencing intra and interspecific genetic resources in coffee. Indeed, community databases federated around a reference genome were initiated in plant genomics with TAIR (2) and Gramene (3) and are proving increasingly important to aggregate, query and retrieve large heterogeneous biological data sets. Various types of plant genomics systems can be distinguished: those based on ENSEMBL system (4) such as Gramene, and others which use and adapt components of the Generic Model Organism Database project (http://www.gmod.org) such as SGN (5) for *Solanaceae* genomes or more recently AIP (6) for Arabidopsis genomes. As previously discussed (7), we opted for GMOD components that are open source, modular, portable and benefiting from a large community support in which we have been involved. Our main strategy

<footnote>*To whom correspondence should be addressed. Tel: +33 4 67 41 61 88; Fax: +33 4 67 41 62 83; Email: alexis.dereeper@ird.fr
Correspondence may also be addressed to Gaëtan Droc. Tel: +33 4 67 61 75 56; Fax: +33 4 67 61 56 05; Email: gaetan.droc@cirad.fr</footnote>
in implementing this hub was to exploit, whenever possible, interconnected generic software solutions to establish a reliable working environment for scientists interested in coffee biology and related topics. The hub architecture is based on the Content Management System (CMS) Drupal with the Tripal module (8) that interacts with the Chado database (9). Together with our Chado controller (10) and Artemis (11), it forms the core of this Community Annotation System (CAS). In fact, similar strategies have been adopted by other information systems (GDR (12), CottonGen (13), Banana Genome Hub (7)) as this facilitates the integration of other GMOD components underlying the hub, such as GBrowse (14), JBrowse (15,16), BioMart (17), Pathway Tools (18), CMAP (19) and Galaxy (20). We also plugged in-house tools developed by the South Green bioinformatics platform (http://www.southgreen.fr) such as GreenPhylDB (21) and SNiPlay (22) (Supplementary Figure S1).

# HUB CONTENT

## Genomics data

*Whole genome sequence data.* The sequence assembly and scaffold anchoring to linkage groups of the genetic map give a total of 569.9 Mb (471.3 without N) divided into 11 pseudomolecules and a sequence made of non-anchored scaffolds randomly concatenated.

*Gene models.* As described by Denoeud et al. (1), automatic gene prediction was performed using the Gaze combiner (23), an integrative gene finding software that combines several evidences such as *ab initio* predictions (Geneid, SNAP and FGenesH), mapping of different protein sequence sets (Blat then Genewise), Expressed Sequence Tags (ESTs) and full-length cDNAs (Blat then Est2genome), as well as RNA-Seq reads from Solexa/Illumina technology (Gmorse) (24). The 25 574 predicted genes obtained from the Gaze combiner have been further annotated to include similarities to proteins of plant model or closely related species, namely tomato or grape, and *in silico* assignment of InterPro protein domains, GO terms and EC (Enzyme Commission) numbers providing information on probable pathways (Table 1).

*Transposable elements.* Transposable elements (TEs) were predicted and classified using the REPET package (25) but also using other tools for LTR retrotransposon (LTR_STRUC), MITEs (MITE hunter) and SINEs. We developed a specific expert procedure (manuscript in preparation) and kept 448,845 TE classified in retrotransposons and DNA transposons superfamilies.

## Transcriptomic data

We collected and inserted the different sources of transcriptomics data that have been provided for protein-coding gene annotation reported by Denoeud et al. (1). This includes unigenes, cDNAs (*C. canephora* and *C. arabica*) as well as RNA-Seq reads from different tissues (root, stamen, pistil, leaf, stem and flower) of different *C. canephora* accessions. For these, clean reads were aligned both on the

reference genome using MapSplice (26) and on the Coding DNA Sequence (CDS) using Burrows-Wheeler Aligner (BWA) (27) for transcript level estimate, normalized expression level (RPKM). In addition, statistical values obtained by other differential expression studies from microarray (28,29) and RNA-Seq (30) experiments were also made available through the CGH. A table summarizing the different sources of cDNAs and RNASeq can be accessed at http://coffee-genome.org/coffeacanephora.

## SNP polymorphisms and genotyping data

A total of 386 560 Single Nucleotide Polymorphisms (SNPs) were identified by Illumina transcriptome sequencing of seven *C. canephora* genotypes, selected to represent the main genetic groups identified in a previous genetic diversity study (31). RNA-Seq reads were aligned directly to the protein-coding sequences using the BWA aligner (27) and SNP discovery was then undertaken with the GATK package (32) using the UnifiedGenotyper module to obtain a list of SNPs and allelic data. Resulting variants were then annotated by SnpEff (33) and provided to biologists in the Hub via JBrowse and SNiPlay (see Managing SNP polymorphisms and genotyping data). Overall, for the seven accessions, 292 125 variants were found to be genotyped without missing data and matching 18 762 genes, including 155 103 non-synonymous mutations (53%).

## Genetic map and molecular markers

Molecular markers and genetic maps are stored in the MoccaDB database (34) and can be viewed and compared using the CMap (19) embedded in the Hub. Four *Coffea* genetic maps are currently recorded, including the high-density *C. canephora* consensus map combining Simple Sequence Repeat (SSR), Restriction Fragment Length Polymorphism (RFLP) and SNP (including RADseq) markers used for the anchoring of scaffolds. For the latter, among the 3230 loci distributed on 11 linkage groups, 2564 markers have been anchored and located on scaffolds and were thus cross-linked between CMap and JBrowse. Additional information (genetic diversity data, primers...) can be accessed from CMap which redirects to MoccaDB.

## Gene families and metabolic pathways

*Gene families.* The Coffee Genome Hub enables comparison of gene families within the *Viridiplantae*. Protein-coding sequences were clustered with 36 other plant species and 24 359 (95%) sequences were classified into 4543 clusters (BLASTP 1e-05 and MCL I = 1.2). Approximately 57% of these clusters were functionally annotated based on the curated catalog of gene families available in GreenPhylDB (21). Best Blast Mutual Hits (BBMH) and phylogenetic analyses were subsequently performed and the resulting approximately 1700 phylogenetic trees as well as homology relationships were made available in the Hub.

Overall, the *Coffea canephora* gene family distribution is consistent compared to the other gene families in plants. Using the InterPro domain distribution tool, 63 transcription factors were studied based the transcription-associated
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Overview of the Coffee Genome Hub (A) A gene search is performed with the of SAM dependent carboxyl methyltransferase (IPR005299) InterPro family identifier. The result page returns a list of genes with graphical display on the chromosomes. (B) The gene report summarizes all the data available for a gene and links to additional resources: (C) Gene family—here the distribution of the Sam Dependent Carboxyl Methyltransferase gene family (GP000195) of coffee illustrates its abundance in plants, (D) JBrowse centered on the region of selected gene (±10 kb) and (E) Pathways tools (e.g. biosynthesis of the caffeine).</figcaption></figure>

protein (TAP) classification rules (35) leading to the identi-
fication of RWP-RK expansion (1). It is interesting to note
that there are 211 coffee-specific clusters, ranging from 2 to
18 paralogs and 116 phylum-specific families in the Aster-
ids sharing at least one sequence in common between C.
canephora, Solanum tuberosum and Solanum lypercosicum.
We identified an over-representation of the NB-ARC super-
family (724 sequences bearing the NB-ARC InterPro signa-
ture IPR002182). In addition, a high number of gene copies
(47 sequences) were also detected for the Sam Dependent
Carboxyl Methyltransferase family that is involved in caf-
feine synthesis (Figure 1).

*Metabolic pathways.* Coffee enzymes and metabolic path-
ways were predicted in the same manner as in the MusaCyc

(7) using respectively PRIAM (36) and Pathway Tools. The percentage of CoffeeCyc enzymes and transporters predicted was 30.6% (7597 enzymes and 231 transporters for a total of 25 574 polypeptides), against 33.5% for AraCyc v18.1 (8884 enzymes and 312 transporters for a total of 27 416), 29.7% for GrapeCyc v4.0 (7559 enzymes and 256 transporters for a total of 26 346) and 24.1% for SolCyc v3.2 (8033 enzymes and 344 transporters for a total of 34 729). The number of pathways was 330 for coffee, 521 for Arabidopsis, 456 for tomato and 485 for grapevine.

Whole genome duplications and synteny blocks

The coffee, tomato and grape proteomes were compared
using BLASTP (e-value 1e-20) and the 10 best hits were
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Transcriptomics data exploration using the Coffee Genome Hub. (A) JBrowse displays alignments of RNA-Seq reads to the genome and allows for each gene a graphical bar representation of RPKM expression values. (B) Heatmap representation of expression values using a user-defined list of genes. (C) Differential expression values (log2ratio, P-value) can be searched by comparison between samples/conditions, and then intersected between studies.</figcaption></figure>

retained. Chromosome segments of genomes containing at least 10 orthologous genes were considered as syntenic regions. A local version of the Plant Genome Duplication Database (37) was implemented in the CGH with a dynamic dot plot allowing the display of syntenic regions and further access to the list of orthologous gene pairs.

The analysis of the paralogous relationships within the reference genome revealed that nearly all chromosome segments were duplicated in three copies. These duplications observed in the coffee genome originate from the ancestral triplication of eudicots (38).

The comparison of coffee and grape genomes showed that globally, one coffee segment corresponds to three grape

segments and that one grape segment corresponds to three coffee segments. The comparison with the tomato genome (which has been subjected to an additional triplication event (39)) showed that one tomato segment can correspond to two or three coffee segments and that one coffee segment can correspond to up to six tomato segments. These comparative analyses with grape and tomato genomes support that no additional Whole Genome Duplication (WGD) events have occurred in the coffee genome following the ancestral triplication of eudicots. All comparisons can be easily visualized at http://coffee-genome.org/syntenic.dotplot. On the CGH dot plots, gene pairs of the synteny blocks are
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Management of SNP polymorphisms in the Coffee Genome Hub. (A) Users can retrieve polymorphic positions based on a subset of genotypes and a subset of genes. The database outputs SNPs together with annotations, minor allele frequency and genotypic data. (B) Connection with JBrowse allows visualizing and browsing the genomic location of the selected SNP. (C) Resulting SNPs can be sent to our tool that calculates and displays the distribution of SNP density on the genome or to a SNP-based distance tree analysis.</figcaption></figure>

<table><caption>**Table 1.** The Hub content. Number of entries by data types as of September 2014</caption><thead><tr><th>Data type</th><th>Number of entries</th></tr></thead><tbody><tr><td>Genome</td><td>1</td></tr><tr><td>Genes</td><td>25 574</td></tr><tr><td>Similarities with proteome of related species (Arabidopsis thaliana, Solanum tuberosum, Vitis vinifera) + Gentianales proteins of Uniprot</td><td>155 283</td></tr><tr><td>Similarities with other Uniprot proteins</td><td>7 498 085</td></tr><tr><td>Similarities with Coffee ESTs</td><td>524 675</td></tr><tr><td>Transposable Elements</td><td>448 845</td></tr><tr><td>Bac Ends sequences</td><td>68 542 (BstYI) + 68 928 (HindIII)</td></tr><tr><td>SSRs</td><td>4 949 134</td></tr><tr><td>SNPs</td><td>386 560</td></tr><tr><td>Anchored genetic markers</td><td>2564</td></tr><tr><td>Expression studies (RNASeq samples)</td><td>10</td></tr><tr><td>Expression studies (microarray samples)</td><td>18</td></tr><tr><td>Gene families</td><td>4543</td></tr><tr><td>Metabolic Pathways</td><td>330</td></tr><tr><td>Synteny relationships (number of syntenic segments)</td><td>566 (Coffee-Coffee)<br>960 (Coffee-Grape)<br>1409 (Coffee-Tomato)</td></tr></tbody></table>
painted in different colors according to the seven ancestral core eudicot chromosomes (1).

## Downloads

Assembly of pseudomolecules as well as their structural and functional annotations are available in FASTA and in Generic File Format (GFF3) formats respectively at http://coffee-genome.org/download.

## TOOLS AND FACILITIES

### Advanced search

Simple or advanced search modes are implemented into the hub. Genes can be searched by keyword, locus, InterPro domain, EC number, location or by Gene Ontology identifier (Figure 1A). The results are displayed both as a dynamic table that summarizes information on the corresponding search and graphically on chromosomes within predicted ancestral blocks as explained in Whole genome duplications and synteny blocks. The output table can also be downloaded as file (FASTA or Excel) or can be sent directly (FASTA file) to Galaxy, wherein further analyses can be performed.

### Primer designer and primer blaster: automatic design and validation of primers

Primer Designer was designed to help users build primers that are specific to intended polymerase chain reaction experiments. Primer3 (40) is used to generate the candidate primer pairs for a given template sequence. Another specific tool, Primer Blaster, was designed to test the specificity of any primer pair on the coffee genome by using BLAST. Users can provide or download primers sequences as fasta files, making sure they are in the same order in both cases. As a result, a table displays all tested primers with primer name and positions, location on chromosome, amplicon size and number of hits on the coffee genome. If the primer pair is really specific, this status is clearly marked as 'ok' in the table.

### Genomic viewers

**JBrowse.** The CGH integrates the interactive Ajax-based genome browser JBrowse v1.11.3 (15,16). We decided to incorporate JBrowse as a complementary browser to the traditional CGI-based genome browser (GBrowse) to speed up the performance of genome browsing and improve interactivity. This next-generation genome browser, built with JavaScript and HTML5, allows users to easily navigate and explore the *C. canephora* genome sequence and annotation data over the web. Anticipating for new genomes or high volume data in the future, we also opted for JBrowse for easy addition of new tracks of information.

Users can select tracks and view various genomic features located on the reference genome, such as gene models, transposable elements and repeats, Blast matches of ESTs, SNPs and putative orthologous genes from other model plant species.

Most of the features displayed (CDSs, TEs, ESTs) are clickable and will link to a window detailing information about the selected feature. Some features are thus directly cross-linked to related specific databases (e.g. Gene report provided by Tripal, markers toward MoccaDB, SNPs toward SNiPlay).

In addition, JBrowse shows RNA-Seq raw read alignments to the genome directly from Binary Alignment/Map (BAM) files. This can be of considerable help in assessing validity of predicted transcripts and for improvement of structural annotation. Furthermore, this functionality allows users to visualize and analyse the distribution of read alignments from one or several samples and to quickly see in what proportion reads are aligned to the genome.

Examples of visualization in JBrowse are shown in Figures 1, 2 and 3.

**Chromosome viewer.** To complement the Genome Browsers, we developed a chromosome viewer based on the Highcharts API (http://api.highcharts.com) that provides a quick overview of feature density along the chromosomes. This viewer includes the visualization of the non-uniform distribution of genes (introns and exons) and transposable elements (Copia, Gypsy, LINEs and DNA transposons), as well as SNP density tracks. Users have the ability to interactively zoom in on regions of interest, and to switch to JBrowse by a simple mouse left click on any point of the graph (a popup displays a link toward JBrowse, centered on the region of the selected genomic position (±20 kb)). Output images can be exported in PNG or SVG formats.

This tool is effective in displaying variation in genome structure and, more generally, any other kind of positional relationships between genomic intervals. It allows a quick preview of feature distribution, comparing it with other related species or highlighting large-scale rearrangements or co-localizations of small-scale events. Such a whole-genome overview can potentially act as an entry point for in-depth investigations. This function is especially favorable in Cofee in which both diploid (*C. canephora*) and allotetraploid (*C. arabica*) species exist, as one could envisage comparing genomic structure between *Coffea* species and thus identify potential structural rearrangements during species evolution. The modular structure renders it versatile: additional genomic features can easily be added by providing a tabular file, thus displaying in a sliding window all the values corresponding to the given feature for each genome interval.

### Exploring gene expression level

Transcriptomics data can be explored at three levels:

(i) The JBrowse feature allows users to view raw read alignments to the genome and display the relative abundance of transcript fragments for the different tissues or conditions (Figure 2A)

(ii) A web form enables users to load a list of genes and to dynamically generate a heatmap image showing the expression values for the different tissues, thus enabling comparisons between conditions (Figure 2B).

(iii) Based on differential expression studies (see Transcriptional data), the CGH allows searching for differen-
tially regulated genes respecting a minimum log 2-
fold expression ratio or statistical *P*-values, by com-
paring two experimental conditions (Figure 2C). Re-
sulting genes can be retrieved as a list or be graphically
summarized depending on their genomic location or
their Gene Ontology (GO) functions. Furthermore, the
database offers the possibility to compare between sev-
eral pairs of conditions either from the same or from
different studies. By clicking ‘intersect’ and selecting an
additional pair of conditions, users may combine the
filtering parameters and retrieve common genes that
have been found differentially expressed (up- or down-
regulated genes) in several analyses. Numbers of genes
are displayed in a Venn diagram.

Managing SNP polymorphisms and genotyping data

It is possible to filter SNPs and INDELs, retrieve
polymorphic positions using a selected subset of
accessions/genotypes and of genes or chromosomes,
and thus discriminate between intra- and inter-genotype
SNPs (Figure 3A). Therefore, using a combination of
successive queries, a customizable subset of variants
respecting specific conditions (Minor Allele Frequency,
non-synonymous effect,% missing data or minimum depth
of coverage) can be generated. Notably, this is an efficient
way to distinguish between homoeoSNP and allelic SNPs,
and thus compare SNPs observed in polyploid species such
as C. arabica to those present in diploid species such as C.
canephora or C. eugenioides, as reported in some of our
studies (41,42).

Users can export genotyping data in various formats (e.g. hmapap, fasta) or send them out for a specialized analysis such as GWAS, diversity analysis or SNP-based distance tree analysis (Figure 3C).

In addition, SNP positions and their associated genotyping information can also be accessed from JBrowse viewer, which can display feature data directly from Variant Call Format (VCF) files (Figure 3B). This allows users to visualize and browse the genomic location of the SNPs. Mutual links between the SNiPlay database and GBrowse/JBrowse are available so that users can easily switch aspects of the investigation to interesting SNP polymorphisms.

FUTURE DIRECTIONS: WHAT ELSE?

The coffee community has now a central portal for the
management of the coffee genome with a comprehensive
set of related data sets. Our main objective is to propose a
community-curated information system with high-quality
gene annotations that will support further genomic stud-
ies. To reach this point, we need to promote synergies with
other groups working on coffee and include new data types
as soon as they become available.

The Coffee Genome Hub is part of a dynamic strategy for data integration and interoperability of bioinformatics applications federated around a community annotation system to ensure the accuracy of genome annotations. In that respect, future directions will rely on JBrowse with the added possibility to manage both private and public data depending on the authenticated user. A synchronization of

the flat files and the genomic features stored in Chado will
be necessary to reflect manual curation of the automatic
structural annotation.

SUPPLEMENTARY DATA

Supplementary Data are available at NAR Online.

ACKNOWLEDGEMENT

The Coffee Genome Hub is supported by the South Green
Bioinformatics Platform (http://www.southgreen.fr/). We
thank all bioinformaticians of the South Green network
for participating to the development or addition of original
tools in the platform. We also thank Bertrand Pitollat who
maintains the computing facilities (web servers and HPC)
with so many system requirements. We also acknowledge
with thanks the Coffee International Research Consortium
to complete the sequencing project of *Coffea canephora*, and
more generally the coffee research community for providing
data and for their feedback. We are grateful to all the peo-
ple that interact regularly with us for the improvements of
the various systems described in this publication and help
us to maintain and further develop useful resources for the
community. Finally, we thank Philippe Chatelet and Rachel
Chase for their help with English editing.

FUNDING

Funding for open access charge: Centre de coopération internationale en recherche agronomique pour le développement; Institut de Recherche pour le Développement.

*Conflict of interest statement.* None declared.

REFERENCES

1. Denoeud,F., Carretero-Paulet,L., Dereeper,A., Droc,G., Guyot,R., Pietrella,M., Zheng,C., Alberti,A., Anthony,F., Aprea,G. et al. (2014) The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. *Science*, **345**, 1181–1184.

2. Lamesch,P., Berardini,T.Z., Li,D., Swarbreck,D., Wilks,C., Sasidharan,R., Muller,R., Dreher,K., Alexander,D.L., Garcia-Hernandez,M. et al. (2012) The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. *Nucleic Acids Res.*, **40**, D1202–D1210.

3. Monaco,M.K., Stein,J., Naithani,S., Wei,S., Dharmawardhana,P., Kumari,S., Amarasinghe,V., Youens-Clark,K., Thomason,J., Preece,J. et al. (2014) Gramene 2013: comparative plant genomics resources. *Nucleic Acids Res.*, **42**, D1193–D1199.

4. Flicek,P., Amode,M.R., Barrell,D., Beal,K., Billis,K., Brent,S., Carvalho-Silva,D., Clapham,P., Coates,G., Fitzgerald,S. et al. (2014) Ensembl 2014. *Nucleic Acids Res.*, **42**, D749–D755.

5. Bombarely,A., Menda,N., Tecle,I.Y., Buels,R.M., Strickler,S., Fischer-York,T., Pujar,A., Leto,J., Gosselin,J., Mueller,L.A. et al. (2011) The Sol Genomics Network (solgenomics.net): growing tomatoes using Perl. *Nucleic Acids Res.*, **39**, D1149–D1155.

6. Consortium,T.I.A.I. (2012) Taking the next step: building an Arabidopsis information portal. *Plant Cell Online*, **24**, 2248–2256.

7. Droc,G., Larivière,D., Guignon,V., Yahiaoui,N., This,D., Garsmeur,O., Dereeper,A., Hamelin,C., Argout,X., Dufayard,J.-F. et al. (2013) The Banana Genome Hub. *Database*, **2013**, bat035.

8. Ficklin,S.P., Sanderson,L.-A., Cheng,C.-H., Staton,M.E., Lee,T., Cho,I.-H., Jung,S., Bett,K.E., Main,D. et al. (2011) Tripal: a construction toolkit for online genome databases. *Database*, **2011**, bar044.
9. Mungall, C.J., Emmert, D.B. and FlyBase Consortium (2007) A Chado case study: an ontology-based modular schema for representing genome-associated biological information. *Bioinformatics*, **23**, i337-i346.

10. Guignon, V., Droc, G., Alaux, M., Baurens, F.-C., Garsmeur, O., Poiron, C., Carver, T., Rouard, M. and Bocs, S. (2012) Chado Controller: advanced annotation management with a community annotation system. *Bioinformatics*, **28**, 1054-1056.

11. Carver, T., Berriman, M., Tivey, A., Patel, C., Böhme, U., Barrell, B.G., Parkhill, J. and Rajandream, M.-A. (2008) Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. *Bioinformatics*, **24**, 2672-2676.

12. Jung, S., Ficklin, S.P., Lee, T., Cheng, C.-H., Blenda, A., Zheng, P., Yu, J., Bombarely, A., Cho, I., Ru, S. et al. (2014) The Genome Database for Rosaceae (GDR): year 10 update. *Nucleic Acids Res.*, **42**, D1237–D1244.

13. Yu, J., Jung, S., Cheng, C.-H., Ficklin, S.P., Lee, T., Zheng, P., Jones, D., Percy, R.G. and Main, D. (2014) CottonGen: a genomics, genetics and breeding database for cotton research. *Nucleic Acids Res.*, **42**, D1229–D1236.

14. Stein, L.D. (2013) Using GBrowse 2.0 to visualize and share next-generation sequence data. *Brief. Bioinform.*, **14**, 162–171.

15. Skinner, M.E., Uzilov, A.V., Stein, L.D., Mungall, C.J. and Holmes, I.H. (2009) JBrowse: a next-generation genome browser. *Genome Res.*, **19**, 1630–1638.

16. Westesson, O., Skinner, M. and Holmes, I. (2013) Visualizing next-generation sequencing data with JBrowse. *Brief. Bioinform.*, **14**, 172–177.

17. Kasprzyk, A. (2011) BioMart: driving a paradigm change in biological data management. *Database*, **2011**, bar049.

18. Karp, P.D., Paley, S.M., Krummenacker, M., Latendresse, M., Dale, J.M., Lee, T.J., Kaipa, P., Gilham, F., Spaulding, A., Popescu, L. et al. (2010) Pathway Tools version 13.0: integrated software for pathway/genome informatics and systems biology. *Brief. Bioinform.*, **11**, 40–79.

19. Youens-Clark, K., Faga, B., Yap, I.V., Stein, L. and Ware, D. (2009) CMap 1.01: a comparative mapping application for the internet. *Bioinformatics*, **25**, 3040–3042.

20. Giardine, B., Riemer, C., Hardison, R.C., Burhans, R., Elnitski, L., Shah, P., Zhang, Y., Blankenberg, D., Albert, I., Taylor, J. et al. (2005) Galaxy: a platform for interactive large-scale genome analysis. *Genome Res.*, **15**, 1451–1455.

21. Rouard, M., Guignon, V., Aluome, C., Laporte, M.-A., Droc, G., Walde, C., Zmaske, C.M., Périn, C. and Conte, M.G. (2011) GreenPhylDB v2.0: comparative and functional genomics in plants. *Nucleic Acids Res.*, **39**, D1095–D1102.

22. Dereeper, A., Nicolas, S., Le Cunff, L., Bacilieri, R., Doligez, A., Peros, J.-P., Ruiz, M. and This, P. (2011) SNiPlay: a web-based tool for detection, management and analysis of SNPs. Application to grapevine diversity projects. *BMC Bioinformatics*, **12**, 134.

23. Howe, K.L., Chothia, T. and Durbin, R. (2002) GAZE: a generic framework for the integration of gene-prediction data by dynamic programming. *Genome Res.*, **12**, 1418–1427.

24. Denoeud, F., Aury, J.-M., Da Silva, C., Noel, B., Rogier, O., Delledonne, M., Morgante, M., Valle, G., Wincker, P., Scarpelli, C. et al. (2008) Annotating genomes with massive-scale RNA sequencing. *Genome Biol.*, **9**, R175.

25. Flutre, T., Duprat, E., Feuillet, C. and Quesneville, H. (2011) Considering transposable element diversification in de novo annotation approaches. *PLoS One*, **6**, e16526.

26. Wang, K., Singh, D., Zeng, Z., Coleman, S.J., Huang, Y., Savich, G.L., He, X., Mieczkowski, P., Grimm, S.A., Perou, C.M. et al. (2010) MapSplice: accurate mapping of RNA-seq reads for splice junction discovery. *Nucleic Acids Res.*, **38**, e178.

27. Li, H. and Durbin, R. (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics*, **25**, 1754–1760.

28. Bardil, A., de Almeida, J.D., Combes, M.C., Lashermes, P. and Bertrand, B. (2011) Genomic expression dominance in the natural allopolyploid Coffea arabica is massively affected by growth temperature. *New Phytol.*, **192**, 760–774.

29. Privat, I., Bardil, A., Gomez, A.B., Severac, D., Dantec, C., Fuentes, I., Mueller, L., Joët, T., Pot, D., Foucrier, S. et al. (2011) The «PUCE CAFE» Project: the first 15K coffee microarray, a new tool for discovering candidate genes correlated to agronomic and quality traits. *BMC Genomics*, **12**, 5.

30. Bertrand, B., Stefan, L., Pirrotta, M., Monchaud, D., Bodio, E., Richard, P., Le Gendre, P., Warmerdam, E., de Jager, M.H., Groothuis, G.M.M. et al. (2014) Caffeine-based gold(I) N-heterocyclic carbenes as possible anticancer agents: synthesis and biological properties. *Inorg. Chem.*, **53**, 2296–2303.

31. Leroy, T., De Bellis, F., Legnate, H., Musoli, P., Kalonji, A., Loor Solórzano, R.G. and Cubry, P. (2014) Developing core collections to optimize the management and the exploitation of diversity of the coffee Coffea canephora. *Genetica*, **142**, 185–199.

32. McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., Garimella, K., Altshuler, D., Gabriel, S., Daly, M. et al. (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.*, **20**, 1297–1303.

33. Cingolani, P., Platts, A., Wang, L.L., Coon, M., Nguyen, T., Wang, L., Land, S.J., Lu, X. and Ruden, D.M. (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. *Fly (Austin)*, **6**, 80–92.

34. Plechakova, O., Tranchant-Dubreuil, C., Benedet, F., Couderc, M., Tinaut, A., Viader, V., Block, P.D., Hamon, P., Campa, C., de Kochko, A. et al. (2009) MoccaDB - an integrative database for functional, comparative and diversity studies in the Rubiaceae family. *BMC Plant Biol.*, **9**, 123.

35. Lang, D., Weiche, B., Timmerhaus, G., Richardt, S., Riano-Pachon, D.M., Correa, L.G.G., Reski, R., Mueller-Roeber, B. and Rensing, S.A. (2010) Genome-wide phylogenetic comparative analysis of plant transcriptional regulation: a timeline of loss, gain, expansion, and correlation with complexity. *Genome Biol. Evol.*, **2**, 488–503.

36. Claudel-Renard, C., Chevalet, C., Faraut, T. and Kahn, D. (2003) Enzyme-specific profiles for genome annotation: PRIAM. *Nucleic Acids Res.*, **31**, 6633–6639.

37. Lee, T.-H., Tang, H., Wang, X. and Paterson, A.H. (2013) PGDD: a database of gene and genome duplication in plants. *Nucleic Acids Res.*, **41**, D1152–D1158.

38. Jiao, Y., Leebens-Mack, J., Ayyampalayam, S., Bowers, J.E., McKain, M.R., McNeal, J., Rolf, M., Ruzicka, D.R., Wafula, E., Wickett, N.J. et al. (2012) A genome triplication associated with early diversification of the core eudicots. *Genome Biol.*, **13**, R3.

39. Consortium, T.T.G. (2012) The tomato genome sequence provides insights into fleshy fruit evolution. *Nature*, **485**, 635–641.

40. Koressaar, T. and Remm, M. (2007) Enhancements and modifications of primer design program Primer3. *Bioinformatics*, **23**, 1289–1291.

41. Combes, M.-C., Dereeper, A., Severac, D., Bertrand, B. and Lashermes, P. (2013) Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures. *New Phytol.*, **200**, 251–260.

42. Lashermes, P., Combes, M.-C., Hueber, Y., Severac, D. and Dereeper, A. (2014) Genome rearrangements derived from homeologous recombination following allopolyploidy speciation in coffee. *Plant J. Cell Mol. Biol.*, **78**, 674–685.