A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

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.

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.

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.

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.

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

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.

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations.

Thanks to genome-scale diversity data, present-day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up-to-date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self-fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype-environment correlations and five designed to detect adaptive differentiation (based on F(ST) or similar measures). We show that genotype-environment correlation methods have substantially more power to detect selection than differentiation-based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype-environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.

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.

Biodiversity, evolution and adaptation of cultivated crops.

The human diet depends on very few crops. Current diversity in these crops is the result of a long interaction between farmers and cultivated plants, and their environment. Man largely shaped crop biodiversity from the domestication period 12,000 B.P. to the development of improved varieties during the last century. We illustrate this process through a detailed analysis of the domestication and early diffusion of maize. In smallholder agricultural systems, farmers still have a major impact on crop diversity today. We review several examples of the major impact of man on current diversity. Finally, biodiversity is considered to be an asset for adaptation to current environmental changes. We describe the evolution of pearl millet in West Africa, where average rainfall has decreased over the last forty years. Diversity in cultivated varieties has certainly helped this crop to adapt to climate variation.

Selection for earlier flowering crop associated with climatic variations in the Sahel.

Climate changes will have an impact on food production and will require costly adaptive responses. Adapting to a changing environment will be particularly challenging in sub-Saharan Africa where climate change is expected to have a major impact. However, one important phenomenon that is often overlooked and is poorly documented is the ability of agro-systems to rapidly adapt to environmental variations. Such an adaptation could proceed by the adoption of new varieties or by the adaptation of varieties to a changing environment. In this study, we analyzed these two processes in one of the driest agro-ecosystems in Africa, the Sahel. We performed a detailed study in Niger where pearl millet is the main crop and covers 65% of the cultivated area. To assess how the agro-system is responding to recent recurrent drought, we analyzed samples of pearl millet landraces collected in the same villages in 1976 and 2003 throughout the entire cultivated area of Niger. We studied phenological and morphological differences in the 1976 and 2003 collections by comparing them over three cropping seasons in a common garden experiment. We found no major changes in the main cultivated varieties or in their genetic diversity. However, we observed a significant shift in adaptive traits. Compared to the 1976 samples, samples collected in 2003 displayed a shorter lifecycle, and a reduction in plant and spike size. We also found that an early flowering allele at the PHYC locus increased in frequency between 1976 and 2003. The increase exceeded the effect of drift and sampling, suggesting a direct effect of selection for earliness on this gene. We conclude that recurrent drought can lead to selection for earlier flowering in a major Sahelian crop. Surprisingly, these results suggest that diffusion of crop varieties is not the main driver of short term adaptation to climatic variation.

Divergent expression patterns of miR164 and CUP-SHAPED COTYLEDON genes in palms and other monocots: implication for the evolution of meristem function in angiosperms.

In order to understand how the morphology of plant species has diversified over time, it is necessary to decipher how the underlying developmental programs have evolved. The regulatory network controlling shoot meristem activity is likely to have played an important role in morphological diversification and useful insights can be gained by comparing monocots and eudicots. These two distinct monophyletic groups of angiosperms diverged 130 Ma and are characterized by important differences in their morphology. Several studies of eudicot species have revealed a conserved role for NAM and CUC3 genes in meristem functioning and pattern formation through the definition of morphogenetic boundaries during development. In this study, we show that NAM- and CUC3-related genes are conserved in palms and grasses, their diversification having predated the radiation of monocots and eudicots. Moreover, the NAM-miR164 posttranscriptional regulatory module is also conserved in palm species. However, in contrast to the CUC3-related genes, which share a similar expression pattern between the two angiosperm groups, the expression domain of the NAM-miR164 module differs between monocot and eudicot species. In our studies of spatial expression patterns, we compared existing eudicot data with novel results from our work using two palm species (date palm and oil palm) and two members of the Poaceae (rice and millet). In addition to contrasting results obtained at the gene expression level, major differences were also observed between eudicot and monocot NAM-related genes in the occurrence of putative cis-regulatory elements in their promoter sequences. Overall, our results suggest that although NAM- and CUC3-related proteins are functionally equivalent between monocots and eudicots, evolutionary radiation has resulted in heterotopy through alterations in the expression domain of the NAM-miR164 regulatory module.

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.

The role of regulatory genes during maize domestication: evidence from nucleotide polymorphism and gene expression.

We investigated DNA sequence variation in 72 candidate genes in maize landraces and the wild ancestor of maize, teosinte. The candidate genes were chosen because they exhibit very low sequence diversity among maize inbreds and have sequence homology to known regulatory genes. We observed signatures of selection in 17 candidate genes, indicating that they were potential targets of artificial selection during domestication. In addition, 21 candidate genes were identified as potential targets of natural selection in teosinte. A comparison of the proportion of selected genes between our regulatory genes and genes unfiltered for their potential function (but also with very low sequence diversity among maize inbreds) provided some weak evidence that regulatory genes are overrepresented among selected genes. We detected no significant association between the positions of genes identified as potential targets of selection during domestication and quantitative trait loci (QTL) responsible for maize domestication traits. However, a subset of these genes, those identified by sequence homology as kinase/phosphatase genes, significantly cluster with the domestication QTL. We also analyzed expression profiles of genes in distinct maize tissues and observed that domestication genes are expressed on average at a significantly higher level than neutral genes in reproductive organs, including kernels.

A weak effect of background selection on trinucleotide microsatellites in maize.

Artificial selection during the domestication of maize is thought to have been predominantly positive and to have had little effect on the surrounding neutral diversity because linkage disequilibrium breaks down rapidly when physical distance increases. However, the degree to which indirect selection has shaped neutral diversity in the maize genome during domestication remains unclear. In this study, we investigate the relationship between local recombination rate and neutral polymorphism in maize and in teosinte using both sequence and microsatellite data. To quantify diversity, we estimate 3 parameters expected to differentially reflect the effects of indirect selection and mutation. We find no general correlation between diversity and recombination, indicating that indirect selection has had no genome-wide impact on maize diversity. However, we detect a weak correlation between heterozygosity and recombination for trinucleotide microsatellites deviating from the stepwise mutation model and located within genes (rho = 0.32, P < 0.03). This result can be explained by a background selection hypothesis. The fact that the same correlation is not confirmed for nucleotide diversity suggests that the strength of purifying selection at or near this class of microsatellites is higher than for nucleotide mutations.

Maize association population: a high-resolution platform for quantitative trait locus dissection.

Crop improvement and the dissection of complex genetic traits require germplasm diversity. Although this necessary phenotypic variability exists in diverse maize, most research is conducted using a small subset of inbred lines. An association population of 302 lines is now available--a valuable research tool that captures a large proportion of the alleles in cultivated maize. Provided that appropriate statistical models correcting for population structure are included, this tool can be used in association analyses to provide high-resolution evaluation of multiple alleles. This study describes the population structure of the 302 lines, and investigates the relationship between population structure and various measures of phenotypic and breeding value. On average, our estimates of population structure account for 9.3% of phenotypic variation, roughly equivalent to a major quantitative trait locus (QTL), with a high of 35%. Inclusion of population structure in association models is critical to meaningful analyses. This new association population has the potential to identify QTL with small effects, which will aid in dissecting complex traits and in planning future projects to exploit the rich diversity present in maize.