Rapid evolutionary dynamics in a 2.8-Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Physico-Chemical and Rheological Trait-Based Identification of Indian Wheat Varieties Suitable for Different End-Uses

India has increased its wheat production phenomenally in the last two decades and it now has a buffer stock of 9.7 million tonnes. However, despite the release of several wheat cultivars, the end-use quality traits of Indian wheat varieties have not been explored in-depth to determine the increasing demand of the domestic processing industry as well as export. In this study, 55 wheat genotypes including 47 released varieties, and 8 genetic stocks were grown along with 10 Australian varieties grown during cropping seasons: 2019-2020 and 2020-2021 and diversity in different physiochemical and rheological traits was evaluated. They showed considerable diversity in all the quality traits studied. However, very few genotypes could be found suitable for any one end-use. Five genotypes were found to possess four to five traits for superior bread-making quality. Two varieties and three advanced breeding lines had up to four good chapati quality traits. None of the released varieties investigated had suitable traits for biscuit making; however, two breeding lines possessed requisite quality traits suitable for biscuit making. It is, therefore, concluded that systematic breeding efforts are required to develop genotypes that bring together the most important quality traits in a single genotype to be suitable for domestic industry as well as for export.

Intraspecific divergence of diploid grass Aegilopscomosa is associated with structural chromosome changes

Aegilopscomosa Smith in Sibthorp et Smith, 1806 is diploid grass with MM genome constitution occurring mainly in Greece. Two morphologically distinct subspecies - Ae.c.comosa Chennaveeraiah, 1960 and Ae.c.heldreichii (Holzmann ex Boissier) Eig, 1929 are discriminated within Ae.comosa, however, genetic and karyotypic bases of their divergence are not fully understood. We used Fluorescence in situ hybridization (FISH) with repetitive DNA probes and electrophoretic analysis of gliadins to characterize the genome and karyotype of Ae.comosa to assess the level of their genetic diversity and uncover mechanisms leading to radiation of subspecies. We show that two subspecies differ in size and morphology of chromosomes 3M and 6M, which can be due to reciprocal translocation. Subspecies also differ in the amount and distribution of microsatellite and satellite DNA sequences, the number and position of minor NORs, especially on 3M and 6M, and gliadin spectra mainly in the a-zone. Frequent occurrence of hybrids can be caused by open pollination, which, along with genetic heterogeneity of accessions and, probably, the lack of geographic or genetic barrier between the subspecies, may contribute to extremely broad intraspecific variation of GAA_n and gliadin patterns in Ae.comosa, which are usually not observed in endemic plant species.

Genome-Wide Identification and Expression Analysis of Cysteine-Rich Polycomb-like Protein (CPP) Gene Family in Tomato

The cysteine-rich polycomb-like protein (CPP) gene family is a class of transcription factors containing conserved cysteine-rich CRC structural domains that is involved in the regulation of plant growth and stress tolerance to adversity. Relative to other gene families, the CPP gene family has not received sufficient attention. In this study, six SlCPPs were identified for the first time using the most recent genome-wide identification data of tomato. Subsequently, a phylogenetic analysis classified SlCPPs into four subfamilies. The analysis of cis-acting elements in the promoter indicates that SlCPPs are involved in plant growth and development and also stress response. We present for the first time the prediction of the tertiary structure of these SlCPPs proteins using the AlphaFold2 artificial intelligence system developed by the DeepMind team. Transcriptome data analysis showed that SlCPPs were differentially expressed in different tissues. Gene expression profiling showed that all SlCPPs except SlCPP5 were up-regulated under drought stress; SlCPP2, SlCPP3 and SlCPP4 were up-regulated under cold stress; SlCPP2 and SlCPP5 were up-regulated under salt stress; all SlCPPs were up-regulated under inoculation with Cladosporium fulvum; and SlCPP1, SlCPP3, and SlCPP4 were up-regulated under inoculation with Stemphylium lycopersici. We performed a virus-induced gene silencing experiment on SlCPP3, and the results indicated that SlCPP3 was involved in the response to drought stress. Finally, we predicted the interaction network of the key gene SlCPP3, and there was an interaction relationship between SlCPP3 and 10 genes, such as RBR1 and MSI1. The positive outcome showed that SlCPPs responded to environmental stress. This study provides a theoretical and empirical basis for the response mechanisms of tomato in abiotic stresses.

Characterization of Glutenin Genes in Bread Wheat by Third-Generation RNA Sequencing and the Development of a Glu-1Dx5 Marker Specific for the Extra Cysteine Residue

High-molecular-weight glutenin subunits (HMW-GS) and low-molecular-weight glutenin subunits (LMW-GS) in a mature grain play important roles in the formation of a glutenin macropolymer and gluten quality. To characterize the expressed glutenin genes of the bread wheat variety Xinmai 26 during seed development, a total of 18 full-length transcripts were obtained by the newly emerged third-generation RNA sequencing of the PacBio Sequel II platform, including 5 transcripts of HMW-GS genes and 13 transcripts of LMW-GS genes (8 intact genes and 5 pseudogenes). Combined with the patterns of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS), allelic types of the obtained glutenin genes were, respectively, determined, wherein molecular characterization deduced by transcript1528 (1Dx5) and transcript907 (Glu-A3c) indicated their great influence on dough quality. In addition, a specific functional marker dCAPS5 was developed for the single-nucleotide substitution at position 353 of the 1Dx5 subunit, which was further intensively compared with the other proposed markers to efficiently utilize the 1Dx5 subunit with the extra cysteine residue. This study provides an efficient method to accurately identify and utilize glutenin genes in bread wheat, which is helpful in understanding the contributions of glutenin genes to wheat quality.

Genome-Wide Identification, Characterization and Expression Pattern Analysis of the γ-Gliadin Gene Family in the Durum Wheat (Triticum durum Desf.) Cultivar Svevo

Very recently, the genome of the modern durum wheat cv. Svevo was fully sequenced, and its assembly is publicly available. So, we exploited the opportunity to carry out an in-depth study for the systematic characterization of the γ-gliadin gene family in the cv. Svevo by combining a bioinformatic approach with transcript and protein analysis. We found that the γ-gliadin family consists of nine genes that include seven functional genes and two pseudogenes. Three genes, Gli-γ1a, Gli-γ3a and Gli-γ4a, and the pseudogene Gli-γ2a* mapped on the A genome, whereas the remaining four genes, Gli-γ1b, Gli-γ2b, Gli-γ3b and Gli-γ5b, and the pseudogene Gli-γ4b* mapped on the B genome. The functional γ-gliadins presented all six domains and eight-cysteine residues typical of γ-gliadins. The Gli-γ1b also presented an additional cysteine that could possibly have a role in the formation of the gluten network through binding to HMW glutenins. The γ-gliadins from the A and B genome differed in their celiac disease (CD) epitope content and composition, with the γ-gliadins from the B genome showing the highest frequency of CD epitopes. In all the cases, almost all the CD epitopes clustered in the central region of the γ-gliadin proteins. Transcript analysis during seed development revealed that all the functional γ-gliadin genes were expressed with a similar pattern, although significant differences in the transcript levels were observed among individual genes that were sometimes more than 60-fold. A progressive accumulation of the γ-gliadin fraction was observed in the ripening seeds that reached 34% of the total gliadin fraction at harvest maturity. We believe that the insights generated in the present study could aid further studies on gliadin protein functions and future breeding programs aimed at the selection of new healthier durum wheat genotypes.

Single-parent expression drives dynamic gene expression complementation in maize hybrids

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Genomic and functional genomics analyses of gluten proteins and prospect for simultaneous improvement of end-use and health-related traits in wheat

KEY MESSAGE

Recent genomic and functional genomics analyses have substantially improved the understanding on gluten proteins, which are important determinants of wheat grain quality traits. The new insights obtained and the availability of precise, versatile and high-throughput genome editing technologies will accelerate simultaneous improvement of wheat end-use and health-related traits. Being a major staple food crop in the world, wheat provides an indispensable source of dietary energy and nutrients to the human population. As worldwide population grows and living standards rise in both developed and developing countries, the demand for wheat with high quality attributes increases globally. However, efficient breeding of high-quality wheat depends on critically the knowledge on gluten proteins, which mainly include several families of prolamin proteins specifically accumulated in the endospermic tissues of grains. Although gluten proteins have been studied for many decades, efficient manipulation of these proteins for simultaneous enhancement of end-use and health-related traits has been difficult because of high complexities in their expression, function and genetic variation. However, recent genomic and functional genomics analyses have substantially improved the understanding on gluten proteins. Therefore, the main objective of this review is to summarize the genomic and functional genomics information obtained in the last 10 years on gluten protein chromosome loci and genes and the cis- and trans-factors regulating their expression in the grains, as well as the efforts in elucidating the involvement of gluten proteins in several wheat sensitivities affecting genetically susceptible human individuals. The new insights gathered, plus the availability of precise, versatile and high-throughput genome editing technologies, promise to speed up the concurrent improvement of wheat end-use and health-related traits and the development of high-quality cultivars for different consumption needs.

Robust non-syntenic gene expression patterns in diverse maize hybrids during root development

Distantly related maize (Zea mays L.) inbred lines exhibit an exceptional degree of structural genomic diversity, which is probably unique among plants. This study systematically investigated the developmental and genotype-dependent regulation of the primary root transcriptomes of a genetically diverse panel of maize F1-hybrids and their parental inbred lines. While we observed substantial transcriptomic changes during primary root development, we demonstrated that hybrid-associated gene expression patterns, including differential, non-additive, and allele-specific transcriptome profiles, are particularly robust to these developmental fluctuations. For instance, differentially expressed genes with preferential expression in hybrids were highly conserved during development in comparison to their parental counterparts. Similarly, in hybrids a major proportion of non-additively expressed genes with expression levels between the parental values were particularly conserved during development. Importantly, in these expression patterns non-syntenic genes that evolved after the separation of the maize and sorghum lineages were systemically enriched. Furthermore, non-syntenic genes were substantially linked to the conservation of all surveyed gene expression patterns during primary root development. Among all F1-hybrids, between ~40% of the non-syntenic genes with unexpected allelic expression ratios and ~60% of the non-syntenic differentially and non-additively expressed genes were conserved and therefore robust to developmental changes. Hence, the enrichment of non-syntenic genes during primary root development might be involved in the developmental adaptation of maize roots and thus the superior performance of hybrids.

Outlook for coeliac disease patients: towards bread wheat with hypoimmunogenic gluten by gene editing of α- and γ-gliadin gene families

BACKGROUND

Wheat grains contain gluten proteins, which harbour immunogenic epitopes that trigger Coeliac disease in 1-2% of the human population. Wheat varieties or accessions containing only safe gluten have not been identified and conventional breeding alone struggles to achieve such a goal, as the epitopes occur in gluten proteins encoded by five multigene families, these genes are partly located in tandem arrays, and bread wheat is allohexaploid. Gluten immunogenicity can be reduced by modification or deletion of epitopes. Mutagenesis technologies, including CRISPR/Cas9, provide a route to obtain bread wheat containing gluten proteins with fewer immunogenic epitopes.

RESULTS

In this study, we analysed the genetic diversity of over 600 α- and γ-gliadin gene sequences to design six sgRNA sequences on relatively conserved domains that we identified near coeliac disease epitopes. They were combined in four CRISPR/Cas9 constructs to target the α- or γ-gliadins, or both simultaneously, in the hexaploid bread wheat cultivar Fielder. We compared the results with those obtained with random mutagenesis in cultivar Paragon by γ-irradiation. For this, Acid-PAGE was used to identify T1 grains with altered gliadin protein profiles compared to the wild-type endosperm. We first optimised the interpretation of Acid-PAGE gels using Chinese Spring deletion lines. We then analysed the changes generated in 360 Paragon γ-irradiated lines and in 117 Fielder CRISPR/Cas9 lines. Similar gliadin profile alterations, with missing protein bands, could be observed in grains produced by both methods.

CONCLUSIONS

The results demonstrate the feasibility and efficacy of using CRISPR/Cas9 to simultaneously edit multiple genes in the large α- and γ-gliadin gene families in polyploid bread wheat. Additional methods, generating genomics and proteomics data, will be necessary to determine the exact nature of the mutations generated with both methods.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA

BACKGROUND

Omega-5 gliadins are a group of highly repetitive gluten proteins in wheat flour encoded on the 1B chromosome of hexaploid wheat. These proteins are the major sensitizing allergens in a severe form of food allergy called wheat-dependent exercise-induced anaphylaxis (WDEIA). The elimination of omega-5 gliadins from wheat flour through biotechnology or breeding approaches could reduce the immunogenic potential and adverse health effects of the flour.

RESULTS

A mutant line missing low-molecular weight glutenin subunits encoded at the Glu-B3 locus was selected previously from a doubled haploid population generated from two Korean wheat cultivars. Analysis of flour from the mutant line by 2-dimensional gel electrophoresis coupled with tandem mass spectrometry revealed that the omega-5 gliadins and several gamma gliadins encoded by the closely linked Gli-B1 locus were also missing as a result of a deletion of at least 5.8 Mb of chromosome 1B. Two-dimensional immunoblot analysis of flour proteins using sera from WDEIA patients showed reduced IgE reactivity in the mutant relative to the parental lines due to the absence of the major omega-5 gliadins. However, two minor proteins showed strong reactivity to patient sera in both the parental and the mutant lines and also reacted with a monoclonal antibody against omega-5 gliadin. Analysis of the two minor reactive proteins by mass spectrometry revealed that both proteins correspond to omega-5 gliadin genes encoded on chromosome 1D that were thought previously to be pseudogenes.

CONCLUSIONS

While breeding approaches can be used to reduce the levels of the highly immunogenic omega-5 gliadins in wheat flour, these approaches are complicated by the genetic linkage of different classes of gluten protein genes and the finding that omega-5 gliadins may be encoded on more than one chromosome. The work illustrates the importance of detailed knowledge about the genomic regions harboring the major gluten protein genes in individual wheat cultivars for future efforts aimed at reducing the immunogenic potential of wheat flour.

Dynamic Evolution of α-Gliadin Prolamin Gene Family in Homeologous Genomes of Hexaploid Wheat

Wheat Gli-2 loci encode complex groups of α-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of α-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous α-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 α-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of α-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of α-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of α-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.

Single-Parent Expression Is a General Mechanism Driving Extensive Complementation of Non-syntenic Genes in Maize Hybrids

Maize (Zea mays L.) displays an exceptional degree of structural genomic diversity [1, 2]. In addition, variation in gene expression further contributes to the extraordinary phenotypic diversity and plasticity of maize. This study provides a systematic investigation on how distantly related homozygous maize inbred lines affect the transcriptomic plasticity of their highly heterozygous F₁ hybrids. The classical dominance model of heterosis explains the superiority of hybrid plants by the complementation of deleterious parental alleles by superior alleles of the second parent at many loci [3]. Genes active in one inbred line but inactive in another represent an extreme instance of allelic diversity defined as single-parent expression [4]. We observed on average ∼1,000 such genes in all inbred line combinations during primary root development. These genes consistently displayed expression complementation (i.e., activity) in their hybrid progeny. Consequently, extreme expression complementation is a general mechanism that results on average in ∼600 additionally active genes and their encoded biological functions in hybrids. The modern maize genome is complemented by a set of non-syntenic genes, which emerged after the separation of the maize and sorghum lineages and lack syntenic orthologs in any other grass species [5]. We demonstrated that non-syntenic genes are the driving force of gene expression complementation in hybrids. Among those, the highly diversified families of bZIP and bHLH transcription factors [6] are systematically overrepresented. In summary, extreme gene expression complementation extensively shapes the transcriptomic plasticity of maize hybrids and might therefore be one factor controlling the developmental plasticity of hybrids.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

Genome sequence of the progenitor of the wheat D genome Aegilops tauschii

A combination of advanced sequencing and mapping techniques is used to produce a reference genome of Aegilops tauschii, progenitor of the wheat D genome, providing a valuable resource for comparative genetic studies.

Sequencing the genomes of crops plants provides useful resources for crop improvement and breeding. Jan Dvořák, Katrien Devos, Steven Salzberg and colleagues report a reference genome for Aegilops tauschii, the diploid progenitor of the D genome of hexaploid wheat. They use a combination of ordered-clone genome sequencing, whole-genome shotgun sequencing and BioNano optical genome mapping to assemble this large and highly repetitive genome. This provides a useful resource for comparative genomics studies of wheat.

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat¹ (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat^2,3,4. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence⁵. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

New insights into structural organization and gene duplication in a 1.75-Mb genomic region harboring the α-gliadin gene family in Aegilops tauschii, the source of wheat D genome

Among the wheat prolamins important for its end-use traits, α-gliadins are the most abundant, and are also a major cause of food-related allergies and intolerances. Previous studies of various wheat species estimated that between 25 and 150 α-gliadin genes reside in the Gli-2 locus regions. To better understand the evolution of this complex gene family, the DNA sequence of a 1.75-Mb genomic region spanning the Gli-2 locus was analyzed in the diploid grass, Aegilops tauschii, the ancestral source of D genome in hexaploid bread wheat. Comparison with orthologous regions from rice, sorghum, and Brachypodium revealed rapid and dynamic changes only occurring to the Ae. tauschii Gli-2 region, including insertions of high numbers of non-syntenic genes and a high rate of tandem gene duplications, the latter of which have given rise to 12 copies of α-gliadin genes clustered within a 550-kb region. Among them, five copies have undergone pseudogenization by various mutation events. Insights into the evolutionary relationship of the duplicated α-gliadin genes were obtained from their genomic organization, transcription patterns, transposable element insertions and phylogenetic analyses. An ancestral glutamate-like receptor (GLR) gene encoding putative amino acid sensor in all four grass species has duplicated only in Ae. tauschii and generated three more copies that are interspersed with the α-gliadin genes. Phylogenetic inference and different gene expression patterns support functional divergence of the Ae. tauschii GLR copies after duplication. Our results suggest that the duplicates of α-gliadin and GLR genes have likely taken different evolutionary paths; conservation for the former and neofunctionalization for the latter.

Genome-Wide Association and Prediction of Grain and Semolina Quality Traits in Durum Wheat Breeding Populations

Grain yield and semolina quality traits are essential selection criteria in durum wheat breeding. However, high phenotypic screening costs limit selection to relatively few breeding lines in late generations. This selection paradigm confers relatively low selection efficiency due to the advancement of undesirable lines into expensive yield trials for grain yield and quality trait testing. Marker-aided selection can enhance selection efficiency, especially for traits that are difficult or costly to phenotype. The aim of this study was to identify major quality trait quantitative trait loci (QTL) for marker-assisted selection (MAS) and to explore potential application of genomic selection (GS) in a durum wheat breeding program. In this study, genome-wide association mapping was conducted for five quality traits using 1184 lines from the North Dakota State University (NDSU) durum wheat breeding program. Several QTL associated with test weight, semolina color, and gluten strength were identified. Genomic selection models were developed and forward prediction accuracies of 0.27 to 0.66 were obtained for the five quality traits. Our results show the potential for grain and semolina quality traits to be selected more efficiently through MAS and GS with further refinement. Considerable opportunity exists to extend these techniques to other traits such as grain yield and agronomic characteristics, further improving breeding efficiency in durum cultivar development.

Genome-wide analysis of complex wheat gliadins, the dominant carriers of celiac disease epitopes

Gliadins, specified by six compound chromosomal loci (Gli-A1/B1/D1 and Gli-A2/B2/D2) in hexaploid bread wheat, are the dominant carriers of celiac disease (CD) epitopes. Because of their complexity, genome-wide characterization of gliadins is a strong challenge. Here, we approached this challenge by combining transcriptomic, proteomic and bioinformatic investigations. Through third-generation RNA sequencing, full-length transcripts were identified for 52 gliadin genes in the bread wheat cultivar Xiaoyan 81. Of them, 42 were active and predicted to encode 25 α-, 11 γ-, one δ- and five ω-gliadins. Comparative proteomic analysis between Xiaoyan 81 and six newly-developed mutants each lacking one Gli locus indicated the accumulation of 38 gliadins in the mature grains. A novel group of α-gliadins (the CSTT group) was recognized to contain very few or no CD epitopes. The δ-gliadins identified here or previously did not carry CD epitopes. Finally, the mutant lacking Gli-D2 showed significant reductions in the most celiac-toxic α-gliadins and derivative CD epitopes. The insights and resources generated here should aid further studies on gliadin functions in CD and the breeding of healthier wheat.