Complex evolutionary history of coffees revealed by full plastid genomes and 28,800 nuclear SNP analyses, with particular emphasis on Coffea canephora (Robusta coffee).

For decades coffees were associated with the genus Coffea. In 2011, the closely related genus Psilanthus was subsumed into Coffea. However, results obtained in 2017-based on 28,800 nuclear SNPs-indicated that there is not substantial phylogenetic support for this incorporation. In addition, a recent study of 16 plastid full-genome sequences highlighted an incongruous placement of Coffea canephora (Robusta coffee) between maternal and nuclear trees. In this study, similar global features of the plastid genomes of Psilanthus and Coffea are observed. In agreement with morphological and physiological traits, the nuclear phylogenetic tree clearly separates Psilanthus from Coffea (with exception to C. rhamnifolia, closer to Psilanthus than to Coffea). In contrast, the maternal molecular tree was incongruent with both morphological and nuclear differentiation, with four main clades observed, two of which include both Psilanthus and Coffea species, and two with either Psilanthus or Coffea species. Interestingly, Coffea and Psilanthus taxa sampled in West and Central Africa are members of the same group. Several mechanisms such as the retention of ancestral polymorphisms due to incomplete lineage sorting, hybridization leading to homoploidy (without chromosome doubling) and alloploidy (for C. arabica) are involved in the evolutionary history of the coffee species. While sharing similar morphological characteristics, the genetic relationships within C. canephora have shown that some populations are well differentiated and genetically isolated. Given the position of its closely-related species, we may also consider C. canephora to be undergoing a long process of speciation with an intermediate step of (sub-)speciation.

Phenotypic diversity assessment within a major ex situ collection of wild endemic coffees in Madagascar.

BACKGROUND AND AIMS: Like other clades, the Coffea genus is highly diversified on the island of Madagascar. The 66 endemic species have colonized various environments and consequently exhibit a wide diversity of morphological, functional and phenological features and reproductive strategies. The trends of interspecific trait variation, which stems from interactions between genetically defined species and their environment, still needed to be addressed for Malagasy coffee trees. METHODS: Data acquisition was done in the most comprehensive ex situ collection of Madagascan wild Coffea. The structure of endemic wild coffees maintained in an ex situ collection was explored in terms of morphological, phenological and functional traits. The environmental (natural habitat) effect was assessed on traits in species from distinct natural habitats. Phylogenetic signal (Pagel's λ, Blomberg's K) was used to quantify trait proximities among species according to their phylogenetic relatedness. KEY RESULTS: Despite the lack of environmental difference in the ex situ collection, widely diverging phenotypes were observed. Phylogenetic signal was found to vary greatly across and even within trait categories. The highest values were exhibited by the ratio of internode mass to leaf mass, the length of the maturation phase and leaf dry matter content (ratio of dry leaf mass to fresh leaf mass). By contrast, traits weakly linked to phylogeny were either constrained by the original natural environment (leaf size) or under selective pressures (phenological traits). CONCLUSIONS: This study gives insight into complex patterns of trait variability found in an ex situ collection, and underlines the opportunities offered by living ex situ collections for research characterizing phenotypic variation.

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

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

Active transposable elements recover species boundaries and geographic structure in Madagascan coffee species.

The completion of the genome assembly for the economically important coffee plant Coffea canephora (Rubiaceae) has allowed the use of bioinformatic tools to identify and characterize a diverse array of transposable elements (TEs), which can be used in evolutionary studies of the genus. An overview of the copy number and location within the C. canephora genome of four TEs is presented. These are tested for their use as molecular markers to unravel the evolutionary history of the Millotii Complex, a group of six wild coffee (Coffea) species native to Madagascar. Two TEs from the Gypsy superfamily successfully recovered some species boundaries and geographic structure among samples, whereas a TE from the Copia superfamily did not. Notably, species occurring in evergreen moist forests of eastern and southeastern Madagascar were divergent with respect to species in other habitats and regions. Our results suggest that the peak of transpositional activity of the Gypsy and Copia TEs occurred, respectively, before and after the speciation events of the tested Madagascan species. We conclude that the utilization of active TEs has considerable potential to unravel the evolutionary history and delimitation of closely related Coffea species. However, the selection of TE needs to be experimentally tested, since each element has its own evolutionary history. Different TEs with similar copy number in a given species can render different dendrograms; thus copy number is not a good selection criterion to attain phylogenetic resolution.

Partial sequencing reveals the transposable element composition of Coffea genomes and provides evidence for distinct evolutionary stories.

The Coffea genus, 124 described species, has a natural distribution spreading from inter-tropical Africa, to Western Indian Ocean Islands, India, Asia and up to Australasia. Two cultivated species, C. arabica and C. canephora, are intensively studied while, the breeding potential and the genome composition of all the wild species remained poorly uncharacterized. Here, we report the characterization and comparison of the highly repeated transposable elements content of 11 Coffea species representatives of the natural biogeographic distribution. A total of 994 Mb from 454 reads were produced with a genome coverage ranging between 3.2 and 15.7 %. The analyses showed that highly repeated transposable elements, mainly LTR retrotransposons (LTR-RT), represent between 32 and 53 % of Coffea genomes depending on their biogeographic location and genome size. Species from West and Central Africa (Eucoffea) contained the highest LTR-RT content but with no strong variation relative to their genome size. At the opposite, for the insular species (Mascarocoffea), a strong variation of LTR-RT was observed suggesting differential dynamics of these elements in this group. Two LTR-RT lineages, SIRE and Del were clearly differentially accumulated between African and insular species, suggesting these lineages were associated to the genome divergence of Coffea species in Africa. Altogether, the information obtained in this study improves our knowledge and brings new data on the composition, the evolution and the divergence of wild Coffea genomes.

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

BACKGROUND AND AIMS: The coffee genus (Coffea) comprises 124 species, and is indigenous to the Old World Tropics. Due to its immense economic importance, Coffea has been the focus of numerous genetic diversity studies, but despite this effort it remains insufficiently studied. In this study the genetic diversity and genetic structure of Coffea across Africa and the Indian Ocean islands is investigated. METHODS: Genetic data were produced using 13 polymorphic nuclear microsatellite markers (simple sequence repeats, SSRs), including seven expressed sequence tag-SSRs, and the data were analysed using model- and non-model-based methods. The study includes a total of 728 individuals from 60 species. KEY RESULTS: Across Africa and the Indian Ocean islands Coffea comprises a closely related group of species with an overall pattern of genotypes running from west to east. Genetic structure was identified in accordance with pre-determined geographical regions and phylogenetic groups. There is a good relationship between morpho-taxonomic species delimitations and genetic units. Genetic diversity in African and Indian Ocean Coffea is high in terms of number of alleles detected, and Madagascar appears to represent a place of significant diversification in terms of allelic richness and species diversity. CONCLUSIONS: Cross-species SSR transferability in African and Indian Ocean islands Coffea was very efficient. On the basis of the number of private alleles, diversification in East Africa and the Indian Ocean islands appears to be more recent than in West and West-Central Africa, although this general trend is complicated in Africa by the position of species belonging to lineages connecting the main geographical regions. The general pattern of phylogeography is not in agreement with an overall east to west (Mascarene, Madagascar, East Africa, West Africa) increase in genome size, the high proportion of shared alleles between the four regions or the high numbers of exclusive shared alleles between pairs or triplets of regions.

A survey of mangiferin and hydroxycinnamic acid ester accumulation in coffee (Coffea) leaves: biological implications and uses.

BACKGROUND AND AIMS: The phenolic composition of Coffea leaves has barely been studied, and therefore this study conducts the first detailed survey, focusing on mangiferin and hydroxycinnamic acid esters (HCEs). METHODS: Using HPLC, including a new technique allowing quantification of feruloylquinic acid together with mangiferin, and histochemical methods, mangiferin content and tissue localization were compared in leaves and fruits of C. pseudozanguebariae, C. arabica and C. canephora. The HCE and mangiferin content of leaves was evaluated for 23 species native to Africa or Madagascar. Using various statistical methods, data were assessed in relation to distribution, ecology, phylogeny and use. KEY RESULTS: Seven of the 23 species accumulated mangiferin in their leaves. Mangiferin leaf-accumulating species also contain mangiferin in the fruits, but only in the outer (sporophytic) parts. In both leaves and fruit, mangiferin accumulation decreases with ageing. A relationship between mangiferin accumulation and UV levels is posited, owing to localization with photosynthetic tissues, and systematic distribution in high altitude clades and species with high altitude representatives. Analyses of mangiferin and HCE content showed that there are significant differences between species, and that samples can be grouped into species, with few exceptions. These data also provide independent support for various Coffea lineages, as proposed by molecular phylogenetic analyses. Sampling of the hybrids C. arabica and C. heterocalyx cf. indicates that mangiferin and HCE accumulation may be under independent parental influence. CONCLUSIONS: This survey of the phenolic composition in Coffea leaves shows that mangiferin and HCE accumulation corresponds to lineage recognition and species delimitation, respectively. Knowledge of the spectrum of phenolic accumulation within species and populations could be of considerable significance for adaptation to specific environments. The potential health benefits of coffee-leaf tea, and beverages and masticatory products made from the fleshy parts of Coffea fruits, are supported by our phenolic quantification.

The absence of the caffeine synthase gene is involved in the naturally decaffeinated status of Coffea humblotiana, a wild species from Comoro archipelago.

Caffeine is the most consumed alkaloid stimulant in the world. It is synthesized through the activity of three known N-methyltransferase proteins. Here we are reporting on the 422-Mb chromosome-level assembly of the Coffea humblotiana genome, a wild and endangered, naturally caffeine-free, species from the Comoro archipelago. We predicted 32,874 genes and anchored 88.7% of the sequence onto the 11 chromosomes. Comparative analyses with the African Robusta coffee genome (C. canephora) revealed an extensive genome conservation, despite an estimated 11 million years of divergence and a broad diversity of genome sizes within the Coffea genus. In this genome, the absence of caffeine is likely due to the absence of the caffeine synthase gene which converts theobromine into caffeine through an illegitimate recombination mechanism. These findings pave the way for further characterization of caffeine-free species in the Coffea genus and will guide research towards naturally-decaffeinated coffee drinks for consumers.

WCSdb: a database of wild Coffea species.

Coffee is a beverage enjoyed by millions of people worldwide and an important commodity for millions of people. Beside the two cultivated species (Coffea arabica and Coffea canephora), the 139 wild coffee species/taxa belonging to the Coffea genus are largely unknown to coffee scientists and breeders although these species may be crucial for future coffee crop development to face climate changes. Here we present the Wild Coffee Species database (WCSdb) hosted by Pl@ntNet platform (http://publish.plantnet-project.org/project/wildcofdb_en), providing information for 141 coffee species/taxa, for which 84 contain a photo gallery and 82 contain sequencing data (genotyping-by-sequencing, chloroplast or whole genome sequences). The objective of this database is to better understand and characterize the species (identification, morphology, biochemical compounds, genetic diversity and sequence data) in order to better protect and promote them. DATABASE URL: http://publish.plantnet-project.org/project/wildcofdb_en.

Isolation and genetic mapping of a Coffea canephora phenylalanine ammonia-lyase gene (CcPAL1) and its involvement in the accumulation of caffeoyl quinic acids.

Biosynthesis of caffeoylquinic acids occurs via the phenylpropanoid pathway in which the phenylalanine ammonia-lyase (PAL) acts as a key-control enzyme. A full-length cDNA (pF6), corresponding to a PAL gene (CcPAL1), was isolated by screening a Coffea canephora fruit cDNA library and its corresponding genomic sequence was characterized. Amplification of total DNA from seven Coffea species revealed differences in intronic length. This interspecific polymorphism was used to locate the gene on a genetic map established for a backcross progeny between Coffea pseudozanguebariae and C. dewevrei. The CcPAL1 gene was found on the same linkage group, but genetically independent, as a caffeoyl-coenzyme A-O-methyltransferase gene, another gene intervening in the phenylpropanoid pathway. In the same backcross, a lower caffeoylquinic acid content was observed in seeds harvested from plants harbouring the C. pseudozanguebariae CcPAL1 allele. Involvement of the CcPAL1 allelic form in the differential accumulation of caffeoylquinic acids in coffee green beans is then discussed.