Genome sequences of three Aegilops species of the section Sitopsis reveal phylogenetic relationships and provide resources for wheat improvement.

Aegilops is a close relative of wheat (Triticum spp.), and Aegilops species in the section Sitopsis represent a rich reservoir of genetic diversity for the improvement of wheat. To understand their diversity and advance their utilization, we produced whole-genome assemblies of Aegilops longissima and Aegilops speltoides. Whole-genome comparative analysis, along with the recently sequenced Aegilops sharonensis genome, showed that the Ae. longissima and Ae. sharonensis genomes are highly similar and are most closely related to the wheat D subgenome. By contrast, the Ae. speltoides genome is more closely related to the B subgenome. Haplotype block analysis supported the idea that Ae. speltoides genome is closest to the wheat B subgenome, and highlighted variable and similar genomic regions between the three Aegilops species and wheat. Genome-wide analysis of nucleotide-binding leucine-rich repeat (NLR) genes revealed species-specific and lineage-specific NLR genes and variants, demonstrating the potential of Aegilops genomes for wheat improvement.

New insights into the dispersion history and adaptive evolution of taxon Aegilops tauschii in China.

Aegilops tauschii, the wild progenitor of wheat D-genome and a valuable germplasm for wheat improvement, has a wide natural distribution from eastern Turkey to China. However, the phylogenetic relationship and dispersion history of Ae. tauschii in China has not been scientifically clarified. In this study, we genotyped 208 accessions (with 104 in China) using ddRAD sequencing and 55K SNP array, and classified the population into six sublineages. Three possible spreading routes or events were identified, resulting in specific distribution patterns, with four sublineages found in Xinjiang, one in Qinghai, two in Shaanxi and one in Henan. We also established the correlation of SNP-based, karyotype-based and spike-morphology-based techniques to demonstrate the internal classification of Ae. tauschii, and developed consensus dataset with 1245 putative accessions by merging data previously published. Our analysis suggested that eight inter-lineage accessions could be assigned to the putative Lineage 3 and these accessions would help to conserve the genetic diversity of the species. By developing the consensus phylogenetic relationships of Ae. tauschii, our work validated the hypothesis on the dispersal history of Ae. tauschii in China, and contributed to the efficient and comprehensive germplasm-mining of the species.

Wild emmer introgression alters root-to-shoot growth dynamics in durum wheat in response to water stress.

Water deficit during the early vegetative growth stages of wheat (Triticum) can limit shoot growth and ultimately impact grain productivity. Introducing diversity in wheat cultivars to enhance the range of phenotypic responses to water limitations during vegetative growth can provide potential avenues for mitigating subsequent yield losses. We tested this hypothesis in an elite durum wheat background by introducing a series of introgressions from a wild emmer (Triticum turgidum ssp. dicoccoides) wheat. Wild emmer populations harbor rich phenotypic diversity for drought-adaptive traits. To determine the effect of these introgressions on vegetative growth under water-limited conditions, we used image-based phenotyping to catalog divergent growth responses to water stress ranging from high plasticity to high stability. One of the introgression lines exhibited a significant shift in root-to-shoot ratio in response to water stress. We characterized this shift by combining genetic analysis and root transcriptome profiling to identify candidate genes (including a root-specific kinase) that may be linked to the root-to-shoot carbon reallocation under water stress. Our results highlight the potential of introducing functional diversity into elite durum wheat for enhancing the range of water stress adaptation.

The Independent Domestication of Timopheev's Wheat: Insights from Haplotype Analysis of the Brittle rachis 1 (BTR1-A) Gene.

Triticum turgidum and T. timopheevii are two tetraploid wheat species sharing T. urartu as a common ancestor, and domesticated accessions from both of these allopolyploids exhibit nonbrittle rachis (i.e., nonshattering spikes). We previously described the loss-of-function mutations in the Brittle Rachis 1 genes BTR1-A and BTR1-B in the A and B subgenomes, respectively, that are responsible for this most visible domestication trait in T. turgidum. Resequencing of a large panel of wild and domesticated T. turgidum accessions subsequently led to the identification of the two progenitor haplotypes of the btr1-A and btr1-B domesticated alleles. Here, we extended the haplotype analysis to other T. turgidum subspecies and to the BTR1 homologues in the related T. timopheevii species. Our results showed that all the domesticated wheat subspecies within T. turgidum share common BTR1-A and BTR1-B haplotypes, confirming their common origin. In T. timopheevii, however, we identified a novel loss-of-function btr1-A allele underlying a partially brittle spike phenotype. This novel recessive allele appeared fixed within the pool of domesticated Timopheev's wheat but was also carried by one wild timopheevii accession exhibiting partial brittleness. The promoter region for BTR1-B could not be amplified in any T. timopheevii accessions with any T. turgidum primer combination, exemplifying the gene-level distance between the two species. Altogether, our results support the concept of independent domestication processes for the two polyploid, wheat-related species.

Genome-Wide Mapping of Loci for Adult-Plant Resistance to Stripe Rust in Durum Wheat Svevo Using the 90K SNP Array.

Stripe rust is a foliar disease in wheat caused by Puccinia striiformis f. tritici. The best way to protect wheat from this disease is by growing resistant cultivars. Tetraploid wheat can serve as a good source of valuable genetic diversity for various traits. Here, we report the mapping of nine stripe rust resistance quantitative trait loci (QTL) effective against P. striiformis f. tritici in China and Israel. We used recombinant inbred lines (RILs) developed from a cross between the durum wheat cultivar Svevo and Triticum dicoccoides accession Zavitan. By genotyping the RIL population of 137 lines using the wheat 90K single-nucleotide polymorphism array, we mapped an adult-plant resistance locus QYrsv.swust-1BL.1, the most effective QTL, within a 0.75-centimorgan region in T. turgidum subsp. durum 'Svevo' on chromosome arm 1BL, corresponding to the region of 670.7 to 671.5 Mb on the Chinese Spring chromosome arm 1BL. Of the other eight minor-effect stripe rust QTL, seven were from Svevo and mapped on chromosomes 1A, 1B, 2B, 3A, 4A, and 5A, and one was from Zavitan and mapped on chromosome 2A. Several QTL with epistatic effects were identified as well. The markers linked to the resistance QTL can be useful in marker-assisted selection for incorporation of these resistance QTL into both durum and common wheat cultivars.

The Brittle Rachis Trait in Species Belonging to the Triticeae and Its Controlling Genes Btr1 and Btr2.

In many non-cultivated angiosperm species, seed dispersal is facilitated by the shattering of the seed head at maturity; in the Triticeae tribe, to which several of the world's most important cereals belong, shattering takes the form of a disarticulation of the rachis. The products of the genes Btr1 and Btr2 are both required for disarticulation to occur above the rachis nodes within the genera Hordeum (barley) and Triticum/Aegilops (wheat). Here, it has been shown that both Btr1 and Btr2 are specific to the Triticeae tribe, although likely paralogs (Btr1-like and Btr2-like) are carried by the family Poaceae including Triticeae. Aegilops tauschii (the donor of the bread wheat D genome) lacks a copy of Btr1 and disarticulation in this species occurs below, rather than above the rachis node; thus, the product of Btr1 appears to be required for disarticulation to occur above the rachis node.

Correction to: Recombination between homoeologous chromosomes induced in durum wheat by the Aegilops speltoides Su1-Ph1 suppressor.

Unfortunately, the 9th author name was incorrectly published in the original publication. The complete correct name is given below.

Recombination between homoeologous chromosomes induced in durum wheat by the Aegilops speltoides Su1-Ph1 suppressor.

KEY MESSAGE: Su1-Ph1, which we previously introgressed into wheat from Aegilops speltoides, is a potent suppressor of Ph1 and a valuable tool for gene introgression in tetraploid wheat. We previously introgressed Su1-Ph1, a suppressor of the wheat Ph1 gene, from Aegilops speltoides into durum wheat cv Langdon (LDN). Here, we evaluated the utility of the introgressed suppressor for inducing introgression of alien germplasm into durum wheat. We built LDN plants heterozygous for Su1-Ph1 that simultaneously contained a single LDN chromosome 5B and a single Ae. searsii chromosome 5Sse, which targeted them for recombination. We genotyped 28 BC1F1 and 84 F2 progeny with the wheat 90-K Illumina single-nucleotide polymorphism assay and detected extensive recombination between the two chromosomes, which we confirmed by non-denaturing fluorescence in situ hybridization (ND-FISH). We constructed BC1F1 and F2 genetic maps that were 65.31 and 63.71 cM long, respectively. Recombination rates between the 5B and 5Sse chromosomes were double the expected rate computed from their meiotic pairing, which we attributed to selection against aneuploid gametes. Recombination rate between 5B and 5Sse was depressed compared to that between 5B chromosomes in the proximal region of the long arm. We integrated ND-FISH signals into the genetic map and constructed a physical map, which we used to map a 172,188,453-bp Ph1 region. Despite the location of the region in a low-recombination region of the 5B chromosome, we detected three crossovers in it. Our data show that Su1-Ph1 is a valuable tool for gene introgression and gene mapping based on recombination between homoeologous chromosomes in wheat.

Wheat domestication in light of haplotype analyses of the Brittle rachis 1 genes (BTR1-A and BTR1-B).

Wheat domestication was a milestone in the rise of agrarian societies in the Fertile Crescent. As opposed to the freely dispersing seeds of its tetraploid progenitor wild emmer, the hallmark trait of domesticated wheat is intact, harvestable spikes. During domestication, wheat acquired recessive loss-of-function mutations in the Brittle Rachis 1 genes, both in the A genome (BTR1-A) and B genome (BTR1-B). In this study, we probe the geographical provenances of these mutations via haplotype analyses of a collection of wild and domesticated accessions. Our results show that the precursor of the domesticated haplotype of BTR1-A was detected in 32% of the wild accessions gathered throughout the Levant, from central Israel to central Turkey. In contrast, the precursor of the domesticated haplotype of BTR1-B, which carries a distinct 11 bp deletion in the promoter region, was found in only 10% of the tested wild accessions, all from the Southern Levant. Moreover, we identified of a single wild emmer accession in Southern Levant that carries the progenitor haplotypes for both BTR1-A and BTR1-B genes. These observations suggest that at least part of the emmer domestication process occurred in Southern Levant, contrary to the widely held view that the northern part of the Fertile Crescent was the center of wheat domestication.

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.

GNI-A1 mediates trade-off between grain number and grain weight in tetraploid wheat.

KEY MESSAGE: Wild emmer allele of GNI-A1 ease competition among developing grains through the suppression of floret fertility and increase grain weight in tetraploid wheat. Grain yield is a highly polygenic trait determined by the number of grains per unit area, as well as by grain weight. In wheat, grain number and grain weight are usually negatively correlated. Yet, the genetic basis underlying trade-off between the two is mostly unknown. Here, we fine-mapped a grain weight QTL using wild emmer introgressions in a durum wheat background and showed that grain weight is associated with the GNI-A1 gene, a regulator of floret fertility. In-depth characterization of grain number and grain weight indicated that suppression of distal florets by the wild emmer GNI-A1 allele increases weight of proximal grains in basal and central spikelets due to alteration in assimilate distribution. Re-sequencing of GNI-A1 in tetraploid wheat demonstrated the rich allelic repertoire of the wild emmer gene pool, including a rare allele which was present in two gene copies and contained a nonsynonymous mutation in the C-terminus of the protein. Using an F2 population generated from a cross between wild emmer accessions Zavitan, which carries the rare allele, and TTD140, we demonstrated that this unique polymorphism is associated with grain weight, independent of grain number. Moreover, we showed, for the first time, that GNI-A1 proteins are transcriptional activators and that selection targeted compromised activity of the protein. Our findings expand the knowledge of the genetic basis underlying trade-off between key yield components and may contribute to breeding efforts for enhanced grain yield.

Improved Genome Sequence of Wild Emmer Wheat Zavitan with the Aid of Optical Maps.

Wild emmer (Triticum turgidum ssp. dicoccoides) is the progenitor of all modern cultivated tetraploid wheat. Its genome is large (> 10 Gb) and contains over 80% repeated sequences. The successful whole-genome-shotgun assembly of the wild emmer (accession Zavitan) genome sequence (WEW_v1.0) was an important milestone for wheat genomics. In an effort to improve this assembly, an optical map of accession Zavitan was constructed using Bionano Direct Label and Stain (DLS) technology. The map spanned 10.4 Gb. This map and another map produced earlier by us with the Bionano's Nick Label Repair and Stain (NLRS) technology were used to improve the current wild emmer assembly. The WEW_v1.0 assembly consisted of 151,912 scaffolds. Of them, 3,102 could be confidently aligned on the optical maps. Forty-seven were chimeric. They were disjoined and new scaffolds were assembled with the aid of the optical maps. The total number of scaffolds was reduced from 151,912 to 149,252 and N50 increased from 6.96 Mb to 72.63 Mb. Of the 149,252 scaffolds, 485 scaffolds, which accounted for 97% of the total genome length, were aligned and oriented on genetic maps, and new WEW_v2.0 pseudomolecules were constructed. The new pseudomolecules included 333 scaffolds (68.51 Mb) which were originally unassigned, 226 scaffolds (554.84 Mb) were placed into new locations, and 332 scaffolds (394.83 Mb) were re-oriented. The improved wild emmer genome assembly is an important resource for understanding genomic modification that occurred by domestication.

Genome Based Meta-QTL Analysis of Grain Weight in Tetraploid Wheat Identifies Rare Alleles of GRF4 Associated with Larger Grains.

The domestication and subsequent genetic improvement of wheat led to the development of large-seeded cultivated wheat species relative to their smaller-seeded wild progenitors. While increased grain weight (GW) continues to be an important goal of many wheat breeding programs, few genes underlying this trait have been identified despite an abundance of studies reporting quantitative trait loci (QTL) for GW. Here we perform a QTL analysis for GW using a population of recombinant inbred lines (RILs) derived from the cross between wild emmer wheat accession 'Zavitan' and durum wheat variety 'Svevo'. Identified QTLs in this population were anchored to the recent Zavitan reference genome, along with previously published QTLs for GW in tetraploid wheat. This genome-based, meta-QTL analysis enabled the identification of a locus on chromosome 6A whose introgression from wild wheat positively affects GW. The locus was validated using an introgression line carrying the 6A GW QTL region from Zavitan in a Svevo background, resulting in >8% increase in GW compared to Svevo. Using the reference sequence for the 6A QTL region, we identified a wheat ortholog to OsGRF4, a rice gene previously associated with GW. The coding sequence of this gene (TtGRF4-A) contains four single nucleotide polymorphisms (SNPs) between Zavitan and Svevo, one of which reveals the Zavitan allele to be rare in a core collection of wild emmer and completely absent from the domesticated emmer genepool. Similarly, another wild emmer accession (G18-16) was found to carry a rare allele of TtGRF4-A that also positively affects GW and is characterized by a unique SNP absent from the entire core collection. These results exemplify the rich genetic diversity of wild wheat, posit TtGRF4-A as a candidate gene underlying the 6A GW QTL, and suggest that the natural Zavitan and G18-16 alleles of TtGRF4-A have potential to increase wheat yields in breeding programs.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family.

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

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.

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.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides.

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene Yr26.

KEY MESSAGE: NGS-assisted super pooling emerging as powerful tool to accelerate gene mapping and haplotype association analysis within target region uncovering specific linkage SNPs or alleles for marker-assisted gene pyramiding. Conventional gene mapping methods to identify genes associated with important agronomic traits require significant amounts of financial support and time. Here, a single nucleotide polymorphism (SNP)-based mapping approach, RNA-Seq and SNP array assisted super pooling analysis, was used for rapid mining of a candidate genomic region for stripe rust resistance gene Yr26 that has been widely used in wheat breeding programs in China. Large DNA and RNA super-pools were genotyped by Wheat SNP Array and sequenced by Illumina HiSeq, respectively. Hundreds of thousands of SNPs were identified and then filtered by multiple filtering criteria. Among selected SNPs, over 900 were found within an overlapping interval of less than 30 Mb as the Yr26 candidate genomic region in the centromeric region of chromosome arm 1BL. The 235 chromosome-specific SNPs were converted into KASP assays to validate the Yr26 interval in different genetic populations. Using a high-resolution mapping population (> 30,000 gametes), we confined Yr26 to a 0.003-cM interval. The Yr26 target region was anchored to the common wheat IWGSC RefSeq v1.0 and wild emmer WEWSeq v.1.0 sequences, from which 488 and 454 kb fragments were obtained. Several candidate genes were identified in the target genomic region, but there was no typical resistance gene in either genome region. Haplotype analysis identified specific SNPs linked to Yr26 and developed robust and breeder-friendly KASP markers. This integration strategy can be applied to accelerate generating many markers closely linked to target genes/QTL for a trait of interest in wheat and other polyploid species.

A glycosyl transferase family 43 protein involved in xylan biosynthesis is associated with straw digestibility in Brachypodium distachyon.

The recalcitrance of secondary plant cell walls to digestion constrains biomass use for the production of sustainable bioproducts and for animal feed. We screened a population of Brachypodium recombinant inbred lines (RILs) for cell wall digestibility using commercial cellulases and detected a quantitative trait locus (QTL) associated with this trait. Examination of the chromosomal region associated with this QTL revealed a candidate gene that encodes a putative glycosyl transferase family (GT) 43 protein, orthologue of IRX14 in Arabidopsis, and hence predicted to be involved in the biosynthesis of xylan. Arabinoxylans form the major matrix polysaccharides in cell walls of grasses, such as Brachypodium. The parental lines of the RIL population carry alternative nonsynonymous polymorphisms in the BdGT43A gene, which were inherited in the RIL progeny in a manner compatible with a causative role in the variation in straw digestibility. In order to validate the implied role of our candidate gene in affecting straw digestibility, we used RNA interference to lower the expression levels of the BdGT43A gene in Brachypodium. The biomass of the silenced lines showed higher digestibility supporting a causative role of the BdGT43A gene, suggesting that it might form a good target for improving straw digestibility in crops.