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.

A quantitative theory for genomic offset statistics

Genomic offset statistics predict the maladaptation of populations to rapid habitat alteration based on association of genotypes with environmental variation. Despite substantial evidence for empirical validity, genomic offset statistics have well-identified limitations, and lack a theory that would facilitate interpretations of predicted values. Here, we clarified the theoretical relationships between genomic offset statistics and unobserved fitness traits controlled by environmentally selected loci, and proposed a geometric measure to predict fitness after rapid change in local environment. The predictions of our theory were verified in computer simulations and in empirical data on African pearl millet (Cenchrus americanus) obtained from a common garden experiment. Our results proposed a unified perspective on genomic offset statistics, and provided a theoretical foundation necessary when considering their potential application in conservation management in the face of environmental change.

Crop domestication as a step towards reproductive isolation

Speciation, Darwin's mystery of mysteries, is a continuous process that results in genomic divergence accompanied by the gradual increment of reproductive barriers between lineages. Since the beginning of research on the genetics of speciation, several questions have emerged such as: What are the genetic bases of incompatibilities? How many loci are necessary to prevent hybridization and how are they distributed along genomes? Can speciation occur despite gene flow and how common is ecological speciation? Early stages of divergence are key to understand the ecology and genetics of speciation, and semi-isolated species where hybrids can still be produced are particularly relevant This article is protected by copyright. All rights reserved.

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 varietial 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.

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.

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.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

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.

Genetic control of rhizosheath formation in pearl millet

The rhizosheath, the layer of soil that adheres strongly to roots, influences water and nutrients acquisition. Pearl millet is a cereal crop that plays a major role for food security in arid regions of sub-Saharan Africa and India. We previously showed that root-adhering soil mass is a heritable trait in pearl millet and that it correlates with changes in rhizosphere microbiota structure and functions. Here, we studied the correlation between root-adhering soil mass and root hair development, root architecture, and symbiosis with arbuscular mycorrhizal fungi and we analysed the genetic control of this trait using genome wide association (GWAS) combined with bulk segregant analysis and gene expression studies. Root-adhering soil mass was weakly correlated only to root hairs traits in pearl millet. Twelve QTLs for rhizosheath formation were identified by GWAS. Bulk segregant analysis on a biparental population validated five of these QTLs. Combining genetics with a comparison of global gene expression in the root tip of contrasted inbred lines revealed candidate genes that might control rhizosheath formation in pearl millet. Our study indicates that rhizosheath formation is under complex genetic control in pearl millet and suggests that it is mainly regulated by root exudation.

Cultivated and wild pearl millet display contrasting patterns of abundance and co-occurrence in their root mycobiome

Fungal communities associated with roots play a key role in nutrient uptake and in mitigating the abiotic and biotic stress of their host. In this study, we characterized the roots mycobiome of wild and cultivated pearl millet [Pennisetum glaucum (L.) R. Br., synonym: Cenchrus americanus (L.) Morrone] in three agro-ecological areas of Senegal following a rainfall gradient. We hypothesized that wild pearl millet could serve as a reservoir of endophytes for cultivated pearl millet. We therefore analyzed the soil factors influencing fungal community structure and whether cultivated and wild millet shared the same fungal communities. The fungal communities associated with pearl millet were significantly structured according to sites and plant type (wild vs cultivated). Besides, soil pH and phosphorus were the main factors influencing the fungal community structure. We observed a higher fungal diversity in cultivated compared to wild pearl millet. Interestingly, we detected higher relative abundance of putative pathotrophs, especially plant pathogen, in cultivated than in wild millet in semi-arid and semi-humid zones, and higher relative abundance of saprotrophs in wild millet in arid and semi-humid zones. A network analysis based on taxa co-occurrence patterns in the core mycobiome revealed that cultivated millet and wild relatives had dissimilar groups of hub taxa. The identification of the core mycobiome and hub taxa of cultivated and wild pearl millet could be an important step in developing microbiome engineering approaches for more sustainable management practices in pearl millet agroecosystems.

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.

Trees and their seed networks: The social dynamics of urban fruit trees and implications for genetic diversity

Trees are a traditional component of urban spaces where they provide ecosystem services critical to urban wellbeing. In the Tropics, urban trees’ seed origins have rarely been characterized. Yet, understanding the social dynamics linked to tree planting is critical given their influence on the distribution of associated genetic diversity. This study examines elements of these dynamics (seed exchange networks) in an emblematic indigenous fruit tree species from Central Africa, the African plum tree (Dacryodes edulis, Burseraceae), within the urban context of Yaoundé. We further evaluate the consequences of these social dynamics on the distribution of the genetic diversity of the species in the city. Urban trees were planted predominantly using seeds sourced from outside the city, resulting in a level of genetic diversity as high in Yaoundé as in a whole region of production of the species. Debating the different drivers that foster the genetic diversity in planted urban trees, the study argues that cities and urban dwellers can unconsciously act as effective guardians of indigenous tree genetic diversity.

GWAS unveils features between early- and late-flowering pearl millets

Background

Pearl millet, a nutritious food for around 100 million people in Africa and India, displays extensive genetic diversity and a high degree of admixture with wild relatives. Two major morphotypes can be distinguished in Senegal: early-flowering Souna and late-flowering Sanio. Phenotypic variabilities related to flowering time play an important role in the adaptation of pearl millet to climate variability. A better understanding of the genetic makeup of these variabilities would make it possible to breed pearl millet to suit regions with different climates. The aim of this study was to characterize the genetic basis of these phenotypic differences.

Results

We defined a core collection that captures most of the diversity of cultivated pearl millets in Senegal and includes 60 early-flowering Souna and 31 late-flowering Sanio morphotypes. Sixteen agro-morphological traits were evaluated in the panel in the 2016 and 2017 rainy seasons. Phenological and phenotypic traits related with yield, flowering time, and biomass helped differentiate early- and late-flowering morphotypes. Further, using genotyping-by-sequencing (GBS), 21,663 single nucleotide polymorphisms (SNPs) markers with more than 5% of minor allele frequencies were discovered. Sparse non-negative matrix factorization (sNMF) analysis confirmed the genetic structure in two gene pools associated with differences in flowering time. Two chromosomal regions on linkage groups (LG 3) (~ 89.7 Mb) and (LG 6) (~ 68.1 Mb) differentiated two clusters among the early-flowering Souna. A genome-wide association study (GWAS) was used to link phenotypic variation to the SNPs, and 18 genes were linked to flowering time, plant height, tillering, and biomass (P-value < 2.3E-06).

Conclusions

The diversity of early- and late-flowering pearl millet morphotypes in Senegal was captured using a heuristic approach. Key phenological and phenotypic traits, SNPs, and candidate genes underlying flowering time, tillering, biomass yield and plant height of pearl millet were identified. Chromosome rearrangements in LG3 and LG6 were inferred as a source of variation in early-flowering morphotypes. Using candidate genes underlying these features between pearl millet morphotypes will be of paramount importance in breeding for resilience to climatic variability.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07198-2.

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.

Replacement of local crops with alternative varieties adapted to future conditions may improve food security under climate change. Here the authors apply landscape genomics and ensemble climate modelling to pearl millet in West Africa, supporting the potential of transfrontier assisted seed exchange.

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.

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.

Fonio millet genome unlocks African orphan crop diversity for agriculture in a changing climate

Sustainable food production in the context of climate change necessitates diversification of agriculture and a more efficient utilization of plant genetic resources. Fonio millet (Digitaria exilis) is an orphan African cereal crop with a great potential for dryland agriculture. Here, we establish high-quality genomic resources to facilitate fonio improvement through molecular breeding. These include a chromosome-scale reference assembly and deep re-sequencing of 183 cultivated and wild Digitaria accessions, enabling insights into genetic diversity, population structure, and domestication. Fonio diversity is shaped by climatic, geographic, and ethnolinguistic factors. Two genes associated with seed size and shattering showed signatures of selection. Most known domestication genes from other cereal models however have not experienced strong selection in fonio, providing direct targets to rapidly improve this crop for agriculture in hot and dry environments.

Fonio millet genome unlocks African orphan crop diversity for agriculture in a changing climate

Sustainable food production in the context of climate change necessitates diversification of agriculture and a more efficient utilization of plant genetic resources. Fonio millet (Digitaria exilis) is an orphan African cereal crop with a great potential for dryland agriculture. Here, we establish high-quality genomic resources to facilitate fonio improvement through molecular breeding. These include a chromosome-scale reference assembly and deep re-sequencing of 183 cultivated and wild Digitaria accessions, enabling insights into genetic diversity, population structure, and domestication. Fonio diversity is shaped by climatic, geographic, and ethnolinguistic factors. Two genes associated with seed size and shattering showed signatures of selection. Most known domestication genes from other cereal models however have not experienced strong selection in fonio, providing direct targets to rapidly improve this crop for agriculture in hot and dry environments.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

Yam genomics supports West Africa as a major cradle of crop domestication

While there has been progress in our understanding of the origin and history of agriculture in sub-Saharan Africa, a unified perspective is still lacking on where and how major crops were domesticated in the region. Here, we investigated the domestication of African yam (Dioscorea rotundata), a key crop in early African agriculture. Using whole-genome resequencing and statistical models, we show that cultivated yam was domesticated from a forest species. We infer that the expansion of African yam agriculture started in the Niger River basin. This result, alongside with the origins of African rice and pearl millet, supports the hypothesis that the vicinity of the Niger River was a major cradle of African agriculture.

Long-fragment targeted capture for long-read sequencing of plastomes

PREMISE

Third-generation sequencing methods generate significantly longer reads than those produced using alternative sequencing methods. This provides increased possibilities for the study of biodiversity, phylogeography, and population genetics. We developed a protocol for in-solution enrichment hybridization capture of long DNA fragments applicable to complete plastid genomes.

METHODS AND RESULTS

The protocol uses cost-effective in-house probes developed via long-range PCR and was used in six non-model monocot species (Poaceae: African rice, pearl millet, fonio; and three palm species). DNA was extracted from fresh and silica gel-dried leaves. Our protocol successfully captured long-read plastome fragments (3151 bp median on average), with an enrichment rate ranging from 15% to 98%. DNA extracted from silica gel-dried leaves led to low-quality plastome assemblies when compared to DNA extracted from fresh tissue.

CONCLUSIONS

Our protocol could also be generalized to capture long sequences from specific nuclear fragments.

Adaptive Introgression: An Untapped Evolutionary Mechanism for Crop Adaptation

Global environmental changes strongly impact wild and domesticated species biology and their associated ecosystem services. For crops, global warming has led to significant changes in terms of phenology and/or yield. To respond to the agricultural challenges of this century, there is a strong need for harnessing the genetic variability of crops and adapting them to new conditions. Gene flow, from either the same species or a different species, may be an immediate primary source to widen genetic diversity and adaptions to various environments. When the incorporation of a foreign variant leads to an increase of the fitness of the recipient pool, it is referred to as "adaptive introgression". Crop species are excellent case studies of this phenomenon since their genetic variability has been considerably reduced over space and time but most of them continue exchanging genetic material with their wild relatives. In this paper, we review studies of adaptive introgression, presenting methodological approaches and challenges to detecting it. We pay particular attention to the potential of this evolutionary mechanism for the adaptation of crops. Furthermore, we discuss the importance of farmers' knowledge and practices in shaping wild-to-crop gene flow. Finally, we argue that screening the wild introgression already existing in the cultivated gene pool may be an effective strategy for uncovering wild diversity relevant for crop adaptation to current environmental changes and for informing new breeding directions.

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.

Response to early drought stress and identification of QTLs controlling biomass production under drought in pearl millet

Pearl millet plays a major role in food security in arid and semi-arid areas of Africa and India. However, it lags behind the other cereal crops in terms of genetic improvement. The recent sequencing of its genome opens the way to the use of modern genomic tools for breeding. Our study aimed at identifying genetic components involved in early drought stress tolerance as a first step toward the development of improved pearl millet varieties or hybrids. A panel of 188 inbred lines from West Africa was phenotyped under early drought stress and well-irrigated conditions. We found a strong impact of drought stress on yield components. This impact was variable between inbred lines. We then performed an association analysis with a total of 392,493 SNPs identified using Genotyping-by-Sequencing (GBS). Correcting for genetic relatedness, genome wide association study identified QTLs for biomass production in early drought stress conditions and for stay-green trait. In particular, genes involved in the sirohaem and wax biosynthesis pathways were found to co-locate with two of these QTLs. Our results might contribute to breed pearl millet lines with improved yield under drought stress.

Metabarcoding by capture using a single COI probe (MCSP) to identify and quantify fish species in ichthyoplankton swarms

The ability to determine the composition and relative frequencies of fish species in large ichthyoplankton swarms could have extremely important ecological applications However, this task is currently hampered by methodological limitations. We proposed a new method for Amazonian species based on hybridization capture of the COI gene DNA from a distant species (Danio rerio), absent from our study area (the Amazon basin). The COI sequence of this species is approximately equidistant from all COI of Amazonian species available. By using this sequence as probe we successfully facilitated the simultaneous identification of fish larvae belonging to the order Siluriformes and to the Characiformes represented in our ichthyoplankton samples. Species relative frequencies, estimated by the number of reads, showed almost perfect correlations with true frequencies estimated by a Sanger approach, allowing the development of a quantitative approach. We also proposed a further improvement to a previous protocol, which enables lowering the sequencing effort by 40 times. This new Metabarcoding by Capture using a Single Probe (MCSP) methodology could have important implications for ecology, fisheries management and conservation in fish biodiversity hotspots worldwide. Our approach could easily be extended to other plant and animal taxa.

Fonio (Digitaria exilis (Kippist) Stapf): A Socially Embedded Cereal for Food and Nutrition Security in Senegal

Agricultural diversification with neglected and underutilized species is a viable way to sustainably increase the productivity of agrosystems. Understanding the social, cultural, and ecological roles of these species is crucial for their promotion. White fonio (Digitaria exilis), a neglected cereal endemic to West African Sahelian countries, is recognized as a crop for the future due to its cultural, nutritional, and economic values. In this study, we described fonio farming systems in Senegal through an ethnobotanical approach. As expected in family farming systems, farmers largely practiced diversified subsistence agriculture on small plots and relied on local seed exchange networks. The importance of fonio varied among agroecological zones, ethnic groups, and gender. In the Groundnut Basin, where agriculture is more mechanized, late-maturing landraces of fonio are cultivated as a cash crop rather than a staple crop. However, in southern Senegal, where food shortages are recurrent, fonio is more widespread and the cultivation of both early- and late-maturing landraces contributes to local food security. These differences also reflect the cultural status of fonio: different among ethnic groups and more important for women than for men. Finally, the regressive dynamics of fonio cultivation in most regions emphasizes the need to develop integrated conservation and promotion strategies that take into account the diversity of social and agronomical roles of the plant.

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.

No Excess of Cis-Regulatory Variation Associated with Intraspecific Selection in Wild Pearl Millet (Cenchrus americanus)

Several studies suggest that cis-regulatory mutations are the favorite target of evolutionary changes, one reason being that cis-regulatory mutations might have fewer deleterious pleiotropic effects than protein-coding mutations. A review of the process also suggests that this bias towards adaptive cis-regulatory variation might be less pronounced at the intraspecific level compared with the interspecific level. In this study, we assessed the contribution of cis-regulatory variation to adaptation at the intraspecific level using populations of wild pearl millet (Cenchrus americanus ssp. monodii) sampled along an environmental gradient in Niger. From RNA sequencing of hybrids to assess allele-specific expression, we identified genes with cis-regulatory divergence between two parental accessions collected in contrasted environmental conditions. This revealed that ∼15% of transcribed genes showed cis-regulatory variation. Intersecting the gene set exhibiting cis-regulatory variation with the gene set identified as targets of selection revealed no excess of cis-acting mutations among the selected genes. We additionally found no excess of cis-regulatory variation among genes associated with adaptive traits. As our approach relied on methods identifying mainly genes submitted to strong selection pressure or with high phenotypic effect, the contribution of cis-regulatory changes to soft selection or polygenic adaptive traits remains to be tested. However our results favor the hypothesis that enrichment of adaptive cis-regulatory divergence builds up over time. For short evolutionary time-scales, cis-acting mutations are not predominantly involved in adaptive evolution associated with strong selective signal.

Molecular basis of African yam domestication: analyses of selection point to root development, starch biosynthesis, and photosynthesis related genes

BACKGROUND

After cereals, root and tuber crops are the main source of starch in the human diet. Starch biosynthesis was certainly a significant target for selection during the domestication of these crops. But domestication of these root and tubers crops is also associated with gigantism of storage organs and changes of habitat.

RESULTS

We studied here, the molecular basis of domestication in African yam, Dioscorea rotundata. The genomic diversity in the cultivated species is roughly 30% less important than its wild relatives. Two percent of all the genes studied showed evidences of selection. Two genes associated with the earliest stages of starch biosynthesis and storage, the sucrose synthase 4 and the sucrose-phosphate synthase 1 showed evidence of selection. An adventitious root development gene, a SCARECROW-LIKE gene was also selected during yam domestication. Significant selection for genes associated with photosynthesis and phototropism were associated with wild to cultivated change of habitat. If the wild species grow as vines in the shade of their tree tutors, cultivated yam grows in full light in open fields.

CONCLUSIONS

Major rewiring of aerial development and adaptation for efficient photosynthesis in full light characterized yam domestication.

The Empirical Distribution of Singletons for Geographic Samples of DNA Sequences

Rare variants are important for drawing inference about past demographic events in a species history. A singleton is a rare variant for which genetic variation is carried by a unique chromosome in a sample. How singletons are distributed across geographic space provides a local measure of genetic diversity that can be measured at the individual level. Here, we define the empirical distribution of singletons in a sample of chromosomes as the proportion of the total number of singletons that each chromosome carries, and we present a theoretical background for studying this distribution. Next, we use computer simulations to evaluate the potential for the empirical distribution of singletons to provide a description of genetic diversity across geographic space. In a Bayesian framework, we show that the empirical distribution of singletons leads to accurate estimates of the geographic origin of range expansions. We apply the Bayesian approach to estimating the origin of the cultivated plant species Pennisetum glaucum [L.] R. Br. (pearl millet) in Africa, and find support for range expansion having started from Northern Mali. Overall, we report that the empirical distribution of singletons is a useful measure to analyze results of sequencing projects based on large scale sampling of individuals across geographic space.

Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments

Draft genome, 994 re-sequenced lines and GWAS for yield-traits provide a resource of genetics and genomics tools for pearl millet researchers and breeders.

Pearl millet [Cenchrus americanus (L.) Morrone] is a staple food for more than 90 million farmers in arid and semi-arid regions of sub-Saharan Africa, India and South Asia. We report the ∼1.79 Gb draft whole genome sequence of reference genotype Tift 23D₂B₁-P1-P5, which contains an estimated 38,579 genes. We highlight the substantial enrichment for wax biosynthesis genes, which may contribute to heat and drought tolerance in this crop. We resequenced and analyzed 994 pearl millet lines, enabling insights into population structure, genetic diversity and domestication. We use these resequencing data to establish marker trait associations for genomic selection, to define heterotic pools, and to predict hybrid performance. We believe that these resources should empower researchers and breeders to improve this important staple crop.

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.

Pearl Millet Genetic Traits Shape Rhizobacterial Diversity and Modulate Rhizosphere Aggregation

Root exudation contributes to soil carbon allocation and also to microbial C and energy supply, which subsequently impacts soil aggregation around roots. Biologically-driven soil structural formation is an important driver of soil fertility. Plant genetic determinants of exudation and more generally of factors promoting rhizosphere soil aggregation are largely unknown. Here, we characterized rhizosphere aggregation in a panel of 86 pearl millet inbred lines using a ratio of root-adhering soil dry mass per root tissue dry mass (RAS/RT). This ratio showed significant variations between lines, with a roughly 2-fold amplitude between lowest and highest average values. For 9 lines with contrasting aggregation properties, we then compared the bacterial diversity and composition in root-adhering soil. Bacterial α-diversity metrics increased with the "RAS/RT ratio." Regarding taxonomic composition, the Rhizobiales were stimulated in lines showing high aggregation level whereas Bacillales were more abundant in lines with low ratio. 184 strains of cultivable exopolysaccharides-producing bacteria have been isolated from the rhizosphere of some lines, including members from Rhizobiales and Bacillales. However, at this stage, we could not find a correlation between abundance of EPS-producing species in bacterial communities and the ratio RAS/RT. These results illustrated the impact of cereals genetic trait variation on soil physical properties and microbial diversity. This opens the possibility of considering plant breeding to help management of soil carbon content and physical characteristics through carbon rhizodeposition in soil.

Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin

Local people's perceptions of cultivated and wild agrobiodiversity, as well as their management of hybridization are still understudied in Amazonia. Here we analyze domesticated treegourd (Crescentia cujete), whose versatile fruits have technological, symbolic, and medicinal uses. A wild relative (C. amazonica) of the cultivated species grows spontaneously in Amazonian flooded forests. We demonstrated, using whole chloroplast sequences and nuclear microsatellites, that the two species are strongly differentiated. Nonetheless, they hybridize readily throughout Amazonia and the proportions of admixture correlate with fruit size variation of cultivated trees. New morphotypes arise from hybridization, which are recognized by people and named as local varieties. Small hybrid fruits are used to make the important symbolic rattle (maracá), suggesting that management of hybrid trees is an ancient human practice in Amazonia. Effective conservation of Amazonian agrobiodiversity needs to incorporate this interaction between wild and cultivated populations that is managed by smallholder families. Beyond treegourd, our study clearly shows that hybridization plays an important role in tree crop phenotypic diversification and that the integration of molecular analyses and farmers’ perceptions of diversity help disentangle crop domestication history.

New Genetic Insights into Pearl Millet Diversity As Revealed by Characterization of Early- and Late-Flowering Landraces from Senegal

Pearl millet (Pennisetum glaucum (L.) R. Br.) is a staple food and a drought-tolerant cereal well adapted to Sub-Saharan Africa agro-ecosystems. An important diversity of pearl millet landraces has been widely conserved by farmers and therefore could help copping with climate changes and contribute to future food security. Hence, characterizing its genetic diversity and population structure can contribute to better assist breeding programs for a sustainable agricultural productivity enhancement. Toward this goal, a comprehensive panel of 404 accessions were used that correspond to 12 improved varieties, 306 early flowering and 86 late-flowering cultivated landraces from Senegal. Twelve highly polymorphic SSR markers were used to study diversity and population structure. Two genes, PgMADS11 and PgPHYC, were genotyped to assess their association to flowering phenotypic difference in landraces. Results indicate a large diversity and untapped potential of Senegalese pearl millet germplasm as well as a genetic differentiation between early- and late-flowering landraces. Further, a fine-scale genetic difference of PgPHYC and PgMADS11 (SNP and indel, respectively) and co-variation of their alleles with flowering time were found among landraces. These findings highlight new genetic insights of pearl millet useful to define heterotic populations for breeding, genomic association panel, or crosses for trait-specific mapping.

Diversity of Algerian oases date palm (Phoenix dactylifera L., Arecaceae): Heterozygote excess and cryptic structure suggest farmer management had a major impact on diversity

Date palm (Phoenix dactyliferaL.) is the mainstay of oasis agriculture in the Saharan region. It is cultivated in a large part of the Mediterranean coastal area of the Sahara and in most isolated oases in the Algerian desert. We sampled 10 oases in Algeria to understand the structure of date palm diversity from the coastal area to a very isolated desert location. We used 18 microsatellite markers and a chloroplast minisatellite to characterize 414 individual palm trees corresponding to 114 named varieties. We found a significant negative inbreeding coefficient, suggesting active farmer selection for heterozygous individuals. Three distinct genetic clusters were identified, a ubiquitous set of varieties found across the different oases, and two clusters, one of which was specific to the northern area, and the other to the drier southern area of the Algerian Sahara. The ubiquitous cluster presented very striking chloroplast diversity, signing the frequency of haplotypes found in Saudi Arabia, the most eastern part of the date palm range. Exchanges of Middle Eastern and Algerian date palms are known to have occurred and could have led to the introduction of this particular chlorotype. However, Algerian nuclear diversity was not of eastern origin. Our study strongly suggests that the peculiar chloroplastic diversity of date palm is maintained by farmers and could originate from date palms introduced from the Middle East a long time ago, which since then, hasbeen strongly introgressed. This study illustrates the complex structure of date palm diversity in Algerian oases and the role of farmers in shaping such cryptic diversity.

Myosin XI is associated with fitness and adaptation to aridity in wild pearl millet

Phenotypic changes in plants can be observed along many environmental gradients and are determined by both environmental and genetic factors. The identification of alleles associated with phenotypic variations is a rapidly developing area of research. We studied the genetic basis of phenotypic variations in 11 populations of wild pearl millet (Pennisetum glaucum) on two North-South aridity gradients, one in Niger and one in Mali. Most of the 11 phenotypic traits assessed in a common garden experiment varied between the populations studied. Moreover, the size of the inflorescence, the number of flowers and aboveground dry mass co-varied positively with a decrease in rainfall. To decipher the genetic basis of these phenotypes, we used an association mapping strategy with a mixed model. We found two SNPs on the same myosin XI contig significantly associated with variations in the average number of flowers. Both the allele frequency of the two SNPs and the average number of flowers co-varied with the rainfall gradient on the two gradients. Interestingly, this gene was also a target of selection during domestication. The Myosin XI gene is thus a good candidate for fitness-related adaptation in wild populations.

miR2118-triggered phased siRNAs are differentially expressed during the panicle development of wild and domesticated African rice species

BACKGROUND

Rice exhibits a wide range of panicle structures. To explain these variations, much emphasis has been placed on changes in transcriptional regulation, but no large-scale study has yet reported on changes in small RNA regulation in the various rice species. To evaluate this aspect, we performed deep sequencing and expression profiling of small RNAs from two closely related species with contrasting panicle development: the cultivated African rice Oryza glaberrima and its wild relative Oryza barthii.

RESULTS

Our RNA-seq analysis revealed a dramatic difference between the two species in the 21 nucleotide small RNA population, corresponding mainly to miR2118-triggered phased siRNAs. A detailed expression profiling during the panicle development of O. glaberrima and O. barthii using qRT-PCRs and in situ hybridization, confirmed a delayed expression of the phased siRNAs as well as their lncRNA precursors and regulators (miR2118 and MEL1 gene) in O. glaberrima compared to O. barthii. We provide evidence that the 21-nt phasiRNA pathway in rice is associated with male-gametogenesis but is initiated in spikelet meristems.

CONCLUSION

Differential expression of the miR2118-triggered 21-nt phasiRNA pathway between the two African rice species reflects differential rates of determinate fate acquisition of panicle meristems between the two species.

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.

Chloroplast sequence of treegourd (Crescentia cujete, Bignoniaceae) to study phylogeography and domestication

PREMISE OF THE STUDY

Crescentia cujete (Bignoniaceae) fruit rinds are traditionally used for storage vessels and handicrafts. We assembled its chloroplast genome and identified single-nucleotide polymorphisms (SNPs).

METHODS AND RESULTS

Using a genome skimming approach, the whole chloroplast of C. cujete was assembled using 3,106,928 sequence reads of 150 bp. The chloroplast is 154,662 bp in length, structurally divided into a large single copy region (84,788 bp), a small single copy region (18,299 bp), and two inverted repeat regions (51,575 bp) with 88 genes annotated. By resequencing the whole chloroplast, we identified 66 SNPs in C. cujete (N = 30) and 68 SNPs in C. amazonica (N = 6). Nucleotide diversity was estimated at 1.1 × 10^-3 and 3.5 × 10^-3 for C. cujete and C. amazonica, respectively.

CONCLUSIONS

This broadened C. cujete genetic toolkit will be important to study the origin, domestication, diversity, and phylogeography of treegourds in the Neotropics.

A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives

We produced a unique large data set of reference transcriptomes to obtain new knowledge about the evolution of plant genomes and crop domestication. For this purpose, we validated a RNA-Seq data assembly protocol to perform comparative population genomics. For the validation, we assessed and compared the quality of de novo Illumina short-read assemblies using data from two crops for which an annotated reference genome was available, namely grapevine and sorghum. We used the same protocol for the release of 26 new transcriptomes of crop plants and wild relatives, including still understudied crops such as yam, pearl millet and fonio. The species list has a wide taxonomic representation with the inclusion of 15 monocots and 11 eudicots. All contigs were annotated using BLAST, prot4EST and Blast2GO. A strong originality of the data set is that each crop is associated with close relative species, which will permit whole-genome comparative evolutionary studies between crops and their wild-related species. This large resource will thus serve research communities working on both crops and model organisms. All the data are available at http://arcad-bioinformatics.southgreen.fr/.

Spatial and Temporal Variation in Selection of Genes Associated with Pearl Millet Varietal Quantitative Traits In situ

Ongoing global climate changes imply new challenges for agriculture. Whether plants and crops can adapt to such rapid changes is still a widely debated question. We previously showed adaptation in the form of earlier flowering in pearl millet at the scale of a whole country over three decades. However, this analysis did not deal with variability of year to year selection. To understand and possibly manage plant and crop adaptation, we need more knowledge of how selection acts in situ. Is selection gradual, abrupt, and does it vary in space and over time? In the present study, we tracked the evolution of allele frequency in two genes associated with pearl millet phenotypic variation in situ. We sampled 17 populations of cultivated pearl millet over a period of 2 years. We tracked changes in allele frequencies in these populations by genotyping more than seven thousand individuals. We demonstrate that several allele frequencies changes are compatible with selection, by correcting allele frequency changes associated with genetic drift. We found marked variation in allele frequencies from year to year, suggesting a variable selection effect in space and over time. We estimated the strength of selection associated with variations in allele frequency. Our results suggest that the polymorphism maintained at the genes we studied is partially explained by the spatial and temporal variability of selection. In response to environmental changes, traditional pearl millet varieties could rapidly adapt thanks to this available functional variability.