Genomic selection to improve husk tightness based on genomic molecular markers in maize

Introduction

The husk tightness (HTI) in maize plays a crucial role in regulating the water content of ears during the maturity stage, thereby influencing the quality of mechanical grain harvesting in China. Genomic selection (GS), which employs molecular markers, offers a promising approach for identifying and selecting inbred lines with the desired HTI trait in maize breeding. However, the effectiveness of GS is contingent upon various factors, including the genetic architecture of breeding populations, sequencing platforms, and statistical models.

Methods

An association panel of maize inbred lines was grown across three sites over two years, divided into four subgroups. GS analysis for HTI prediction was performed using marker data from three sequencing platforms and six marker densities with six statistical methods.

Results

The findings indicate that a loosely attached husk can aid in the dissipation of water from kernels in temperate maize germplasms across most environments but not nessarily for tropical-origin maize. Considering the balance between GS prediction accuracy and breeding cost, the optimal prediction strategy is the rrBLUP model, the 50K sequencing platform, a 30% proportion of the test population, and a marker density of r2=0.1. Additionally, selecting a specific SS subgroup for sampling the testing set significantly enhances the predictive capacity for husk tightness.

Discussion

The determination of the optimal GS prediction strategy for HTI provides an economically feasible reference for the practice of molecular breeding. It also serves as a reference method for GS breeding of other agronomic traits.

Alternative scoring methods of fusarium head blight resistance for genomic assisted breeding

Fusarium head blight (FHB) is a fungal disease of wheat (Triticum aestivum.L) that causes yield losses and produces mycotoxins which could easily exceed the limits of the EU regulations. Resistance to FHB has a complex genetic architecture and accurate evaluation in breeding programs is key to selecting resistant varieties. The Area Under the Disease Progress Curve (AUDPC) is one of the commonly metric used as a standard methodology to score FHB. Although efficient, AUDPC requires significant costs in phenotyping to cover the entire disease development pattern. Here, we show that there are more efficient alternatives to AUDPC (angle, growing degree days to reach 50% FHB severity, and FHB maximum variance) that reduce the number of field assessments required and allow for fair comparisons between unbalanced evaluations across trials. Furthermore, we found that the evaluation method that captures the maximum variance in FHB severity across plots is the most optimal approach for scoring FHB. In addition, results obtained on experimental data were validated on a simulated experiment where the disease progress curve was modeled as a sigmoid curve with known parameters and assessment protocols were fully controlled. Results show that alternative metrics tested in this study captured key components of quantitative plant resistance. Moreover, the new metrics could be a starting point for more accurate methods for measuring FHB in the field. For example, the optimal interval for FHB evaluation could be predicted using prior knowledge from historical weather data and FHB scores from previous trials. Finally, the evaluation methods presented in this study can reduce the FHB phenotyping burden in plant breeding with minimal losses on signal detection, resulting in a response variable available to use in data-driven analysis such as genome-wide association studies or genomic selection.

Selection strategies to introgress water deficit tolerance derived from Solanum galapagense accession LA1141 into cultivated tomato

Crop wild relatives have been used as a source of genetic diversity for over one hundred years. The wild tomato relative Solanum galapagense accession LA1141 demonstrates the ability to tolerate deficit irrigation, making it a potential resource for crop improvement. Accessing traits from LA1141 through introgression may improve the response of cultivated tomatoes grown in water-limited environments. Canopy temperature is a proxy for physiological traits which are challenging to measure efficiently and may be related to water deficit tolerance. We optimized phenotypic evaluation based on variance partitioning and further show that objective phenotyping methods coupled with genomic prediction lead to gain under selection for water deficit tolerance. The objectives of this work were to improve phenotyping workflows for measuring canopy temperature, mapping quantitative trait loci (QTLs) from LA1141 that contribute to water deficit tolerance and comparing selection strategies. The phenotypic variance attributed to genetic causes for canopy temperature was higher when estimated from thermal images relative to estimates based on an infrared thermometer. Composite interval mapping using BC₂S₃ families, genotyped with single nucleotide polymorphisms, suggested that accession LA1141 contributed alleles that lower canopy temperature and increase plant turgor under water deficit. QTLs for lower canopy temperature were mapped to chromosomes 1 and 6 and explained between 6.6 and 9.5% of the total phenotypic variance. QTLs for higher leaf turgor were detected on chromosomes 5 and 7 and explained between 6.8 and 9.1% of the variance. We advanced tolerant BC₂S₃ families to the BC₂S₅ generation using selection indices based on phenotypic values and genomic estimated breeding values (GEBVs). Phenotypic, genomic, and combined selection strategies demonstrated gain under selection and improved performance compared to randomly advanced BC₂S₅ progenies. Leaf turgor, canopy temperature, stomatal conductance, and vapor pressure deficit (VPD) were evaluated and compared in BC₂S₅ progenies grown under deficit irrigation. Progenies co-selected for phenotypic values and GEBVs wilted less, had significantly lower canopy temperature, higher stomatal conductance, and lower VPD than randomly advanced lines. The fruit size of water deficit tolerant selections was small compared to the recurrent parent. However, lines with acceptable yield, canopy width, and quality parameters were recovered. These results suggest that we can create selection indices to improve water deficit tolerance in a recurrent parent background, and additional crossing and evaluation are warranted.

Cotton Breeding in Australia: Meeting the Challenges of the 21st Century

The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program is the sole breeding effort for cotton in Australia, developing high performing cultivars for the local industry which is worth∼AU$3 billion per annum. The program is supported by Cotton Breeding Australia, a Joint Venture between CSIRO and the program's commercial partner, Cotton Seed Distributors Ltd. (CSD). While the Australian industry is the focus, CSIRO cultivars have global impact in North America, South America, and Europe. The program is unique compared with many other public and commercial breeding programs because it focuses on diverse and integrated research with commercial outcomes. It represents the full research pipeline, supporting extensive long-term fundamental molecular research; native and genetically modified (GM) trait development; germplasm enhancement focused on yield and fiber quality improvements; integration of third-party GM traits; all culminating in the release of new commercial cultivars. This review presents evidence of past breeding successes and outlines current breeding efforts, in the areas of yield and fiber quality improvement, as well as the development of germplasm that is resistant to pests, diseases and abiotic stressors. The success of the program is based on the development of superior germplasm largely through field phenotyping, together with strong commercial partnerships with CSD and Bayer CropScience. These relationships assist in having a shared focus and ensuring commercial impact is maintained, while also providing access to markets, traits, and technology. The historical successes, current foci and future requirements of the CSIRO cotton breeding program have been used to develop a framework designed to augment our breeding system for the future. This will focus on utilizing emerging technologies from the genome to phenome, as well as a panomics approach with data management and integration to develop, test and incorporate new technologies into a breeding program. In addition to streamlining the breeding pipeline for increased genetic gain, this technology will increase the speed of trait and marker identification for use in genome editing, genomic selection and molecular assisted breeding, ultimately producing novel germplasm that will meet the coming challenges of the 21st Century.

Genomic Selection for Wheat Blast in a Diversity Panel, Breeding Panel and Full-Sibs Panel

Wheat blast is an emerging threat to wheat production, due to its recent migration to South Asia and Sub-Saharan Africa. Because genomic selection (GS) has emerged as a promising breeding strategy, the key objective of this study was to evaluate it for wheat blast phenotyped at precision phenotyping platforms in Quirusillas (Bolivia), Okinawa (Bolivia) and Jashore (Bangladesh) using three panels: (i) a diversity panel comprising 172 diverse spring wheat genotypes, (ii) a breeding panel comprising 248 elite breeding lines, and (iii) a full-sibs panel comprising 298 full-sibs. We evaluated two genomic prediction models (the genomic best linear unbiased prediction or GBLUP model and the Bayes B model) and compared the genomic prediction accuracies with accuracies from a fixed effects model (with selected blast-associated markers as fixed effects), a GBLUP + fixed effects model and a pedigree relationships-based model (ABLUP). On average, across all the panels and environments analyzed, the GBLUP + fixed effects model (0.63 ± 0.13) and the fixed effects model (0.62 ± 0.13) gave the highest prediction accuracies, followed by the Bayes B (0.59 ± 0.11), GBLUP (0.55 ± 0.1), and ABLUP (0.48 ± 0.06) models. The high prediction accuracies from the fixed effects model resulted from the markers tagging the 2NS translocation that had a large effect on blast in all the panels. This implies that in environments where the 2NS translocation-based blast resistance is effective, genotyping one to few markers tagging the translocation is sufficient to predict the blast response and genome-wide markers may not be needed. We also observed that marker-assisted selection (MAS) based on a few blast-associated markers outperformed GS as it selected the highest mean percentage (88.5%) of lines also selected by phenotypic selection and discarded the highest mean percentage of lines (91.8%) also discarded by phenotypic selection, across all panels. In conclusion, while this study demonstrates that MAS might be a powerful strategy to select for the 2NS translocation-based blast resistance, we emphasize that further efforts to use genomic tools to identify non-2NS translocation-based blast resistance are critical.

Increased Predictive Accuracy of Multi-Environment Genomic Prediction Model for Yield and Related Traits in Spring Wheat (Triticum aestivum L.)

Genomic selection (GS) has the potential to improve the selection gain for complex traits in crop breeding programs from resource-poor countries. The GS model performance in multi-environment (ME) trials was assessed for 141 advanced breeding lines under four field environments via cross-predictions. We compared prediction accuracy (PA) of two GS models with or without accounting for the environmental variation on four quantitative traits of significant importance, i.e., grain yield (GRYLD), thousand-grain weight, days to heading, and days to maturity, under North and Central Indian conditions. For each trait, we generated PA using the following two different ME cross-validation (CV) schemes representing actual breeding scenarios: (1) predicting untested lines in tested environments through the ME model (ME_CV1) and (2) predicting tested lines in untested environments through the ME model (ME_CV2). The ME predictions were compared with the baseline single-environment (SE) GS model (SE_CV1) representing a breeding scenario, where relationships and interactions are not leveraged across environments. Our results suggested that the ME models provide a clear advantage over SE models in terms of robust trait predictions. Both ME models provided 2-3 times higher prediction accuracies for all four traits across the four tested environments, highlighting the importance of accounting environmental variance in GS models. While the improvement in PA from SE to ME models was significant, the CV1 and CV2 schemes did not show any clear differences within ME, indicating the ME model was able to predict the untested environments and lines equally well. Overall, our results provide an important insight into the impact of environmental variation on GS in smaller breeding programs where these programs can potentially increase the rate of genetic gain by leveraging the ME wheat breeding trials.

An Assessment of the Factors Influencing the Prediction Accuracy of Genomic Prediction Models Across Multiple Environments

The effects of climate change create formidable challenges for breeders striving to produce sufficient food quantities in rapidly changing environments. It is therefore critical to investigate the ability of multi-environment genomic prediction (GP) models to predict genomic estimated breeding values (GEBVs) in extreme environments. Exploration of the impact of training set composition on the accuracy of such GEBVs is also essential. Accordingly, we examined the influence of the number of training environments and the use of environmental covariates (ECs) in GS models on four subsets of n = 500 lines of the soybean nested association mapping (SoyNAM) panel grown in nine environments in the US-North Central Region. The ensuing analyses provided insights into the influence of both of these factors for predicting grain yield in the most and the least extreme of these environments. We found that only a subset of the available environments was needed to obtain the highest observed prediction accuracies. The inclusion of ECs in the GP model did not substantially increase prediction accuracies relative to competing models, and instead more often resulted in negative prediction accuracies. Combined with the overall low prediction accuracies for grain yield in the most extreme environment, our findings highlight weaknesses in current GP approaches for prediction in extreme environments, and point to specific areas on which to focus future research efforts.

GWAS-Assisted Genomic Prediction to Predict Resistance to Septoria Tritici Blotch in Nordic Winter Wheat at Seedling Stage

Septoria tritici blotch (STB) disease caused by Zymoseptoria tritici is one of the most damaging diseases of wheat causing significant yield losses worldwide. Identification and employment of resistant germplasm is the most cost-effective method to control STB. In this study, we characterized seedling stage resistance to STB in 175 winter wheat landraces and old cultivars of Nordic origin. The study revealed significant (p < 0.05) phenotypic differences in STB severity in the germplasm. Genome-wide association analysis (GWAS) using five different algorithms identified ten significant markers on five chromosomes. Six markers were localized within a region of 2 cM that contained seven candidate genes on chromosome 1B. Genomic prediction (GP) analysis resulted in a model with an accuracy of 0.47. To further improve the prediction efficiency, significant markers identified by GWAS were included as fixed effects in the GP model. Depending on the number of fixed effect markers, the prediction accuracy improved from 0.47 (without fixed effects) to 0.62 (all non-redundant GWAS markers as fixed effects), respectively. The resistant genotypes and single-nucleotide polymorphism (SNP) markers identified in the present study will serve as a valuable resource for future breeding for STB resistance in wheat. The results also highlight the benefits of integrating GWAS with GP to further improve the accuracy of GP.

Genomic selection on shelling percentage and other traits for maize

Genomic selection (GS) is the one of the new method for molecular marker-assisted selection (MAS) that can improve selection efficiency and thereby accelerate selective breeding progress. In the present study, we used the exotic germplasm LK1 to improve the shelling percentage of Qi319 by GS. Genome-wide marker effects for each trait were estimated based on the performance of the testcross and SNP data for F₂ progenies in the training population. The accuracy of genomic predictions was estimated as the correlation between marker-predicted genotypic values and phenotypic values of the testcrosses for each trait in the validation population. Our study result indicated that selection response for shell percentage was 33.7%, which is greater than those for grain yield, kernel number per ear, or grain moisture at harvest. Selection response for tassel branch number and weight per 100 kernels was greater than 60%. The Higher trait heritability resulted in better prediction efficiency; Prediction accuracy increased with the training population size; Prediction efficiency did not differ significantly between SNP densities of 1000 bp and 55,000 bp. The results of the present research project will provide a basis for genome-wide selection technology in maize breeding, and lay the groundwork for the application of GS to germplasms that are useful in China.

Genetic Correlation, Genome-Wide Association and Genomic Prediction of Portable NIRS Predicted Carotenoids in Cassava Roots

Random forests (RF) was used to correlate spectral responses to known wet chemistry carotenoid concentrations including total carotenoid content (TCC), all-trans β-carotene (ATBC), violaxanthin (VIO), lutein (LUT), 15-cis beta-carotene (15CBC), 13-cis beta-carotene (13CBC), alpha-carotene (AC), 9-cis beta-carotene (9CBC), and phytoene (PHY) from laboratory analysis of 173 cassava root samples in Columbia. The cross-validated correlations between the actual and estimated carotenoid values using RF ranged from 0.62 in PHY to 0.97 in ATBC. The developed models were used to evaluate the carotenoids of 594 cassava clones with spectral information collected across three locations in a national breeding program (NRCRI, Umudike), Nigeria. Both populations contained cassava clones characterized as white and yellow. The NRCRI evaluated phenotypes were used to assess the genetic correlations, conduct genome-wide association studies (GWAS), and genomic predictions. Estimates of genetic correlation showed various levels of the relationship among the carotenoids. The associations between TCC and the individual carotenoids were all significant (P < 0.001) with high positive values (r > 0.75, except in LUT and PHY where r < 0.3). The GWAS revealed significant genomic regions on chromosomes 1, 2, 4, 13, 14, and 15 associated with variation in at least one of the carotenoids. One of the identified candidate genes, phytoene synthase (PSY) has been widely reported for variation in TCC in cassava. On average, genomic prediction accuracies from the single-trait genomic best linear unbiased prediction (GBLUP) and RF as well as from a multiple-trait GBLUP model ranged from ∼0.2 in LUT and PHY to 0.52 in TCC. The multiple-trait GBLUP model gave slightly higher accuracies than the single trait GBLUP and RF models. This study is one of the initial attempts in understanding the genetic basis of individual carotenoids and demonstrates the usefulness of NIRS in cassava improvement.

Quantitative Genetics and Genomics Converge to Accelerate Forest Tree Breeding

Forest tree breeding has been successful at delivering genetically improved material for multiple traits based on recurrent cycles of selection, mating, and testing. However, long breeding cycles, late flowering, variable juvenile-mature correlations, emerging pests and diseases, climate, and market changes, all pose formidable challenges. Genetic dissection approaches such as quantitative trait mapping and association genetics have been fruitless to effectively drive operational marker-assisted selection (MAS) in forest trees, largely because of the complex multifactorial inheritance of most, if not all traits of interest. The convergence of high-throughput genomics and quantitative genetics has established two new paradigms that are changing contemporary tree breeding dogmas. Genomic selection (GS) uses large number of genome-wide markers to predict complex phenotypes. It has the potential to accelerate breeding cycles, increase selection intensity and improve the accuracy of breeding values. Realized genomic relationships matrices, on the other hand, provide innovations in genetic parameters' estimation and breeding approaches by tracking the variation arising from random Mendelian segregation in pedigrees. In light of a recent flow of promising experimental results, here we briefly review the main concepts, analytical tools and remaining challenges that currently underlie the application of genomics data to tree breeding. With easy and cost-effective genotyping, we are now at the brink of extensive adoption of GS in tree breeding. Areas for future GS research include optimizing strategies for updating prediction models, adding validated functional genomics data to improve prediction accuracy, and integrating genomic and multi-environment data for forecasting the performance of genetic material in untested sites or under changing climate scenarios. The buildup of phenotypic and genome-wide data across large-scale breeding populations and advances in computational prediction of discrete genomic features should also provide opportunities to enhance the application of genomics to tree breeding.

Genetics of Resistance and Pathogenicity in the Maize/Setosphaeria turcica Pathosystem and Implications for Breeding

Northern corn leaf blight (NCLB), the most devastating leaf pathogen in maize (Zea mays L.), is caused by the heterothallic ascomycete Setosphaeria turcica. The pathogen population shows an extremely high genetic diversity in tropical and subtropical regions. Varietal resistance is the most efficient technique to control NCLB. Host resistance can be qualitative based on race-specific Ht genes or quantitative controlled by many genes with small effects. Quantitative resistance is moderately to highly effective and should be more durable combatting all races of the pathogen. Quantitative resistance must, however, be analyzed in many environments (= location × year combinations) to select stable resistances. In the tropical and subtropical environments, quantitative resistance is the preferred option to manage NCLB epidemics. Resistance level can be increased in practical breeding programs by several recurrent selection cycles based on disease severity rating and/or by genomic selection. This review aims to address two important aspects of the NCLB pathosystem: the genetics of the fungus S. turcica and the modes of inheritance of the host plant maize, including successful breeding strategies regarding NCLB resistance. Both drivers of this pathosystem, pathogen, and host, must be taken into account to result in more durable resistance.

Genomics-assisted breeding in fruit trees

Recent advancements in genomic analysis technologies have opened up new avenues to promote the efficiency of plant breeding. Novel genomics-based approaches for plant breeding and genetics research, such as genome-wide association studies (GWAS) and genomic selection (GS), are useful, especially in fruit tree breeding. The breeding of fruit trees is hindered by their long generation time, large plant size, long juvenile phase, and the necessity to wait for the physiological maturity of the plant to assess the marketable product (fruit). In this article, we describe the potential of genomics-assisted breeding, which uses these novel genomics-based approaches, to break through these barriers in conventional fruit tree breeding. We first introduce the molecular marker systems and whole-genome sequence data that are available for fruit tree breeding. Next we introduce the statistical methods for biparental linkage and quantitative trait locus (QTL) mapping as well as GWAS and GS. We then review QTL mapping, GWAS, and GS studies conducted on fruit trees. We also review novel technologies for rapid generation advancement. Finally, we note the future prospects of genomics-assisted fruit tree breeding and problems that need to be overcome in the breeding.

Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding

Marker-assisted selection (MAS) refers to the use of molecular markers to assist phenotypic selections in crop improvement. Several types of molecular markers, such as single nucleotide polymorphism (SNP), have been identified and effectively used in plant breeding. The application of next-generation sequencing (NGS) technologies has led to remarkable advances in whole genome sequencing, which provides ultra-throughput sequences to revolutionize plant genotyping and breeding. To further broaden NGS usages to large crop genomes such as maize and wheat, genotyping-by-sequencing (GBS) has been developed and applied in sequencing multiplexed samples that combine molecular marker discovery and genotyping. GBS is a novel application of NGS protocols for discovering and genotyping SNPs in crop genomes and populations. The GBS approach includes the digestion of genomic DNA with restriction enzymes followed by the ligation of barcode adapter, PCR amplification and sequencing of the amplified DNA pool on a single lane of flow cells. Bioinformatic pipelines are needed to analyze and interpret GBS datasets. As an ultimate MAS tool and a cost-effective technique, GBS has been successfully used in implementing genome-wide association study (GWAS), genomic diversity study, genetic linkage analysis, molecular marker discovery and genomic selection under a large scale of plant breeding programs.

Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding

Marker-assisted selection (MAS) refers to the use of molecular markers to assist phenotypic selections in crop improvement. Several types of molecular markers, such as single nucleotide polymorphism (SNP), have been identified and effectively used in plant breeding. The application of next-generation sequencing (NGS) technologies has led to remarkable advances in whole genome sequencing, which provides ultra-throughput sequences to revolutionize plant genotyping and breeding. To further broaden NGS usages to large crop genomes such as maize and wheat, genotyping-by-sequencing (GBS) has been developed and applied in sequencing multiplexed samples that combine molecular marker discovery and genotyping. GBS is a novel application of NGS protocols for discovering and genotyping SNPs in crop genomes and populations. The GBS approach includes the digestion of genomic DNA with restriction enzymes followed by the ligation of barcode adapter, PCR amplification and sequencing of the amplified DNA pool on a single lane of flow cells. Bioinformatic pipelines are needed to analyze and interpret GBS datasets. As an ultimate MAS tool and a cost-effective technique, GBS has been successfully used in implementing genome-wide association study (GWAS), genomic diversity study, genetic linkage analysis, molecular marker discovery and genomic selection under a large scale of plant breeding programs.

Potential assessment of genome-wide association study and genomic selection in Japanese pear Pyrus pyrifolia

Although the potential of marker-assisted selection (MAS) in fruit tree breeding has been reported, bi-parental QTL mapping before MAS has hindered the introduction of MAS to fruit tree breeding programs. Genome-wide association studies (GWAS) are an alternative to bi-parental QTL mapping in long-lived perennials. Selection based on genomic predictions of breeding values (genomic selection: GS) is another alternative for MAS. This study examined the potential of GWAS and GS in pear breeding with 76 Japanese pear cultivars to detect significant associations of 162 markers with nine agronomic traits. We applied multilocus Bayesian models accounting for ordinal categorical phenotypes for GWAS and GS model training. Significant associations were detected at harvest time, black spot resistance and the number of spurs and two of the associations were closely linked to known loci. Genome-wide predictions for GS were accurate at the highest level (0.75) in harvest time, at medium levels (0.38-0.61) in resistance to black spot, firmness of flesh, fruit shape in longitudinal section, fruit size, acid content and number of spurs and at low levels (<0.2) in all soluble solid content and vigor of tree. Results suggest the potential of GWAS and GS for use in future breeding programs in Japanese pear.