Genome analysis in Avena sativa reveals hidden breeding barriers and opportunities for oat improvement

A pangenome reveals LTR repeat dynamics as a major driver of genome evolution in Chenopodium

The genus Chenopodium L. is characterized by its wide geographic distribution and ecological adaptability. Species such as quinoa (Chenopodium quinoa Willd.) have served as domesticated staple crops for centuries. Wild Chenopodium species exhibit diverse niche adaptations and are important genetic reservoirs for beneficial agronomic traits, including disease resistance and climate hardiness. To harness the potential of the wild taxa for crop improvement, we developed a Chenopodium pangenome through the assembly and comparative analyses of 12 Chenopodium species that encompass the eight known genome types (A-H). Six of the species are new chromosome-scale assemblies, and many are polyploids; thus, a total of 20 genomes were included in the pangenome analyses. We show that the genomes vary dramatically in size with the D genome being the smallest (∼370 Mb) and the B genome being the largest (∼700 Mb) and that genome size was correlated with independent expansions of the Copia and Gypsy LTR retrotransposon families, suggesting that transposable elements have played a critical role in the evolution of the Chenopodium genomes. We annotated a total of 33,457 pan-Chenopodium gene families, of which ∼65% were classified as shell (2% private). Phylogenetic analysis clarified the evolutionary relationships among the genome lineages, notably resolving the taxonomic placement of the F genome while highlighting the uniqueness of the A genome in the Western Hemisphere. These genomic resources are particularly important for understanding the secondary and tertiary gene pools available for the improvement of the domesticated chenopods while furthering our understanding of the evolution and complexity within the genus.

Genomic strategies to facilitate breeding for increased β-Glucan content in oat (Avena sativa L.)

BACKGROUND

Hexaploid oat (Avena sativa L.) is a commercially important cereal crop due to its soluble dietary fiber β-glucan, a hemicellulose known to prevent cardio-vascular diseases. To maximize health benefits associated with the consumption of oat-based food products, breeding efforts have aimed at increasing the β-glucan content in oat groats. However, progress has been limited. To accelerate oat breeding efforts, we leveraged existing breeding datasets (1,230 breeding lines from South Dakota State University oat breeding program grown in multiple environments between 2015 and 2022) to conduct a genome-wide association study (GWAS) to increase our understanding of the genetic control of beta-glucan content in oats and to compare strategies to implement genomic selection (GS) to increase genetic gain for β-glucan content in oat.

RESULTS

Large variation for β-glucan content was observed with values ranging between 3.02 and 7.24%. An independent GWAS was performed for each breeding panel in each environment and identified 22 loci distributed over fourteen oat chromosomes significantly associated with β-glucan content. Comparison based on physical position showed that 12 out of 22 loci coincided with previously identified β-glucan QTLs, and three loci are in the vicinity of cellulose synthesis genes, Cellulose synthase-like (Csl). To perform a GWAS analysis across all breeding datasets, the β-glucan content of each breeding line was predicted for each of the 26 environments. The overall GWAS identified 73 loci, of which 15 coincided with loci identified for individual environments and 37 coincided with previously reported β-glucan QTLs not identified when performing the GWAS in single years. In addition, 21 novel loci were identified that were not reported in the previous studies. The proposed approach increased our ability to detect significantly associated markers. The comparison of multiple GS scenarios indicated that using a specific set of markers as a fixed effect in GS models did not increase the prediction accuracy. However, the use of multi-environment data in the training population resulted in an increase in prediction accuracy (0.61-0.72) as compared to single-year (0.28-0.48) data. The use of USDA-SoyWheOatBar-3 K genotyping array data resulted in a similar level of prediction accuracy as did genotyping-by-sequencing data.

CONCLUSION

This study identified and confirmed the location of multiple loci associated with β-glucan content. The proposed genomic strategies significantly increase both our ability to detect significant markers in GWAS and the accuracy of genomic predictions. The findings of this study can be useful to accelerate the genetic improvement of β-glucan content and other traits.

Phenotypic characterization and candidate gene analysis of a short kernel and brassinosteroid insensitive mutant from hexaploid oat (Avena sativa)

In an ethyl methanesulfonate oat (Avena sativa) mutant population we have found a mutant with striking differences to the wild-type (WT) cv. Belinda. We phenotyped the mutant and compared it to the WT. The mutant was crossed to the WT and mapping-by-sequencing was performed on a pool of F2 individuals sharing the mutant phenotype, and variants were called. The impacts of the variants on genes present in the reference genome annotation were estimated. The mutant allele frequency distribution was combined with expression data to identify which among the affected genes was likely to cause the observed phenotype. A brassinosteroid sensitivity assay was performed to validate one of the identified candidates. A literature search was performed to identify homologs of genes known to be involved in seed shape from other species. The mutant had short kernels, compact spikelets, altered plant architecture, and was found to be insensitive to brassinosteroids when compared to the WT. The segregation of WT and mutant phenotypes in the F2 population was indicative of a recessive mutation of a single locus. The causal mutation was found to be one of 123 single-nucleotide polymorphisms (SNPs) spanning the entire chromosome 3A, with further filtering narrowing this down to six candidate genes. In-depth analysis of these candidate genes and the brassinosteroid sensitivity assay suggest that a Pro303Leu substitution in AVESA.00010b.r2.3AG0419820.1 could be the causal mutation of the short kernel mutant phenotype. We identified 298 oat proteins belonging to orthogroups of previously published seed shape genes, with AVESA.00010b.r2.3AG0419820.1 being the only of these affected by a SNP in the mutant. The AVESA.00010b.r2.3AG0419820.1 candidate is functionally annotated as a GSK3/SHAGGY-like kinase with homologs in Arabidopsis, wheat, barley, rice, and maize, with several of these proteins having known mutants giving rise to brassinosteroid insensitivity and shorter seeds. The substitution in AVESA.00010b.r2.3AG0419820.1 affects a residue with a known gain-of function substitution in Arabidopsis BRASSINOSTEROID-INSENSITIVE2. We propose a gain-of-function mutation in AVESA.00010b.r2.3AG0419820.1 as the most likely cause of the observed phenotype, and name the gene AsGSK2.1. The findings presented here provide potential targets for oat breeders, and a step on the way towards understanding brassinosteroid signaling, seed shape and nutrition in oats.

Genome-Wide Association Study Uncovers Genomic Regions Associated with Coleoptile Length in a Worldwide Collection of Oat

The length of coleoptile is crucial for determining the sowing depth of oats in low-precipitation regions, which is significant for oat breeding programs. In this study, a diverse panel of 243 oat accessions was used to explore coleoptile length in two independent experiments. The panel exhibited significant variation in coleoptile length, ranging from 4.66 to 8.76 cm. Accessions from Africa, America, and the Mediterranean region displayed longer coleoptile lengths than those from Asia and Europe. Genome-wide association studies (GWASs) using 26,196 SNPs identified 34 SNPs, representing 32 quantitative trait loci (QTLs) significantly associated with coleoptile length. Among these QTLs, six were consistently detected in both experiments, explaining 6.43% to 10.07% of the phenotypic variation. The favorable alleles at these stable loci additively increased coleoptile length, offering insights for pyramid breeding. Gene Ontology (GO) analysis of the 350 candidate genes underlying the six stable QTLs revealed significant enrichment in cell development-related processes. Several phytochrome-related genes, including auxin transporter-like protein 1 and cytochrome P450 proteins, were found within these QTLs. Further validation of these loci will enhance our understanding of coleoptile length regulation. This study provides new insights into the genetic architecture of coleoptile length in oats.

Genome-wide identification and functional analysis of the cellulose synthase-like gene superfamily in common oat (Avena sativa L.)

Hemicelluloses constitute approximately one-third of the plant cell wall and can be used as a dietary fiber and food additive, and as raw materials for biofuels. Although genes involved in hemicelluloses synthesis have been investigated in some model plants, no comprehensive analysis has been conducted in common oat at present. In this study, we identified and systematically analyzed the cellulose synthase-like gene (Csl) family members in common oat and investigated them using various bioinformatics tools. The results showed that there are 76 members of the oat Csl gene family distributed on 17 chromosomes, and phylogenetic analysis indicated that the 76 Csl genes belong to the CslA, CslC, CslD, CslE, CslF, CslH, and CslJ subfamilies. A total of 14 classes of cis-acting elements were identified in the promoter regions, including hormone response, light response, cell development, and defense stress elements. The collinearity analysis identified 28 pairs of segmentally duplicated genes, most of which were found on chromosomes 2D and 6A. Expression pattern analysis showed that oat Csl genes display strong tissue-specific expression; of the 76 Csl genes, 33 were significantly up-regulated in stems and 30 were up-regulated in immature seeds. The expression of most members of the AsCsl gene family is repressed by abiotic stress, while the expression of some members is up-regulated by light. Immunoelectron microscopy shows that the product of AsCsl61, a member of CslF subfamily, mediates (1,3; 1,4)-β-D-glucan synthesis in transgenic Arabidopsis. These findings provide a fundamental understanding of the structural, functional, and evolutionary features of the oat Csl genes and may contribute to our general understanding of hemicellulose biosynthesis. Moreover, this information will be helpful in designing experiments for genetic manipulation of mixed-linkage glucan (MLG) synthesis with the goal of quality improvement in oat.

Characterization of the Metallothionein Gene Family in Avena sativa L. and the Gene Expression during Seed Germination and Heavy Metal Stress

Metallothioneins (MTs) are a family of small proteins rich in cysteine residues. The sulfhydryl group of metallothioneins can bind to metal ions, maintaining metal homeostasis and protecting the cells from damage caused by toxic heavy metals. Moreover, MTs can function as reactive oxygen species scavengers since cysteine thiols undergo reversible and irreversible oxidation. Here, we identified 21 metallothionein genes (AsMTs) in the oat (Avena sativa L.) genome, which were divided into four types depending on the amino acid sequences of putative proteins encoded by identified genes. Analysis of promoter sequences showed that MTs might respond to a variety of stimuli, including biotic and abiotic stresses and phytohormones. The results of qRT-PCR showed that all four types of AsMTs are differentially expressed during the first 48 hours of seed germination. Moreover, stress induced by the application of zinc, cadmium, and a mixture of zinc and cadmium affects the expression of oat MTs variously depending on the MT type, indicating that AsMT1-4 fulfil different roles in plant cells.

Historical Selection, Adaptation Signatures, and Ambiguity of Introgressions in Wheat

Wheat was one of the crops domesticated in the Fertile Crescent region approximately 10,000 years ago. Despite undergoing recent polyploidization, hull-to-free-thresh transition events, and domestication bottlenecks, wheat is now grown in over 130 countries and accounts for a quarter of the world's cereal production. The main reason for its widespread success is its broad genetic diversity that allows it to thrive in different environments. To trace historical selection and hybridization signatures, genome scans were performed on two datasets: approximately 113K SNPs from 921 predominantly bread wheat accessions and approximately 110K SNPs from about 400 wheat accessions representing all ploidy levels. To identify environmental factors associated with the loci, a genome-environment association (GEA) was also performed. The genome scans on both datasets identified a highly differentiated region on chromosome 4A where accessions in the first dataset were dichotomized into a group (n = 691), comprising nearly all cultivars, wild emmer, and most landraces, and a second group (n = 230), dominated by landraces and spelt accessions. The grouping of cultivars is likely linked to their potential ancestor, bread wheat cv. Norin-10. The 4A region harbored important genes involved in adaptations to environmental conditions. The GEA detected loci associated with latitude and temperature. The genetic signatures detected in this study provide insight into the historical selection and hybridization events in the wheat genome that shaped its current genetic structure and facilitated its success in a wide spectrum of environmental conditions. The genome scans and GEA approaches applied in this study can help in screening the germplasm housed in gene banks for breeding, and for conservation purposes.

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

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

Genome-Wide Association Mapping of QTL Underlying Groat Protein Content of a Diverse Panel of Oat Accessions

Groat protein content (GPC) is a key quality trait attribute in oat. Understanding the variation of GPC in oat germplasms and identifying genomic regions associated with GPC are essential for improving this trait. In this study, the GPC of 174 diverse oat accessions was evaluated in three field trials. The results showed a wide variation in GPC, ranging from 6.97% to 22.24% in this panel. Hulless oats displayed a significantly higher GPC compared to hulled oats across all environments. A GWAS analysis was performed based on 38,313 high-quality SNPs, which detected 27 non-redundant QTLs with 41 SNPs significantly associated with GPC. Two QTLs on chromosome 6C (QTL16) and 4D (QTL11) were consistently detected in multiple environments, with QTL16 being the most significant and explaining the highest proportion of the phenotypical variation in all tested environments except in CZ20. Haplotype analysis showed that the favorable haplotypes for GPC are more prevalent in hulless oats. These findings provide a foundation for future efforts to incorporate favorable alleles into new cultivars through introgression, fine mapping, and cloning of promising QTLs.

Non-B-form DNA tends to form in centromeric regions and has undergone changes in polyploid oat subgenomes

Centromeres are the specialized regions of the chromosomes that direct faithful chromosome segregation during cell division. Despite their functional conservation, centromeres display features of rapidly evolving DNA and wide evolutionary diversity in size and organization. Previous work found that the noncanonical B-form DNA structures are abundant in the centromeres of several eukaryotic species with a possible implication for centromere specification. Thus far, systematic studies into the organization and function of non-B-form DNA in plants remain scarce. Here, we applied the oat system to investigate the role of non-B-form DNA in centromeres. We conducted chromatin immunoprecipitation sequencing using an antibody to the centromere-specific histone H3 variant (CENH3); this accurately positioned oat centromeres with different ploidy levels and identified a series of centromere-specific sequences including minisatellites and retrotransposons. To define genetic characteristics of oat centromeres, we surveyed the repeat sequences and found that dyad symmetries were abundant in oat centromeres and were predicted to form non-B-DNA structures in vivo. These structures including bent DNA, slipped DNA, Z-DNA, G-quadruplexes, and R-loops were prone to form within CENH3-binding regions. Dynamic conformational changes of predicted non-B-DNA occurred during the evolution from diploid to tetraploid to hexaploid oat. Furthermore, we applied the single-molecule technique of AFM and DNA:RNA immunoprecipitation with deep sequencing to validate R-loop enrichment in oat centromeres. Centromeric retrotransposons exhibited strong associations with R-loop formation. Taken together, our study elucidates the fundamental character of non-B-form DNA in the oat genome and reveals its potential role in centromeres.

Advancing understanding of oat phenology for crop adaptation

Oat (Avena sativa) is an annual cereal grown for forage, fodder and grain. Seasonal flowering behaviour, or phenology, is a key contributor to the success of oat as a crop. As a species, oat is a vernalization-responsive long-day plant that flowers after winter as days lengthen in spring. Variation in both vernalization and daylength requirements broadens adaptation of oat and has been used to breed modern cultivars with seasonal flowering behaviours suited to different regions, sowing dates and farming practices. This review examines the importance of variation in oat phenology for crop adaptation. Strategies to advance understanding of the genetic basis of oat phenology are then outlined. These include the potential to transfer knowledge from related temperate cereals, particularly wheat (Triticum aestivum) and barley (Hordeum vulgare), to provide insights into the potential molecular basis of variation in oat phenology. Approaches that use emerging genomic resources to directly investigate the molecular basis of oat phenology are also described, including application of high-resolution genome-wide diversity surveys to map genes linked to variation in flowering behaviour. The need to resolve the contribution of individual phenology genes to crop performance by developing oat genetic resources, such as near-isogenic lines, is emphasised. Finally, ways that deeper knowledge of oat phenology can be applied to breed improved varieties and to inform on-farm decision-making are outlined.

Archetypes of inflorescence: genome-wide association networks of panicle morphometric, growth, and disease variables in a multiparent oat population

There is limited information regarding the morphometric relationships of panicle traits in oat (Avena sativa) and their contribution to phenology and growth, physiology, and pathology traits important for yield. To model panicle growth and development and identify genomic regions associated with corresponding traits, 10 diverse spring oat mapping populations (n = 2,993) were evaluated in the field and 9 genotyped via genotyping-by-sequencing. Representative panicles from all progeny individuals, parents, and check lines were scanned, and images were analyzed using manual and automated techniques, resulting in over 60 unique panicle, rachis, and spikelet variables. Spatial modeling and days to heading were used to account for environmental and phenological variances, respectively. Panicle variables were intercorrelated, providing reproducible archetypal and growth models. Notably, adult plant resistance for oat crown rust was most prominent for taller, stiff stalked plants having a more open panicle structure. Within and among family variance for panicle traits reflected the moderate-to-high heritability and mutual genome-wide associations (hotspots) with numerous high-effect loci. Candidate genes and potential breeding applications are discussed. This work adds to the growing genetic resources for oat and provides a unique perspective on the genetic basis of panicle architecture in cereal crops.

The mosaic oat genome gives insights into a uniquely healthy cereal crop

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia^1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

Assembly of the hexaploid oat genome and its diploid and tetraploid relatives clarifies the evolutionary history of oat and allows mapping of genes for agronomic traits.