The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

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.

Glutaredoxin regulation of primary root growth is associated with early drought stress tolerance in pearl millet.

Seedling root traits impact plant establishment under challenging environments. Pearl millet is one of the most heat and drought tolerant cereal crops that provides a vital food source across the sub-Saharan Sahel region. Pearl millet's early root system features a single fast-growing primary root which we hypothesize is an adaptation to the Sahelian climate. Using crop modeling, we demonstrate that early drought stress is an important constraint in agrosystems in the Sahel where pearl millet was domesticated. Furthermore, we show that increased pearl millet primary root growth is correlated with increased early water stress tolerance in field conditions. Genetics including genome-wide association study and quantitative trait loci (QTL) approaches identify genomic regions controlling this key root trait. Combining gene expression data, re-sequencing and re-annotation of one of these genomic regions identified a glutaredoxin-encoding gene PgGRXC9 as the candidate stress resilience root growth regulator. Functional characterization of its closest Arabidopsis homolog AtROXY19 revealed a novel role for this glutaredoxin (GRX) gene clade in regulating cell elongation. In summary, our study suggests a conserved function for GRX genes in conferring root cell elongation and enhancing resilience of pearl millet to its Sahelian environment.

Pearl millet is a staple food for over 90 million people living in regions of Africa and India that typically experience high temperatures and little rainfall. It was domesticated about 4,500 years ago in the Sahel region of West Africa and is one of the most heat and drought tolerant cereal crops worldwide. In most plants, organs known as roots absorb water and essential nutrients from the soil. Young pearl millet plants develop a fast-growing primary root, but it is unclear how this unique feature helps the crop to grow in hot and dry conditions. Using weather data collected from the Sahel over a 20-year period, Fuente, Grondin et al. predicted by modelling that early drought stress is the major factor limiting pearl millet growth and yield in this region. Field experiments found that plants with primary roots that grow faster within soil were better at tolerating early drought than those with slower growing roots. Further work using genetic approaches revealed that a gene known as PgGRXC9 promotes the growth of the primary root. To better understand how this gene works, the team examined a very similar gene in a well-studied model plant known as Arabidopsis. This suggested that PgGRXC9 helps the primary root to grow by stimulating cell elongation within the root. Since it is well adapted to dry conditions, pearl millet is expected to play an important role in helping agriculture adjust to climate change. The findings of Fuente, Grondin et al. may be used by plant breeders to create more resilient and productive varieties of pearl millet.

Genome-wide association analysis identifies natural allelic variants associated with panicle architecture variation in African rice, Oryza glaberrima Steud.

African rice (Oryza glaberrima Steud), a short-day cereal crop closely related to Asian rice (Oryza sativa L.), has been cultivated in Sub-Saharan Africa for ∼ 3,000 years. Although less cultivated globally, it is a valuable genetic resource in creating high-yielding cultivars that are better adapted to diverse biotic and abiotic stresses. While inflorescence architecture, a key trait for rice grain yield improvement, has been extensively studied in Asian rice, the morphological and genetic determinants of this complex trait are less understood in African rice. In this study, using a previously developed association panel of 162 O. glaberrima accessions and new SNP variants characterized through mapping to a new version of the O. glaberrima reference genome, we conducted a genome-wide association study of four major morphological panicle traits. We have found a total of 41 stable genomic regions that are significantly associated with these traits, of which 13 co-localized with previously identified QTLs in O. sativa populations and 28 were unique for this association panel. Additionally, we found a genomic region of interest on chromosome 3 that was associated with the number of spikelets and primary and secondary branches. Within this region was localized the O. sativa ortholog of the PHYTOCHROME B gene (Oglab_006903/OgPHYB). Haplotype analysis revealed the occurrence of natural sequence variants at the OgPHYB locus associated with panicle architecture variation through modulation of the flowering time phenotype, whereas no equivalent alleles were found in O. sativa. The identification in this study of genomic regions specific to O. glaberrima indicates panicle-related intra-specific genetic variation in this species, increasing our understanding of the underlying molecular processes governing panicle architecture. Identified candidate genes and major haplotypes may facilitate the breeding of new African rice cultivars with preferred panicle traits.

Chlorophyll fluorescence-based high-throughput phenotyping facilitates the genetic dissection of photosynthetic heat tolerance in African (Oryza glaberrima) and Asian (Oryza sativa) rice.

Rising temperatures and extreme heat events threaten rice production. Half of the global population relies on rice for basic nutrition, and therefore developing heat-tolerant rice is essential. During vegetative development, reduced photosynthetic rates can limit growth and the capacity to store soluble carbohydrates. The photosystem II (PSII) complex is a particularly heat-labile component of photosynthesis. We have developed a high-throughput chlorophyll fluorescence-based screen for photosynthetic heat tolerance capable of screening hundreds of plants daily. Through measuring the response of maximum PSII efficiency to increasing temperature, this platform generates data for modelling the PSII-temperature relationship in large populations in a small amount of time. Coefficients from these models (photosynthetic heat tolerance traits) demonstrated high heritabilities across African (Oryza glaberrima) and Asian (Oryza sativa, Bengal Assam Aus Panel) rice diversity sets, highlighting valuable genetic variation accessible for breeding. Genome-wide association studies were performed across both species for these traits, representing the first documented attempt to characterize the genetic basis of photosynthetic heat tolerance in any species to date. A total of 133 candidate genes were highlighted. These were significantly enriched with genes whose predicted roles suggested influence on PSII activity and the response to stress. We discuss the most promising candidates for improving photosynthetic heat tolerance in rice.

An improved assembly of the pearl millet reference genome using Oxford Nanopore long reads and optical mapping.

Pearl millet (Pennisetum glaucum (L.)) R. Br. syn. Cenchrus americanus (L.) Morrone) is an important crop in South Asia and sub-Saharan Africa which contributes to ensuring food security. Its genome has an estimated size of 1.76 Gb and displays a high level of repetitiveness above 80%. A first assembly was previously obtained for the Tift 23D2B1-P1-P5 cultivar genotype using short-read sequencing technologies. This assembly is, however, incomplete and fragmented with around 200 Mb unplaced on chromosomes. We report here an improved quality assembly of the pearl millet Tift 23D2B1-P1-P5 cultivar genotype obtained with an approach combining Oxford Nanopore long reads and Bionano Genomics optical maps. This strategy allowed us to add around 200 Mb at the chromosome-level assembly. Moreover, we strongly improved continuity in the order of the contigs and scaffolds within the chromosomes, particularly in the centromeric regions. Notably, we added more than 100 Mb around the centromeric region on chromosome 7. This new assembly also displayed a higher gene completeness with a complete BUSCO score of 98.4% using the Poales database. This more complete and higher quality assembly of the Tift 23D2B1-P1-P5 genotype now available to the community will help in the development of research on the role of structural variants and more broadly in genomics studies and the breeding of pearl millet.

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.

Redox-related gene expression and sugar accumulation patterns are altered in the edible inflorescence produced by the cultivated form of pacaya palm (Chamaedorea tepejilote).

BACKGROUND AND AIMS: The pacaya palm is a dioecious neotropical palm species that is exploited in Latin America for its male inflorescence, which is edible when immature. It is cultivated, in a non-intensive manner, in Guatemala, where a morphotype occurs that produces much larger, more highly branched inflorescences compared with wild palms. We sought to identify molecular factors underlying this phenotypic divergence, which is likely to be a product of domestication. METHODS: We performed RNA-seq-based studies on immature pacaya palm male inflorescences in order to identify genes that might be directly or indirectly affected in their expression in relation to domestication. We also measured the accumulation of a range of soluble sugar molecules to provide information on the biochemical status of the two different types of material. KEY RESULTS: A total of 408 genes were found to display significantly different expression levels between the wild and cultivated morphotypes. Three different functional categories were found to be enriched in the gene set that was upregulated in the cultivated morphotype: redox balance; secondary metabolism; and transport. Several sugars were found to accumulate at higher levels in inflorescences of the cultivated morphotype, in particular myo-inositol, fructose and glucose. CONCLUSIONS: The observed upregulation of redox-related genes in the cultivated morphotype is corroborated by the observation of higher myo-inositol accumulation, which has been shown to be associated with enhanced scavenging of reactive oxygen species in other plants and which may affect meristem activity.

Aquaporins are main contributors to root hydraulic conductivity in pearl millet [Pennisetum glaucum (L) R. Br.].

Pearl millet is a key cereal for food security in arid and semi-arid regions but its yield is increasingly threatened by water stress. Physiological mechanisms relating to conservation of soil water or increased water use efficiency can alleviate that stress. Aquaporins (AQP) are water channels that mediate root water transport, thereby influencing plant hydraulics, transpiration and soil water conservation. However, AQP remain largely uncharacterized in pearl millet. Here, we studied AQP function in root water transport in two pearl millet lines contrasting for water use efficiency (WUE). We observed that these lines also contrasted for root hydraulic conductivity (Lpr) and AQP contribution to Lpr. The line with lower WUE showed significantly higher AQP contribution to Lpr. To investigate AQP isoforms contributing to Lpr, we developed genomic approaches to first identify the entire AQP family in pearl millet and secondly, characterize the plasma membrane intrinsic proteins (PIP) gene expression profile. We identified and annotated 33 AQP genes in pearl millet, among which ten encoded PIP isoforms. PgPIP1-3 and PgPIP1-4 were significantly more expressed in the line showing lower WUE, higher Lpr and higher AQP contribution to Lpr. Overall, our study suggests that the PIP1 AQP family are the main regulators of Lpr in pearl millet and may possibly be associated with mechanisms associated to whole plant water use. This study paves the way for further investigations on AQP functions in pearl millet hydraulics and adaptation to environmental stresses.

Genome Wide Association Study Pinpoints Key Agronomic QTLs in African Rice Oryza glaberrima.

BACKGROUND: African rice, Oryza glaberrima, is an invaluable resource for rice cultivation and for the improvement of biotic and abiotic resistance properties. Since its domestication in the inner Niger delta ca. 2500 years BP, African rice has colonized a variety of ecologically and climatically diverse regions. However, little is known about the genetic basis of quantitative traits and adaptive variation of agricultural interest for this species. RESULTS: Using a reference set of 163 fully re-sequenced accessions, we report the results of a Genome Wide Association Study carried out for African rice. We investigated a diverse panel of traits, including flowering date, panicle architecture and resistance to Rice yellow mottle virus. For this, we devised a pipeline using complementary statistical association methods. First, using flowering time as a target trait, we found several association peaks, one of which co-localised with a well described gene in the Asian rice flowering pathway, OsGi, and identified new genomic regions that would deserve more study. Then we applied our pipeline to panicle- and resistance-related traits, highlighting some interesting genomic regions and candidate genes. Lastly, using a high-resolution climate database, we performed an association analysis based on climatic variables, searching for genomic regions that might be involved in adaptation to climatic variations. CONCLUSION: Our results collectively provide insights into the extent to which adaptive variation is governed by sequence diversity within the O. glaberrima genome, paving the way for in-depth studies of the genetic basis of traits of interest that might be useful to the rice breeding community.

MooSciTIC: Training of trainers in West African research and higher education.

The MooSciTIC project is a capacity-building initiative targeting West African research scientists and higher education teachers. The project aimed to improve the self-reliance of researchers and upgrade research practices by providing on-site summer schools on trans-disciplinary topics such as scientific writing, communication, and integrity. Here, we explain how this program was designed and implemented and share the positive responses from our trainees, hoping to inspire similar initiatives.

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.

The Rise and Fall of African Rice Cultivation Revealed by Analysis of 246 New Genomes.

African rice (Oryza glaberrima) was domesticated independently from Asian rice. The geographical origin of its domestication remains elusive. Using 246 new whole-genome sequences, we inferred the cradle of its domestication to be in the Inner Niger Delta. Domestication was preceded by a sharp decline of most wild populations that started more than 10,000 years ago. The wild population collapse occurred during the drying of the Sahara. This finding supports the hypothesis that depletion of wild resources in the Sahara triggered African rice domestication. African rice cultivation strongly expanded 2,000 years ago. During the last 5 centuries, a sharp decline of its cultivation coincided with the introduction of Asian rice in Africa. A gene, PROG1, associated with an erect plant architecture phenotype, showed convergent selection in two rice cultivated species, Oryza glaberrima from Africa and Oryza sativa from Asia. In contrast, a shattering gene, SH5, showed selection signature during African rice domestication, but not during Asian rice domestication. Overall, our genomic data revealed a complex history of African rice domestication influenced by important climatic changes in the Saharan area, by the expansion of African agricultural society, and by recent replacement by another domesticated species.

De Novo Assemblies of Three Oryza glaberrima Accessions Provide First Insights about Pan-Genome of African Rices.

Oryza glaberrima is one of the two cultivated species of rice, and harbors various interesting agronomic traits, especially in biotic and abiotic resistance, compared with its Asian cousin O. sativa. A previous reference genome was published but newer studies highlighted some missing parts. Moreover, global species diversity is known nowadays to be represented by more than one single individual. For that purpose, we sequenced, assembled and annotated de novo three different cultivars from O. glaberrima. After validating our assemblies, we were able to better solve complex regions than the previous assembly and to provide a first insight in pan-genomic divergence between individuals. The three assemblies shown large common regions, but almost 25% of the genome present collinearity breakpoints or are even individual specific.

LTRclassifier: A website for fast structural LTR retrotransposons classification in plants.

Automatic classification of LTR retrotransposons is a big challenge in the area of massive genomics. Many tools were developed to detect them but automatic classification is somehow challenging. Here we propose a simple approach, LTRclassifier, based on HMM recognition followed by BLAST analyses (i) to classify plant LTR retrotransposons in their respective superfamily, and (ii) to provide automatically a basic functional annotation of these elements. The method was tested on various TE databases, and shown to be robust and fast. This tool is available as a web service implemented at IRD bioinformatics facility, http://LTRclassifier.ird.fr/.

TOGGLE: toolbox for generic NGS analyses.

BACKGROUND: The explosion of NGS (Next Generation Sequencing) sequence data requires a huge effort in Bioinformatics methods and analyses. The creation of dedicated, robust and reliable pipelines able to handle dozens of samples from raw FASTQ data to relevant biological data is a time-consuming task in all projects relying on NGS. To address this, we created a generic and modular toolbox for developing such pipelines. RESULTS: TOGGLE (TOolbox for Generic nGs anaLysEs) is a suite of tools able to design pipelines that manage large sets of NGS softwares and utilities. Moreover, TOGGLE offers an easy way to manipulate the various options of the different softwares through the pipelines in using a single basic configuration file, which can be changed for each assay without having to change the code itself. We also describe one implementation of TOGGLE in a complete analysis pipeline designed for SNP discovery for large sets of genomic data, ready to use in different environments (from a single machine to HPC clusters). CONCLUSION: TOGGLE speeds up the creation of robust pipelines with reliable log tracking and data flow, for a large range of analyses. Moreover, it enables Biologists to concentrate on the biological relevance of results, and change the experimental conditions easily. The whole code and test data are available at https://github.com/SouthGreenPlatform/TOGGLE .

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.

In silico mining of microsatellites in coding sequences of the date palm (Arecaceae) genome, characterization, and transferability.

PREMISE OF THE STUDY: To complement existing sets of primarily dinucleotide microsatellite loci from noncoding sequences of date palm, we developed primers for tri- and hexanucleotide microsatellite loci identified within genes. Due to their conserved genomic locations, the primers should be useful in other palm taxa, and their utility was tested in seven other Phoenix species and in Chamaerops, Livistona, and Hyphaene. • METHODS AND RESULTS: Tandem repeat motifs of 3-6 bp were searched using a simple sequence repeat (SSR)-pipeline package in coding portions of the date palm draft genome sequence. Fifteen loci produced highly consistent amplification, intraspecific polymorphisms, and stepwise mutation patterns. • CONCLUSIONS: These microsatellite loci showed sufficient levels of variability and transferability to make them useful for population genetic, selection signature, and interspecific gene flow studies in Phoenix and other Coryphoideae genera.

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.