Genome scans capture key adaptation and historical hybridization signatures in tetraploid wheat.

Tetraploid wheats (Triticum turgidum L.), including durum wheat (T. turgidum ssp. durum (Desf.) Husn.), are important crops with high nutritional and cultural values. However, their production is constrained by sensitivity to environmental conditions. In search of adaptive genetic signatures tracing historical selection and hybridization events, we performed genome scans on two datasets: (1) Durum Global Diversity Panel comprising a total of 442 tetraploid wheat and wild progenitor accessions including durum landraces (n = 286), domesticated emmer (T. turgidum ssp. dicoccum (Schrank) Thell.; n = 103) and wild emmer (T. turgidum ssp. dicoccoides (Korn. ex Asch. & Graebn.) Thell.; n = 53) wheats genotyped using the 90K single nucleotide polymorphism (SNP) array, and (2) a second dataset comprising a total 121 accessions of nine T. turgidum subspecies including wild emmer genotyped with >100 M SNPs from whole-genome resequencing. The genome scan on the first dataset detected six outlier loci on chromosomes 1A, 1B, 3A (n = 2), 6A, and 7A. These loci harbored important genes for adaptation to abiotic stresses, phenological responses, such as seed dormancy, circadian clock, flowering time, and key yield-related traits, including pleiotropic genes, such as HAT1, KUODA1, CBL1, and ZFN1. The scan on the second dataset captured a highly differentiated region on chromosome 2B that shows significant differentiation between two groups: one group consists of Georgian (T. turgidum ssp. paleocolchicum A. Love & D. Love) and Persian (T. turgidum ssp. carthlicum (Nevski) A. Love & D. Love) wheat accessions, while the other group comprises all the remaining tetraploids including wild emmer. This is consistent with a previously reported introgression in this genomic region from T. timopheevii Zhuk. which naturally cohabit in the Georgian and neighboring areas. This region harbored several adaptive genes, including the thermomorphogenesis gene PIF4, which confers temperature-resilient disease resistance and regulates other biological processes. Genome scans can be used to fast-track germplasm housed in gene banks and in situ; which helps to identify environmentally resilient accessions for breeding and/or to prioritize them for conservation.

Multi-locus genome-wide association studies reveal the genetic architecture of Fusarium head blight resistance in durum wheat.

Durum wheat is more susceptible to Fusarium head blight (FHB) than other types or classes of wheat. The disease is one of the most devastating in wheat; it reduces yield and end-use quality and contaminates the grain with fungal mycotoxins such as deoxynivalenol (DON). A panel of 265 Canadian and European durum wheat cultivars, as well as breeding and experimental lines, were tested in artificially inoculated field environments (2019-2022, inclusive) and two greenhouse trials (2019 and 2020). The trials were assessed for FHB severity and incidence, visual rating index, Fusarium-damaged kernels, DON accumulation, anthesis or heading date, maturity date, and plant height. In addition, yellow pigment and protein content were analyzed for the 2020 field season. To capture loci underlying FHB resistance and related traits, GWAS was performed using single-locus and several multi-locus models, employing 13,504 SNPs. Thirty-one QTL significantly associated with one or more FHB-related traits were identified, of which nine were consistent across environments and associated with multiple FHB-related traits. Although many of the QTL were identified in regions previously reported to affect FHB, the QTL QFhb-3B.2, associated with FHB severity, incidence, and DON accumulation, appears to be novel. We developed KASP markers for six FHB-associated QTL that were consistently detected across multiple environments and validated them on the Global Durum Panel (GDP). Analysis of allelic diversity and the frequencies of these revealed that the lines in the GDP harbor between zero and six resistance alleles. This study provides a comprehensive assessment of the genetic basis of FHB resistance and DON accumulation in durum wheat. Accessions with multiple favorable alleles were identified and will be useful genetic resources to improve FHB resistance in durum breeding programs through marker-assisted recurrent selection and gene stacking.

Historical Selection, Adaptation Signatures, and Ambiguity of Introgressions in Wheat.

Wheat was one of the crops domesticated in the Fertile Crescent region approximately 10,000 years ago. Despite undergoing recent polyploidization, hull-to-free-thresh transition events, and domestication bottlenecks, wheat is now grown in over 130 countries and accounts for a quarter of the world's cereal production. The main reason for its widespread success is its broad genetic diversity that allows it to thrive in different environments. To trace historical selection and hybridization signatures, genome scans were performed on two datasets: approximately 113K SNPs from 921 predominantly bread wheat accessions and approximately 110K SNPs from about 400 wheat accessions representing all ploidy levels. To identify environmental factors associated with the loci, a genome-environment association (GEA) was also performed. The genome scans on both datasets identified a highly differentiated region on chromosome 4A where accessions in the first dataset were dichotomized into a group (n = 691), comprising nearly all cultivars, wild emmer, and most landraces, and a second group (n = 230), dominated by landraces and spelt accessions. The grouping of cultivars is likely linked to their potential ancestor, bread wheat cv. Norin-10. The 4A region harbored important genes involved in adaptations to environmental conditions. The GEA detected loci associated with latitude and temperature. The genetic signatures detected in this study provide insight into the historical selection and hybridization events in the wheat genome that shaped its current genetic structure and facilitated its success in a wide spectrum of environmental conditions. The genome scans and GEA approaches applied in this study can help in screening the germplasm housed in gene banks for breeding, and for conservation purposes.

Loci harboring genes with important role in drought and related abiotic stress responses in flax revealed by multiple GWAS models.

KEY MESSAGE: QTNs associated with drought tolerance traits and indices were identified in a flax mini-core collection through multiple GWAS models and phenotyping at multiple locations under irrigated and non-irrigated field conditions. Drought is a critical phenomenon challenging today's agricultural sector. Crop varieties adapted to moisture deficit are becoming vital. Flax can be greatly affected by limiting moisture conditions, especially during the early development and reproductive stages. Here, a mini-core collection comprising genotypes from more than 20 major growing countries was evaluated for 11 drought-related traits in irrigated and non-irrigated fields for 3 years. Heritability of the traits ranged from 44.7 to 86%. Six of the 11 traits showed significant phenotypic difference between irrigated and non-irrigated conditions. A genome-wide association study (GWAS) was performed for these six traits and their corresponding stress indices based on 106 genotypes and 12,316 single nucleotide polymorphisms (SNPs) using six multi-locus and one single-locus models. The SNPs were then assigned to 8050 linkage disequilibrium (LD) blocks to which a restricted two-stage multi-locus multi-allele GWAS was applied. A total of 144 quantitative trait nucleotides (QTNs) and 13 LD blocks were associated with at least one trait or stress index. Of these, 16 explained more than 15% of the genetic variance. Most large-effect QTN loci harbored gene(s) previously predicted to play role(s) in the associated traits. Genes mediating responses to abiotic stresses resided at loci associated with stress indices. Flax genes Lus10009480 and Lus10030150 that are predicted to encode WAX INDUCER1 and STRESS-ASSOCIATED PROTEIN (SAP), respectively, are among the important candidates detected. Accessions with multiple favorable alleles outperformed others for grain yield, thousand seed weight and fiber/biomass in non-irrigated conditions, suggesting their potential usefulness in breeding and genomic selection.

Genome-wide signatures in flax pinpoint to adaptive evolution along its ecological gradient.

BACKGROUND: Flax is one of the eight founder crops of agriculture. It is believed to have been domesticated as a long-day plant that has since spread to survive in a wide range of eco-geographic regions extending from the warm Indian subcontinent to the low latitude east African highlands and to the cool and high-latitude Eurasia. Understanding the genetic basis underlying its adaptation and selection events throughout its dispersion is essential to develop cultivars adapted to local environmental conditions. METHODS: Here we detected genetic signatures of local adaptation and selection events of flax based on 385 accessions from all major flax growing regions of the world using genome scan methods and three genomic datasets: (1) a genome-wide dataset of more than 275K single nucleotide polymorphisms (SNPs), (2) a filtered dataset of 23K SNPs with minor allele frequency >10% and, (3) a 34K exon-derived SNP dataset. RESULTS: Principal component (PC) and fixation index (FS⁢T)-based genome scans yielded consistent outlier SNP loci on chromosomes 1, 8, 9 and 12. Additional loci on chromosomes 3, 7, 8, 10, 11, 13 and 14 were detected using both the PC and FS⁢T methods in two of the three datasets. A genome-environment association (GEA) analysis using the 23K dataset and the first PC of cropping season temperature, day-length and latitude identified significant SNPs on chromosomes 3, 7, 9 and 13. CONCLUSIONS: Most of the loci detected by the three methods harbored relevant genes for local adaptation, including some that play roles in day-length, light and other biotic and abiotic stresses responses. Such genetic signatures may help to select pre-breeding materials potentially adapted to specific growing niches prior to field performance trials. Given the current low genotyping cost and freely available environmental data, the genome scans along with GEA can readily provide opportunity to sort out materials suitable to various environmental conditions from large set of germplasm in gene banks and/or in situ, thereby assisting the breeding and genetic conservation efforts.

The Complex Genetic Architecture of Early Root and Shoot Traits in Flax Revealed by Genome-Wide Association Analyses.

Roots are fundamental organs for water and nutrient uptake as well as for signal transduction in response to biotic and abiotic stresses. Flax has a shallow tap root system that relies mostly on top soil nutrient and moisture resources. The crop can easily be outcompeted by weeds or other crops in intercropping systems, especially in moisture deficit conditions. However, there is a wide range of variation among genotypes in terms of performance under scarce resources such as moisture limitation. Here we phenotyped 15 root, two shoot traits and shoot to root dry weight ratio on 115 flax accessions grown in a hydroponic pouch system and performed a genome-wide association study (GWAS) based on seven different models to identify quantitative trait loci underlying these traits. Significant variation among genotypes was observed for the two shoot and 12 of the 14 root traits. Shoot dry weight was correlated with root network volume, length, surface area, and root dry weight (r > 0.5, P < 0.001) but not significantly correlated with root depth (r = 0.033, P > 0.05). The seven GWAS models detected a total of 228 quantitative trait nucleotides (QTNs) for 16 traits. Most loci, defined by an interval of 100 kb up and downstream of the QTNs, harbored genes known to play role(s) in root and shoot development, suggesting them as candidates. Examples of candidate genes linked to root network QTNs included genes encoding GRAS transcription factors, mitogen-activated protein kinases, and auxin related lateral organ boundary proteins while QTN loci for shoot dry weight harbored genes involved in photomorphogenesis and plant immunity. These results provide insights into the genetic bases of early shoot and root development traits in flax that could be capitalized upon to improve its root architecture, particularly in view of better withstanding water limiting conditions during the cropping season.

The genetic structure of flax illustrates environmental and anthropogenic selections that gave rise to its eco-geographical adaptation.

Flax, one of the eight founder crops of agriculture, has been cultivated for its oil and/or fiber for millennia. Understanding genetic divergence and geographic origins of germplasm in line with their cultivation history and ecological adaptation are essential for conservation and breeding. Here we performed a genome-wide assessment based on more than 51,000 single nucleotide polymorphic sites defining 383 flax accessions from a core collection representing 37 flax growing countries. Population structure analysis resulted in a total of 12 populations that were pooled into four major groups: Temperate, South Asian, Abyssinian and Mediterranean. The vast majority (n = 335) belonged to the Temperate group that comprised eight populations including one dominated by fiber flax. Genetic variation between fiber and oil morphotypes was less pronounced than variation within morphotypes. The genetic variation among groups and populations was attributed in part to eco-geographic and anthropogenic factors. Genetic signatures indicated loci under strong selection by environmental factors such as day length. A high concentration of private haplotypes were observed in the South Asian, Mediterranean and Abyssinian populations despite their low genotype representation, hinting at the long history of the crop in these regions. The addition of genotypes from these three regions would enrich the core collection by capturing a wider genetic breadth for breeding and conservation.