Haitian coffee agroforestry systems harbor complex arabica variety mixtures and under-recognized genetic diversity.

Though facing significant challenges, coffee (Coffea arabica) grown in Haitian agroforestry systems are important contributors to rural livelihoods and provide several ecosystem services. However, little is known about their genetic diversity and the variety mixtures used. In light of this, there is a need to characterize Haitian coffee diversity to help inform revitalization of this sector. We sampled 28 diverse farms in historically important coffee growing regions of northern and southern Haiti. We performed KASP-genotyping of SNP markers and HiPlex multiplex amplicon sequencing for haplotype calling on our samples, as well as several Ethiopian and commercial accessions from international collections. This allowed us to assign Haitian samples to varietal groups. Our analyses revealed considerable genetic diversity in Haitian farms, higher in fact than many farmers realized. Notably, genetic structure analyses revealed the presence of clusters related to Typica, Bourbon, and Catimor groups, another group that was not represented in our reference accession panel, and several admixed individuals. Across the study areas, we found both mixed-variety farms and monovarietal farms with the historical and traditional Typica variety. This study is, to our knowledge, the first to genetically characterize Haitian C. arabica variety mixtures, and report the limited cultivation of C. canephora (Robusta coffee) in the study area. Our results show that some coffee farms are repositories of historical, widely-abandoned varieties while others are generators of new diversity through genetic mixing.

Genome-Wide Admixture Mapping Identifies Wild Ancestry-of-Origin Segments in Cultivated Robusta Coffee.

Humans have had a major influence on the dissemination of crops beyond their native range, thereby offering new hybridization opportunities. Characterizing admixed genomes with mosaic origins generates valuable insight into the adaptive history of crops and the impact on current varietal diversity. We applied the ELAI tool-an efficient local ancestry inference method based on a two-layer hidden Markov model to track segments of wild origin in cultivated accessions in the case of multiway admixtures. Source populations-which may actually be limited and partially admixed-must be generally specified when using such inference models. We thus developed a framework to identify local ancestry with admixed source populations. Using sequencing data for wild and cultivated Coffea canephora (commonly called Robusta), our approach was found to be highly efficient and accurate on simulated hybrids. Application of the method to assess elite Robusta varieties from Vietnam led to the identification of an accession derived from a likely backcross between two genetic groups from the Congo Basin and the western coastal region of Central Africa. Admixtures resulting from crop hybridization and diffusion could thus lead to the generation of elite high-yielding varieties. Our methods should be widely applicable to gain insight into the role of hybridization during plant and animal evolutionary history.

Intraspecific variation in growth response to drought stress across geographic locations and genetic groups in Coffea canephora.

Uganda lies within the drier end of the natural distribution range of Coffea canephora and contains unexplored genetic material that could be drought-adapted and useful for developing climate-resilient varieties. Using water treatment: (i) ample and (ii) restricted-water, the response of 148 genotypes were studied comprising wild, feral and cultivated C. canephora. Biomass allocation, standing leaf area and leaf area growth data were collected. Linear mixed effect models and PCA were used to the analyze effect of water treatment on genotypes from different: (i) cultivation status, (ii) genetic groups and (iii) locations. We also assessed the relationship between drought tolerance for relative growth rate in leaf area (RGRA), total number of leaves (TNL), total leaf area (TLA) and total leaf dry weight (TLDW) of genotypes at final harvest. Restricted-water reduced RGRA across genetic groups (3.2-32.5%) and locations (7.1-36.7%) but not cultivation status. For TNL, TLA and TLDW, genotypes that performed well in ample-water performed worse under restricted-water, indicating growth-tolerance trade-off. Drought tolerance in RGRA and TNL were negatively correlated with wetness index suggesting some degree of adaptation to local climate. Findings indicate a growth-tolerance trade-off within this tropical tree species and drought tolerance of Uganda's C. canephora is somewhat associated with local climate.

Ecological and genomic vulnerability to climate change across native populations of Robusta coffee (Coffea canephora).

The assessment of population vulnerability under climate change is crucial for planning conservation as well as for ensuring food security. Coffea canephora is, in its native habitat, an understorey tree that is mainly distributed in the lowland rainforests of tropical Africa. Also known as Robusta, its commercial value constitutes a significant revenue for many human populations in tropical countries. Comparing ecological and genomic vulnerabilities within the species' native range can provide valuable insights about habitat loss and the species' adaptive potential, allowing to identify genotypes that may act as a resource for varietal improvement. By applying species distribution models, we assessed ecological vulnerability as the decrease in climatic suitability under future climatic conditions from 492 occurrences. We then quantified genomic vulnerability (or risk of maladaptation) as the allelic composition change required to keep pace with predicted climate change. Genomic vulnerability was estimated from genomic environmental correlations throughout the native range. Suitable habitat was predicted to diminish to half its size by 2050, with populations near coastlines and around the Congo River being the most vulnerable. Whole-genome sequencing revealed 165 candidate SNPs associated with climatic adaptation in C. canephora, which were located in genes involved in plant response to biotic and abiotic stressors. Genomic vulnerability was higher for populations in West Africa and in the region at the border between DRC and Uganda. Despite an overall low correlation between genomic and ecological vulnerability at broad scale, these two components of vulnerability overlap spatially in ways that may become damaging. Genomic vulnerability was estimated to be 23% higher in populations where habitat will be lost in 2050 compared to regions where habitat will remain suitable. These results highlight how ecological and genomic vulnerabilities are relevant when planning on how to cope with climate change regarding an economically important species.

Adaptive potential of Coffea canephora from Uganda in response to climate change.

Understanding vulnerabilities of plant populations to climate change could help preserve their biodiversity and reveal new elite parents for future breeding programmes. To this end, landscape genomics is a useful approach for assessing putative adaptations to future climatic conditions, especially in long-lived species such as trees. We conducted a population genomics study of 207 Coffea canephora trees from seven forests along different climate gradients in Uganda. For this, we sequenced 323 candidate genes involved in key metabolic and defence pathways in coffee. Seventy-one single nucleotide polymorphisms (SNPs) were found to be significantly associated with bioclimatic variables, and were thereby considered as putatively adaptive loci. These SNPs were linked to key candidate genes, including transcription factors, like DREB-like and MYB family genes controlling plant responses to abiotic stresses, as well as other genes of organoleptic interest, such as the DXMT gene involved in caffeine biosynthesis and a putative pest repellent. These climate-associated genetic markers were used to compute genetic offsets, predicting population responses to future climatic conditions based on local climate change forecasts. Using these measures of maladaptation to future conditions, substantial levels of genetic differentiation between present and future diversity were estimated for all populations and scenarios considered. The populations from the forests Zoka and Budongo, in the northernmost zone of Uganda, appeared to have the lowest genetic offsets under all predicted climate change patterns, while populations from Kalangala and Mabira, in the Lake Victoria region, exhibited the highest genetic offsets. The potential of these findings in terms of ex situ conservation strategies are discussed.

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.

McSwan: A joint site frequency spectrum method to detect and date selective sweeps across multiple population genomes.

Inferring the mode and tempo of natural selection helps further our understanding of adaptation to past environmental changes. Here, we introduce McSwan, a method to detect and date past and recent natural selection events in the case of a hard sweep. The method is based on the comparison of site frequency spectra obtained under various demographic models that include selection. McSwan demonstrated high power (high sensitivity and specificity) in capturing hard selective sweep events without requiring haplotype phasing. It performed slightly better than SweeD when the recent effective population size was low and the genomic region was small. We then applied our method to a European (CEU) and an African (LWK) human re-sequencing data set. Most hard sweeps were detected in the CEU population (96%). Moreover, hard sweeps in the African population were estimated to have occurred further back in time (mode: 43,625 years BP) compared to those of Europeans (mode: 24,850 years BP). Most of the estimated ages of hard sweeps in Europeans were associated with the Last Glacial Maximum and were enriched in immunity-associated genes.

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.

Two disjunct Pleistocene populations and anisotropic postglacial expansion shaped the current genetic structure of the relict plant Amborella trichopoda.

Past climate fluctuations shaped the population dynamics of organisms in space and time, and have impacted their present intra-specific genetic structure. Demo-genetic modelling allows inferring the way past demographic and migration dynamics have determined this structure. Amborella trichopoda is an emblematic relict plant endemic to New Caledonia, widely distributed in the understory of non-ultramafic rainforests. We assessed the influence of the last glacial climates on the demographic history and the paleo-distribution of 12 Amborella populations covering the whole current distribution. We performed coalescent genetic modelling of these dynamics, based on both whole-genome resequencing and microsatellite genotyping data. We found that the two main genetic groups of Amborella were shaped by the divergence of two ancestral populations during the last glacial maximum. From 12,800 years BP, the South ancestral population has expanded 6.3-fold while the size of the North population has remained stable. Recent asymmetric gene flow between the groups further contributed to the phylogeographical pattern. Spatially explicit coalescent modelling allowed us to estimate the location of ancestral populations with good accuracy (< 22 km) and provided indications regarding the mid-elevation pathways that facilitated post-glacial expansion.

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 Amborella vacuolar processing enzyme family.

Most vacuolar proteins are synthesized on rough endoplasmic reticulum as proprotein precursors and then transported to the vacuoles, where they are converted into their respective mature forms by vacuolar processing enzymes (VPEs). In the case of the seed storage proteins, this process is of major importance, as it conditions the establishment of vigorous seedlings. Toward the goal of identifying proteome signatures that could be associated with the origin and early diversification of angiosperms, we previously characterized the 11S-legumin-type seed storage proteins from Amborella trichopoda, a rainforest shrub endemic to New Caledonia that is also the probable sister to all other angiosperms (Amborella Genome Project, 2013). In the present study, proteomic and genomic approaches were used to characterize the VPE family in this species. Three genes were found to encode VPEs in the Amborella's genome. Phylogenetic analyses showed that the Amborella sequences grouped within two major clades of angiosperm VPEs, indicating that the duplication that generated the ancestors of these clades occurred before the most recent common ancestor of living angiosperms. A further important duplication within the VPE family appears to have occurred in common ancestor of the core eudicots, while many more recent duplications have also occurred in specific taxa, including both Arabidopsis thaliana and Amborella. An analysis of natural genetic variation for each of the three Amborella VPE genes revealed the absence of selective forces acting on intronic and exonic single-nucleotide polymorphisms among several natural Amborella populations in New Caledonia.

The coffee genome hub: a resource for coffee genomes.

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

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis.

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia.

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range-wide genetic diversity based on nuclear microsatellite data and detected four well-differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21,000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra- and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.

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.

Microsatellite markers for Amborella (Amborellaceae), a monotypic genus endemic to New Caledonia.

PREMISE OF THE STUDY: Informative markers are required for assessing the diversity of Amborella trichopoda, the only species of its order, endemic to New Caledonia and considered to be the sister species to all flowering plants. Therefore, expressed sequence tag (EST)-based microsatellite markers were developed. • METHODS AND RESULTS: Fifty-five microsatellite loci were characterized in 14896 putative unigenes, which were generated by assembling A. trichopoda ESTs from the public sequence database. Seventeen markers revealed polymorphism in 80 adult shrubs from three populations. The number of alleles per locus ranged from two to 12, with a total of 132 alleles scored. The mean expected heterozygosity per population ranged from 0.336 to 0.567. • CONCLUSIONS: These markers offer an appropriate amount of variation to investigate genetic diversity structure, gene flow, and other conservation issues.

Ancestral synteny shared between distantly-related plant species from the asterid (Coffea canephora and Solanum Sp.) and rosid (Vitis vinifera) clades.

BACKGROUND: Coffee trees (Rubiaceae) and tomato (Solanaceae) belong to the Asterid clade, while grapevine (Vitaceae) belongs to the Rosid clade. Coffee and tomato separated from grapevine 125 million years ago, while coffee and tomato diverged 83-89 million years ago. These long periods of divergent evolution should have permitted the genomes to reorganize significantly. So far, very few comparative mappings have been performed between very distantly related species belonging to different clades. We report the first multiple comparison between species from Asterid and Rosid clades, to examine both macro-and microsynteny relationships. RESULTS: Thanks to a set of 867 COSII markers, macrosynteny was detected between coffee, tomato and grapevine. While coffee and tomato genomes share 318 orthologous markers and 27 conserved syntenic segments (CSSs), coffee and grapevine also share a similar number of syntenic markers and CSSs: 299 and 29 respectively. Despite large genome macrostructure reorganization, several large chromosome segments showed outstanding macrosynteny shedding new insights into chromosome evolution between Asterids and Rosids. We also analyzed a sequence of 174 kb containing the ovate gene, conserved in a syntenic block between coffee, tomato and grapevine that showed a high-level of microstructure conservation. A higher level of conservation was observed between coffee and grapevine, both woody and long life-cycle plants, than between coffee and tomato. Out of 16 coffee genes of this syntenic segment, 7 and 14 showed complete synteny between coffee and tomato or grapevine, respectively. CONCLUSIONS: These results show that significant conservation is found between distantly related species from the Asterid (Coffea canephora and Solanum sp.) and Rosid (Vitis vinifera) clades, at the genome macrostructure and microstructure levels. At the ovate locus, conservation did not decline in relation to increasing phylogenetic distance, suggesting that the time factor alone does not explain divergences. Our results are considerably useful for syntenic studies between supposedly remote species for the isolation of important genes for agronomy.

Site-Specific Insertion Polymorphism of the MITE Alex-1 in the Genus Coffea Suggests Interspecific Gene Flow.

Miniature Inverted-repeat Transposable Elements (MITEs) are small nonautonomous class-II transposable elements distributed throughout eukaryotic genomes. We identified a novel family of MITEs (named Alex) in the Coffea canephora genome often associated with expressed sequences. The Alex-1 element is inserted in an intron of a gene at the CcEIN4 locus. Its mobility was demonstrated by sequencing the insertion site in C. canephora accessions and Coffea species. Analysis of the insertion polymorphism of Alex-1 at this locus in Coffea species and in C. canephora showed that there was no relationship between the geographical distribution of the species, their phylogenetic relationships, and insertion polymorphism. The intraspecific distribution of C. canephora revealed an original situation within the E diversity group. These results suggest possibly greater gene flow between species than previously thought. This MITE family will enable the study of the C. canephora genome evolution, phylogenetic relationships, and possible gene flows within the Coffea genus.

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

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

Two novel Ty1-copia retrotransposons isolated from coffee trees can effectively reveal evolutionary relationships in the Coffea genus (Rubiaceae).

In the study, we developed new markers for phylogenetic relationships and intraspecies differentiation in Coffea. Nana and Divo, two novel Ty1-copia LTR-retrotransposon families, were isolated through C. canephora BAC clone sequencing. Nana- and Divo-based markers were used to test their: (1) ability to resolve recent phylogenetic relationships; (2) efficiency in detecting intra-species differentiation. Sequence-specific amplification polymorphism (SSAP), retrotransposon-microsatellite amplified polymorphism (REMAP) and retrotransposon-based insertion polymorphism (RBIP) approaches were applied to 182 accessions (31 Coffea species and one Psilanthus accession). Nana- and Divo-based markers revealed contrasted transpositional histories. At the BAC clone locus, RBIP results on C. canephora demonstrated that Nana insertion took place prior to C. canephora differentiation, while Divo insertion occurred after differentiation. Combined SSAP and REMAP data showed that Nana could resolve Coffea lineages, while Divo was efficient at a lower taxonomic level. The combined results indicated that the retrotransposon-based markers were useful in highlighting Coffea genetic diversity and the chronological pattern of speciation/differentiation events. Ongoing complete sequencing of the C. canephora genome will soon enable exhaustive identification of LTR-RTN families, as well as more precise in-depth analyses on contributions to genome size variation and Coffea evolution.