Correction: Single Marker and Haplotype-Based Association Analysis of Semolina and Pasta Colour in Elite Durum Wheat Breeding Lines Using a High-Density Consensus Map

Integrated Assays of Genome-Wide Association Study, Multi-Omics Co-Localization, and Machine Learning Associated Calcium Signaling Genes with Oilseed Rape Resistance to Sclerotinia sclerotiorum

Sclerotinia sclerotiorum (Ss) is one of the most devastating fungal pathogens, causing huge yield loss in multiple economically important crops including oilseed rape. Plant resistance to Ss pertains to quantitative disease resistance (QDR) controlled by multiple minor genes. Genome-wide identification of genes involved in QDR to Ss is yet to be conducted. In this study, we integrated several assays including genome-wide association study (GWAS), multi-omics co-localization, and machine learning prediction to identify, on a genome-wide scale, genes involved in the oilseed rape QDR to Ss. Employing GWAS and multi-omics co-localization, we identified seven resistance-associated loci (RALs) associated with oilseed rape resistance to Ss. Furthermore, we developed a machine learning algorithm and named it Integrative Multi-Omics Analysis and Machine Learning for Target Gene Prediction (iMAP), which integrates multi-omics data to rapidly predict disease resistance-related genes within a broad chromosomal region. Through iMAP based on the identified RALs, we revealed multiple calcium signaling genes related to the QDR to Ss. Population-level analysis of selective sweeps and haplotypes of variants confirmed the positive selection of the predicted calcium signaling genes during evolution. Overall, this study has developed an algorithm that integrates multi-omics data and machine learning methods, providing a powerful tool for predicting target genes associated with specific traits. Furthermore, it makes a basis for further understanding the role and mechanisms of calcium signaling genes in the QDR to Ss.

Identification of superior and rare haplotypes to optimize branch number in soybean

KEY MESSAGE

Using the integrated approach in the present study, we identified eleven significant SNPs, seven stable QTLs and 20 candidate genes associated with branch number in soybean. Branch number is a key yield-related quantitative trait that directly affects the number of pods and seeds per soybean plant. In this study, an integrated approach with a genome-wide association study (GWAS) and haplotype and candidate gene analyses was used to determine the detailed genetic basis of branch number across a diverse set of soybean accessions. The GWAS revealed a total of eleven SNPs significantly associated with branch number across three environments using the five GWAS models. Based on the consistency of the SNP detection in multiple GWAS models and environments, seven genomic regions within the physical distance of ± 202.4 kb were delineated as stable QTLs. Of these QTLs, six QTLs were novel, viz., qBN7, qBN13, qBN16, qBN18, qBN19 and qBN20, whereas the remaining one, viz., qBN12, has been previously reported. Moreover, 11 haplotype blocks, viz., Hap4, Hap7, Hap12, Hap13A, Hap13B, Hap16, Hap17, Hap18, Hap19A, Hap19B and Hap20, were identified on nine different chromosomes. Haplotype allele number across the identified haplotype blocks varies from two to five, and different branch number phenotype is regulated by these alleles ranging from the lowest to highest through intermediate branching. Furthermore, 20 genes were identified underlying the genomic region of ± 202.4 kb of the identified SNPs as putative candidates; and six of them showed significant differential expression patterns among the soybean cultivars possessing contrasting branch number, which might be the potential candidates regulating branch number in soybean. The findings of this study can assist the soybean breeding programs for developing cultivars with desirable branch numbers.

Multi-model genome-wide association studies for appearance quality in rice

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

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

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

Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat

Terminal drought stress poses a big challenge to sustain wheat grain production in rain-fed environments. This study aimed to utilize the genetically diverse pre-breeding lines for identification of genomic regions associated with agro-physiological traits at terminal stage drought stress in wheat. A total of 339 pre-breeding lines panel derived from three-way crosses of 'exotics × elite × elite' lines were evaluated in field conditions at Obregon, Mexico for two years under well irrigated as well as drought stress environments. Drought stress was imposed at flowering by skipping the irrigations at pre and post anthesis stage. Results revealed that drought significantly reduced grain yield (Y), spike length (SL), number of grains spikes-1 (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased. Population structure analysis in this panel uncovered three sub-populations. Genome wide linkage disequilibrium (LD) decay was observed at 2.5 centimorgan (cM). The haplotypes-based genome wide association study (GWAS) identified significant associations of Y, SL, and TKW on three chromosomes; 4A (HB10.7), 2D (HB6.10) and 3B (HB8.12), respectively. Likewise, associations on chromosomes 6B (HB17.1) and 3A (HB7.11) were found for NGS while on chromosome 3A (HB7.12) for KA. The genomic analysis information generated in the study can be efficiently utilized to improve Y and/or related parameters under terminal stage drought stress through marker-assisted breeding.