Jean-François Rami

Plant breeding is crucial to develop varieties able to cope with climate change and support food and feed value chains. Genomic prediction (GP) has been a major step in increasing their efficiency and recently, phenomic prediction (PP) has gained attention as a promising complementary approach to GP, potentially further increasing this efficiency. Factors impacting PP are not fully clarified. Thus, we studied the impacts of spectra preprocessing, prediction methods, population structure, training set size, near infrared reflectance spectroscopy (NIRS) acquisition environment, and wavelength selection on a large multi‐parental sorghum population including 2498 BC1F3:5 families from 29 crosses with a strong population structure. Using 51,545 single nucleotide polymorphisms and 1154 NIRS features, we show that PP can reach predictive abilities (PAs) similar to GP, that it is less affected by population structure, and can reach its maximal PA with smaller training sets than GP, but its performances are trait dependent. We also show that NIRS can be acquired in a reference environment to perform prediction in other environments and that it is possible to randomly select wavelengths to perform predictions. Finally, we show that spectra preprocessing and statistical methods have an inconsistent impact on PA. Our study confirms that PP is a relevant trait prediction method that deserves attention to optimize breeding schemes. The main challenges for the future will be to better understand the information contained in the spectra and disentangle their genetic and proxy components to optimize the use of PP in breeding programs.

While the green revolution adapted a handful of crops to homogenous and high-input industrialized agriculture, much of the global population still relies on local food production from low-input smallholder farms that grow highly variable crop cultivars. The high diversity of the grain and bioenergy crop sorghum, and many other crops that were not homogenized during the green revolution, not only provides the raw materials for breeders to make substantial gains in cultivar improvement, but also constrains breeding efforts due to highly specialized locally adapted plant phenotypes. Here, we construct a 33-member pangenome and identify trait-associated variants in 1,988 cultivars and landraces. We then apply these resources to explore the complex interplay between historical contingency, ongoing adaptation, and the potential for future gains through climate-aware genome-enabled breeding. Specifically, our analyses conclusively demonstrate that multiple nested, deeply diverged, and previously uncharacterized structural variants in the domestication gene SHATTERING1 distinguish the previously established multicentric origin of sorghum. We then apply landscape genomics tests to reveal how gene flow, adaptation, and secondary contact created the complex genetic mosaic in current global breeding networks. Further analysis of climate-gene associations highlights candidate loci underlying adaptation, including the biosynthetic gene cluster for the cyanogenic glucoside dhurrin. Combined, the pangenome-informed variants developed here will enable both trait discovery and subsequent marker assays to accelerate breeding and provide a framework for similar applications in other diverse and non-model crops.

Plant breeding efficiency is crucial to develop varieties able to cope with climate change and support food and feed value chains. Genomic prediction (GP) has been a major step in increasing this efficiency and is now routinely used in breeding programs. Recently, phenomic prediction (PP) has gained attention as a promising complementary approach to GP, further increasing the breeding programs’ efficiency. Factors impacting the predictive ability (PA) of PP have been studied on many species but are not fully clarified. In this context, we studied the impacts of spectra pre-processing, prediction methods, population structure, training set size, NIRS acquisition environment and wavelength selection on a large multi-parental sorghum population including 2498 genotypes. Our results show that PP can compete with GP, that it is less affected by population structure, and can reach its maximal PA with smaller training sets than GP, but its performances are trait dependant. We also show that NIRS can be acquired in a reference environment to perform prediction in other environments and that it is possible to randomly select as little as 10 wavelengths to perform predictions. Finally, we show that spectra pre-processing, and statistical methods have a limited and unclear impact on PA. Our study confirms that PP is a relevant trait prediction method that deserves attention to optimize breeding schemes. The main challenges for the future will be to better understand the information contained in the spectra and disentangle their genetic and proxy components to optimize the use of PP in breeding programs. Key message Phenomic prediction is promising for sorghum breeding. Geneticists’ methods may not be suited to optimally extract spectral information.

Motivation: Genotyping of bi-parental populations can be performed with low-coverage next-generation sequencing (LC-NGS). This allows the creation of highly saturated genetic maps at reasonable cost, precisely localized recombination breakpoints (i.e., the crossovers), and minimized mapping intervals for quantitative-trait locus analysis. The main issues with these low-coverage genotyping methods are (1) poor performance at heterozygous loci, (2) high percentage of missing data, (3) local errors due to erroneous mapping of sequencing reads and reference genome mistakes, and (4) global, technical errors inherent to NGS itself. Recent methods like Tassel-FSFHap or LB-Impute are excellent at addressing issues 1 and 2, but nonetheless perform poorly when issues 3 and 4 are persistent in a dataset (i.e., “noisy” data). Here, we present a new algorithm for imputation of LC-NGS data that eliminates the need of complex pre-filtering of noisy data, accurately types heterozygous chromosomal regions, precisely estimates crossover positions, corrects erroneous data, and imputes missing data. The imputation of genotypes and recombination breakpoints is based on maximum-likelihood estimation. We compare its performance with Tassel-FSFHap and LB-Impute using simulated data and two real datasets. NOISYmputer is consistently more efficient than the two other software tested and reaches average breakpoint precision of 99.9% and average recall of 99.6% on illumina simulated dataset. NOISYmputer consistently provides precise map size estimations when applied to real datasets while alternative tools may exhibit errors ranging from 3 to 1845 times the real size of the chromosomes in centimorgans. Furthermore, the algorithm is not only highly effective in terms of precision and recall but is also particularly economical in its use of RAM and computation time, being much faster than Hidden Markov Model methods. Availability: NOISYmputer and its source code are available as a multiplatform (Linux, macOS, Windows) Java executable at the URL https://gitlab.cirad.fr/noisymputer/noisymputerstandalone/-/tree/1.0.0-RELEASE?reftype=tags.

Population growth and the impacts of climate change are placing increasing pressure on global agriculture and breeding programmes. Recent advancements in phenotyping techniques, genotyping technologies, and predictive modelling are accelerating genetic gains in breeding programmes, helping researchers and breeders develop improved crops more efficiently. However, these advancements have also led to an overwhelming torrent of fragmented data, creating significant challenges in data integration and management. To address this issue, the Breeding Application Programming Interface (BrAPI) project was established as a standardized data model for breeding data. BrAPI is an international, community-driven effort that facilitates interoperability among databases and tools, improving the sharing and interpretation of breeding-related data. This open-source standard is software-agnostic and can be used by anyone interested in breeding, phenotyping, germplasm, genotyping, and agronomy data management. This manuscript provides an overview of the BrAPI project, highlighting the significant progress made in the development of the data standard and the expansion of its community. It also presents a showcase of the wide variety of BrAPI-compatible tools that have been built to enhance breeding and research activities, demonstrating how the project is advancing agricultural innovation and data management practices.

Motivation Bi-parental populations genotyping can be performed with low-coverage next-generation sequencing (LC-NGS). This allows the creation of highly saturated genetic maps at reasonable cost, precisely localized recombination breakpoints, and minimized mapping intervals for quantitative-trait locus analysis. The main issues with these low-coverage genotyping methods are (1) poor performance at heterozygous loci, (2) a high percentage of missing data, (3) local errors due to erroneous mapping of sequencing reads and reference genome mistakes, and (4) global, technical errors inherent to NGS itself. Recent methods like Tassel-FSFHap or LB-Impute are excellent at addressing issues 1 and 2, but nonetheless perform poorly when issues 3 and 4 are persistent in a dataset (i.e., “noisy” data). Here, we present an algorithm for imputation of LC-NGS data that eliminates the need of complex pre-filtering of noisy data, accurately types heterozygous chromosomal regions, precisely estimates crossover positions, corrects erroneous data, and imputes missing data. The imputation of genotypes and recombination breakpoints is based on maximum-likelihood estimation. We compare its performance with Tassel-FSFHap and LB-Impute using simulated data and two real datasets. Furthermore, the algorithm is much faster than the Hidden Markov-Chains method. Availability NOISYmputer is available as a multiplatform (Linux, macOS, Windows) Java executable at the URL https://gitlab.cirad.fr/noisymputer/noisymputerstandalone/-/tree/1.0.0-RELEASE?ref_type=tags. The source code is available at the same URL.

Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.

Cultivated peanut or groundnut (Arachis hypogaea L.) is a grain legume grown in many developing countries by smallholder farmers for food, feed, and/or income. The speciation of the cultivated species, that involved polyploidization followed by domestication, greatly reduced its variability at the DNA level. Mobilizing peanut diversity is a prerequisite for any breeding program for overcoming the main constraints that plague production and for increasing yield in farmer fields. In this study, the Groundnut Improvement Network for Africa assembled a collection of 1,049 peanut breeding lines, varieties, and landraces from 9 countries in Africa. The collection was genotyped with the Axiom_Arachis2 48K SNP array and 8,229 polymorphic single nucleotide polymorphism (SNP) markers were used to analyze the genetic structure of this collection and quantify the level of genetic diversity in each breeding program. A supervised model was developed using dapc to unambiguously assign 542, 35, and 172 genotypes to the Spanish, Valencia, and Virginia market types, respectively. Distance-based clustering of the collection showed a clear grouping structure according to subspecies and market types, with 73% of the genotypes classified as fastigiata and 27% as hypogaea subspecies. Using STRUCTURE, the global structuration was confirmed and showed that, at a minimum membership of 0.8, 76% of the varieties that were not assigned by dapc were actually admixed. This was particularly the case of most of the genotype of the Valencia subgroup that exhibited admixed genetic heritage. The results also showed that the geographic origin (i.e. East, Southern, and West Africa) did not strongly explain the genetic structure. The gene diversity managed by each breeding program, measured by the expected heterozygosity, ranged from 0.25 to 0.39, with the Niger breeding program having the lowest diversity mainly because only lines that belong to the fastigiata subspecies are used in this program. Finally, we developed a core collection composed of 300 accessions based on breeding traits and genetic diversity. This collection, which is composed of 205 genotypes of fastigiata subspecies (158 Spanish and 47 Valencia) and 95 genotypes of hypogaea subspecies (all Virginia), improves the genetic diversity of each individual breeding program and is, therefore, a unique resource for allele mining and breeding.

Quantitative Trait Loci (QTL) mapping has been thoroughly used in peanut genetics and breeding in spite of the narrow genetic diversity and the segmental tetraploid nature of the cultivated species. QTL mapping is helpful for identifying the genomic regions that contribute to traits, for estimating the extent of variation and the genetic action (i.e., additive, dominant, or epistatic) underlying this variation, and for pinpointing genetic correlations between traits. The aim of this paper is to review the recently published studies on QTL mapping with a particular emphasis on mapping populations used as well as traits related to kernel quality. We found that several populations have been used for QTL mapping including interspecific populations developed from crosses between synthetic tetraploids and elite varieties. Those populations allowed the broadening of the genetic base of cultivated peanut and helped with the mapping of QTL and identifying beneficial wild alleles for economically important traits. Furthermore, only a few studies reported QTL related to kernel quality. The main quality traits for which QTL have been mapped include oil and protein content as well as fatty acid compositions. QTL for other agronomic traits have also been reported. Among the 1261 QTL reported in this review, and extracted from the most relevant studies on QTL mapping in peanut, 413 (~33%) were related to kernel quality showing the importance of quality in peanut genetics and breeding. Exploiting the QTL information could accelerate breeding to develop highly nutritious superior cultivars in the face of climate change.

Peanuts (Arachis hypogaea L.) are an allotetraploid grain legume mainly cultivated by poor farmers in Africa, in degraded soil and with low input systems. Further understanding nodulation genetic mechanisms could be a relevant option to facilitate the improvement of yield and lift up soil without synthetic fertilizers. We used a subset of 83 chromosome segment substitution lines (CSSLs) derived from the cross between a wild synthetic tetraploid AiAd (Arachis ipaensis × Arachis duranensis)^4× and the cultivated variety Fleur11, and evaluated them for traits related to BNF under shade-house conditions. Three treatments were tested: without nitrogen; with nitrogen; and without nitrogen, but with added0 Bradyrhizobium vignae strain ISRA400. The leaf chlorophyll content and total biomass were used as surrogate traits for BNF. We found significant variations for both traits specially linked to BNF, and four QTLs (quantitative trait loci) were consistently mapped. At all QTLs, the wild alleles decreased the value of the trait, indicating a negative effect on BNF. A detailed characterization of the lines carrying those QTLs in controlled conditions showed that the QTLs affected the nitrogen fixation efficiency, nodule colonization, and development. Our results provide new insights into peanut nodulation mechanisms and could be used to target BNF traits in peanut breeding programs.