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

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Patterns of genomic variation reveal a single evolutionary origin of the wild allotetraploid Mimulus sookensis

Polyploidy occurs across the tree of life and is especially common in plants. Because newly formed cytotypes are often incompatible with their progenitors, polyploidy is also said to trigger "instantaneous" speciation. If a polyploid can self-fertilize or reproduce asexually, it is even possible for one individual to produce an entirely new lineage, but how often this scenario occurs is unclear. Here, we investigate the evolutionary history of the wild allotetraploid Mimulus sookensis, which was formed through hybridization between self-compatible, diploid species in the Mimulus guttatus complex. We generate a chromosome-scale reference assembly for M. sookensis and define its distinct subgenomes. Despite previous reports suggesting multiple origins of this highly selfing polyploid, we discover patterns of population genomic variation that provide unambiguous support for a single origin. One M. sookensis subgenome is clearly derived from the selfer Mimulus nasutus, which organellar variation suggests is the maternal progenitor. The ancestor of the other subgenome is less certain, but it shares variation with both Mimulus decorus and M. guttatus, two outcrossing diploids with geographic ranges that overlap broadly with M. sookensis. This study establishes M. sookensis as an example of instantaneous speciation, likely facilitated by the polyploid's predisposition to self-fertilize.

Genome-wide analysis of MYB transcription factor family and AsMYB1R subfamily contribution to ROS homeostasis regulation in Avena sativa under PEG-induced drought stress

BACKGROUND

The myeloblastosis (MYB) transcription factor (TF) family is one of the largest and most important TF families in plants, playing an important role in a life cycle and abiotic stress.

RESULTS

In this study, 268 Avena sativa MYB (AsMYB) TFs from Avena sativa were identified and named according to their order of location on the chromosomes, respectively. Phylogenetic analysis of the AsMYB and Arabidopsis MYB proteins were performed to determine their homology, the AsMYB1R proteins were classified into 5 subgroups, and the AsMYB2R proteins were classified into 34 subgroups. The conserved domains and gene structure were highly conserved among the subgroups. Eight differentially expressed AsMYB genes were screened in the transcriptome of transcriptional data and validated through RT-qPCR. Three genes in AsMYB2R subgroup, which are related to the shortened growth period, stomatal closure, and nutrient and water transport by PEG-induced drought stress, were investigated in more details. The AsMYB1R subgroup genes LHY and REV 1, together with GST, regulate ROS homeostasis to ensure ROS signal transduction and scavenge excess ROS to avoid oxidative damage.

CONCLUSION

The results of this study confirmed that the AsMYB TFs family is involved in the homeostatic regulation of ROS under drought stress. This lays the foundation for further investigating the involvement of the AsMYB TFs family in regulating A. sativa drought response mechanisms.

Hijacking a rapid and scalable metagenomic method reveals subgenome dynamics and evolution in polyploid plants

Premise

The genomes of polyploid plants archive the evolutionary events leading to their present forms. However, plant polyploid genomes present numerous hurdles to the genome comparison algorithms for classification of polyploid types and exploring genome dynamics.

Methods

Here, the problem of intra- and inter-genome comparison for examining polyploid genomes is reframed as a metagenomic problem, enabling the use of the rapid and scalable MinHashing approach. To determine how types of polyploidy are described by this metagenomic approach, plant genomes were examined from across the polyploid spectrum for both k-mer composition and frequency with a range of k-mer sizes. In this approach, no subgenome-specific k-mers are identified; rather, whole-chromosome k-mer subspaces were utilized.

Results

Given chromosome-scale genome assemblies with sufficient subgenome-specific repetitive element content, literature-verified subgenomic and genomic evolutionary relationships were revealed, including distinguishing auto- from allopolyploidy and putative progenitor genome assignment. The sequences responsible were the rapidly evolving landscape of transposable elements. An investigation into the MinHashing parameters revealed that the downsampled k-mer space (genomic signatures) produced excellent approximations of sequence similarity. Furthermore, the clustering approach used for comparison of the genomic signatures is scrutinized to ensure applicability of the metagenomics-based method.

Discussion

The easily implementable and highly computationally efficient MinHashing-based sequence comparison strategy enables comparative subgenomics and genomics for large and complex polyploid plant genomes. Such comparisons provide evidence for polyploidy-type subgenomic assignments. In cases where subgenome-specific repeat signal may not be adequate given a chromosomes' global k-mer profile, alternative methods that are more specific but more computationally complex outperform this approach.

Formation of lipid-derived volatile products through lipoxygenase (LOX)- and hydroperoxide lyase (HPL)- mediated pathway in oat, barley and soy bean

The aim of this study was to explore the formation of volatile lipid oxidation products by the lipoxygenase (LOX)-hydroperoxide lyase (HPL)-mediated pathway in oat, barley and soy bean. LOX activity was found only in barley and soy bean samples, but the lipase and HPL activity was detected in all samples. HPL showed particularly high activity with 13-hydroperoxides, while the activity was quite low when using 9-hydroperoxides, especially in the oat and barley. The optimum pH for HPL in different samples was similar, i.e., pH 6-7. In this condition, the volatile compounds formed dramatically with aldehydes and furans as the dominant products. Furthermore, a remarkable enzymatic degradation of lipids occurred during the preparation of food models with highly refined rapeseed oil (RO) and rapeseed oil fatty acid (ROFA) emulsions, where the ROFAs were more prone to oxidation than RO. This study shows the significance of lipid-degrading enzymes in plant-food flavour formation.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Assessing myBaits Target Capture Sequencing Methodology Using Short-Read Sequencing for Variant Detection in Oat Genomics and Breeding

The integration of target capture systems with next-generation sequencing has emerged as an efficient tool for exploring specific genetic regions with a high resolution and facilitating the rapid discovery of novel alleles. Despite these advancements, the application of targeted sequencing methodologies, such as the myBaits technology, in polyploid oat species remains relatively unexplored. In this study, we utilized the myBaits target capture method offered by Daicel Arbor Biosciences to detect variants and assess their reliability for variant detection in oat genomics and breeding. Ten oat genotypes were carefully chosen for targeted sequencing, focusing on specific regions on chromosome 2A to detect variants. The selected region harbors 98 genes. Precisely designed baits targeting the genes within these regions were employed for the target capture sequencing. We employed various mappers and variant callers to identify variants. After the identification of variants, we focused on the variants identified via all variants callers to assess the applicability of the myBaits sequencing methodology in oat breeding. In our efforts to validate the identified variants, we focused on two SNPs, one deletion and one insertion identified via all variant callers in the genotypes KF-318 and NOS 819111-70 but absent in the remaining eight genotypes. The Sanger sequencing of targeted SNPs failed to reproduce target capture data obtained through the myBaits technology. Similarly, the validation of deletion and insertion variants via high-resolution melting (HRM) curve analysis also failed to reproduce target capture data, again suggesting limitations in the reliability of the myBaits target capture sequencing using short-read sequencing for variant detection in the oat genome. This study shed light on the importance of exercising caution when employing the myBaits target capture strategy for variant detection in oats. This study provides valuable insights for breeders seeking to advance oat breeding efforts and marker development using myBaits target capture sequencing, emphasizing the significance of methodological sequencing considerations in oat genomics research.

Genome Variability in Artificial Allopolyploid Hybrids of Avena sativa L. and Avena macrostachya Balansa ex Coss. et Durieu Based on Marker Sequences of Satellite DNA and the ITS1-5.8S rDNA Region

Artificial hybrids between cultivated Avena species and wild Avena macrostachya that possess genes for resistance to biotic and abiotic stresses can be important for oat breeding. For the first time, a comprehensive study of genomes of artificial fertile hybrids Avena sativa × Avena macrostachya and their parental species was carried out based on the chromosome FISH mapping of satellite DNA sequences (satDNAs) and also analysis of intragenomic polymorphism in the 18S-ITS1-5.8S rDNA region, using NGS data. Chromosome distribution patterns of marker satDNAs allowed us to identify all chromosomes in the studied karyotypes, determine their subgenomic affiliation, and detect several chromosome rearrangements. Based on the obtained cytogenomic data, we revealed differences between two A. macrostachya subgenomes and demonstrated that only one of them was inherited in the studied octoploid hybrids. Ribotype analyses showed that the second major ribotype of A. macrostachya was species-specific and was not represented in rDNA pools of the octoploids, which could be related to the allopolyploid origin of this species. Our results indicate that the use of marker satDNAs in cytogenomic studies can provide important data on genomic relationships within Avena allopolyploid species and hybrids, and also expand the potential for interspecific crosses for breeding.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Differential expression and global analysis of miR156/SQUAMOSA promoter binding-like proteins (SPL) module in oat

SQUAMOSA promoter binding-like proteins (SPLs) are important transcription factors that influence growth phase transition and reproduction in plants. SPLs are targeted by miR156 but the SPL/miR156 module is completely unknown in oat. We identified 28 oat SPL genes (AsSPLs) distributed across all 21 oat chromosomes except for 4C and 6D. The oat- SPL gene family represented six of eight SPL phylogenetic groups, with no AsSPLs in groups 3 and 7. A novel oat miR156 (AsmiR156) family with 21 precursors divided into 7 groups was characterized. A total of 16 AsSPLs were found to be targeted by AsmiR156. Intriguingly, AsSPL3s showed high transcript abundance during early inflorescence (GS-54), as compared to the lower abundance of AsmiR156, indicating their role in reproductive development. Unravelling the SPL/miR156 regulatory hub and alterations in expression patterns of AsSPLs could provide an essential toolbox for genetic improvement in the cultivated oat.

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.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Determining the Impact of Genotype × Environment on Oat Protein Isolate Composition Using HPLC and LC-MS Techniques

The effect of genotype and environment on oat protein composition was analyzed through size exclusion-high-performance liquid chromatography (SE-HPLC) and liquid chromatography-mass spectrometry (LC-MS) to characterize oat protein isolate (OPI) extracted from three genotypes grown at three locations in the Canadian Prairies. SE-HPLC identified four fractions in OPI, including polymeric globulins, avenins, glutelins, and albumins, and smaller proteins. The protein composition was dependent on the environment, rather than the genotype. The proteins identified through LC-MS were grouped into eight categories, including globulins, prolamins/avenins, glutelins, enzymes/albumins, enzyme inhibitors, heat shock proteins, grain softness proteins, and allergenic proteins. Three main globulin protein types were also identified, including the P14812|SSG2-12S seed storage globulin, the Q6UJY8_TRITU-globulin, and the M7ZQM3_TRIUA-Globulin-1 S. Principal component analysis indicated that samples from Manitoba showed a positive association with the M7ZQM3_TRIUA-Globulin-1 S allele and Q6UJY8_TRITU-globulin, while samples from Alberta and Saskatchewan had a negative association with them. The results show that the influence of G × E on oat protein fractions and their relative composition is crucial to understanding genotypes' behavior in response to different environments.

Genome-Wide Identification of NAC Family Genes in Oat and Functional Characterization of AsNAC109 in Abiotic Stress Tolerance

The plant-specific NAC gene family is one of the largest transcription factor families, participating in plant growth regulation and stress response. Despite extensive characterization in various plants, our knowledge of the NAC family in oat is lacking. Herein, we identified 333 NAC genes from the latest release of the common oat genome. We provide a comprehensive overview of the oat NAC gene family, covering gene structure, chromosomal localization, phylogenetic characteristics, conserved motif compositions, and gene duplications. AsNAC gene expression in different tissues and the response to various abiotic stresses were characterized using RT-qPCR. The main driver of oat NAC gene family expansion was identified as segmental duplication using collinearity analysis. In addition, the functions of AsNAC109 in regulating abiotic stress tolerance in Arabidopsis were clarified. This is the first genome-wide investigation of the NAC gene family in cultivated oat, which provided a unique resource for subsequent research to elucidate the mechanisms responsible for oat stress tolerance and provides valuable clues for the improvement of stress resistance in cultivated oat.

Efficacy of oats in dyslipidemia: a systematic review and meta-analysis

Background: In recent years, oats' effect on lowering serum cholesterol has been recognized. However, no systematic reviews summarized the effect of daily consumption of oat-based products on serum lipids in patients with dyslipidemia. Methods: We searched eight databases and two clinical trial registries from inception to July 31, 2023. We included randomized controlled trials (RCTs) evaluating the efficacy of oat-based products (≥4 weeks) on lipid levels or cardiovascular events in patients with dyslipidemia. Two authors independently screened articles, extracted data and assessed the risk of bias of included studies with Cochrane risk-of-bias tool 2.0. We used STATA 17.0 to conduct meta-analysis and Grading of Recommendations Assessment, Development and Evaluation (GRADE) to assess the certainty of evidence. Results: We finally included 17 eligible trials with 1731 subjects. The oat intervention varied from oat β-glucan-based products to oat bran-based products and wholegrain oat. Overall, the risk of bias of included trials was high or some concerns were noted because of the inadequate randomization, allocation concealment, and inappropriate data analysis method. Compared to the placebo or usual diet, one study indicated that oat-based products have no significant difference in major cardiovascular events. Pooled estimates showed that oat-based products may result in a large reduction in LDL-C (WMD, -0.24 mmol L-1; 95% CI: -0.33, -0.15) (moderate certainty) and TC (WMD, -0.32 mmol L-1; 95% CI: -0.48, -0.17) (moderate certainty). Compared to other diets (mainly other cereals), oat-based products probably reduce the level of LDL-C (WMD, -0.17 mmol L-1; 95% CI: -0.25, -0.08) (moderate certainty) and TC (WMD, -0.21 mmol L-1; 95% CI: -0.30, -0.12) (moderate certainty). Both groups showed that oat-based products had little effect on HDL-C and TG (moderate certainty). Oat-related adverse events were mostly gastrointestinal such as diarrhea, nausea, and flatulence being the most prevalent. Conclusions: Oat-based products may reduce TC and LDL-C, but have little effect on TG, HDL-C, and major cardiovascular events in patients with dyslipidemia.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

The complex polyploid genome architecture of sugarcane

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

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.

An improved chromosome-level genome assembly of perennial ryegrass (Lolium perenne L.)

This work is an update and extension of the previously published article "Ultralong Oxford Nanopore Reads Enable the Development of a Reference-Grade Perennial Ryegrass Genome Assembly" by Frei et al. The published genome assembly of the doubled haploid perennial ryegrass (Lolium perenne L.) genotype Kyuss (Kyuss v1.0) marked a milestone for forage grass research and breeding. However, order and orientation errors may exist in the pseudo-chromosomes of Kyuss, since barley (Hordeum vulgare L.), which diverged 30 million years ago from perennial ryegrass, was used as the reference to scaffold Kyuss. To correct for structural errors possibly present in the published Kyuss assembly, we de novo assembled the genome again and generated 50-fold coverage high-throughput chromosome conformation capture (Hi-C) data to assist pseudo-chromosome construction. The resulting new chromosome-level assembly Kyuss v2.0 showed improved quality with high contiguity (contig N50 = 120 Mb), high completeness (total BUSCO score = 99%), high base-level accuracy (QV = 50), and correct pseudo-chromosome structure (validated by Hi-C contact map). This new assembly will serve as a better reference genome for Lolium spp. and greatly benefit the forage and turf grass research community.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-Wide Analysis of the Oat (Avena sativa) HSP90 Gene Family Reveals Its Identification, Evolution, and Response to Abiotic Stress

Oats (Avena sativa) are an important cereal crop and cool-season forage worldwide. Heat shock protein 90 (HSP90) is a protein ubiquitously expressed in response to heat stress in almost all plants. To date, the HSP90 gene family has not been comprehensively reported in oats. Herein, we have identified twenty HSP90 genes in oats and elucidated their evolutionary pathways and responses to five abiotic stresses. The gene structure and motif analyses demonstrated consistency across the phylogenetic tree branches, and the groups exhibited relative structural conservation. Additionally, we identified ten pairs of segmentally duplicated genes in oats. Interspecies synteny analysis and orthologous gene identification indicated that oats share a significant number of orthologous genes with their ancestral species; this implies that the expansion of the oat HSP90 gene family may have occurred through oat polyploidization and large fragment duplication. The analysis of cis-acting elements revealed their influential role in the expression pattern of HSP90 genes under abiotic stresses. Analysis of oat gene expression under high-temperature, salt, cadmium (Cd), polyethylene glycol (PEG), and abscisic acid (ABA) stresses demonstrated that most AsHSP90 genes were significantly up-regulated by heat stress, particularly AsHSP90-7, AsHSP90-8, and AsHSP90-9. This study offers new insights into the amplification and evolutionary processes of the AsHSP90 protein, as well as its potential role in response to abiotic stresses. Furthermore, it lays the groundwork for understanding oat adaptation to abiotic stress, contributing to research and applications in plant breeding.

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.

Genome-wide analysis and expression of the aquaporin gene family in Avena sativa L

Background

Oat (Avena sativa L.) belongs to the early maturity grass subfamily of the Gramineae subfamily oats (Avena) and has excellent characteristics, such as tolerance to barrenness, salt, cold, and drought. Aquaporin (AQP) proteins belong to the major intrinsic protein (MIP) superfamily, are widely involved in plant growth and development, and play an important role in abiotic stress responses. To date, previous studies have not identified or analyzed the AsAQP gene family system, and functional studies of oat AQP genes in response to drought, cold, and salt stress have not been performed.

Methods

In this study, AQP genes (AsAQP) were identified from the oat genome, and various bioinformatics data on the AQP gene family, gene structure, gene replication, promoters and regulatory networks were analyzed. Quantitative real-time PCR technology was used to verify the expression patterns of the AQP gene family in different oat tissues under different abiotic stresses.

Results

In this study, a total of 45 AQP genes (AsAQP) were identified from the oat reference genome. According to a phylogenetic analysis, 45 AsAQP were divided into 4 subfamilies (PIP, SIP, NIP, and TIP). Among the 45 AsAQP, 23 proteins had interactions, and among these, 5AG0000633.1 had the largest number of interacting proteins. The 20 AsAQP genes were expressed in all tissues, and their expression varied greatly among different tissues and organs. All 20 AsAQP genes responded to salt, drought and cold stress. The NIP subfamily 6Ag0000836.1 gene was significantly upregulated under different abiotic stresses and could be further verified as a key candidate gene.

Conclusion

The findings of this study provide a comprehensive list of members and their sequence characteristics of the AsAQP protein family, laying a solid theoretical foundation for further functional analysis of AsAQP in oats. This research also offers valuable reference for the creation of stress-tolerant oat varieties through genetic engineering techniques.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

Volatile metabolomics reveals the characteristics of the unique flavor substances in oats

Oats is a cereal well known for its high nutritional value and unique flavor. This study investigated the metabolomics data from oats, wheat, and barley using broadly targeted GC-MS metabonomic techniques. A total of 437 volatile organic compounds (VOCs) were identified, of which 414 were shared metabolites, with three metabolites unique to oats. Three hundred and seven differentially accumulated metabolites (DAMs) were screened from all the comparison groups, of which 27 metabolites were shared by oats and barley, and 121 shared by oats and wheat. Terpenoids and esters were the key metabolites determining the differences in flavor. A KEGG analysis indicated that the alpha-linolenic acid and phenylalanine pathways were the most significant metabolic pathways. The 42 DAMs found may be the main substances leading to the flavor differences between the different varieties. Overall, this study reveals the main reasons for the unique flavor of oats through metabolomic evidence.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

From Petri Dish to Field: Plant Tissue Culture and Genetic Engineering of Oats for Improved Agricultural Outcomes

Oats (Avena sativa) hold immense economic and nutritional value as a versatile crop. They have long been recognized as an exceptional choice for human consumption and animal feed. Oats' unique components, including proteins, starches, and β-glucans, have led to its widespread use in various food products such as bread, noodles, flakes, and milk. The popularity of oat milk as a vegan alternative to dairy milk has soared due to the increasing number of vegetarians/vegans and growing environmental awareness. Oat milk offers a sustainable option with reduced greenhouse gas emissions during its production, rendering it an appropriate choice for individuals who are lactose-intolerant or have dairy allergies. To ensure improved adaptability and enhanced nutrition, the development of new oat varieties is crucial, considering factors like cultivation, climate, and growing conditions. Plant cell culture plays a crucial role in both traditional and contemporary breeding methods. In classical breeding, plant cell culture facilitates the rapid production of double haploid plants, which can be employed to accelerate the breeding process. In modern breeding methods, it enables genetic manipulation and precise genome editing at the cellular level. This review delves into the importance of oats and their diverse applications, highlighting the advantages of plant cell culture in both classical and modern breeding methods. Specifically, it provides an overview of plant tissue culture, encompassing genetic transformation, haploid technology, protoplast technology, and genome editing.

Transcriptomic response for revealing the molecular mechanism of oat flowering under different photoperiods

Proper flowering is essential for the reproduction of all kinds of plants. Oat is an important cereal and forage crop; however, its cultivation is limited because it is a long-day plant. The molecular mechanism by which oats respond to different photoperiods is still unclear. In this study, oat plants were treated under long-day and short-day photoperiods for 10 days, 15 days, 20 days, 25 days, 30 days, 40 days and 50 days, respectively. Under the long-day treatment, oats entered the reproductive stage, while oats remained vegetative under the short-day treatment. Forty-two samples were subjected to RNA-Seq to compare the gene expression patterns of oat under long- and short-day photoperiods. A total of 634-5,974 differentially expressed genes (DEGs) were identified for each time point, while the floral organ primordium differentiation stage showed the largest number of DEGs, and the spikelet differentiation stage showed the smallest number. Gene Ontology (GO) analysis showed that the plant hormone signaling transduction and hormone metabolism processes significantly changed in the photoperiod regulation of flowering time in oat. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Mapman analysis revealed that the DEGs were mainly concentrated in the circadian rhythm, protein antenna pathways and sucrose metabolism process. Additionally, transcription factors (TFs) involved in various flowering pathways were explored. Combining all this information, we established a molecular model of oat flowering induced by a long-day photoperiod. Taken together, the long-day photoperiod has a large effect at both the morphological and transcriptomic levels, and these responses ultimately promote flowering in oat. Our findings expand the understanding of oat as a long-day plant, and the explored genes could be used in molecular breeding to help break its cultivation limitations in the future.

Transcriptional and protein structural characterization of homogentisate phytyltransferase genes in barley, wheat, and oat

BACKGROUND

Homogentisate phytyltransferase (HPT) is the critical enzyme for the biosynthesis of tocopherols (vitamin E), which are the major lipid-soluble antioxidants and help plants adapt to various stress conditions. HPT is generally strictly conserved in various plant genomes; however, a divergent lineage HPT2 was identified recently in some Triticeae species. The molecular function and transcriptional profiles of HPT2 remain to be characterized.

RESULTS

In this study, we performed comprehensive transcriptome data mining of HPT1 and HPT2 in different tissues and stages of barley (Hordeum vulgare), wheat (Triticum aestivum), and oat (Avena sativa), followed by qRT-PCR experiments on HPT1 and HPT2 in different tissues of barley and wheat. We found that the common HPT1 genes (HvHPT1, TaHPT1s, and AsHPT1s) displayed a conserved transcriptional pattern in the three target species and were universally transcribed in various tissues, with a notable preference in leaf. In contrast, HPT2 genes (HvHPT2, TaHPT2, and AsHPT2) were specifically transcribed in spike (developmentally up-regulated) and shoot apex tissues, displaying a divergent tissue-specific pattern. Cis-regulatory elements prediction in the promoter region identified common factors related to light-, plant hormone-, low temperature-, drought- and defense- responses in both HPT1s and HPT2s. We observed the transcriptional up-regulation of HvHPT1 and HvHPT2 under various stress conditions, supporting their conserved function in environmental adaption. We detected a clear, relaxed selection pressure in the HPT2 lineage, consistent with the predicted evolution pattern following gene duplication. Protein structural modelling and substrate docking analyses identified putative catalytic amino acid residues for HvHPT1 and HvHPT2, which are strictly conserved and consistent with their function in vitamin E biosynthesis.

CONCLUSIONS

We confirmed the presence of two lineages of HPT in Triticeae and Aveninae, including hexaploid oat, and characterized their transcriptional profiles based on transcriptome and qRT-PCR data. HPT1s were ubiquitously transcribed in various tissues, whilst HPT2s were highly expressed in specific stages and tissue. The active transcription of HPT2s, together with its conserved cis-elements and protein structural features, support HPT2s' role in tocopherol production in Triticeae. This study is the first protein structural analysis on the membrane-bound plant HPTs and provides valuable insights into its catalytic mechanism.

Genome-Wide Investigation and Expression Analysis of the Catalase Gene Family in Oat Plants (Avena sativa L.)

Through the degradation of reactive oxygen species (ROS), different antioxidant enzymes, such as catalase (CAT), defend organisms against oxidative stress. These enzymes are crucial to numerous biological functions, like plant development and defense against several biotic and abiotic stresses. However, despite the major economic importance of Avena sativa around the globe, little is known about the CAT gene's structure and organization in this crop. Thus, a genome-wide investigation of the CAT gene family in oat plants has been carried out to characterize the potential roles of those genes under different stressors. Bioinformatic approaches were used in this study to predict the AvCAT gene's structure, secondary and tertiary protein structures, physicochemical properties, phylogenetic tree, and expression profiling under diverse developmental and biological conditions. A local Saudi oat variety (AlShinen) was used in this work. Here, ten AvCAT genes that belong to three groups (Groups I-III) were identified. All identified CATs harbor the two conserved domains (pfam00199 and pfam06628), a heme-binding domain, and a catalase activity motif. Moreover, identified AvCAT proteins were located in different compartments in the cell, such as the peroxisome, mitochondrion, and cytoplasm. By analyzing their promoters, different cis-elements were identified as being related to plant development, maturation, and response to different environmental stresses. Gene expression analysis revealed that three different AvCAT genes belonging to three different subgroups showed noticeable modifications in response to various stresses, such as mannitol, salt, and ABA. As far as we know, this is the first report describing the genome-wide analysis of the oat catalase gene family, and these data will help further study the roles of catalase genes during stress responses, leading to crop improvement.

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.

The dynamics of touch-responsive gene expression in cereals

Wind, rain, herbivores, obstacles, neighbouring plants, etc. provide important mechanical cues to steer plant growth and survival. Mechanostimulation to stimulate yield and stress resistance of crops is of significant research interest, yet a molecular understanding of transcriptional responses to touch is largely absent in cereals. To address this, we performed whole-genome transcriptomics following mechanostimulation of wheat, barley, and the recent genome-sequenced oat. The largest transcriptome changes occurred ±25 min after touching, with most of the genes being upregulated. While most genes returned to basal expression level by 1-2 h in oat, many genes retained high expression even 4 h post-treatment in barley and wheat. Functional categories such as transcription factors, kinases, phytohormones, and Ca2+ regulation were affected. In addition, cell wall-related genes involved in (hemi)cellulose, lignin, suberin, and callose biosynthesis were touch-responsive, providing molecular insight into mechanically induced changes in cell wall composition. Furthermore, several cereal-specific transcriptomic footprints were identified that were not observed in Arabidopsis. In oat and barley, we found evidence for systemic spreading of touch-induced signalling. Finally, we provide evidence that both the jasmonic acid-dependent and the jasmonic acid-independent pathways underlie touch-signalling in cereals, providing a detailed framework and marker genes for further study of (a)biotic stress responses in cereals.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum

BACKGROUND

Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation.

RESULTS

Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum.

CONCLUSIONS

Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

Triticeae crop genome biology: an endless frontier

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Effect of Dynamic High-Pressure Microfluidizer on Physicochemical and Microstructural Properties of Whole-Grain Oat Pulp

By avoiding the filtration step and utilizing the whole components of oats, the highest utilization rate of raw materials, improving the nutritional value of products and reducing environmental pollution, can be achieved in the production of whole-grain oat drinks. This study innovatively introduced a dynamic high-pressure microfluidizer (DHPM) into the processing of whole-grain oat pulp, which aimed to achieve the efficient crushing, homogenizing and emulsification of starch, dietary fiber and other substances. Due to DHPM processing, the instability index and slope value were reduced, whereas the β-glucan content, soluble protein content and soluble dietary fiber content were increased. In the samples treated with a pressure of 120 MPa and 150 MPa, 59% and 67% more β-glucan content was released, respectively. The soluble dietary fiber content in the samples treated with a pressure of 120 MPa and 150 MPa was increased by 44.8% and 43.2%, respectively, compared with the sample treated with a pressure of 0 MPa. From the perspective of the relative stability of the sample and nutrient enhancement, the processing pressure of 120 MPa was a good choice. In addition, DHPM processing effectively reduced the average particle size and the relaxation time of the water molecules of whole-grain oat pulp, whereas it increased the apparent viscosity of whole-grain oat pulp; all of the above changes alleviated the gravitational subsidence of particles to a certain extent, and thus the overall stability of the system was improved. Furthermore, CLSM and AFM showed that the samples OM-120 and OM-150 had a more uniform and stable structural system as a whole. This study could provide theoretical guidance for the development of a whole-grain oat drink with improved quality and consistency.

Genetic Analysis and Physical Mapping of Oat Adult Plant Resistance Loci Against Puccinia coronata f. sp. avenae

Six quantitative trait loci (QTLs) for adult plant resistance against oat crown rust (Puccinia coronata f. sp. avenae) were identified from mapping three recombinant inbred populations. Using genotyping-by-sequencing with markers called against the OT3098 v1 reference genome, the QTLs were mapped on six different chromosomes: Chr1D, Chr4D, Chr5A, Chr5D, Chr7A, and Chr7C. Composite interval mapping with marker cofactor selection showed that the phenotypic variance explained by all identified QTLs for coefficient of infection range from 12.2 to 46.9%, whereas heritability estimates ranged from 0.11 to 0.38. The significant regions were narrowed down to intervals of 3.9 to 25 cM, equivalent to physical distances of 11 to 133 Mb. At least two flanking single-nucleotide polymorphism markers were identified within 10 cM of each QTL that could be used in marker-assisted introgression, pyramiding, and selection. The additive effects of the QTLs in each population were determined using single-nucleotide polymorphism haplotype data, which showed a significantly lower coefficient of infection in lines homozygous for the resistant alleles. Analysis of pairwise linkage disequilibrium also revealed high correlation of markers and presence of linkage blocks in the significant regions. To further facilitate marker-assisted breeding, polymerase chain reaction allelic competitive extension (PACE) markers for the adult plant resistance loci were developed. Putative candidate genes were also identified in each of the significant regions, which include resistance gene analogs that encode for kinases, ligases, and predicted receptors of avirulence proteins from pathogens.

Genome-wide analysis and expression profiling of glyoxalase gene families in oat (Avena sativa) indicate their responses to abiotic stress during seed germination

Abiotic stresses have deleterious effects on seed germination and seedling establishment, leading to significant crop yield losses. Adverse environmental conditions can cause the accumulation of methylglyoxal (MG) within plant cells, which can negatively impact plant growth and development. The glyoxalase system, which consists of the glutathione (GSH)-dependent enzymes glyoxalase I (GLX1) and glyoxalase II (GLX2), as well as the GSH-independent glyoxalase III (GLX3 or DJ-1), plays a crucial role in detoxifying MG. However, genome-wide analysis of glyoxalase genes has not been performed for one of the agricultural important species, oat (Avena sativa). This study identified a total of 26 AsGLX1 genes, including 8 genes encoding Ni²⁺-dependent GLX1s and 2 genes encoding Zn²⁺-dependent GLX1s. Additionally, 14 AsGLX2 genes were identified, of which 3 genes encoded proteins with both lactamase B and hydroxyacylglutathione hydrolase C-terminal domains and potential catalytic activity, and 15 AsGLX3 genes encoding proteins containing double DJ-1 domains. The domain architecture of the three gene families strongly correlates with the clades observed in the phylogenetic trees. The AsGLX1, AsGLX2, and AsGLX3 genes were evenly distributed in the A, C, and D subgenomes, and gene duplication of AsGLX1 and AsGLX3 genes resulted from tandem duplications. Besides the core cis-elements, hormone responsive elements dominated the promoter regions of the glyoxalase genes, and stress responsive elements were also frequently observed. The subcellular localization of glyoxalases was predicted to be primarily in the cytoplasm, chloroplasts, and mitochondria, with a few presents in the nucleus, which is consistent with their tissue-specific expression. The highest expression levels were observed in leaves and seeds, indicating that these genes may play important roles in maintaining leaf function and ensuring seed vigor. Moreover, based on in silico predication and expression pattern analysis, AsGLX1-7A, AsGLX2-5D, AsDJ-1-5D, AsGLX1-3D2, and AsGLX1-2A were suggested as promising candidate genes for improving stress resistance or seed vigor in oat. Overall, the identification and analysis of the glyoxalase gene families in this study can provide new strategies for improving oat stress resistance and seed vigor.

Genome-Wide Identification of NAC Gene Family and Expression Analysis under Abiotic Stresses in Avena sativa

In this study, a total of 177 NAC members were identified in Avena sativa, located on 21 chromosomes. Phylogenetic analysis showed that AsNAC proteins could be divided into seven subfamilies (I-VII), and that proteins in the same subfamily have similar protein motifs. Gene structure analysis found that NAC introns ranged from 1 to 17. Cis-element analysis of the promoter indicated that the gene family may have stress-related elements and growth regulation elements. Through qRT-PCR experiments, we speculated that AsNACs genes can respond to abiotic stresses such as cold, freezing, salt, and saline alkali. This study provides a theoretical basis for further exploring the function of the NAC gene family in A. sativa.

The mitochondrial genome of the diploid oat Avena longiglumis

BACKGROUND

Avena longiglumis Durieu (2n = 2x = 14) is a wild relative of cultivated oat (Avena sativa, 2n = 6x = 42) with good agronomic and nutritional traits. The plant mitochondrial genome has a complex organization and carries genetic traits of value in exploiting genetic resources, not least male sterility alleles used to generate F₁ hybrid seeds. Therefore, we aim to complement the chromosomal-level nuclear and chloroplast genome assemblies of A. longiglumis with the complete assembly of the mitochondrial genome (mitogenome) based on Illumina and ONT long reads, comparing its structure with Poaceae species.

RESULTS

The complete mitochondrial genome of A. longiglumis can be represented by one master circular genome being 548,445 bp long with a GC content of 44.05%. It can be represented by linear or circular DNA molecules (isoforms or contigs), with multiple alternative configurations mediated by long (4,100-31,235 bp) and medium (144-792 bp) size repeats. Thirty-five unique protein-coding genes, three unique rRNA genes, and 11 unique tRNA genes are identified. The mitogenome is rich in duplications (up to 233 kb long) and multiple tandem or simple sequence repeats, together accounting for more than 42.5% of the total length. We identify homologous sequences between the mitochondrial, plastid and nuclear genomes, including the exchange of eight plastid-derived tRNA genes, and nuclear-derived retroelement fragments. At least 85% of the mitogenome is duplicated in the A. longiglumis nuclear genome. We identify 269 RNA editing sites in mitochondrial protein-coding genes including stop codons truncating ccmFC transcripts.

CONCLUSIONS

Comparative analysis with Poaceae species reveals the dynamic and ongoing evolutionary changes in mitochondrial genome structure and gene content. The complete mitochondrial genome of A. longiglumis completes the last link of the oat reference genome and lays the foundation for oat breeding and exploiting the biodiversity in the genus.

Genome-wide-association study and transcriptome analysis reveal the genetic basis controlling the formation of leaf wax in Brassica napus

Cuticular wax protects plants from various biotic and abiotic stresses. However, the genetic network of wax biosynthesis and the environmental factors influencing leaf wax production in rapeseed (Brassica napus) remains unclear. Here, we demonstrated the role of leaf wax in the resistance to Sclerotinia infection in rapeseed. We found that leaves grown under high light intensity had higher expression of genes involved in wax biosynthesis, and produced more wax on the leaf surface, compared with those grown under low light conditions. Genome-wide association study (GWAS) identified 89 single nucleotide polymorphisms significantly associated with leaf wax coverage. A cross-analysis between GWAS and differentially expressed genes (DEGs) in the leaf epidermis of the accessions with contrasting differences in wax content revealed 17 candidate genes that control this variation in rapeseed. Selective sweep analysis combined with DEG analysis unveiled 510 candidate genes with significant selective signatures. From the candidate genes, we selected BnaA02.LOX4, a putative lipoxygenase, and BnaCnn.CER1, BnaA02.CER3, BnaC02.CER3, and BnaA01.CER4 (ECERIFERUM1-4) that were putatively responsible for wax biosynthesis, to analyse the allelic forms and haplotypes corresponding to high or low leaf wax coverage. These data enrich our knowledge about wax formation, and provide a gene pool for breeding an ideal leaf wax content in rapeseed.

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.

Avenanthramides, Distinctive Hydroxycinnamoyl Conjugates of Oat, Avena sativa L.: An Update on the Biosynthesis, Chemistry, and Bioactivities

Avenanthramides are a group of N-cinnamoylanthranilic acids (phenolic alkaloid compounds) that are produced in oat plants as phytoalexins, in response to pathogen attack and elicitation. The enzyme catalysing the cinnamamide-generating reaction is hydroxycinnamoyl-CoA: hydroxyanthranilate N-hydroxycinnamoyltransferase (HHT, a member of the super family of BAHD acyltransferases). HHT from oat appears to have a narrow range of substrate usage, with preferred use of 5-hydroxyanthranilic acid (and to a lesser extent, other hydroxylated and methoxylated derivatives) as acceptor molecules, but is able to use both substituted cinnamoyl-CoA and avenalumoyl-CoA thioesters as donor molecules. Avenanthramides thus combine carbon skeletons from both the stress-inducible shikimic acid and phenylpropanoid pathways. These features contribute to the chemical characteristics of avenanthramides as multifunctional plant defence compounds, as antimicrobial agents and anti-oxidants. Although avenanthramides are naturally and uniquely synthesised in oat plants, these molecules also exhibit medicinal and pharmaceutical uses important for human health, prompting research into utilisation of biotechnology to enhance agriculture and value-added production.

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.

Genetic mapping of the powdery mildew resistance gene Pm7 on oat chromosome 5D

Three independent experiments with different genetic backgrounds mapped the resistance gene Pm7 in the oat genome to the distal part of the long arm of chromosome 5D. Resistance of oat to Blumeria graminis DC. f. sp. avenae is an important breeding goal in Central and Western Europe. In this study, the position of the effective and widely used resistance gene Pm7 in the oat genome was determined based on three independent experiments with different genetic backgrounds: genome-wide association mapping in a diverse set of inbred oat lines and binary phenotype mapping in two bi-parental populations. Powdery mildew resistance was assessed in the field as well as by detached leaf tests in the laboratory. Genotyping-by-sequencing was conducted to establish comprehensive genetic fingerprints for subsequent genetic mapping experiments. All three mapping approaches located the gene to the distal part of the long arm of chromosome 5D in the hexaploid oat genome sequences of OT3098 and 'Sang.' Markers from this region were homologous to a region of chromosome 2Ce of the C-genome species, Avena eriantha, the donor of Pm7, which appears to be the ancestral source of a translocated region on the hexaploid chromosome 5D.

Changes in physio-biochemical parameters and expression of metallothioneins in Avena sativa L. in response to drought

Drought is one of the major threats to food security. Among several mechanisms involved in plant stress tolerance, one protein family-the plant metallothioneins (MTs)-shows great promise for enhancing drought resistance. Plant metallothioneins in oat (Avena sativa L.) have not yet been deeply analysed, and the literature lacks a comprehensive study of the whole family of plant MTs in response to drought. In this study, we showed that the number and nature of cis-elements linked with stress response in promoters of AsMTs1-3 differed depending on the MT type. Drought stress in oat plants caused an increase in the expression of AsMT2 and AsMT3 and a decrease in the expression of AsMT1 compared to well-watered plants. Moreover, the low values of relative water content, water use efficiency, net photosynthesis (P_N), transpiration (E), stomatal conductance (g_s), chlorophyll a, and carotenoid were accompanied by high levels of electrolyte leakage, internal CO₂ concentration (C_i) and abscisic acid content, and high activity of antioxidants enzymes in plants under drought stress. The present study puts forward the idea that AsMTs are crucial for oat response to drought stress not only by regulating antioxidant activity but also by changing the plant water regime and photosynthesis. Our results support the hypothesis that structural differences among types of plant MTs reflect their diversified physiological roles.

Biocontrol Effect of Clonostachys rosea on Fusarium graminearum Infection and Mycotoxin Detoxification in Oat (Avena sativa)

Oat (Avena sativa) is susceptible to Fusarium head blight (FHB). The quality of oat grain is threatened by the accumulation of mycotoxins, particularly the trichothecene deoxynivalenol (DON), which also acts as a virulence factor for the main pathogen Fusarium graminearum. The plant can defend itself, e.g., by DON detoxification by UGT-glycosyltransferases (UTGs) and accumulation of PR-proteins, even though these mechanisms do not deliver effective levels of resistance. We studied the ability of the fungal biocontrol agent (BCA) Clonostachys rosea to reduce FHB and mycotoxin accumulation. Greenhouse trials showed that C. rosea-inoculation of oat spikelets at anthesis 3 days prior to F. graminearum inoculation reduced both the amount of Fusarium DNA (79%) and DON level (80%) in mature oat kernels substantially. DON applied to C. rosea-treated spikelets resulted in higher conversion of DON to DON-3-Glc than in mock treated plants. Moreover, there was a significant enhancement of expression of two oat UGT-glycosyltransferase genes in C. rosea-treated oat. In addition, C. rosea treatment activated expression of genes encoding four PR-proteins and a WRKY23-like transcription factor, suggesting that C. rosea may induce resistance in oat. Thus, C. rosea IK726 has strong potential to be used as a BCA against FHB in oat as it inhibits F. graminearum infection effectively, whilst detoxifying DON mycotoxin rapidly.

Development of a set of novel binary expression vectors for plant gene function analysis and genetic transformation

With the advent of multiple omics and Genome-Wide Association Studies (GWAS) technology, genome-scale functional analysis of candidate genes is to be conducted in diverse plant species. Construction of plant binary expression vectors is the prerequisite for gene function analysis. Therefore, it is of significance to develop a set of plant binary expression vectors with highly efficient, inexpensive, and convenient cloning method, and easy-to-use in screening of positive recombinant in Escherichia coli. In this study, we developed a set of plant binary expression vectors, termed pBTR vectors, based on Golden Gate cloning using BsaI restriction site. Foreign DNA fragment of interest (FDI) can be cloned into the destination pBTR by one-step digestion-ligation reaction in a single tube, and even the FDI contains internal BsaI site(s). Markedly, in one digestion-ligation reaction, multiple FDIs (exemplified by cloning four soybean Glyma.02g025400, Glyma.05g201700, Glyma.06g165700, and Glyma.17g095000 genes) can be cloned into the pBTR vector to generate multiple corresponding expression constructs (each expression vector carrying an FDI). In addition, the pBTR vectors carry the visual marker, a brightness monomeric red fluorescent protein mScarlet-I, that can be observed with the unaided eye in screening of positive recombinants without the use of additional reagents/equipment. The reliability of the pBTR vectors was validated in plants by overexpression of AtMyb75/PAP1 in tomato and GUSPlus in soybean roots via Agrobacterium rhizogenes-mediated transformation, promoter activity analysis of AtGCSpro in Arabidopsis via A. tumefaciens-mediated transformation, and protein subcellular localization of the Vitis vinifera VvCEB1_opt in tobacco, respectively. These results demonstrated that the pBTR vectors can be used in analysis of gene (over)expression, promoter activity, and protein subcellular localization. These vectors will contribute to speeding up gene function analysis and the process of plant molecular breeding.