Genome sequence of the progenitor of the wheat D genome Aegilops tauschii.

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

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.

Dynamic evolution of resistance gene analogs in the orthologous genomic regions of powdery mildew resistance gene MlIW170 in Triticum dicoccoides and Aegilops tauschii.

KEY MESSAGE: Rapid evolution of powdery mildew resistance gene MlIW170 orthologous genomic regions in wheat subgenomes. Wheat is one of the most important staple grain crops in the world and also an excellent model for plant ploidy evolution research with different ploidy levels from diploid to hexaploid. Powdery mildew disease caused by Blumeria graminis f.sp. tritici can result in significant loss in both grain yield and quality in wheat. In this study, the wheat powdery mildew resistance gene MlIW170 locus located at the Triticum dicoccoides chromosome 2B short arm was further characterized by constructing and sequencing a BAC-based physical map contig covering a 0.3 cM genetic distance region (880 kb) and developing additional markers to delineate the resistance gene within a 0.16 cM genetic interval (372 kb). Comparative analyses of the T. dicoccoides 2BS region with the orthologous Aegilops tauschii 2DS region showed great gene colinearity, including the structure organization of both types of RGA1/2-like and RPS2-like resistance genes. Comparative analyses with the orthologous regions from Brachypodium and rice genomes revealed considerable dynamic evolutionary changes that have re-shaped this MlIW170 region in the wheat genome, resulting in a high number of non-syntenic genes including resistance-related genes. This result might reflect the rapid evolution in R-gene regions. Phylogenetic analysis on these resistance-related gene sequences indicated the duplication of these genes in the MlIW170 region, occurred before the separation of the wheat B and D genomes.

High-density genetic linkage map construction and QTL mapping of grain shape and size in the wheat population Yanda1817 × Beinong6.

High-density genetic linkage maps are necessary for precisely mapping quantitative trait loci (QTLs) controlling grain shape and size in wheat. By applying the Infinium iSelect 9K SNP assay, we have constructed a high-density genetic linkage map with 269 F 8 recombinant inbred lines (RILs) developed between a Chinese cornerstone wheat breeding parental line Yanda1817 and a high-yielding line Beinong6. The map contains 2431 SNPs and 128 SSR & EST-SSR markers in a total coverage of 3213.2 cM with an average interval of 1.26 cM per marker. Eighty-eight QTLs for thousand-grain weight (TGW), grain length (GL), grain width (GW) and grain thickness (GT) were detected in nine ecological environments (Beijing, Shijiazhuang and Kaifeng) during five years between 2010-2014 by inclusive composite interval mapping (ICIM) (LOD ≥ 2.5). Among which, 17 QTLs for TGW were mapped on chromosomes 1A, 1B, 2A, 2B, 3A, 3B, 3D, 4A, 4D, 5A, 5B and 6B with phenotypic variations ranging from 2.62% to 12.08%. Four stable QTLs for TGW could be detected in five and seven environments, respectively. Thirty-two QTLs for GL were mapped on chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 4D, 5A, 5B, 6B, 7A and 7B, with phenotypic variations ranging from 2.62% to 44.39%. QGl.cau-2A.2 can be detected in all the environments with the largest phenotypic variations, indicating that it is a major and stable QTL. For GW, 12 QTLs were identified with phenotypic variations range from 3.69% to 12.30%. We found 27 QTLs for GT with phenotypic variations ranged from 2.55% to 36.42%. In particular, QTL QGt.cau-5A.1 with phenotypic variations of 6.82-23.59% was detected in all the nine environments. Moreover, pleiotropic effects were detected for several QTL loci responsible for grain shape and size that could serve as target regions for fine mapping and marker assisted selection in wheat breeding programs.

Comparative genetic mapping and genomic region collinearity analysis of the powdery mildew resistance gene Pm41.

KEY MESSAGE: By applying comparative genomics analyses, a high-density genetic linkage map narrowed the powdery mildew resistance gene Pm41 originating from wild emmer in a sub-centimorgan genetic interval. Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici, results in large yield losses worldwide. A high-density genetic linkage map of the powdery mildew resistance gene Pm41, originating from wild emmer (Triticum turgidum var. dicoccoides) and previously mapped to the distal region of chromosome 3BL bin 0.63-1.00, was constructed using an F5:6 recombinant inbred line population derived from a cross of durum wheat cultivar Langdon and wild emmer accession IW2. By applying comparative genomics analyses, 19 polymorphic sequence-tagged site markers were developed and integrated into the Pm41 genetic linkage map. Ultimately, Pm41 was mapped in a 0.6 cM genetic interval flanked by markers XWGGC1505 and XWGGC1507, which correspond to 11.7, 19.2, and 24.9 kb orthologous genomic regions in Brachypodium, rice, and sorghum, respectively. The XWGGC1506 marker co-segregated with Pm41 and could be served as a starting point for chromosome landing and map-based cloning as well as marker-assisted selection of Pm41. Detailed comparative genomics analysis of the markers flanking the Pm41 locus in wheat and the putative orthologous genes in Brachypodium, rice, and sorghum suggests that the gene order is highly conserved between rice and sorghum. However, intra-chromosome inversions and re-arrangements are evident in the wheat and Brachypodium genomic regions, and gene duplications are also present in the orthologous genomic regions of Pm41 in wheat, indicating that the Brachypodium gene model can provide more useful information for wheat marker development.

Fine physical and genetic mapping of powdery mildew resistance gene MlIW172 originating from wild emmer (Triticum dicoccoides).

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most important wheat diseases in the world. In this study, a single dominant powdery mildew resistance gene MlIW172 was identified in the IW172 wild emmer accession and mapped to the distal region of chromosome arm 7AL (bin7AL-16-0.86-0.90) via molecular marker analysis. MlIW172 was closely linked with the RFLP probe Xpsr680-derived STS marker Xmag2185 and the EST markers BE405531 and BE637476. This suggested that MlIW172 might be allelic to the Pm1 locus or a new locus closely linked to Pm1. By screening genomic BAC library of durum wheat cv. Langdon and 7AL-specific BAC library of hexaploid wheat cv. Chinese Spring, and after analyzing genome scaffolds of Triticum urartu containing the marker sequences, additional markers were developed to construct a fine genetic linkage map on the MlIW172 locus region and to delineate the resistance gene within a 0.48 cM interval. Comparative genetics analyses using ESTs and RFLP probe sequences flanking the MlIW172 region against other grass species revealed a general co-linearity in this region with the orthologous genomic regions of rice chromosome 6, Brachypodium chromosome 1, and sorghum chromosome 10. However, orthologous resistance gene-like RGA sequences were only present in wheat and Brachypodium. The BAC contigs and sequence scaffolds that we have developed provide a framework for the physical mapping and map-based cloning of MlIW172.

Comparative high-resolution mapping of the wax inhibitors Iw1 and Iw2 in hexaploid wheat.

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (W1 and W2) and non-glaucousness loci (Iw1 and Iw2). The non-glaucousness (Iw) loci act as inhibitors of the glaucousness loci (W). High-resolution comparative genetic linkage maps of the wax inhibitors Iw1 originating from Triticum dicoccoides, and Iw2 from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic sequences corresponding to the syntenic regions of the Iw loci in wheat. Eleven Iw1 and eight Iw2 linked EST markers were developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iw1 locus mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iw1 locus on chromosome 2BS. Both Iw1 and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical mapping and map-based cloning of the wax inhibitors in wheat.