The drought-adapted shrub guayule (Parthenium argentatum) produces rubber, a natural product of major commercial importance, and two co-products with potential industrial use: terpene resin and the carbohydrate fructan. The rubber content of guayule plants subjected to water stress is higher compared to that of well-irrigated plants, a fact consistently reported in guayule field evaluations. To better understand how drought influences rubber biosynthesis at the molecular level, a comprehensive transcriptome database was built from drought-stressed guayule stem tissues using de novo RNA-seq and genome-guided assembly, followed by annotation and expression analysis. Despite having higher rubber content, most rubber biosynthesis related genes were down-regulated in drought-stressed guayule, compared to well-irrigated plants, suggesting post-transcriptional effects may regulate drought-induced rubber accumulation. On the other hand, terpene resin biosynthesis genes were unevenly affected by water stress, implying unique environmental influences over transcriptional control of different terpene compounds or classes. Finally, drought induced expression of fructan catabolism genes in guayule and significantly suppressed these fructan biosynthesis genes. It appears then, that in guayule cultivation, irrigation levels might be calibrated in such a regime to enable tunable accumulation of rubber, resin and fructan.

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and a valuable genetic resource for wheat improvement. Several reference-quality genome sequences have been reported for A. tauschii accession AL8/78. A new genome sequence assembly (Aet v6.0) built from long Pacific Biosciences HiFi reads and employing an optical genome map constructed with a new technology is reported here for this accession. The N50 contig length of 31.81 Mb greatly exceeded that of the previous AL8/78 genome sequence assembly (Aet v5.0). Of 1,254 super-scaffolds, 92, comprising 98% of the total super-scaffold length, were anchored on a high-resolution genetic map, and pseudomolecules were assembled. The number of gaps in the pseudomolecules was reduced from 52,910 in Aet v5.0 to 351 in Aet v6.0. Gene models were transferred from the Aet v5.0 assembly into the Aet v6.0 assembly. A total of 40,447 putative orthologous gene pairs were identified between the Aet v6.0 and Chinese Spring wheat IWGSC RefSer v2.1 D-subgenome pseudomolecules. Orthologous gene pairs were used to compare the structure of the A. tauschii and wheat D-subgenome pseudomolecules. A total of 223 structural differences were identified. They included 44 large differences in sequence orientation and 25 differences in sequence location. A technique for discriminating between assembly errors and real structural variation between closely related genomes is suggested.

Aegilops caudata L. is a rich source of resistance genes to major wheat pathogens and pests, yet its complex genome structure and lack of a reference sequence have hindered genetic analysis, gene discovery, and introgression. Here, we report a high-quality genome sequence assembly of Ae. caudata accession S740-69, integrated with optical mapping to precisely delineate introgressed segments from structurally altered genomic regions in wheat-Ae. caudata translocation lines carrying a novel stem rust resistance gene. Resequencing and phenotyping of 95 diverse accessions, coupled with k-mer-based association mapping, enabled population-level identification of several novel loci for resistance to wheat pests and diseases, such as greenbug and stem rust. Virus-induced gene silencing confirmed Aecau6C01G127270 (SrAect1) as the first functionally characterized stem rust resistance gene from Ae. caudata. This new genomic framework provides a foundation for systematic mining of resistance genes in this species and their efficient introgression into wheat. ### Competing Interest Statement The authors have declared no competing interest. United States Department of Agriculture, Agricultural Research Service (USDA-ARS), 2030-21430-015-000D

This work reports the draft genome of Agrobacterium fabrum strain 1D1416. The assembled genome is composed of a 2,837,379-bp circular chromosome, a 2,043,296-bp linear chromosome, a 519,735-bp AT1 plasmid, a 188,396-bp AT2 plasmid, and a 196,706-bp Ti virulence plasmid. The nondisarmed strain produces gall-like structures in citrus tissue.

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 R-390 and T-392 residues of TGB1(ND) (ND18 strain) and the G(196) and K-197 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.

This work reports the draft genome sequence of Agrobacterium fabrum strain 1D1104. The assembled genome is composed of a 2,774,783-bp circular chromosome, a 2,110,112-bp linear chromosome, an AT plasmid of 133,577 bp, and four unassembled contigs of 5,389,544 bp, 42,391 bp, 41,768 bp, and 35,476 bp.

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.

Powdery mildew, caused by Blumeria graminis f. sp. tritici ( Bgt ), is one of the most destructive diseases that pose a great threat to wheat production. Wheat landraces represent a rich source of powdery mildew resistance. Here, we report the map-based cloning of powdery mildew resistance gene Pm24 from Chinese wheat landrace Hulutou. It encodes a tandem kinase protein (TKP) with putative kinase-pseudokinase domains, designated WHEAT TANDEM KINASE 3 (WTK3). The resistance function of Pm24 was validated by transgenic assay, independent mutants, and allelic association analyses. Haplotype analysis revealed that a rare 6-bp natural deletion of lysine-glycine codons, endemic to wheat landraces of Shaanxi Province, China, in the kinase I domain (Kin I) of WTK3 is critical for the resistance function. Transgenic assay of WTK3 chimeric variants revealed that only the specific two amino acid deletion, rather than any of the single or more amino acid deletions, in the Kin I of WTK3 is responsible for gaining the resistance function of WTK3 against the Bgt fungus.

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.

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically-contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants, and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

This work reports the draft genome sequence of Agrobacterium tumefaciens strain 1D1526. The assembled genome is composed of a 2,881,823-bp circular chromosome, a 2,235,711-bp linear chromosome, and a 44,582-bp unassembled contig.

Sharp eyespot, mainly caused by the soil-borne fungus Rhizoctonia cerealis, affects wheat（Triticum aestivum L.） production worldwide. In this study, we isolated Ta CML36 gene encoding a wheat calmodulin-like protein, and studied its defense role in protection against R. cerealis. Transcription of Ta CML36 was significantly elevated by both R. cerealis infection and exogenous ethylene treatment. Transcription was higher in resistant wheat lines than in susceptible ones. There were copies of Ta CML36 on chromosomes 5 A, 5 B, and 5 D. The Ta CML36 protein is composed of 183 amino acids and contains two calcium-binding EFhand domains. Subcellular localization assays in wheat indicated that Ta CML36 localizes in both the cytoplasm and nucleus. Virus-induced gene silencing and disease assessment indicated that compared to the controls, Ta CML36-silenced wheat plants displayed significantly reduced resistance to R. cerealis and had greater fungal biomass, suggesting that knockdown of Ta CML36 impaired host resistance. Knockdown of Ta CML36 also significantly repressed expression of pathogenesis-related genes such as Chitinase 1,PDF35, and PR17 C, the ethylene response factor-encoding gene Ta PIE1, and ethylene biosynthesis gene ACO2. Collectively, our results suggest that Ta CML36 positively participates in the innate immune response to R. cerealis infection by modulating expression of defense-associated genes possibly in the ethylene signaling pathway.

This work reports the draft genome sequence of Agrobacterium fabrum strain 1D159 (also known as ATCC strain 27912). The assembled genome is composed of a 2,861,352-bp circular chromosome, a 2,058,040-bp linear chromosome, a 519,735-bp AT plasmid, and the 223,394-bp Ti virulence plasmid. The wild nondisarmed strain produces small gall-like structures in citrus.

α-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.

This work reports the draft genome of Serratia sp. 1D1416. The assembled genome contains a 5,552,016-bp circular chromosome. The strain was discovered in a mixed culture from a gall isolated from Euonymus japonicas.

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.

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.8Mb 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.

Background Guayule ( Parthenium argentatum A. Gray) is a rubber-producing desert shrub native to Mexico and the United States. Guayule represents an alternative to Hevea brasiliensis as a source for commercial natural rubber. The efficient application of modern molecular/genetic tools to guayule improvement requires characterization of its genome. Results The 1.6 Gb guayule genome was sequenced, assembled and annotated. The final 1.5 Gb assembly, while fragmented (N 50 = 22 kb), maps > 95% of the shotgun reads and is essentially complete. Approximately 40,000 transcribed, protein encoding genes were annotated on the assembly. Further characterization of this genome revealed 15 families of small, microsatellite-associated, transposable elements (TEs) with unexpected chromosomal distribution profiles. These SaTar ( Sa tellite Tar geted) elements, which are non-autonomous Mu- like elements (MULEs), were frequently observed in multimeric linear arrays of unrelated individual elements within which no individual element is interrupted by another. This uniformly non-nested TE multimer architecture has not been previously described in either eukaryotic or prokaryotic genomes. Five families of similarly distributed non-autonomous MULEs (microsatellite associated, modularly assembled) were characterized in the rice genome. Families of TEs with similar structures and distribution profiles were identified in sorghum and citrus. Conclusion The sequencing and assembly of the guayule genome provides a foundation for application of current crop improvement technologies to this plant. In addition, characterization of this genome revealed SaTar elements with distribution profiles unique among TEs. Satar targeting appears based on an alternative MULE recombination mechanism with the potential to impact gene evolution.

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.

Gliadins are a major component of wheat seed proteins. However, the complex homoeologous Gli-2 loci (Gli-A2, -B2 and -D2) that encode the α-gliadins in commercial wheat are still poorly understood. Here we analyzed the Gli-D2 locus of Xiaoyan 81 (Xy81), a winter wheat cultivar. A total of 421.091 kb of the Gli-D2 sequence was assembled from sequencing multiple bacterial artificial clones, and 10 α-gliadin genes were annotated. Comparative genomic analysis showed that Xy81 carried only eight of the α-gliadin genes of the D genome donor Aegilops tauschii, with two of them each experiencing a tandem duplication. A mutant line lacking Gli-D2 (DLGliD2) consistently exhibited better breadmaking quality and dough functionalities than its progenitor Xy81, but without penalties in other agronomic traits. It also had an elevated lysine content in the grains. Transcriptome analysis verified the lack of Gli-D2 α-gliadin gene expression in DLGliD2. Furthermore, the transcript and protein levels of protein disulfide isomerase were both upregulated in DLGliD2 grains. Consistent with this finding, DLGliD2 had increased disulfide content in the flour. Our work sheds light on the structure and function of Gli-D2 in commercial wheat, and suggests that the removal of Gli-D2 and the gliadins specified by it is likely to be useful for simultaneously enhancing the end-use and health-related traits of common wheat. Because gliadins and homologous proteins are widely present in grass species, the strategy and information reported here may be broadly useful for improving the quality traits of diverse cereal crops.