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.

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.

Common microgeographical selection patterns revealed in four European conifers.

Microgeographical adaptation occurs when the effects of directional selection persist despite gene flow. Traits and genetic loci under selection can then show adaptive divergence, against the backdrop of little differentiation at other traits or loci. How common such events are and how strong the selection is that underlies them remain open questions. Here, we discovered and analysed microgeographical patterns of genomic divergence in four European and Mediterranean conifers with widely differing life-history traits and ecological requirements (Abies alba MIll., Cedrus atlantica [Endl.] Manetti, Pinus halepensis Mill. and Pinus pinaster Aiton) by screening pairs from geographically close forest stands sampled along steep ecological gradients. We inferred patterns of genomic divergence by applying a combination of divergence outlier detection methods, demographic modelling, Approximate Bayesian Computation inferences and genomic annotation to genomic data. Surprisingly for such small geographical scales, we showed that selection is strong in all species but generally affects different loci in each. A clear signature of selection was systematically detected on a fraction of the genome, of the order of 0.1%-1% of the loci depending on the species. The novel modelling method we designed for estimating selection coefficients showed that the microgeographical selection coefficient scaled by population size (Ns) was 2-30. Our results convincingly suggest that selection maintains within-population diversity at microgeographical scales in spatially heterogeneous environments. Such genetic diversity is likely to be a major reservoir of adaptive potential, helping populations to adapt under fluctuating environmental conditions.

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.

Physiological and genetic control of transpiration efficiency in African rice, Oryza glaberrima Steud.

Improving crop water use efficiency, the amount of carbon assimilated as biomass per unit of water used by a plant, is of major importance as water for agriculture becomes scarcer. In rice, the genetic bases of transpiration efficiency, the derivation of water use efficiency at the whole-plant scale, and its putative component trait transpiration restriction under high evaporative demand remain unknown. These traits were measured in 2019 in a panel of 147 African rice (Oryza glaberrima) genotypes known to be potential sources of tolerance genes to biotic and abiotic stresses. Our results reveal that higher transpiration efficiency is associated with transpiration restriction in African rice. Detailed measurements in a subset of highly contrasted genotypes in terms of biomass accumulation and transpiration confirmed these associations and suggested that root to shoot ratio played an important role in transpiration restriction. Genome wide association studies identified marker-trait associations for transpiration response to evaporative demand, transpiration efficiency, and its residuals, with links to genes involved in water transport and cell wall patterning. Our data suggest that root-shoot partitioning is an important component of transpiration restriction that has a positive effect on transpiration efficiency in African rice. Both traits are heritable and define targets for breeding rice with improved water use strategies.

Interactions between microenvironment, selection and genetic architecture drive multiscale adaptation in a simulation experiment.

When environmental conditions differ both within and among populations, multiscale adaptation results from processes at both scales and interference across scales. We hypothesize that within-population environmental heterogeneity influences the chance of success of migration events, both within and among populations, and maintains within-population adaptive differentiation. We used a simulation approach to analyse the joint effects of environmental heterogeneity patterns, selection intensity and number of QTL controlling a selected trait on local adaptation in a hierarchical metapopulation design. We show the general effects of within-population environmental heterogeneity: (i) it increases occupancy rate at the margins of distribution ranges, under extreme environments and high levels of selection; (ii) it increases the adaptation lag in all environments; (iii) it impacts the genetic variance in each environment, depending on the ratio of within- to between-populations environmental heterogeneity; (iv) it reduces the selection-induced erosion of adaptive gene diversity. Most often, the smaller the number of QTL involved, the stronger are these effects. We also show that both within- and between-populations phenotypic differentiation (Q ST ) mainly results from covariance of QTL effects rather than QTL differentiation (F STq ), that within-population QTL differentiation is negligible, and that stronger divergent selection is required to produce adaptive differentiation within populations than among populations. With a high number of QTL, when the difference between environments within populations exceeds the smallest difference between environments across populations, high levels of within-population differentiation can be reached, reducing differentiation among populations. Our study stresses the need to account for within-population environmental heterogeneity when investigating local adaptation.

Genetic diversity and population structure in Vitis species illustrate phylogeographic patterns in eastern North America.

Geographical distribution and diversity of current plant species have been strongly shaped by climatic oscillations during the Quaternary. Analysing the resulting divergence among species and differentiation within species is crucial to understand the evolution of taxa like the Vitis genus, which provides very useful genetic resources for grapevine improvement and might reveal original recolonization patterns due to growth habit and dispersal mode. Here, we studied the genetic structure in natural populations of three species from eastern North America: Vitis aestivalis, V. cinerea and V. riparia using different marker types. Vitis aestivalis and V. cinerea showed higher diversity than V. riparia. The two former species are less differentiated, confirming an earlier divergence of V. riparia. V. aestivalis and V. riparia exhibited different genetic groups on both sides of the Appalachian Mountains that could mirror different recolonization routes from southern refugia. Genetic structure was stronger in V. cinerea, for which two varieties (var. berlandieri and var. cinerea) are morphologically recognized. Our results confirm this distinction and suggest the existence of three other lineages within var. cinerea. These discontinuities appear linked to adaptation of var. berlandieri to dry and limy areas of Texas and partially to the Mississippi River Valley. Rapid range expansions from refugia upon climate warming are also suggested by the low linkage disequilibrium values observed. Furthermore, large variation for downy mildew resistance was observed in the three species. Our findings appeared consistent with the vegetation history of eastern North America.

QTL detection for growth and latex production in a full-sib rubber tree population cultivated under suboptimal climate conditions.

BACKGROUND: Rubber tree is cultivated in mainly Southeast Asia and is by far the most significant source of natural rubber production worldwide. However, the genetic architecture underlying the primary agronomic traits of this crop has not been widely characterized. This study aimed to identify quantitative trait loci (QTLs) associated with growth and latex production using a biparental population established in suboptimal growth conditions in Brazil. RESULTS: A full-sib population composed of 251 individuals was developed from crossing two high-producing Asiatic rubber tree cultivars, PR 255 and PB 217. This mapping population was genotyped with microsatellite markers from enriched genomic libraries or transcriptome datasets and single-nucleotide polymorphism (SNP) markers, leading to construction of a saturated multipoint integrated genetic map containing 354 microsatellite and 151 SNP markers. Height and circumference measurements repeated over a six-year period and registration of cumulative latex production during six consecutive months on the same individuals allowed in-depth characterization of the genetic values of several growth traits and precocious latex production. Growth traits, circumference and height, were overall positively correlated, whereas latex production was not correlated or even negatively correlated with growth traits. A total of 86 distinct QTLs were identified, most of which were detected for only one trait. Among these QTLs, 15 were linked to more than one phenotypic trait (up to 4 traits simultaneously). Latex production and circumference increments during the last wintering period were associated with the highest numbers of identified QTLs (eleven and nine, respectively), jointly explaining the most significantly observed phenotypic variances (44.1% and 44.4%, respectively). The most important QTL for latex production, located on linkage group 16, had an additive effect of the male parent PB 217 and corresponded to a QTL at the same position detected in a previous study carried out in Thailand for the biparental population RRIM 600 x PB 217. CONCLUSIONS: Our results identified a set of significant QTLs for rubber tree, showing that the performance of modern Asiatic cultivars can still be improved and paving the way for further marker-assisted selection, which could accelerate breeding programs.

Developing core collections to optimize the management and the exploitation of diversity of the coffee Coffea canephora.

The management of diversity for conservation and breeding is of great importance for all plant species and is particularly true in perennial species, such as the coffee Coffea canephora. This species exhibits a large genetic and phenotypic diversity with six different diversity groups. Large field collections are available in the Ivory Coast, Uganda and other Asian, American and African countries but are very expensive and time consuming to establish and maintain in large areas. We propose to improve coffee germplasm management through the construction of genetic core collections derived from a set of 565 accessions that are characterized with 13 microsatellite markers. Core collections of 12, 24 and 48 accessions were defined using two methods aimed to maximize the allelic diversity (Maximization strategy) or genetic distance (Maximum-Length Sub-Tree method). A composite core collection of 77 accessions is proposed for both objectives of an optimal management of diversity and breeding. This core collection presents a gene diversity value of 0.8 and exhibits the totality of the major alleles (i.e., 184) that are present in the initial set. The seven proposed core collections constitute a valuable tool for diversity management and a foundation for breeding programs. The use of these collections for collection management in research centers and breeding perspectives for coffee improvement are discussed.

Molecular characterisation of Pinus sylvestris (L.) in Ireland at the western limit of the species distribution.

BACKGROUND: Scots pine (Pinus sylvestris L.) underwent significant population declines across much of northwest Europe during the mid-to-late Holocene and was thought to have become extirpated in Ireland from about 400 AD. However, most extant populations are plantations reintroduced from Scotland. Others are naturalised therefrom and one in Western Ireland is a putative relict. In this paper, Scots pine in Ireland are genetically described for the first time. RESULTS: Using two mitochondrial (mtDNA) loci, eight chloroplast (cpSSR) and 18 nuclear (nSSR) loci, the genetic composition and diversity of 19 Irish Scots pine populations is described and compared to other European populations. All trees sampled in Ireland were fixed for mitotype a, which is the most common across northwest Europe. By contrast, cpSSR (HCP = 0.967) and nSSR (He = 0.540) variation was high, and comparable with estimates for other regions across the species range. Differentiation at both sets of loci were similarly low (cpSSR FST = 0.019; nSSR FST = 0.018), but populations from continental Europe were significantly differentiated from all Irish populations based on nSSR variation. CONCLUSIONS: All Irish Scots pine are likely part of a common Irish-Scottish gene pool which diverged from continental Scots pine following post-glacial recolonisation. A high genetic diversity and an absence of evidence of inbreeding suggests the regional decline of Scots pine did not critically reduce allelic variation. The post-glacial relationship between Irish and Scottish pine is discussed, and a suggestion from recent palaeoecological work that reintroduced Scots pine be managed as a native species is now further supported by genetic data.

Genomic evidence of genuine wild versus admixed olive populations evolving in the same natural environments in western Mediterranean Basin.

Crop-to-wild gene flow is a mechanism process widely documented, both in plants and animals. This can have positive or negative impacts on the evolution of admixed populations in natural environments, yet the phenomenon is still misunderstood in long-lived woody species, contrary to short-lived crops. Wild olive Olea europaea L. occurs in the same eco-geographical range as domesticated olive, i.e. the Mediterranean Basin (MB). Moreover, it is an allogamous and anemophilous species whose seeds are disseminated by birds, i.e. factors that drive gene flow between crops and their wild relatives. Here we investigated the genetic structure of western MB wild olive populations in natural environments assuming a homogenous gene pool with limited impact of cultivated alleles, as previously suggested. We used a target sequencing method based on annotated genes from the Farga reference genome to analyze 27 western MB olive tree populations sampled in natural environments in France, Spain and Morocco. We also target sequenced cultivated olive tree accessions from the Worldwide Olive Germplasm Bank of Marrakech and Porquerolles and from an eastern MB wild olive tree population. We combined PCA, sNMF, pairwise FST and TreeMix and clearly identified genuine wild olive trees throughout their natural distribution range along a north-south gradient including, for the first time, in southern France. However, contrary to our assumption, we highlighted more admixed than genuine wild olive trees. Our results raise questions regarding the admixed population evolution pattern in this environment, which might be facilitated by crop-to-wild gene flow.

An initial assessment of linkage disequilibrium (LD) in coffee trees: LD patterns in groups of Coffea canephora Pierre using microsatellite analysis.

BACKGROUND: A reciprocal recurrent selection program has been under way for the Coffea canephora coffee tree for approximately thirty years in the Ivory Coast. Association genetics would help to speed up this program by more rapidly selecting zones of interest in the genome. However, prior to any such studies, the linkage disequilibrium (LD) needs to be assessed between the markers on the genome. These data are essential for guiding association studies. RESULTS: This article describes the first results of an LD assessment in a coffee tree species. Guinean and Congolese breeding populations of C. canephora have been used for this work, with the goal of identifying ways of using these populations in association genetics. We identified changes in the LD along the genome within the different C. canephora diversity groups. In the different diversity groups studied, the LD was variable. Some diversity groups displayed disequilibria over long distances (up to 25 cM), whereas others had disequilibria not exceeding 1 cM. We also discovered a fine structure within the Guinean group. CONCLUSIONS: Given these results, association studies can be used within the species C. canephora. The coffee recurrent selection scheme being implemented in the Ivory Coast can thus be optimized. Lastly, our results could be used to improve C. arabica because one of its parents is closely related to C. canephora.

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

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

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

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

Impact of Seed Origin and Genetic Drift of Improved Rice Variety IR841 in Benin.

BACKGROUND: Rice production is important for food security in Benin, with a national production largely dominated by the cultivation of the aromatic ecotype IR841, by far the most appreciated by Beninese consumers. However, to maintain agronomical qualities of a given cultivar, the origin and quality of seeds are mandatory and at the heart of the maintenance of yield and of market requests. Following this idea, the objective of the current study was thus to investigate the genetic purity of IR841 varieties collected from rice producers across 22 villages in Benin in relation to their agronomical performances. RESULTS: For this, agromorphological evaluation of 72 accessions based on 13 quantitative descriptors followed by genotyping with the Illumina Infinium rice 7 K SNP array of 9 accessions was carried out in the presence of 2 controls. Agromorphologic as well as genetic and phylogenetic analyses revealed two groups, the first one Okouta97, Koum47, Nana30, Man118, Ang1 and control sample IR841-2) was characterized by seed accessions provided by the formal seed system, while the second (Koum53, Tchaka41 and Koud46) comprising seeds from local markets or from previous harvests and showing a depression in agronomic performances. CONCLUSION: We showed that IR841 seed purity is mandatory for the completion of agronomical performance, and that the farmers' choice of seeds must be guided and informed to ensure sustainability and food security.

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.

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.

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.

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.