An online database for einkorn wheat to aid in gene discovery and functional genomics studies.

Diploid A-genome wheat (einkorn wheat) presents a nutrition-rich option as an ancient grain crop and a resource for the improvement of bread wheat against abiotic and biotic stresses. Realizing the importance of this wheat species, reference-level assemblies of two einkorn wheat accessions were generated (wild and domesticated). This work reports an einkorn genome database that provides an interface to the cereals research community to perform comparative genomics, applied genetics and breeding research. It features queries for annotated genes, the use of a recent genome browser release, and the ability to search for sequence alignments using a modern BLAST interface. Other features include a comparison of reference einkorn assemblies with other wheat cultivars through genomic synteny visualization and an alignment visualization tool for BLAST results. Altogether, this resource will help wheat research and breeding. Database URL  https://wheat.pw.usda.gov/GG3/pangenome.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool.

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5Uu and 5Mg) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5Mg, 5Uu, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

Toward reducing the immunogenic potential of wheat flour: identification and characterization of wheat lines missing omega-5 gliadins encoded by the 1D chromosome.

KEY MESSAGE: Eleven wheat lines that are missing genes for the 1D-encoded omega-5 gliadins will facilitate breeding efforts to reduce the immunogenic potential of wheat flour for patients susceptible to wheat allergy. Efforts to reduce the levels of allergens in wheat flour that cause wheat-dependent exercise-induced anaphylaxis are complicated by the presence of genes encoding omega-5 gliadins on both chromosomes 1B and 1D of hexaploid wheat. In this study, we screened 665 wheat germplasm samples using gene specific DNA markers for omega-5 gliadins encoded by the genes on 1D chromosome that were obtained from the reference wheat Chinese Spring. Eleven wheat lines missing the PCR product corresponding to 1D omega-5 gliadin gene sequences were identified. Two of the lines contained the 1BL·1RS translocation. Relative quantification of gene copy numbers by qPCR revealed that copy numbers of 1D omega-5 gliadins in the other nine lines were comparable to those in 1D null lines of Chinese Spring, while copy numbers of 1B omega-5 gliadins were like those of Chinese Spring. 2-D immunoblot analysis of total flour proteins from the selected lines using a specific monoclonal antibody against the N-terminal sequence of omega-5 gliadin showed no reactivity in regions of the blots containing previously identified 1D omega-5 gliadins. Interestingly, RP-UPLC analysis of the gliadin fractions of the selected lines indicated that the expression of omega-1,2 gliadins was also significantly reduced in seven of the lines, implying that 1D omega-5 gliadin and 1D omega-1,2 gliadin genes are tightly linked on the Gli-D1 loci of chromosome 1D. Wheat lines missing the omega-5 gliadins encoded by the genes on 1D chromosome should be useful in future breeding efforts to reduce the immunogenic potential of wheat flour.

The CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein.

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R390 and T392 residues of TGB1ND (ND18 strain) and the G196 and K197 residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

CRISPR-Cas9 Gene Editing of the Sal1 Gene Family in Wheat.

The highly conserved Sal1 encodes a bifunctional enzyme with inositol polyphosphate-1-phosphatase and 3' (2'), 5'-bisphosphate nucleotidase activity and has been shown to alter abiotic stress tolerance in plants when disrupted. Precise gene editing techniques were used to generate Sal1 mutants in hexaploid bread wheat. The CRISPR (Clustered Regulatory Interspaced Short Palindromic Repeats) Cas9 system with three guide RNAs (gRNAs) was used to inactivate six Sal1 homologous genes within the Bobwhite wheat genome. The resulting mutant wheat plants with all their Sal1 genes disabled had slimmer stems, had a modest reduction in biomass and senesced more slowly in water limiting conditions, but did not exhibit improved yield under drought conditions. Our results show that multiplexed gRNAs enabled effective targeted gene editing of the Sal1 gene family in hexaploid wheat. These Sal1 mutant wheat plants will be a resource for further research studying the function of this gene family in wheat.

A Rapid, Reliable RP-UPLC Method for Large-Scale Analysis of Wheat HMW-GS Alleles.

High-molecular-weight glutenin subunits (HMW-GS) account for only 10% of total wheat storage proteins, but play an important role in the processing quality of wheat flour. Therefore, identifying HMW-GS alleles associated with good end-use quality provides important information for wheat breeders. To rapidly, accurately and reproducibly identify HMW-GS, we established an optimized reversed-phase ultra-performance liquid chromatography (RP-UPLC) method. Separation parameters were optimized using an ACQUITY UPLC Protein BEH C4 column and stepwise ACN gradient, and the separation patterns and retention times (RTs) of 22 subunits were comparatively analyzed in 16 standard wheat cultivars. All HMW-GS proteins were well separated within about 5.5 min, and all analyses were complete within 12 min. We distinguished the 16 subunits based on RT, although three subunits in 1Bx (1Bx7/1Bx7OE and 1Bx17) and three subunits in 1By (1By8*, 1By9 and 1By15) had overlapping RTs; these were differentiated by SDS-PAGE. To distinguish 1Bx7 and 1Bx7OE, which differ in protein abundance, RP-UPLC was combined with PCR analysis of DNA junction markers. The optimized method was successfully applied to determine HMW-GS alleles in a large collection of bread wheat germplasm (1787 lines). This protocol is an appropriate option for selecting lines harboring favorable HMW-GS alleles in wheat breeding.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation.

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Proteomic Determination of Low-Molecular-Weight Glutenin Subunit Composition in Aroona Near-Isogenic Lines and Standard Wheat Cultivars.

The low-molecular weight glutenin subunit (LMW-GS) composition of wheat (Triticum aestivum) flour has important effects on end-use quality. However, assessing the contributions of each LMW-GS to flour quality remains challenging because of the complex LMW-GS composition and allelic variation among wheat cultivars. Therefore, accurate and reliable determination of LMW-GS alleles in germplasm remains an important challenge for wheat breeding. In this study, we used an optimized reversed-phase HPLC method and proteomics approach comprising 2-D gels coupled with liquid chromatography-tandem mass spectrometry (MS/MS) to discriminate individual LMW-GSs corresponding to alleles encoded by the Glu-A3, Glu-B3, and Glu-D3 loci in the 'Aroona' cultivar and 12 'Aroona' near-isogenic lines (ARILs), which contain unique LMW-GS alleles in the same genetic background. The LMW-GS separation patterns for 'Aroona' and ARILs on chromatograms and 2-D gels were consistent with those from a set of 10 standard wheat cultivars for Glu-3. Furthermore, 12 previously uncharacterized spots in 'Aroona' and ARILs were excised from 2-D gels, digested with chymotrypsin, and subjected to MS/MS. We identified their gene haplotypes and created a 2-D gel map of LMW-GS alleles in the germplasm for breeding and screening for desirable LMW-GS alleles for wheat quality improvement.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly.

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

Characterization of Pm68, a new powdery mildew resistance gene on chromosome 2BS of Greek durum wheat TRI 1796.

KEY MESSAGE: New powdery mildew resistance gene Pm68 was found in the terminal region of chromosome 2BS of Greek durum wheat TRI 1796. The co-segregated molecular markers could be used for MAS. Durum wheat (Triticum turgidum L. var. durum Desf.) is not only an important cereal crop for pasta making, but also a genetic resource for common wheat improvement. In the present study, a Greek durum wheat TRI 1796 was found to confer high resistance to all 22 tested isolates of Blumeria graminis f. sp. tritici (Bgt). Inheritance study on the F1 plants and the F2 population derived from the cross TRI 1796/PI 584832 revealed that the resistance in TRI 1796 was controlled by a single dominant gene, herein designated Pm68. Using the bulked segregant RNA-Seq (BSR-Seq) analysis combined with molecular analysis, Pm68 was mapped to the terminal part of the short arm of chromosome 2B and flanked by markers Xdw04 and Xdw12/Xdw13 with genetic distances of 0.22 cM each. According to the reference genome of durum wheat cv. Svevo, the corresponding physical region spanned the Pm68 locus was about 1.78-Mb, in which a number of disease resistance-related genes were annotated. This study reports the new powdery mildew resistance gene Pm68 that would be a valuable resource for improvement of both common wheat and durum wheat. The co-segregated markers (Xdw05-Xdw11) developed here would be useful tools for marker-assisted selection (MAS) in breeding.

Comparison of MALDI-TOF-MS and RP-HPLC as Rapid Screening Methods for Wheat Lines With Altered Gliadin Compositions.

The wheat gliadins are a complex group of flour proteins that can trigger celiac disease and serious food allergies. As a result, mutation breeding and biotechnology approaches are being used to develop new wheat lines with reduced immunogenic potential. Key to these efforts is the development of rapid, high-throughput methods that can be used as a first step in selecting lines with altered gliadin contents. In this paper, we optimized matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and reversed-phase high-performance liquid chromatography (RP-HPLC) methods for the separation of gliadins from Triticum aestivum cv. Chinese Spring (CS). We evaluated the quality of the resulting profiles using the complete set of gliadin gene sequences recently obtained from this cultivar as well as a set of aneuploid lines in CS. The gliadins were resolved into 13 peaks by MALDI-TOF-MS. α- or γ-gliadins that contain abundant celiac disease epitopes and are likely targets for efforts to reduce the immunogenicity of flour were found in several peaks. However, other peaks contained multiple α- and γ-gliadins, including one peak with as many as 12 different gliadins. In comparison, separation of proteins by RP-HPLC yielded 28 gliadin peaks, including 13 peaks containing α-gliadins and eight peaks containing γ-gliadins. While the separation of α- and γ-gliadins gliadins achieved by RP-HPLC was better than that achieved by MALDI-TOF-MS, it was not possible to link peaks with individual protein sequences. Both MALDI-TOF-MS and RP-HPLC provided adequate separation of ω-gliadins. While MALDI-TOF-MS is faster and could prove useful in studies that target specific gliadins, RP-HPLC is an effective method that can be applied more broadly to detect changes in gliadin composition.

A high-quality reference genome sequence of Salvia miltiorrhiza provides insights into tanshinone synthesis in its red rhizomes.

Salvia miltiorrhiza Bunge, also known as red sage or Danshen, is an important traditional Chinese medicine (TCM) that has been used for thousands of years to treat cardiovascular and other diseases. It is also considered an important model TCM plant. Here, a high-quality reference genome of S. miltiorrhiza was generated by combining PacBio long-read sequencing and chromatin interaction mapping (Hi-C) technologies, resulting in the chromosome-scale assembly of a 594.75-Mb genome sequence with a contig N50 of 2.70 Mb. This assembly shows the highest level of continuity for a Danshen genome generated thus far. The S. miltiorrhiza genome contained 32,483 protein-coding genes, with a repetitive DNA content of approximately 64.84%. The high percentage of young LTRs suggests that multiple TE transposition bursts occurred recently in S. miltiorrhiza. Genes unique to secondary metabolism pathways were expanded in the S. miltiorrhiza genome. A new CYP450 gene cluster was identified in the phloem of red roots where active components were synthesized. This reference genome sequence will facilitate future studies aimed at the elucidation of the secondary metabolism synthesis pathway and the genetic improvement of S. miltiorrhiza.

Exploiting the reference genome sequence of hexaploid wheat: a proteomic study of flour proteins from the cultivar Chinese Spring.

Although the economic value of wheat flour is determined by the complement of gluten proteins, these proteins have been challenging to study because of the complexity of the major protein groups and the tremendous sequence diversity among wheat cultivars. The completion of a high-quality wheat genome sequence from the reference wheat Chinese Spring recently facilitated the assembly and annotation of a complete set of gluten protein genes from a single cultivar, making it possible to link individual proteins in the flour to specific gene sequences. In a proteomic analysis of total wheat flour protein from Chinese Spring using quantitative two-dimensional gel electrophoresis combined with tandem mass spectrometry, gliadins or low-molecular-weight glutenin subunits were identified as the predominant proteins in 72 protein spots. Individual spots were associated with 40 of 56 Chinese Spring gene sequences, including 16 of 26 alpha gliadins, 10 of 11 gamma gliadins, six of seven omega gliadins, one of two delta gliadins, and nine of ten LMW-GS. Most genes that were not associated with protein spots were either expressed at low levels in endosperm or encoded proteins with high similarity to other proteins. A wide range of protein accumulation levels were observed and discrepancies between transcript levels and protein levels were noted. This work together with similar studies using other commercial cultivars should provide new insight into the molecular basis of wheat flour quality and allergenic potential.

GrainGenes: centralized small grain resources and digital platform for geneticists and breeders.

GrainGenes (https://wheat.pw.usda.gov or https://graingenes.org) is an international centralized repository for curated, peer-reviewed datasets useful to researchers working on wheat, barley, rye and oat. GrainGenes manages genomic, genetic, germplasm and phenotypic datasets through a dynamically generated web interface for facilitated data discovery. Since 1992, GrainGenes has served geneticists and breeders in both the public and private sectors on six continents. Recently, several new datasets were curated into the database along with new tools for analysis. The GrainGenes homepage was enhanced by making it more visually intuitive and by adding links to commonly used pages. Several genome assemblies and genomic tracks are displayed through the genome browsers at GrainGenes, including the Triticum aestivum (bread wheat) cv. 'Chinese Spring' IWGSC RefSeq v1.0 genome assembly, the Aegilops tauschii (D genome progenitor) Aet v4.0 genome assembly, the Triticum turgidum ssp. dicoccoides (wild emmer wheat) cv. 'Zavitan' WEWSeq v.1.0 genome assembly, a T. aestivum (bread wheat) pangenome, the Hordeum vulgare (barley) cv. 'Morex' IBSC genome assembly, the Secale cereale (rye) select 'Lo7' assembly, a partial hexaploid Avena sativa (oat) assembly and the Triticum durum cv. 'Svevo' (durum wheat) RefSeq Release 1.0 assembly. New genetic maps and markers were added and can be displayed through CMAP. Quantitative trait loci, genetic maps and genes from the Wheat Gene Catalogue are indexed and linked through the Wheat Information System (WheatIS) portal. Training videos were created to help users query and reach the data they need. GSP (Genome Specific Primers) and PIECE2 (Plant Intron Exon Comparison and Evolution) tools were implemented and are available to use. As more small grains reference sequences become available, GrainGenes will play an increasingly vital role in helping researchers improve crops.

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.

OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species.

OrthoVenn is a powerful web platform for the comparison and analysis of whole-genome orthologous clusters. Here we present an updated version, OrthoVenn2, which provides new features that facilitate the comparative analysis of orthologous clusters among up to 12 species. Additionally, this update offers improvements to data visualization and interpretation, including an occurrence pattern table for interrogating the overlap of each orthologous group for the queried species. Within the occurrence table, the functional annotations and summaries of the disjunctions and intersections of clusters between the chosen species can be displayed through an interactive Venn diagram. To facilitate a broader range of comparisons, a larger number of species, including vertebrates, metazoa, protists, fungi, plants and bacteria, have been added in OrthoVenn2. Finally, a stand-alone version is available to perform large dataset comparisons and to visualize results locally without limitation of species number. In summary, OrthoVenn2 is an efficient and user-friendly web server freely accessible at https://orthovenn2.bioinfotoolkits.net.

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.

Reassessment of the evolution of wheat chromosomes 4A, 5A, and 7B.

KEY MESSAGE: Comparison of genome sequences of wild emmer wheat and Aegilops tauschii suggests a novel scenario of the evolution of rearranged wheat chromosomes 4A, 5A, and 7B. Past research suggested that wheat chromosome 4A was subjected to a reciprocal translocation T(4AL;5AL)1 that occurred in the diploid progenitor of the wheat A subgenome and to three major rearrangements that occurred in polyploid wheat: pericentric inversion Inv(4AS;4AL)1, paracentric inversion Inv(4AL;4AL)1, and reciprocal translocation T(4AL;7BS)1. Gene collinearity along the pseudomolecules of tetraploid wild emmer wheat (Triticum turgidum ssp. dicoccoides, subgenomes AABB) and diploid Aegilops tauschii (genomes DD) was employed to confirm these rearrangements and to analyze the breakpoints. The exchange of distal regions of chromosome arms 4AS and 4AL due to pericentric inversion Inv(4AS;4AL)1 was detected, and breakpoints were validated with an optical Bionano genome map. Both breakpoints contained satellite DNA. The breakpoints of reciprocal translocation T(4AL;7BS)1 were also found. However, the breakpoints that generated paracentric inversion Inv(4AL;4AL)1 appeared to be collocated with the 4AL breakpoints that had produced Inv(4AS;4AL)1 and T(4AL;7BS)1. Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 either originated sequentially, and Inv(4AL;4AL)1 was produced by recurrent chromosome breaks at the same breakpoints that generated Inv(4AS;4AL)1 and T(4AL;7BS)1, or Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 originated simultaneously. We prefer the latter hypothesis since it makes fewer assumptions about the sequence of events that produced these chromosome rearrangements.

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.

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size.

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.