A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat

Identification of candidate genes, regions and markers for pre-harvest sprouting resistance in wheat (Triticum aestivum L.)

BACKGROUND

Pre-harvest sprouting (PHS) of wheat grain leads to a reduction in grain yield and quality. The availability of markers for marker-assisted selection (MAS) of PHS resistance will serve to enhance breeding selection and advancement of lines for cultivar development. The aim of this study was to identify candidate regions and develop molecular markers for PHS resistance in wheat. This was achieved via high density mapping of single nucleotide polymorphism (SNP) markers from an Illumina 90 K Infinium Custom Beadchip in a doubled haploid (DH) population derived from a RL4452/'AC Domain' cross and subsequent detection of quantitative trait loci (QTL) for PHS related traits (falling number [FN], germination index [GI] and sprouting index [SI]). SNP marker sequences flanking QTL were used to locate colinear regions in Brachypodium and rice, and identify genic markers associated with PHS resistance that can be utilized for MAS in wheat.

RESULTS

A linkage map spanning 2569.4 cM was constructed with a total of 12,201 SNP, simple sequence repeat (SSR), diversity arrays technology (DArT) and expressed sequence tag (EST) markers. QTL analyses using Multiple Interval Mapping (MIM) identified four QTL for PHS resistance traits on chromosomes 3B, 4A, 7B and 7D. Sequences of SNPs flanking these QTL were subject to a BLASTN search on the International Wheat Genome Sequencing Consortium (IWGSC) database (http://wheat-urgi.versailles.inra.fr/Seq-Repository). Best survey sequence hits were subject to a BLASTN search on Gramene (www.gramene.org) against both Brachypodium and rice databases, and candidate genes and regions for PHS resistance were identified. A total of 18 SNP flanking sequences on chromosomes 3B, 4A, 7B and 7D were converted to KASP markers and validated with matching genotype calls of Infinium SNP data.

CONCLUSIONS

Our study identified candidate genes involved in abscissic acid (ABA) and gibberellin (GA) metabolism, and flowering time in four genomic regions of Brachypodium and rice respectively, in addition to 18 KASP markers for PHS resistance in wheat. These markers can be deployed in future genetic studies of PHS resistance and might also be useful in the evaluation of PHS in germplasm and breeding material.

The Na(+) transporter, TaHKT1;5-D, limits shoot Na(+) accumulation in bread wheat

Bread wheat (Triticum aestivum L.) has a major salt tolerance locus, Kna1, responsible for the maintenance of a high cytosolic K(+) /Na(+) ratio in the leaves of salt stressed plants. The Kna1 locus encompasses a large DNA fragment, the distal 14% of chromosome 4DL. Limited recombination has been observed at this locus making it difficult to map genetically and identify the causal gene. Here, we decipher the function of TaHKT1;5-D, a candidate gene underlying the Kna1 locus. Transport studies using the heterologous expression systems Saccharomyces cerevisiae and Xenopus laevis oocytes indicated that TaHKT1;5-D is a Na(+) -selective transporter. Transient expression in Arabidopsis thaliana mesophyll protoplasts and in situ polymerase chain reaction indicated that TaHKT1;5-D is localised on the plasma membrane in the wheat root stele. RNA interference-induced silencing decreased the expression of TaHKT1;5-D in transgenic bread wheat lines which led to an increase in the Na(+) concentration in the leaves. This indicates that TaHKT1;5-D retrieves Na(+) from the xylem vessels in the root and has an important role in restricting the transport of Na(+) from the root to the leaves in bread wheat. Thus, TaHKT1;5-D confers the essential salinity tolerance mechanism in bread wheat associated with the Kna1 locus via shoot Na(+) exclusion and is critical in maintaining a high K(+) /Na(+) ratio in the leaves. These findings show there is potential to increase the salinity tolerance of bread wheat by manipulation of HKT1;5 genes.

Extensive pericentric rearrangements in the bread wheat (Triticum aestivum L.) genotype "Chinese Spring" revealed from chromosome shotgun sequence data

The bread wheat (Triticum aestivum L.) genotype "Chinese Spring" ("CS") is the reference base in wheat genetics and genomics. Pericentric rearrangements in this genotype were systematically assessed by analyzing homoeoloci for a set of nonredundant genes from Brachypodium distachyon, Triticum urartu, and Aegilops tauschii in the CS chromosome shotgun sequence obtained from individual chromosome arms flow-sorted from CS aneuploid lines. Based on patterns of their homoeologous arm locations, 551 genes indicated the presence of pericentric inversions in at least 10 of the 21 chromosomes. Available data from deletion bin-mapped expressed sequence tags and genetic mapping in wheat indicated that all inversions had breakpoints in the low-recombinant gene-poor pericentromeric regions. The large number of putative intrachromosomal rearrangements suggests the presence of extensive structural differences among the three subgenomes, at least some of which likely occurred during the production of the aneuploid lines of this hexaploid wheat genotype. These differences could have significant implications in wheat genome research where comparative approaches are used such as in ordering and orientating sequence contigs and in gene cloning.

Utilization of deletion bins to anchor and order sequences along the wheat 7B chromosome

A total of 3,671 sequence contigs and scaffolds were mapped to deletion bins on wheat chromosome 7B providing a foundation for developing high-resolution integrated physical map for this chromosome. Bread wheat (Triticum aestivum L.) has a large, complex and highly repetitive genome which is challenging to assemble into high quality pseudo-chromosomes. As part of the international effort to sequence the hexaploid bread wheat genome by the international wheat genome sequencing consortium (IWGSC) we are focused on assembling a reference sequence for chromosome 7B. The successful completion of the reference chromosome sequence is highly dependent on the integration of genetic and physical maps. To aid the integration of these two types of maps, we have constructed a high-density deletion bin map of chromosome 7B. Using the 270 K Nimblegen comparative genomic hybridization (CGH) array on a set of cv. Chinese spring deletion lines, a total of 3,671 sequence contigs and scaffolds (~7.8 % of chromosome 7B physical length) were mapped into nine deletion bins. Our method of genotyping deletions on chromosome 7B relied on a model-based clustering algorithm (Mclust) to accurately predict the presence or absence of a given genomic sequence in a deletion line. The bin mapping results were validated using three different approaches, viz. (a) PCR-based amplification of randomly selected bin mapped sequences (b) comparison with previously mapped ESTs and (c) comparison with a 7B genetic map developed in the present study. Validation of the bin mapping results suggested a high accuracy of the assignment of 7B sequence contigs and scaffolds to the 7B deletion bins.

A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome

An ordered draft sequence of the 17-gigabase hexaploid bread wheat (Triticum aestivum) genome has been produced by sequencing isolated chromosome arms. We have annotated 124,201 gene loci distributed nearly evenly across the homeologous chromosomes and subgenomes. Comparative gene analysis of wheat subgenomes and extant diploid and tetraploid wheat relatives showed that high sequence similarity and structural conservation are retained, with limited gene loss, after polyploidization. However, across the genomes there was evidence of dynamic gene gain, loss, and duplication since the divergence of the wheat lineages. A high degree of transcriptional autonomy and no global dominance was found for the subgenomes. These insights into the genome biology of a polyploid crop provide a springboard for faster gene isolation, rapid genetic marker development, and precise breeding to meet the needs of increasing food demand worldwide.

Structural and functional partitioning of bread wheat chromosome 3B

We produced a reference sequence of the 1-gigabase chromosome 3B of hexaploid bread wheat. By sequencing 8452 bacterial artificial chromosomes in pools, we assembled a sequence of 774 megabases carrying 5326 protein-coding genes, 1938 pseudogenes, and 85% of transposable elements. The distribution of structural and functional features along the chromosome revealed partitioning correlated with meiotic recombination. Comparative analyses indicated high wheat-specific inter- and intrachromosomal gene duplication activities that are potential sources of variability for adaption. In addition to providing a better understanding of the organization, function, and evolution of a large and polyploid genome, the availability of a high-quality sequence anchored to genetic maps will accelerate the identification of genes underlying important agronomic traits.

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.

Development of T. aestivum L.-H. californicum alien chromosome lines and assignment of homoeologous groups of Hordeum californicum chromosomes

Hordeum californicum (2n = 2x = 14, HH) is resistant to several wheat diseases and tolerant to lower nitrogen. In this study, a molecular karyotype of H. californicum chromosomes in the Triticum aestivum L. cv. Chinese Spring (CS)-H. californicum amphidiploid (2n = 6x = 56, AABBDDHH) was established. By genomic in situ hybridization (GISH) and multicolor fluorescent in situ hybridization (FISH) using repetitive DNA clones (pTa71, pTa794 and pSc119.2) as probes, the H. californicum chromosomes could be differentiated from each other and from the wheat chromosomes unequivocally. Based on molecular karyotype and marker analyses, 12 wheat-alien chromosome lines, including four disomic addition lines (DAH1, DAH3, DAH5 and DAH6), five telosomic addition lines (MtH7L, MtH1S, MtH1L, DtH6S and DtH6L), one multiple addition line involving H. californicum chromosome H2, one disomic substitution line (DSH4) and one translocation line (TH7S/1BL), were identified from the progenies derived from the crosses of CS-H. californicum amphidiploid with common wheat varieties. A total of 482 EST (expressed sequence tag) or SSR (simple sequence repeat) markers specific for individual H. californicum chromosomes were identified, and 47, 50, 45, 49, 21, 51 and 40 markers were assigned to chromosomes H1, H2, H3, H4, H5, H6 and H7, respectively. According to the chromosome allocation of these markers, chromosomes H2, H3, H4, H5, and H7 of H. californicum have relationship with wheat homoeologous groups 5, 2, 6, 3, and 1, and hence could be designated as 5H(c), 2H(c), 6H(c), 3H(c) and 1H(c), respectively. The chromosomes H1 and H6 were designated as 7H(c) and 4H(c), respectively, by referring to SSR markers located on rye chromosomes.

A novel genetic map of wheat: utility for mapping QTL for yield under different nitrogen treatments

Background

Common wheat (Triticum aestivum L.) is one of the most important food crops worldwide. Wheat varieties that maintain yield (YD) under moderate or even intense nitrogen (N) deficiency can adapt to low input management systems. A detailed genetic map is necessary for both wheat molecular breeding and genomics research. In this study, an F_6:7 recombinant inbred line population comprising 188 lines was used to construct a novel genetic map and subsequently to detect quantitative trait loci (QTL) for YD and response to N stress.

Results

A genetic map consisting of 591 loci distributed across 21 wheat chromosomes was constructed. The map spanned 3930.7 cM, with one marker per 6.7 cM on average. Genomic simple sequence repeat (g-SSR), expressed sequence tag-derived microsatellite (e-SSR), diversity arrays technology (DArT), sequence-tagged sites (STS), sequence-related amplified polymorphism (SRAP), and inter-simple sequence repeat (ISSR) molecular markers were included in the map. The linear relationships between loci found in the present map and in previously compiled physical maps were presented, which were generally in accordance. Information on the genetic and physical positions and allele sizes (when possible) of 17 DArT, 50 e-SSR, 44 SRAP, five ISSR, and two morphological markers is reported here for the first time. Seven segregation distortion regions (SDR) were identified on chromosomes 1B, 3BL, 4AL, 6AS, 6AL, 6BL, and 7B. A total of 22 and 12 QTLs for YD and yield difference between the value (YDDV) under HN and the value under LN were identified, respectively. Of these, QYd-4B-2 and QYddv-4B, two major stable QTL, shared support interval with alleles from KN9204 increasing YD in LN and decreasing YDDV. We probe into the use of these QTLs in wheat breeding programs. Moreover, factors affecting the SDR and total map length are discussed in depth.

Conclusions

This novel map may facilitate the use of novel markers in wheat molecular breeding programs and genomics research. Moreover, QTLs for YD and YDDV provide useful markers for wheat molecular breeding programs designed to increase yield potential under N stress.

Patterns of homoeologous gene expression shown by RNA sequencing in hexaploid bread wheat

BACKGROUND

Bread wheat (Triticum aestivum) has a large, complex and hexaploid genome consisting of A, B and D homoeologous chromosome sets. Therefore each wheat gene potentially exists as a trio of A, B and D homoeoloci, each of which may contribute differentially to wheat phenotypes. We describe a novel approach combining wheat cytogenetic resources (chromosome substitution 'nullisomic-tetrasomic' lines) with next generation deep sequencing of gene transcripts (RNA-Seq), to directly and accurately identify homoeologue-specific single nucleotide variants and quantify the relative contribution of individual homoeoloci to gene expression.

RESULTS

We discover, based on a sample comprising ~5-10% of the total wheat gene content, that at least 45% of wheat genes are expressed from all three distinct homoeoloci. Most of these genes show strikingly biased expression patterns in which expression is dominated by a single homoeolocus. The remaining ~55% of wheat genes are expressed from either one or two homoeoloci only, through a combination of extensive transcriptional silencing and homoeolocus loss.

CONCLUSIONS

We conclude that wheat is tending towards functional diploidy, through a variety of mechanisms causing single homoeoloci to become the predominant source of gene transcripts. This discovery has profound consequences for wheat breeding and our understanding of wheat evolution.

SNP Discovery for mapping alien introgressions in wheat

Background

Monitoring alien introgressions in crop plants is difficult due to the lack of genetic and molecular mapping information on the wild crop relatives. The tertiary gene pool of wheat is a very important source of genetic variability for wheat improvement against biotic and abiotic stresses. By exploring the 5M^g short arm (5M^gS) of Aegilops geniculata, we can apply chromosome genomics for the discovery of SNP markers and their use for monitoring alien introgressions in wheat (Triticum aestivum L).

Results

The short arm of chromosome 5M^g of Ae. geniculata Roth (syn. Ae. ovata L.; 2n = 4x = 28, U^gU^gM^gM^g) was flow-sorted from a wheat line in which it is maintained as a telocentric chromosome. DNA of the sorted arm was amplified and sequenced using an Illumina Hiseq 2000 with ~45x coverage. The sequence data was used for SNP discovery against wheat homoeologous group-5 assemblies. A total of 2,178 unique, 5M^gS-specific SNPs were discovered. Randomly selected samples of 59 5M^gS-specific SNPs were tested (44 by KASPar assay and 15 by Sanger sequencing) and 84% were validated. Of the selected SNPs, 97% mapped to a chromosome 5M^g addition to wheat (the source of t5M^gS), and 94% to 5M^g introgressed from a different accession of Ae. geniculata substituting for chromosome 5D of wheat. The validated SNPs also identified chromosome segments of 5M^gS origin in a set of T5D-5M^g translocation lines; eight SNPs (25%) mapped to TA5601 [T5DL · 5DS-5M^gS(0.75)] and three (8%) to TA5602 [T5DL · 5DS-5M^gS (0.95)]. SNPs (gsnp_5ms83 and gsnp_5ms94), tagging chromosome T5DL · 5DS-5M^gS(0.95) with the smallest introgression carrying resistance to leaf rust (Lr57) and stripe rust (Yr40), were validated in two released germplasm lines with Lr57 and Yr40 genes.

Conclusion

This approach should be widely applicable for the identification of species/genome-specific SNPs. The development of a large number of SNP markers will facilitate the precise introgression and monitoring of alien segments in crop breeding programs and further enable mapping and cloning novel genes from the wild relatives of crop plants.

Fine mapping and metabolic and physiological characterization of the glume glaucousness inhibitor locus Iw3 derived from wild wheat

This research provided the first view of metabolic and physiological effect of a tissue-specific glaucousness inhibitor in wheat and laid foundation for map-based cloning of the Iw3 locus. Cuticular wax constitutes the outermost layer of plant skin, and its composition greatly impacts plant appearance and plant-environment interaction. Epicuticular wax in the upper part of adult wheat plants can form the glaucousness, which is associated with drought tolerance. In this research, we characterized a glume-specific glaucousness inhibitor, Iw3, by fine mapping, physiological, and molecular approaches. Iw3 inhibits glaucousness formation by altering wax composition. Compared to the wild type, Iw3 eliminated β-diketone, reduced 47 % primary alcohols, but increased aldehyde 400-fold and alkanes fivefold, which led to 30 % reduction of total glume wax load. Loss of the glaucousness increased cuticle permeability, suggesting an important role in drought sensitivity. Genetically, the glaucousness-inhibiting effect by Iw3 is partially dominant in a dosage-dependent manner. We localized the Iw3 locus within a 0.13-cM interval delimited by marker loci Xpsp3000 and XWL3096. Of the 53 wax genes assayed, we detected transcription changes in nine genes by Iw3, downregulation of Cer4-1 and upregulation of other five Cer4 and three KCS homologs. All these results provided initial insights into Iw3-mediated regulation of wax metabolism and paved way for in-depth characterization of the Iw3 locus and the glaucousness-related β-diketone pathway.

Genetic rearrangements of six wheat-agropyron cristatum 6P addition lines revealed by molecular markers

Agropyron cristatum (L.) Gaertn. (2n = 4x = 28, PPPP) not only is cultivated as pasture fodder but also could provide many desirable genes for wheat improvement. It is critical to obtain common wheat–A. cristatum alien disomic addition lines to locate the desired genes on the P genome chromosomes. Comparative analysis of the homoeologous relationships between the P genome chromosome and wheat genome chromosomes is a key step in transferring different desirable genes into common wheat and producing the desired alien translocation line while compensating for the loss of wheat chromatin. In this study, six common wheat–A. cristatum disomic addition lines were produced and analyzed by phenotypic examination, genomic in situ hybridization (GISH), SSR markers from the ABD genomes and STS markers from the P genome. Comparative maps, six in total, were generated and demonstrated that all six addition lines belonged to homoeologous group 6. However, chromosome 6P had undergone obvious rearrangements in different addition lines compared with the wheat chromosome, indicating that to obtain a genetic compensating alien translocation line, one should recombine alien chromosomal regions with homoeologous wheat chromosomes. Indeed, these addition lines were classified into four types based on the comparative mapping: 6P_I, 6P_II, 6P_III, and 6P_IV. The different types of chromosome 6P possessed different desirable genes. For example, the 6P_I type, containing three addition lines, carried genes conferring high numbers of kernels per spike and resistance to powdery mildew, important traits for wheat improvement. These results may prove valuable for promoting the development of conventional chromosome engineering techniques toward molecular chromosome engineering.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

Gene expression analysis in the roots of salt-stressed wheat and the cytogenetic derivatives of wheat combined with the salt-tolerant wheatgrass, Lophopyrum elongatum

Using microarray analysis, we identified regulatory and signaling-related genes with differential expression in three genotypes with varying degrees of salt tolerance, Triticum aestivum , the amphiploid, and the wheat substitution line DS3E(3A). Lophopyrum elongatum is among one of the most salt-tolerant members of the Triticeae; important genetic stocks developed from crosses between wheat and L. elongatum provide a unique opportunity to compare gene expression in response to salt stress between these highly related species. The octaploid amphiploid contains the entire genome of T. aestivum and L. elongatum, and the disomic substitution line DS3E(3A) has chromosome 3A of wheat replaced by chromosome 3E of L. elongatum. In this study, microarray analysis was used to characterize gene expression profiles in the roots of three genotypes, Triticum aestivum, the octaploid amphiploid, and the wheat DS3E(3A) substitution line, in response to salt stress. We first examined changes in gene expression in wheat over a time course of 3 days of salt stress, and then compared changes in gene expression in wheat, the T. aestivum × L. elongatum amphiploid, and in the DS3E(3A) substitution line after 3 days of salt stress. In the time course experiment, 237 genes had 1.5 fold or greater change at least one out of three time points assayed in the experiment. The comparison between the three genotypes revealed 304 genes with significant differences in changes of expression between the genotypes. Forty-two of these genes had at least a twofold change in expression in response to salt treatment; 18 of these genes have signaling or regulatory function. Genes with significant differences in induction or repression between genotypes included transcription factors, protein kinases, ubiquitin ligases and genes related to phospholipid signaling.

A physical map of the short arm of wheat chromosome 1A

Bread wheat (Triticum aestivum) has a large and highly repetitive genome which poses major technical challenges for its study. To aid map-based cloning and future genome sequencing projects, we constructed a BAC-based physical map of the short arm of wheat chromosome 1A (1AS). From the assembly of 25,918 high information content (HICF) fingerprints from a 1AS-specific BAC library, 715 physical contigs were produced that cover almost 99% of the estimated size of the chromosome arm. The 3,414 BAC clones constituting the minimum tiling path were end-sequenced. Using a gene microarray containing ∼40 K NCBI UniGene EST clusters, PCR marker screening and BAC end sequences, we arranged 160 physical contigs (97 Mb or 35.3% of the chromosome arm) in a virtual order based on synteny with Brachypodium, rice and sorghum. BAC end sequences and information from microarray hybridisation was used to anchor 3.8 Mbp of Illumina sequences from flow-sorted chromosome 1AS to BAC contigs. Comparison of genetic and synteny-based physical maps indicated that ∼50% of all genetic recombination is confined to 14% of the physical length of the chromosome arm in the distal region. The 1AS physical map provides a framework for future genetic mapping projects as well as the basis for complete sequencing of chromosome arm 1AS.

Reticulate evolution of the rye genome

Rye (Secale cereale) is closely related to wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to its large genome (~8 Gb) and its regional importance, genome analysis of rye has lagged behind other cereals. Here, we established a virtual linear gene order model (genome zipper) comprising 22,426 or 72% of the detected set of 31,008 rye genes. This was achieved by high-throughput transcript mapping, chromosome survey sequencing, and integration of conserved synteny information of three sequenced model grass genomes (Brachypodium distachyon, rice [Oryza sativa], and sorghum [Sorghum bicolor]). This enabled a genome-wide high-density comparative analysis of rye/barley/model grass genome synteny. Seventeen conserved syntenic linkage blocks making up the rye and barley genomes were defined in comparison to model grass genomes. Six major translocations shaped the modern rye genome in comparison to a putative Triticeae ancestral genome. Strikingly dissimilar conserved syntenic gene content, gene sequence diversity signatures, and phylogenetic networks were found for individual rye syntenic blocks. This indicates that introgressive hybridizations (diploid or polyploidy hybrid speciation) and/or a series of whole-genome or chromosome duplications played a role in rye speciation and genome evolution.

Evaluation of genome sequencing quality in selected plant species using expressed sequence tags

BACKGROUND

With the completion of genome sequencing projects for more than 30 plant species, large volumes of genome sequences have been produced and stored in online databases. Advancements in sequencing technologies have reduced the cost and time of whole genome sequencing enabling more and more plants to be subjected to genome sequencing. Despite this, genome sequence qualities of multiple plants have not been evaluated.

METHODOLOGY/PRINCIPAL FINDING

Integrity and accuracy were calculated to evaluate the genome sequence quality of 32 plants. The integrity of a genome sequence is presented by the ratio of chromosome size and genome size (or between scaffold size and genome size), which ranged from 55.31% to nearly 100%. The accuracy of genome sequence was presented by the ratio between matched EST and selected ESTs where 52.93% ∼ 98.28% and 89.02% ∼ 98.85% of the randomly selected clean ESTs could be mapped to chromosome and scaffold sequences, respectively. According to the integrity, accuracy and other analysis of each plant species, thirteen plant species were divided into four levels. Arabidopsis thaliana, Oryza sativa and Zea mays had the highest quality, followed by Brachypodium distachyon, Populus trichocarpa, Vitis vinifera and Glycine max, Sorghum bicolor, Solanum lycopersicum and Fragaria vesca, and Lotus japonicus, Medicago truncatula and Malus × domestica in that order. Assembling the scaffold sequences into chromosome sequences should be the primary task for the remaining nineteen species. Low GC content and repeat DNA influences genome sequence assembly.

CONCLUSION

The quality of plant genome sequences was found to be lower than envisaged and thus the rapid development of genome sequencing projects as well as research on bioinformatics tools and the algorithms of genome sequence assembly should provide increased processing and correction of genome sequences that have already been published.

Chromosomal location of wheat genes of the carotenoid biosynthetic pathway and evidence for a catalase gene on chromosome 7A functionally associated with flour b* colour variation

Knowledge of molecular and genetic mechanisms controlling wheat grain quality characteristics is significant for improving flour for end-product functionality. Flour b* colour is an important quality trait for breeding wheat varieties to produce grain for specific market requirements. The degree of flour yellowness is due to the accumulation of carotenoids in grain, particularly lutein. Flour b* is under polygenic control and quantitative trait loci (QTL) have frequently been reported on chromosome 7AL. Analysis of carotenoid genes showed that phytoene synthase (PSY) co-located to the QTL on 7AL but other genes at this locus are also thought to contribute flour b* colour variation. This study used the wheat genome survey sequence and identified the chromosomal location of all wheat carotenoid genes, but none other than PSY were located on 7AL and, therefore, other genes may control flour b* colour variation including oxidative genes that degrade carotenoids. An investigation of EST bin mapped to 7AL identified a gene encoding a catalase enzyme (Cat3-A1) that was phylogenetically related to other plant class III enzymes, co-located to the QTL for flour b* colour variation on 7AL in three mapping populations and expressed during seed development. Therefore, Cat3-A1 was functionally associated with flour b* colour variation. Catalase acts upon hydrogen peroxide as a substrate and it was postulated that Cat3-A1 alleles control varying degrees of bleaching action on lutein in developing wheat grain. Markers for Cat3-A1 developed in this study can be used in conjunction with other candidate gene markers including phytoene synthase and lycopene-ε-cylase to develop a molecular signature for selecting lines with specific flour b* colour values in wheat breeding.

A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

Genetic association of OPR genes with resistance to Hessian fly in hexaploid wheat

BACKGROUND

Hessian fly (Mayetiola destructor) is one of the most destructive pests of wheat. The genes encoding 12-oxo-phytodienoic acid reductase (OPR) and lipoxygenase (LOX) play critical roles in insect resistance pathways in higher plants, but little is known about genes controlling resistance to Hessian fly in wheat.

RESULTS

In this study, 154 F6:8 recombinant inbred lines (RILs) generated from a cross between two cultivars, 'Jagger' and '2174' of hexaploid wheat (2n = 6 × =42; AABBDD), were used to map genes associated with resistance to Hessian fly. Two QTLs were identified. The first one was a major QTL on chromosome 1A (QHf.osu-1A), which explained 70% of the total phenotypic variation. The resistant allele at this locus in cultivar 2174 could be orthologous to one or more of the previously mapped resistance genes (H9, H10, H11, H16, and H17) in tetraploid wheat. The second QTL was a minor QTL on chromosome 2A (QHf.osu-2A), which accounted for 18% of the total phenotypic variation. The resistant allele at this locus in 2174 is collinear to an Yr17-containing-fragment translocated from chromosome 2N of Triticum ventricosum (2n = 4 × =28; DDNN) in Jagger. Genetic mapping results showed that two OPR genes, TaOPR1-A and TaOPR2-A, were tightly associated with QHf.osu-1A and QHf.osu-2A, respectively. Another OPR gene and three LOX genes were mapped but not associated with Hessian fly resistance in the segregating population.

CONCLUSIONS

This study has located two major QTLs/genes in bread wheat that can be directly used in wheat breeding programs and has also provided insights for the genetic association and disassociation of Hessian fly resistance with OPR and LOX genes in wheat.

Systems genetics of environmental response in the mature wheat embryo

Quantitative phenotypic traits are influenced by genetic and environmental variables as well as the interaction between the two. Underlying genetic × environment interaction is the influence that the surrounding environment exerts on gene expression. Perturbation of gene expression by environmental factors manifests itself in alterations to gene co-expression networks and ultimately in phenotypic plasticity. Comparative gene co-expression networks have been used to uncover biological mechanisms that differentiate tissues or other biological factors. In this study, we have extended consensus and differential Weighted Gene Co-Expression Network Analysis to compare the influence of different growing environments on gene co-expression in the mature wheat (Triticum aestivum) embryo. This network approach was combined with mapping of individual gene expression QTL to examine the genetic control of environmentally static and variable gene expression. The approach is useful for gene expression experiments containing multiple environments and allowed for the identification of specific gene co-expression modules responsive to environmental factors. This procedure identified conserved coregulation of gene expression between environments related to basic developmental and cellular functions, including protein localization and catabolism, vesicle composition/trafficking, Golgi transport, and polysaccharide metabolism among others. Environmentally unique modules were found to contain genes with predicted functions in responding to abiotic and biotic environmental variables. These findings represent the first report using consensus and differential Weighted Gene Co-expression Network Analysis to characterize the influence of environment on coordinated transcriptional regulation.

A non-additive interaction in a single locus causes a very short root phenotype in wheat

Non-additive allelic interactions underlie over dominant and under dominant inheritance, which explain positive and negative heterosis. These heteroses are often observed in the aboveground traits, but rarely reported in root. We identified a very short root (VSR) phenotype in the F1 hybrid between the common wheat (Triticum aestivum L.) landrace Chinese Spring and synthetic wheat accession TA4152-71. When germinated in tap water, primary roots of the parental lines reached ~15 cm 10 days after germination, but those of the F1 hybrid were ~3 cm long. Selfing populations segregated at a 1 (long-root) to 1 (short-root) ratio, indicating that VSR is controlled by a non-additive interaction between two alleles in a single gene locus, designated as Vsr1. Genome mapping localized the Vsr1 locus in a 3.8-cM interval delimited by markers XWL954 and XWL2506 on chromosome arm 5DL. When planted in vermiculite with supplemental fertilizer, the F1 hybrid had normal root growth, virtually identical to the parental lines, but the advanced backcrossing populations segregated for VSR, indicating that the F1 VSR expression was suppressed by interactions between other genes in the parental background and the vermiculite conditions. Preliminary physiological analyses showed that the VSR suppression is independent of light status but related to potassium homeostasis. Phenotyping additional hybrids between common wheat and synthetics revealed a high VSR frequency and their segregation data suggested more Vsr loci involved. Because the VSR plants can be regularly maintained and readily phenotyped at the early developmental stage, it provides a model for studies of non-additive interactions in wheat.

Physical mapping integrated with syntenic analysis to characterize the gene space of the long arm of wheat chromosome 1A

Background

Bread wheat (Triticum aestivum L.) is one of the most important crops worldwide and its production faces pressing challenges, the solution of which demands genome information. However, the large, highly repetitive hexaploid wheat genome has been considered intractable to standard sequencing approaches. Therefore the International Wheat Genome Sequencing Consortium (IWGSC) proposes to map and sequence the genome on a chromosome-by-chromosome basis.

Methodology/Principal Findings

We have constructed a physical map of the long arm of bread wheat chromosome 1A using chromosome-specific BAC libraries by High Information Content Fingerprinting (HICF). Two alternative methods (FPC and LTC) were used to assemble the fingerprints into a high-resolution physical map of the chromosome arm. A total of 365 molecular markers were added to the map, in addition to 1122 putative unique transcripts that were identified by microarray hybridization. The final map consists of 1180 FPC-based or 583 LTC-based contigs.

Conclusions/Significance

The physical map presented here marks an important step forward in mapping of hexaploid bread wheat. The map is orders of magnitude more detailed than previously available maps of this chromosome, and the assignment of over a thousand putative expressed gene sequences to specific map locations will greatly assist future functional studies. This map will be an essential tool for future sequencing of and positional cloning within chromosome 1A.

PmX: a recessive powdery mildew resistance gene at the Pm4 locus identified in wheat landrace Xiaohongpi

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most devastating foliar diseases of wheat and imposes a constant challenge on wheat breeders. Xiaohongpi, a Chinese landrace of wheat (Triticum aestivum L.), shows resistance to powdery mildew during the entire growth stage in the field and under controlled conditions. The F1 plants from cross of the powdery mildew susceptible cultivar Yangmai158 with Xiaohongpi were susceptible to isolate Bgt19, the locally most prevalent Bgt isolate. In the derived F2 population and F3 progenies, the resistance segregation deviated significantly from the one-gene Mendelian ratio. However, marker analysis indicated that only one recessive gene conferred the resistance, which co-segregated with Xsts-bcd1231 that showed co-segregation with Pm4a in different studies. Allelism test indicated that this recessive resistance gene, designated as pmX, is either allelic or tightly linked to Pm4a. The pmX gene was different from Pm4 alleles in resistance spectrum. Examination of the genotype frequencies at pmX and the linked marker loci in the F2 population showed that a genetic variation favoring the transmission of Xiaohongpi alleles could be the cause of deviated segregation. Mapping of the pmX-linked markers using Chinese Spring deletion lines indicated that it resides in the 0.85-1.00 bin of chromosome 2AL.

Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation

About 8,000 years ago in the Fertile Crescent, a spontaneous hybridization of the wild diploid grass Aegilops tauschii (2n = 14; DD) with the cultivated tetraploid wheat Triticum turgidum (2n = 4x = 28; AABB) resulted in hexaploid wheat (T. aestivum; 2n = 6x = 42; AABBDD). Wheat has since become a primary staple crop worldwide as a result of its enhanced adaptability to a wide range of climates and improved grain quality for the production of baker's flour. Here we describe sequencing the Ae. tauschii genome and obtaining a roughly 90-fold depth of short reads from libraries with various insert sizes, to gain a better understanding of this genetically complex plant. The assembled scaffolds represented 83.4% of the genome, of which 65.9% comprised transposable elements. We generated comprehensive RNA-Seq data and used it to identify 43,150 protein-coding genes, of which 30,697 (71.1%) were uniquely anchored to chromosomes with an integrated high-density genetic map. Whole-genome analysis revealed gene family expansion in Ae. tauschii of agronomically relevant gene families that were associated with disease resistance, abiotic stress tolerance and grain quality. This draft genome sequence provides insight into the environmental adaptation of bread wheat and can aid in defining the large and complicated genomes of wheat species.

Draft genome of the wheat A-genome progenitor Triticum urartu

Bread wheat (Triticum aestivum, AABBDD) is one of the most widely cultivated and consumed food crops in the world. However, the complex polyploid nature of its genome makes genetic and functional analyses extremely challenging. The A genome, as a basic genome of bread wheat and other polyploid wheats, for example, T. turgidum (AABB), T. timopheevii (AAGG) and T. zhukovskyi (AAGGA(m)A(m)), is central to wheat evolution, domestication and genetic improvement. The progenitor species of the A genome is the diploid wild einkorn wheat T. urartu, which resembles cultivated wheat more extensively than do Aegilops speltoides (the ancestor of the B genome) and Ae. tauschii (the donor of the D genome), especially in the morphology and development of spike and seed. Here we present the generation, assembly and analysis of a whole-genome shotgun draft sequence of the T. urartu genome. We identified protein-coding gene models, performed genome structure analyses and assessed its utility for analysing agronomically important genes and for developing molecular markers. Our T. urartu genome assembly provides a diploid reference for analysis of polyploid wheat genomes and is a valuable resource for the genetic improvement of wheat.

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

Fine mapping of wheat stripe rust resistance gene Yr26 based on collinearity of wheat with Brachypodium distachyon and rice

The Yr26 gene, conferring resistance to all currently important races of Puccinia striiformis f. sp. tritici (Pst) in China, was previously mapped to wheat chromosome deletion bin C-1BL-6-0.32 with low-density markers. In this study, collinearity of wheat to Brachypodium distachyon and rice was used to develop markers to saturate the chromosomal region containing the Yr26 locus, and a total of 2,341 F2 plants and 551 F2∶3 progenies derived from Avocet S×92R137 were used to develop a fine map of Yr26. Wheat expressed sequence tags (ESTs) located in deletion bin C-1BL-6-0.32 were used to develop sequence tagged site (STS) markers. The EST-STS markers flanking Yr26 were used to identify collinear regions of the rice and B. distachyon genomes. Wheat ESTs with significant similarities in the two collinear regions were selected to develop conserved markers for fine mapping of Yr26. Thirty-one markers were mapped to the Yr26 region, and six of them cosegregated with the resistance gene. Marker orders were highly conserved between rice and B. distachyon, but some rearrangements were observed between rice and wheat. Two flanking markers (CON-4 and CON-12) further narrowed the genomic region containing Yr26 to a 1.92 Mb region in B. distachyon chromosome 3 and a 1.17 Mb region in rice chromosome 10, and two putative resistance gene analogs were identified in the collinear region of B. distachyon. The markers developed in this study provide a potential target site for further map-based cloning of Yr26 and should be useful in marker assisted selection for pyramiding the gene with other resistance genes.

Wheat centromeric retrotransposons: the new ones take a major role in centromeric structure

The physical map of the hexaploid wheat chromosome 3B was screened using centromeric DNA probes. A 1.1-Mb region showing the highest number of positive bacterial artificial chromosome (BAC) clones was fully sequenced and annotated, revealing that 96% of the DNA consisted of transposable elements, mainly long terminal repeat (LTR) retrotransposons (88%). Estimation of the insertion time of the transposable elements revealed that CRW (also called Cereba) and Quinta are the youngest elements at the centromeres of common wheat (Triticum spp.) and its diploid ancestors, with Quinta being younger than CRW in both diploid and hexaploid wheats. Chromatin immunoprecipitation experiments revealed that both CRW and Quinta families are targeted by the centromere-specific histone H3 variant CENH3. Immuno colocalization of retroelements and CENH3 antibody indicated that a higher proportion of Quinta than CRWs was associated with CENH3, although CRWs were more abundant. Long arrays of satellite repeats were also identified in the wheat centromere regions, but they lost the ability to bind with CENH3. In addition to transposons, two functional genes and one pseudogene were identified. The gene density in the centromere appeared to be between three and four times lower than the average gene density of chromosome 3B. Comparisons with related grasses also indicated a loss of microcollinearity in this region. Finally, comparison of centromeric sequences of Aegilops tauschii (DD), Triticum boeoticum (AA) and hexaploid wheat revealed that the centromeres in both the polyploids and diploids are still undergoing dynamic changes, and that the new CRWs and Quintas may have undertaken the core role in kinetochore formation.

Fine mapping of leaf rust resistance gene LrZH84 using expressed sequence tag and sequence-tagged site markers, and allelism with other genes on wheat chromosome 1B

Zhou 8425B, possessing the leaf rust resistance gene LrZH84, is an elite wheat (Triticum aestivum) parental line in the Yellow-Huai Valley region of China. In the present study, 2,086 F(2) plants derived from Zhou 8425B/Chinese Spring were used for fine mapping of LrZH84 with expressed sequence tag (EST) and sequence-tagged site (STS) markers. Seventy inter-simple sequence repeat EST and STS markers on 1BL were used to screen the two parents and resistant and susceptible bulks; those polymorphic were used to analyze the entire F(2) population. Three EST markers (BF474863, BE497107, and CD373538) were closely linked to LrZH84, with genetic distances of 0.7, 0.7, and 1.7 cM, respectively. STS marker Hbsf-1 was developed from the sequences of polymerase chain reaction fragments amplified from EST marker BF474863. LrZH84 was 8.19 cM proximal to Lr44, but may be allelic to LrXi and LrG98 although they showed different reactions with some Puccinia triticina pathotypes.

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.

Wheat CBF gene family: identification of polymorphisms in the CBF coding sequence

Expression of cold-regulated genes needed for protection against freezing stress is mediated, in part, by the CBF transcription factor family. Previous studies with temperate cereals suggested that the CBF gene family in wheat was large, and that CBF genes were at the base of an important low temperature tolerance trait. Therefore, the goal of our study was to identify the CBF repertoire in the freezing-tolerant hexaploid wheat cultivar Norstar, and then to examine if the coding region of CBF genes in two spring cultivars contain polymorphisms that could affect the protein sequence and structure. Our analyses reveal that hexaploid wheat contains a complex CBF family consisting of at least 65 CBF genes of which 60 are known to be expressed in the cultivar Norstar. They represent 27 paralogous genes with 1-3 homeologous copies for the A, B, and D genomes. The cultivar Norstar contains two pseudogenes and at least 24 additional proteins having sequences and (or) structures that deviate from the consensus in the conserved AP2 DNA-binding and (or) C-terminal activation-domains. This suggests that in cultivars such as Norstar, low temperature tolerance may be increased through breeding of additional optimal alleles. The examination of the CBF repertoire present in the two spring cultivars, Chinese Spring and Manitou, reveals that they have additional polymorphisms affecting conserved positions in these domains. Understanding the effects of these polymorphisms will provide additional information for the selection of optimum CBF alleles in Triticeae breeding programs.

Analysis of the bread wheat genome using whole-genome shotgun sequencing

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

Discovery of high-confidence single nucleotide polymorphisms from large-scale de novo analysis of leaf transcripts of Aegilops tauschii, a wild wheat progenitor

Construction of high-resolution genetic maps is important for genetic and genomic research, as well as for molecular breeding. Single nucleotide polymorphisms (SNPs) are the predominant class of genetic variation and can be used as molecular markers. Aegilops tauschii, the D-genome donor of common wheat, is considered a valuable genetic resource for wheat improvement. Our previous study implied that Ae. tauschii accessions can be genealogically divided into two major lineages. In this study, the transcriptome of two Ae. tauschii accessions from each lineage, lineage 1 (L1) and 2 (L2), was sequenced, yielding 9435 SNPs and 739 insertion/deletion polymorphisms (indels) after de novo assembly of the reads. Based on 36 contig sequences, 31 SNPs and six indels were validated on 20 diverse Ae. tauschii accessions. Because almost all of the SNP markers were polymorphic between L1 and L2, and the D-genome donor of common wheat is presumed to belong to L2, these markers are available for D-genome typing in crosses between common wheat varieties and L1-derived synthetic wheat. Due to the conserved synteny between wheat and barley chromosomes, the high-density expressed sequence tag barley map and the hypothetical gene order in barley can be applied to develop markers on target chromosomal regions in wheat.

Clusters of genes encoding fructan biosynthesizing enzymes in wheat and barley

Fructans are soluble carbohydrates with health benefits and possible roles in plant adaptation. Fructan biosynthetic genes were isolated using comparative genomics and physical mapping followed by BAC sequencing in barley. Genes encoding sucrose:sucrose 1-fructosyltransferase (1-SST), fructan:fructan 1-fructosyltransferase (1-FFT) and sucrose:fructan 6-fructosyltransferase (6-SFT) were clustered together with multiple copies of vacuolar invertase genes and a transposable element on two barley BAC. Intron-exon structures of the genes were similar. Phylogenetic analysis of the fructosyltransferases and invertases in the Poaceae showed that the fructan biosynthetic genes may have evolved from vacuolar invertases. Quantitative real-time PCR was performed using leaf RNA extracted from three wheat cultivars grown under different conditions. The 1-SST, 1-FFT and 6-SFT genes had correlated expression patterns in our wheat experiment and in existing barley transcriptome database. Single nucleotide polymorphism (SNP) markers were developed and successfully mapped to a major QTL region affecting wheat grain fructan accumulation in two independent wheat populations. The alleles controlling high- and low- fructan in parental lines were also found to be associated in fructan production in a diverse set of 128 wheat lines. To the authors' knowledge, this is the first report on the mapping and sequencing of a fructan biosynthetic gene cluster and in particular, the isolation of a novel 1-FFT gene from barley.

Expressed sequence tag-PCR markers for identification of alien barley chromosome 2H in wheat

We developed EST-PCR markers specific to barley chromosome 2H, for the purpose of effectively tracing alien chromosomes or chromosome parts in the wheat genetic background. The target alien chromosome 2H confers high resistance to pre-harvest sprouting, which is a worldwide natural disaster in wheat. A total of 120 primer pairs were selected by combining the wheat group 2 chromosomes of the EST database and the genome sequences of the new model plant Brachypodium distachyon. Seventy-seven of 120 primer pairs were polymorphic and 31 of 120 primer pairs were monomorphic between a set of wheat-barley chromosome 2H disomic addition/substitution lines and their parents by agarose gel electrophoresis and polyacrylamide gel electrophoresis. Thirty of 77 polymorphic primer pairs including primer pair P120 derived from the basi gene were chromosome 2H-specific. These markers are expected to be valuable in screening of wheat-barley chromosome 2H recombination lines and pre-harvest sprouting resistant varieties.

Genetic analysis and molecular mapping of a new powdery mildew resistant gene Pm46 in common wheat

Powdery mildew (PM), caused by Blumeria graminis f. sp. tritici (Bgt), has become a serious disease and caused severe yield losses in the wheat production worldwide. Resistance gene(s) in wheat cultivars can be quickly overcome by newly evolved pathogen races when these genes are employed for long time or in a large area. It is urgent to search for new sources of resistance to be used in wheat breeding. Tabasco is a German resistant cultivar and a new source of resistance gene(s) to PM. An F(2) population was developed from a cross between Tabasco and a Chinese susceptible cultivar Ningnuo 1. Infection types in 472 F(2) plants and 436 F(2-3) families were evaluated by inoculating plants with isolate Bgt19. Results showed that a single dominant gene, designed Pm46, controlled powdery mildew resistance in Tabasco. This gene was located to the short arm of chromosome 5D (5DS) and flanked by simple sequence repeat markers Xgwm205 and Xcfd81 at 18.9 cM apart. Because another resistance gene Pm2 was also located on 5DS, 15 Bgt isolates were used to inoculate Tabasco and Ulka/8*Cc (Pm2 carrier). The results showed that Tabasco was highly resistant to all of the 15 isolates tested, while Ulka/8*Cc was susceptible to 4 of the isolates, suggesting that Tabasco may carry resistant gene(s) different from Pm2 gene in Ulka/8*Cc. To test the allelism between Pm46 and Pm2, an F(2) population between Tabasco and Ulka/8*Cc was developed. Isolate Bgt2, avirulent to both parents, was used to evaluate the F(2) population and two susceptible plants were identified from 536 progenies with F(2) plants. This result indicated that Pm46 is not allelic to Pm2. Therefore, Pm46 is a new gene for PM resistance identified in this study.

Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.

Construction of whole genome radiation hybrid panels and map of chromosome 5A of wheat using asymmetric somatic hybridization

To explore the feasibility of constructing a whole genome radiation hybrid (WGRH) map in plant species with large genomes, asymmetric somatic hybridization between wheat (Triticum aestivum L.) and Bupleurum scorzonerifolium Willd. was performed. The protoplasts of wheat were irradiated with ultraviolet light (UV) and gamma-ray and rescued by protoplast fusion using B. scorzonerifolium as the recipient. Assessment of SSR markers showed that the radiation hybrids have the average marker retention frequency of 15.5%. Two RH panels (RHPWI and RHPWII) that contained 92 and 184 radiation hybrids, respectively, were developed and used for mapping of 68 SSR markers in chromosome 5A of wheat. A total of 1557 and 2034 breaks were detected in each panel. The RH map of chromosome 5A based on RHPWII was constructed. The distance of the comprehensive map was 2103 cR and the approximate resolution was estimated to be ∼501.6 kb/break. The RH panels evaluated in this study enabled us to order the ESTs in a single deletion bin or in the multiple bins cross the chromosome. These results demonstrated that RH mapping via protoplast fusion is feasible at the whole genome level for mapping purposes in wheat and the potential value of this mapping approach for the plant species with large genomes.

Genome-wide SNP identification in multiple morphotypes of allohexaploid tall fescue (Festuca arundinacea Schreb)

BACKGROUND

Single nucleotide polymorphisms (SNPs) provide essential tools for the advancement of research in plant genomics, and the development of SNP resources for many species has been accelerated by the capabilities of second-generation sequencing technologies. The current study aimed to develop and use a novel bioinformatic pipeline to generate a comprehensive collection of SNP markers within the agriculturally important pasture grass tall fescue; an outbreeding allopolyploid species displaying three distinct morphotypes: Continental, Mediterranean and rhizomatous.

RESULTS

A bioinformatic pipeline was developed that successfully identified SNPs within genotypes from distinct tall fescue morphotypes, following the sequencing of 414 polymerase chain reaction (PCR) - generated amplicons using 454 GS FLX technology. Equivalent amplicon sets were derived from representative genotypes of each morphotype, including six Continental, five Mediterranean and one rhizomatous. A total of 8,584 and 2,292 SNPs were identified with high confidence within the Continental and Mediterranean morphotypes respectively. The success of the bioinformatic approach was demonstrated through validation (at a rate of 70%) of a subset of 141 SNPs using both SNaPshot™ and GoldenGate™ assay chemistries. Furthermore, the quantitative genotyping capability of the GoldenGate™ assay revealed that approximately 30% of the putative SNPs were accessible to co-dominant scoring, despite the hexaploid genome structure. The sub-genome-specific origin of each SNP validated from Continental tall fescue was predicted using a phylogenetic approach based on comparison with orthologous sequences from predicted progenitor species.

CONCLUSIONS

Using the appropriate bioinformatic approach, amplicon resequencing based on 454 GS FLX technology is an effective method for the identification of polymorphic SNPs within the genomes of Continental and Mediterranean tall fescue. The GoldenGate™ assay is capable of high-throughput co-dominant SNP allele detection, and minimises the problems associated with SNP genotyping in a polyploid by effectively reducing the complexity to a diploid system. This SNP collection may now be refined and used in applications such as cultivar identification, genetic linkage map construction, genome-wide association studies and genomic selection in tall fescue. The bioinformatic pipeline described here represents an effective general method for SNP discovery within outbreeding allopolyploid species.

High-density mapping and marker development for the powdery mildew resistance gene PmAS846 derived from wild emmer wheat (Triticum turgidum var. dicoccoides)

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is an important foliar disease of wheat worldwide. The dominant powdery mildew resistance gene PmAS846 was transferred to the hexaploid wheat lines N9134 and N9738 from wild emmer wheat (Triticum dicoccoides) in 1995, and it is still one of the most effective resistance genes in China. A high resolution genetic map for PmAS846 locus was constructed using two F(2) populations and corresponding F(2:3) families developed from the crosses of N9134/Shaanyou 225 and N9738/Huixianhong. Synteny between wheat and Brachypodium distachyon and rice was used to develop closely linked molecular markers to reduce the genetic interval around PmAS846. Twenty-six expressed sequence tag-derived markers were mapped to the PmAS846 locus. Five markers co-segregated with PmAS846 in the F(2) population of N9134/Shaanyou 225. PmAS846 was physically located to wheat chromosome 5BL bin 0.75-0.76 within a gene-rich region. The markers order is conserved between wheat and Brachypodium distachyon, but rearrangements are present in rice. Two markers, BJ261635 and CJ840011 flanked PmAS846 and narrowed PmAS846 to a region that is collinear with 197 and 112 kb genomic regions on Brachypodium chromosome 4 and rice chromosome 9, respectively. The genes located on the corresponding homologous regions in Brachypodium, rice and barley could be considered for further marker saturation and identification of potential candidate genes for PmAS846. The markers co-segregating with PmAS846 provide a potential target site for positional cloning of PmAS846, and can be used for marker-assisted selection of this gene.

Identification and validation of a major quantitative trait locus for slow-rusting resistance to stripe rust in wheat

Stripe (yellow) rust, caused by Puccinia striiformis Westend. f. sp. tritici Eriks (Pst), is one of the most important wheat (Triticum aestivum L.) diseases and causes significant yield losses. A recombinant inbred (RI) population derived from a cross between Yanzhan 1 and Xichang 76-9 cultivars was evaluated for resistance to wheat stripe rust strain CYR32 at both the seedling and adult plant stages. Four resistance quantitative trait loci (QTLs) were detected in this population, in which the major one, designated as Yrq1, was mapped on chromosome 2DS. The strategy of using the Brachypodium distachyon genome, wheat expressed sequence tags and a draft DNA sequences (scaffolds) of the D-genome (Aegilops tauschii Coss.) for the development of simple sequence repeat (SSR) markers was successfully used to identify 147 SSRs in hexaploid wheat. Of the 19 polymorphic SSRs in the RI population, 17 SSRs were mapped in the homeologous group 2 chromosomes near Yrq1 region and eight SSRs were genetically mapped in the 2.7 cM region of Yrq1, providing abundant DNA markers for fine-mapping of Yrq1 and marker-assisted selection in wheat breeding program. The effectiveness of Yrq1 was validated in an independent population, indicating that this resistance QTL can be successfully transferred into a susceptible cultivar for improvement of stripe rust resistance.