Species-level ichthyoplankton dynamics for 97 fishes in two major river basins of the Amazon using quantitative metabarcoding.

The Amazon basin holds the world's largest freshwater fish diversity. Information on the intensity and timing of reproductive ecology of Amazonian fish is scant. We use a metabarcoding method by capture using a single probe to quantify species-level ichthyoplankton dynamics. We sampled the Marañón and the Ucayali rivers in Peru monthly for 2 years. We identified 97 species that spawned mainly during the flood start, the flood end or the receding periods, although some species had spawning activity in more than one period. This information was new for 40 of the species in the Amazon basin and 80 species in Peru. Most species ceased spawning for a month during a strong hydrological anomaly in January 2016, demonstrating the rapidity with which they react to environmental modifications during the breeding season. We also document another unreported event in the Amazon basin, the inverse phenology of species belonging to one genus (Triportheus). Overall larval flow in the Marañón was more than twice that of the Ucayali, including for most commercial species (between two and 20 times higher), whereas the Ucayali accounts for ~80% of the fisheries landings in the region. Our results are discussed in the light of the main anthropogenic threats to fishes, hydropower dam construction and the Hidrovía Amazónica, and should serve as a pre-impact baseline.

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.

Transferability, development of simple sequence repeat (SSR) markers and application to the analysis of genetic diversity and population structure of the African fan palm (Borassus aethiopum Mart.) in Benin.

BACKGROUND: In Sub-Saharan Africa, Borassus aethiopum Mart. (African fan palm) is an important non-timber forest product-providing palm that faces multiple anthropogenic threats to its genetic diversity. However, this species is so far under-studied, which prevents its sustainable development as a resource. The present work is a first attempt at characterizing the genetic diversity and population structure of B. aethiopum across nine collection sites spanning the three climatic regions of Benin, West Africa, through the use of microsatellite markers. RESULTS: During a first phase we relied on the reported transferability of primers developed in other palm species. We find that, in disagreement with previously published results, only 22.5% of the markers tested enable amplification of B. aethiopum DNA and polymorphism detection is very low. In a second phase, we generated a B. aethiopum-specific genomic dataset through high-throughput sequencing and used it for the de novo detection of microsatellite loci. Among the primer pairs targeting these, 11 detected polymorphisms and were further used for analyzing genetic diversity. Across the nine sites, expected heterozygosity (He) ranges from 0.263 to 0.451 with an overall average of 0.354, showing a low genetic diversity. Analysis of molecular variance (AMOVA) shows that within-site variation accounts for 53% of the genetic variation. Accordingly, the low number of migrants and positive values of the fixation index (F) in sites from both the Central (Sudano-Guinean) and the Southern (Guinean) climatic regions suggest limited gene flow between sites. The global correlation between genetic and geographic distances is weak; however, our clustering analyses indicate that B. aethiopum palms from Savè (Center) are genetically more similar to those from the North than to samples from other Central sites. CONCLUSIONS: In the light of our results, we discuss the use of inter-species transfer vs. de novo development of microsatellite markers in genetic diversity analyses targeting under-studied species, and suggest future applications for our molecular resources. We propose that, while prominent short-range pollen and seed dispersal in Benin explain most of our results, gene flux between the Central and Northern regions, as a result of animal and/or human migrations, might underlie the Savè discrepancy.

Microsatellite markers development for Indonesian nutmeg (Myristica fragrans Houtt.) and transferability to other Myristicaceae spp.

Myristica fragrans (Myristicaceae) is a tropical evergreen tree that yields the two famous spices: nutmeg and mace. Despite its socio-economic importance, the spatial distribution of its genetic diversity is barely documented. In this aim, 48 nuclear microsatellite markers were isolated of which 14 were polymorphic in M. fragrans. Number of alleles per locus ranged from 2 to 6. The level of observed heterozygosity ranged from 0.038 to 0.929 across loci. Transferability of these microsatellites in other Myristica species (M. fatua, M. argentea, and M. crassipes) and Myristicaceae species (Horsfieldia palauensis) was tested and successful. These new microsatellites will be useful for future investigation on genetic diversity and population structure of M. fragrans and phylogenetically-related species.

A single amino acid substitution (H451Y) in Leishmania calcium-dependent kinase SCAMK confers high tolerance and resistance to antimony.

BACKGROUND: For almost a century, antimonials have remained the first-line drugs for the treatment of leishmaniasis. However, little is known about their mode of action and clinical resistance mechanisms. OBJECTIVES: We have previously shown that Leishmania nicotinamidase (PNC1) is an essential enzyme for parasite NAD+ homeostasis and virulence in vivo. Here, we found that parasites lacking the pnc1 gene (Δpnc1) are hypersusceptible to the active form of antimony (SbIII) and used these mutant parasites to better understand antimony's mode of action and the mechanisms leading to resistance. METHODS: SbIII-resistant WT and Δpnc1 parasites were selected in vitro by a stepwise selection method. NAD(H)/NADP(H) dosages and quantitative RT-PCR experiments were performed to explain the susceptibility differences observed between strains. WGS and a marker-free CRISPR/Cas9 base-editing approach were used to identify and validate the role of a new resistance mutation. RESULTS: NAD+-depleted Δpnc1 parasites were highly susceptible to SbIII and this phenotype could be rescued by NAD+ precursor or trypanothione precursor supplementation. Δpnc1 parasites could become resistant to SbIII by an unknown mechanism. WGS revealed a unique amino acid substitution (H451Y) in an EF-hand domain of an orphan calcium-dependent kinase, recently named SCAMK. When introduced into a WT reference strain by base editing, the H451Y mutation allowed Leishmania parasites to survive at extreme concentrations of SbIII, potentiating the rapid emergence of resistant parasites. CONCLUSIONS: These results establish that Leishmania SCAMK is a new central hub of antimony's mode of action and resistance development, and uncover the importance of drug tolerance mutations in the evolution of parasite drug resistance.

Identification of a Hypervirulent Pathotype of Rice yellow mottle virus: A Threat to Genetic Resistance Deployment in West-Central Africa.

Rice yellow mottle virus (RYMV) causes high losses to rice production in Africa. Several sources of varietal high resistance are available but the emergence of virulent pathotypes that are able to overcome one or two resistance alleles can sometimes occur. Both resistance spectra and viral adaptability have to be taken into account to develop sustainable rice breeding strategies against RYMV. In this study, we extended previous resistance spectrum analyses by testing the rymv1-4 and rymv1-5 alleles that are carried by the rice accessions Tog5438 and Tog5674, respectively, against isolates that are representative of RYMV genetic and pathogenic diversity. Our study revealed a hypervirulent pathotype, named thereafter pathotype T', that is able to overcome all known sources of high resistance. This pathotype, which is spatially localized in West-Central Africa, appears to be more abundant than previously suspected. To better understand the adaptive processes of pathotype T', molecular determinants of resistance breakdown were identified via Sanger sequencing and validated through directed mutagenesis of an infectious clone. These analyses confirmed the key role of convergent nonsynonymous substitutions in the central part of the viral genome-linked protein to overcome RYMV1-mediated resistance. In addition, deep-sequencing analyses revealed that resistance breakdown does not always coincide with fixed mutations. Actually, virulence mutations that are present in a small proportion of the virus population can be sufficient for resistance breakdown. Considering the spatial distribution of RYMV strains in Africa and their ability to overcome the RYMV resistance genes and alleles, we established a resistance-breaking risk map to optimize strategies for the deployment of sustainable and resistant rice lines in Africa.

Genome scan reveals selection acting on genes linked to stress response in wild pearl millet.

Uncovering genomic regions involved in adaption is a major goal in evolutionary biology. High-throughput sequencing now makes it possible to tackle this challenge in nonmodel species. Yet, despite the increasing number of methods targeted to specifically detect genomic footprints of selection, the complex demography of natural populations often causes high rates of false positive in gene discoveries. The aim of this study was to identify climate adaptations in wild pearl millet populations, Cenchrus americanus ssp. monodii. We focused on two climate gradients, one in Mali and one in Niger. We used a two-step strategy to limit false-positive outliers. First, we considered gradients as biological replicates and performed RNA sequencing of four populations at the extremities. We combined four methods-three based on differentiation among populations and one based on diversity patterns within populations-to identify outlier SNPs from a set of 87 218 high-quality SNPs. Among 11 155 contigs of pearl millet reference transcriptome, 540 exhibited selection signals as evidenced by at least one of the four methods. In a second step, we genotyped 762 samples in 11 additional populations distributed along the gradients using SNPs from the detected contigs and random SNPs as control. We further assessed selection on this large data set using a differentiation-based method and a method based on correlations with environmental variables based. Four contigs displayed consistent signatures between the four extreme and 11 additional populations, two of which were linked to abiotic and biotic stress responses.

Association mapping, patterns of linkage disequilibrium and selection in the vicinity of the PHYTOCHROME C gene in pearl millet.

KEY MESSAGE: Linkage analysis confirmed the association in the region of PHYC in pearl millet. The comparison of genes found in this region suggests that PHYC is the best candidate. Major efforts are currently underway to dissect the phenotype-genotype relationship in plants and animals using existing populations. This method exploits historical recombinations accumulated in these populations. However, linkage disequilibrium sometimes extends over a relatively long distance, particularly in genomic regions containing polymorphisms that have been targets for selection. In this case, many genes in the region could be statistically associated with the trait shaped by the selected polymorphism. Statistical analyses could help in identifying the best candidate genes into such a region where an association is found. In a previous study, we proposed that a fragment of the PHYTOCHROME C gene (PHYC) is associated with flowering time and morphological variations in pearl millet. In the present study, we first performed linkage analyses using three pearl millet F2 families to confirm the presence of a QTL in the vicinity of PHYC. We then analyzed a wider genomic region of ~100 kb around PHYC to pinpoint the gene that best explains the association with the trait in this region. A panel of 90 pearl millet inbred lines was used to assess the association. We used a Markov chain Monte Carlo approach to compare 75 markers distributed along this 100-kb region. We found the best candidate markers on the PHYC gene. Signatures of selection in this region were assessed in an independent data set and pointed to the same gene. These results foster confidence in the likely role of PHYC in phenotypic variation and encourage the development of functional studies.

Association studies including genotype by environment interactions: prospects and limits.

BACKGROUND: Association mapping studies offer great promise to identify polymorphisms associated with phenotypes and for understanding the genetic basis of quantitative trait variation. To date, almost all association mapping studies based on structured plant populations examined the main effects of genetic factors on the trait but did not deal with interactions between genetic factors and environment. In this paper, we propose a methodological prospect of mixed linear models to analyze genotype by environment interaction effects using association mapping designs. First, we simulated datasets to assess the power of linear mixed models to detect interaction effects. This simulation was based on two association panels composed of 90 inbreds (pearl millet) and 277 inbreds (maize). RESULTS: Based on the simulation approach, we reported the impact of effect size, environmental variation, allele frequency, trait heritability, and sample size on the power to detect the main effects of genetic loci and diverse effect of interactions implying these loci. Interaction effects specified in the model included SNP by environment interaction, ancestry by environment interaction, SNP by ancestry interaction and three way interactions. The method was finally used on real datasets from field experiments conducted on the two considered panels. We showed two types of interactions effects contributing to genotype by environment interactions in maize: SNP by environment interaction and ancestry by environment interaction. This last interaction suggests differential response at the population level in function of the environment. CONCLUSIONS: Our results suggested the suitability of mixed models for the detection of diverse interaction effects. The need of samples larger than that commonly used in current plant association studies is strongly emphasized to ensure rigorous model selection and powerful interaction assessment. The use of ancestry interaction component brought valuable information complementary to other available approaches.

Evolutionary history of pearl millet (Pennisetum glaucum [L.] R. Br.) and selection on flowering genes since its domestication.

The plant domestication process is associated with considerable modifications of plant phenotype. The identification of the genetic basis of this adaptation is of great interest for evolutionary biology. One of the methods used to identify such genes is the detection of signatures of selection. However, domestication is generally associated with major demographic effects. It is therefore crucial to disentangle the effects of demography and selection on diversity. In this study, we investigated selection in a flowering time pathway during domestication of pearl millet. We first used a random set of 20 genes to model pearl millet domestication using approximate Bayesian computation. This analysis showed that a model with exponential growth and wild-cultivated gene flow was well supported by our data set. Under this model, the domestication date of pearl millet is estimated at around 4,800 years ago. We assessed selection in 15 pearl millet DNA sequences homologous to flowering time genes and showed that these genes underwent selection more frequently than expected. We highlighted significant signatures of selection in six pearl millet flowering time genes associated with domestication or improvement of pearl millet. Moreover, higher deviations from neutrality were found for circadian clock-associated genes. Our study provides new insights into the domestication process of pearl millet and shows that a category of genes of the flowering pathway were preferentially selected during pearl millet domestication.

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.

Unveiling biogeographical patterns of the ichthyofauna in the Tuichi basin, a biodiversity hotspot in the Bolivian Amazon, using environmental DNA.

To date, more than 2400 valid fish species have been recorded in the Amazon basin. However, some regions remain poorly documented. This is the case in the Beni basin and in particular in one of its main sub-basins, the Tuichi, an Andean foothills rivers flowing through the Madidi National Park in the Bolivian Amazonia. The knowledge of its ichthyological diversity is, however, essential for the management and protection of aquatic ecosystems, which are threatened by the development of infrastructures (dams, factories and cities), mining and deforestation. Environmental DNA (eDNA) has been relatively little used so far in the Amazon basin. We sampled eDNA from water in 34 sites in lakes and rivers in the Beni basin including 22 sites in the Tuichi sub-basin, during the dry season. To assess the biogeographical patterns of the amazonian ichthyofauna, we implemented a metabarcoding approach using two pairs of specific primers designed and developed in our laboratory to amplify two partially overlapping CO1 fragments, one of 185bp and another of 285bp. We detected 252 fish taxa (207 at species level) among which 57 are newly identified for the Beni watershed. Species compositions are significantly different between lakes and rivers but also between rivers according to their hydrographic rank and altitude. Furthermore, the diversity patterns are related to the different hydro-ecoregions through which the Tuichi flows. The eDNA approach makes it possible to identify and complete the inventory of the ichthyofauna in this still poorly documented Amazon basin. However, taxonomic identification remains constrained by the lack of reference barcodes in public databases and does not allow the assignment of all OTUs. Our results can be taken into account in conservation and management strategies and could serve as a baseline for future studies, including on other Andean tributaries.

Bactris gasipaes Kunth var. gasipaes complete plastome and phylogenetic analysis.

Bactris gasipaes var. gasipaes (Arecaceae, Palmae) is an economically and socially important plant species for populations across tropical South and Central America. It has been domesticated from its wild variety, B. gasipaes var. chichagui, since pre-Columbian times. In this study, we sequenced the plastome of the cultivated variety, B. gasipaes Kunth var. gasipaes and compared it with the published plastome of the wild variety. The chloroplast sequence obtained was 156,580 bp. The cultivated chloroplast sequence was conserved compared to the wild type sequence with 99.8% of nucleotide identity. We did, however, identify multiple Single Nucleotide Variants (SNVs), insertions, microsatellites and a resolved region of missing nucleotides. A SNV in one of the core barcode markers (matK) was detected between the wild and cultivated accessions. Phylogenetic analysis was carried out across the Arecaceae family and compared to previous reports, resulting in an identical topology. This study is a step forward in understanding the genome evolution of this species.

Genomic footprints of selection in early-and late-flowering pearl millet landraces.

Pearl millet is among the top three-cereal production in one of the most climate vulnerable regions, sub-Saharan Africa. Its Sahelian origin makes it adapted to grow in poor sandy soils under low soil water regimes. Pearl millet is thus considered today as one of the most interesting crops to face the global warming. Flowering time, a trait highly correlated with latitude, is one of the key traits that could be modulated to face future global changes. West African pearl millet landraces, can be grouped into early- (EF) and late-flowering (LF) varieties, each flowering group playing a specific role in the functioning and resilience of Sahelian smallholders. The aim of this study was thus to detect genes linked to flowering but also linked to relevant traits within each flowering group. We thus investigated genomic and phenotypic diversity in 109 pearl millet landrace accessions, i.e., 66 early-flowering and 43 late-flowering, grown in the groundnut basin, the first area of rainfed agriculture in Senegal dominated by dry cereals (millet, maize, and sorghum) and legumes (groundnuts, cowpeas). We were able to confirm the role of PhyC gene in pearl millet flowering and identify several other genes that appear to be as much as important, such as FSR12 and HAC1. HAC1 and two other genes appear to be part of QTLs previously identified and deserve further investigation. At the same time, we were able to highlight a several genes and variants that could contribute to the improvement of pearl millet yield, especially since their impact was demonstrated across flowering cycles.

Genetic basis of pearl millet adaptation along an environmental gradient investigated by a combination of genome scan and association mapping.

Identifying the molecular bases of adaptation is a key issue in evolutionary biology. Genome scan is an efficient approach for identifying important molecular variation involved in adaptation. Association mapping also offers an opportunity to gain insight into genotype-phenotype relationships. Using these two approaches coupled with environmental data should help to come up with a refined picture of the evolutionary process underlying adaptation. In this study, we first conducted a selection scan analysis on a transcription factor gene family. We focused on the MADS-box gene family, a gene family which plays a crucial role in vegetative and flower development. Twenty-one pearl millet populations were sampled along an environmental gradient in West Africa. We identified one gene, i.e. PgMADS11, using Bayesian analysis to detect selection signatures. Polymorphism at this gene was also associated with flowering time variation in an association mapping framework. Finally, we found that PgMADS11 allele frequencies were closely associated with annual rainfall. Overall, we determined an efficient way to detect functional polymorphisms associated with climate variation in non-model plants by combining genome scan and association mapping. These results should help monitor the impact of recent climatic changes on plant adaptation.

Association studies identify natural variation at PHYC linked to flowering time and morphological variation in pearl millet.

The identification of genes selected during and after plant domestication is an important research topic to enhance knowledge on adaptative evolution. Adaptation to different climates was a key factor in the spread of domesticated crops. We conducted a study to identify genes responsible for these adaptations in pearl millet and developed an association framework to identify genetic variations associated with the phenotype in this species. A set of 90 inbred lines genotyped using microsatellite loci and AFLP markers was used. The population structure was assessed using two different Bayesian approaches that allow inbreeding or not. Association studies were performed using a linear mixed model considering both the population structure and familial relationships between inbred lines. We assessed the ability of the method to limit the number of false positive associations on the basis of the two different Bayesian methods, the number of populations considered and different morphological traits while also assessing the power of the methodology to detect given additive effects. Finally, we applied this methodology to a set of eight pearl millet genes homologous to cereal flowering pathway genes. We found significant associations between several polymorphisms of the pearl millet PHYC gene and flowering time, spike length, and stem diameter in the inbred line panel. To validate this association, we performed a second association analysis in a different set of pearl millet individuals from Niger. We confirmed a significant association between genetic variation in this gene and these characters.

Evidence of diploidy in the wild Amerindian yam, a putative progenitor of the endangered species Dioscorea trifida (Dioscoreaceae).

The fundamental question about Dioscorea trifida (Dioscoreaceae), the most important Amerindian yam, that remains unresolved concerns its evolutionary origin, since no wild relative has been reported. In this paper we report the existence of D. trifida's wild relative for the first time. The diploidy of wild D. trifida (2n = 40) is clearly demonstrated by flow cytometry, chromosome counts, and microsatellite pattern analysis, whereas the cultivated form was previously shown to be autotetraploid (2n = 80). In the coastal region where the wild and cultivated forms are sympatric, tetraploid and triploid cytotypes coexist within the same populations. In the sites where the wild and cultivated forms are allopatric, the wild diploid cytotype predominates. AFLP (amplified fragment length polymorphism) analyses gave an initial idea of the position of the wild forms in relation to the cultivated forms. All the wild and cultivated types form a monophyletic group structured into two major subgroups corresponding to the tetraploid cytotype of the cultivated form and the diploid cytotype of the wild form. The triploid cytotypes of the wild form are in an intermediary position. Wild accessions are grouped on the basis of their geographic origin. The data presented in this paper are significant for the effective breeding and conservation of D. trifida and to assess its genetic diversity and population structure for the general understanding of the evolution and domestication of the species.

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.

Abandonment of pearl millet cropping and homogenization of its diversity over a 40 year period in Senegal.

Cultivated diversity is considered an insurance against major climatic variability. However, since the 1980s, several studies have shown that climate variability and agricultural changes may already have locally eroded crop genetic diversity. We studied pearl millet diversity in Senegal through a comparison of pearl millet landraces collected 40 years apart. We found that more than 20% of villages visited in 1976 had stopped growing pearl millet. Despite this, its overall genetic diversity has been maintained but differentiation between early- and late-flowering accessions has been reduced. We also found stronger crop-to-wild gene flow than wild-to-crop gene flow and that wild-to-crop gene flow was weaker in 2016 than in 1976. In conclusion, our results highlight genetic homogenization in Senegal. This homogenization within cultivated pearl millet and between wild and cultivated forms is a key factor in genetic erosion and it is often overlooked. Improved assessment and conservation strategies are needed to promote and conserve both wild and cultivated pearl millet diversity.

Pearl millet genomic vulnerability to climate change in West Africa highlights the need for regional collaboration.

Climate change is already affecting agro-ecosystems and threatening food security by reducing crop productivity and increasing harvest uncertainty. Mobilizing crop diversity could be an efficient way to mitigate its impact. We test this hypothesis in pearl millet, a nutritious staple cereal cultivated in arid and low-fertility soils in sub-Saharan Africa. We analyze the genomic diversity of 173 landraces collected in West Africa together with an extensive climate dataset composed of metrics of agronomic importance. Mapping the pearl millet genomic vulnerability at the 2050 horizon based on the current genomic-climate relationships, we identify the northern edge of the current areas of cultivation of both early and late flowering varieties as being the most vulnerable to climate change. We predict that the most vulnerable areas will benefit from using landraces that already grow in equivalent climate conditions today. However, such seed-exchange scenarios will require long distance and trans-frontier assisted migrations. Leveraging genetic diversity as a climate mitigation strategy in West Africa will thus require regional collaboration.