Identification of a major dominant gene for race-nonspecific tan spot resistance in wild emmer wheat.

KEY MESSAGE: A single dominant gene found in tetraploid and hexaploid wheat controls broad-spectrum race-nonspecific resistance to the foliar disease tan spot caused by Pyrenophora tritici-repentis. Tan spot is an important foliar disease of durum and common wheat caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis. Genetic studies in common wheat have shown that pathogen-produced necrotrophic effectors interact with host genes in an inverse gene-for-gene manner to cause disease, but quantitative trait loci (QTLs) with broad race-nonspecific resistance also exist. Less work has been done to understand the genetics of tan spot interactions in durum wheat. Here, we evaluated a set of Langdon durum-wild emmer (Triticum turgidum ssp. dicoccoides) disomic chromosome substitution lines for reaction to four P. tritici-repentis isolates representing races 1, 2, 3, and 5 to identify wild emmer chromosomes potentially containing tan spot resistance genes. Chromosome 3B from the wild emmer accession IsraelA rendered the tan spot-susceptible durum cultivar Langdon resistant to all four fungal isolates. Genetic analysis indicated that a single dominant gene, designated Tsr7, governed resistance. Detailed mapping experiments showed that the Tsr7 locus is likely the same as the race-nonspecific QTL previously identified in the hexaploid wheat cultivars BR34 and Penawawa. Four user-friendly SNP-based semi-thermal asymmetric reverse PCR (STARP) markers cosegregated with Tsr7 and should be useful for marker-assisted selection of resistance. In addition to 3B, other wild emmer chromosomes contributed moderate levels of tan spot resistance, and, as has been shown previously for tetraploid wheat, the Tsn1-Ptr ToxA interaction was not associated with susceptibility. This is the first report of a major dominant gene governing resistance to tan spot in tetraploid wheat.

Durum wheat genome highlights past domestication signatures and future improvement targets.

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Identification and Validation of a New Source of Low Grain Cadmium Accumulation in Durum Wheat.

Cadmium (Cd) is a heavy metal that has no known biological function and is toxic for many living organisms. The maximum level of Cd concentration allowed in the international market for wheat grain is 0.2 mg kg-1 Because phenotyping for Cd uptake is expensive and time consuming, molecular markers associated with genes conferring low Cd uptake would expedite selection and lead to the development of durum cultivars with reduced Cd concentrations. Here, we identified single nucleotide polymorphisms (SNPs) associated with a novel low Cd uptake locus in the durum experimental line D041735, which has hexaploid common wheat in its pedigree. Genetic analysis revealed a single major QTL for Cd uptake on chromosome arm 5BL within a 0.3 cM interval flanked by SNP markers. Analysis of the intervening sequence revealed a gene with homology to an aluminum-induced protein as a candidate gene. Validation and allelism tests revealed that the low Cd uptake gene identified in this study is different from the closely linked Cdu1-B gene, which also resides on 5BL. This study therefore showed that the durum experimental line D041735 contains a novel low Cd uptake gene that was likely acquired from hexaploid wheat.

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.

Molecular cloning and characterization of two novel genes from hexaploid wheat that encode double PR-1 domains coupled with a receptor-like protein kinase.

Hexaploid wheat (Triticum aestivum L.) contains at least 23 TaPr-1 genes encoding the group 1 pathogenesis-related (PR-1) proteins as identified in our previous work. Here, we report the cloning and characterization of TaPr-1-rk1 and TaPr-1-rk2, two novel genes closely related to the wheat PR-1 family. The two TaPr-1-rk genes are located on homoeologous chromosomes 3D and 3A, respectively, and each contains a large open reading frame (7385 or 6060 bp) that is interrupted by seven introns and subjected to alternative splicing (AS) with five or six isoforms of mRNA transcripts. The deduced full-length TaPR-1-RK1 and TaPR-1-RK2 proteins (95% identity) contain two repeat PR-1 domains, the second of which is fused via a transmembrane helix to a serine/threonine kinase catalytic (STKc) domain characteristic of receptor-like protein kinases. Phylogenetic analysis indicated that the two PR-1 domains of the TaPR-1-RK proteins form sister clades with their homologues identified in other monocot plants and are well separated from stand-alone PR-1 proteins, whereas the STKc domains may have originated from cysteine-rich receptor-like kinases (CRKs). Reverse-transcriptase-PCR analysis revealed that the TaPr-1-rk genes are predominantly expressed in wheat leaves and their expression levels are elevated in response to pathogen attack, such as infection by barley stripe mosaic virus (BSMV), and also to stress conditions, most obviously, to soil salinity. This is the first report of PR-1-CRK hybrid proteins in wheat. The data may shed new insights into the function/evolutionary origin of the PR-1 family and the STKc-mediated defense/stress response pathways in plants.

Molecular and comparative mapping of genes governing spike compactness from wild emmer wheat.

The development and morphology of the wheat spike is important because the spike is where reproduction occurs and it holds the grains until harvest. Therefore, genes that influence spike morphology are of interest from both theoretical and practical stand points. When substituted for the native chromosome 2A in the tetraploid Langdon (LDN) durum wheat background, the Triticum turgidum ssp. dicoccoides chromosome 2A from accession IsraelA confers a short, compact spike with fewer spikelets per spike compared to LDN. Molecular mapping and quantitative trait loci (QTL) analysis of these traits in a homozygous recombinant population derived from LDN × the chromosome 2A substitution line (LDNIsA-2A) indicated that the number of spikelets per spike and spike length were controlled by linked, but different, loci on the long arm of 2A. A QTL explaining most of the variation for spike compactness coincided with the QTL for spike length. Comparative mapping indicated that the QTL for number of spikelets per spike overlapped with a previously mapped QTL for Fusarium head blight susceptibility. The genes governing spike length and compactness were not orthologous to either sog or C, genes known to confer compact spikes in diploid and hexaploid wheat, respectively. Mapping and sequence analysis indicated that the gene governing spike length and compactness derived from wild emmer could be an ortholog of the barley Cly1/Zeo gene, which research indicates is an AP2-like gene pleiotropically affecting cleistogamy, flowering time, and rachis internode length. This work provides researchers with knowledge of new genetic loci and associated markers that may be useful for manipulating spike morphology in durum wheat.

Map-based analysis of the tenacious glume gene Tg-B1 of wild emmer and its role in wheat domestication.

The domestication of wheat was instrumental in spawning the civilization of humankind, and it occurred through genetic mutations that gave rise to types with non-fragile rachises, soft glumes, and free-threshing seed. Wild emmer (Triticum turgidum ssp. dicoccoides), the tetraploid AB-genome progenitor of domesticated wheat has genes that confer tenacious glumes (Tg) that underwent genetic mutations to give rise to free-threshing wheat. Here, we evaluated disomic substitution lines involving chromosomes 2A and 2B of wild emmer accessions substituted for homologous chromosomes in tetraploid and hexaploid backgrounds. The results suggested that both chromosomes 2A and 2B of wild emmer possess genes that inhibit threshability. A population of recombinant inbred lines derived from the tetraploid durum wheat variety Langdon crossed with a Langdon - T. turgidum ssp. dicoccoides accession PI 481521 chromosome 2B disomic substitution line was used to develop a genetic linkage map of 2B, evaluate the genetics of threshability, and map the gene derived from PI 481521 that inhibited threshability. A 2BS linkage map comprised of 58 markers was developed, and markers delineated the gene to a 2.3 cM interval. Comparative analysis with maps containing the tenacious glume gene Tg-D1 on chromosome arm 2 DS from Aegilops tauschii, the D genome progenitor of hexaploid wheat, revealed that the gene inhibiting threshability in wild emmer was homoeologous to Tg-D1 and therefore designated Tg-B1. Comparative analysis with rice and Brachypodium distachyon indicated a high level of divergence and poorly conserved colinearity, particularly near the Tg-B1 locus. These results provide a foundation for further studies involving Tg-B1, which, together with Tg-D1, had profound influences on wheat domestication.

Dimerization and protease resistance: new insight into the function of PR-1.

The group 1 pathogenesis-related (PR-1) proteins have long been considered hallmarks of hypersensitive response/defense pathways in plants, but their biochemical functions are still obscure despite resolution of the NMR/X-ray structures of several PR-1-like proteins, including P14a (the prototype PR-1). We report here the characterization of two basic PR-1 proteins (PR-1-1 and PR-1-5) recently identified from hexaploid wheat (Triticum aestivum). Both proteins were expressed in Pichia pastoris as a single major species of ∼15 kDa. Sequence identity of the expressed PR-1 proteins was verified by MALDI-TOF/TOF analysis. Accumulation of the native PR-1-5 protein in pathogen-challenged wheat was confirmed by protein gel blot analysis. Low-temperature SDS-PAGE and yeast two-hybrid assays revealed that PR-1-1 exists primarily as a monomer whereas PR-1-5 forms homodimers. Both PR-1 proteins are resistant to proteases compared to bovine serum albumin, but PR-1-1 shows resistance mainly to subtilisin and protease K (serine proteases) whereas PR-1-5 shows resistance to subtilisin, protease K and papain (a cysteine protease). Site-specific mutations at the five putative active sites in the PR-1 domain all affected dimerization, with the mutations at Glu-72 and Glu-102 (in the PR-1-5 numeration) also diminishing protease resistance. Sequence analysis revealed that the Glu-72 and Glu-102 residues are located in motif-like sequences that are conserved in both PR-1 and the human apoptosis-related caspase proteins. These findings prompt us to examine the function of PR-1 for a role in protease-mediated programmed cell death pathways in plants.

Comparative analysis of syntenic genes in grass genomes reveals accelerated rates of gene structure and coding sequence evolution in polyploid wheat.

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%-25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

The cysteine rich necrotrophic effector SnTox1 produced by Stagonospora nodorum triggers susceptibility of wheat lines harboring Snn1.

The wheat pathogen Stagonospora nodorum produces multiple necrotrophic effectors (also called host-selective toxins) that promote disease by interacting with corresponding host sensitivity gene products. SnTox1 was the first necrotrophic effector identified in S. nodorum, and was shown to induce necrosis on wheat lines carrying Snn1. Here, we report the molecular cloning and validation of SnTox1 as well as the preliminary characterization of the mechanism underlying the SnTox1-Snn1 interaction which leads to susceptibility. SnTox1 was identified using bioinformatics tools and verified by heterologous expression in Pichia pastoris. SnTox1 encodes a 117 amino acid protein with the first 17 amino acids predicted as a signal peptide, and strikingly, the mature protein contains 16 cysteine residues, a common feature for some avirulence effectors. The transformation of SnTox1 into an avirulent S. nodorum isolate was sufficient to make the strain pathogenic. Additionally, the deletion of SnTox1 in virulent isolates rendered the SnTox1 mutated strains avirulent on the Snn1 differential wheat line. SnTox1 was present in 85% of a global collection of S. nodorum isolates. We identified a total of 11 protein isoforms and found evidence for strong diversifying selection operating on SnTox1. The SnTox1-Snn1 interaction results in an oxidative burst, DNA laddering, and pathogenesis related (PR) gene expression, all hallmarks of a defense response. In the absence of light, the development of SnTox1-induced necrosis and disease symptoms were completely blocked. By comparing the infection processes of a GFP-tagged avirulent isolate and the same isolate transformed with SnTox1, we conclude that SnTox1 may play a critical role during fungal penetration. This research further demonstrates that necrotrophic fungal pathogens utilize small effector proteins to exploit plant resistance pathways for their colonization, which provides important insights into the molecular basis of the wheat-S. nodorum interaction, an emerging model for necrotrophic pathosystems.

Duplication and partitioning in evolution and function of homoeologous Q loci governing domestication characters in polyploid wheat.

The Q gene encodes an AP2-like transcription factor that played an important role in domestication of polyploid wheat. The chromosome 5A Q alleles (5AQ and 5Aq) have been well studied, but much less is known about the q alleles on wheat homoeologous chromosomes 5B (5Bq) and 5D (5Dq). We investigated the organization, evolution, and function of the Q/q homoeoalleles in hexaploid wheat (Triticum aestivum L.). Q/q gene sequences are highly conserved within and among the A, B, and D genomes of hexaploid wheat, the A and B genomes of tetraploid wheat, and the A, S, and D genomes of the diploid progenitors, but the intergenic regions of the Q/q locus are highly divergent among homoeologous genomes. Duplication of the q gene 5.8 Mya was likely followed by selective loss of one of the copies from the A genome progenitor and the other copy from the B, D, and S genomes. A recent V(329)-to-I mutation in the A lineage is correlated with the Q phenotype. The 5Bq homoeoalleles became a pseudogene after allotetraploidization. Expression analysis indicated that the homoeoalleles are coregulated in a complex manner. Combined phenotypic and expression analysis indicated that, whereas 5AQ plays a major role in conferring domestication-related traits, 5Dq contributes directly and 5Bq indirectly to suppression of the speltoid phenotype. The evolution of the Q/q loci in polyploid wheat resulted in the hyperfunctionalization of 5AQ, pseudogenization of 5Bq, and subfunctionalization of 5Dq, all contributing to the domestication traits.

A unique wheat disease resistance-like gene governs effector-triggered susceptibility to necrotrophic pathogens.

Plant disease resistance is often conferred by genes with nucleotide binding site (NBS) and leucine-rich repeat (LRR) or serine/threonine protein kinase (S/TPK) domains. Much less is known about mechanisms of susceptibility, particularly to necrotrophic fungal pathogens. The pathogens that cause the diseases tan spot and Stagonospora nodorum blotch on wheat produce effectors (host-selective toxins) that induce susceptibility in wheat lines harboring corresponding toxin sensitivity genes. The effector ToxA is produced by both pathogens, and sensitivity to ToxA is governed by the Tsn1 gene on wheat chromosome arm 5BL. Here, we report the cloning of Tsn1, which was found to have disease resistance gene-like features, including S/TPK and NBS-LRR domains. Mutagenesis revealed that all three domains are required for ToxA sensitivity, and hence disease susceptibility. Tsn1 is unique to ToxA-sensitive genotypes, and insensitive genotypes are null. Sequencing and phylogenetic analysis indicated that Tsn1 arose in the B-genome diploid progenitor of polyploid wheat through a gene-fusion event that gave rise to its unique structure. Although Tsn1 is necessary to mediate ToxA recognition, yeast two-hybrid experiments suggested that the Tsn1 protein does not interact directly with ToxA. Tsn1 transcription is tightly regulated by the circadian clock and light, providing further evidence that Tsn1-ToxA interactions are associated with photosynthesis pathways. This work suggests that these necrotrophic pathogens may thrive by subverting the resistance mechanisms acquired by plants to combat other pathogens.

SnTox3 acts in effector triggered susceptibility to induce disease on wheat carrying the Snn3 gene.

The necrotrophic fungus Stagonospora nodorum produces multiple proteinaceous host-selective toxins (HSTs) which act in effector triggered susceptibility. Here, we report the molecular cloning and functional characterization of the SnTox3-encoding gene, designated SnTox3, as well as the initial characterization of the SnTox3 protein. SnTox3 is a 693 bp intron-free gene with little obvious homology to other known genes. The predicted immature SnTox3 protein is 25.8 kDa in size. A 20 amino acid signal sequence as well as a possible pro sequence are predicted. Six cysteine residues are predicted to form disulfide bonds and are shown to be important for SnTox3 activity. Using heterologous expression in Pichia pastoris and transformation into an avirulent S. nodorum isolate, we show that SnTox3 encodes the SnTox3 protein and that SnTox3 interacts with the wheat susceptibility gene Snn3. In addition, the avirulent S. nodorum isolate transformed with SnTox3 was virulent on host lines expressing the Snn3 gene. SnTox3-disrupted mutants were deficient in the production of SnTox3 and avirulent on the Snn3 differential wheat line BG220. An analysis of genetic diversity revealed that SnTox3 is present in 60.1% of a worldwide collection of 923 isolates and occurs as eleven nucleotide haplotypes resulting in four amino acid haplotypes. The cloning of SnTox3 provides a fundamental tool for the investigation of the S. nodorum-wheat interaction, as well as vital information for the general characterization of necrotroph-plant interactions.

Identification and characterization of segregation distortion loci along chromosome 5B in tetraploid wheat.

Segregation distortion genes are widespread in plants and animals and function by their effect on competition among gametes for preferential fertilization. In this study, we evaluated the segregation distortion of molecular markers in multiple reciprocal backcross populations derived from unique cytogenetic stocks involving the durum cultivar Langdon (LDN) and wild emmer accessions that allowed us to study the effects of chromosome 5B in isolation. No segregation distortion of female gametes was observed, but three populations developed to analyze segregation of male gametes had genomic regions containing markers with skewed segregation ratios. One region of distortion was due to preferential transmission of LDN alleles over wild emmer alleles through male gametes. Another region required the presence of LDN 5B chromosomes in the female for preferential fertilization by male gametes harboring LDN alleles indicating that the corresponding genes in the female gametes can govern genes affecting segregation distortion of male gametes. A third region of distortion was the result of preferential transmission of wild emmer alleles over LDN alleles through male gametes. These results indicate the existence of different distorter/meiotic drive elements among different genotypes and show that distortion factors along wheat chromosome 5B differ in chromosomal location as well as underlying mechanisms.