Harnessing translational research in wheat for climate resilience.

Despite being the world's most widely grown crop, research investments in wheat (Triticum aestivum and Triticum durum) fall behind those in other staple crops. Current yield gains will not meet 2050 needs, and climate stresses compound this challenge. However, there is good evidence that heat and drought resilience can be boosted through translating promising ideas into novel breeding technologies using powerful new tools in genetics and remote sensing, for example. Such technologies can also be applied to identify climate resilience traits from among the vast and largely untapped reserve of wheat genetic resources in collections worldwide. This review describes multi-pronged research opportunities at the focus of the Heat and Drought Wheat Improvement Consortium (coordinated by CIMMYT), which together create a pipeline to boost heat and drought resilience, specifically: improving crop design targets using big data approaches; developing phenomic tools for field-based screening and research; applying genomic technologies to elucidate the bases of climate resilience traits; and applying these outputs in developing next-generation breeding methods. The global impact of these outputs will be validated through the International Wheat Improvement Network, a global germplasm development and testing system that contributes key productivity traits to approximately half of the global wheat-growing area.

Mapping for Adult-Plant Resistance Against Septoria Tritici Blotch in a Common Wheat Line Murga.

Septoria tritici blotch (STB) is a major foliar disease globally that is notorious for quickly developing fungicide resistance, making host resistance an indispensable component in mitigating STB. The International Maize and Wheat Improvement Center (CIMMYT) wheat line Murga is well known for its high, durable, and broad-spectrum resistance against STB infection. This study aimed to investigate the resistance mechanism of Murga to facilitate its utilization in breeding. A recombinant inbred line population was derived from a cross between Murga and STB-susceptible line Huirivis#1, comprising 297 progenies. The population was evaluated for adult-plant STB resistance in Toluca, Mexico (from 2017 to 2019), and in La Estanzuela, Uruguay (from 2016 to 2018). Genotyping was performed with the DArTseq platform. Quantitative trait locus (QTL) mapping indicated a major and stable QTL on chromosome 3DL, explaining a phenotypic variation for STB of 41.2 to 62.5% in Mexico and 27.5 to 40.3% in Uruguay. This QTL was regarded as Stb16 based on the comparison of its physical position, the possible origin from synthetic wheat, and its broad-spectrum resistance. Additional QTLs with minor effects were identified on chromosomes 2B, 2D, 3A, 3B, and 5B. The QTL on 5BS was significant in four of the six environments and must be new. Murga was the resistant donor for all QTLs except for those on 2B and 3A. Being an elite breeding line, Stb16 carrier Murga could be used as a promising STB resistance donor. Rational employment of Stb16 could contribute to STB management yet avoid the rapid emergence of Stb16-virulent isolates.

Genomic analysis of Spanish wheat landraces reveals their variability and potential for breeding.

BACKGROUND: One of the main goals of the plant breeding in the twenty-first century is the development of crop cultivars that can maintain current yields in unfavorable environments. Landraces that have been grown under varying local conditions include genetic diversity that will be essential to achieve this objective. The Center of Plant Genetic Resources of the Spanish Institute for Agriculture Research maintains a broad collection of wheat landraces. These accessions, which are locally adapted to diverse eco-climatic conditions, represent highly valuable materials for breeding. However, their efficient use requires an exhaustive genetic characterization. The overall aim of this study was to assess the diversity and population structure of a selected set of 380 Spanish landraces and 52 reference varieties of bread and durum wheat by high-throughput genotyping. RESULTS: The DArTseq GBS approach generated 10 K SNPs and 40 K high-quality DArT markers, which were located against the currently available bread and durum wheat reference genomes. The markers with known locations were distributed across all chromosomes with relatively well-balanced genome-wide coverage. The genetic analysis showed that the Spanish wheat landraces were clustered in different groups, thus representing genetic pools providing a range of allelic variation. The subspecies had a major impact on the population structure of the durum wheat landraces, with three distinct clusters that corresponded to subsp. durum, turgidum and dicoccon being identified. The population structure of bread wheat landraces was mainly biased by geographic origin. CONCLUSIONS: The results showed broader genetic diversity in the landraces compared to a reference set that included commercial varieties, and higher divergence between the landraces and the reference set in durum wheat than in bread wheat. The analyses revealed genomic regions whose patterns of variation were markedly different in the landraces and reference varieties, indicating loci that have been under selection during crop improvement, which could help to target breeding efforts. The results obtained from this work will provide a basis for future genome-wide association studies.

The genome of Eucalyptus grandis.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Quantitative Trait Loci Mapping for Spot Blotch Resistance in Two Biparental Mapping Populations of Bread Wheat.

Spot blotch (SB), caused by Bipolaris sorokiniana, is a major fungal disease of wheat in South Asia and South America. Two biparental mapping populations with 232 F2:7 progenies each were generated, with CIMMYT breeding lines CASCABEL and KATH as resistant parents and CIANO T79 as the common susceptible parent. The two populations were evaluated for field SB resistance in CIMMYT's Agua Fria station for three consecutive cropping seasons, with artificial inoculation. Genotyping was done with the DArTseq platform and approximately 1,500 high quality and nonredundant markers were used for quantitative trait loci (QTL) mapping. In both populations, a major QTL was found on chromosome 5A in the Vrn-A1 region, explaining phenotypic variations of 13.5 to 25.9%, which turned up to be less- or nonsignificant when days to heading and plant height were used as covariates in the analysis, implying a disease escape mechanism. Another major QTL was located on chromosome 5B in CASCABEL, accounting for 8.9 to 21.4% of phenotypic variation. Minor QTL were found on 4A and 4B in CASCABEL; 1B, 4B, and 4D in KATH; and 1B, 2B, and 4B in CIANO T79. Through an analysis of QTL projection onto the IWGSC Chinese Spring reference genome, the 5B QTL in CASCABEL was mapped in the Sb2 region, delimited by the single nucleotide polymorphism marker wsnp_Ku_c50354_55979952 and the simple sequence repeat marker gwm213, with a physical distance of about 14 Mb to the Tsn1 locus.

Genome-Wide Genetic Diversity and Population Structure of Tunisian Durum Wheat Landraces Based on DArTseq Technology.

Tunisia, being part of the secondary center of diversity for durum wheat, has rich unexploited landraces that are being continuously lost and replaced by high yielding modern cultivars. This study aimed to investigate the genetic diversity and population structure of 196 durum wheat lines issued from landraces collected from Tunisia using Diversity Array Technology sequencing (DArTseq) and to understand possible ways of introduction in comparing them to landraces from surrounding countries. A total of 16,148 polymorphic DArTseq markers covering equally the A and B genomes were effective to assess the genetic diversity and to classify the accessions. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish five distinct groups that matched well with the farmer's variety nomenclature. Interestingly, Mahmoudi and Biskri landraces constitute the same gene pool while Jenah Zarzoura constitutes a completely different group. Analysis of molecular variance (AMOVA) showed that the genetic variation was among rather than within the landraces. DAPC analysis of the Tunisian, Mediterranean and West Asian landraces confirmed our previous population structure and showed a genetic similarity between the Tunisian and the North African landraces with the exception of Jenah Zarzoura being the most distant. The genomic characterization of the Tunisian collection will enhance their conservation and sustainable use.

Resistance to Spot Blotch in Two Mapping Populations of Common Wheat Is Controlled by Multiple QTL of Minor Effects.

Spot blotch (SB) is an important fungal disease of wheat in South Asia and South America. Host resistance is regarded as an economical and environmentally friendly approach of controlling SB, and the inheritance of resistance is mostly quantitative. In order to gain a better understanding on the SB resistance mechanism in CIMMYT germplasm, two bi-parental mapping populations were generated, both comprising 232 F2:7 progenies. Elite CIMMYT breeding lines, BARTAI and WUYA, were used as resistant parents, whereas CIANO T79 was used as susceptible parent in both populations. The two populations were evaluated for field SB resistance at CIMMYT's Agua Fria station for three consecutive years, from the 2012⁻2013 to 2014⁻2015 cropping seasons. Phenological traits like plant height (PH) and days to heading (DH) were also determined. Genotyping was performed using the DArTSeq genotyping-by-sequencing (GBS) platform, and a few D-genome specific SNPs and those for phenological traits were integrated for analysis. The most prominent quantitative trait locus (QTL) in both populations was found on chromosome 5AL at the Vrn-A1 locus, explaining phenotypic variations of 7⁻27%. Minor QTL were found on chromosomes 1B, 3A, 3B, 4B, 4D, 5B and 6D in BARTAI and on chromosomes 1B, 2A, 2D and 4B in WUYA, whereas minor QTL contributed by CIANO T79 were identified on chromosome 1B, 1D, 3A, 4B and 7A. In summary, resistance to SB in the two mapping populations was controlled by multiple minor QTL, with strong influence from Vrn-A1.

Genome-wide association analysis of tan spot disease resistance in durum wheat accessions from Tunisia.

Background: Tunisia harbors a rich collection of unexploited durum wheat landraces (Triticum durum ssp. durum) that have been gradually replaced by elite cultivars since the 1970s. These landraces represent an important potential source for broadening the genetic background of elite durum wheat cultivars and for the introgression of novel genes for key traits, including disease resistance, into these cultivars. Methods: In this study, single nucleotide polymorphism (SNP) markers were used to investigate the genetic diversity and population structure of a core collection of 235 durum wheat accessions consisting mainly of landraces. The high phenotypic and genetic diversity of the fungal pathogen Pyrenophora tritici-repentis (cause of tan spot disease of wheat) in Tunisia allowed the assessment of the accessions for tan spot resistance at the adult plant stage under field conditions over three cropping seasons. A genome-wide association study (GWAS) was performed using a 90k SNP array. Results: Bayesian population structure analysis with 9191 polymorphic SNP markers classified the accessions into two groups, where groups 1 and 2 included 49.79% and 31.49% of the accessions, respectively, while the remaining 18.72% were admixtures. Principal coordinate analysis, the unweighted pair group method with arithmetic mean and the neighbor-joining method clustered the accessions into three to five groups. Analysis of molecular variance indicated that 76% of the genetic variation was among individuals and 23% was between individuals. Genome-wide association analyses identified 26 SNPs associated with tan spot resistance and explained between 8.1% to 20.2% of the phenotypic variation. The SNPs were located on chromosomes 1B (1 SNP), 2B (4 SNPs), 3A (2 SNPs), 3B (2 SNPs), 4A (2 SNPs), 4B (1 SNP), 5A (2 SNPs), 5B (4 SNPs), 6A (5 SNPs), 6B (2 SNPs), and 7B (1 SNP). Four markers, one on each of chromosomes 1B, and 5A, and two on 5B, coincided with previously reported SNPs for tan spot resistance, while the remaining SNPs were either novel markers or closely related to previously reported SNPs. Eight durum wheat accessions were identified as possible novel sources of tan spot resistance that could be introgressed into elite cultivars. Conclusion: The results highlighted the significance of chromosomes 2B, 5B, and 6A as genomic regions associated with tan spot resistance.

A reference linkage map for Eucalyptus.

BACKGROUND: Genetic linkage maps are invaluable resources in plant research. They provide a key tool for many genetic applications including: mapping quantitative trait loci (QTL); comparative mapping; identifying unlinked (i.e. independent) DNA markers for fingerprinting, population genetics and phylogenetics; assisting genome sequence assembly; relating physical and recombination distances along the genome and map-based cloning of genes. Eucalypts are the dominant tree species in most Australian ecosystems and of economic importance globally as plantation trees. The genome sequence of E. grandis has recently been released providing unprecedented opportunities for genetic and genomic research in the genus. A robust reference linkage map containing sequence-based molecular markers is needed to capitalise on this resource. Several high density linkage maps have recently been constructed for the main commercial forestry species in the genus (E. grandis, E. urophylla and E. globulus) using sequenced Diversity Arrays Technology (DArT) and microsatellite markers. To provide a single reference linkage map for eucalypts a composite map was produced through the integration of data from seven independent mapping experiments (1950 individuals) using a marker-merging method. RESULTS: The composite map totalled 1107 cM and contained 4101 markers; comprising 3880 DArT, 213 microsatellite and eight candidate genes. Eighty-one DArT markers were mapped to two or more linkage groups, resulting in the 4101 markers being mapped to 4191 map positions. Approximately 13% of DArT markers mapped to identical map positions, thus the composite map contained 3634 unique loci at an average interval of 0.31 cM. CONCLUSION: The composite map represents the most saturated linkage map yet produced in Eucalyptus. As the majority of DArT markers contained on the map have been sequenced, the map provides a direct link to the E. grandis genome sequence and will serve as an important reference for progressing eucalypt research.

Genomic selection for growth and wood quality in Eucalyptus: capturing the missing heritability and accelerating breeding for complex traits in forest trees.

• Genomic selection (GS) is expected to cause a paradigm shift in tree breeding by improving its speed and efficiency. By fitting all the genome-wide markers concurrently, GS can capture most of the 'missing heritability' of complex traits that quantitative trait locus (QTL) and association mapping classically fail to explain. Experimental support of GS is now required. • The effectiveness of GS was assessed in two unrelated Eucalyptus breeding populations with contrasting effective population sizes (N(e) = 11 and 51) genotyped with > 3000 DArT markers. Prediction models were developed for tree circumference and height growth, wood specific gravity and pulp yield using random regression best linear unbiased predictor (BLUP). • Accuracies of GS varied between 0.55 and 0.88, matching the accuracies achieved by conventional phenotypic selection. Substantial proportions (74-97%) of trait heritability were captured by fitting all genome-wide markers simultaneously. Genomic regions explaining trait variation largely coincided between populations, although GS models predicted poorly across populations, likely as a result of variable patterns of linkage disequilibrium, inconsistent allelic effects and genotype × environment interaction. • GS brings a new perspective to the understanding of quantitative trait variation in forest trees and provides a revolutionary tool for applied tree improvement. Nevertheless population-specific predictive models will likely drive the initial applications of GS in forest tree breeding.

Emmer Wheat Eco-Geographic and Genomic Congruence Shapes Phenotypic Performance under Mediterranean Climate.

Emmer wheat (Triticum turgidum ssp. dicoccum) is one of the world's oldest domesticated crops, and it harbors a potentially rich reservoir of agronomic and nutritional quality trait variations. The growing global demand for plant-based health-food niche markets has promoted new commercial interest in ancient grains, including Emmer wheat. Although T. dicoccum can also perform well under harsh environments, its cultivation along the Mediterranean agro-ecosystems is sparse. Here, we analyze a unique tetraploid wheat collection (n = 121) representing a wide geographic range of Emmer accessions, using 9897 DArTseq markers and on-field phenotypic characterization to quantify the extent of diversity among populations and the interactions between eco-geographic, genetic, and phenotypic attributes. Population genomic inferences based on the DArTseq data indicated that the collection could be split into four distinguished clusters in accordance with their eco-geographic origin although significant phenotypic variation was observed within clusters. Superior early vegetative vigor, shorter plant height, and early phenology were observed among emmer wheat accessions from Ethiopia compared to accessions from northern regions. This adaptive advantage highlights the potential of emmer wheat as an exotic germplasm for wheat improvement through breeding. The direct integration of such germplasm into conventional or organic farming agro-systems under the Mediterranean basin climate is also discussed.

Exploring the Genetic Diversity and Population Structure of Wheat Landrace Population Conserved at ICARDA Genebank.

Landraces are considered a valuable source of potential genetic diversity that could be used in the selection process in any plant breeding program. Here, we assembled a population of 600 bread wheat landraces collected from eight different countries, conserved at the ICARDA's genebank, and evaluated the genetic diversity and the population structure of the landraces using single nucleotide polymorphism (SNP) markers. A total of 11,830 high-quality SNPs distributed across the genomes A (40.5%), B (45.9%), and D (13.6%) were used for the final analysis. The population structure analysis was evaluated using the model-based method (STRUCTURE) and distance-based methods [discriminant analysis of principal components (DAPC) and principal component analysis (PCA)]. The STRUCTURE method grouped the landraces into two major clusters, with the landraces from Syria and Turkey forming two clusters with high proportions of admixture, whereas the DAPC and PCA analysis grouped the population into three subpopulations mostly according to the geographical information of the landraces, i.e., Syria, Iran, and Turkey with admixture. The analysis of molecular variance revealed that the majority of the variation was due to genetic differences within the populations as compared with between subpopulations, and it was the same for both the cluster-based and distance-based methods. Genetic distance analysis was also studied to estimate the differences between the landraces from different countries, and it was observed that the maximum genetic distance (0.389) was between the landraces from Spain and Palestine, whereas the minimum genetic distance (0.013) was observed between the landraces from Syria and Turkey. It was concluded from the study that the model-based methods (DAPC and PCA) could dissect the population structure more precisely when compared with the STRUCTURE method. The population structure and genetic diversity analysis of the bread wheat landraces presented here highlight the complex genetic architecture of the landraces native to the Fertile Crescent region. The results of this study provide useful information for the genetic improvement of hexaploid wheat and facilitate the use of landraces in wheat breeding programs.

Worldwide Selection Footprints for Drought and Heat in Bread Wheat (Triticum aestivum L.).

Genome-environment Associations (GEA) or Environmental Genome-Wide Association scans (EnvGWAS) have been poorly applied for studying the genomics of adaptive traits in bread wheat landraces (Triticum aestivum L.). We analyzed 990 landraces and seven climatic variables (mean temperature, maximum temperature, precipitation, precipitation seasonality, heat index of mean temperature, heat index of maximum temperature, and drought index) in GEA using the FarmCPU approach with GAPIT. Historical temperature and precipitation values were obtained as monthly averages from 1970 to 2000. Based on 26,064 high-quality SNP loci, landraces were classified into ten subpopulations exhibiting high genetic differentiation. The GEA identified 59 SNPs and nearly 89 protein-encoding genes involved in the response processes to abiotic stress. Genes related to biosynthesis and signaling are mainly mediated by auxins, abscisic acid (ABA), ethylene (ET), salicylic acid (SA), and jasmonates (JA), which are known to operate together in modulation responses to heat stress and drought in plants. In addition, we identified some proteins associated with the response and tolerance to stress by high temperatures, water deficit, and cell wall functions. The results provide candidate regions for selection aimed to improve drought and heat tolerance in bread wheat and provide insights into the genetic mechanisms involved in adaptation to extreme environments.

Genome-Wide Association Study for Resistance to Tan Spot in Synthetic Hexaploid Wheat.

Synthetic hexaploid wheat (SHW) has shown effective resistance to a diversity of diseases and insects, including tan spot, which is caused by Pyrenophora tritici-repentis, being an important foliar disease that can attack all types of wheat and several grasses. In this study, 443 SHW plants were evaluated for their resistance to tan spot under controlled environmental conditions. Additionally, a genome-wide association study was conducted by genotyping all entries with the DArTSeq technology to identify marker-trait associations for tan spot resistance. Of the 443 SHW plants, 233 showed resistant and 183 moderately resistant reactions, and only 27 were moderately susceptible or susceptible to tan spot. Durum wheat (DW) parents of the SHW showed moderately susceptible to susceptible reactions. A total of 30 significant marker-trait associations were found on chromosomes 1B (4 markers), 1D (1 marker), 2A (1 marker), 2D (2 markers), 3A (4 markers), 3D (3 markers), 4B (1 marker), 5A (4 markers), 6A (6 markers), 6B (1 marker) and 7D (3 markers). Increased resistance in the SHW in comparison to the DW parents, along with the significant association of resistance with the A and B genome, supported the concept of activating epistasis interaction across the three wheat genomes. Candidate genes coding for F-box and cytochrome P450 proteins that play significant roles in biotic stress resistance were identified for the significant markers. The identified resistant SHW lines can be deployed in wheat breeding for tan spot resistance.

Genome-Wide Association Study for Spot Blotch Resistance in Synthetic Hexaploid Wheat.

Spot blotch (SB) caused by Bipolaris sorokiniana (Sacc.) Shoem is a destructive fungal disease affecting wheat and many other crops. Synthetic hexaploid wheat (SHW) offers opportunities to explore new resistance genes for SB for introgression into elite bread wheat. The objectives of our study were to evaluate a collection of 441 SHWs for resistance to SB and to identify potential new genomic regions associated with the disease. The panel exhibited high SB resistance, with 250 accessions showing resistance and 161 showing moderate resistance reactions. A genome-wide association study (GWAS) revealed a total of 41 significant marker-trait associations for resistance to SB, being located on chromosomes 1B, 1D, 2A, 2B, 2D, 3A, 3B, 3D, 4A, 4D, 5A, 5D, 6D, 7A, and 7D; yet none of them exhibited a major phenotypic effect. In addition, a partial least squares regression was conducted to validate the marker-trait associations, and 15 markers were found to be most important for SB resistance in the panel. To our knowledge, this is the first GWAS to investigate SB resistance in SHW that identified markers and resistant SHW lines to be utilized in wheat breeding.

State of ex situ conservation of landrace groups of 25 major crops.

Crop landraces have unique local agroecological and societal functions and offer important genetic resources for plant breeding. Recognition of the value of landrace diversity and concern about its erosion on farms have led to sustained efforts to establish ex situ collections worldwide. The degree to which these efforts have succeeded in conserving landraces has not been comprehensively assessed. Here we modelled the potential distributions of eco-geographically distinguishable groups of landraces of 25 cereal, pulse and starchy root/tuber/fruit crops within their geographic regions of diversity. We then analysed the extent to which these landrace groups are represented in genebank collections, using geographic and ecological coverage metrics as a proxy for genetic diversity. We find that ex situ conservation of landrace groups is currently moderately comprehensive on average, with substantial variation among crops; a mean of 63% ± 12.6% of distributions is currently represented in genebanks. Breadfruit, bananas and plantains, lentils, common beans, chickpeas, barley and bread wheat landrace groups are among the most fully represented, whereas the largest conservation gaps persist for pearl millet, yams, finger millet, groundnut, potatoes and peas. Geographic regions prioritized for further collection of landrace groups for ex situ conservation include South Asia, the Mediterranean and West Asia, Mesoamerica, sub-Saharan Africa, the Andean mountains of South America and Central to East Asia. With further progress to fill these gaps, a high degree of representation of landrace group diversity in genebanks is feasible globally, thus fulfilling international targets for their ex situ conservation.

Population genetic analysis and phylogeny reconstruction in Eucalyptus (Myrtaceae) using high-throughput, genome-wide genotyping.

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species.

Labelling Selective Sweeps Used in Durum Wheat Breeding from a Diverse and Structured Panel of Landraces and Cultivars.

A panel of 387 durum wheat genotypes including Mediterranean landraces and modern cultivars was characterized with 46,161 diversity arrays technology (DArTseq) markers. Analysis of population structure uncovered the existence of five subpopulations (SP) related to the pattern of migration of durum wheat from the domestication area to the west of the Mediterranean basin (SPs 1, 2, and 3) and further improved germplasm (SPs 4 and 5). The total genetic diversity (HT) was 0.40 with a genetic differentiation (GST) of 0.08 and a mean gene flow among SPs of 6.02. The lowest gene flow was detected between SP 1 (presumably the ancient genetic pool of the panel) and SPs 4 and 5. However, gene flow from SP 2 to modern cultivars was much higher. The highest gene flow was detected between SP 3 (western Mediterranean germplasm) and SP 5 (North American and European cultivars). A genome wide association study (GWAS) approach using the top ten eigenvectors as phenotypic data revealed the presence of 89 selective sweeps, represented as quantitative trait loci (QTL) hotspots, widely distributed across the durum wheat genome. A principal component analysis (PCoA) using 147 markers with -log10p > 5 identified three regions located on chromosomes 2A, 2B and 3A as the main drivers for differentiation of Mediterranean landraces. Gene flow between SPs offers clues regarding the putative use of Mediterranean old durum germplasm by the breeding programs represented in the structure analysis. EigenGWAS identified selective sweeps among landraces and modern cultivars. The analysis of the corresponding genomic regions in the 'Zavitan', 'Svevo' and 'Chinese Spring' genomes discovered the presence of important functional genes including Ppd, Vrn, Rht, and gene models involved in important biological processes including LRR-RLK, MADS-box, NAC, and F-box.

Genetic Diversity and Population Structure Analysis of Triticum aestivum L. Landrace Panel from Afghanistan.

Landraces are a potential source of genetic diversity and provide useful genetic resources to cope with the current and future challenges in crop breeding. Afghanistan is located close to the centre of origin of hexaploid wheat. Therefore, understanding the population structure and genetic diversity of Afghan wheat landraces is of enormous importance in breeding programmes for the development of high-yielding cultivars as well as broadening the genetic base of bread wheat. Here, a panel of 363 bread wheat landraces collected from seven north and north-eastern provinces of Afghanistan were evaluated for population structure and genetic diversity using single nucleotide polymorphic markers (SNPs). The genotyping-by-sequencing of studied landraces after quality control provided 4897 high-quality SNPs distributed across the genomes A (33.75%), B (38.73%), and D (27.50%). The population structure analysis was carried out by two methods using model-based STRUCTURE analysis and cluster-based discriminant analysis of principal components (DAPC). The analysis of molecular variance showed a higher proportion of variation within the sub-populations compared with the variation observed as a whole between sub-populations. STRUCTURE and DAPC analysis grouped the majority of the landraces from Badakhshan and Takhar together in one cluster and the landraces from Baghlan and Kunduz in a second cluster, which is in accordance with the micro-climatic conditions prevalent within the north-eastern agro-ecological zone. Genetic distance analysis was also studied to identify differences among the Afghan regions; the strongest correlation was observed for the Badakhshan and Takhar (0.003), whereas Samangan and Konarha (0.399) showed the highest genetic distance. The population structure and genetic diversity analysis highlighted the complex genetic variation present in the landraces which were highly correlated to the geographic origin and micro-climatic conditions within the agro-climatic zones of the landraces. The higher proportions of admixture could be attributed to historical unsupervised exchanges of seeds between the farmers of the central and north-eastern provinces of Afghanistan. The results of this study will provide useful information for genetic improvement in wheat and is essential for association mapping and genomic prediction studies to identify novel sources for resistance to abiotic and biotic stresses.

Comparison of Genomic Prediction Methods for Yellow, Stem, and Leaf Rust Resistance in Wheat Landraces from Afghanistan.

Wheat rust diseases, including yellow rust (Yr; also known as stripe rust) caused by Puccinia striiformis Westend. f. sp. tritici, leaf rust (Lr) caused by Puccinia triticina Eriks. and stem rust (Sr) caused by Puccinia graminis Pres f. sp. tritici are major threats to wheat production all around the globe. Durable resistance to wheat rust diseases can be achieved through genomic-assisted prediction of resistant accessions to increase genetic gain per unit time. Genomic prediction (GP) is a promising technology that uses genomic markers to estimate genomic-assisted breeding values (GBEVs) for selecting resistant plant genotypes and accumulating favorable alleles for adult plant resistance (APR) to wheat rust diseases. To evaluate GP we compared the predictive ability of nine different parametric, semi-parametric and Bayesian models including Genomic Unbiased Linear Prediction (GBLUP), Ridge Regression (RR), Least Absolute Shrinkage and Selection Operator (LASSO), Elastic Net (EN), Bayesian Ridge Regression (BRR), Bayesian A (BA), Bayesian B (BB), Bayesian C (BC) and Reproducing Kernel Hilbert Spacing model (RKHS) to estimate GEBV's for APR to yellow, leaf and stem rust of wheat in a panel of 363 bread wheat landraces of Afghanistan origin. Based on five-fold cross validation the mean predictive abilities were 0.33, 0.30, 0.38, and 0.33 for Yr (2016), Yr (2017), Lr, and Sr, respectively. No single model outperformed the rest of the models for all traits. LASSO and EN showed the lowest predictive ability in four of the five traits. GBLUP and RR gave similar predictive abilities, whereas Bayesian models were not significantly different from each other as well. We also investigated the effect of the number of genotypes and the markers used in the analysis on the predictive ability of the GP model. The predictive ability was highest with 1000 markers and there was a linear trend in the predictive ability and the size of the training population. The results of the study are encouraging, confirming the feasibility of GP to be effectively applied in breeding programs for resistance to all three wheat rust diseases.