Mining Centuries Old In situ Conserved Turkish Wheat Landraces for Grain Yield and Stripe Rust Resistance Genes

Genomic wide association study and selective sweep analysis identify genes associated with improved yield under drought in Turkish winter wheat germplasm

A panel comprising of 84 Turkish winter wheat landraces (LR) and 73 modern varieties (MV) was analyzed with genome wide association study (GWAS) to identify genes/genomic regions associated with increased yield under favorable and drought conditions. In addition, selective sweep analysis was conducted to detect signatures of selection in the winter wheat genome driving the differentiation between LR and MV, to gather an understanding of genomic regions linked to adaptation and yield improvement. The panel was genotyped with 25 K wheat SNP array and phenotyped for agronomic traits for two growing seasons (2018 and 2019) in Konya, Turkey. Year 2018 was treated as drought environment due to very low precipitation prior to heading whereas year 2019 was considered as a favorable season. GWAS conducted with SNPs and haplotype blocks using mixed linear model identified 18 genomic regions in the vicinities of known genes i.e., TaERF3-3A, TaERF3-3B, DEP1-5A, FRIZZY PANICLE-2D, TaSnRK23-1A, TaAGL6-A, TaARF12-2A, TaARF12-2B, WAPO1, TaSPL16-7D, TaTGW6-A1, KAT-2B, TaOGT1, TaSPL21-6B, TaSBEIb, trs1/WFZP-A, TaCwi-A1-2A and TaPIN1-7A associated with grain yield (GY) and yield related traits. Haplotype-based GWAS identified five haplotype blocks (H1A-42, H2A-71, H4A-48, H7B-123 and H7B-124), with the favorable haplotypes showing a yield increase of > 700 kg/ha in the drought season. SNP-based GWAS, detected only one larger effect genomic region on chromosome 7B, in common with haplotype-based GWAS. On an average, the percentage variation (PV) explained by haplotypes was 8.0% higher than PV explained by SNPs for all the investigated traits. Selective sweep analysis detected 39 signatures of selection between LR and MV of which 15 were within proximity of known functional genes controlling flowering (PRR-A1, PPR-D1, TaHd1-6B), GY and GY components (TaSus2-2B, TaGS2-B1, AG1-1A/WAG1-1A, DUO-A1, DUO-B1, AG2-3A/WAG2-3A, TaLAX1, TaSnRK210-4A, FBP, TaLAX1, TaPIL1 and AP3-1-7A/WPA3-7A) and 10 regions underlying various transcription factors and regulatory genes. The study outcomes contribute to utilization of LR in breeding winter wheat.

Genome-wide association mapping in elite winter wheat breeding for yield improvement

Increased grain yield (GY) is the primary breeding target of wheat breeders. We performed the genome-wide association study (GWAS) on 168 elite winter wheat lines from an ongoing breeding program to identify the main determinants of grain yield. Sequencing of Diversity Array Technology fragments (DArTseq) resulted in 19,350 single-nucleotide polymorphism (SNP) and presence-absence variation (PAV) markers. We identified 15 main genomic regions located in ten wheat chromosomes (1B, 2B, 2D, 3A, 3D, 5A, 5B, 6A, 6B, and 7B) that explained from 7.9 to 20.3% of the variation in grain yield and 13.3% of the yield stability. Loci identified in the reduced genepool are important for wheat improvement using marker-assisted selection. We found marker-trait associations between three genes involved in starch biosynthesis and grain yield. Two starch synthase genes (TraesCS2B03G1238800 and TraesCS2D03G1048800) and a sucrose synthase gene (TraesCS3D03G0024300) were found in regions of QGy.rut-2B.2, QGy.rut-2D.1, and QGy.rut-3D, respectively. These loci and other significantly associated SNP markers found in this study can be used for pyramiding favorable alleles in high-yielding varieties or to improve the accuracy of prediction in genomic selection.

Genome-wide association study of resistance to PstS2 and Warrior races of Puccinia striiformis f. sp. tritici (stripe rust) in bread wheat landraces

Stripe or yellow rust, caused by Puccinia striiformis Westend. f. sp. tritici is a major threat to bread wheat production worldwide. The breakdown in resistance of certain major genes and newly emerging aggressive races of stripe rusts pose serious concerns in all main wheat growing areas of the world. To identify new sources of resistance and associated QTL for effective utilization in future breeding programs an association mapping (AM) panel comprising of 600 bread wheat landraces collected from eight different countries conserved at ICARDA gene bank were evaluated for seedling and adult plant resistance against the PstS2 and Warrior races of stripe rust at the Regional Cereal Rust Research Center (RCRRC), Izmir, Turkey during 2016, 2018 and 2019. A set of 25,169 informative SNP markers covering the whole genome were used to examine the population structure, linkage disequilibrium and marker-trait associations in the AM panel. The genome-wide association study (GWAS) was carried out using a Mixed Linear Model (MLM). We identified 47 SNP markers across 19 chromosomes with significant SNP-trait associations for both seedling stage and adult plant resistance. The threshold of significance for all SNP-trait associations was determined by the false discovery rate (q) ≤ 0.05. Three genomic regions (QYr.1D_APR, QYr.3A_seedling and QYr.7D_seedling) identified in this study do not correspond to previously reported Yr genes or QTL, suggesting new genomic regions for stripe rust resistance.

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.

Genome-wide association analysis of Mexican bread wheat landraces for resistance to yellow and stem rust

Deploying under-utilized landraces in wheat breeding has been advocated to accelerate genetic gains in current era of genomics assisted breeding. Mexican bread wheat landraces (Creole wheats) represent an important resource for the discovery of novel alleles including disease resistance. A core set of 1,098 Mexican landraces was subjected to multi-location testing for rust diseases in India, Mexico and Kenya. The landrace core set showed a continuous variation for yellow (YR) and stem rust (SR) disease severity. Principal component analysis differentiated Mexican landraces into three groups based on their respective collection sites. Linkage disequilibrium (LD) decay varied from 10 to 32 Mb across chromosomes with an averge of 23Mb across whole genome. Genome-wide association analysis revealed marker-trait associations for YR resistance in India and Mexico as well as for SR resistance in Kenya. In addition, significant additive-additive interaction effects were observed for both YR and SR resistance including genomic regions on chromosomes 1BL and 3BS, which co-locate with pleiotropic genes Yr29/Lr46/Sr58/Pm39/Ltn2 and Sr2/Yr30/Lr27, respectively. Study reports novel genomic associations for YR (chromosomes 1AL, 2BS, and 3BL) and SR (chromosomes 2AL, 4DS, and 5DS). The novel findings in Creole wheat landraces can be efficiently utilized for the wheat genetic improvement.

Genome-Wide Association Study of Tan Spot Resistance in a Hexaploid Wheat Collection From Kazakhstan

Tan spot, caused by Pyrenophora tritici-repentis, is a serious foliar disease of wheat in Kazakhstan with reported yield losses as high as 50% during epidemic years. Here, we report the evaluation of a collection of 191 hexaploid spring and winter wheat lines for tan spot resistance and its underlying genetic architecture using genome-wide association study (GWAS). Our wheat collection comprised candidate varieties from Kazakhstan, Russia, and CIMMYT. It was genotyped using the DArTseq technology and phenotyped for resistance to tan spot at seedling and adult plant stages in Kazakhstan. DArTseq SNPs revealed high genetic diversity (average polymorphic information content = 0.33) in the panel and genome-wide linkage disequilibrium decay at 22 Mb (threshold r² = 0.1). Principal component analysis revealed a clear separation of Eurasian germplasm from CIMMYT and IWWIP lines. GWAS identified 34 marker-trait associations (MTA) for resistance to tan spot and the amount of phenotypic variation explained by these MTA ranged from 4% to 13.7%. Our results suggest the existence of novel valuable resistant alleles on chromosomes 3BS, and 5DL and 6AL for resistance to Race 1 and Race 5, respectively, in addition to known genes tsn1 and tsc2. On chromosome 6AL, a genomic region spanning 3 Mb was identified conferring resistance to both Race 1 and Race 5. Epistatic interaction of associated loci was revealed on chromosomes 1B, 5B, 7B, 5A, and 6A contributing to additional variation of 3.2–11.7%. Twenty-five lines with the best allele combinations of SNPs associated with resistance to both races have been identified as candidates for future variety release and breeding. The results of the present study will be further validated in other independent genetic backgrounds to be able to use markers in breeding.

Genome-wide association mapping reveals potential novel loci controlling stripe rust resistance in a Chinese wheat landrace diversity panel from the southern autumn-sown spring wheat zone

BACKGROUND

Stripe rust, caused by the fungal pathogen Puccinia striiformis f. sp. tritici (Pst), is a serious foliar disease of wheat. Identification of novel stripe rust resistance genes and cultivation of resistant cultivars are considered to be the most effective approaches to control this disease. In this study, we evaluated the infection type (IT), disease severity (DS) and area under the disease progress curve (AUDPC) of 143 Chinese wheat landrace accessions for stripe rust resistance. Assessments were undertaken in five environments at the adult-plant stage with Pst mixture races under field conditions. In addition, IT was assessed at the seedling stage with two prevalent Pst races (CYR32 and CYR34) under a controlled greenhouse environment.

RESULTS

Seventeen accessions showed stable high-level resistance to stripe rust across all environments in the field tests. Four accessions showed resistance to the Pst races CYR32 and CYR34 at the seedling stage. Combining phenotypic data from the field and greenhouse trials with 6404 markers that covered the entire genome, we detected 17 quantitative trait loci (QTL) on 11 chromosomes for IT associated with seedling resistance and 15 QTL on seven chromosomes for IT, final disease severity (FDS) or AUDPC associated with adult-plant resistance. Four stable QTL detected on four chromosomes, which explained 9.99-23.30% of the phenotypic variation, were simultaneously associated with seedling and adult-plant resistance. Integrating a linkage map of stripe rust resistance in wheat, 27 QTL overlapped with previously reported genes or QTL, whereas four and one QTL conferring seedling and adult-plant resistance, respectively, were mapped distantly from previously reported stripe rust resistance genes or QTL and thus may be novel resistance loci.

CONCLUSIONS

Our results provided an integrated overview of stripe rust resistance resources in a wheat landrace diversity panel from the southern autumn-sown spring wheat zone of China. The identified resistant accessions and resistance loci will be useful in the ongoing effort to develop new wheat cultivars with strong resistance to stripe rust.

Haplotype-Based, Genome-Wide Association Study Reveals Stable Genomic Regions for Grain Yield in CIMMYT Spring Bread Wheat

We untangled key regions of the genetic architecture of grain yield (GY) in CIMMYT spring bread wheat by conducting a haplotype-based, genome-wide association study (GWAS), together with an investigation of epistatic interactions using seven large sets of elite yield trials (EYTs) consisting of a total of 6,461 advanced breeding lines. These lines were phenotyped under irrigated and stress environments in seven growing seasons (2011-2018) and genotyped with genotyping-by-sequencing markers. Genome-wide 519 haplotype blocks were constructed, using a linkage disequilibrium-based approach covering 14,036 Mb in the wheat genome. Haplotype-based GWAS identified 7, 4, 10, and 15 stable (significant in three or more EYTs) associations in irrigated (I), mild drought (MD), severe drought (SD), and heat stress (HS) testing environments, respectively. Considering all EYTs and the four testing environments together, 30 stable associations were deciphered with seven hotspots identified on chromosomes 1A, 1B, 2B, 4A, 5B, 6B, and 7B, where multiple haplotype blocks were associated with GY. Epistatic interactions contributed significantly to the genetic architecture of GY, explaining variation of 3.5-21.1%, 3.7-14.7%, 3.5-20.6%, and 4.4- 23.1% in I, MD, SD, and HS environments, respectively. Our results revealed the intricate genetic architecture of GY, controlled by both main and epistatic effects. The importance of these results for practical applications in the CIMMYT breeding program is discussed.

A novel approach to the characterization of old wheat (Triticum aestivum L.) varieties by complex rheological analysis

BACKGROUND

Lines of the internationally recognized old Hungarian Bánkúti 1201 variety are important genetic resources for breeding programmes. Their protein composition and gluten dependent technological traits have been comprehensively studied, however, little information is available about their carbohydrate dependent viscous properties. The aim of this work was to obtain comprehensive rheological characterization of all sublines of Bánkúti 1201 maintained at Martonvásár and to investigate their variability if the carbohydrate dependent viscous behaviour was also included in the analyses.

RESULTS

The majority of the lines reflected the famously good mixing quality of Bánkúti, however, much higher diversity of pasting behaviour was detected. Cluster analysis of the Mixolab data was performed resulting in four sample groups. Since several lines of similar mixing properties had significantly different pasting characteristics, it was assumed that classification was mainly based on the viscous properties. From each cluster two to three representative samples were selected for wider examination using conventional testing methods. These results also supported the higher variability of pasting behaviour of the lines, which can be critical for end product quality. The members of the second cluster can be highlighted due to their waxy wheat like behaviour.

CONCLUSIONS

Possible reasons for the great variability of pasting behaviour could be the compositional and structural differences of starch and other carbohydrates (e.g. arabinoxylans). Complex rheological characterization and study of molecular background can provide information about important traits from the point of view of technology and product development, which are unknown in the case of old wheat varieties and landraces. © 2020 Society of Chemical Industry.

Incorporating Genome-Wide Association Mapping Results Into Genomic Prediction Models for Grain Yield and Yield Stability in CIMMYT Spring Bread Wheat

Untangling the genetic architecture of grain yield (GY) and yield stability is an important determining factor to optimize genomics-assisted selection strategies in wheat. We conducted in-depth investigation on the above using a large set of advanced bread wheat lines (4,302), which were genotyped with genotyping-by-sequencing markers and phenotyped under contrasting (irrigated and stress) environments. Haplotypes-based genome-wide-association study (GWAS) identified 58 associations with GY and 15 with superiority index Pi (measure of stability). Sixteen associations with GY were "environment-specific" with two on chromosomes 3B and 6B with the large effects and 8 associations were consistent across environments and trials. For Pi, 8 associations were from chromosomes 4B and 7B, indicating 'hot spot' regions for stability. Epistatic interactions contributed to an additional 5-9% variation on average. We further explored whether integrating consistent and robust associations identified in GWAS as fixed effects in prediction models improves prediction accuracy. For GY, the model accounting for the haplotype-based GWAS loci as fixed effects led to up to 9-10% increase in prediction accuracy, whereas for Pi this approach did not provide any advantage. This is the first report of integrating genetic architecture of GY and yield stability into prediction models in wheat.

Haplotypes-based genetic analysis: benefits and challenges

The increasing availability of Single Nucleotide Polymorphisms (SNPs) discovered by Next Generation Sequencing will enable a range of new genetic analyses in crops, which was not possible before. Concomitantly, researchers will face the challenge of handling large data sets at the whole-genome level. By grouping thousands of SNPs into a few hundred haplotype blocks, complexity of the data can be reduced with fewer statistical tests and a lower probability of spurious associations. Owing to the strong genome structure present in breeding lines of most crops, the deployment of haplotypes could be a powerful complement to improve efficiency of marker-assisted and genomic selection. This review describes in brief the commonly used approaches to construct haplotype blocks and some examples in animals and crops are cited where haplotype-based dissection of traits were proven beneficial. Some important considerations and facts while working with haplotypes in crops are reviewed at the end.

Validation of Candidate Gene-Based Markers and Identification of Novel Loci for Thousand-Grain Weight in Spring Bread Wheat

Increased thousand-grain weight (TGW) is an important breeding target for indirectly improving grain yield (GY). Fourteen reported candidate genes known to enhance TGW were evaluated in two independent and existing datasets of wheat at CIMMYT, the Elite Yield Trial (EYT) from 2015 to 2016 (EYT2015-16) and the Wheat Association Mapping Initiative (WAMI) panel, to study their allele effects on TGW and to verify their suitability for marker-assisted selection. Of these, significant associations were detected for only one gene (TaGs3-D1) in the EYT2015-16 and two genes (TaTGW6 and TaSus1) in WAMI. The reported favorable alleles of TaGs3-D1 and TaTGW6 genes decreased TGW in the datasets. A haplotype-based genome wide association study was implemented to identify the genetic determinants of TGW on a large set of CIMMYT germplasm (4,302 lines comprising five EYTs), which identified 15 haplotype blocks to be significantly associated with TGW. Four of them, identified on chromosomes 4A, 6A, and 7A, were associated with TGW in at least three EYTs. The locus on chromosome 6A (Hap-6A-13) had the largest effect on TGW and additionally GY with increases of up to 2.60 g and 258 kg/ha, respectively. Discovery of novel TGW loci described in our study expands the opportunities for developing diagnostic markers and for multi-gene pyramiding to derive new allele combinations for enhanced TGW and GY in CIMMYT wheat.

Towards the Genomic Basis of Local Adaptation in Landraces

Landraces are key elements of agricultural biodiversity that have long been considered a source of useful traits. Their importance goes beyond subsistence agriculture and the essential need to preserve genetic diversity, because landraces are farmer-developed populations that are often adapted to environmental conditions of significance to tackle environmental concerns. It is therefore increasingly important to identify adaptive traits in crop landraces and understand their molecular basis. This knowledge is potentially useful for promoting more sustainable agricultural techniques, reducing the environmental impact of high-input cropping systems, and diminishing the vulnerability of agriculture to global climate change. In this review, we present an overview of the opportunities and limitations offered by landraces’ genomics. We discuss how rapid advances in DNA sequencing techniques, plant phenotyping, and recombinant DNA-based biotechnology encourage both the identification and the validation of the genomic signature of local adaptation in crop landraces. The integration of ‘omics’ sciences, molecular population genetics, and field studies can provide information inaccessible with earlier technological tools. Although empirical knowledge on the genetic and genomic basis of local adaptation is still fragmented, it is predicted that genomic scans for adaptation will unlock an intraspecific molecular diversity that may be different from that of modern varieties.

Joining smallholder farmers' traditional knowledge with metric traits to select better varieties of Ethiopian wheat

Smallholder farming communities face highly variable climatic conditions that threaten locally adapted, low-input agriculture. The benefits of modern crop breeding may fail to reach their fields when broadly adapted genetic materials do not address local requirements. To date, participatory methods only scratched the surface of the exploitability of farmers’ traditional knowledge in breeding. In this study, 30 smallholder farmers in each of two locations in Ethiopia provided quantitative evaluations of earliness, spike morphology, tillering capacity and overall quality on 400 wheat genotypes, mostly traditional varieties, yielding altogether 192,000 data points. Metric measurements of ten agronomic traits were simultaneously collected, allowing to systematically break down farmers’ preferences on quantitative phenotypes. Results showed that the relative importance of wheat traits differed by gender and location. Farmer traits were variously contributed by metric traits, and could only partially be explained by them. Eventually, farmer trait values were used to produce a ranking of the 400 wheat varieties identifying the trait combinations most desired by farmers. The study scale and methods lead to a better understanding of the quantitative basis of Ethiopian smallholder farmer preference in wheat, broadening the discussion for the future of local, sustainable breeding efforts accommodating farmers’ knowledge.