Genetic Diversity, Linkage Disequilibrium and Population Structure of Bulgarian Bread Wheat Assessed by Genome-Wide Distributed SNP Markers: From Old Germplasm to Semi-Dwarf Cultivars

Assessment of genetic structure and trait associations of Watkins wheat landraces under Egyptian field conditions

Background

Wheat landraces represent a reservoir of genetic diversity that can support wheat improvement through breeding. A core panel of 300 Watkins wheat landraces, as well as 16 non-Watkins landraces and elite wheat cultivars, was grown during the 2020-2021 and 2021-2022 seasons at four Agricultural Research Stations in Egypt, Gemmiza, Nubaria, Sakha, and Sids, to evaluate the core panel for agromorphological and yield-related traits. The genetic population structure within these genotypes were assessed using 35,143 single nucleotide polymorphisms (SNPs).

Results

Cluster analyses using Discriminant Analysis of Principal Components (DAPC) and k-means revealed three clusters with moderate genetic differentiation and population structure, possibly due to wheat breeding systems and geographical isolation. The best ancestry was k = 4, but k = 2 and k = 3 were also significant. A genome-wide association study (GWAS) identified clustered marker trait associations (MTAs) linked to thousand kernel weight on chromosome 5A, plant height on chromosomes 3B and 1D, days to heading on chromosomes 2A, 4B, 5B and 1D, and plant maturity on chromosomes 3A, 2B, and 6B. In the future, these MTAs can be used to accelerate the incorporation of beneficial alleles into locally adapted germplasm through marker-assisted selection. Gene enrichment analysis identified key genes within these loci, including Reduced height-1 (Rht-A1) and stress-related genes.

Conclusion

These findings underscore significant genetic connections and the involvement of crucial biological pathways.

Exploring Novel Genomic Loci and Candidate Genes Associated with Plant Height in Bulgarian Bread Wheat via Multi-Model GWAS

In the context of crop breeding, plant height (PH) plays a pivotal role in determining straw and grain yield. Although extensive research has explored the genetic control of PH in wheat, there remains an opportunity for further advancements by integrating genomics with growth-related phenomics. Our study utilizes the latest genome-wide association scan (GWAS) techniques to unravel the genetic basis of temporal variation in PH across 179 Bulgarian bread wheat accessions, including landraces, tall historical, and semi-dwarf modern varieties. A GWAS was performed with phenotypic data from three growing seasons, the calculated best linear unbiased estimators, and the leveraging genotypic information from the 25K Infinium iSelect array, using three statistical methods (MLM, FarmCPU, and BLINK). Twenty-five quantitative trait loci (QTL) associated with PH were identified across fourteen chromosomes, encompassing 21 environmentally stable quantitative trait nucleotides (QTNs), and four haplotype blocks. Certain loci (17) on chromosomes 1A, 1B, 1D, 2A, 2D, 3A, 3B, 4A, 5B, 5D, and 6A remain unlinked to any known Rht (Reduced height) genes, QTL, or GWAS loci associated with PH, and represent novel regions of potential breeding significance. Notably, these loci exhibit varying effects on PH, contribute significantly to natural variance, and are expressed during seedling to reproductive stages. The haplotype block on chromosome 6A contains five QTN loci associated with reduced height and two loci promoting height. This configuration suggests a substantial impact on natural variation and holds promise for accurate marker-assisted selection. The potentially novel genomic regions harbor putative candidate gene coding for glutamine synthetase, gibberellin 2-oxidase, auxin response factor, ethylene-responsive transcription factor, and nitric oxide synthase; cell cycle-related genes, encoding cyclin, regulator of chromosome condensation (RCC1) protein, katanin p60 ATPase-containing subunit, and expansins; genes implicated in stem mechanical strength and defense mechanisms, as well as gene regulators such as transcription factors and protein kinases. These findings enrich the pool of semi-dwarfing gene resources, providing the potential to further optimize PH, improve lodging resistance, and achieve higher grain yields in bread wheat.

Effectiveness and Genetic Control of Trichoderma spp. as a Biological Control of Wheat Powdery Mildew Disease

Wheat powdery mildew (WPM) is one of the most devasting diseases that affects wheat yield worldwide. Few efforts have been made to control such a serious disease. An effective way to control WPM is urgently needed. Biological control is an effective way to control plant diseases worldwide. In this study, the efficiency of three different Trichoderma spp. in controlling WPM at the seedling growth stage was tested using 35 highly diverse wheat genotypes. Highly significant differences were found in WPM resistance among the four treatments, confirming the efficiency of Trichoderma in controlling WPM. Of the three species, T. asperellum T34 (T34) was the most effective species in controlling WPM, as it reduced the symptoms by 50.56%. A set of 196 wheat genotypes was used to identify the genetic control of the WPM resistance induced by T34. A total of 39, 27, and 18 gene models were identified to contain the significant markers under Pm, T34, and the improvement in powdery mildew resistance due to T34 (T34_improvement) conditions. Furthermore, no gene model was common between T34 and Pm, suggesting the presence of completely different genetic systems controlling the resistance under T34 and Pm. The functional annotation and biological process pathways of the detected gene models confirm their association with the normal and induced resistance. This study, for the first time, confirms the efficiency of T34 in controlling WPM and provides a deep understanding of the genetic control of induced and normal resistance to WPM.

Identifying the physiological traits associated with DNA marker using genome wide association in wheat under heat stress

Heat stress poses a significant environmental challenge that profoundly impacts wheat productivity. It disrupts vital physiological processes such as photosynthesis, by impeding the functionality of the photosynthetic apparatus and compromising plasma membrane stability, thereby detrimentally affecting grain development in wheat. The scarcity of identified marker trait associations pertinent to thermotolerance presents a formidable obstacle in the development of marker-assisted selection strategies against heat stress. To address this, wheat accessions were systematically exposed to both normal and heat stress conditions and phenotypic data were collected on physiological traits including proline content, canopy temperature depression, cell membrane injury, photosynthetic rate, transpiration rate (at vegetative and reproductive stage and 'stay-green'. Principal component analysis elucidated the most significant contributors being proline content, transpiration rate, and canopy temperature depression, which exhibited a synergistic relationship with grain yield. Remarkably, cluster analysis delineated the wheat accessions into four discrete groups based on physiological attributes. Moreover, to explore the relationship between physiological traits and DNA markers, 158 wheat accessions were genotyped with 186 SSRs. Allelic frequency and polymorphic information content value were found to be highest on genome A (4.94 and 0.688), chromosome 1A (5.00 and 0.712), and marker Xgwm44 (13.0 and 0.916). Population structure, principal coordinate analysis and cluster analysis also partitioned the wheat accessions into four subpopulations based on genotypic data, highlighting their genetic homogeneity. Population diversity and presence of linkage disequilibrium established the suitability of population for association mapping. Additionally, linkage disequilibrium decay was most pronounced within a 15-20 cM region on chromosome 1A. Association mapping revealed highly significant marker trait associations at Bonferroni correction P < 0.00027. Markers Xwmc418 (located on chromosome 3D) and Xgwm233 (chromosome 7A) demonstrated associations with transpiration rate, while marker Xgwm494 (chromosome 3A) exhibited an association with photosynthetic rates at both vegetative and reproductive stages under heat stress conditions. Additionally, markers Xwmc201 (chromosome 6A) and Xcfa2129 (chromosome 1A) displayed robust associations with canopy temperature depression, while markers Xbarc163 (chromosome 4B) and Xbarc49 (chromosome 5A) were strongly associated with cell membrane injury at both stages. Notably, marker Xbarc49 (chromosome 5A) exhibited a significant association with the 'stay-green' trait under heat stress conditions. These results offers the potential utility in marker-assisted selection, gene pyramiding and genomic selection models to predict performance of wheat accession under heat stress conditions.

Genome-wide scanning to identify and validate single nucleotide polymorphism markers associated with drought tolerance in spring wheat seedlings

Unlike other growth stages of wheat, very few studies on drought tolerance have been done at the seedling stage, and this is due to the complexity and sensitivity of this stage to drought stress resulting from climate change. As a result, the drought tolerance of wheat seedlings is poorly understood and very few genes associated with drought tolerance at this stage were identified. To address this challenge, a set of 172 spring wheat genotypes representing 20 different countries was evaluated under drought stress at the seedling stage. Drought stress was applied on all tested genotypes by water withholding for 13 days. Two types of traits, namely morphological and physiological traits were scored on the leaves of all tested genotypes. Genome-wide association study (GWAS) is one of the effective genetic analysis methods that was used to identify target single nucleotide polymorphism (SNP) markers and candidate genes for later use in marker-assisted selection. The tested plant materials were genotyped using 25k Infinium iSelect array (25K) (herein after it will be identified as 25K) (for 172 genotypes) and genotyping-by-sequencing (GBS) (for 103 genotypes), respectively. The results of genotyping revealed 21,093 25K and 11,362 GBS-SNPs, which were used to perform GWAS analysis for all scored traits. The results of GWAS revealed that 131 and 55 significant SNPs were controlling morphological and physiological traits, respectively. Moreover, a total of eight and seven SNP markers were found to be associated with more than one morphological and physiological trait under drought stress, respectively. Remarkably, 10 significant SNPs found in this study were previously reported for their association with drought tolerance in wheat. Out of the 10 validated SNP markers, four SNPs were associated with drought at the seedling stage, while the remaining six SNPs were associated with drought stress at the reproductive stage. Moreover, the results of gene enrichment revealed 18 and six pathways as highly significant biological and molecular pathways, respectively. The selection based on drought-tolerant alleles revealed 15 genotypes with the highest number of different drought-tolerant alleles. These genotypes can be used as candidate parents in future breeding programs to produce highly drought-tolerant genotypes with high genetic diversity. Our findings in this study provide novel markers and useful information on the genetic basis of drought tolerance at early growth stages.

GWAS in a Collection of Bulgarian Old and Modern Bread Wheat Accessions Uncovers Novel Genomic Loci for Grain Protein Content and Thousand Kernel Weight

Genetic enhancement of grain production and quality is a priority in wheat breeding projects. In this study, we assessed two key agronomic traits-grain protein content (GPC) and thousand kernel weight (TKW)-across 179 Bulgarian contemporary and historic varieties and landraces across three growing seasons. Significant phenotypic variation existed for both traits among genotypes and seasons, and no discernible difference was evident between the old and modern accessions. To understand the genetic basis of the traits, we conducted a genome-wide association study with MLM using phenotypic data from the crop seasons, best linear unbiased estimators, and genotypic data from the 25K Infinium iSelect array. As a result, we detected 16 quantitative trait nucleotides (QTNs) associated with GPC and 15 associated with TKW, all of which passed the false discovery rate threshold. Seven loci favorably influenced GPC, resulting in an increase of 1.4% to 8.1%, while four loci had a positive impact on TKW with increases ranging from 1.9% to 8.4%. While some loci confirmed previously published associations, four QTNs linked to GPC on chromosomes 2A, 7A, and 7B, as well as two QTNs related to TKW on chromosomes 1B and 6A, may represent novel associations. Annotations for proteins involved in the senescence-associated nutrient remobilization and in the following buildup of resources required for seed germination have been found for selected putative candidate genes. These include genes coding for storage proteins, cysteine proteases, cellulose-synthase, alpha-amylase, transcriptional regulators, and F-box and RWP-RK family proteins. Our findings highlight promising genomic regions for targeted breeding programs aimed at improving grain yield and protein content.

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.

Novel genomic regions on chromosome 5B controlling wheat powdery mildew seedling resistance under Egyptian conditions

Wheat powdery mildew (PM) causes significant yield losses worldwide. None of the Egyptian wheat cultivars was detected to be highly resistant to such a severe disease. Therefore, a diverse spring wheat panel was evaluated for PM seedling resistance using different Bgt conidiospores collected from Egyptian fields in two growing seasons. The evaluation was done in two separate experiments. Highly significant differences were found between the two experiments suggesting the presence of different isolates populations. Highly significant differences were found among the tested genotypes confirming the ability to improve PM resistance using the recent panel. Genome-wide association study (GWAS) was done for each experiment separately and a total of 71 significant markers located within 36 gene models were identified. The majority of these markers are located on chromosome 5B. Haplotype block analysis identified seven blocks containing the significant markers on chromosome 5B. Five gene models were identified on the short arm of the chromosome. Gene enrichment analysis identified five and seven pathways based on the biological process and molecular functions respectively for the detected gene models. All these pathways are associated with disease resistance in wheat. The genomic regions on 5B seem to be novel regions that are associated with PM resistance under Egyptian conditions. Selection of superior genotypes was done and Grecian genotypes seem to be a good source for improving PM resistance under Egyptian conditions.

Identification of genomic regions associated with cereal cyst nematode (Heterodera avenae Woll.) resistance in spring and winter wheat

Cereal cyst nematode (CCN) is a major threat to cereal crop production globally including wheat (Triticum aestivum L.). In the present study, single-locus and multi-locus models of Genome-Wide Association Study (GWAS) were used to find marker trait associations (MTAs) against CCN (Heterodera avenae) in wheat. In total, 180 wheat accessions (100 spring and 80 winter types) were screened against H. avenae in two independent years (2018/2019 "Environment 1" and 2019/2020 "Environment 2") under controlled conditions. A set of 12,908 SNP markers were used to perform the GWAS. Altogether, 11 significant MTAs, with threshold value of -log10 (p-values) ≥ 3.0, were detected using 180 wheat accessions under combined environment (CE). A novel MTA (wsnp_Ex_c53387_56641291) was detected under all environments (E1, E2 and CE) and considered to be stable MTA. Among the identified 11 MTAs, eight were novel and three were co-localized with previously known genes/QTLs/MTAs. In total, 13 putative candidate genes showing differential expression in roots, and known to be involved in plant defense mechanisms were reported. These MTAs could help us to identify resistance alleles from new sources, which could be used to identify wheat varieties with enhanced CCN resistance.

Utilizing Genomics to Characterize the Common Oat Gene Pool—The Story of More than a Century of Polish Breeding

This study was undertaken to investigate the diversity and population structure of 487 oat accessions, including breeding lines from the ongoing programs of the three largest Polish breeding companies, along with modern and historical Polish and foreign cultivars. The analysis was based on 7411 DArTseq-derived SNPs distributed among three sub-genomes (A, C, and D). The heterogeneity of the studied material was very low, as only cultivars and advanced breeding lines were examined. Principal component analysis (PCA), principal coordinate analysis (PCoA), and cluster and STRUCTURE analyses found congruent results, which show that most of the examined cultivars and materials from Polish breeding programs formed major gene pools, that only some accessions derived from Strzelce Plant Breeding, and that foreign cultivars were outside of the main group. During the 120 year oat breeding process, only 67 alleles from the old gene pool were lost and replaced by 67 new alleles. The obtained results indicate that no erosion of genetic diversity was observed within the Polish native oat gene pool. Moreover, current oat breeding programs have introduced 673 new alleles into the gene pool relative to historical cultivars. The analysis also showed that most of the changes in relation to historical cultivars occurred within the A sub-genome with emphasis on chromosome 6A. The targeted changes were the rarest in the C sub-genome. This study showed that Polish oat breeding based mainly on traditional breeding methods—although focused on improving traits typical to this crop, i.e., enhancing the grain yield and quality and improving adaptability—did not significantly narrow the oat gene pool and in fact produced cultivars that are not only competitive in the European market but are also reservoirs of new alleles that were not found in the analyzed foreign materials.

In-Depth Understanding of the Genetic Control of Stripe Rust Resistance (Puccinia striiformis f. sp. tritici) Induced in Wheat (Triticum aestivum) by Trichoderma asperellum T34

Wheat stripe rust (caused by Puccinia striiformis f. tritici Erikss.) causes severe yield losses worldwide. Due to the continuous appearance of new stripe rust races, resistance has been broken in most of the highly resistant genotypes in Egypt and worldwide. Therefore, looking for new ways to resist such a severe disease is urgently needed. Trichoderma asperellum strain T34 has been known as an effective bioagent against many crop diseases. It exists naturally in Egyptian fields. Therefore, in our study, the effectiveness of strain T34 was tested as a bioagent against wheat stripe rust. For this purpose, 198 spring wheat genotypes were tested for their resistance against two different P. striiformis f. tritici populations collected from the Egyptian fields. The most highly aggressive P. striiformis f. tritici population was used to test the effectiveness of strain T34. Highly significant differences were found between strain T34 and stripe rust, suggesting the effectiveness of strain T34 in stripe rust resistance. A genome-wide association study identified 48 gene models controlling resistance under normal conditions and 46 gene models controlling strain T34-induced resistance. Of these gene models, only one common gene model was found, suggesting the presence of two different genetic systems controlling resistance under each condition. The pathways of the biological processes were investigated under both conditions. This study provided in-depth understanding of genetic control and, hence, will accelerate the future of wheat breeding programs for stripe rust resistance.

Genetic Variants Associated with Long-Terminal Repeats Can Diagnostically Classify Cannabis Varieties

Cannabis sativa (Cannabis) has recently been legalized in multiple countries globally for either its recreational or medicinal use. This, in turn, has led to a marked increase in the number of Cannabis varieties available for use in either market. However, little information currently exists on the genetic distinction between adopted varieties. Such fundamental knowledge is of considerable value and underpins the accelerated development of both a nascent pharmaceutical industry and the commercial recreational market. Therefore, in this study, we sought to assess genetic diversity across 10 Cannabis varieties by undertaking a reduced representation shotgun sequencing approach on 83 individual plants to identify variations which could be used to resolve the genetic structure of the assessed population. Such an approach also allowed for the identification of the genetic features putatively associated with the production of secondary metabolites in Cannabis. Initial analysis identified 3608 variants across the assessed population with phylogenetic analysis of this data subsequently enabling the confident grouping of each variety into distinct subpopulations. Within our dataset, the most diagnostically informative single nucleotide polymorphisms (SNPs) were determined to be associated with the long-terminal repeat (LTRs) class of retroelements, with 172 such SNPs used to fully resolve the genetic structure of the assessed population. These 172 SNPs could be used to design a targeted resequencing panel, which we propose could be used to rapidly screen different Cannabis plants to determine genetic relationships, as well as to provide a more robust, scientific classification of Cannabis varieties as the field moves into the pharmaceutical sphere.

Growth responses and genetic variation among highly ecologically diverse spring wheat genotypes grown under seawater stress

Most of the freshwaters worldwide are used for agriculture. Freshwater sources are expected to decline and will not suffice to support the food production needed for the growing population. Therefore, growing crops with seawater might constitute a solution. However, very little work has been done on the effect of seawater stress on wheat, an important cereal crop. The present study aimed to determine whether particular wheat genotypes provided better resistance to seawater stress. A set of 80 highly diverse spring wheat genotypes collected from different countries in Europe, Asia, Africa, North and South America was exposed to 50% seawater stress at the early growth stage. Four seeding shoot and root traits were scored for all genotypes. High genetic variations were found among all genotypes for the epicotyl length (EL), hypocotyl length (HL), number of radicles (NOR), and fresh weight (FW). Eight genotypes with high-performance scores of seedling traits were selected. The correlation analyses revealed highly significant correlations among all traits scored in this study. The strongest correlation was found between the NOR and the other seeding traits. Thus, the NOR might be an important adaptive trait for seawater tolerance. The genetic diversity among all genotypes was investigated based on genetic distance. A wide range of genetic distances among all genotypes was found. There was also a great genetic distance among the eight selected genotypes. In particular, the genetic distance between ATRI 5310 (France) and the other seven genotypes was the greatest. Such high genetic diversity might be utilized to select highly divergent genotypes for crossing in a future breeding program. The present study provides very useful information on the presence of different genetic resources in wheat for seawater tolerance.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

Diversity Assessment of the Montenegrin Maize Landrace Gene Pool Maintained in Two Gene Banks

Due to the loss of agro-biodiversity, there is a strong effort to find apparent and efficient mechanisms for the conservation and sustainable use of genetic diversity. A joint monitoring of the diversity and collections structure of the Montenegrin maize landraces conserved in the Serbian (MRIZPGB) and Montenegrin (MGB) gene banks has been conducted in order to improve the composition of the collections and to identify and eliminate possible redundancy. Based on a separate analysis of white- and yellow-orange maize landraces, it can be concluded that the diversity and evolution of distinct maize landraces grown and collected in Montenegro have been simultaneously shaped by both environmental (i.e., natural selection) and socially driven factors (farmers’ selection, migration and colonization processes of the human population). Although it has been determined that the authenticity and variability of the Montenegrin maize landraces gene pool have largely been preserved in the MRIZPGB collection, a significant amount of redundancy was observed. The obtained results will contribute to the cost-efficient conservation of the maize gene pool in the Montenegrin and Serbian gene banks. The recognized and well-preserved original variability of the MRIZPGB and MGB Montenegrin gene pool represents a valuable source for pre-breeding activities on broadening the white and flint maize breeding programmes under temperate conditions.