Domestication evolution,genetics and genomics in wheat

Domestication of plants and animals is the major factor underlying human civilization and is a gigantic evolutionary experiment of adaptation and speciation, generating incipient species. Wheat is one of the most important grain crops in the world, and consists mainly of two types: the hexaploid bread wheat (Triticum aestivum) accounting for about 95% of world wheat production, and the tetraploid durum wheat (T. durum) accounting for the other 5%. In this review, we summarize and discuss research on wheat domestication, mainly focusing on recent findings in genetics and genomics studies. T. aestivum originated from a cross between domesticated emmer wheat T. dicoccum and the goat grass Aegilops tauschii, most probably in the south and west of the Caspian Sea about 9,000 years ago. Wild emmer wheat has the same genome formula as durum wheat and has contributed two genomes to bread wheat, and is central to wheat domestication. Domestication has genetically not only transformed the brittle rachis, tenacious glume and non-free threshability, but also modified yield and yield components in wheat. Wheat domestication involves a limited number of chromosome regions, or domestication syndrome factors, though many relevant quantitative trait loci have been detected. On completion of the genome sequencing of diploid wild wheat (T. urartu or Ae. tauschii), domestication syndrome factors and other relevant genes could be isolated, and effects of wheat domestication could be determined. The achievements of domestication genetics and robust research programs in Triticeae genomics are of greatly help in conservation and exploitation of wheat germplasm and genetic improvement of wheat cultivars.     

Patterns and impacting factors of gene evolutionary rate between wild and cultivated emmer wheat (Triticum turgidum)

Tetraploid emmer wheat (Triticum turgidum L., BBAA) is the founder progenitor of bread wheat, providing the valuable genetic resource and gene pool for wheat improvement. However, the evolutionary trajectory of tetraploid wheat, especially the evolutionary fate of different types of genes has not been well studied. In this study, the rate of non-synonymous substitution (d_N) and synonymous substitution (d_S) was calculated by comparing the orthologs between the wild emmer and cultivated durum wheat at the whole genome and subgenome levels to obtain the positively selected genes (PSGs) and negatively selected genes (NSGs). Then, mutation rate, gene length, exon number, GC content, codon bias, and expression level were comprehensively investigated and compared between the PSGs and NSGs. Within both wild emmer and cultivated durum wheat, PSGs between A and B subgenome displayed shorter gene and exon lengths as well as fewer exon numbers compared with NSGs, whereas from wild emmer to cultivated durum wheat, PSGs showed longer gene length and more exon numbers. Furthermore, PSGs displayed much higher expression levels and stronger codon usage bias, but lower genetic diversity compared with NSGs. Finally, two PSGs TdER1-6B, and TdLC7-2A, were found to play the crucial roles in regulating grain width and plant height of tetraploid wheat, respectively. This study systematically investigated the evolutionary, structural, and functional difference between PSGs and NSGs in tetraploid wheat, which will contribute to a better understanding of the selective mode and evolutionary trajectory during wheat domestication and evolution.     

Identification and genetic characterization of an <Emphasis Type="Italic">Aegilops tauschii</Emphasis> ortholog of the wheat leaf rust disease resistance gene <Emphasis Type="Italic">Lr1</Emphasis>

Aegilops tauschii (goat grass) is the progenitor of the D genome in hexaploid bread wheat. We have screened more than 200 Ae. tauschii accessions for resistance against leaf rust (Puccinia triticina) isolates, which are avirulent on the leaf rust resistance gene Lr1. Approximately 3.5% of the Ae. tauschii accessions displayed the same low infection type as the tester line Thatcher Lr1. The accession Tr.t. 213, which showed resistance after artificial infection with Lr1 isolates both in Mexico and in Switzerland, was chosen for further analysis. Genetic analysis showed that the resistance in this accession is controlled by a single dominant gene, which mapped at the same chromosomal position as Lr1 in wheat. It was delimited in a 1.3-cM region between the restriction fragment length polymorphism (RFLP) markers ABC718 and PSR567 on chromosome 5DL of Ae. tauschii. The gene was more tightly linked to PSR567 (0.47 cM) than to ABC718 (0.79 cM). These results indicate that the resistance gene in Ae. tauschii accession Tr.t. 213 is an ortholog of the leaf rust resistance gene Lr1 of bread wheat, suggesting that Lr1 originally evolved in diploid goat grass and was introgressed into the wheat D genome during or after domestication of hexaploid wheat. Compared to hexaploid wheat, higher marker polymorphism and recombination frequencies were observed in the region of the Lr1 ortholog in Ae. tauschii. The identification of Lr1^Ae, the orthologous gene of wheat Lr1, in Ae. tauschii will allow map-based cloning of Lr1 from this genetically simpler, diploid genome.Hong-Qing Ling and Jiwen Qiu have contributed equally to this work     

Wheat genetic diversity trends during domestication and breeding

It has been claimed that plant breeding reduces genetic diversity in elite germplasm which could seriously jeopardize the continued ability to improve crops. The main objective of this study was to examine the loss of genetic diversity in spring bread wheat during (1) its domestication, (2) the change from traditional landrace cultivars (LCs) to modern breeding varieties, and (3) 50 years of international breeding. We studied 253 CIMMYT or CIMMYT-related modern wheat cultivars, LCs, and Triticum tauschii accessions, the D-genome donor of wheat, with 90 simple sequence repeat (SSR) markers dispersed across the wheat genome. A loss of genetic diversity was observed from T. tauschii to the LCs, and from the LCs to the elite breeding germplasm. Wheats genetic diversity was narrowed from 1950 to 1989, but was enhanced from 1990 to 1997. Our results indicate that breeders averted the narrowing of the wheat germplasm base and subsequently increased the genetic diversity through the introgression of novel materials. The LCs and T. tauschii contain numerous unique alleles that were absent in modern spring bread wheat cultivars. Consequently, both the LCs and T. tauschii represent useful sources for broadening the genetic base of elite wheat breeding germplasm.     

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome‐based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow‐sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA^Lys and tRNA^Met species, while the nonrepetitive assembly reveals tRNA^Ala species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome‐specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome‐level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.     

Fine genetic mapping of greenbug aphid-resistance gene <Emphasis Type="Italic">Gb3</Emphasis> in <Emphasis Type="Italic">Aegilops tauschii</Emphasis>

The greenbug, Schizaphis graminum (Rondani), is an important aphid pest of small grain crops especially wheat (Triticum aestivum L., 2n = 6x = 42, genomes AABBDD) in many parts of the world. The greenbug-resistance gene Gb3 originated from Aegilops tauschii Coss. (2n = 2x = 14, genome D^tD^t) has shown consistent and durable resistance against prevailing greenbug biotypes in wheat fields. We previously mapped Gb3 in a recombination-rich, telomeric bin of wheat chromosome arm 7DL. In this study, high-resolution genetic mapping was carried out using an F_2:3 segregating population derived from two Ae. tauschii accessions, the resistant PI 268210 (original donor of Gb3 in the hexaploid wheat germplasm line ‘Largo’) and susceptible AL8/78. Molecular markers were developed by exploring bin-mapped wheat RFLPs, SSRs, ESTs and the Ae. tauschii physical map (BAC contigs). Wheat EST and Ae. tauschii BAC end sequences located in the deletion bin 7DL3-0.82–1.00 were used to design STS (sequence tagged site) or CAPS (Cleaved Amplified Polymorphic Sequence) markers. Forty-five PCR-based markers were developed and mapped to the chromosomal region spanning the Gb3 locus. The greenbug-resistance gene Gb3 now was delimited in an interval of 1.1 cM by two molecular markers (HI067J6-R and HI009B3-R). This localized high-resolution genetic map with markers closely linked to Gb3 lays a solid foundation for map based cloning of Gb3 and marker-assisted selection of this gene in wheat breeding.     

Development and QTL assessment of Triticum aestivum–Aegilops tauschii introgression lines

A set of 84 bread wheat lines, each containing a single homozygous introgression of the Aegilops tauschii genome was produced in the ‘Chinese Spring’ background via backcrossing of the D-genome chromosome substitution lines ‘Chinese Spring’/Sears’s ‘Synthetic 6x’ with the recurrent parent and subsequent selfing. The development of the lines was accompanied by microsatellite marker assisted selection. With the exception of three telomeric regions at chromosomes 1DL, 4DL and 7DS, and a region of less than 24 cM on the chromosome arm 3DL, the genome of Ae. tauschii is fully represented in these lines. The newly developed lines were used for the discovery of morphological and agronomical quantitative trait loci (QTLs) from the wild species. Fifty-two introgression lines were grown in the field and evaluated for six traits including flowering time, plant height, ear length, spikelet number, fertility and grain weight per ear. Seventeen significant QTLs were detected, Ae. tauschii contributed favourable alleles at nine loci influencing five traits. The whole set of 84 homozygous lines provides a tool for further testing the effects and stability of the detected QTLs and for the evaluation of new traits.     

Comparative analysis of the D genome-encoded high-molecular weight subunits of glutenin

Four genes encoding novel 1Dx-type high-molecular weight (HMW) subunits were amplified by polymerase chain reaction, two each from Aegilops tauschii and bread wheat Triticum aestivum. The two subunits from Ae. tauschii (1Dx2.1^t and 1Dx2^t) were both very similar in sequence to subunit 1Dx2 from bread wheat. In contrast, the two novel bread wheat subunits (1Dx2.2 and 1Dx2.2*) differed from subunit 1Dx2 in having different internally duplicated regions (of 132 and 186 amino acid, respectively) within their repetitive domains. These duplicated sequences were located adjacent to the regions from which they had been duplicated and had complete intact repeat motifs at each end. The implications of these results for HMW subunit evolution and wheat quality improvement are discussed.     

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.     

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole‐genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co‐cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low‐density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info .     

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

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

Sequencing of 15 622 gene‐bearing BACs clarifies the gene‐dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole‐genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene‐containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical‐mapped gene‐bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene‐enriched BACs and are characterized by high recombination rates, there are also gene‐dense regions with suppressed recombination. We made use of published map‐anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D‐genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley–Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map‐based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene‐dense but low recombination is particularly relevant.     

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.     

A high-density genetic linkage map of Aegilops tauschii, the D-genome progenitor of bread wheat

 Aegilops tauschii is the diploid D-genome progenitor of bread wheat (Triticum aestivum L. em Thell, 2n=6x=42, AABBDD). A genetic linkage map of the Ae. tauschii genome was constructed, composed of 546 loci. One hundred and thirty two loci (24%) gave distorted segregation ratios. Sixty nine probes (13%) detected multiple copies in the genome. One hundred and twenty three of the 157 markers shared between the Ae. tauschii genetic and T. aestivum physical maps were colinear. The discrepancy in the order of five markers on the Ae. tauschii 3DS genetic map versus the T. aestivum 3D physical map indicated a possible inversion. Further work is needed to verify the discrepancies in the order of markers on the 4D, 5D and 7D Ae. tauschii genetic maps versus the physical and genetic maps of T. aestivum. Using common markers, 164 agronomically important genes were assigned to specific regions on Ae. tauschii linkage, and T. aestivum physical, maps. This information may be useful for map-based cloning and marker-assisted plant breeding. Received: 23 March 1998 / Accepted: 27 October 1998     

Physical mapping and identification of a candidate for the leaf rust resistance gene <Emphasis Type="Italic">Lr1</Emphasis> of wheat

Lr1 is a dominant leaf rust resistance gene located on chromosome 5DL of bread wheat and the wild species Aegilops tauschii. In this study, three polymorphic markers (WR001, WR002, and WR003) were developed from resistance gene analogs (RGAs) clustering around the Lr1 locus. Using these and other markers, Lr1 was mapped to a genetic interval of 0.79 cM in Ae. tauschii and 0.075 cM in wheat. The CAPS marker WR003, derived from LR1RGA1, co-segregated with Lr1 in both mapping populations of wheat and Ae. tauschii. For isolation of Lr1, two genomic BAC libraries (from Ae. tauschii and hexaploid wheat) were screened using the tightly flanking marker PSR567F and a set of nested primers derived from the conserved region of the RGA sequences. Approximately 400 kb BAC contig spanning the Lr1 locus was constructed. The LR1RGA1 encoding a CC-NBS-leucine-rich repeat (LRR) type of protein was the only one of the four RGAs at the Lr1 locus, which co-segregated with leaf rust resistance. Therefore, it represents a very good candidate for Lr1. The allelic sequences of LR1RGA1 from resistant and susceptible lines revealed a divergent DNA sequence block of ∼605 bp encoding the LRR repeats 9–15, whereas the rest of the sequences were mostly identical. Within this sequence block, the 48 non-synonymous changes resulted in 44 amino acid differences. This indicates that LR1RGA1 likely evolved through one or more recombination or gene conversion events with unknown genes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development of simple sequence repeat markers specific for the Lr34 resistance region of wheat using sequence information from rice and Aegilops tauschii

Hexaploid wheat (Triticum aestivum L.) originated about 8,000 years ago from the hybridization of tetraploid wheat with diploid Aegilops tauschii Coss. containing the D-genome. Thus, the bread wheat D-genome is evolutionary young and shows a low degree of polymorphism in the bread wheat gene pool. To increase marker density around the durable leaf rust resistance gene Lr34 located on chromosome 7DS, we used molecular information from the orthologous region in rice. Wheat expressed sequence tags (wESTs) were identified by homology with the rice genes in the interval of interest, but were monomorphic in the ‘Arina’ × ‘Forno’ mapping population. To derive new polymorphic markers, bacterial artificial chromosome (BAC) clones representing a total physical size of ∼1 Mb and belonging to four contigs were isolated from Ae. tauschii by hybridization screening with wheat ESTs. Several BAC clones were low-pass sequenced, resulting in a total of ∼560 kb of sequence. Ten microsatellite sequences were found, and three of them were polymorphic in our population and were genetically mapped close to Lr34. Comparative analysis of marker order revealed a large inversion between the rice genome and the wheat D-genome. The SWM10 microsatellite is closely linked to Lr34 and has the same allele in the three independent sources of Lr34: ‘Frontana’, ‘Chinese Spring’, and ‘Forno’, as well in most of the genotypes containing Lr34. Therefore, SWM10 is a highly useful marker to assist selection for Lr34 in breeding programs worldwide.     

Molecular characterization and genomic organization of low molecular weight glutenin subunit genes at the Glu-3 loci in hexaploid wheat (Triticum aestivum L.)

In this study, we report on the molecular characterization and genomic organization of the low molecular weight glutenin subunit (LMW-GS) gene family in hexaploid wheat (Triticum aestivum L.). Eighty-two positive BAC clones were identified to contain LMW-GS genes from the hexaploid wheat ‘Glenlea’ BAC library via filter hybridization and PCR validation. Twelve unique LMW glutenin genes and seven pseudogenes were isolated from these positive BAC clones by primer-template mismatch PCR and subsequent primer walking using hemi-nested touchdown PCR. These genes were sequenced and each consisted of a single-open reading frame (ORF) and untranslated 5′ and 3′ flanking regions. All 12 LMW glutenin subunits contained eight cysteine residues. The LMW-m-type subunits are the most abundant in hexaploid wheat. Of the 12 LMW-GS, 1, 2 and 9 are i-type, s-type and m-type, respectively. The phylogenetic analysis suggested that the LMW-i type gene showed greater differences to LMW-s and LMW-m-type genes, which, in turn, were more closely related to one another. On the basis of their N-terminal sequences, they were classified into nine groups. Fingerprinting of the 82 BAC clones indicated 30 BAC clones assembled into eight contigs, while the remaining clones were singletons. BAC end sequencing of the 82 clones revealed that long terminal repeat (LTR) retrotransposons were abundant in the Glu-3 regions. The average physical distance between two adjacent LMW-GS genes was estimated to be 81 kb. Most of LMW-GS genes are located in the d-genome, suggesting that the Glu-D3 locus is much larger than the Glu-B3 locus and Glu-A3 locus. Alignments of sequences indicated that the same type (starting with the same N-terminal sequence) LMW-GS genes were highly conserved in the homologous genomes between hexaploid wheat and its donors such as durum wheat and T. tauschii. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.