The coffee genome hub: a resource for coffee genomes

Genomic-assisted breeding for climate-smart coffee

Coffee is a universal beverage that drives a multi-industry market on a global basis. Today, the sustainability of coffee production is threatened by accelerated climate changes. In this work, we propose the implementation of genomic-assisted breeding for climate-smart coffee in Coffea canephora. This species is adapted to higher temperatures and is more resilient to biotic and abiotic stresses. After evaluating two populations, over multiple harvests, and under severe drought weather condition, we dissected the genetic architecture of yield, disease resistance, and quality-related traits. By integrating genome-wide association studies and diallel analyses, our contribution is four-fold: (i) we identified a set of molecular markers with major effects associated with disease resistance and post-harvest traits, while yield and plant architecture presented a polygenic background; (ii) we demonstrated the relevance of nonadditive gene actions and projected hybrid vigor when genotypes from different geographically botanical groups are crossed; (iii) we computed medium-to-large heritability values for most of the traits, representing potential for fast genetic progress; and (iv) we provided a first step toward implementing molecular breeding to accelerate improvements in C. canephora. Altogether, this work is a blueprint for how quantitative genetics and genomics can assist coffee breeding and support the supply chain in the face of the current global changes.

Metabolic Pathway Reconstruction Indicates the Presence of Important Medicinal Compounds in Coffea Such as L-DOPA

The use of transcriptomic data to make inferences about plant metabolomes is a useful tool to help the discovery of important compounds in the available biodiversity. To unveil previously undiscovered metabolites of Coffea, of phytotherapeutic and economic value, we employed 24 RNAseq libraries. These libraries were sequenced from leaves exposed to a diverse range of environmental conditions. Subsequently, the data were meticulously processed to create models of putative metabolic networks, which shed light on the production of potential natural compounds of significant interest. Then, we selected one of the predicted compounds, the L-3,4-dihydroxyphenylalanine (L-DOPA), to be analyzed by LC-MS/MS using three biological replicates of flowers, leaves, and fruits from Coffea arabica and Coffea canephora. We were able to identify metabolic pathways responsible for producing several compounds of economic importance. One of the identified pathways involved in isoquinoline alkaloid biosynthesis was found to be active and producing L-DOPA, which is a common product of POLYPHENOL OXIDASES (PPOs, EC 1.14.18.1 and EC 1.10.3.1). We show that coffee plants are a natural source of L-DOPA, a widely used medicine for treatment of the human neurodegenerative condition called Parkinson's disease. In addition, dozens of other compounds with medicinal significance were predicted as potential natural coffee products. By further refining analytical chemistry techniques, it will be possible to enhance the characterization of coffee metabolites, enabling a deeper understanding of their properties and potential applications in medicine.

Genome assemblies and comparison of two Neotropical spiral gingers: Costus pulverulentus and C. lasius

The spiral gingers (Costus L.) are a pantropical genus of herbaceous perennial monocots; the Neotropical clade of Costus radiated rapidly in the past few million years into over 60 species. The Neotropical spiral gingers have a rich history of evolutionary and ecological research that can motivate and inform modern genetic investigations. Here, we present the first 2 chromosome-level genome assemblies in the genus, for C. pulverulentus and C. lasius, and briefly compare their synteny. We assembled the C. pulverulentus genome from a combination of short-read data, Chicago and Dovetail Hi-C chromatin-proximity sequencing, and alignment with a linkage map. We annotated the genome by mapping a C. pulverulentus transcriptome and querying mapped transcripts against a protein database. We assembled the C. lasius genome with Pacific Biosciences HiFi long reads and alignment to the C. pulverulentus genome. These 2 assemblies are the first published genomes for non-cultivated tropical plants. These genomes solidify the spiral gingers as a model system and will facilitate research on the poorly understood genetic basis of tropical plant diversification.

Predicción bioinformática de proteínas NBS-LRR en el genoma de Coffea arabica

Gracias al acceso al genoma completo de Coffea arabica y el Desarrollo de multiples herramientas de bioinformartica que permite la búsqueda de genes de resistencia de plantas (R-genes), ha sido posible implementar estas estrategias en programas de mejora genética. En las plantas, los R-genes codifican proteínas involucradas en mecanismos de defensa contra patógenos. Los genes con dominios tipo Nucleotide-Binding-Site Leucine-Rich-Repeat (NBS-LRR) forman la familia de R-genes de plantas más grande. El objetivo de este estudio fue identificar genes de proteínas NBS-LRR en el genoma de C. arabica utilizando un enfoque bioinformático. Identificamos motivos conservados de R-genes de C. arabica relacionados con genes similares encontrados en Coffea canephora y Coffea eugenoides, dos especies evolutivas relacionadas con C. arabica. Los resultados de estos análisis revelaron proteínas con origen evolutivo provenientes de dicotiledóneo ancestrales, así como proteínas de resistencia específicas del género Coffea. Además, todas las secuencias de los R-genes de C. arabica mostraron una gran similitud con proteína CNL de Arabidopsis thaliana. Finalmente, la presencia de motivos altamente conservados, la distribución cromosómica y las relaciones filogenéticas de los R-genes de C. arabica muestran procesos de coevolución con patógenos adaptados, demostrando de esta manera la importancia del estudio de estos genes en la inmunidad del café. Palabras clave: Café, NBS-LRR, Proteínas de Resistencia, Bioinformática.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Genome-Wide Analysis, Modeling, and Identification of Amino Acid Binding Motifs Suggest the Involvement of GH3 Genes during Somatic Embryogenesis of Coffea canephora

Auxin plays a central role in growth and plant development. To maintain auxin homeostasis, biological processes such as biosynthesis, transport, degradation, and reversible conjugation are essential. The Gretchen Hagen 3 (GH3) family genes codify for the enzymes that esterify indole-3-acetic acid (IAA) to various amino acids, which is a key process in the induction of somatic embryogenesis (SE). The GH3 family is one of the principal families of early response to auxin genes, exhibiting IAA-amido synthetase activity to maintain optimal levels of free auxin in the cell. In this study, we carried out a systematic identification of the GH3 gene family in the genome of Coffea canephora, determining a total of 18 CcGH3 genes. Analysis of the genetic structures and phylogenetic relationships of CcGH3 genes with GH3 genes from other plant species revealed that they could be clustered in two major categories with groups 1 and 2 of the GH3 family of Arabidopsis. We analyzed the transcriptome expression profiles of the 18 CcGH3 genes using RNA-Seq analysis-based data and qRT-PCR during the different points of somatic embryogenesis induction. Furthermore, the endogenous quantification of free and conjugated indole-3-acetic acid (IAA) suggests that the various members of the CcGH3 genes play a crucial role during the embryogenic process of C. canephora. Three-dimensional modeling of the selected CcGH3 proteins showed that they consist of two domains: an extensive N-terminal domain and a smaller C-terminal domain. All proteins analyzed in the present study shared a unique conserved structural topology. Additionally, we identified conserved regions that could function to bind nucleotides and specific amino acids for the conjugation of IAA during SE in C. canephora. These results provide a better understanding of the C. canephora GH3 gene family for further exploration and possible genetic manipulation.

An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis)

Tea, coffee, and cocoa are the three most popular nonalcoholic beverages in the world and have extremely high economic and cultural value. The genomes of four tea plant varieties have recently been sequenced, but there is some debate regarding the characterization of a whole-genome duplication (WGD) event in tea plants. Whether the WGD in the tea plant is shared with other plants in order Ericales and how it contributed to tea plant evolution remained unanswered. Here we re-analyzed the tea plant genome and provided evidence that tea experienced only WGD event after the core-eudicot whole-genome triplication (WGT) event. This WGD was shared by the Polemonioids-Primuloids-Core Ericales (PPC) sections, encompassing at least 17 families in the order Ericales. In addition, our study identified eight pairs of duplicated genes in the catechins biosynthesis pathway, four pairs of duplicated genes in the theanine biosynthesis pathway, and one pair of genes in the caffeine biosynthesis pathway, which were expanded and retained following this WGD. Nearly all these gene pairs were expressed in tea plants, implying the contribution of the WGD. This study shows that in addition to the role of the recent tandem gene duplication in the accumulation of tea flavor-related genes, the WGD may have been another main factor driving the evolution of tea flavor.

Genetic diversity of native and cultivated Ugandan Robusta coffee (Coffea canephora Pierre ex A. Froehner): Climate influences, breeding potential and diversity conservation

Wild genetic resources and their ability to adapt to environmental change are critically important in light of the projected climate change, while constituting the foundation of agricultural sustainability. To address the expected negative effects of climate change on Robusta coffee trees (Coffea canephora), collecting missions were conducted to explore its current native distribution in Uganda over a broad climatic range. Wild material from seven forests could thus be collected. We used 19 microsatellite (SSR) markers to assess genetic diversity and structure of this material as well as material from two ex-situ collections and a feral population. The Ugandan C. canephora diversity was then positioned relative to the species’ global diversity structure. Twenty-two climatic variables were used to explore variations in climatic zones across the sampled forests. Overall, Uganda’s native C. canephora diversity differs from other known genetic groups of this species. In northwestern (NW) Uganda, four distinct genetic clusters were distinguished being from Zoka, Budongo, Itwara and Kibale forests A large southern-central (SC) cluster included Malabigambo, Mabira, and Kalangala forest accessions, as well as feral and cultivated accessions, suggesting similarity in genetic origin and strong gene flow between wild and cultivated compartments. We also confirmed the introduction of Congolese varieties into the SC region where most Robusta coffee production takes place. Identified populations occurred in divergent environmental conditions and 12 environmental variables significantly explained 16.3% of the total allelic variation across populations. The substantial genetic variation within and between Ugandan populations with different climatic envelopes might contain adaptive diversity to cope with climate change. The accessions that we collected have substantially enriched the diversity hosted in the Ugandan collections and thus contribute to ex situ conservation of this vital genetic resource. However, there is an urgent need to develop strategies to enhance complementary in-situ conservation of Coffea canephora in native forests in northwestern Uganda.

The NanDeSyn database for Nannochloropsis systems and synthetic biology

Nannochloropsis species, unicellular industrial oleaginous microalgae, are model organisms for microalgal systems and synthetic biology. To facilitate community-based annotation and mining of the rapidly accumulating functional genomics resources, we have initiated an international consortium and present a comprehensive multi-omics resource database named Nannochloropsis Design and Synthesis (NanDeSyn; http://nandesyn.single-cell.cn). Via the Tripal toolkit, it features user-friendly interfaces hosting genomic resources with gene annotations and transcriptomic and proteomic data for six Nannochloropsis species, including two updated genomes of Nannochloropsis oceanica IMET1 and Nannochloropsis salina CCMP1776. Toolboxes for search, Blast, synteny view, enrichment analysis, metabolic pathway analysis, a genome browser, etc. are also included. In addition, functional validation of genes is indicated based on phenotypes of mutants and relevant bibliography. Furthermore, epigenomic resources are also incorporated, especially for sequencing of small RNAs including microRNAs and circular RNAs. Such comprehensive and integrated landscapes of Nannochloropsis genomics and epigenomics will promote and accelerate community efforts in systems and synthetic biology of these industrially important microalgae.

Genome-Wide Screening and Characterization of Non-Coding RNAs in Coffea canephora

Coffea canephora grains are highly traded commodities worldwide. Non-coding RNAs (ncRNAs) are transcriptional products involved in genome regulation, environmental responses, and plant development. There is not an extensive genome-wide analysis that uncovers the ncRNA portion of the C. canephora genome. This study aimed to provide a curated characterization of six ncRNA classes in the Coffea canephora genome. For this purpose, we employed a combination of similarity-based and structural-based computational approaches with stringent curation. Candidate ncRNA loci had expression evidence analyzed using sRNA-seq libraries. We identified 7455 ncRNA loci (6976 with transcriptional evidence) in the C. canephora genome. This comprised of total 115 snRNAs, 1031 snoRNAs, 92 miRNA precursors, 602 tRNAs, 72 rRNAs, and 5064 lncRNAs. For miRNAs, we identified 159 putative high-confidence targets. This study was the most extensive genomic catalog of curated ncRNAs in the Coffea genus. This data might help elaborating more robust hypotheses in future comparative genomic studies as well as gene regulation and genome dynamics, helping to understand the molecular basis of domestication, environmental adaptation, resistance to pests and diseases, and coffee productivity.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a ‘SyntenyViewer’ to view genome synteny between different cucurbit species and an ‘RNA-Seq’ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Transcriptome analysis of the induction of somatic embryogenesis in Coffea canephora and the participation of ARF and Aux/IAA genes

Background

Somatic embryogenesis (SE) is a useful biotechnological tool to study the morpho-physiological, biochemical and molecular processes during the development of Coffea canephora. Plant growth regulators (PGR) play a key role during cell differentiation in SE. The Auxin-response-factor (ARF) and Auxin/Indole-3-acetic acid (Aux/IAA) are fundamental components involved in the signaling of the IAA. The IAA signaling pathway activates or represses the expression of genes responsive to auxins during the embryogenic transition of the somatic cells. The growing development of new generation sequencing technologies (NGS), as well as bioinformatics tools, has allowed us to broaden the landscape of SE study of various plant species and identify the genes directly involved.

Methods

Analysis of transcriptome expression profiles of the C. canephora genome and the identification of a particular set of differentially expressed genes (DEG) during SE are described in this study.

Results

A total of eight ARF and seven Aux/IAA differentially expressed genes were identified during the different stages of the SE induction process. The quantitative expression analysis showed that ARF18 and ARF5 genes are highly expressed after 21 days of the SE induction, while Aux/IAA7 and Aux/IAA12 genes are repressed.

Discussion

The results of this study allow a better understanding of the genes involved in the auxin signaling pathway as well as their expression profiles during the SE process.

Development and evaluation of a genome-wide Coffee 8.5K SNP array and its application for high-density genetic mapping and for investigating the origin of Coffea arabica L

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single-nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high-density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora-derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high-density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.

Homeologous regulation of Frigida-like genes provides insights on reproductive development and somatic embryogenesis in the allotetraploid Coffea arabica

Coffea arabica is an allotetraploid of high economic importance. C. arabica transcriptome is a combination of the transcripts of two parental genomes (C. eugenioides and C. canephora) that gave rise to the homeologous genes of the species. Previous studies have reported the transcriptional dynamics of C. arabica. In these reports, the ancestry of homeologous genes was identified and the overall regulation of homeologous differential expression (HDE) was explored. One of these genes is part of the FRIGIDA-like family (FRL), which includes the Arabidopsis thaliana flowering-time regulation protein, FRIGIDA (FRI). As nonfunctional FRI proteins give rise to rapid-cycling summer annual ecotypes instead of vernalization-responsive winter-annuals, allelic variation in FRI can modulate flowering time in A. thaliana. Using bioinformatics, genomic analysis, and the evaluation of gene expression of homeologs, we characterized the FRL gene family in C. arabica. Our findings indicate that C. arabica expresses 10 FRL homeologs, and that, throughout flower and fruit development, these genes are differentially transcribed. Strikingly, in addition to confirming the expression of FRL genes during zygotic embryogenesis, we detected FRL expression during direct somatic embryogenesis, a novel finding regarding the FRL gene family. The HDE profile of FRL genes suggests an intertwined homeologous gene regulation. Furthermore, we observed that FLC gene of C. arabica has an expression profile similar to that of CaFRL genes.

Rice Galaxy: an open resource for plant science

BACKGROUND

Rice molecular genetics, breeding, genetic diversity, and allied research (such as rice-pathogen interaction) have adopted sequencing technologies and high-density genotyping platforms for genome variation analysis and gene discovery. Germplasm collections representing rice diversity, improved varieties, and elite breeding materials are accessible through rice gene banks for use in research and breeding, with many having genome sequences and high-density genotype data available. Combining phenotypic and genotypic information on these accessions enables genome-wide association analysis, which is driving quantitative trait loci discovery and molecular marker development. Comparative sequence analyses across quantitative trait loci regions facilitate the discovery of novel alleles. Analyses involving DNA sequences and large genotyping matrices for thousands of samples, however, pose a challenge to non-computer savvy rice researchers.

FINDINGS

The Rice Galaxy resource has shared datasets that include high-density genotypes from the 3,000 Rice Genomes project and sequences with corresponding annotations from 9 published rice genomes. The Rice Galaxy web server and deployment installer includes tools for designing single-nucleotide polymorphism assays, analyzing genome-wide association studies, population diversity, rice-bacterial pathogen diagnostics, and a suite of published genomic prediction methods. A prototype Rice Galaxy compliant to Open Access, Open Data, and Findable, Accessible, Interoperable, and Reproducible principles is also presented.

CONCLUSIONS

Rice Galaxy is a freely available resource that empowers the plant research community to perform state-of-the-art analyses and utilize publicly available big datasets for both fundamental and applied science.

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Glycosyltransferase Family 61 in Liliopsida (Monocot): The Story of a Gene Family Expansion

Plant cell walls play a fundamental role in several plant traits and also influence crop use as livestock nutrition or biofuel production. The Glycosyltransferase family 61 (GT61) is involved in the synthesis of cell wall xylans. In grasses (Poaceae), a copy number expansion was reported for the GT61 family, and raised the question of the evolutionary history of this gene family in a broader taxonomic context. A phylogenetic study was performed on GT61 members from 13 species representing the major angiosperm clades, in order to classify the genes, reconstruct the evolutionary history of this gene family and study its expansion in monocots. Four orthogroups (OG) were identified in angiosperms with two of them displaying a copy number expansion in monocots. These copy number expansions resulted from both tandem and segmental duplications during the genome evolution of monocot lineages. Positive selection footprints were detected on the ancestral branch leading to one of the orthogroups suggesting that the gene number expansion was accompanied by functional diversification, at least partially. We propose an OG-based classification framework for the GT61 genes at different taxonomic levels of the angiosperm useful for any further functional or translational biology study.

Long-read sequencing of the coffee bean transcriptome reveals the diversity of full-length transcripts

Polyploidization contributes to the complexity of gene expression, resulting in numerous related but different transcripts. This study explored the transcriptome diversity and complexity of the tetraploid Arabica coffee (Coffea arabica) bean. Long-read sequencing (LRS) by Pacbio Isoform sequencing (Iso-seq) was used to obtain full-length transcripts without the difficulty and uncertainty of assembly required for reads from short-read technologies. The tetraploid transcriptome was annotated and compared with data from the sub-genome progenitors. Caffeine and sucrose genes were targeted for case analysis. An isoform-level tetraploid coffee bean reference transcriptome with 95 995 distinct transcripts (average 3236 bp) was obtained. A total of 88 715 sequences (92.42%) were annotated with BLASTx against NCBI non-redundant plant proteins, including 34 719 high-quality annotations. Further BLASTn analysis against NCBI non-redundant nucleotide sequences, Coffea canephora coding sequences with UTR, C. arabica ESTs, and Rfam resulted in 1213 sequences without hits, were potential novel genes in coffee. Longer UTRs were captured, especially in the 5΄UTRs, facilitating the identification of upstream open reading frames. The LRS also revealed more and longer transcript variants in key caffeine and sucrose metabolism genes from this polyploid genome. Long sequences (>10 kilo base) were poorly annotated. LRS technology shows the limitation of previous studies. It provides an important tool to produce a reference transcriptome including more of the diversity of full-length transcripts to help understand the biology and support the genetic improvement of polyploid species such as coffee.

Inpactor, Integrated and Parallel Analyzer and Classifier of LTR Retrotransposons and Its Application for Pineapple LTR Retrotransposons Diversity and Dynamics

One particular class of Transposable Elements (TEs), called Long Terminal Repeats (LTRs), retrotransposons, comprises the most abundant mobile elements in plant genomes. Their copy number can vary from several hundreds to up to a few million copies per genome, deeply affecting genome organization and function. The detailed classification of LTR retrotransposons is an essential step to precisely understand their effect at the genome level, but remains challenging in large-sized genomes, requiring the use of optimized bioinformatics tools that can take advantage of supercomputers. Here, we propose a new tool: Inpactor, a parallel and scalable pipeline designed to classify LTR retrotransposons, to identify autonomous and non-autonomous elements, to perform RT-based phylogenetic trees and to analyze their insertion times using High Performance Computing (HPC) techniques. Inpactor was tested on the classification and annotation of LTR retrotransposons in pineapple, a recently-sequenced genome. The pineapple genome assembly comprises 44% of transposable elements, of which 23% were classified as LTR retrotransposons. Exceptionally, 16.4% of the pineapple genome assembly corresponded to only one lineage of the Gypsy superfamily: Del, suggesting that this particular lineage has undergone a significant increase in its copy numbers. As demonstrated for the pineapple genome, Inpactor provides comprehensive data of LTR retrotransposons’ classification and dynamics, allowing a fine understanding of their contribution to genome structure and evolution. Inpactor is available at https://github.com/simonorozcoarias/Inpactor.

MGIS: managing banana (Musa spp.) genetic resources information and high-throughput genotyping data

Unraveling the genetic diversity held in genebanks on a large scale is underway, due to advances in Next-generation sequence (NGS) based technologies that produce high-density genetic markers for a large number of samples at low cost. Genebank users should be in a position to identify and select germplasm from the global genepool based on a combination of passport, genotypic and phenotypic data. To facilitate this, a new generation of information systems is being designed to efficiently handle data and link it with other external resources such as genome or breeding databases. The Musa Germplasm Information System (MGIS), the database for global ex situ-held banana genetic resources, has been developed to address those needs in a user-friendly way. In developing MGIS, we selected a generic database schema (Chado), the robust content management system Drupal for the user interface, and Tripal, a set of Drupal modules which links the Chado schema to Drupal. MGIS allows germplasm collection examination, accession browsing, advanced search functions, and germplasm orders. Additionally, we developed unique graphical interfaces to compare accessions and to explore them based on their taxonomic information. Accession-based data has been enriched with publications, genotyping studies and associated genotyping datasets reporting on germplasm use. Finally, an interoperability layer has been implemented to facilitate the link with complementary databases like the Banana Genome Hub and the MusaBase breeding database.

Database URL:https://www.crop-diversity.org/mgis/

Genome-wide association study reveals candidate genes influencing lipids and diterpenes contents in Coffea arabica L

Lipids, including the diterpenes cafestol and kahweol, are key compounds that contribute to the quality of coffee beverages. We determined total lipid content and cafestol and kahweol concentrations in green beans and genotyped 107 Coffea arabica accessions, including wild genotypes from the historical FAO collection from Ethiopia. A genome-wide association study was performed to identify genomic regions associated with lipid, cafestol and kahweol contents and cafestol/kahweol ratio. Using the diploid Coffea canephora genome as a reference, we identified 6,696 SNPs. Population structure analyses suggested the presence of two to three groups (K = 2 and K = 3) corresponding to the east and west sides of the Great Rift Valley and an additional group formed by wild accessions collected in western forests. We identified 5 SNPs associated with lipid content, 4 with cafestol, 3 with kahweol and 9 with cafestol/kahweol ratio. Most of these SNPs are located inside or near candidate genes related to metabolic pathways of these chemical compounds in coffee beans. In addition, three trait-associated SNPs showed evidence of directional selection among cultivated and wild coffee accessions. Our results also confirm a great allelic richness in wild accessions from Ethiopia, especially in accessions originating from forests in the west side of the Great Rift Valley.

Shaping the evolutionary tree of green plants: evidence from the GST family

Glutathione-S-transferases (GSTs) are encoded by genes belonging to a wide ubiquitous family in aerobic species and catalyze the conjugation of electrophilic substrates to glutathione (GSH). GSTs are divided in different classes, both in plants and animals. In plants, GSTs function in several pathways, including those related to secondary metabolites biosynthesis, hormone homeostasis, defense from pathogens and allow the prevention and detoxification of damage from heavy metals and herbicides. 1107 GST protein sequences from 20 different plant species with sequenced genomes were analyzed. Our analysis assigns 666 unclassified GSTs proteins to specific classes, remarking the wide heterogeneity of this gene family. Moreover, we highlighted the presence of further subclasses within each class. Regarding the class GST-Tau, one possible subclass appears to be present in all the Tau members of ancestor plant species. Moreover, the results highlight the presence of members of the Tau class in Marchantiophytes and confirm previous observations on the absence of GST-Tau in Bryophytes and green algae. These results support the hypothesis regarding the paraphyletic origin of Bryophytes, but also suggest that Marchantiophytes may be on the same branch leading to superior plants, depicting an alternative model for green plants evolution.

Transcripts of pectin-degrading enzymes and isolation of complete cDNA sequence of a pectate lyase gene induced by coffee white stem borer (Xylotrechus quadripes) in the bark tissue of Coffea canephora (robusta coffee)

Of the two commercially cultivated coffee (Coffea) species, C. arabica (arabica) is highly susceptible and C. canephora (robusta) is highly resistant to the insect pest Xylotrechus quadripes (Coleoptera: Cerambycidae), commonly known as coffee white stem borer (CWSB). We constructed a forward-subtracted cDNA library by Suppression Subtractive Hybridization (SSH) from robusta bark tissue for profiling genes induced by CWSB infestation. Among the 265 unigenes of the SSH EST library, 7 unigenes (5 contigs and 2 singletons) matching different pectin-degrading enzymes were discovered. These ESTs matched one pectate lyase, three polygalacturonases, and one pectin acetylesterase gene. Quantitative real-time PCR (qRT-PCR) revealed that CWSB infestation strongly induces the pectate lyase gene at 72 h. Complete cDNA sequence of the pectate lyase gene was obtained through 3' and 5' RACE reactions. It was a 1595 bp long sequence that included full CDS and both UTRs. Against C. canephora genome sequences in Coffee Genome Hub database ( http://coffee-genome.org/ ), it had 22 matches to different pectate lyase genes mapped on 9 of the 11 pseudochromosomes, the top match being Cc07_g00190 Pectate lyase. In NCBI database, it matched pectate lyase sequences of several plants. Apart from C. canephora, the closest pectate lyase matches were from Sesamum indicum and Nicotiana tabacum. The pectinolytic enzymes discovered here are thought to play a role in the production of oligogalacturonides (OGs) which act as Damage-Associated Molecular Pattern (DAMP) signals eliciting innate immunity in plants. The pectate lyase gene, induced by CWSB infestation, along with other endogenous pectinolytic enzymes and CWSB-specific elicitors, may be involved in triggering basal defense responses to protect the CWSB-damaged tissue against pathogens, as well as to contain CWSB in robusta.

Genome-wide identification, classification and transcriptional analysis of nitrate and ammonium transporters in Coffea

Nitrogen (N) is quantitatively the main nutrient required by coffee plants, with acquisition mainly by the roots and mostly exported to coffee beans. Nitrate (NO3-) and ammonium (NH4+) are the most important inorganic sources for N uptake. Several N transporters encoded by different gene families mediate the uptake of these compounds. They have an important role in source preference for N uptake in the root system. In this study, we performed a genome-wide analysis, including in silico expression and phylogenetic analyses of AMT1, AMT2, NRT1/PTR, and NRT2 transporters in the recently sequenced Coffea canephora genome. We analyzed the expression of six selected transporters in Coffea arabica roots submitted to N deficiency. N source preference was also analyzed in C. arabica using isotopes. C. canephora N transporters follow the patterns observed for most eudicots, where each member of the AMT and NRT families has a particular role in N mobilization, and where some of these are modulated by N deficiency. Despite the prevalence of putative nitrate transporters in the Coffea genome, ammonium was the preferential inorganic N source for N-starved C. arabica roots. This data provides an important basis for fundamental and applied studies to depict molecular mechanisms involved in N uptake in coffee trees.

Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis

Coffea arabica L. is an important crop in several developing countries. Despite its economic importance, minimal transcriptome data are available for fruit tissues, especially during fruit development where several compounds related to coffee quality are produced. To understand the molecular aspects related to coffee fruit and grain development, we report a large-scale transcriptome analysis of leaf, flower and perisperm fruit tissue development. Illumina sequencing yielded 41,881,572 high-quality filtered reads. De novo assembly generated 65,364 unigenes with an average length of 1,264 bp. A total of 24,548 unigenes were annotated as protein coding genes, including 12,560 full-length sequences. In the annotation process, we identified nine candidate genes related to the biosynthesis of raffinose family oligossacarides (RFOs). These sugars confer osmoprotection and are accumulated during initial fruit development. Four genes from this pathway had their transcriptional pattern validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Furthermore, we identified ~24,000 putative target sites for microRNAs (miRNAs) and 134 putative transcriptionally active transposable elements (TE) sequences in our dataset. This C. arabica transcriptomic atlas provides an important step for identifying candidate genes related to several coffee metabolic pathways, especially those related to fruit chemical composition and therefore beverage quality. Our results are the starting point for enhancing our knowledge about the coffee genes that are transcribed during the flowering and initial fruit development stages.

Transcriptome Analysis of Hamelia patens (Rubiaceae) Anthers Reveals Candidate Genes for Tapetum and Pollen Wall Development

Studies of the anther transcriptome on non-model plants without a known genome are surprisingly scarce. RNA-Seq and digital gene expression (DGE) profiling provides a comprehensive approach to identify candidate genes contributing to developmental processes in non-model species. Here we built a transcriptome library of developing anthers of Hamelia patens and analyzed DGE profiles from each stage to identify genes that regulate tapetum and pollen development. In total 7,720 putative differentially expressed genes across four anther stages were identified. The number of putative stage-specific genes was: 776 at microspore mother cell stage, 807 at tetrad stage, 322 at uninucleate microspore stage, and the highest number (1,864) at bicellular pollen stage. GO enrichment analysis revealed 243 differentially expressed and 108 stage-specific genes that are potentially related to tapetum development, sporopollenin synthesis, and pollen wall. The number of expressed genes, their function and expression profiles were all significantly correlated with anther developmental processes. Overall comparisons of anther and pollen transcriptomes with those of rice and Arabidopsis together with the expression profiles of homologs of known anther-expressed genes, revealed conserved patterns and also divergence. The divergence may reflect taxon-specific differences in gene expression, the use RNA-seq as a more sensitive methodology, variation in tissue composition and sampling strategies. Given the lack of genomic sequence, this study succeeded in assigning putative identity to a significant proportion of anther-expressed genes and genes relevant to tapetum and pollen development in H. patens. The anther transcriptome revealed a molecular distinction between developmental stages, serving as a resource to unravel the functions of genes involved in anther development in H. patens and informing the analysis of other members of the Rubiaceae.

Evolutionary Analyses of GRAS Transcription Factors in Angiosperms

GRAS transcription factors (TFs) play critical roles in plant growth and development such as gibberellin and mycorrhizal signaling. Proteins belonging to this gene family contain a typical GRAS domain in the C-terminal sequence, whereas the N-terminal region is highly variable. Although, GRAS genes have been characterized in a number of plant species, their classification is still not completely resolved. Based on a panel of eight representative species of angiosperms, we identified 29 orthologous groups or orthogroups (OGs) for the GRAS gene family, suggesting that at least 29 ancestor genes were present in the angiosperm lineage before the "Amborella" evolutionary split. Interestingly, some taxonomic groups were missing members of one or more OGs. The gene number expansion usually observed in transcription factors was not observed in GRAS while the genome triplication ancestral to the eudicots (γ hexaploidization event) was detectable in a limited number of GRAS orthogroups. We also found conserved OG-specific motifs in the variable N-terminal region. Finally, we could regroup OGs in 17 subfamilies for which names were homogenized based on a literature review and described 5 new subfamilies (DLT, RAD1, RAM1, SCLA, and SCLB). This study establishes a consistent framework for the classification of GRAS members in angiosperm species, and thereby a tool to correctly establish the orthologous relationships of GRAS genes in most of the food crops in order to facilitate any subsequent functional analyses in the GRAS gene family. The multi-fasta file containing all the sequences used in our study could be used as database to perform diagnostic BLASTp to classify GRAS genes from other non-model species.

Advances in genomics for the improvement of quality in coffee

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

Shift in precipitation regime promotes interspecific hybridization of introduced Coffea species

The frequency of plant species introductions has increased in a highly connected world, modifying species distribution patterns to include areas outside their natural ranges. These introductions provide the opportunity to gain new insight into the importance of flowering phenology as a component of adaptation to a new environment. Three Coffea species, C. arabica, C. canephora (Robusta), and C. liberica, native to intertropical Africa have been introduced to New Caledonia. On this archipelago, a secondary contact zone has been characterized where these species coexist, persist, and hybridize spontaneously. We investigated the impact of environmental changes undergone by each species following its introduction in New Caledonia on flowering phenology and overcoming reproductive barriers between sister species. We developed species distribution models and compared both environmental envelopes and climatic niches between native and introduced hybrid zones. Flowering phenology was monitored in a population in the hybrid zone along with temperature and precipitation sequences recorded at a nearby weather station. The extent and nature of hybridization events were characterized using chloroplast and nuclear microsatellite markers. The three Coffea species encountered weak environmental suitability compared to their native ranges when introduced to New Caledonia, especially C. arabica and C. canephora. The niche of the New Caledonia hybrid zone was significantly different from all three species' native niches based on identity tests (I Similarity and D Schoener's Similarity Indexes). This area appeared to exhibit intermediate conditions between the native conditions of the three species for temperature-related variables and divergent conditions for precipitation-related ones. Flowering pattern in these Coffea species was shown to have a strong genetic component that determined the time between the triggering rain and anthesis (flower opening), specific to each species. However, a precipitation regime different from those in Africa was directly involved in generating partial flowering overlap between species and thus in allowing hybridization and interspecific gene flow. Interspecific hybrids accounted for 4% of the mature individuals in the sympatric population and occurred between each pair of species with various level of introgression. Adaptation to new environmental conditions following introduction of Coffea species to New Caledonia has resulted in a secondary contact between three related species, which would not have happened in their native ranges, leading to hybridization and gene flow.

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

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

Identification of Putative Molecular Markers Associated with Root Traits in Coffea canephora Pierre ex Froehner

Coffea canephora exhibit poor root system and are very sensitive to drought stress that affects growth and production. Deeper root system has been largely empirical as better avoidance to soil water limitation in drought condition. The present study aimed to identify molecular markers linked to high root types in Coffea canephora using molecular markers. Contrasting parents, L1 valley with low root and S.3334 with high root type, were crossed, and 134 F1 individuals were phenotyped for root and associated physiological traits (29 traits) and genotyped with 41 of the 320 RAPD and 9 of the 55 SSR polymorphic primers. Single marker analysis was deployed for detecting the association of markers linked to root associated traits by SAS software. There were 13 putative RAPD markers associated with root traits such as root length, secondary roots, root dry weight, and root to shoot ratio, in which root length associated marker OPS1850 showed high phenotypic variance of 6.86%. Two microsatellite markers linked to root length (CPCM13400) and root to shoot ratio (CM211300). Besides, 25 markers were associated with more than one trait and few of the markers were associated with positively related physiological traits and can be used in marker assisted trait selection.