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.

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.

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

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.

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

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 5S^se, which targeted them for recombination. We genotyped 28 BC₁F₁ and 84 F₂ 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 BC₁F₁ and F₂ genetic maps that were 65.31 and 63.71 cM long, respectively. Recombination rates between the 5B and 5S^se chromosomes were double the expected rate computed from their meiotic pairing, which we attributed to selection against aneuploid gametes. Recombination rate between 5B and 5S^se 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.

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.

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.

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

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 F₂ 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.

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.

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.

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.

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

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.

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.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

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.

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

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.

Exploring the metabolic variation between domesticated and wild tetraploid wheat genotypes in response to corn leaf aphid infestation

Infestation of Triticum (wheat) plants by their pest Rhopalosiphum maidis (corn leaf aphid) causes severe vegetative damage. Despite the agro-economic importance of wheat, the metabolic diversity of Triticum turgidum (tetraploid wheat) in response to aphid attack has not been sufficiently addressed. In this study, we compared the metabolic diversity of two tetraploid wheat genotypes, domesticated and wild emmer. The plants were grown in a control growth room and infested with aphids for 96 h. Our untargeted metabolic analysis performed on plants with and without aphids revealed massive differences between the two genotypes. The targeted metabolic analysis highlighted the differences in the biosynthesis of phytohormones. The aphid progeny was lower in the cultivated durum wheat than in the wild emmer wheat. Overall, these observations emphasize the potential of using the natural diversity of wheat species to better understand the metabolic responses to pest damage.

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.

Wild emmer genome architecture and diversity elucidate wheat evolution and domestication

Wheat (Triticum spp.) is one of the founder crops that likely drove the Neolithic transition to sedentary agrarian societies in the Fertile Crescent more than 10,000 years ago. Identifying genetic modifications underlying wheat's domestication requires knowledge about the genome of its allo-tetraploid progenitor, wild emmer (T. turgidum ssp. dicoccoides). We report a 10.1-gigabase assembly of the 14 chromosomes of wild tetraploid wheat, as well as analyses of gene content, genome architecture, and genetic diversity. With this fully assembled polyploid wheat genome, we identified the causal mutations in Brittle Rachis 1 (TtBtr1) genes controlling shattering, a key domestication trait. A study of genomic diversity among wild and domesticated accessions revealed genomic regions bearing the signature of selection under domestication. This reference assembly will serve as a resource for accelerating the genome-assisted improvement of modern wheat varieties.

On the Origin of the Non-brittle Rachis Trait of Domesticated Einkorn Wheat

Einkorn and emmer wheat together with barley were among the first cereals domesticated by humans more than 10,000 years ago, long before durum or bread wheat originated. Domesticated einkorn wheat differs from its wild progenitor in basic morphological characters such as the grain dispersal system. This study identified the Non-brittle rachis 1 (btr1) and Non-brittle rachis 2 (btr2) in einkorn as homologous to barley. Re-sequencing of the Btr1 and Btr2 in a collection of 53 lines showed that a single non-synonymous amino acid substitution (alanine to threonine) at position 119 at btr1, is responsible for the non-brittle rachis trait in domesticated einkorn. Tracing this haplotype variation back to wild einkorn samples provides further evidence that the einkorn progenitor came from the Northern Levant. We show that the geographical origin of domesticated haplotype coincides with the non-brittle domesticated barley haplotypes, which suggest the non-brittle rachis phenotypes of einkorn and barley were fixed in same geographic area in today's South-east Turkey.

Introgression of the Aegilops speltoides Su1-Ph1 Suppressor into Wheat

Meiotic pairing between homoeologous chromosomes in polyploid wheat is inhibited by the Ph1 locus on the long arm of chromosome 5 in the B genome. Aegilops speltoides (genomes SS), the closest relative of the progenitor of the wheat B genome, is polymorphic for genetic suppression of Ph1. Using this polymorphism, two major suppressor loci, Su1-Ph1 and Su2-Ph1, have been mapped in Ae. speltoides. Su1-Ph1 is located in the distal, high-recombination region of the long arm of the Ae. speltoides chromosome 3S. Its location and tight linkage to marker Xpsr1205-3S makes Su1-Ph1 a suitable target for introgression into wheat. Here, Xpsr1205-3S was introgressed into hexaploid bread wheat cv. Chinese Spring (CS) and from there into tetraploid durum wheat cv. Langdon (LDN). Sequential fluorescence in situ hybridization and genomic in situ hybridization showed that an Ae. speltoides segment with Xpsr1205-3S replaced the distal end of the long arm of chromosome 3A. In the CS genetic background, the chromosome induced homoeologous chromosome pairing in interspecific hybrids with Ae. peregrina but not in progenies from crosses involving alien disomic substitution lines. In the LDN genetic background, the chromosome induced homoeologous chromosome pairing in both interspecific hybrids and progenies from crosses involving alien disomic substitution lines. We conclude that the recombined chromosome harbors Su1-Ph1 but its expression requires expression of complementary gene that is present in LDN but absent in CS. We suggest that it is unlikely that Su1-Ph1 and ZIP4-1, a paralog of Ph1 located on wheat chromosomes 3A and 3B and Ae. tauschii chromosome 3D, are equivalent. The utility of Su1-Ph1 for induction of recombination between homoeologous chromosomes in wheat is illustrated.

A High-Density Genetic Map of Wild Emmer Wheat from the Karaca Dağ Region Provides New Evidence on the Structure and Evolution of Wheat Chromosomes

Wild emmer (Triticum turgidum ssp. dicoccoides) is a progenitor of all cultivated wheat grown today. It has been hypothesized that emmer was domesticated in the Karaca Dağ region in southeastern Turkey. A total of 445 recombinant inbred lines of T. turgidum ssp. durum cv. 'Langdon' x wild emmer accession PI 428082 from this region was developed and genotyped with the Illumina 90K single nucleotide polymorphism Infinium assay. A genetic map comprising 2,650 segregating markers was constructed. The order of the segregating markers and an additional 8,264 co-segregating markers in the Aegilops tauschii reference genome sequence was used to compare synteny of the tetraploid wheat with the Brachypodium distachyon, rice, and sorghum. These comparisons revealed the presence of 15 structural chromosome rearrangements, in addition to the already known 4A-5A-7B rearrangements. The most common type was an intra-chromosomal translocation in which the translocated segment was short and was translocated only a short distance along the chromosome. A large reciprocal translocation, one small non-reciprocal translocation, and three large and one small paracentric inversions were also discovered. The use of inversions for a phylogeny reconstruction in the Triticum-Aegilops alliance was illustrated. The genetic map was inconsistent with the current model of evolution of the rearranged chromosomes 4A-5A-7B. Genetic diversity in the rearranged chromosome 4A showed that the rearrangements might have been contemporary with wild emmer speciation. A selective sweep was found in the centromeric region of chromosome 4A in Karaca Dağ wild emmer but not in 4A of T. aestivum. The absence of diversity from a large portion of chromosome 4A of wild emmer, believed to be ancestral to all domesticated wheat, is puzzling.

Haplotype Analysis of the Pre-harvest Sprouting Resistance Locus Phs-A1 Reveals a Causal Role of TaMKK3-A in Global Germplasm

Pre-harvest sprouting (PHS) is an important cause of quality loss in many cereal crops and is particularly prevalent and damaging in wheat. Resistance to PHS is therefore a valuable target trait in many breeding programs. The Phs-A1 locus on wheat chromosome arm 4AL has been consistently shown to account for a significant proportion of natural variation to PHS in diverse mapping populations. However, the deployment of sprouting resistance is confounded by the fact that different candidate genes, including the tandem duplicated Plasma Membrane 19 (PM19) genes and the mitogen-activated protein kinase kinase 3 (TaMKK3-A) gene, have been proposed to underlie Phs-A1. To further define the Phs-A1 locus, we constructed a physical map across this interval in hexaploid and tetraploid wheat. We established close proximity of the proposed candidate genes which are located within a 1.2 Mb interval. Genetic characterization of diverse germplasm used in previous genetic mapping studies suggests that TaMKK3-A, and not PM19, is the major gene underlying the Phs-A1 effect in European, North American, Australian and Asian germplasm. We identified the non-dormant TaMKK3-A allele at low frequencies within the A-genome diploid progenitor Triticum urartu genepool, and show an increase in the allele frequency in modern varieties. In United Kingdom varieties, the frequency of the dormant TaMKK3-A allele was significantly higher in bread-making quality varieties compared to feed and biscuit-making cultivars. Analysis of exome capture data from 58 diverse hexaploid wheat accessions identified fourteen haplotypes across the extended Phs-A1 locus and four haplotypes for TaMKK3-A. Analysis of these haplotypes in a collection of United Kingdom and Australian cultivars revealed distinct major dormant and non-dormant Phs-A1 haplotypes in each country, which were either rare or absent in the opposing germplasm set. The diagnostic markers and haplotype information reported in the study will help inform the choice of germplasm and breeding strategies for the deployment of Phs-A1 resistance into breeding germplasm.

The Solanum tuberosum KST1 partial promoter as a tool for guard cell expression in multiple plant species

To date, guard cell promoters have been examined in only a few species, primarily annual dicots. A partial segment of the potato (Solanum tuberosum) KST1 promoter (KST1 partial promoter, KST1ppro) has previously been shown to confer guard cell expression in potato, tomato (Solanum lycopersicum), citrus [Troyer citrange (C. sinensis×Poncirus trifoliata)], and Arabidopsis (Arabidopsis thaliana). Here, we describe an extensive analysis of the expression pattern of KST1ppro in eight (previously reported, as well as new) species from five different angiosperm families, including the Solanaceae and the Cucurbitaceae, Arabidopsis, the monocot barley (Hordeum vulgare), and two perennial species: grapevine (Vitis vinifera) and citrus. Using confocal imaging and three-dimensional movies, we demonstrate that KST1ppro drives guard cell expression in all of these species, making it the first dicot-originated guard cell promoter shown to be active in a monocot and the first promoter reported to confer guard cell expression in barley and cucumber (Cucumis sativus). The results presented here indicate that KST1ppro can be used to drive constitutive guard cell expression in monocots and dicots and in both annual and perennial plants. In addition, we show that the KST1ppro is active in guard cells shortly after the symmetric division of the guard mother cell and generates stable expression in mature guard cells. This allows us to follow the spatial and temporal distribution of stomata in cotyledons and true leaves.

Acceleration of leaf senescence is slowed down in transgenic barley plants deficient in the DNA/RNA-binding protein WHIRLY1

WHIRLY1 in barley was isolated as a potential regulator of the senescence-associated gene HvS40. In order to investigate whether the plastid-nucleus-located DNA/RNA-binding protein WHIRLY1 plays a role in regulation of leaf senescence, primary foliage leaves from transgenic barley plants with an RNAi-mediated knockdown of the WHIRLY1 gene were characterized by typical senescence parameters, namely pigment contents, function and composition of the photosynthetic apparatus, as well as expression of selected genes known to be either down- or up-regulated during leaf senescence. When the plants were grown at low light intensity, senescence progression was similar between wild-type and RNAi-W1 plants. Likewise, dark-induced senescence of detached leaves was not affected by reduction of WHIRLY1. When plants were grown at high light intensity, however, senescence was induced prematurely in wild-type plants but was delayed in RNAi-W1 plants. This result suggests that WHIRLY1 plays a role in light sensing and/or stress communication between chloroplasts and the nucleus.

QTLs for uniform grain dimensions and germination selected during wheat domestication are co-located on chromosome 4B

A major locus on the long arm of wheat chromosome 4B controls within-spikelet variation in both grain size and seed dormancy, the latter an important survival mechanism likely eliminated from wild wheat during domestication. Seed dormancy can increase the probability of survival of at least some progeny under unstable environmental conditions. In wild emmer wheat, only one of the two grains in a spikelet germinates during the first rainy season following maturation; and this within-plant variation in seed dormancy is associated with both grain dimension differences and position within the spikelet. Here, in addition to characterizing these associations, we elucidate the genetic mechanism controlling differential grain dimensions and dormancy within wild tetraploid wheat spikelets using phenotypic data from a wild emmer × durum wheat population and a high-density genetic map. We show that in wild emmer, the lower grain within the spikelet is about 30 % smaller and more dormant than the larger, upper grain that germinates usually within 3 days. We identify a major locus on the long arm of chromosome 4B that explains >40 % of the observed variation in grain dimensions and seed dormancy within spikelets. This locus, designated QGD-4BL, is validated using an independent set of wild emmer × durum wheat genetic stocks. The domesticated variant of this novel locus on chromosome 4B, likely fixed during the process of wheat domestication, favors spikelets with seeds of uniform size and synchronous germination. The identification of locus QGD-4BL enhances our knowledge of the genetic basis of the domestication syndrome of one of our most important crops.

A Metabolic Gene Cluster in the Wheat W1 and the Barley Cer-cqu Loci Determines β-Diketone Biosynthesis and Glaucousness

The glaucous appearance of wheat (Triticum aestivum) and barley (Hordeum vulgare) plants, that is the light bluish-gray look of flag leaf, stem, and spike surfaces, results from deposition of cuticular β-diketone wax on their surfaces; this phenotype is associated with high yield, especially under drought conditions. Despite extensive genetic and biochemical characterization, the molecular genetic basis underlying the biosynthesis of β-diketones remains unclear. Here, we discovered that the wheat W1 locus contains a metabolic gene cluster mediating β-diketone biosynthesis. The cluster comprises genes encoding proteins of several families including type-III polyketide synthases, hydrolases, and cytochrome P450s related to known fatty acid hydroxylases. The cluster region was identified in both genetic and physical maps of glaucous and glossy tetraploid wheat, demonstrating entirely different haplotypes in these accessions. Complementary evidence obtained through gene silencing in planta and heterologous expression in bacteria supports a model for a β-diketone biosynthesis pathway involving members of these three protein families. Mutations in homologous genes were identified in the barley eceriferum mutants defective in β-diketone biosynthesis, demonstrating a gene cluster also in the β-diketone biosynthesis Cer-cqu locus in barley. Hence, our findings open new opportunities to breed major cereal crops for surface features that impact yield and stress response.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Regulation of Zn and Fe transporters by the GPC1 gene during early wheat monocarpic senescence

BACKGROUND

During wheat senescence, leaf components are degraded in a coordinated manner, releasing amino acids and micronutrients which are subsequently transported to the developing grain. We have previously shown that the simultaneous downregulation of Grain Protein Content (GPC) transcription factors, GPC1 and GPC2, greatly delays senescence and disrupts nutrient remobilization, and therefore provide a valuable entry point to identify genes involved in micronutrient transport to the wheat grain.

RESULTS

We generated loss-of-function mutations for GPC1 and GPC2 in tetraploid wheat and showed in field trials that gpc1 mutants exhibit significant delays in senescence and reductions in grain Zn and Fe content, but that mutations in GPC2 had no significant effect on these traits. An RNA-seq study of these mutants at different time points showed a larger proportion of senescence-regulated genes among the GPC1 (64%) than among the GPC2 (37%) regulated genes. Combined, the two GPC genes regulate a subset (21.2%) of the senescence-regulated genes, 76.1% of which are upregulated at 12 days after anthesis, before the appearance of any visible signs of senescence. Taken together, these results demonstrate that GPC1 is a key regulator of nutrient remobilization which acts predominantly during the early stages of senescence. Genes upregulated at this stage include transporters from the ZIP and YSL gene families, which facilitate Zn and Fe export from the cytoplasm to the phloem, and genes involved in the biosynthesis of chelators that facilitate the phloem-based transport of these nutrients to the grains.

CONCLUSIONS

This study provides an overview of the transport mechanisms activated in the wheat flag leaf during monocarpic senescence. It also identifies promising targets to improve nutrient remobilization to the wheat grain, which can help mitigate Zn and Fe deficiencies that afflict many regions of the developing world.

Introgression of leaf rust and stripe rust resistance from Sharon goatgrass (Aegilops sharonensis Eig) into bread wheat (Triticum aestivum L.)

Leaf rust and stripe rust are devastating wheat diseases, causing significant yield losses in many regions of the world. The use of resistant varieties is the most efficient way to protect wheat crops from these diseases. Sharon goatgrass (Aegilops sharonensis or AES), which is a diploid wild relative of wheat, exhibits a high frequency of leaf and stripe rust resistance. We used the resistant AES accession TH548 and induced homoeologous recombination by the ph1b allele to obtain resistant wheat recombinant lines carrying AES chromosome segments in the genetic background of the spring wheat cultivar Galil. The gametocidal effect from AES was overcome by using an "anti-gametocidal" wheat mutant. These recombinant lines were found resistant to highly virulent races of the leaf and stripe rust pathogens in Israel and the United States. Molecular DArT analysis of the different recombinant lines revealed different lengths of AES segments on wheat chromosome 6B, which indicates the location of both resistance genes.

Senescence, nutrient remobilization, and yield in wheat and barley

Cereals including wheat and barley are of primary importance to ensure food security for the 21st century. A combination of lab- and field-based approaches has led to a considerably improved understanding of the importance of organ and particularly of whole-plant (monocarpic) senescence for wheat and barley yield and quality. A delicate balance between senescence timing, grain nutrient content, nutrient-use efficiency, and yield needs to be considered to (further) improve cereal varieties for a given environment and end use. The recent characterization of the Gpc-1 (NAM-1) genes in wheat and barley demonstrates the interdependence of these traits. Lines or varieties with functional Gpc-1 genes demonstrate earlier senescence and enhanced grain protein and micronutrient content but, depending on the environment, somewhat reduced yields. A major effort is needed to dissect regulatory networks centred on additional wheat and barley transcription factors and signalling pathways influencing the senescence process. Similarly, while important molecular details of nutrient (particularly nitrogen) remobilization from senescing organs to developing grains have been identified, important knowledge gaps remain. The genes coding for the major proteases involved in senescence-associated plastidial protein degradation are largely unknown. Membrane transport proteins involved in the different transport steps occurring between senescing organ (such as leaf mesophyll) cells and protein bodies in the endosperm of developing grains remain to be identified or further characterized. Existing data suggest that an improved understanding of all these steps will reveal additional, important targets for continued cereal improvement.

Functional characterization of GPC-1 genes in hexaploid wheat

In wheat, monocarpic senescence is a tightly regulated process during which nitrogen (N) and micronutrients stored pre-anthesis are remobilized from vegetative tissues to the developing grains. Recently, a close connection between senescence and remobilization was shown through the map-based cloning of the GPC (grain protein content) gene in wheat. GPC-B1 encodes a NAC transcription factor associated with earlier senescence and increased grain protein, iron and zinc content, and is deleted or non-functional in most commercial wheat varieties. In the current research, we identified 'loss of function' ethyl methanesulfonate mutants for the two GPC-B1 homoeologous genes; GPC-A1 and GPC-D1, in a hexaploid wheat mutant population. The single gpc-a1 and gpc-d1 mutants, the double gpc-1 mutant and control lines were grown under field conditions at four locations and were characterized for senescence, GPC, micronutrients and yield parameters. Our results show a significant delay in senescence in both the gpc-a1 and gpc-d1 single mutants and an even stronger effect in the gpc-1 double mutant in all the environments tested in this study. The accumulation of total N in the developing grains showed a similar increase in the control and gpc-1 plants until 25 days after anthesis (DAA) but at 41 and 60 DAA the control plants had higher grain N content than the gpc-1 mutants. At maturity, GPC in all mutants was significantly lower than in control plants while grain weight was unaffected. These results demonstrate that the GPC-A1 and GPC-D1 genes have a redundant function and play a major role in the regulation of monocarpic senescence and nutrient remobilization in wheat.

Increased copy number at the HvFT1 locus is associated with accelerated flowering time in barley

A precise regulation of flowering time is critical for plant reproductive success, and therefore, a better understanding of the natural variation in genes regulating the initiation of the reproductive phase is required to develop well-adapted varieties. In both monocot and dicot species, the FLOWERING LOCUS T (FT) is a central integrator of seasonal signals perceived by the leaves. The encoded mobile protein (florigen) is transmitted to the apical meristem where it induces flowering. The FT homolog in barley (Hordeum vulgare L.), designated HvFT1, was shown to correspond to the vernalization locus VRN-H3, and natural alleles for spring and winter growth habit were identified. In this study, we demonstrate that the HvFT1 allele present in the barley genetic stock (BGS213) associated with a dominant spring growth habit carries at least four identical copies of HvFT1, whereas most barley varieties have a single copy. Increased copy number is associated with earlier transcriptional up-regulation of HvFT1 and a spring growth habit. This allele is epistatic to winter alleles for VRN-H1 and VRN-H2. Among accessions with one HvFT1 copy, haplotype differences in the HvFT1 promoter and first intron are also associated with differences in flowering time, which are modulated by genetic background. These different HvFT1 alleles can be used to develop barley varieties adapted to different or changing environments. Our results, together with studies of other wheat and barley flowering genes, show that copy number variation plays an important role in the regulation of developmental processes in the temperate cereals.

Identification and characterization of a novel powdery mildew resistance gene PmG3M derived from wild emmer wheat, Triticum dicoccoides

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most important wheat diseases worldwide. Wild emmer wheat, Triticum turgidum ssp. dicoccoides, the tetraploid ancestor (AABB) of domesticated bread and durum wheat, harbors many important alleles for resistance to various diseases, including powdery mildew. In the current study, two tetraploid wheat mapping populations, derived from a cross between durum wheat (cv. Langdon) and wild emmer wheat (accession G-305-3M), were used to identify and map a novel powdery mildew resistance gene. Wild emmer accession G-305-3M was resistant to all 47 Bgt isolates tested, from Israel and Switzerland. Segregation ratios of F(2) progenies and F(6) recombinant inbred line (RIL) mapping populations, in their reactions to inoculation with Bgt, revealed a Mendelian pattern (3:1 and 1:1, respectively), indicating the role of a single dominant gene derived from T. dicoccoides accession G-305-3M. This gene, temporarily designated PmG3M, was mapped on chromosome 6BL and physically assigned to chromosome deletion bin 6BL-0.70-1.00. The F(2) mapping population was used to construct a genetic map of the PmG3M gene region consisted of six simple sequence repeats (SSR), 11 resistance gene analog (RGA), and two target region amplification polymorphism (TRAP) markers. A second map, constructed based on the F(6) RIL population, using a set of skeleton SSR markers, confirmed the order of loci and distances obtained for the F(2) population. The discovery and mapping of this novel powdery mildew resistance gene emphasize the importance of the wild emmer wheat gene pool as a source for crop improvement.

Divergent functions of orthologous NAC transcription factors in wheat and rice

The wheat GPC-B1 gene located on chromosome 6B is an early regulator of senescence and affects remobilization of protein and minerals to the grain. GPC-B1 is a NAC transcription factor and has a paralogous copy on chromosome 2B in wheat, GPC-B2. The closest rice homolog to both wheat GPC genes is Os07g37920 which is located on rice chromosome 7 and is colinear with GPC-B2. Since rice is a diploid species with a sequenced genome, we initiated the study of Os07g37920 to develop a simpler model to study senescence and mineral remobilization in cereals. We developed eleven independent RNA interference transgenic rice lines (Os07g37920-RNAi) and 10 over-expressing transgenic lines (Os07g37920-OE), but none of them showed differences in senescence. Transgenic Os07g37920-RNAi rice plants had reduced proportions of viable pollen grains and were male-sterile, but were able to produce seeds by cross pollination. Analysis of the flower morphology of the transgenic rice plants showed that anthers failed to dehisce. Transgenic Os07g37920-OE lines showed no sterility or anther dehiscence problems. Os07g37920 transcript levels were higher in stamens compared to leaves and significantly reduced in the transgenic Os07g37920-RNAi plants. Wheat GPC genes showed the opposite transcription profile (higher transcript levels in leaves than in flowers) and plants carrying knock-out mutations of all GPC-1 and GPC-2 genes exhibited delayed senescence but normal anther dehiscence and fertility. These results indicate a functional divergence of the homologous wheat and rice NAC genes and suggest the need for separate studies of the function and targets of these transcription factors in wheat and rice.

Effect of the down-regulation of the high Grain Protein Content (GPC) genes on the wheat transcriptome during monocarpic senescence

BACKGROUND

Increasing the nutrient concentration of wheat grains is important to ameliorate nutritional deficiencies in many parts of the world. Proteins and nutrients in the wheat grain are largely derived from the remobilization of degraded leaf molecules during monocarpic senescence. The down-regulation of the NAC transcription factor Grain Protein Content (GPC) in transgenic wheat plants delays senescence (>3 weeks) and reduces the concentration of protein, Zn and Fe in the grain (>30%), linking senescence and nutrient remobilization.Based on the early and rapid up-regulation of GPC in wheat flag leaves after anthesis, we hypothesized that this transcription factor is an early regulator of monocarpic senescence. To test this hypothesis, we used high-throughput mRNA-seq technologies to characterize the effect of the GPC down-regulation on the wheat flag-leaf transcriptome 12 days after anthesis. At this early stage of senescence GPC transcript levels are significantly lower in transgenic GPC-RNAi plants than in the wild type, but there are still no visible phenotypic differences between genotypes.

RESULTS

We generated 1.4 million 454 reads from early senescing flag leaves (average ~350 nt) and assembled 1.2 million into 30,497 contigs that were used as a reference to map 145 million Illumina reads from three wild type and four GPC-RNAi plants. Following normalization and statistical testing, we identified a set of 691 genes differentially regulated by GPC (431 ≥ 2-fold change). Transcript level ratios between transgenic and wild type plants showed a high correlation (R = 0.83) between qRT-PCR and Illumina results, providing independent validation of the mRNA-seq approach. A set of differentially expressed genes were analyzed across an early senescence time-course.

CONCLUSIONS

Monocarpic senescence is an active process characterized by large-scale changes in gene expression which begins considerably before the appearance of visual symptoms of senescence. The mRNA-seq approach used here was able to detect small differences in transcript levels during the early stages of senescence. This resulted in an extensive list of GPC-regulated genes, which includes transporters, hormone regulated genes, and transcription factors. These GPC-regulated genes, particularly those up-regulated during senescence, provide valuable entry points to dissect the early stages of monocarpic senescence and nutrient remobilization in wheat.

Wheat flowering repressor VRN2 and promoter CO2 compete for interactions with NUCLEAR FACTOR-Y complexes

The transition from vegetative to reproductive development in the temperate cereals is mainly regulated by seasonal cues including vernalization (determined mainly by VRN1 and VRN2 genes) and photoperiod (determined mainly by PPD1 and CO2 genes). The wheat VRN3 gene, which is similar to Arabidopsis FT, plays a central role in the integration of the competing signals from these two pathways. Under long days, VRN3 transcription is down-regulated by VRN2, a unique flowering repressor in cereals, and up-regulated by CO2. Overexpression of VRN3 overcomes VRN2 repression and promotes VRN1 transcription and flowering initiation. Using yeast two- and three-hybrid assays we show here that the CCT domains present in VRN2 and CO2 proteins interact with the same subset of eight NF-Y proteins, and that these CCT proteins compete with NF-YA for interactions with NF-YB proteins. We have confirmed all these interactions in vitro, and the interactions between VRN2 and two of the three NF-YB proteins were further confirmed in planta. In addition, we show that mutations in the CCT domain of VRN2 that eliminate the vernalization requirement in winter wheat also reduce the strength of the interactions between VRN2 and NF-Y proteins, and the ability of VRN2 to compete with CO2. Taken together, our results suggest that the interactions between CCT and NF-Y proteins play an important role in the integration of the vernalization and photoperiod seasonal signals, and provide a flexible combinatorial system to integrate multiple developmental and environmental signals in the regulation of flowering initiation in the temperate cereals.

Regulation of freezing tolerance and flowering in temperate cereals: the VRN-1 connection

In winter wheat (Triticum spp.) and barley (Hordeum vulgare) varieties, long exposures to nonfreezing cold temperatures accelerate flowering time (vernalization) and improve freezing tolerance (cold acclimation). However, when plants initiate their reproductive development, freezing tolerance decreases, suggesting a connection between the two processes. To better understand this connection, we used two diploid wheat (Triticum monococcum) mutants, maintained vegetative phase (mvp), that carry deletions encompassing VRN-1, the major vernalization gene in temperate cereals. Homozygous mvp/mvp plants never flower, whereas plants carrying at least one functional VRN-1 copy (Mvp/-) exhibit normal flowering and high transcript levels of VRN-1 under long days. The Mvp/- plants showed reduced freezing tolerance and reduced transcript levels of several cold-induced C-REPEAT BINDING FACTOR transcription factors and COLD REGULATED genes (COR) relative to the mvp/mvp plants. Diploid wheat accessions with mutations in the VRN-1 promoter, resulting in high transcript levels under both long and short days, showed a significant down-regulation of COR14b under long days but not under short days. Taken together, these studies suggest that VRN-1 is required for the initiation of the regulatory cascade that down-regulates the cold acclimation pathway but that additional genes regulated by long days are required for the down-regulation of the COR genes. In addition, our results show that allelic variation in VRN-1 is sufficient to determine differences in freezing tolerance, suggesting that quantitative trait loci for freezing tolerance previously mapped on this chromosome region are likely a pleiotropic effect of VRN-1 rather than the effect of a separate closely linked locus (FROST RESISTANCE-1), as proposed in early freezing tolerance studies.

Small RNAs, DNA methylation and transposable elements in wheat

Background

More than 80% of the wheat genome is composed of transposable elements (TEs). Since active TEs can move to different locations and potentially impose a significant mutational load, their expression is suppressed in the genome via small non-coding RNAs (sRNAs). sRNAs guide silencing of TEs at the transcriptional (mainly 24-nt sRNAs) and post-transcriptional (mainly 21-nt sRNAs) levels. In this study, we report the distribution of these two types of sRNAs among the different classes of wheat TEs, the regions targeted within the TEs, and their impact on the methylation patterns of the targeted regions.

Results

We constructed an sRNA library from hexaploid wheat and developed a database that included our library and three other publicly available sRNA libraries from wheat. For five completely-sequenced wheat BAC contigs, most perfectly matching sRNAs represented TE sequences, suggesting that a large fraction of the wheat sRNAs originated from TEs. An analysis of all wheat TEs present in the Triticeae Repeat Sequence database showed that sRNA abundance was correlated with the estimated number of TEs within each class. Most of the sRNAs perfectly matching miniature inverted repeat transposable elements (MITEs) belonged to the 21-nt class and were mainly targeted to the terminal inverted repeats (TIRs). In contrast, most of the sRNAs matching class I and class II TEs belonged to the 24-nt class and were mainly targeted to the long terminal repeats (LTRs) in the class I TEs and to the terminal repeats in CACTA transposons. An analysis of the mutation frequency in potentially methylated sites revealed a three-fold increase in TE mutation frequency relative to intron and untranslated genic regions. This increase is consistent with wheat TEs being preferentially methylated, likely by sRNA targeting.

Conclusions

Our study examines the wheat epigenome in relation to known TEs. sRNA-directed transcriptional and post-transcriptional silencing plays important roles in the short-term suppression of TEs in the wheat genome, whereas DNA methylation and increased mutation rates may provide a long-term mechanism to inactivate TEs.

Characterization of the maintained vegetative phase deletions from diploid wheat and their effect on VRN2 and FT transcript levels

Allelic differences at the VRN1 (AP1/CAL/FRU), VRN2 (ZCCT) and VRN3 (FT) vernalization genes affect flowering time in wheat. The two maintained vegetative phase (mvp) mutants from Triticum monococcum L., previously reported as carrying a single gene (VRN1) deletion, are incapable of flowering. In this study, we show that both mvp lines have larger deletions that include the genes AGLG1, CYS, PHYC, VRN1 and possibly others. The original mvp deletions were generated in lines that lack the VRN2 gene. Therefore, to study the effect of the mvp deletions on the regulation of VRN2 we generated populations segregating for both genes simultaneously. The two mvp deletions co-segregated with the non-flowering phenotype, but surprisingly, the lines homozygous for the mvp mutations showed reduced transcript levels of both VRN2 and FT relative to the wild type. The VRN1 deletion is an unlikely cause of the down-regulation of VRN2 since VRN2 transcript levels are higher in the fall, before VRN1 is expressed, and are down-regulated by VRN1. Since both VRN2 and FT are regulated by light and photoperiod, their down-regulation in the mvp mutants might be related to the deletion of the PHYC photoreceptor. However, alternative hypotheses including combinations of other genes deleted in the mvp mutants cannot be ruled out. Until the specific gene(s) responsible for the down-regulation of VRN2 and FT and the non-flowering phenotype are precisely identified, it is premature to use these results to postulate alternative flowering models.

Vrn-D4 is a vernalization gene located on the centromeric region of chromosome 5D in hexaploid wheat

Natural variation in wheat requirement of long exposures to cold temperatures to accelerate flowering (vernalization) is mainly controlled by the Vrn-1, Vrn-2, Vrn-3, and Vrn-4 loci. The first three loci have been well characterized, but limited information is available for Vrn-4. So far, natural variation for Vrn-4 has been detected only in the D genome (Vrn-D4), and genetic stocks for this gene are available in Triple Dirk (TDF, hereafter). We detected heterogeneity in the Vrn-1 alleles present in different TDF stocks, which may explain inconsistencies among previous studies. A correct TDF seed stock from Japan carrying recessive vrn-A1, vrn-B1, and vrn-D1 alleles was crossed with three different winter cultivars to generate F(2) mapping populations. Most of the variation in flowering time in these three populations was controlled by a single locus, Vrn-D4, which was mapped within a 1.8 cM interval flanked by markers Xcfd78 and Xbarc205 in the centromeric region of chromosome 5D. A factorial ANOVA for heading time using Vrn-D4 alleles and vernalization as factors showed a significant interaction (P < 0.0001), which confirmed that the Vrn-D4 effect on flowering time is modulated by vernalization. Comparison of the different Triple Dirk stocks revealed that Vrn-B1, Vrn-D1, and Vrn-D4 all have a small residual response to vernalization, but Vrn-D4 differs from the other two in its response to short vernalization periods. The precise mapping and characterization of Vrn-D4 presented here represent a first step toward the positional cloning of this gene.

Regulation of flowering in temperate cereals

Long exposure to cold (vernalization) accelerates flowering in winter cereals, a process regulated by the VRN1 (approximately AP1), VRN2, and VRN3 (approximately FT) vernalization genes. Flowering during the fall is prevented by the VRN2 downregulation of VRN3 and low VRN1 transcription. Vernalization induces VRN1, which is followed by the downregulation of VRN2, thereby releasing VRN3. In the longer days of spring, photoperiod genes PPD1 and CO upregulate VRN3, which induces VRN1 above the threshold levels required for flowering initiation. VRN3 transcription is modulated through interactions involving CCT-domain proteins and HAP2/HAP3/HAP5 complexes coded by multiple genes. The vast number of HAP-CCT combinations can provide the flexibility required for integrating seasonal cues and different stress signals in the regulation of the transition to flowering.

A kinase-START gene confers temperature-dependent resistance to wheat stripe rust

Stripe rust is a devastating fungal disease that afflicts wheat in many regions of the world. New races of Puccinia striiformis, the pathogen responsible for this disease, have overcome most of the known race-specific resistance genes. We report the map-based cloning of the gene Yr36 (WKS1), which confers resistance to a broad spectrum of stripe rust races at relatively high temperatures (25 degrees to 35 degrees C). This gene includes a kinase and a putative START lipid-binding domain. Five independent mutations and transgenic complementation confirmed that both domains are necessary to confer resistance. Yr36 is present in wild wheat but is absent in modern pasta and bread wheat varieties, and therefore it can now be used to improve resistance to stripe rust in a broad set of varieties.

Genetic and molecular characterization of the VRN2 loci in tetraploid wheat

Winter wheat (Triticum spp.) varieties require long exposures to low temperatures to flower, a process called vernalization. The VRN2 locus includes two completely linked zinc finger-CCT domain genes (ZCCT1 and ZCCT2) that act as flowering repressors down-regulated during vernalization. Deletions or mutations in these two genes result in the elimination of the vernalization requirement in diploid wheat (Triticum monococcum). However, natural allelic variation in these genes has not been described so far in polyploid wheat (tetraploid Triticum turgidum and hexaploid Triticum aestivum). A tetraploid wheat population segregating for both VRN-A2 and VRN-B2 loci facilitated the characterization of different alleles. Comparisons between functional and nonfunctional alleles revealed that both ZCCT1 and ZCCT2 genes are able to confer vernalization requirement and that different ZCCT genes are functional in different genomes. ZCCT1 and ZCCT2 proteins from nonfunctional vrn2 alleles have mutations at arginine amino acids at position 16, 35, or 39 of the CCT domain. These positions are conserved between CCT and HEME ACTIVATOR PROTEIN2 (HAP2) proteins, supporting a model in which the action of CCT domains is mediated by their interactions with HAP2/HAP3/HAP5 complexes. This study also revealed natural variation in gene copy number, including a duplication of the functional ZCCT-B2 gene and deletions or duplications of the complete VRN-B2 locus. Allelic variation at the VRN-B2 locus was associated with a partially dominant effect, which suggests that variation in the number of functional ZCCT genes can be used to expand allelic diversity for heading time in polyploid wheat and, hopefully, improve its adaptation to different environments.

Identification of a novel gene (Hsdr4) involved in water-stress tolerance in wild barley

Drought is one of the most severe stresses limiting plant growth and yield. Genes involved in water stress tolerance of wild barley (Hordeum spontaneoum), the progenitor of cultivated barley, were investigated using genotypes contrasting in their response to water stress. Gene expression profiles of water-stress tolerant vs. water-stress sensitive wild barley genotypes, under severe dehydration stress applied at the seedling stage, were compared using cDNA-AFLP analysis. Of the 1100 transcript-derived fragments (TDFs) amplified about 70 displayed differential expression between control and stress conditions. Eleven of them showed clear difference (up- or down-regulation) between tolerant and susceptible genotypes. These TDFs were isolated, sequenced and tested by RT-PCR. The differential expression of seven TDFs was confirmed by RT-PCR, and TDF-4 was selected as a promising candidate gene for water-stress tolerance. The corresponding gene, designated Hsdr4 (Hordeum spontaneum dehydration-responsive), was sequenced and the transcribed and flanking regions were determined. The deduced amino acid sequence has similarity to the rice Rho-GTPase-activating protein-like with a Sec14 p-like lipid-binding domain. Analysis of Hsdr4 promoter region that was isolated by screening a barley BAC library, revealed a new putative miniature inverted-repeat transposable element (MITE), and several potential stress-related binding sites for transcription factors (MYC, MYB, LTRE, and GT-1), suggesting a role of the Hsdr4 gene in plant tolerance to dehydration stress. Furthermore, the Hsdr4 gene was mapped using wild barley mapping population to the long arm of chromosome 3H between markers EBmac541 and EBmag705, within a region that previously was shown to affect osmotic adaptation in barley.

A NAC Gene regulating senescence improves grain protein, zinc, and iron content in wheat

Enhancing the nutritional value of food crops is a means of improving human nutrition and health. We report here the positional cloning of Gpc-B1, a wheat quantitative trait locus associated with increased grain protein, zinc, and iron content. The ancestral wild wheat allele encodes a NAC transcription factor (NAM-B1) that accelerates senescence and increases nutrient remobilization from leaves to developing grains, whereas modern wheat varieties carry a nonfunctional NAM-B1 allele. Reduction in RNA levels of the multiple NAM homologs by RNA interference delayed senescence by more than 3 weeks and reduced wheat grain protein, zinc, and iron content by more than 30%.

Physical map of the wheat high-grain protein content gene Gpc-B1 and development of a high-throughput molecular marker

Grain protein content (GPC) is important for human nutrition and has a strong influence on pasta and bread quality. A quantitative trait locus, derived from a Triticum turgidum ssp. dicoccoides accession (DIC), with an average increase in GPC of 14 g kg(-1) was mapped on chromosome 6BS. Using the wheat-rice colinearity, a high-density map of the wheat region was developed and the quantitative trait locus was mapped as a simple Mendelian locus designated Gpc-B1. A physical map of approx. 250 kb of the Gpc-B1 region was developed using a tetraploid wheat bacterial artificial chromosome library. The constructed physical map included the two Gpc-B1 flanking markers and one potential candidate gene from the colinear rice region completely linked to Gpc-B1. The relationship between physical and genetic distances and the feasibility of isolating genes by positional cloning in wheat are discussed. A high-throughput codominant marker, Xuhw89, was developed. A 4-bp deletion present in the DIC allele was absent in a collection of 117 cultivated tetraploid and hexaploid wheat germplasm, suggesting that this marker will be useful to incorporate the high GPC allele from the DIC accession studied here into commercial wheat varieties.

High-temperature adult-plant (HTAP) stripe rust resistance gene Yr36 from Triticum turgidum ssp. dicoccoides is closely linked to the grain protein content locus Gpc-B1

Several new races of the stripe rust pathogen have become frequent throughout the wheat growing regions of the United States since 2000. These new races are virulent to most of the wheat seedling resistance genes limiting the resistance sources that can be used to combat this pathogen. High-temperature adult-plant (HTAP) stripe rust resistance has proven to be more durable than seedling resistance due to its non-race-specific nature, but its use is limited by the lack of mapping information. We report here the identification of a new HTAP resistance gene from Triticum turgidum ssp. dicoccoides (DIC) designated as Yr36. Lines carrying this gene were susceptible to almost all the stripe rust pathogen races tested at the seedling stage but showed adult-plant resistance to the prevalent races in California when tested at high diurnal temperatures. Isogenic lines for this gene were developed by six backcross generations. Field tests in two locations showed increased levels of field resistance to stripe rust and increased yields in isogenic lines carrying the Yr36 gene compared to those without the gene. Recombinant substitution lines of chromosome 6B from DIC in the isogenic background of durum cv. Langdon were used to map the Yr36 gene on the short arm of chromosome 6B completely linked to Xbarc101, and within a 2-cM interval defined by PCR-based markers Xucw71 and Xbarc136. Flanking locus Xucw71 is also closely linked to the grain protein content locus Gpc-B1 (0.3-cM). Marker-assisted selection strategies are presented to improve stripe rust resistance and simultaneously select for high or low Gpc-B1 alleles.