Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Identification of the causal mutation in early heading mutant of bread wheat (Triticum aestivum L.) using MutMap approach.

In bread wheat (Triticum aestivum L.), fine-tuning the heading time is essential to maximize grain yield. Photoperiod-1 (Ppd-1) and VERNALIZATION 1 (Vrn-1) are major genes affecting photoperiod sensitivity and vernalization requirements, respectively. These genes have predominantly governed heading timing. However, Ppd-1 and Vrn-1 significantly impact heading dates, necessitating another gene that can slightly modify heading dates for fine-tuning. In this study, we developed an early heading mutant from the ethyl methanesulfonate-mutagenized population of the Japanese winter wheat cultivar "Kitahonami." MutMap analysis identified a nonsense mutation in the clock component gene Wheat PHYTOCLOCK 1/LUX ARRHYTHMO (WPCL-D1) as the probable SNP responsible for the early heading mutant on chromosome 3D. Segregation analysis using F2 and F3 populations confirmed that plants carrying the wpcl-D1 allele headed significantly earlier than those with the functional WPCL-D1. The early heading mutant exhibited increased expression levels of Ppd-1 and circadian clock genes, such as WPCL1 and LATE ELONGATED HYPOCOTYL (LHY). Notably, the transcript accumulation levels of Ppd-A1 and Ppd-D1 were influenced by the copy number of the functional WPCL1 gene. These results suggest that a loss-of-function mutation in WPCL-D1 is the causal mutation for the early heading phenotype. Adjusting the functional copy number of WPCL1 will be beneficial in fine-tuning of heading dates. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01478-5.

Structural Dynamics of E6AP E3 Ligase HECT Domain and Involvement of a Flexible Hinge Loop in the Ubiquitin Chain Synthesis Mechanism.

Ubiquitin (Ub) ligases E3 are important factors in selecting target proteins for ubiquitination and determining the type of polyubiquitin chains on the target proteins. In the HECT (homologous to E6AP C-terminus)-type E3 ligases, the HECT domain is composed of an N-lobe and a C-lobe that are connected by a flexible hinge loop. The large conformational rearrangement of the HECT domain via the flexible hinge loop is essential for the HECT-type E3-mediated Ub transfer from E2 to a target protein. However, detailed insights into the structural dynamics of the HECT domain remain unclear. Here, we provide the first direct demonstration of the structural dynamics of the HECT domain using high-speed atomic force microscopy at the nanoscale. We also found that the flexibility of the hinge loop has a great impact not only on its structural dynamics but also on the formation mechanism of free Ub chains.

Natural variations of wheat EARLY FLOWERING 3 highlight their contributions to local adaptation through fine-tuning of heading time.

KEY MESSAGE: We identified a large chromosomal deletion containing TaELF-B3 that confers early flowering in wheat. This allele has been preferred in recent wheat breeding in Japan to adapt to the environment. Heading at the appropriate time in each cultivation region can greatly contribute to stabilizing and maximizing yield. Vrn-1 and Ppd-1 are known as the major genes for vernalization requirement and photoperiod sensitivity in wheat. Genotype combinations of Vrn-1 and Ppd-1 can explain the variation in heading time. However, the genes that can explain the remaining variations in heading time are largely unknown. In this study, we aimed to identify the genes conferring early heading using doubled haploid lines derived from Japanese wheat varieties. Quantitative trait locus (QTL) analysis revealed a significant QTL on the long arm of chromosome 1B in multiple growing seasons. Genome sequencing using Illumina short reads and Pacbio HiFi reads revealed a large deletion of a ~ 500 kb region containing TaELF-B3, an orthologue of Arabidopsis clock gene EARLY FLOWERING 3 (ELF3). Plants with the deleted allele of TaELF-B3 (ΔTaELF-B3 allele) headed earlier only under short-day vernalization conditions. Higher expression levels of clock- and clock-output genes, such as Ppd-1 and TaGI, were observed in plants with the ΔTaELF-B3 allele. These results suggest that the deletion of TaELF-B3 causes early heading. Of the TaELF-3 homoeoalleles conferring early heading, the ΔTaELF-B3 allele showed the greatest effect on the early heading phenotype in Japan. The higher allele frequency of the ΔTaELF-B3 allele in western Japan suggests that the ΔTaELF-B3 allele was preferred during recent breeding to adapt to the environment. TaELF-3 homoeologs will help to expand the cultivated area by fine-tuning the optimal timing of heading in each environment.

Mutations within the miR172 target site of wheat AP2 homoeologs regulate lodicule size and rachis internode length.

Closed fertilization in flowers, or cleistogamy, reduces the risk of fungal infection in Triticeae crops. In barley (Hordeum vulgare), cleistogamy is determined by a single recessive gene, cly1, which results from a single nucleotide polymorphism within the microRNA172 target site of the Apetala2 (AP2) transcription factor gene. The recessive cly1 allele negatively regulates the development of lodicules, keeping florets closed at anthesis. However, cleistogamy is not evident in hexaploid wheat (Triticum aestivum) cultivars. This study aimed at identifying mutations in wheat AP2 orthologs by ethyl methane sulfonate-induced mutagenesis and high-resolution melt analysis. Although flowers of AP2 mutants induced in the A and D genomes opened at anthesis, their lodicule size was significantly smaller, especially in the direction of depth, than that of wild-type plants. One of the mutants that carried a nucleotide replacement in AP2 from the D genome produced a compact spike caused by a substantial decrease in rachis internode length, analogous to the barley dense spike. Cleistogamous hexaploid wheat might be generated by combining effective mutant alleles of AP2-homoeologous genes.

Genome sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants.

BACKGROUND: Gamma-irradiated mutants of Triticum aestivum L., hexaploid wheat, provide novel and agriculturally important traits and are used as breeding materials. However, the identification of causative genomic regions of mutant phenotypes is challenging because of the large and complicated genome of hexaploid wheat. Recently, the combined use of high-quality reference genome sequences of common wheat and cost-effective resequencing technologies has made it possible to evaluate genome-wide polymorphisms, even in complex genomes. RESULTS: To investigate whether the genome sequencing approach can effectively detect structural variations, such as deletions, frequently caused by gamma irradiation, we selected a grain-hardness mutant from the gamma-irradiated population of Japanese elite wheat cultivar "Kitahonami." The Hardness (Ha) locus, including the puroindoline protein-encoding genes Pina-D1 and Pinb-D1 on the short arm of chromosome 5D, primarily regulates the grain hardness variation in common wheat. We performed short-read genome sequencing of wild-type and grain-hardness mutant plants, and subsequently aligned their short reads to the reference genome of the wheat cultivar "Chinese Spring." Genome-wide comparisons of depth-of-coverage between wild-type and mutant strains detected ~ 130 Mbp deletion on the short arm of chromosome 5D in the mutant genome. Molecular markers for this deletion were applied to the progeny populations generated by a cross between the wild-type and the mutant. A large deletion in the region including the Ha locus was associated with the mutant phenotype, indicating that the genome sequencing is a powerful and efficient approach for detecting a deletion marker of a gamma-irradiated mutant phenotype. In addition, we investigated a pre-harvest sprouting tolerance mutant and identified a 67.8 Mbp deletion on chromosome 3B where Viviparous-B1 and GRAS family transcription factors are located. Co-dominant markers designed to detect the deletion-polymorphism confirmed the association with low germination rate, leading to pre-harvest sprouting tolerance. CONCLUSIONS: Short read-based genome sequencing of gamma-irradiated mutants facilitates the identification of large deletions linked to mutant phenotypes when combined with segregation analyses in progeny populations. This method allows effective application of mutants with agriculturally important traits in breeding using marker-assisted selection.

Allelic variations of Vrn-1 and Ppd-1 genes in Japanese wheat varieties reveal the genotype-environment interaction for heading time.

The timing of heading is largely affected by environmental conditions. In wheat, Vrn-1 and Ppd-1 have been identified as the major genes involved in vernalization requirement and photoperiod sensitivity, respectively. To compare the effects of Vrn-1 and Ppd-1 alleles on heading time under different environments, we genotyped Vrn-1 and Ppd-1 homoeologues and measured the heading time at Morioka, Tsukuba and Chikugo in Japan for two growing seasons. A total of 128 Japanese and six foreign varieties, classified into four populations based on the 519 genome-wide SNPs, were used for analysis. Varieties with the spring alleles (Vrn-D1a or Vrn-D1b) at the Vrn-D1 locus and insensitive allele (Hapl-I) at the Ppd-D1 locus were found in earlier heading varieties. The effects of Vrn-D1 and Ppd-D1 on heading time were stronger than those of the other Vrn-1 and Ppd-1 homoeologues. Analysis of variance revealed that heading time was significantly affected by the genotype-environment interactions. Some Vrn-1 and Ppd-1 alleles conferred earlier or later heading in specific environments, indicating that the effect of both alleles on the timing of heading depends on the environment. Information on Vrn-1 and Ppd-1 alleles, together with heading time in various environments, provide useful information for wheat breeding.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Structural features of two major nucleolar organizer regions (NORs), Nor-B1 and Nor-B2, and chromosome-specific rRNA gene expression in wheat.

The reference genome sequence of wheat 'Chinese Spring' (CS) is now available (IWGSC RefSeq v1.0), but the core sequences defining the nucleolar organizer regions (NORs) have not been characterized. We estimated that the total copy number of the rDNA units in the wheat genome is 11 160, of which 30.5%, 60.9% and 8.6% are located on Nor-B1 (1B), Nor-B2 (6B) and other NORs, respectively. The total length of the NORs is estimated to be 100 Mb, corresponding to approximately 10% of the unassembled portion of the genome not represented in RefSeq v1.0. Four subtypes (S1-S4) of the rDNA units were identified based on differences within the 3' external transcribed spacer regions in Nor-B1 and Nor-B2, and quantitative PCR indicated locus-specific variation in rDNA subtype contents. Expression analyses of rDNA subtypes revealed that S1 was predominantly expressed and S2 weakly expressed, in contrast to the relative abundance of rDNA subtypes in the wheat genome. These results suggest a regulation mechanism of differential rDNA expression based on sequence differences. S3 expression increased in the ditelosomic lines Dt1BL and Dt6BL, suggesting that S3 is subjected to chromosome-mediated silencing. Structural differences were detected in the regions surrounding the NOR among homoeologous chromosomes of groups 1 and 6. The adjacent regions distal to the major NORs were expanded compared with their homoeologous counterparts, and the gene density of these expanded regions was relatively low. We provide evidence that these regions are likely to be important for autoregulation of the associated major NORs as well as silencing of minor NORs.

Ubiquitin chain specificities of E6AP E3 ligase and its HECT domain.

Ubiquitination of target proteins is accomplished by isopeptide bond formation between the carboxy group of the C-terminal glycine (Gly) residue of ubiquitin (Ub) and the ɛ-amino group of lysine (Lys) on the target proteins. The formation of an isopeptide bond between Ubs that gives rise to a poly-Ub chain on the target proteins and the types of poly-Ub chains formed depend on which of the seven Lys residues or N-terminal methionine (Met) residue on Ub is used for chain elongation. To understand the linkage specificity mechanism of Ub chains on E3, the previous study established an assay to monitor the formation of a free diubiquitin chain (Ub2 chain synthesis assay) by HECT type E3 ligase. In this study, we investigated Ub2 chain specificity using E6AP HECT domain. We here demonstrate the importance of the N-terminal domain of full length E6AP for Ub2 chain specificity.

Characterization of a mini core collection of Japanese wheat varieties using single-nucleotide polymorphisms generated by genotyping-by-sequencing.

A core collection of Japanese wheat varieties (JWC) consisting of 96 accessions was established based on their passport data and breeding pedigrees. To clarify the molecular basis of the JWC collection, genome-wide single-nucleotide polymorphism (SNP) genotyping was performed using the genotyping-by-sequencing (GBS) approach. Phylogenetic tree and population structure analyses using these SNP data revealed the genetic diversity and relationships among the JWC accessions, classifying them into four groups; "varieties in the Hokkaido area", "modern varieties in the northeast part of Japan", "modern varieties in the southwest part of Japan" and "classical varieties including landraces". This clustering closely reflected the history of wheat breeding in Japan. Furthermore, to demonstrate the utility of the JWC collection, we performed a genome-wide association study (GWAS) for three traits, namely, "days to heading in autumn sowing", "days to heading in spring sowing" and "culm length". We found significantly associated SNP markers with each trait, and some of these were closely linked to known major genes for heading date or culm length on the genetic map. Our study indicates that this JWC collection is a useful set of germplasm for basic and applied research aimed at understanding and utilizing the genetic diversity among Japanese wheat varieties.

A high-resolution physical map integrating an anchored chromosome with the BAC physical maps of wheat chromosome 6B.

BACKGROUND: A complete genome sequence is an essential tool for the genetic improvement of wheat. Because the wheat genome is large, highly repetitive and complex due to its allohexaploid nature, the International Wheat Genome Sequencing Consortium (IWGSC) chose a strategy that involves constructing bacterial artificial chromosome (BAC)-based physical maps of individual chromosomes and performing BAC-by-BAC sequencing. Here, we report the construction of a physical map of chromosome 6B with the goal of revealing the structural features of the third largest chromosome in wheat. RESULTS: We assembled 689 informative BAC contigs (hereafter reffered to as contigs) representing 91% of the entire physical length of wheat chromosome 6B. The contigs were integrated into a radiation hybrid (RH) map of chromosome 6B, with one linkage group consisting of 448 loci with 653 markers. The order and direction of 480 contigs, corresponding to 87% of the total length of 6B, were determined. We also characterized the contigs that contained a part of the nucleolus organizer region or centromere based on their positions on the RH map and the assembled BAC clone sequences. Analysis of the virtual gene order along 6B using the information collected for the integrated map revealed the presence of several chromosomal rearrangements, indicating evolutionary events that occurred on chromosome 6B. CONCLUSIONS: We constructed a reliable physical map of chromosome 6B, enabling us to analyze its genomic structure and evolutionary progression. More importantly, the physical map should provide a high-quality and map-based reference sequence that will serve as a resource for wheat chromosome 6B.

Characterisation of the wheat (Triticum aestivum L.) transcriptome by de novo assembly for the discovery of phosphate starvation-responsive genes: gene expression in Pi-stressed wheat.

BACKGROUND: Phosphorus (P) is an essential macronutrient for plant growth and development. To modulate their P homeostasis, plants must balance P uptake, mobilisation, and partitioning to various organs. Despite the worldwide importance of wheat as a cultivated food crop, molecular mechanisms associated with phosphate (Pi) starvation in wheat remain unclear. To elucidate these mechanisms, we used RNA-Seq methods to generate transcriptome profiles of the wheat variety 'Chinese Spring' responding to 10 days of Pi starvation. RESULTS: We carried out de novo assembly on 73.8 million high-quality reads generated from RNA-Seq libraries. We then constructed a transcript dataset containing 29,617 non-redundant wheat transcripts, comprising 15,047 contigs and 14,570 non-redundant full-length cDNAs from the TriFLDB database. When compared with barley full-length cDNAs, 10,656 of the 15,047 contigs were unalignable, suggesting that many might be distinct from barley transcripts. The average expression level of the contigs was lower than that of the known cDNAs, implying that these contigs included transcripts that were rarely represented in the full-length cDNA library. Within the non-redundant transcript set, we identified 892-2,833 responsive transcripts in roots and shoots, corresponding on average to 23.4% of the contigs not covered by cDNAs in TriFLDB under Pi starvation. The relative expression level of the wheat IPS1 (Induced by Phosphate Starvation 1) homologue, TaIPS1, was 341-fold higher in roots and 13-fold higher in shoots; this finding was further confirmed by qRT-PCR analysis. A comparative analysis of the wheat- and rice-responsive transcripts for orthologous genes under Pi-starvation revealed commonly upregulated transcripts, most of which appeared to be involved in a general response to Pi starvation, namely, an IPS1-mediated signalling cascade and its downstream functions such as Pi remobilisation, Pi uptake, and changes in Pi metabolism. CONCLUSIONS: Our transcriptome profiles demonstrated the impact of Pi starvation on global gene expression in wheat. This study revealed that enhancement of the Pi-mediated signalling cascade using IPS1 is a potent adaptation mechanism to Pi starvation that is conserved in both wheat and rice and validated the effectiveness of using short-read next-generation sequencing data for wheat transcriptome analysis in the absence of reference genome information.

Segregation distortion caused by weak hybrid necrosis in recombinant inbred lines of common wheat.

Segregation distortion of molecular markers is closely related to hybrid incompatibility in progeny from intraspecific crosses. Recent reports in higher plants have demonstrated that hybrid sterility results in segregation distortion at the causal gene regions in progeny of intraspecific crosses. Ne1 and Ne2 complementary loci are known to control hybrid necrosis in intraspecific crosses of common wheat cultivars. Here, we examine the effect of a weak necrosis allele Ne1(w) on the segregation ratio of molecular markers in recombinant inbred lines (RILs) of common wheat. Some RILs showed accelerated cell death in the leaves at the heading stage due to the epistatic interaction between two quantitative trait loci (QTL) on chromosomes 5B and 2B. Chromosomal localization of these QTL corresponding to Ne1(w) and Ne2 showed distorted segregation ratios of assigned markers having oppositely biased direction. Although the Ne1(w) and Ne2 interaction had no obvious effect on seed fertility, Ne1(w) reduced completion of grain development under the Ne2-homozygous background. This reduction might be one of causes that induces segregation distortion in the 5B and 2B chromosomal regions of RILs. The present study demonstrated that weak hybrid necrosis has limited phenotypic effects; it causes segregation distortion in progeny from intraspecific crosses.

A major quantitative trait locus for cold-responsive gene expression is linked to frost-resistance gene Fr-A2 in common wheat.

Low temperature induces expression of Cor (cold-responsive)/Lea (late embryogenesis-abundant) gene family members through C-repeat binding factor (CBF) transcription factors in common wheat. However, the relationship between the genetic loci controlling cold-responsive gene expression and freezing tolerance is unclear. In expression quantitative trait locus (eQTL) analysis, accumulated transcripts of Cor/Lea and CBF genes were quantified in recombinant inbred lines derived from a cross between two common wheat cultivars with different levels of freezing tolerance. Four eQTLs controlling five cold-responsive genes were found, and the major eQTL with the greatest effect was located on the long arm of chromosome 5A. At least the 1D and 5A eQTLs played important roles in development of freezing tolerance in common wheat. The chromosomal location of the 5A eQTL, controlling four cold-responsive genes, coincided with a region homoeologous to a frost-tolerance locus (Fr-A (m) 2) reported as a CBF cluster region in einkorn wheat. The 5A eQTL plays a significant role through Cor/Lea gene expression in cold acclimation of wheat. In addition, our results suggest that one or more CBF copies at the Fr-2 region positively regulate other copies, which might amplify the positive effects of the CBF cluster on downstream Cor/Lea gene activation.

bex-db: Bioinformatics workbench for comprehensive analysis of barley-expressed genes.

Barley (Hordeum vulgare) is one of the world's most important cereal crops. Although its large and complex genome has held back barley genomics for quite a while, the whole genome sequence was released in 2012 by the International Barley Genome Sequencing Consortium (IBSC). Moreover, more than 30,000 barley full-length cDNAs (FLcDNAs) are now available in the public domain. Here we present the Barley Gene Expression Database (bex-db: http://barleyflc.dna.affrc.go.jp/bexdb/index.html) as a repository of transcriptome data including the sequences and the expression profiles of barley genes resulting from microarray analysis. In addition to FLcDNA sequences, bex-db also contains partial sequences of more than 309,000 novel expressed sequence tags (ESTs). Users can browse the data via keyword, sequence homology and expression profile search options. A genome browser was also developed to display the chromosomal locations of barley FLcDNAs and wheat (Triticum aestivum) transcripts as well as Aegilops tauschii gene models on the IBSC genome sequence for future comparative analysis of orthologs among Triticeae species. The bex-db should provide a useful resource for further genomics studies and development of genome-based tools to enhance the progress of the genetic improvement of cereal crops.

Holographic polymer-dispersed liquid crystal memory for optically reconfigurable gate array using subwavelength grating mask.

Holographic polymer-dispersed liquid crystal (HPDLC) memory formed by a subwavelength grating (SWG) mask is presented for new optical information processing. The SWG structure in a photomask is formed on the SiO(2) plate using the anisotropic reactive ion etching technique. The configuration contexts for optically reconfigurable gate arrays (ORGAs) are stored in the HPDLC memory by polarization modulation property based on the form birefringence of the SWG plate. The configuration context pattern in the HPDLC memory is reconstructed to write it for the ORGAs under parallel programmability.

Overexpression of the pathogen-inducible wheat TaWRKY45 gene confers disease resistance to multiple fungi in transgenic wheat plants.

Recently we cloned and characterized the gene for the wheat transcription factor TaWRKY45 and showed that TaWRKY45 was upregulated in response to benzothiadiazole (BTH) and Fusarium head blight (FHB) and that its overexpression conferred enhanced resistance against F. graminearum. To characterize the functional role of TaWRKY45 in the disease resistance of wheat, in the present study we conducted expression analyses of TaWRKY45 with inoculations of powdery mildew and leaf rust and evaluated TaWRKY45-overexpressing wheat plants for resistance to these diseases. TaWRKY45 was upregulated in response to infections with Blumeria graminis, a causal fungus for powdery mildew, and Puccinia triticina, a causal fungus for leaf rust. Constitutive overexpression of the TaWRKY45 transgene conferred enhanced resistance against these two fungi on transgenic wheat plants grown under greenhouse conditions. However, the expression of two resistance-related genes, Pm3 and Lr34, was not induced by the inoculation with powdery mildew in TaWRKY45-overexpressing wheat plants. These results suggest that TaWRKY45 is involved in the defense responses for multiple fungal diseases in wheat but that resistance involving TaWRKY45 differs from at least Pm3 and/or Lr34-related resistance. Our present and previous studies indicate that TaWRKY45 may be potentially utilized to improve a wide range of disease resistance in wheat.

Multiple Wheat Genomes Reveal Novel Gli-2 Sublocus Location and Variation of Celiac Disease Epitopes in Duplicated α-Gliadin Genes.

The seed protein α-gliadin is a major component of wheat flour and causes gluten-related diseases. However, due to the complexity of this multigene family with a genome structure composed of dozens of copies derived from tandem and genome duplications, little was known about the variation between accessions, and thus little effort has been made to explicitly target α-gliadin for bread wheat breeding. Here, we analyzed genomic variation in α-gliadins across 11 recently published chromosome-scale assemblies of hexaploid wheat, with validation using long-read data. We unexpectedly found that the Gli-B2 locus is not a single contiguous locus but is composed of two subloci, suggesting the possibility of recombination between the two during breeding. We confirmed that the number of immunogenic epitopes among 11 accessions varied. The D subgenome of a European spelt line also contained epitopes, in agreement with its hybridization history. Evolutionary analysis identified amino acid sites under diversifying selection, suggesting their functional importance. The analysis opens the way for improved grain quality and safety through wheat breeding.

Identification of quantitative trait loci for ABA responsiveness at the seedling stage associated with ABA-regulated gene expression in common wheat.

Responsiveness to abscisic acid (ABA) during vegetative growth plays an important role in regulating adaptive responses to various environmental conditions, including activation of a number of ABA-responsive genes. However, the relationship between gene expression and responsiveness to ABA at the seedling stage has not been well studied in wheat. In the present study, quantitative trait locus (QTL) analysis for ABA responsiveness at the seedling stage was performed using recombinant inbred lines derived from a cross between common wheat cultivars showing different ABA responsiveness. Five QTLs were found to be significant, located on chromosomes 1B, 2A, 3A, 6D and 7B. The QTL with the greatest effect was located on chromosome 6D and explained 11.12% of the variance in ABA responsiveness. The other QTLs each accounted for approximately 5-8% of the phenotypic variation. Expression analyses of three ABA-responsive Cor/Lea genes, Wdhn13, Wrab15 and Wrab17, showed that allelic differences in QTLs on chromosomes 2A, 6D and 7B influenced expression of these genes in seedlings treated with ABA. The 3A QTL appeared to be involved in the regulatory system of Wdhn13 and Wrab15, but not Wrab17. The effects of the 2A and 6D QTLs on gene expression were relatively large. The combination of alleles at the QTLs resulted in an additive or synergistic effect on Cor/Lea expression. These results indicate that the QTLs influencing ABA responsiveness are associated with ABA-regulated gene expression and suggest that the QTL on chromosome 6D with the largest effect acts as a key regulator of ABA responses including seedling growth arrest and gene expression during the vegetative stage.