GWAS: Fast-forwarding gene identification and characterization in temperate Cereals: lessons from Barley - A review

Application of SVR-Mediated GWAS for Identification of Durable Genetic Regions Associated with Soybean Seed Quality Traits

Soybean (Glycine max L.) is an important food-grade strategic crop worldwide because of its high seed protein and oil contents. Due to the negative correlation between seed protein and oil percentage, there is a dire need to detect reliable quantitative trait loci (QTL) underlying these traits in order to be used in marker-assisted selection (MAS) programs. Genome-wide association study (GWAS) is one of the most common genetic approaches that is regularly used for detecting QTL associated with quantitative traits. However, the current approaches are mainly focused on estimating the main effects of QTL, and, therefore, a substantial statistical improvement in GWAS is required to detect associated QTL considering their interactions with other QTL as well. This study aimed to compare the support vector regression (SVR) algorithm as a common machine learning method to fixed and random model circulating probability unification (FarmCPU), a common conventional GWAS method in detecting relevant QTL associated with soybean seed quality traits such as protein, oil, and 100-seed weight using 227 soybean genotypes. The results showed a significant negative correlation between soybean seed protein and oil concentrations, with heritability values of 0.69 and 0.67, respectively. In addition, SVR-mediated GWAS was able to identify more relevant QTL underlying the target traits than the FarmCPU method. Our findings demonstrate the potential use of machine learning algorithms in GWAS to detect durable QTL associated with soybean seed quality traits suitable for genomic-based breeding approaches. This study provides new insights into improving the accuracy and efficiency of GWAS and highlights the significance of using advanced computational methods in crop breeding research.

Genetic dissection of protein and starch during wheat grain development using QTL mapping and GWAS

Protein, starch, and their components are important for wheat grain yield and end-products, which are affected by wheat grain development. Therefore, QTL mapping and a genome-wide association study (GWAS) of grain protein content (GPC), glutenin macropolymer content (GMP), amylopectin content (GApC), and amylose content (GAsC) were performed on wheat grain development at 7, 14, 21, and 28 days after anthesis (DAA) in two environments using a recombinant inbred line (RIL) population of 256 stable lines and a panel of 205 wheat accessions. A total of 29 unconditional QTLs, 13 conditional QTLs, 99 unconditional marker-trait associations (MTAs), and 14 conditional MTAs significantly associated (p < 10^-4) with four quality traits were found to be distributed on 15 chromosomes, with the phenotypic variation explained (PVE) ranging from 5.35% to 39.86%. Among these genomic variations, three major QTLs [QGPC3B, QGPC2A, and QGPC(S3|S2)3B] and SNP clusters on the 3A and 6B chromosomes were detected for GPC, and the SNP TA005876-0602 was stably expressed during the three periods in the natural population. The QGMP3B locus was detected five times in three developmental stages in two environments with 5.89%-33.62% PVE, and SNP clusters for GMP content were found on the 3A and 3B chromosomes. For GApC, the QGApC3B.1 locus had the highest PVE of 25.69%, and SNP clusters were found on chromosomes 4A, 4B, 5B, 6B, and 7B. Four major QTLs of GAsC were detected at 21 and 28 DAA. Most interestingly, both QTL mapping and GWAS analysis indicated that four chromosomes (3B, 4A, 6B, and 7A) were mainly involved in the development of protein, GMP, amylopectin, and amylose synthesis. Of these, the wPt-5870-wPt-3620 marker interval on chromosome 3B seemed to be most important because it played an important role in the synthesis of GMP and amylopectin before 7 DAA, in the synthesis of protein and GMP from 14 to 21 DAA, and in the development of GApC and GAsC from 21 to 28 DAA. Using the annotation information of IWGSC Chinese Spring RefSeq v1.1 genome assembly, we predicted 28 and 69 candidate genes for major loci from QTL mapping and GWAS, respectively. Most of them have multiple effects on protein and starch synthesis during grain development. These results provide new insights and information for the potential regulatory network between grain protein and starch synthesis.

Genetic analysis of protein content and oil content in soybean by genome-wide association study

Soybean seed protein content (PC) and oil content (OC) have important economic value. Detecting the loci/gene related to PC and OC is important for the marker-assisted selection (MAS) breeding of soybean. To detect the stable and new loci for PC and OC, a total of 320 soybean accessions collected from the major soybean-growing countries were used to conduct a genome-wide association study (GWAS) by resequencing. The PC ranged from 37.8% to 46.5% with an average of 41.1% and the OC ranged from 16.7% to 22.6% with an average of 21.0%. In total, 23 and 29 loci were identified, explaining 3.4%-15.4% and 5.1%-16.3% of the phenotypic variations for PC and OC, respectively. Of these, eight and five loci for PC and OC, respectively, overlapped previously reported loci and the other 15 and 24 loci were newly identified. In addition, nine candidate genes were identified, which are known to be involved in protein and oil biosynthesis/metabolism, including lipid transport and metabolism, signal transduction, and plant development pathway. These results uncover the genetic basis of soybean protein and oil biosynthesis and could be used to accelerate the progress in enhancing soybean PC and OC.

Genotypic variation of tocopherol content in a representative genetic population and genome-wide association study on tocopherol in rapeseed (Brassica napus)

Tocopherols (Tocs) are a kind of lipid-soluble substance required for the normal physiological function of mammals, particularly their antioxidant capacity. Rapeseed (Brassica napus) oil is an important source of exogenous Tocs. However, the genotypic differences in the total Toc contents, the Toc composition in the seeds, and the molecular markers associated with the seed Toc remain largely unknown. Here, we selected 290 rapeseed accessions based on the resequencing of 991 genomes in a worldwide collection of rapeseed germplasm. The contents of the four Toc isoforms, namely, α-, β-, γ-, and δ-Tocs, were also measured. Results show that the total Toc content and the γ-/α-Toc ratio varied greatly across the accessions from 85.34 to 387.00 mg/mg and 0.65 to 5.03, respectively. Furthermore, we conducted genome-wide association studies on the Tocs, which identified 28 and 73 single nucleotide polymorphisms significantly associated with the variation of total Toc content and γ-/α-Toc ratio, respectively. Bna.C02.VTE4, a putative orthologue of Arabidopsis VITAMIN E DEFICIENT 4, was tightly associated with the γ-/α-Toc ratio. This study recommends specific genetic materials with particularly high total Toc and/or low γ-/α-Toc ratio and the molecular markers and haplotypes associated with these quality traits for rapeseed breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01394-0.

Effectiveness of TaDreb-B1 and 1-FEH w3 KASP Markers in Spring and Winter Wheat Populations for Marker-Assisted Selection to Improve Drought Tolerance

Due to the advances in DNA markers, kompetitive allele-specific PCR (KASP) markers could accelerate breeding programs and genetically improve drought tolerance. Two previously reported KASP markers, TaDreb-B1 and 1-FEH w3, were investigated in this study for the marker-assisted selection (MAS) of drought tolerance. Two highly diverse spring and winter wheat populations were genotyped using these two KASP markers. The same populations were evaluated for drought tolerance at seedling (drought stress) and reproductive (normal and drought stress) growth stages. The single-marker analysis revealed a high significant association between the target allele of 1-FEH w3 and drought susceptibility in the spring population, while the marker-trait association was not significant in the winter population. The TaDreb-B1 marker did not have any highly significant association with seedling traits, except the sum of leaf wilting in the spring population. For field experiments, SMA revealed very few negative and significant associations between the target allele of the two markers and yield traits under both conditions. The results of this study revealed that the use of TaDreb-B1 provided better consistency in improving drought tolerance than 1-FEH w3.

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.

A novel soybean malectin-like receptor kinase-encoding gene, GmMLRK1, provides resistance to soybean mosaic virus

Soybean mosaic virus (SMV) severely damages soybean [Glycine max (L.) Merr.] yield and seed quality. Moreover, the underlying genetic determinants of resistance to SMV remain largely unknown. Here, we performed a genome-wide association study (GWAS) of SMV resistance in a panel of 219 diverse soybean accessions across four environments and identified a new resistance-related gene, GmMLRK1, at the major resistance locus Rsv4 on chromosome 2. GmMLRK1 encodes a malectin-like receptor kinase (RK) that was induced earlier and to a greater degree in leaves of the SMV-resistant cultivar Kefeng No. 1 than in those of the susceptible cultivar Nannong 1138-2 after inoculation. We demonstrated that soybean plants overexpressing GmMLRK1 show broad-spectrum resistance to both strains SC7 and SC3 on the basis of reduced viral accumulation, increased reactive oxygen species production, and local cell death associated with the hypersensitive response. In contrast, GmMLRK1 knockout mutants were more susceptible to both pathotypes. Haplotype analysis revealed the presence of five haplotypes (H1-H5) within the soybean population, and only H1 provided SMV resistance, which was independent of its tightly linked SMV resistance gene RNase-H at the same locus. These results report a novel gene that adds new understanding of SMV resistance and can be used for breeding resistant soybean accessions.

Genome-wide survey identified superior and rare haplotypes for plant height in the north-eastern soybean germplasm of China

The proper and efficient utilization of natural genetic diversity can significantly impact crop improvements. Plant height is a quantitative trait governing the plant type as well as the yield and quality of soybean. Here, we used a combined approach including a genome-wide association study (GWAS) and haplotype and candidate gene analyses to explore the genetic basis of plant height in diverse natural soybean populations. For the GWAS analysis, we used the whole-genome resequencing data of 196 diverse soybean cultivars collected from different accumulated temperature zones of north-eastern China to detect the significant single-nucleotide polymorphisms (SNPs) associated with plant height across three environments (E1, E2, and E3). A total of 33 SNPs distributed on four chromosomes, viz., Chr.02, Chr.04, Chr.06, and Chr.19, were identified to be significantly associated with plant height across the three environments. Among them, 23 were consistently detected in two or more environments and the remaining 10 were identified in only one environment. Interestingly, all the significant SNPs detected on the respective chromosomes fell within the physical interval of linkage disequilibrium (LD) decay (± 38.9 kb). Hence, these genomic regions were considered to be four quantitative trait loci (QTLs), viz., qPH2, qPH4, qPH6, and qPH19, regulating plant height. Moreover, the genomic region flanking all significant SNPs on four chromosomes exhibited strong LD. These significant SNPs thus formed four haplotype blocks, viz., Hap-2, Hap-4, Hap-6, and Hap-19. The number of haplotype alleles underlying each block varied from four to six, and these alleles regulate the different phenotypes of plant height ranging from dwarf to extra-tall heights. Nine candidate genes were identified within the four haplotype blocks, and these genes were considered putative candidates regulating soybean plant height. Hence, these stable QTLs, superior haplotypes, and candidate genes (after proper validation) can be deployed for the development of soybean cultivars with desirable plant heights.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01363-7.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

Nutritional Genomic Approach for Improving Grain Protein Content in Wheat

Grain protein content (GPC) is a key aspect of grain quality, a major determinant of the flour functional properties and grain nutritional value of bread wheat. Exploiting diverse germplasms to identify genes for improving crop performance and grain nutritional quality is needed to enhance food security. Here, we evaluated GPC in a panel of 255 Triticum aestivum L. accessions from 27 countries. GPC determined in seeds from three consecutive crop seasons varied from 8.6 to 16.4% (11.3% on average). Significant natural phenotypic variation in GPC among genotypes and seasons was detected. The population was evaluated for the presence of the trait-linked single nucleotide polymorphism (SNP) markers via a genome-wide association study (GWAS). GWAS analysis conducted with calculated best linear unbiased estimates (BLUEs) of phenotypic data and 90 K SNP array using the fixed and random model circulating probability unification (FarmCPU) model identified seven significant genomic regions harboring GPC-associated markers on chromosomes 1D, 3A, 3B, 3D, 4B and 5A, of which those on 3A and 3B shared associated SNPs with at least one crop season. The verified SNP-GPC associations provide new promising genomic signals on 3A (SNPs: Excalibur_c13709_2568 and wsnp_Ku_c7811_13387117) and 3B (SNP: BS00062734_51) underlying protein improvement in wheat. Based on the linkage disequilibrium for significant SNPs, the most relevant candidate genes within a 4 Mbp-window included genes encoding a subtilisin-like serine protease; amino acid transporters; transcription factors; proteins with post-translational regulatory functions; metabolic proteins involved in the starch, cellulose and fatty acid biosynthesis; protective and structural proteins, and proteins associated with metal ions transport or homeostasis. The availability of molecular markers within or adjacent to the sequences of the detected candidate genes might assist a breeding strategy based on functional markers to improve genetic gains for GPC and nutritional quality in wheat.

Genome-wide association studies reveal novel loci for resistance to groundnut rosette disease in the African core groundnut collection

We identified markers associated with GRD resistance after screening an Africa-wide core collection across three seasons in Uganda Groundnut is cultivated in several African countries where it is a major source of food, feed and income. One of the major constraints to groundnut production in Africa is groundnut rosette disease (GRD), which is caused by a complex of three agents: groundnut rosette assistor luteovirus, groundnut rosette umbravirus and its satellite RNA. Despite several years of breeding for GRD resistance, the genetics of the disease is not fully understood. The objective of the current study was to use the African core collection to establish the level of genetic variation in their response to GRD, and to map genomic regions responsible for the observed resistance. The African groundnut core genotypes were screened across two GRD hotspot locations in Uganda (Nakabango and Serere) for 3 seasons. The Area Under Disease Progress Curve combined with 7523 high quality SNPs were analyzed to establish marker-trait associations (MTAs). Genome-Wide Association Studies based on Enriched Compressed Mixed Linear Model detected 32 MTAs at Nakabango: 21 on chromosome A04, 10 on B04 and 1 on B08. Two of the significant markers were localised on the exons of a putative TIR-NBS-LRR disease resistance gene on chromosome A04. Our results suggest the likely involvement of major genes in the resistance to GRD but will need to be further validated with more comprehensive phenotypic and genotypic datasets. The markers identified in the current study will be developed into routine assays and validated for future genomics-assisted selection for GRD resistance in groundnut.

Genome-wide association analysis of spike and kernel traits in the U.S. hard winter wheat

A better understanding of the genetic control of spike and kernel traits that have higher heritability can help in the development of high-yielding wheat varieties. Here, we identified the marker-trait associations (MTAs) for various spike- and kernel-related traits in winter wheat (Triticum aestivum L.) through genome-wide association studies (GWAS). An association mapping panel comprising 297 hard winter wheat accessions from the U.S. Great Plains was evaluated for eight spike- and kernel-related traits in three different environments. A GWAS using 15,590 single-nucleotide polymorphisms (SNPs) identified a total of 53 MTAs for seven spike- and kernel-related traits, where the highest number of MTAs were identified for spike length (16) followed by the number of spikelets per spike (15) and spikelet density (11). Out of 53 MTAs, 14 were considered to represent stable quantitative trait loci (QTL) as they were identified in multiple environments. Five multi-trait MTAs were identified for various traits including the number of spikelets per spike (NSPS), spikelet density (SD), kernel width (KW), and kernel area (KA) that could facilitate the pyramiding of yield-contributing traits. Further, a significant additive effect of accumulated favorable alleles on the phenotype of four spike-related traits suggested that breeding lines and cultivars with a higher number of favorable alleles could be a valuable resource for breeders to improve yield-related traits. This study improves the understanding of the genetic basis of yield-related traits in hard winter wheat and provides reliable molecular markers that will facilitate marker-assisted selection (MAS) in wheat breeding programs.

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

The worldwide production of eggplant is estimated at about 58 Mt, with China, India and Egypt being the major producing countries. Breeding efforts in the species have mainly focused on increasing productivity, abiotic and biotic tolerance/resistance, shelf-life, the content of health-promoting metabolites in the fruit rather than decreasing the content of anti-nutritional compounds in the fruit. From the literature, we collected information on mapping quantitative trait loci (QTLs) affecting eggplant's traits following a biparental or multi-parent approach as well as genome-wide association (GWA) studies. The positions of QTLs were lifted according to the eggplant reference line (v4.1) and more than 700 QTLs were identified, here organized into 180 quantitative genomic regions (QGRs). Our findings thus provide a tool to: (i) determine the best donor genotypes for specific traits; (ii) narrow down QTL regions affecting a trait by combining information from different populations; (iii) pinpoint potential candidate genes.

A genome-wide association study (GWAS) of the personality constructs in CPAI-2 in Taiwanese Hakka populations

Here in this study we adopted genome-wide association studies (GWAS) to investigate the genetic components of the personality constructs in the Chinese Personality Assessment Inventory 2 (CPAI-2) in Taiwanese Hakka populations, who are likely the descendants of a recent admixture between a group of Chinese immigrants with high emigration intention and a group of the Taiwanese aboriginal population generally without it. A total of 279 qualified participants were examined and genotyped by an Illumina array with 547,644 SNPs to perform the GWAS. Although our sample size is small and that unavoidably limits our statistical power (Type 2 error but not Type 1 error), we still found three genomic regions showing strong association with Enterprise, Diversity, and Logical vs. Affective Orientation, respectively. Multiple genes around the identified regions were reported to be nervous system related, which suggests that genetic variants underlying the certain personalities should indeed exist in the nearby areas. It is likely that the recent immigration and admixture history of the Taiwanese Hakka people created strong linkage disequilibrium between the emigration intention-related genetic variants and their neighboring genetic markers, so that we could identify them despite with only limited statistical power.

Association mapping unravels the genetics controlling seedling drought stress tolerance in winter wheat

Drought is a major constraint in wheat (Triticum aestivum L.) grain yield. The present work aimed to identify quantitative trait nucleotides (QTNs)/ candidate genes influencing drought tolerance-related traits at the seedling stage in 261 accessions of a diverse winter wheat panel. Seeds from three consecutive years were exposed to polyethylene glycol 12% (PEG-6000) and a control treatment (distilled water). The Farm-CPU method was used for the association analysis with 17,093 polymorphic SNPs. PEG treatment reduced shoot length (SL) (-36.3%) and root length (RL) (-11.3%) compared with control treatments, while the coleoptile length (CL) was increased by 11% under drought conditions, suggesting that it might be considered as an indicator of stress-tolerance. Interestingly, we revealed 70 stable QTN across 17 chromosomes. Eight QTNs related to more than one trait were detected on chromosomes 1B, 2A (2), 2B, 2D, 4B, 7A, and 7B and located nearby or inside candidate genes within the linkage disequilibrium (LD) interval. For instance, the QTN on chromosome 2D is located inside the gene TraesCS2D02G133900 that controls the variation of CL_S and SL_C. The allelic variation at the candidate genes showed significant influence on the associated traits, demonstrating their role in controlling the natural variation of multi-traits of drought stress tolerance. The gene expression of these candidate genes under different stress conditions validates their biological role in stress tolerance. Our findings offer insight into understanding the genetic factors and diverse mechanisms in response to water shortage conditions that are important for wheat improvement and adaptation at early developmental stages.

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.

cgSSR marker-based genome-wide association study identified genomic regions for panicle characters and yield in rice (Oryza sativa L.)

BACKGROUND

To improve production efficiency, positive alleles corresponding to yield-related attributes must be accumulated in a single elite background. We designed and used cgSSR markers, which are superior to random SSR markers in genome-wide association study, to identify genomic regions that contribute to panicle characters and grain yield in this study.

RESULTS

As evidenced by the high polymorphic information content value and gene diversity coefficient, the new cgSSR markers were determined to be highly informative. These cgSSR markers were employed to generate genotype data for an association panel evaluated for four panicle characters and grain yield over three seasons. For five traits, 17 significant marker-trait associations on six chromosomes were discovered. The percentage of phenotypic variance that could be explained ranged from 4% to 13%. Unrelated gene-derived markers had a strong association with target traits as well.

CONCLUSION

Trait-associated cgSSR markers derived from corresponding or related genes ensure their utility in direct allele selection, while other linked markers aid in allele selection indirectly by altering the phenotype of interest. Through a marker-assisted breeding approach, these marker-trait associations can be leveraged to accumulate favourable alleles for yield enhancement in rice. © 2022 Society of Chemical Industry.

Digitalizing breeding in plants: A new trend of next-generation breeding based on genomic prediction

As the world's population grows and food needs diversification, the demand for cereals and horticultural crops with beneficial traits increases. In order to meet a variety of demands, suitable cultivars and innovative breeding methods need to be developed. Breeding methods have changed over time following the advance of genetics. With the advent of new sequencing technology in the early 21st century, predictive breeding, such as genomic selection (GS), emerged when large-scale genomic information became available. GS shows good predictive ability for the selection of individuals with traits of interest even for quantitative traits by using various types of the whole genome-scanning markers, breaking away from the limitations of marker-assisted selection (MAS). In the current review, we briefly describe the history of breeding techniques, each breeding method, various statistical models applied to GS and methods to increase the GS efficiency. Consequently, we intend to propose and define the term digital breeding through this review article. Digital breeding is to develop a predictive breeding methods such as GS at a higher level, aiming to minimize human intervention by automatically proceeding breeding design, propagating breeding populations, and to make selections in consideration of various environments, climates, and topography during the breeding process. We also classified the phases of digital breeding based on the technologies and methods applied to each phase. This review paper will provide an understanding and a direction for the final evolution of plant breeding in the future.

Boosting of tau protein aggregation by CD40 and CD48 gene expression in Alzheimer's disease

Neurodegenerative diseases result from the interplay of abnormal gene expression and various pathological factors. Therefore, a disease-specific integrative genetic approach is required to understand the complexities and causes of target diseases. Recent studies have identified the correlation between genes encoding several transmembrane proteins, such as the cluster of differentiation (CD) and Alzheimer's disease (AD) pathogenesis. In this study, CD48 and CD40 gene expression in AD, a neurodegenerative disease, was analyzed to infer this link. Total RNA sequencing was performed using an Alzheimer's disease mouse model brain and blood, and gene expression was determined using a genome-wide association study (GWAS). We observed a marked elevation of CD48 and CD40 genes in Alzheimer's disease. Indeed, the upregulation of both CD48 and CD40 genes was significantly increased in the severe Alzheimer's disease group. With the elevation of CD48 and CD40 genes in Alzheimer's disease, associations of protein levels were also markedly increased in tissues. In addition, overexpression of CD48 and CD40 genes triggered tau aggregation, and co-expression of these genes accelerated aggregation. The nuclear factor kappa B (NF-ĸB) signaling pathway was enriched by CD48 and CD40 gene expression: it was also associated with tau pathology. Our data suggested that the CD48 and CD40 genes are novel AD-related genes, and this approach may be useful as a diagnostic or therapeutic target for the disease.

A Haplotype-Based GWAS Identified Trait-Improving QTL Alleles Controlling Agronomic Traits under Contrasting Nitrogen Fertilization Treatments in the MAGIC Wheat Population WM-800

The multi-parent-advanced-generation-intercross (MAGIC) population WM-800 was developed by intercrossing eight modern winter wheat cultivars to enhance the genetic diversity present in breeding populations. We cultivated WM-800 during two seasons in seven environments under two contrasting nitrogen fertilization treatments. WM-800 lines exhibited highly significant differences between treatments, as well as high heritabilities among the seven agronomic traits studied. The highest-yielding WM-line achieved an average yield increase of 4.40 dt/ha (5.2%) compared to the best founder cultivar Tobak. The subsequent genome-wide-association-study (GWAS), which was based on haplotypes, located QTL for seven agronomic traits including grain yield. In total, 40, 51, and 46 QTL were detected under low, high, and across nitrogen treatments, respectively. For example, the effect of QYLD_3A could be associated with the haplotype allele of cultivar Julius increasing yield by an average of 4.47 dt/ha (5.2%). A novel QTL on chromosome 2B exhibited pleiotropic effects, acting simultaneously on three-grain yield components (ears-per-square-meter, grains-per-ear, and thousand-grain-weight) and plant-height. These effects may be explained by a member of the nitrate-transporter-1 (NRT1)/peptide-family, TaNPF5.34, located 1.05 Mb apart. The WM-800 lines and favorable QTL haplotypes, associated with yield improvements, are currently implemented in wheat breeding programs to develop advanced nitrogen-use efficient wheat cultivars.

Genome-Wide Association Study Revealed SNP Alleles Associated with Seed Size Traits in African Yam Bean (Sphenostylis stenocarpa (Hochst ex. A. Rich.) Harms)

Seed size is an important yield and quality-determining trait in higher plants and is also crucial to their evolutionary fitness. In African yam bean (AYB), seed size varies widely among different accessions. However, the genetic basis of such variation has not been adequately documented. A genome-wide marker-trait association study was conducted to identify genomic regions associated with four seed size traits (seed length, seed width, seed thickness, and 100-seed weight) in a panel of 195 AYB accessions. A total of 5416 SNP markers were generated from the diversity array technology sequence (DArTseq) genotype-by-sequencing (GBS)- approach, in which 2491 SNPs were retained after SNP quality control and used for marker-trait association analysis. Significant phenotypic variation was observed for the traits. Broad-sense heritability ranged from 50.0% (seed width) to 66.4% (seed length). The relationships among the traits were positive and significant. Genome-wide association study (GWAS) using the general linear model (GLM) and the mixed linear model (MLM) approaches identified 12 SNP markers significantly associated with seed size traits across the six test environments. The 12 makers explained 6.5-10.8% of the phenotypic variation. Two markers (29420334|F|0-52:C>G-52:C>G and 29420736|F|0-57:G>T-57:G>T) with pleiotropic effects associated with seed width and seed thickness were found. A candidate gene search identified five significant markers (100026424|F|0-37:C>T-37:C>T, 100041049|F|0-42:G>C-42:G>C, 100034480|F|0-31:C>A-31:C>A, 29420365|F|0-55:C>G-55:C>G, and 29420736|F|0-57:G>T-57:G>T) located close to 43 putative genes whose encoding protein products are known to regulate seed size traits. This study revealed significant makers not previously reported for seed size in AYB and could provide useful information for genomic-assisted breeding in AYB.

Diversity and determinants of bitterness, astringency, and fat content in cultivated Nacional and native Amazonian cocoa accessions from Ecuador

Cocoa (Theobroma cacao L.) is the only tree that can produce cocoa. Cocoa beans are highly sought after by chocolate makers to produce chocolate. Cocoa can be fine aromatic, characterized by floral and fruity notes, or it can be described as standard cocoa with a more pronounced cocoa aroma and bitterness. In this study, the genetic and biochemical determinants of sensorial notes and nonvolatile compounds related to bitterness, astringency, fat content, and protein content will be investigated in two populations: a cultivated modern Nacional population and a population of cocoa accessions collected recently in the Ecuadorian South Amazonia area of origin of the Nacional ancestral variety. For this purpose, a genome-wide association study (GWAS) was carried out on both populations, with results of biochemical compounds evaluated by near-infrared spectroscopy (NIRS) assays and with sensory evaluations. Twenty areas of associations were detected for sensorial data especially bitterness and astringency. Fifty-three areas of associations were detected linked to nonvolatile compounds. A total of 81 candidate genes could be identified in the areas of the association.

Genome-wide association mapping for field and seedling resistance to the emerging Puccinia graminis f. sp. tritici race TTRTF in wheat

Stem rust of wheat (Triticum spp.), caused by Puccinia graminis f. sp. tritici (Pgt), is one of the most impactful wheat diseases because of its threat to global wheat production. While disease mitigation has primarily been achieved through the deployment of resistant wheat varieties, emerging new virulent races continue to pose risks to the crop. For example, races such as Ug99 (TTKSK), TKTTF, and TTRTF have caused epidemics in different wheat growing regions of the world in recent years. A continual search for new and effective sources of resistance is therefore necessary to safeguard wheat production. This study assessed a breeding panel from the Ethiopian Institute of Agricultural Research (EIAR) wheat breeding program for seedling and field plant resistance to TTRTF and reports genomic regions conferring resistance to TTRTF. Trait correlations (r) were medium to strong (range = .38-.71) and heritabilities were moderate (.32-.56). Association analysis for resistance to TTRTF resulted in detection of 20 markers in 11 chromosomes; the marker S1B_175439851 was associated with resistance at both seedling and adult plant stages. Models with two to four QTL combinations reduced seedling and field disease severity by 12-48 and 9-17%, respectively. Genomic prediction for TTRTF resistance resulted in low to moderately-high predictions (mean correlations of .25-.47). Identification of resistant lines and QTL in the EIAR population is expected to assist in selection toward improved resistance to TTRTF. Specifically, the application of genomic selection (GS) in identifying resistant lines in future related breeding populations will further assist breeding efforts against this new stem rust pathogen race.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

Validation of a QTL on Chromosome 1DS Showing a Major Effect on Salt Tolerance in Winter Wheat

Salt stress is one the most destructive abiotic stressors, causing yield losses in wheat worldwide. A prerequisite for improving salt tolerance is the identification of traits for screening genotypes and uncovering causative genes. Two populations of F₃ lines developed from crosses between sensitive and tolerant parents were tested for salt tolerance at the seedling stage. Based on their response, the offspring were classified as salt sensitive and tolerant. Under saline conditions, tolerant genotypes showed lower Na⁺ and proline content but higher K⁺, higher chlorophyll content, higher K⁺/Na⁺ ratio, higher PSII activity levels, and higher photochemical efficiency, and were selected for further molecular analysis. Five stress responsive QTL identified in a previous study were validated in the populations. A QTL on the short arm of chromosome 1D showed large allelic effects in several salt tolerant related traits. An expression analysis of associated candidate genes showed that TraesCS1D02G052200 and TraesCS5B02G368800 had the highest expression in most tissues. Furthermore, qRT-PCR expression analysis revealed that ZIP-7 had higher differential expressions under saline conditions compared to KefC, AtABC8 and 6-SFT. This study provides information on the genetic and molecular basis of salt tolerance that could be useful in development of salt-tolerant wheat varieties.

Identification of Genetic Variations and Candidate Genes Responsible for Stalk Sugar Content and Agronomic Traits in Fresh Corn via GWAS across Multiple Environments

The stem and leaves of fresh corn plants can be used as green silage or can be converted to biofuels, and the stalk sugar content and yield directly determine the application value of fresh corn. To identify the genetic variations and candidate genes responsible for the related traits in fresh corn, the genome-wide scan and genome-wide association analysis (GWAS) were performed. A total of 32 selective regions containing 172 genes were detected between sweet and waxy corns. Using the stalk sugar content and seven other agronomic traits measured in four seasons over two years, the GWAS identified ninety-two significant single nucleotide polymorphisms (SNPs). Most importantly, seven SNPs associated with the stalk sugar content were detected across multiple environments, which could explain 13.68–17.82% of the phenotypic variation. Accessions differing in genotype for certain significant SNPs showed significant variation in the stalk sugar content and other agronomic traits, and the expression levels of six important candidate genes were significantly different between two materials with different stalk sugar content. The genetic variations and candidate genes provide valuable resources for future studies of the molecular mechanism of the stalk sugar content and establish the foundation for molecular marker-assisted breeding of fresh corn.

Modulating physiological and transcriptional regulatory mechanisms for enhanced climate resilience in cereal crops

Climate change adversely affects the yield and productivity of cereal crops, which consequently impacts food security. Therefore, studying stress acclimation, particularly transcriptional patterns and morpho-physiological responses of cereal crops to different stresses, will provide insights into the molecular determinants underlying climate resilience. The availability of advanced tools and approaches has enabled the characterization of plants at morphological, physiological, biochemical, and molecular levels, which will lead to the identification of genomic regions regulating the stress responses at these levels. This will further facilitate using transgenic, breeding, or genome editing approaches to manipulate the identified regions (genes, alleles, or QTLs) to enhance stress resilience. Next-generation sequencing approaches have advanced the identification of causal genes and markers in the genomes through forward or reverse genetics. In this context, the review enumerates the progress of dissecting the molecular mechanisms underlying transcriptional and physiological responses of major cereals to climate-induced stresses. The review systematically discusses different tools and approaches available to study the response of plants to various stresses and identify the molecular determinants regulating stress-resilience. Further, the application of genomics-assisted breeding, transgene-, and targeted editing-based approaches for modulating the genetic determinants for enhanced climate resilience has been elaborated.

Genome-Wide SNP Markers Accelerate Perennial Forest Tree Breeding Rate for Disease Resistance through Marker-Assisted and Genome-Wide Selection

The ecological and economic importance of forest trees is evident and their survival is necessary to provide the raw materials needed for wood and paper industries, to preserve the diversity of associated animal and plant species, to protect water and soil, and to regulate climate. Forest trees are threatened by anthropogenic factors and biotic and abiotic stresses. Various diseases, including those caused by fungal pathogens, are one of the main threats to forest trees that lead to their dieback. Genomics and transcriptomics studies using next-generation sequencing (NGS) methods can help reveal the architecture of resistance to various diseases and exploit natural genetic diversity to select elite genotypes with high resistance to diseases. In the last two decades, QTL mapping studies led to the identification of QTLs related to disease resistance traits and gene families and transcription factors involved in them, including NB-LRR, WRKY, bZIP and MYB. On the other hand, due to the limitation of recombination events in traditional QTL mapping in families derived from bi-parental crosses, genome-wide association studies (GWAS) that are based on linkage disequilibrium (LD) in unstructured populations overcame these limitations and were able to narrow down QTLs to single genes through genotyping of many individuals using high-throughput markers. Association and QTL mapping studies, by identifying markers closely linked to the target trait, are the prerequisite for marker-assisted selection (MAS) and reduce the breeding period in perennial forest trees. The genomic selection (GS) method uses the information on all markers across the whole genome, regardless of their significance for development of a predictive model for the performance of individuals in relation to a specific trait. GS studies also increase gain per unit of time and dramatically increase the speed of breeding programs. This review article is focused on the progress achieved in the field of dissecting forest tree disease resistance architecture through GWAS and QTL mapping studies. Finally, the merit of methods such as GS in accelerating forest tree breeding programs is also discussed.

Elucidating the genetic architecture controlling antioxidant status and ionic balance in barley under salt stress

Association genetic analysis empowered us to identify candidate genes underlying natural variation of morpho-physiological, antioxidants, and grain yield-related traits in barley. Novel intriguing genomic regions were identified and dissected. Salinity stress is one of the abiotic stresses that influence the morpho-physiological, antioxidants, and yield-related traits in crop plants. The plants of a core set of 138 diverse barley accessions were analyzed after exposure to salt stress under field conditions during the reproductive phase. A genome-wide association scan (GWAS) was then conducted using 19,276 single nucleotide polymorphisms (SNPs) to uncover the genetic basis of morpho-physiological and grain-related traits. A wide range of responses to salt stress by the accessions was explored in the current study. GWAS detected 263 significantly associated SNPs with the antioxidants, K⁺/Na⁺ content ratio, and agronomic traits. Five genomic regions harbored interesting putative candidate genes within LD ± 1.2 Mbp. Choromosome 2H harbored many candidate genes associated with the antioxidants ascorbic acid (AsA) and glutathione (GSH), such as superoxide dismutase (SOD), ascorbate peroxidase (APX), and glutathione reductase (GR), under salt stress. Markedly, an A:C SNP at 153,773,211 bp on chromosome 7H is located inside the gene HORVU.MOREX.r3.7HG0676830 (153,772,300-153,774,057 bp) that was annotated as L-gulonolactone oxidase, regulating the natural variation of SOD_S and APX_S. The allelic variation at this SNP reveals a negative selection of accessions carrying the C allele, predominantly found in six-rowed spring landraces originating from Far-, Near-East, and central Asia carrying photoperiod sensitive alleles having lower activity of enzymatic antioxidants. The SNP-trait associations detected in the current study constitute a benchmark for developing molecular selection tools for antioxidant compound selection in barley.

Association mapping and candidate genes for physiological non-destructive traits: Chlorophyll content, canopy temperature, and specific leaf area under normal and saline conditions in wheat

Wheat plants experience substantial physiological adaptation when exposed to salt stress. Identifying such physiological mechanisms and their genetic control is especially important to improve its salt tolerance. In this study, leaf chlorophyll content (CC), leaf canopy temperature (CT), and specific leaf area (SLA) were scored in a set of 153 (103 having the best genotypic data were used for GWAS analysis) highly diverse wheat genotypes under control and salt stress. On average, CC and SLA decreased under salt stress, while the CT average was higher under salt stress compared to the control. CT was negatively and significantly correlated with CC under both conditions, while no correlation was found between SLA and CC and CT together. High genetic variation and broad-sense-heritability estimates were found among genotypes for all traits. The genome wide association study revealed important QTLs for CC under both conditions (10) and SLA under salt stress (four). These QTLs were located on chromosomes 1B, 2B, 2D, 3A, 3B, 5A, 5B, and 7B. All QTLs detected in this study had major effects with R² extending from 20.20% to 30.90%. The analysis of gene annotation revealed three important candidate genes (TraesCS5A02G355900, TraesCS1B02G479100, and TraesCS2D02G509500). These genes are found to be involved in the response to salt stress in wheat with high expression levels under salt stress compared to control based on mining in data bases.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

Genome-wide association mapping for root traits associated with frost tolerance in faba beans using KASP-SNP markers

Frost is an abiotic stress factor that threatens plant development and crop productivity not only in cold regions but also in temperate zones. Roots play an important role in plant growth during frost stress. Therefore, variation in root characteristics could be studied to improve frost tolerance in winter faba bean. The present study aimed to identify the genomic regions that control frost tolerance in a winter faba bean population by focusing on root-related traits. A set of 185 genotypes were tested for frost tolerance under artificial frost growth conditions at −16°C, −18°C, and −19°C in a growth chamber. Frost stress reduced the root-related parameters in all genotypes, with a wide variation among genotypes. A genome-wide association study identified nine novel single-nucleotide polymorphisms that are associated with the root-related traits. The most frost-tolerant genotypes were identified; two genotypes, S_028 and S_220, exhibited remarkable performance under frost stress. Moreover, they harbored all four of the alleles favorable for frost tolerance. Remarkably, two markers showed genetic pleiotropic effects with positive allele effects on root fresh matter and root dry matter. Thus, both genotypes can be implemented in a breeding program to provide the alleles for healthier roots under frost conditions to develop more frost-tolerant varieties, and the two markers can be used to screen large collections to select for frost tolerance. These results may provide novel insights for improving frost tolerance in faba beans and in other legume crops.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

Genotyping by Sequencing Advancements in Barley

Barley is considered an ideal crop to study cereal genetics due to its close relationship with wheat and diploid ancestral genome. It plays a crucial role in reducing risks to global food security posed by climate change. Genetic variations in the traits of interest in crops are vital for their improvement. DNA markers have been widely used to estimate these variations in populations. With the advancements in next-generation sequencing, breeders could access different types of genetic variations within different lines, with single-nucleotide polymorphisms (SNPs) being the most common type. However, genotyping barley with whole genome sequencing (WGS) is challenged by the higher cost and computational demand caused by the large genome size (5.5GB) and a high proportion of repetitive sequences (80%). Genotyping-by-sequencing (GBS) protocols based on restriction enzymes and target enrichment allow a cost-effective SNP discovery by reducing the genome complexity. In general, GBS has opened up new horizons for plant breeding and genetics. Though considered a reliable alternative to WGS, GBS also presents various computational difficulties, but GBS-specific pipelines are designed to overcome these challenges. Moreover, a robust design for GBS can facilitate the imputation to the WGS level of crops with high linkage disequilibrium. The complete exploitation of GBS advancements will pave the way to a better understanding of crop genetics and offer opportunities for the successful improvement of barley and its close relatives.

A genome-wide association study unveils key chromosome regions involved in determining sodium accumulation in wheat under conditions of low potassium supply

Improving nutrient use efficiency is an important objective in modern breeding programs. In this work, we examined potassium utilization efficiency (KUtE) and traits potentially related to it in a formerly genotyped, geographically diverse population of bread wheat (Triticum aestivum) under low potassium supply conditions. Our results unveil the existence of a large variation within the population for the traits examined. A genome-wide association study, based on a single-locus model, identified 15 markers associated with some of those traits. No marker-trait association was found using that tool for KUtE, but the use of a multi-locus approach suggested that additional marker-trait associations may be present, including whole-plant KUtE. Besides, the existence of a significant correlation between KUtE and sodium accumulation in shoots suggests the possibility of pyramiding traits associated with sodium homeostasis to improve this efficiency. In this regard, two discrete regions mapped on the long arm of chromosome 1B (1BLA and 1BLB) were associated with variation in sodium accumulation as detected with the single and multi-locus models used. Further exploration of the potential function of the genes placed in these regions, and their expression patterns, suggested likely candidates for this trait. Among the candidates placed in 1BLA region, we found TraesCS1B02G370500, TraesCS1B02G370600, and TraesCS1B02G370900, coding for putative Calcineurin B like proteins. Region 1BLB contain TraesCS1B02G388900 coding for a kinase and other genes including TraesCS1B02G389700, TraesCS1B02G389800 and TraesCS1B02G389900 coding for Ethylene-responsive transcription factors. The information here provided can be useful in breeding programs aimed to manipulate sodium accumulation through marker-assisted selection.

New Hope for Genome Editing in Cultivated Grasses: CRISPR Variants and Application

With the advent of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein (Cas) mediated genome editing, crop improvement has progressed significantly in recent years. In this genome editing tool, CRISPR-associated Cas nucleases are restricted to their target of DNA by their preferred protospacer adjacent motifs (PAMs). A number of CRISPR-Cas variants have been developed e.g. CRISPR-Cas9, -Cas12a and -Cas12b, with different PAM requirements. In this mini-review, we briefly explain the components of the CRISPR-based genome editing tool for crop improvement. Moreover, we intend to highlight the information on the latest development and breakthrough in CRISPR technology, with a focus on a comparison of major variants (CRISPR-Cas9, -Cas12a, and -Cas12b) to the newly developed CRISPR-SpRY that have nearly PAM-less genome editing ability. Additionally, we briefly explain the application of CRISPR technology in the improvement of cultivated grasses with regard to biotic and abiotic stress tolerance as well as improving the quality and yield.

Genetic Architecture of Powdery Mildew Resistance Revealed by a Genome-Wide Association Study of a Worldwide Collection of Flax (Linum usitatissimum L.)

Powdery mildew is one of the most important diseases of flax and is particularly prejudicial to its yield and oil or fiber quality. This disease, caused by the obligate biotrophic ascomycete Oïdium lini, is progressing in France. Genetic resistance of varieties is critical for the control of this disease, but very few resistance genes have been identified so far. It is therefore necessary to identify new resistance genes to powdery mildew suitable to the local context of pathogenicity. For this purpose, we studied a worldwide diversity panel composed of 311 flax genotypes both phenotyped for resistance to powdery mildew resistance over 2 years of field trials in France and resequenced. Sequence reads were mapped on the CDC Bethune reference genome revealing 1,693,910 high-quality SNPs, further used for both population structure analysis and genome-wide association studies (GWASs). A number of four major genetic groups were identified, separating oil flax accessions from America or Europe and those from Asia or Middle-East and fiber flax accessions originating from Eastern Europe and those from Western Europe. A number of eight QTLs were detected at the false discovery rate threshold of 5%, located on chromosomes 1, 2, 4, 13, and 14. Taking advantage of the moderate linkage disequilibrium present in the flax panel, and using the available genome annotation, we identified potential candidate genes. Our study shows the existence of new resistance alleles against powdery mildew in our diversity panel, of high interest for flax breeding program.

Reproductive-Stage Heat Stress in Cereals: Impact, Plant Responses and Strategies for Tolerance Improvement

Reproductive-stage heat stress (RSHS) poses a major constraint to cereal crop production by damaging main plant reproductive structures and hampering reproductive processes, including pollen and stigma viability, pollination, fertilization, grain setting and grain filling. Despite this well-recognized fact, research on crop heat stress (HS) is relatively recent compared to other abiotic stresses, such as drought and salinity, and in particular, RSHS studies in cereals are considerably few in comparison with seedling-stage and vegetative-stage-centered studies. Meanwhile, climate change-exacerbated HS, independently or synergistically with drought, will have huge implications on crop performance and future global food security. Fortunately, due to their sedentary nature, crop plants have evolved complex and diverse transient and long-term mechanisms to perceive, transduce, respond and adapt to HS at the molecular, cell, physiological and whole plant levels. Therefore, uncovering the molecular and physiological mechanisms governing plant response and tolerance to RSHS facilitates the designing of effective strategies to improve HS tolerance in cereal crops. In this review, we update our understanding of several aspects of RSHS in cereals, particularly impacts on physiological processes and yield; HS signal perception and transduction; and transcriptional regulation by heat shock factors and heat stress-responsive genes. We also discuss the epigenetic, post-translational modification and HS memory mechanisms modulating plant HS tolerance. Moreover, we offer a critical set of strategies (encompassing genomics and plant breeding, transgenesis, omics and agronomy) that could accelerate the development of RSHS-resilient cereal crop cultivars. We underline that a judicious combination of all of these strategies offers the best foot forward in RSHS tolerance improvement in cereals. Further, we highlight critical shortcomings to RSHS tolerance investigations in cereals and propositions for their circumvention, as well as some knowledge gaps, which should guide future research priorities. Overall, our review furthers our understanding of HS tolerance in plants and supports the rational designing of RSHS-tolerant cereal crop cultivars for the warming climate.

Genome-Wide Screening of Broad-Spectrum Resistance to Leaf Rust (Puccinia triticina Eriks) in Spring Wheat (Triticum aestivum L.)

Wheat leaf rust (LR) causes significant yield losses worldwide. In Egypt, resistant cultivars began to lose their efficiency in leaf rust resistance. Therefore, a diverse spring wheat panel was evaluated at the seedling stage to identify new sources of broad-spectrum seedling resistance against the Egyptian Puccinia triticina (Pt) races. In three different experiments, seedling evaluation was done using Pt spores collected from different fields and growing seasons. Highly significant differences were found among experiments confirming the presence of different races population in each experiment. Highly significant differences were found among the tested genotypes confirming the ability to select superior genotypes. Genome-wide association study (GWAS) was conducted for each experiment and a set of 87 markers located within 48 gene models were identified. The identified gene models were associated with disease resistance in wheat. Five gene models were identified to resist all Pt races in at least two experiments and could be identified as stable genes under Egyptian conditions. Ten genotypes from five different countries were stable against all the tested Pt races but showed different degrees of resistance.

Identification of Putative SNP Markers Associated with Resistance to Egyptian Loose Smut Race(s) in Spring Barley

Loose smut (LS) disease is a serious problem that affects barley yield. Breeding of resistant cultivars and identifying new genes controlling LS has received very little attention. Therefore, it is important to understand the genetic basis of LS control in order to genetically improve LS resistance. To address this challenge, a set of 57 highly diverse barley genotypes were inoculated with Egyptian loose smut race(s) and the infected seeds/plants were evaluated in two growing seasons. Loose smut resistance (%) was scored on each genotype. High genetic variation was found among all tested genotypes indicating considerable differences in LS resistance that can be used for breeding. The broad-sense heritability (H²) of LS (0.95) was found. Moreover, genotyping-by-sequencing (GBS) was performed on all genotypes and generated in 16,966 SNP markers which were used for genetic association analysis using single-marker analysis. The analysis identified 27 significant SNPs distributed across all seven chromosomes that were associated with LS resistance. One SNP (S6_17854595) was located within the HORVU6Hr1G010050 gene model that encodes a protein kinase domain-containing protein (similar to the Un8 LS resistance gene, which contains two kinase domains). A TaqMan marker (0751D06 F6/R6) for the Un8 gene was tested in the diverse collection. The results indicated that none of the Egyptian genotypes had the Un8 gene. The result of this study provided new information on the genetic control of LS resistance. Moreover, good resistance genotypes were identified and can be used for breeding cultivars with improved resistance to Egyptian LS.

Genome-Wide Association Mapping of Grain Metal Accumulation in Wheat

Increasing wheat grain yield while ignoring grain quality and metal accumulation can result in metal deficiencies, particularly in countries where bread wheat accounts for the majority of daily dietary regimes. When the accumulation level exceeds a certain threshold, it becomes toxic and causes various diseases. Biofortification is an effective method of ensuring nutritional security. We screened 200 spring wheat advanced lines from the wheat association mapping initiative for Mn, Fe, Cu, Zn, Ni, and Cd concentrations. Interestingly, high-yielding genotypes had high essential metals, such as Mn, Fe, Cu, and Zn, but low levels of toxic metals, such as Ni and Cd. Positive correlations were found between all metals except Ni and Cd, where no correlation was found. We identified 142 significant SNPs, 26 of which had possible pleiotropic effects on two or more metals. Several QTLs co-located with previously mapped QTL for the same or other metals, whereas others were new. Our findings contribute to wheat genetic biofortification through marker-assisted selection, ensuring nutritional security in the long run.

Genomic Regions, Molecular Markers, and Flanking Genes of Metribuzin Tolerance in Wheat (Triticum aestivum L.)

Understanding the genetics of metribuzin (a group C herbicide) tolerance in wheat is vital in developing tolerant cultivars to improve wheat productivity in dryland farming systems. This study investigated metribuzin tolerance in wheat by conducting a Genome-wide Association Studies (GWAS) with a panel of 150 wheat genotypes of diverse genetic backgrounds and genotyped them with the wheat 90 K SNP genotyping assay. The phenotyping was conducted in a temperature-controlled glasshouse at the University of Western Australia (UWA). Genotypes were sprayed with a metribuzin dose of 400 grams of active ingredient (g. a.i.) ha^-1 as pre-emergent in a specialized spraying cabinet and transferred to the glasshouse where the tolerance level of the genotypes was assessed by measuring the relative reduction in chlorophyll content of the leaves. The decrease in chlorophyll content of the treated plants compared to the control was regarded as the phytotoxic effects of metribuzin. GWAS analysis following a mixed linear model revealed 19 genomic regions with significant marker-trait associations (MTAs), including ten on chromosome 6A, three on chromosome 2B, and one on chromosomes 3A, 5B, 6B 6D, 7A, and 7B, respectively. Sequences of the significant markers were blasted against the wheat genome, IWGSC RefSeq V1.0, and candidate genes having annotations related to herbicide tolerance in wheat, especially in pathways reported to be involved in metribuzin tolerance, such as cytochrome P450 pathways and ATP Binding Cassette (ABC) superfamilies, were identified in these genomic regions. These included TraesCS6A01G028800, TraesCS6A02G353700, TraesCS6A01G326200, TraesCS7A02G331000, and TraesCS2B01G465200. These genomic regions were validated on 30 top tolerant and 30 most susceptible genotypes using the five closest SSR makers to the flanked SNPs. Sufficient polymorphism was detected on two markers (wms193 and barc1036) that were found to differentiate between the susceptible and tolerant alleles and a t-test analysis of the phenotypic data shows a significant (value of p < 0.001) difference suggesting that these markers can be used for marker-assisted selection (MAS) in metribuzin studies and wheat breeding programs.

Combined GWAS and QTL mapping revealed candidate genes and SNP network controlling recovery and tolerance traits associated with drought tolerance in seedling winter wheat

To date, very little research on drought tolerance has been conducted at the seedling stage in winter wheat. In this study, two types of traits, namely tolerance and recovery traits, associated with drought tolerance were scored in biparental mapping population (BPP) and association mapping population (A-set). The results of this study revealed no or weak significant correlation between the two types of traits. Based on GWAS and QTL mapping analyses, all QTLs associated with recovery traits were completely different from those associated with tolerance traits except one QTL in each population that was found to be associated with one tolerance trait and one recovery trait. The analysis of SNP and gene networks confirmed the results of combined GWAS and QTL mapping. One SNP marker located on the 2B chromosome (S2B_26494801) was found to be associated with recovery traits in both populations. The results of this study provided new information on understanding and improving drought tolerance in winter wheat.

Genome-wide association mapping and genomic prediction for kernel color traits in intermediate wheatgrass (Thinopyrum intermedium)

BACKGROUND

Intermediate wheatgrass (IWG) is a novel perennial grain crop currently undergoing domestication. It offers important ecosystem benefits while producing grain suitable for human consumption. Several aspects of plant biology and genetic control are yet to be studied in this new crop. To understand trait behavior and genetic characterization of kernel color in IWG breeding germplasm from the University of Minnesota was evaluated for the CIELAB components (L*, a*, b*) and visual differences. Trait values were used in a genome-wide association scan to reveal genomic regions controlling IWG's kernel color. The usability of genomic prediction in predicting kernel color traits was also evaluated using a four-fold cross validation method.

RESULTS

A wide phenotypic variation was observed for all four kernel color traits with pairwise trait correlations ranging from - 0.85 to 0.27. Medium to high estimates of broad sense trait heritabilities were observed and ranged from 0.41 to 0.78. A genome-wide association scan with single SNP markers detected 20 significant marker-trait associations in 9 chromosomes and 23 associations in 10 chromosomes using multi-allelic haplotype blocks. Four of the 20 significant SNP markers and six of the 23 significant haplotype blocks were common between two or more traits. Evaluation of genomic prediction of kernel color traits revealed the visual score to have highest mean predictive ability (r² = 0.53); r² for the CIELAB traits ranged from 0.29-0.33. A search for candidate genes led to detection of seven IWG genes in strong alignment with MYB36 transcription factors from other cereal crops of the Triticeae tribe. Three of these seven IWG genes had moderate similarities with R-A1, R-B1, and R-D1, the three genes that control grain color in wheat.

CONCLUSIONS

We characterized the distribution of kernel color in IWG for the first time, which revealed a broad phenotypic diversity in an elite breeding germplasm. Identification of genetic loci controlling the trait and a proof-of-concept that genomic selection might be useful in selecting genotypes of interest could help accelerate the breeding of this novel crop towards specific end-use.

Identification and Validation of Quantitative Trait Loci for Wheat Dwarf Virus Resistance in Wheat (Triticum spp.)

Wheat dwarf virus (WDV) is transmitted by the leafhopper Psammotettix alienus. As a major pathogen in wheat and other cereals, WDV causes high yield losses in many European countries. Due to climate change, insect-transmitted viruses will become more important and the restrictions in the use of insecticides efficient against P. alienus renders growing of WDV resistant/tolerant varieties the only effective strategy to control WDV. So far, there is little information about the possible sources of resistance and no known information about the genome regions responsible for the resistance. In a screening for WDV resistance using artificial inoculation in gauze houses, a panel of 500 wheat accessions including cultivars, gene bank accessions, and wild relatives of wheat was phenotyped for virus titer, infection rate, as well as plant height and yield parameters relative to healthy controls of the same genotype. Additionally, 85 T. aestivum-Ae. tauschii intogression lines were tested for WDV resistance in the greenhouse. A subset of 250 hexaploid wheat accessions was genotyped with the 15k iSelect SNP Chip. By genome-wide association study (GWAS), the quantitative trait loci (QTL) for partial WDV resistance were identified. Within these studies, one cultivar was identified showing an average infection rate of only 5.7%. By analyzing single seed descent (SSD) and doubled haploid (DH) populations comprising 153 and 314 individuals for WDV resistance and by genotyping these with the 25k iSelect SNP Chip, QTL for yield per plant, thousand-grain weight, and relative virus titer were validated on chromosomes 1B, 2B, 3B, 4B, 4A, 5A, 6A, and 7A. These results will be the basis for marker-assisted selection for WDV resistance to replacing the laborious, time-consuming, and technically challenging phenotyping with WDV bearing leafhoppers.

Discovery of Resistance Genes in Rye by Targeted Long-Read Sequencing and Association Genetics

The majority of released rye cultivars are susceptible to leaf rust because of a low level of resistance in the predominant hybrid rye-breeding gene pools Petkus and Carsten. To discover new sources of leaf rust resistance, we phenotyped a diverse panel of inbred lines from the less prevalent Gülzow germplasm using six distinct isolates of Puccinia recondita f. sp. secalis and found that 55 out of 92 lines were resistant to all isolates. By performing a genome-wide association study using 261,406 informative SNP markers, we identified five resistance-associated QTLs on chromosome arms 1RS, 1RL, 2RL, 5RL and 7RS. To identify candidate Puccinia recondita (Pr) resistance genes in these QTLs, we sequenced the rye nucleotide-binding leucine-rich repeat (NLR) intracellular immune receptor complement using a Triticeae NLR bait-library and PacBio^® long-read single-molecule high-fidelity (HiFi) sequencing. Trait-genotype correlations across 10 resistant and 10 susceptible lines identified four candidate NLR-encoding Pr genes. One of these physically co-localized with molecular markers delimiting Pr3 on chromosome arm 1RS and the top-most resistance-associated QTL in the panel.

Association mapping of autumn-seeded rye (Secale cereale L.) reveals genetic linkages between genes controlling winter hardiness and plant development

Winter field survival (WFS) in autumn-seeded winter cereals is a complex trait associated with low temperature tolerance (LTT), prostrate growth habit (PGH), and final leaf number (FLN). WFS and the three sub-traits were analyzed by a genome-wide association study of 96 rye (Secale cereal L.) genotypes of different origins and winter-hardiness levels. A total of 10,244 single nucleotide polymorphism (SNP) markers were identified by genotyping by sequencing and 259 marker-trait-associations (MTAs; p < 0.01) were revealed by association mapping. The ten most significant SNPs (p < 1.49e−04) associated with WFS corresponded to nine strong candidate genes: Inducer of CBF Expression 1 (ICE1), Cold-regulated 413-Plasma Membrane Protein 1 (COR413-PM1), Ice Recrystallization Inhibition Protein 1 (IRIP1), Jasmonate-resistant 1 (JAR1), BIPP2C1-like protein phosphatase, Chloroplast Unusual Positioning Protein-1 (CHUP1), FRIGIDA-like 4 (FRL4-like) protein, Chalcone Synthase 2 (CHS2), and Phenylalanine Ammonia-lyase 8 (PAL8). Seven of the candidate genes were also significant for one or several of the sub-traits supporting the hypothesis that WFS, LTT, FLN, and PGH are genetically interlinked. The winter-hardy rye genotypes generally carried additional allele variants for the strong candidate genes, which suggested allele diversity was a major contributor to cold acclimation efficiency and consistent high WFS under varying field conditions.