Genomic Prediction for Winter Survival of Lowland Switchgrass in the Northern USA

Mapping quantitative trait loci for biomass yield and yield-related traits in lowland switchgrass (Panicum virgatum L.) multiple populations

Switchgrass can be used as an alternative source for bioenergy production. Many breeding programs focus on the genetic improvement of switchgrass for increasing biomass yield. Quantitative trait loci (QTL) mapping can help to discover marker-trait associations and accelerate the breeding process through marker-assisted selection. To identify significant QTL, this study mapped 7 hybrid populations and one combined of 2 hybrid populations (30-96 F1s) derived from Alamo and Kanlow genotypes. The populations were evaluated for biomass yield, plant height, and crown size in a simulated-sward plot with 2 replications at 2 locations in Tennessee from 2019 to 2021. The populations showed significant genetic variation for the evaluated traits and exhibited transgressive segregation. The 17,251 single nucleotide polymorphisms (SNPs) generated through genotyping-by-sequencing (GBS) were used to construct a linkage map using a fast algorithm for multiple outbred families. The linkage map spanned 1,941 cM with an average interval of 0.11 cM between SNPs. The QTL analysis was performed on evaluated traits for each and across environments (year and location) that identified 5 QTL for biomass yield (logarithm of the odds, LOD 3.12-4.34), 4 QTL for plant height (LOD 3.01-5.64), and 7 QTL for crown size (LOD 3.0-4.46) (P ≤ 0.05). The major QTL for biomass yield, plant height, and crown size resided on chromosomes 8N, 6N, and 8K explained phenotypic variations of 5.6, 5.1, and 6.6%, respectively. SNPs linked to QTL could be incorporated into marker-assisted breeding to maximize the selection gain in switchgrass breeding.

Divergent Metabolic Changes in Rhizomes of Lowland and Upland Switchgrass (Panicum virgatum) from Early Season through Dormancy Onset

High-biomass-yielding southerly adapted switchgrasses (Panicum virgatum L.) frequently suffer from unpredictable winter hardiness at more northerly sites arising from damage to rhizomes that prevent effective spring regrowth. Previously, changes occurring over the growing season in rhizomes sampled from a cold-adapted tetraploid upland cultivar, Summer, demonstrated a role for abscisic acid (ABA), starch accumulation, and transcriptional reprogramming as drivers of dormancy onset and potential keys to rhizome health during winter dormancy. Here, rhizome metabolism of a high-yielding southerly adapted tetraploid switchgrass cultivar, Kanlow-which is a significant source of genetics for yield improvement-was studied over a growing season at a northern site. Metabolite levels and transcript abundances were combined to develop physiological profiles accompanying greening through the onset of dormancy in Kanlow rhizomes. Next, comparisons of the data to rhizome metabolism occurring in the adapted upland cultivar Summer were performed. These data revealed both similarities as well as numerous differences in rhizome metabolism that were indicative of physiological adaptations unique to each cultivar. Similarities included elevated ABA levels and accumulation of starch in rhizomes during dormancy onset. Notable differences were observed in the accumulation of specific metabolites, the expression of genes encoding transcription factors, and several enzymes linked to primary metabolism.

Genomic prediction of switchgrass winter survivorship across diverse lowland populations

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

Mapping and selection of downy mildew resistance in spinach cv. whale by low coverage whole genome sequencing

Spinach (Spinacia oleracea) is a popular leafy vegetable crop and commercial production is centered in California and Arizona in the US. The oomycete Peronospora effusa causes the most important disease in spinach, downy mildew. A total of nineteen races of P. effusa are known, with more than 15 documented in the last three decades, and the regular emergence of new races is continually overcoming the genetic resistance to the pathogen. This study aimed to finely map the downy mildew resistance locus RPF3 in spinach, identify single nucleotide polymorphism (SNP) markers associated with the resistance, refine the candidate genes responsible for the resistance, and evaluate the prediction performance using multiple machine learning genomic prediction (GP) methods. Segregating progeny population developed from a cross of resistant cultivar Whale and susceptible cultivar Viroflay to race 5 of P. effusa was inoculated under greenhouse conditions to determine downy mildew disease response across the panel. The progeny panel and the parents were resequenced at low coverage (1x) to identify genome wide SNP markers. Association analysis was performed using disease response phenotype data and SNP markers in TASSEL, GAPIT, and GENESIS programs and mapped the race 5 resistance loci (RPF3) to 1.25 and 2.73 Mb of Monoe-Viroflay chromosome 3 with the associated SNP in the 1.25 Mb region was 0.9 Kb from the NBS-LRR gene SOV3g001250. The RPF3 locus in the 1.22-1.23 Mb region of Sp75 chromosome 3 is 2.41-3.65 Kb from the gene Spo12821 annotated as NBS-LRR disease resistance protein. This study extended our understanding of the genetic basis of downy mildew resistance in spinach cultivar Whale and mapped the RPF3 resistance loci close to the NBS-LRR gene providing a target to pursue functional validation. Three SNP markers efficiently selected resistance based on multiple genomic selection (GS) models. The results from this study have added new genomic resources, generated an informed basis of the RPF3 locus resistant to spinach downy mildew pathogen, and developed markers and prediction methods to select resistant lines.

Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods

Downy mildew, commercially the most important disease of spinach, is caused by the obligate oomycete Peronospora effusa. In the past two decades, new pathogen races have repeatedly overcome the resistance used in newly released cultivars, urging the need for more durable resistance. Commercial spinach cultivars are bred with major R genes to impart resistance to downy mildew pathogens and are effective against some pathogen races/isolates. This work aimed to evaluate the worldwide USDA spinach germplasm collections and commercial cultivars for resistance to downy mildew pathogen in the field condition under natural inoculum pressure and conduct genome wide association analysis (GWAS) to identify resistance-associated genomic regions (alleles). Another objective was to evaluate the prediction accuracy (PA) using several genomic prediction (GP) methods to assess the potential implementation of genomic selection (GS) to improve spinach breeding for resistance to downy mildew pathogen. More than four hundred diverse spinach genotypes comprising USDA germplasm accessions and commercial cultivars were evaluated for resistance to downy mildew pathogen between 2017-2019 in Salinas Valley, California and Yuma, Arizona. GWAS was performed using single nucleotide polymorphism (SNP) markers identified via whole genome resequencing (WGR) in GAPIT and TASSEL programs; detected 14, 12, 5, and 10 significantly associated SNP markers with the resistance from four tested environments, respectively; and the QTL alleles were detected at the previously reported region of chromosome 3 in three of the four experiments. In parallel, PA was assessed using six GP models and seven unique marker datasets for field resistance to downy mildew pathogen across four tested environments. The results suggest the suitability of GS to improve field resistance to downy mildew pathogen. The QTL, SNP markers, and PA estimates provide new information in spinach breeding to select resistant plants and breeding lines through marker-assisted selection (MAS) and GS, eventually helping to accumulate beneficial alleles for durable disease resistance.

One hundred years into the study of ecotypes, new advances are being made through large-scale field experiments in perennial plant systems

A hundred years after Turesson first clearly described how locally adaptive variation is distributed within species, plant biologists are making major breakthroughs in our understanding of mechanisms underlying adaptation from local populations to the scale of continents. Although the genetics of local adaptation has typically been studied in smaller reciprocal transplant experiments, it is now being evaluated with whole genomes in large-scale networks of common garden experiments with perennial switchgrass and poplar trees. These studies support the hypothesis that a complex combination of loci, both with and without adaptive trade-offs, underlies local adaptation and that hybridization and adaptive introgression play a key role in the evolution of these species. Future studies incorporating high-throughput phenotyping, gene expression, and modeling will be used to predict responses of these species to climate change.

Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm

White rust, caused by Albugo occidentalis, is one of the major yield-limiting diseases of spinach (Spinacia oleracea) in some major commercial production areas, particularly in southern Texas in the United States. The use of host resistance is the most economical and environment-friendly approach to managing white rust in spinach production. The objectives of this study were to conduct a genome-wide associating study (GWAS), to identify single nucleotide polymorphism (SNP) markers associated with white rust resistance in spinach, and to perform genomic prediction (GP) to estimate the prediction accuracy (PA). A GWAS panel of 346 USDA (US Dept. of Agriculture) germplasm accessions was phenotyped for white rust resistance under field conditions and GWAS was performed using 13 235 whole-genome resequencing (WGR) generated SNPs. Nine SNPs, chr2_53 049 132, chr3_58 479 501, chr3_95 114 909, chr4_9 176 069, chr4_17 807 168, chr4_83 938 338, chr4_87 601 768, chr6_1 877 096, and chr6_31 287 118, located on chromosomes 2, 3, 4, and 6 were associated with white rust resistance in this GWAS panel. Four scenarios were tested for PA using Pearson's correlation coefficient (r) between the genomic estimation breeding value (GEBV) and the observed values: (1) different ratios between the training set and testing set (fold), (2) different GP models, (3) different SNP numbers in three different SNP sets, and (4) the use of GWAS-derived significant SNP markers. The results indicated that a 2- to 10-fold difference in the various GP models had similar, although not identical, averaged r values in each SNP set; using GWAS-derived significant SNP markers would increase PA with a high r-value up to 0.84. The SNP markers and the high PA can provide valuable information for breeders to improve spinach by marker-assisted selection (MAS) and genomic selection (GS).

Genomic Prediction of Complex Traits in Forage Plants Species: Perennial Grasses Case

The majority of forage grass species are obligate outbreeders. Their breeding classically consists of an initial selection on spaced plants for highly heritable traits such as disease resistances and heading date, followed by familial selection on swards for forage yield and quality traits. The high level of diversity and heterozygosity, and associated decay of linkage disequilibrium (LD) over very short genomic distances, has hampered the implementation of genomic selection (GS) in these species. However, next generation sequencing technologies in combination with the development of genomic resources have recently facilitated implementation of GS in forage grass species such as perennial ryegrass (Lolium perenne L.), switchgrass (Panicum virgatum L.), and timothy (Phleum pratense L.). Experimental work and simulations have shown that GS can increase significantly the genetic gain per unit of time for traits with different levels of heritability. The main reasons are (1) the possibility to select single plants based on their genomic estimated breeding values (GEBV) for traits measured at sward level, (2) a reduction in the duration of selection cycles, and less importantly (3) an increase in the selection intensity associated with an increase in the genetic variance used for selection. Nevertheless, several factors should be taken into account for the successful implementation of GS in forage grasses. For example, it has been shown that the level of relatedness between the training and the selection population is particularly critical when working with highly structured meta-populations consisting of several genetic groups. A sufficient number of markers should be used to estimate properly the kinship between individuals and to reflect the variability of major QTLs. It is also important that the prediction models are trained for relevant environments when dealing with traits with high genotype × environment interaction (G × E). Finally, in these outbreeding species, measures to reduce inbreeding should be used to counterbalance the high selection intensity that can be achieved in GS.

Subsampling and DNA pooling can increase gains through genomic selection in switchgrass

Genomic selection (GS) can accelerate breeding cycles in perennial crops such as the bioenergy grass switchgrass (Panicum virgatum L.). The sequencing costs of GS can be reduced by pooling DNA samples in the training population (TP), only sequencing TP phenotypic outliers, or pooling candidate population (CP) samples. These strategies were simulated for two traits (spring vigor and anthesis date) in three breeding populations. Sequencing only the outlier 50% of the TP phenotype distribution resulted in a penalty of <5% of the predictive ability, measured using cross-validation. Predictive ability also decreased when sequencing progressively fewer TP DNA pools, but TPs constructed from only two phenotypically contrasting DNA samples retained a mean of >80% predictive ability relative to individual TP sequencing. Novel group testing methods allowed greater than one CP individual to be screened per sequenced DNA sample but resulted in a predictive ability penalty. To determine the impact of reduced sequencing, genetic gain was calculated for seven GS scenarios with variable sequencing budgets. Reduced TP sequencing and most CP pooling methods were superior to individual sequence-based GS when sequencing resources were restricted (2,000 DNA samples per 5-yr cycle). Only one scenario was superior to individual sequencing when sequencing budgets were large (8,000 DNA samples per 5-yr cycle). This study highlights multiple routes for reduced sequencing costs in GS.

Genetic loci associated with winter survivorship in diverse lowland switchgrass populations

High winter mortality limits biomass yield of lowland switchgrass (Panicum virgatum L.) planted in the northern latitudes of North America. Breeding of cold tolerant switchgrass cultivars requires many years due to its perennial growth habit and the unpredictable winter selection pressure that is required to identify winter-hardy individuals. Identification of causal genetic variants for winter survivorship would accelerate the improvement of switchgrass biomass production. The objective of this study was to identify allelic variation associated with winter survivorship in lowland switchgrass populations using bulk segregant analysis (BSA). Twenty-nine lowland switchgrass populations were evaluated for winter survival at two locations in southern Wisconsin and 21 populations with differential winter survivorship were used for BSA. A maximum of 10% of the individuals (8-20) were bulked to create survivor and nonsurvivor DNA pools from each population and location. The DNA pools were evaluated using exome capture sequencing, and allele frequencies were used to conduct statistical tests. The BSA tests revealed nine quatitative trait loci (QTL) from tetraploid populations and seven QTL from octoploid populations. Many QTL were population-specific, but some were identified in multiple populations that originated across a broad geographic landscape. Four QTL (at positions 88 Mb on chromosome 2N, 115 Mb on chromosome 5K, and 1 and 100 Mb on chromosome 9N) were potentially the most useful QTL. Markers associated with winter survivorship in this study can be used to accelerate breeding cycles of lowland switchgrass populations and should lead to improvements in adaptation within USDA hardiness zones 4 and 5.

Mapping freezing tolerance QTL in alfalfa: based on indoor phenotyping

BACKGROUND

Winter freezing temperature impacts alfalfa (Medicago sativa L.) persistence and seasonal yield and can lead to the death of the plant. Understanding the genetic mechanisms of alfalfa freezing tolerance (FT) using high-throughput phenotyping and genotyping is crucial to select suitable germplasm and develop winter-hardy cultivars. Several clones of an alfalfa F1 mapping population (3010 x CW 1010) were tested for FT using a cold chamber. The population was genotyped with SNP markers identified using genotyping-by-sequencing (GBS) and the quantitative trait loci (QTL) associated with FT were mapped on the parent-specific linkage maps. The ultimate goal is to develop non-dormant and winter-hardy alfalfa cultivars that can produce extended growth in the areas where winters are often mild.

RESULTS

Alfalfa FT screening method optimized in this experiment comprises three major steps: clone preparation, acclimation, and freezing test. Twenty clones of each genotype were tested, where 10 samples were treated with freezing temperature, and 10 were used as controls. A moderate positive correlation (r ~ 0.36, P < 0.01) was observed between indoor FT and field-based winter hardiness (WH), suggesting that the indoor FT test is a useful indirect selection method for winter hardiness of alfalfa germplasm. We detected a total of 20 QTL associated with four traits; nine for visual rating-based FT, five for percentage survival (PS), four for treated to control regrowth ratio (RR), and two for treated to control biomass ratio (BR). Some QTL positions overlapped with WH QTL reported previously, suggesting a genetic relationship between FT and WH. Some favorable QTL from the winter-hardy parent (3010) were from the potential genic region for a cold tolerance gene CBF. The BLAST alignment of a CBF sequence of M. truncatula, a close relative of alfalfa, against the alfalfa reference showed that the gene's ortholog resides around 75 Mb on chromosome 6.

CONCLUSIONS

The indoor freezing tolerance selection method reported is useful for alfalfa breeders to accelerate breeding cycles through indirect selection. The QTL and associated markers add to the genomic resources for the research community and can be used in marker-assisted selection (MAS) for alfalfa cold tolerance improvement.

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass

Long-term climate change and periodic environmental extremes threaten food and fuel security¹ and global crop productivity^2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience⁵, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation⁶-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

QTL mapping of winter dormancy and associated traits in two switchgrass pseudo-F1 populations: lowland x lowland and lowland x upland

BACKGROUND

Switchgrass (Panicum virgatum) undergoes winter dormancy by sensing photoperiod and temperature changes. It transitions to winter dormancy in early fall following at the end of reproduction and exits dormancy in the spring. The duration of the growing season affects the accumulation of biomass and yield. In this study, we conducted QTL mapping of winter dormancy measured by fall regrowth height (FRH) and normalized difference vegetation index (NDVI), spring emergence (SE), and flowering date (FD) in two bi-parental pseudo-F1 populations derived from crosses between the lowland AP13 with the lowland B6 (AB) with 285 progenies, and the lowland B6 with the upland VS16 (BV) with 227 progenies.

RESULTS

We identified 18 QTLs for FRH, 18 QTLs for NDVI, 21 QTLs for SE, and 30 QTLs for FD. The percent variance explained by these QTLs ranged between 4.21-23.27% for FRH, 4.47-24.06% for NDVI, 4.35-32.77% for SE, and 4.61-29.74% for FD. A higher number of QTL was discovered in the BV population, suggesting more variants in the lowland x upland population contributing to the expression of seasonal dormancy underlying traits. We identified 9 regions of colocalized QTL with possible pleiotropic gene action. The positive correlation between FRH or NDVI with dry biomass weight suggests that winter dormancy duration could affect switchgrass biomass yield. The medium to high heritability levels of FRH (0.55-0.64 H²) and NDVI (0.30-0.61 H²) implies the possibility of using the traits for indirect selection for biomass yield.

CONCLUSION

Markers found within the significant QTL interval can serve as genomic resources for breeding non-dormant and semi-dormant switchgrass cultivars for the southern regions, where growers can benefit from the longer production season.

Natural variation and genomic prediction of growth, physiological traits, and nitrogen-use efficiency in perennial ryegrass under low-nitrogen stress

Genomic prediction of nitrogen-use efficiency (NUE) has not previously been studied in perennial grass species exposed to low-N stress. Here, we conducted a genomic prediction of physiological traits and NUE in 184 global accessions of perennial ryegrass (Lolium perenne) in response to a normal (7.5 mM) and low (0.75 mM) supply of N. After 21 d of treatment under greenhouse conditions, significant variations in plant height increment (ΔHT), leaf fresh weight (LFW), leaf dry weight (LDW), chlorophyll index (Chl), chlorophyll fluorescence, leaf N and carbon (C) contents, C/N ratio, and NUE were observed in accessions , but to a greater extent under low-N stress. Six genomic prediction models were applied to the data, namely the Bayesian method Bayes C, Bayesian LASSO, Bayesian Ridge Regression, Ridge Regression-Best Linear Unbiased Prediction, Reproducing Kernel Hilbert Spaces, and randomForest. These models produced similar prediction accuracy of traits within the normal or low-N treatments, but the accuracy differed between the two treatments. ΔHT, LFW, LDW, and C were predicted slightly better under normal N with a mean Pearson r-value of 0.26, compared with r=0.22 under low N, while the prediction accuracies for Chl, N, C/N, and NUE were significantly improved under low-N stress with a mean r=0.45, compared with r=0.26 under normal N. The population panel contained three population structures, which generally had no effect on prediction accuracy. The moderate prediction accuracies obtained for N, C, and NUE under low-N stress are promising, and suggest a feasible means by which germplasm might be initially assessed for further detailed studies in breeding programs.

Advancing crested wheatgrass [Agropyron cristatum (L.) Gaertn.] breeding through genotyping-by-sequencing and genomic selection

Crested wheatgrass [Agropyron cristatum (L.) Gaertn.] provides high quality, highly palatable forage for early season grazing. Genetic improvement of crested wheatgrass has been challenged by its complex genome, outcrossing nature, long breeding cycle, and lack of informative molecular markers. Genomic selection (GS) has potential for improving traits of perennial forage species, and genotyping-by-sequencing (GBS) has enabled the development of genome-wide markers in non-model polyploid plants. An attempt was made to explore the utility of GBS and GS in crested wheatgrass breeding. Sequencing and phenotyping 325 genotypes representing 10 diverse breeding lines were performed. Bioinformatics analysis identified 827, 3,616, 14,090 and 46,136 single nucleotide polymorphism markers at 20%, 30%, 40% and 50% missing marker levels, respectively. Four GS models (BayesA, BayesB, BayesCπ, and rrBLUP) were examined for the accuracy of predicting nine agro-morphological and three nutritive value traits. Moderate accuracy (0.20 to 0.32) was obtained for the prediction of heading days, leaf width, plant height, clump diameter, tillers per plant and early spring vigor for genotypes evaluated at Saskatoon, Canada. Similar accuracy (0.29 to 0.35) was obtained for predicting fall regrowth and plant height for genotypes evaluated at Swift Current, Canada. The Bayesian models displayed similar or higher accuracy than rrBLUP. These findings show the feasibility of GS application for a non-model species to advance plant breeding.

Genomic Prediction and Indirect Selection for Grain Yield in US Pacific Northwest Winter Wheat Using Spectral Reflectance Indices from High-Throughput Phenotyping

Secondary traits from high-throughput phenotyping could be used to select for complex target traits to accelerate plant breeding and increase genetic gains. This study aimed to evaluate the potential of using spectral reflectance indices (SRI) for indirect selection of winter-wheat lines with high yield potential and to assess the effects of including secondary traits on the prediction accuracy for yield. A total of five SRIs were measured in a diversity panel, and F5 and doubled haploid wheat breeding populations planted between 2015 and 2018 in Lind and Pullman, WA. The winter-wheat panels were genotyped with 11,089 genotyping-by-sequencing derived markers. Spectral traits showed moderate to high phenotypic and genetic correlations, indicating their potential for indirect selection of lines with high yield potential. Inclusion of correlated spectral traits in genomic prediction models resulted in significant (p < 0.001) improvement in prediction accuracy for yield. Relatedness between training and test populations and heritability were among the principal factors affecting accuracy. Our results demonstrate the potential of using spectral indices as proxy measurements for selecting lines with increased yield potential and for improving prediction accuracy to increase genetic gains for complex traits in US Pacific Northwest winter wheat.