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

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

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.