Development of an expressed sequence tag (EST) resource for wheat (Triticum aestivum L.): EST generation, unigene analysis, probe selection and bioinformatics for a 16,000-locus bin-delineated map

Structural and Functional Characterization at the Molecular Level of the MATE Gene Family in Wheat in Silico

A series of multidrug extransporters known as the multidrug and potentially toxic extrusion (MATE) genes are found in all living things and are crucial for the removal of heavy metal ions, metalloids, exogenous xenobiotics, endogenous secondary metabolites, and other toxic substances from the cells. However, there has only been a small amount of them in silico analysis of the MATE family of genes in plant species. In the current study, the MATE gene family was characterized in silico where two families and seven subfamilies based on their evolutionary relationships were proposed. Plant breeders may use TraesCS1D02G030400, TraesCS4B02G244400, and TraesCS1A02G029900 genes for marker-assisted or transgenic breeding to develop novel cultivars since these genes have been hypothesized from protein-protein interaction study to play a critical role in the transport of toxic chemicals across cells. The exon number varies from 01 to 14. One exon has TraesCS1A02G188100, TraesCS5B02G562500, TraesCS6A02G256400, and TraesCS6D02G384300 genes, while 14 exons have only two genes that are TraesCS6A02G418800 and TraesCS6D02G407900. Biological stress (infestations of disease) affects the expression of most of the MATE genes, with the gene TraesCS5D02G355500 having the highest expression level in the wheat expression browser tool. Using the Grain interpretation search engine tool, it is found that the vast bulk of MATE genes are voiced throughout biotic environmental stresses caused by disease pests, with the genotype TraesCS5B02G326600.1 from family 1 exhibiting the greatest level of expression throughout Fusarium head blight infection by Fusarium graminearum after 4 days of infection. The researchers constructed 39 ternary plots, each with a distinct degree of expression under biotic and abiotic stress settings, and observed that 44% of the triplets have imbalanced outputs (extreme values) due to their higher tissue specificity and increased intensity.

Identification of a Recessive Gene PmQ Conferring Resistance to Powdery Mildew in Wheat Landrace Qingxinmai Using BSR-Seq Analysis

Wheat powdery mildew is caused by Blumeria graminis f. sp. tritici (Bgt), a biotrophic fungal species. It is very important to mine new powdery mildew (Pm) resistance genes for developing resistant wheat cultivars to reduce the deleterious effects of the disease. This study was carried out to characterize the Pm gene in Qingxinmai, a winter wheat landrace from Xinjiang, China. Qingxinmai is resistant to many Bgt isolates collected from different wheat fields in China. F₁, F₂, and F_2:3 generations of the cross between Qingxinmai and powdery mildew susceptible line 041133 were developed. It was confirmed that a single recessive gene, PmQ, conferred the seedling resistance to a Bgt isolate in Qingxinmai. Bulked segregant analysis-RNA-Seq (BSR-Seq) was performed on the bulked homozygous resistant and susceptible F_2:3 families, which detected 57 single nucleotide polymorphism (SNP) variants that were enriched in a 40 Mb genomic interval on chromosome arm 2BL. Based on the flanking sequences of the candidate SNPs extracted from the Chinese Spring reference genome, 485 simple sequence repeat (SSR) markers were designed. Six polymorphic SSR markers, together with nine markers that were anchored on chromosome arm 2BL, were used to construct a genetic linkage map for PmQ. This gene was placed in a 1.4 cM genetic interval between markers Xicsq405 and WGGBH913 corresponding to 4.9 Mb physical region in the Chinese Spring reference genome. PmQ differed from most of the other Pm genes identified on chromosome arm 2BL based on its position and/or origin. However, this gene and Pm63 from an Iranian common wheat landrace were located in a similar genomic region, so they may be allelic.

Fine mapping of the wheat powdery mildew resistance gene Pm52 using comparative genomics analysis and the Chinese Spring reference genomic sequence

A high-resolution genetic linkage map was constructed using the comparative genomics analysis approach and the wheat reference genome, which placed wheat powdery mildew resistance gene Pm52 in a 0.21-cM genetic interval on chromosome arm 2BL. The gene Pm52 confers resistance to powdery mildew and has been previously mapped on chromosome arm 2BL in winter wheat cultivar Liangxing 99. Because of its effectiveness against the disease, this study was initiated to finely map Pm52 using the comparative genomics analysis approach and the wheat reference genomic sequence. Based on the EST sequences that were located in the chromosome region flanking Pm52, four EST-SSR markers were developed, and another nine SSR markers were developed using the comparative genomics technology. These thirteen markers were integrated into a genetic linkage map using an F_2:3 subpopulation of the Liangxing 99 × Zhongzuo 9504 cross. Pm52 was mapped within a 3.2-cM genetic interval in the subpopulation that corresponded to a ~40-Mb genomic interval on chromosome arm 2BL of the Chinese Spring reference genome. The Pm52-flanking markers Xicsl163 and Xicsl62 identified 344 recombinant individuals from 8820 F₂ plants. Nine SSR markers generated from the Chinese Spring genomic interval were incorporated into a high-resolution genetic linkage map, which placed Pm52 in a 0.21-cM genetic interval corresponding to 5.6-Mb genomic region. The constructed high-resolution genetic linkage map will facilitate the map-based cloning of Pm52 and its marker-assisted selection.

Unraveling Molecular and Genetic Studies of Wheat (Triticum aestivum L.) Resistance against Factors Causing Pre-Harvest Sprouting

Pre-harvest sprouting (PHS) is one of the most important factors having adverse effects on yield and grain quality all over the world, particularly in wet harvest conditions. PHS is controlled by both genetic and environmental factors and the interaction of these factors. Breeding varieties with high PHS resistance have important implications for reducing yield loss and improving grain quality. The rapid advancements in the wheat genomic database along with transcriptomic and proteomic technologies have broadened our knowledge for understanding the regulatory mechanism of PHS resistance at transcriptomic and post-transcriptomic levels. In this review, we have described in detail the recent advancements on factors influencing PHS resistance, including grain color, seed dormancy, α-amylase activity, plant hormones (especially abscisic acid and gibberellin), and QTL/genes, which are useful for mining new PHS-resistant genes and developing new molecular markers for multi-gene pyramiding breeding of wheat PHS resistance, and understanding the complicated regulatory mechanism of PHS resistance.

A High-Density Consensus Map of Common Wheat Integrating Four Mapping Populations Scanned by the 90K SNP Array

A high-density consensus map is a powerful tool for gene mapping, cloning and molecular marker-assisted selection in wheat breeding. The objective of this study was to construct a high-density, single nucleotide polymorphism (SNP)-based consensus map of common wheat (Triticum aestivum L.) by integrating genetic maps from four recombinant inbred line populations. The populations were each genotyped using the wheat 90K Infinium iSelect SNP assay. A total of 29,692 SNP markers were mapped on 21 linkage groups corresponding to 21 hexaploid wheat chromosomes, covering 2,906.86 cM, with an overall marker density of 10.21 markers/cM. Compared with the previous maps based on the wheat 90K SNP chip detected 22,736 (76.6%) of the SNPs with consistent chromosomal locations, whereas 1,974 (6.7%) showed different chromosomal locations, and 4,982 (16.8%) were newly mapped. Alignment of the present consensus map and the wheat expressed sequence tags (ESTs) Chromosome Bin Map enabled assignment of 1,221 SNP markers to specific chromosome bins and 819 ESTs were integrated into the consensus map. The marker orders of the consensus map were validated based on physical positions on the wheat genome with Spearman rank correlation coefficients ranging from 0.69 (4D) to 0.97 (1A, 4B, 5B, and 6A), and were also confirmed by comparison with genetic position on the previously 40K SNP consensus map with Spearman rank correlation coefficients ranging from 0.84 (6D) to 0.99 (6A). Chromosomal rearrangements reported previously were confirmed in the present consensus map and new putative rearrangements were identified. In addition, an integrated consensus map was developed through the combination of five published maps with ours, containing 52,607 molecular markers. The consensus map described here provided a high-density SNP marker map and a reliable order of SNPs, representing a step forward in mapping and validation of chromosomal locations of SNPs on the wheat 90K array. Moreover, it can be used as a reference for quantitative trait loci (QTL) mapping to facilitate exploitation of genes and QTL in wheat breeding.

A whole-genome, radiation hybrid mapping resource of hexaploid wheat

Generating a contiguous, ordered reference sequence of a complex genome such as hexaploid wheat (2n = 6x = 42; approximately 17 GB) is a challenging task due to its large, highly repetitive, and allopolyploid genome. In wheat, ordering of whole-genome or hierarchical shotgun sequencing contigs is primarily based on recombination and comparative genomics-based approaches. However, comparative genomics approaches are limited to syntenic inference and recombination is suppressed within the pericentromeric regions of wheat chromosomes, thus, precise ordering of physical maps and sequenced contigs across the whole-genome using these approaches is nearly impossible. We developed a whole-genome radiation hybrid (WGRH) resource and tested it by genotyping a set of 115 randomly selected lines on a high-density single nucleotide polymorphism (SNP) array. At the whole-genome level, 26 299 SNP markers were mapped on the RH panel and provided an average mapping resolution of approximately 248 Kb/cR1500 with a total map length of 6866 cR1500 . The 7296 unique mapping bins provided a five- to eight-fold higher resolution than genetic maps used in similar studies. Most strikingly, the RH map had uniform bin resolution across the entire chromosome(s), including pericentromeric regions. Our research provides a valuable and low-cost resource for anchoring and ordering sequenced BAC and next generation sequencing (NGS) contigs. The WGRH developed for reference wheat line Chinese Spring (CS-WGRH), will be useful for anchoring and ordering sequenced BAC and NGS based contigs for assembling a high-quality, reference sequence of hexaploid wheat. Additionally, this study provides an excellent model for developing similar resources for other polyploid species.

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i.e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects.

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.

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.

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.

dbWFA: a web-based database for functional annotation of Triticum aestivum transcripts

The functional annotation of genes based on sequence homology with genes from model species genomes is time-consuming because it is necessary to mine several unrelated databases. The aim of the present work was to develop a functional annotation database for common wheat Triticum aestivum (L.). The database, named dbWFA, is based on the reference NCBI UniGene set, an expressed gene catalogue built by expressed sequence tag clustering, and on full-length coding sequences retrieved from the TriFLDB database. Information from good-quality heterogeneous sources, including annotations for model plant species Arabidopsis thaliana (L.) Heynh. and Oryza sativa L., was gathered and linked to T. aestivum sequences through BLAST-based homology searches. Even though the complexity of the transcriptome cannot yet be fully appreciated, we developed a tool to easily and promptly obtain information from multiple functional annotation systems (Gene Ontology, MapMan bin codes, MIPS Functional Categories, PlantCyc pathway reactions and TAIR gene families). The use of dbWFA is illustrated here with several query examples. We were able to assign a putative function to 45% of the UniGenes and 81% of the full-length coding sequences from TriFLDB. Moreover, comparison of the annotation of the whole T. aestivum UniGene set along with curated annotations of the two model species assessed the accuracy of the annotation provided by dbWFA. To further illustrate the use of dbWFA, genes specifically expressed during the early cell division or late storage polymer accumulation phases of T. aestivum grain development were identified using a clustering analysis and then annotated using dbWFA. The annotation of these two sets of genes was consistent with previous analyses of T. aestivum grain transcriptomes and proteomes.

Database URL:urgi.versailles.inra.fr/dbWFA/

The gene space in wheat: the complete γ-gliadin gene family from the wheat cultivar Chinese Spring

The complete set of unique γ-gliadin genes is described for the wheat cultivar Chinese Spring using a combination of expressed sequence tag (EST) and Roche 454 DNA sequences. Assemblies of Chinese Spring ESTs yielded 11 different γ-gliadin gene sequences. Two of the sequences encode identical polypeptides and are assumed to be the result of a recent gene duplication. One gene has a 3' coding mutation that changes the reading frame in the final eight codons. A second assembly of Chinese Spring γ-gliadin sequences was generated using Roche 454 total genomic DNA sequences. The 454 assembly confirmed the same 11 active genes as the EST assembly plus two pseudogenes not represented by ESTs. These 13 γ-gliadin sequences represent the complete unique set of γ-gliadin genes for cv Chinese Spring, although not ruled out are additional genes that are exact duplications of these 13 genes. A comparison with the ESTs of two other hexaploid cultivars (Butte 86 and Recital) finds that the most active genes are present in all three cultivars, with exceptions likely due to too few ESTs for detection in Butte 86 and Recital. A comparison of the numbers of ESTs per gene indicates differential levels of expression within the γ-gliadin gene family. Genome assignments were made for 6 of the 13 Chinese Spring γ-gliadin genes, i.e., one assignment from a match to two γ-gliadin genes found within a tetraploid wheat A genome BAC and four genes that match four distinct γ-gliadin sequences assembled from Roche 454 sequences from Aegilops tauschii, the hexaploid wheat D-genome ancestor.

Cloning and characterization of two Argonaute genes in wheat (Triticum aestivum L.)

BACKGROUND

Argonaute proteins are key components of RNA interference (RNAi), playing important roles in RNA-directed gene silencing. Various classes of Argonaute genes have been identified from plants and might be involved in developmental regulation. However, little is known about these genes in wheat (Triticum aestivum).

RESULTS

In this study, two full-length cDNAs of Argonaute were cloned from wheat, designated as TaAGO1b and TaAGO4. The cDNA of TaAGO1b is 3273 bp long and encodes 868 amino acids, with a predicted molecular weight of ~97.78 kDa and pI of 9.29. The 3157-bp TaAGO4 encodes 916 amino acids, with a molecular mass of 102.10 kDa and pI of 9.12. Genomics analysis showed that TaAGO1b and TaAGO4 contain 20 and 18 introns, respectively. Protein structural analysis demonstrated that typical PAZ and PIWI domains were found in both TaAGO1b and TaAGO4. From the highly conserved PIWI domains, we detected conserved Asp-Asp-His (DDH) motifs that function as a catalytic triad and have critical roles during the process of sequence-specific cleavage in the RNAi machinery. Structural modelling indicated that both TaAGOs can fold to a specific α/β structure. Moreover, the three aligned DDH residues are spatially close to each other at the "slicer" site of the PIWI domain. Expression analysis indicated that both genes are ubiquitously expressed in vegetative and reproductive organs, including the root, stem, leaf, anther, ovule, and seed. However, they are differentially expressed in germinating endosperm tissues. We were interested to learn that the two TaAGOs are also differentially expressed in developing wheat plants and that their expression patterns are variously affected by vernalization treatment. Further investigation revealed that they can be induced by cold accumulation during vernalization.

CONCLUSIONS

Two putative wheat Argonaute genes, TaAGO1b and TaAGO4, were cloned. Phylogenetic analysis, prediction of conserved domains and catalytic motifs, and modelling of their protein structures suggested that they encode functional Argonaute proteins. Temporal and spatial expression analyses indicated that these genes are potentially involved in developmental regulation of wheat plants.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.

Large-scale development of functional markers in Brassica species

Numerous quantitative trait loci (QTL) have been detected in Brassica species, but fine-mapping of major QTL has advanced slowly. The development of functional markers can overcome this barrier. We used publicly available PlantGDB-assembled unique transcripts (PUTs) from Brassica species to design 7836 functional simple sequence repeat (SSR) primer pairs. Functional annotation of the PUTs containing SSRs was done by Blast2GO. The PUTs harbouring SSRs were mainly involved with nucleotide or protein binding and enzyme activity, and preferentially functioned in membranes and cytoplasm. Totally, 210 PUT primer pairs were selected to test their polymorphism, stability, and PCR quality. Approximately 70% (147) of the primer pairs resulted in successful amplification with an average polymorphic information content (PIC) value of 0.49. The highest level of polymorphism was dinucleotide repeat SSRs, followed by tri- and mononucleotide repeats. Approximately 60% of the primer pairs showed good transferability among Brassica species. These results show that the development of markers from PUTs is a feasible and simple approach to develop functional SSR markers on a large scale across Brassica species. In addition, these markers can provide a novel alternative that is a putative approach for rapid determination of candidate genes, genetic mapping, genetic diversity analysis, and comparative mapping in Brassica species.

Identification of drought-responsive universal stress proteins in viridiplantae

Genes encoding proteins that contain the universal stress protein (USP) domain are known to provide bacteria, archaea, fungi, protozoa, and plants with the ability to respond to a plethora of environmental stresses. Specifically in plants, drought tolerance is a desirable phenotype. However, limited focused and organized functional genomic datasets exist on drought-responsive plant USP genes to facilitate their characterization. The overall objective of the investigation was to identify diverse plant universal stress proteins and Expressed Sequence Tags (ESTs) responsive to water-deficit stress. We hypothesize that cross-database mining of functional annotations in protein and gene transcript bioinformatics resources would help identify candidate drought-responsive universal stress proteins and transcripts from multiple plant species. Our bioinformatics approach retrieved, mined and integrated comprehensive functional annotation data on 511 protein and 1561 ESTs sequences from 161 viridiplantae taxa. A total of 32 drought-responsive ESTs from 7 plant genera Glycine, Hordeum, Manihot, Medicago, Oryza, Pinus and Triticum were identified. Two Arabidopsis USP genes At3g62550 and At3g53990 that encode ATP-binding motif were up-regulated in a drought microarray dataset. Further, a dataset of 80 simple sequence repeats (SSRs) linked to 20 singletons and 47 transcript assembles was constructed. Integrating the datasets on SSRs and drought-responsive ESTs identified three drought-responsive ESTs from bread wheat (BE604157), soybean (BM887317) and maritime pine (BX682209). The SSR sequence types were CAG, ATA and AT respectively. The datasets from cross-database mining provide organized resources for the characterization of USP genes as useful targets for engineering plant varieties tolerant to unfavorable environmental conditions.

Pathogenesis-associated transcriptional patterns in Triticeae

The Triticeae tribe of the plant Poaceae family contains some of the most important cereal crop plants for nutrition of humans and livestock such as wheat and barley. Despite the agronomical relevance of plant immunity, knowledge on mechanisms of disease or resistance in Triticeae is limited. It is hardly understood what actually stops a microbial invader when restricted by the plant and in how far a susceptible host plant contributes to pathogenesis. Transcriptional reprogramming of the host plant may be involved in both immunity and disease. This paper gives an overview about recent analyses of global pathogenesis-related transcriptional patterns in response of Triticeae to biotrophic or non-biotrophic fungal pathogens and their toxins. It highlights enriched biological functions in association with successful plant defence or disease as well as experiments that successfully translated gene expression data into analysis of gene functions.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

In silico identification of coffee genome expressed sequences potentially associated with resistance to diseases

Sequences potentially associated with coffee resistance to diseases were identified by in silico analyses using the database of the Brazilian Coffee Genome Project (BCGP). Keywords corresponding to plant resistance mechanisms to pathogens identified in the literature were used as baits for data mining. Expressed sequence tags (ESTs) related to each of these keywords were identified with tools available in the BCGP bioinformatics platform. A total of 11,300 ESTs were mined. These ESTs were clustered and formed 979 EST-contigs with similarities to chitinases, kinases, cytochrome P450 and nucleotide binding site-leucine rich repeat (NBS-LRR) proteins, as well as with proteins related to disease resistance, pathogenesis, hypersensitivity response (HR) and plant defense responses to diseases. The 140 EST-contigs identified through the keyword NBS-LRR were classified according to function. This classification allowed association of the predicted products of EST-contigs with biological processes, including host defense and apoptosis, and with molecular functions such as nucleotide binding and signal transducer activity. Fisher's exact test was used to examine the significance of differences in contig expression between libraries representing the responses to biotic stress challenges and other libraries from the BCGP. This analysis revealed seven contigs highly similar to catalase, chitinase, protein with a BURP domain and unknown proteins. The involvement of these coffee proteins in plant responses to disease is discussed.

Development of EST-PCR markers for Thinopyrum intermedium chromosome 2Ai#2 and their application in characterization of novel wheat-grass recombinants

A series of expressed sequence tags-derived polymerase chain reaction (EST-PCR) markers specific to chromosome 2Ai#2 from Thinopyrum intermedium were developed in this study using a new integrative approach. The target alien chromosome confers high resistance to barley yellow dwarf virus (BYDV), which is a severe virus disease in wheat. To generate markers evenly distributed on 2Ai#2, a total of 105 primer pairs were designed based on mapped ESTs from 8 bins of wheat chromosome 2B with intron-prediction by aligning ESTs with genomic sequences of the new model plant Brachypodium distachyon. Eight and seven polymorphic markers on the short arm and the long arm of chromosome 2Ai#2, respectively, were obtained with a polymorphism rate of 14.3%. These chromosome 2Ai#2-specific EST-PCR markers were then used in tracing and exploring the structural variation of the alien chromosome in the population derived from the immature embryo culture of the cross between N452, a 2Ai#2(2D) substitution line, and common wheat CB037. Two centric fusion of translocations involving 2Ai#2 short or long arm with wheat chromosome 2D and some new genetic stocks including telosomes with the alien chromosome short or long arm were identified in the SC(3) generations, which provided basic materials to further study the mechanism of the BYDV resistance. BYDV tests in two field seasons suggest that the BYDV resistance was mainly conferred by the short arm, gene interaction on both arms of the alien chromosome was discussed.

A BAC-based physical map of Brachypodium distachyon and its comparative analysis with rice and wheat

Background

Brachypodium distachyon (Brachypodium) has been recognized as a new model species for comparative and functional genomics of cereal and bioenergy crops because it possesses many biological attributes desirable in a model, such as a small genome size, short stature, self-pollinating habit, and short generation cycle. To maximize the utility of Brachypodium as a model for basic and applied research it is necessary to develop genomic resources for it. A BAC-based physical map is one of them. A physical map will facilitate analysis of genome structure, comparative genomics, and assembly of the entire genome sequence.

Results

A total of 67,151 Brachypodium BAC clones were fingerprinted with the SNaPshot HICF fingerprinting method and a genome-wide physical map of the Brachypodium genome was constructed. The map consisted of 671 contigs and 2,161 clones remained as singletons. The contigs and singletons spanned 414 Mb. A total of 13,970 gene-related sequences were detected in the BAC end sequences (BES). These gene tags aligned 345 contigs with 336 Mb of rice genome sequence, showing that Brachypodium and rice genomes are generally highly colinear. Divergent regions were mainly in the rice centromeric regions. A dot-plot of Brachypodium contigs against the rice genome sequences revealed remnants of the whole-genome duplication caused by paleotetraploidy, which were previously found in rice and sorghum. Brachypodium contigs were anchored to the wheat deletion bin maps with the BES gene-tags, opening the door to Brachypodium-Triticeae comparative genomics.

Conclusion

The construction of the Brachypodium physical map, and its comparison with the rice genome sequence demonstrated the utility of the SNaPshot-HICF method in the construction of BAC-based physical maps. The map represents an important genomic resource for the completion of Brachypodium genome sequence and grass comparative genomics. A draft of the physical map and its comparisons with rice and wheat are available at .

Sequencing over 13 000 expressed sequence tags from six subtractive cDNA libraries of wild and modern wheats following slow drought stress

A deeper understanding of the drought response and genetic improvement of the cultivated crops for better tolerance requires attention because of the complexity of the drought response syndrome and the loss of genetic diversity during domestication. We initially screened about 200 wild emmer wheat genotypes and then focused on 26 of these lines, which led to the selection of two genotypes with contrasting responses to water deficiency. Six subtractive cDNA libraries were constructed, and over 13 000 expressed sequence tags (ESTs) were sequenced using leaf and root tissues of wild emmer wheat genotypes TR39477 (tolerant) and TTD-22 (sensitive), and modern wheat variety Kiziltan drought stressed for 7 d. Clustering and assembly of ESTs resulted in 2376 unique sequences (1159 without hypothetical proteins and no hits), 75% of which were represented only once. At this level of EST sampling, each tissue shared a very low percentage of transcripts (13-26%). The data obtained indicated that the genotypes shared common elements of drought stress as well as distinctly differential expression patterns that might be illustrative of their contrasting ability to tolerate water deficiencies. The new EST data generated here provide a highly diverse and rich source for gene discovery in wheat and other grasses.

Expressed sequence tag analysis and development of gene associated markers in a near-isogenic plant system of Eragrostis curvula

Eragrostis curvula (Schrad.) Nees is a forage grass native to the semiarid regions of Southern Africa, which reproduces mainly by pseudogamous diplosporous apomixis. A collection of ESTs was generated from four cDNA libraries, three of them obtained from panicles of near-isogenic lines with different ploidy levels and reproductive modes, and one obtained from 12 days-old plant leaves. A total of 12,295 high-quality ESTs were clustered and assembled, rendering 8,864 unigenes, including 1,490 contigs and 7,394 singletons, with a genome coverage of 22%. A total of 7,029 (79.11%) unigenes were functionally categorized by BLASTX analysis against sequences deposited in public databases, but only 37.80% could be classified according to Gene Ontology. Sequence comparison against the cereals genes indexes (GI) revealed 50% significant hits. A total of 254 EST-SSRs were detected from 219 singletons and 35 from contigs. Di- and tri- motifs were similarly represented with percentages of 38.95 and 40.16%, respectively. In addition, 190 SNPs and Indels were detected in 18 contigs generated from 3 to 4 libraries. The ESTs and the molecular markers obtained in this study will provide valuable resources for a wide range of applications including gene identification, genetic mapping, cultivar identification, analysis of genetic diversity, phenotype mapping and marker assisted selection.

Pm37, a new broadly effective powdery mildew resistance gene from Triticum timopheevii

Powdery mildew is an important foliar disease in wheat, especially in areas with a cool or maritime climate. A dominant powdery mildew resistance gene transferred to the hexaploid germplasm line NC99BGTAG11 from T. timopheevii subsp. armeniacum was mapped distally on the long arm of chromosome 7A. Differential reactions were observed between the resistance gene in NC99BGTAG11 and the alleles of the Pm1 locus that is also located on chromosome arm 7AL. Observed segregation in F2:3 lines from the cross NC99BGTAG11xAxminster (Pm1a) demonstrate that germplasm line NC99BGTAG11 carries a novel powdery mildew resistance gene, which is now designated as Pm37. This new gene is highly effective against all powdery mildew isolates tested so far. Analyses of the population with molecular markers indicate that Pm37 is located 16 cM proximal to the Pm1 complex. Simple sequence repeat (SSR) markers Xgwm332 and Xwmc790 were located 0.5 cM proximal and distal, respectively, to Pm37. In order to identify new markers in the region, wheat expressed sequence tags (ESTs) located in the distal 10% of 7AL that were orthologous to sequences from chromosome 6 of rice were targeted. The two new EST-derived STS markers were located distal to Pm37 and one marker was closely linked to the Pm1a region. These new markers can be used in marker-assisted selection schemes to develop wheat cultivars with pyramids of powdery mildew resistance genes, including combinations of Pm37 in coupling linkage with alleles of the Pm1 locus.

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences

Due in part to its small genome (approximately 350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at http://brachypodium.pw.usda.gov . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize.

Comparative analyses between lolium/festuca introgression lines and rice reveal the major fraction of functionally annotated gene models is located in recombination-poor/very recombination-poor regions of the genome

Publication of the rice genome sequence has allowed an in-depth analysis of genome organization in a model monocot plant species. This has provided a powerful tool for genome analysis in large-genome unsequenced agriculturally important monocot species such as wheat, barley, rye, Lolium, etc. Previous data have indicated that the majority of genes in large-genome monocots are located toward the ends of chromosomes in gene-rich regions that undergo high frequencies of recombination. Here we demonstrate that a substantial component of the coding sequences in monocots is localized proximally in regions of very low and even negligible recombination frequencies. The implications of our findings are that during domestication of monocot plant species selection has concentrated on genes located in the terminal regions of chromosomes within areas of high recombination frequency. Thus a large proportion of the genetic variation available for selection of superior plant genotypes has not been exploited. In addition our findings raise the possibility of the evolutionary development of large supergene complexes that confer a selective advantage to the individual.

PCR-based landmark unique gene (PLUG) markers effectively assign homoeologous wheat genes to A, B and D genomes

BACKGROUND

EST-PCR markers normally represent specific products from target genes, and are therefore effective tools for genetic analysis. However, because wheat is an allohexaploid plant, PCR products derived from homoeologous genes are often simultaneously amplified. Such products may be easier to differentiate if they include intron sequences, which are more polymorphic than exon sequences. However, genomic sequence data for wheat are limited; therefore it is difficult to predict the location of introns. By using the similarities in gene structures between rice and wheat, we developed a system called PLUG (PCR-based Landmark Unique Gene) to design primers so that PCR products include intron sequences. We then investigated whether products amplified using such primers could serve as markers able to distinguish multiple products derived from homoeologous genes.

RESULTS

The PLUG system consists of the following steps: (1) Single-copy rice genes (Landmark Unique Gene loci; LUGs) exhibiting high degrees of homology to wheat UniGene sequences are extracted; (2) Alignment analysis is carried out using the LUGs and wheat UniGene sequences to predict exon-exon junctions, and LUGs which can be used to design wheat primers flanking introns (TaEST-LUGs) are extracted; and (3) Primers are designed in an interactive manner. From a total of 4,312 TaEST-LUGs, 24 loci were randomly selected and used to design primers. With all of these primer sets, we obtained specific, intron-containing products from the target genes. These markers were assigned to chromosomes using wheat nullisomic-tetrasomic lines. By PCR-RFLP analysis using agarose gel electrophoresis, 19 of the 24 markers were located on at least one chromosome.

CONCLUSION

In the development of wheat EST-PCR markers capable of efficiently sorting products derived from homoeologous genes, it is important to design primers able to amplify products that include intron sequences with insertion/deletion polymorphisms. Using the PLUG system, wheat EST sequences that can be used for marker development are selected based on comparative genomics with rice, and then primer sets flanking intron sequences are prepared in an interactive, semi-automatic manner. Hence, the PLUG system is an effective tool for large-scale marker development.

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.

Analysis of the wheat endosperm transcriptome

Among the cereals, wheat is the most widely grown geographically and is part of the staple diet in much of the world. Understanding how the cereal endosperm develops and functions will help generate better tools to manipulate grain qualities important to end-users. We used a genomics approach to identify and characterize genes that are expressed in the wheat endosperm. We analyzed the 17,949 publicly available wheat endosperm EST sequences to identify genes involved in the biological processes that occur within this tissue. Clustering and assembly of the ESTs resulted in the identification of 6,187 tentative unique genes, 2,358 of which formed contigs and 3,829 remained as singletons. A BLAST similarity search against the NCBI non-redundant sequence database revealed abundant messages for storage proteins, putative defense proteins, and proteins involved in starch and sucrose metabolism. The level of abundance of the putatively identified genes reflects the physiology of the developing endosperm. Half of the identified genes have unknown functions. Approximately 61% of the endosperm ESTs has been tentatively mapped in the hexaploid wheat genome. Using microarrays for global RNA profiling, we identified endosperm genes that are specifically up regulated in the developing grain.

Triticeae genomics: advances in sequence analysis of large genome cereal crops

Whole genome sequencing provides direct access to all genes of an organism and represents an essential step towards a systematic understanding of (crop) plant biology. Wheat and barley, two of the most important crop species worldwide, have two- to five-fold larger genomes than human - too large to be completely sequenced at current costs. Nevertheless, significant progress has been made to unlock the gene contents of these species by sequencing expressed sequence tags (EST) for high-density mapping and as a basis for elucidating gene function on a large scale. Several megabases of genomic (BAC) sequences have been obtained providing a first insight into the complexity of these huge cereal genomes. However, to fully exploit the information of the wheat and barley genomes for crop improvement, sequence analysis of a significantly larger portion of the Triticeae genomes is needed. In this review an overview of the current status of Triticeae genome sequencing and a perspective concerning future developments in cereal structural genomics is provided.

Application of Genomics to Molecular Breeding of Wheat and Barley

In wheat and barley, several generations of selectable molecular markers have been included in the genetic maps; and a large number of qualitative and quantitative traits were located in the genomes, some of which are being routinely selected in marker-assisted breeding programs. In recent years, a large number of expressed sequence tags (ESTs) have been generated for wheat and barley that have been used for development of functional molecular markers, preparation of transcript maps, and construction of cDNA arrays. These functional genomic resources combined together with new approaches such as expression genetics, association mapping, allele mining, and informatics (bioinformatic tools) possess potential to identify genes responsible for a trait and their deployment in practical plant breeding. High costs currently limit the implementation of functional genomics in breeding programs. The potential applications together with some examples as well as challenges for applying genomics research in breeding activities are discussed. Genomics research will continue to enhance the efficiency and precision for crop improvement but will not replace conventional breeding and evaluation methods.

The development of a PCR-based marker for PSY1 from Hordeum chilense, a candidate gene for carotenoid content accumulation in tritordeum seeds

Hexaploid tritordeums are the amphiploids derived from the cross between the wild barley Hordeum chilense and durum wheat. Tritordeums are characterised by higher yellow pigment content in their seeds than their durum wheat progenitors due to certain H. chilense genes located on the α arm of chromosome 7Hch. In this work a candidate gene approach based on the phytoene synthase gene (PSY) was followed to investigate whether PSY1 may be responsible for the high carotenoid content in tritordeum and to develop a diagnostic marker for H. chilense PSY. This gene codes for the first step in the carotenoid biosynthetic pathway. It was first demonstrated that PSY is duplicated in H. chilense, Triticum urartu, and durum wheat (PSY1 and PSY2), and subsequently a diagnostic cleaved amplified polymorphism (CAP) marker able to differentiate between H. chilense and durum wheat PSY1 was developed. Using this CAP marker and a set of H. chilense-common wheat addition lines it was found that PSY1 is located on the α arm of chromosome 7Hch, where the gene(s) for yellow pigment content are located. PSY1 is located on chromosomes 7A and 7B of durum wheat as demonstrated using Langdon substitution lines. Furthermore, synteny between rice and wheat indicates that PSY1 should be located on the long arms of chromosomes 7A and 7B, in agreement with QTL data for yellow pigment content. Together, these results suggest that PSY1 may be a good candidate gene for further work with yellow pigment content in both durum wheat and tritordeum. In addition, the diagnostic CAP marker developed will be used in our breeding program to transfer H. chilense genes to durum wheat, to evaluate their potential for durum wheat improvement.

Use of a large-scale Triticeae expressed sequence tag resource to reveal gene expression profiles in hexaploid wheat (Triticum aestivum L.)

The US Wheat Genome Project, funded by the National Science Foundation, developed the first large public Triticeae expressed sequence tag (EST) resource. Altogether, 116,272 ESTs were produced, comprising 100,674 5' ESTs and 15 598 3' ESTs. These ESTs were derived from 42 cDNA libraries, which were created from hexaploid bread wheat (Triticum aestivum L.) and its close relatives, including diploid wheat (T. monococcum L. and Aegilops speltoides L.), tetraploid wheat (T. turgidum L.), and rye (Secale cereale L.), using tissues collected from various stages of plant growth and development and under diverse regimes of abiotic and biotic stress treatments. ESTs were assembled into 18,876 contigs and 23,034 singletons, or 41,910 wheat unigenes. Over 90% of the contigs contained fewer than 10 EST members, implying that the ESTs represented a diverse selection of genes and that genes expressed at low and moderate to high levels were well sampled. Statistical methods were used to study the correlation of gene expression patterns, based on the ESTs clustered in the 1536 contigs that contained at least 10 5' EST members and thus representing the most abundant genes expressed in wheat. Analysis further identified genes in wheat that were significantly upregulated (p < 0.05) in tissues under various abiotic stresses when compared with control tissues. Though the function annotation cannot be assigned for many of these genes, it is likely that they play a role associated with the stress response. This study predicted the possible functionality for 4% of total wheat unigenes, which leaves the remaining 96% with their functional roles and expression patterns largely unknown. Nonetheless, the EST data generated in this project provide a diverse and rich source for gene discovery in wheat.

Resistance gene analogues of wheat: molecular genetic analysis of ESTs

Using two divergent nucleotide binding site (NBS) regions from wheat sequences of the NBS-LRR (leucine rich repeat) class, we retrieved 211 wheat and barley NBS-containing resistance gene analogue (RGA) expressed sequence tags (ESTs). These ESTs were grouped into 129 gene sequence groups that contained ESTs that were at least 70% identical at the DNA level over at least 200 bp. Probes were obtained for 89 of these RGA families and chromosome locations were determined for 72 of these probes using nullitetrasomic Chinese Spring wheat lines. RFLP analysis of 49 of these RGA probes revealed 65 mappable polymorphic bands in the doubled haploid Cranbrook x Halberd wheat population (C x H). These bands mapped to 49 loci in C x H. RGA loci were detected on all 21 chromosomes using the nullitetrasomic lines and on 18 chromosomes (linkage groups) in the C x H map. This identified a set of potential markers that could be developed further for use in mapping and ultimately cloning NBS-LRR-type disease resistance genes in wheat.

Structural and functional analyses of the wheat genomes based on expressed sequence tags (ESTs) related to abiotic stresses

To gain insights into the structure and function of the wheat (Triticum aestivum L.) genomes, we identified 278 ESTs related to abiotic stress (cold, heat, drought, salinity, and aluminum) from 7671 ESTs previously mapped to wheat chromosomes. Of the 278 abiotic stress related ESTs, 259 (811 loci) were assigned to chromosome deletion bins and analyzed for their distribution pattern among the 7 homoeologous chromosome groups. Distribution of abiotic stress related EST loci were not uniform throughout the different regions of the chromosomes of the 3 wheat genomes. Both the short and long arms of group 4 chromosomes showed a higher number of loci in their distal regions compared with proximal regions. Of the 811 loci, the number of mapped loci on the A, B, and D genomes were 258, 281, and 272, respectively. The highest number of abiotic stress related loci were found in homoeologous chromosome group 2 (142 loci) and the lowest number were found in group 6 (94 loci). When considering the genome-specific ESTs, the B genome showed the highest number of unique ESTs (7 loci), while none were found in the D genome. Similarly, considering homoeologous group-specific ESTs, group 2 showed the highest number with 16 unique ESTs (58 loci), followed by group 4 with 9 unique ESTs (33 loci). Many of the classified proteins fell into the biological process categories associated with metabolism, cell growth, and cell maintenance. Most of the mapped ESTs fell into the category of enzyme activity (28%), followed by binding activity (27%). Enzymes related to abiotic stress such as β-galactosidase, peroxidase, glutathione reductase, and trehalose-6-phosphate synthase were identified. The comparison of stress-responsive ESTs with genomic sequences of rice (Oryza sativa L.) chromosomes revealed the complexities of colinearity. This bin map provides insight into the structural and functional details of wheat genomic regions in relation to abiotic stress.

EST sequencing and phylogenetic analysis of the model grass Brachypodium distachyon

Brachypodium distachyon (Brachypodium) is a temperate grass with the physical and genomic attributes necessary for a model system (small size, rapid generation time, self-fertile, small genome size, diploidy in some accessions). To increase the utility of Brachypodium as a model grass, we sequenced 20,440 expressed sequence tags (ESTs) from five cDNA libraries made from leaves, stems plus leaf sheaths, roots, callus and developing seed heads. The ESTs had an average trimmed length of 650 bp. Blast nucleotide alignments against SwissProt and GenBank non-redundant databases were performed and a total of 99.9% of the ESTs were found to have some similarity to existing protein or nucleotide sequences. Tentative functional classification of 77% of the sequences was possible by association with gene ontology or clusters of orthologous group's index descriptors. To demonstrate the utility of this EST collection for studying cell wall composition, we identified homologs for the genes involved in the biosynthesis of lignin subunits. A subset of the ESTs was used for phylogenetic analysis that reinforced the close relationship of Brachypodium to wheat and barley.

Wheat EST resources for functional genomics of abiotic stress

Background

Wheat is an excellent species to study freezing tolerance and other abiotic stresses. However, the sequence of the wheat genome has not been completely characterized due to its complexity and large size. To circumvent this obstacle and identify genes involved in cold acclimation and associated stresses, a large scale EST sequencing approach was undertaken by the Functional Genomics of Abiotic Stress (FGAS) project.

Results

We generated 73,521 quality-filtered ESTs from eleven cDNA libraries constructed from wheat plants exposed to various abiotic stresses and at different developmental stages. In addition, 196,041 ESTs for which tracefiles were available from the National Science Foundation wheat EST sequencing program and DuPont were also quality-filtered and used in the analysis. Clustering of the combined ESTs with d2_cluster and TGICL yielded a few large clusters containing several thousand ESTs that were refractory to routine clustering techniques. To resolve this problem, the sequence proximity and "bridges" were identified by an e-value distance graph to manually break clusters into smaller groups. Assembly of the resolved ESTs generated a 75,488 unique sequence set (31,580 contigs and 43,908 singletons/singlets). Digital expression analyses indicated that the FGAS dataset is enriched in stress-regulated genes compared to the other public datasets. Over 43% of the unique sequence set was annotated and classified into functional categories according to Gene Ontology.

Conclusion

We have annotated 29,556 different sequences, an almost 5-fold increase in annotated sequences compared to the available wheat public databases. Digital expression analysis combined with gene annotation helped in the identification of several pathways associated with abiotic stress. The genomic resources and knowledge developed by this project will contribute to a better understanding of the different mechanisms that govern stress tolerance in wheat and other cereals.

A global assembly of cotton ESTs

Approximately 185,000 Gossypium EST sequences comprising >94,800,000 nucleotides were amassed from 30 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including drought stress and pathogen challenges. These libraries were derived from allopolyploid cotton (Gossypium hirsutum; A(T) and D(T) genomes) as well as its two diploid progenitors, Gossypium arboreum (A genome) and Gossypium raimondii (D genome). ESTs were assembled using the Program for Assembling and Viewing ESTs (PAVE), resulting in 22,030 contigs and 29,077 singletons (51,107 unigenes). Further comparisons among the singletons and contigs led to recognition of 33,665 exemplar sequences that represent a nonredundant set of putative Gossypium genes containing partial or full-length coding regions and usually one or two UTRs. The assembly, along with their UniProt BLASTX hits, GO annotation, and Pfam analysis results, are freely accessible as a public resource for cotton genomics. Because ESTs from diploid and allotetraploid Gossypium were combined in a single assembly, we were in many cases able to bioinformatically distinguish duplicated genes in allotetraploid cotton and assign them to either the A or D genome. The assembly and associated information provide a framework for future investigation of cotton functional and evolutionary genomics.

Proteome of amyloplasts isolated from developing wheat endosperm presents evidence of broad metabolic capability

By contrast to chloroplasts, our knowledge of amyloplasts--organelles that synthesize and store starch in heterotrophic plant tissues--is in a formative stage. While our understanding of what is considered their primary function, i.e. the biosynthesis and degradation of starch, has increased dramatically in recent years, relatively little is known about other biochemical processes taking place in these organelles. To help fill this gap, a proteomic analysis of amyloplasts isolated from the starchy endosperm of wheat seeds (10 d post-anthesis) has been conducted. The study has led to the identification of 289 proteins that function in a range of processes, including carbohydrate metabolism, cytoskeleton/plastid division, energetics, nitrogen and sulphur metabolism, nucleic acid-related reactions, synthesis of various building blocks, protein-related reactions, transport, signalling, stress, and a variety of other activities grouped under 'miscellaneous'. The function of 12% of the proteins was unknown. The results highlight the role of the amyloplast as a starch-storing organelle that fulfills a spectrum of biosynthetic needs of the parent tissue. When compared with a recent proteomic analysis of whole endosperm, the current study demonstrates the advantage of using isolated organelles in proteomic studies.

Transcriptome analysis of salinity stress responses in common wheat using a 22k oligo-DNA microarray

In this study, we constructed a 22k wheat oligo-DNA microarray. A total of 148,676 expressed sequence tags of common wheat were collected from the database of the Wheat Genomics Consortium of Japan. These were grouped into 34,064 contigs, which were then used to design an oligonucleotide DNA microarray. Following a multistep selection of the sense strand, 21,939 60-mer oligo-DNA probes were selected for attachment on the microarray slide. This 22k oligo-DNA microarray was used to examine the transcriptional response of wheat to salt stress. More than 95% of the probes gave reproducible hybridization signals when targeted with RNAs extracted from salt-treated wheat shoots and roots. With the microarray, we identified 1,811 genes whose expressions changed more than 2-fold in response to salt. These included genes known to mediate response to salt, as well as unknown genes, and they were classified into 12 major groups by hierarchical clustering. These gene expression patterns were also confirmed by real-time reverse transcription-PCR. Many of the genes with unknown function were clustered together with genes known to be involved in response to salt stress. Thus, analysis of gene expression patterns combined with gene ontology should help identify the function of the unknown genes. Also, functional analysis of these wheat genes should provide new insight into the response to salt stress. Finally, these results indicate that the 22k oligo-DNA microarray is a reliable method for monitoring global gene expression patterns in wheat.

Analysis of expressed sequence tags and the identification of associated short tandem repeats in switchgrass

Switchgrass is a large, North American, perennial grass that is being evaluated as a potential energy crop. Expressed sequence tags (ESTs) were generated from four switchgrass cv. "Kanlow" cDNA libraries to create a gene inventory of 7,810 unique gene clusters from a total of 11,990 individual sequences. Blast similarity searches to SwissProt and GenBank non-redundant protein and nucleotide databases were performed and a total of 79% of these unique clusters were found to be similar to existing protein or nucleotide sequences. Tentative functional classification of 61% of the sequences was possible by association with appropriate gene ontology descriptors. Significant differential representation between genes in leaf, stem, crown, and callus libraries was observed for many highly expressed genes The unique gene clusters were screened for the presence of short tandem repeats for further development as microsatellite markers. A total of 334 gene clusters contained repeats representing 3.8% of the ESTs queried.

Developmental changes in the metabolic protein profiles of wheat endosperm

A combined two-dimensional gel electrophoresis-mass spectrometry approach was utilized to identify over 250 proteins of wheat (Triticum aestivum L., cv. Butte 86) starchy endosperm that participate in 13 biochemical processes: ATP interconversion reactions, carbohydrate metabolism, cell division, cytoskeleton, lipid metabolism, nitrogen metabolism, protein synthesis/assembly, protein turnover, signal transduction, protein storage, stress/defense, transcription/translation, and transport. Endosperm protein populations were compared at early (10 days post-anthesis, dpa) and late (36 dpa) stages of grain development. Analysis of protein number and spot volume revealed that carbohydrate metabolism, transcription/translation, and protein synthesis/assembly were the principal endosperm functions at 10 dpa followed by nitrogen metabolism, protein turnover, cytoskeleton, cell division, signal transduction, and lipid metabolism. Carbohydrate metabolism and protein synthesis/assembly were also major functions at 36 dpa, but stress/defense and storage were predominant. The results provide insight into biochemical events taking place during wheat grain development and highlight the value of proteomics in characterizing complex biochemical processes. Further, the proteome maps will facilitate future studies addressing the effects of genetic and environmental factors on the development and quality of wheat grain.

Analysis of expressed sequence tags and the identification of associated short tandem repeats in switchgrass

Switchgrass is a large, North American, perennial grass that is being evaluated as a potential energy crop. Expressed sequence tags (ESTs) were generated from four switchgrass cv. "Kanlow" cDNA libraries to create a gene inventory of 7,810 unique gene clusters from a total of 11,990 individual sequences. Blast similarity searches to SwissProt and GenBank non-redundant protein and nucleotide databases were performed and a total of 79% of these unique clusters were found to be similar to existing protein or nucleotide sequences. Tentative functional classification of 61% of the sequences was possible by association with appropriate gene ontology descriptors. Significant differential representation between genes in leaf, stem, crown, and callus libraries was observed for many highly expressed genes The unique gene clusters were screened for the presence of short tandem repeats for further development as microsatellite markers. A total of 334 gene clusters contained repeats representing 3.8% of the ESTs queried.

Tempos of gene locus deletions and duplications and their relationship to recombination rate during diploid and polyploid evolution in the Aegilops-Triticum alliance

The origin of tetraploid wheat and the divergence of diploid ancestors of wheat A and D genomes were estimated to have occurred 0.36 and 2.7 million years ago, respectively. These estimates and the evolutionary history of 3159 gene loci were used to estimate the rates with which gene loci have been deleted and duplicated during the evolution of wheat diploid ancestors and during the evolution of polyploid wheat. During diploid evolution, the deletion rate was 2.1 x 10(-3) locus(-1) MY(-1) for single-copy loci and 1.0 x 10(-2) locus(-1) MY(-1) for loci in paralogous sets. Loci were duplicated with a rate of 2.9 x 10(-3) locus(-1) MY(-1) during diploid evolution. During polyploid evolution, locus deletion and locus duplication rates were 1.8 x 10(-2) and 1.8 x 10(-3) locus(-1) MY(-1), respectively. Locus deletion and duplication rates correlated positively with the distance of the locus from the centromere and the recombination rate during diploid evolution. The functions of deleted and duplicated loci were inferred to gain insight into the surprisingly high rate of deletions of loci present apparently only once in a genome. The significance of these findings for genome evolution at the diploid and polyploid level is discussed.

Deletion mapping of homoeologous group 6-specific wheat expressed sequence tags

To localize wheat (Triticum aestivum L.) ESTs on chromosomes, 882 homoeologous group 6-specific ESTs were identified by physically mapping 7965 singletons from 37 cDNA libraries on 146 chromosome, arm, and sub-arm aneuploid and deletion stocks. The 882 ESTs were physically mapped to 25 regions (bins) flanked by 23 deletion breakpoints. Of the 5154 restriction fragments detected by 882 ESTs, 2043 (loci) were localized to group 6 chromosomes and 806 were mapped on other chromosome groups. The number of loci mapped was greatest on chromosome 6B and least on 6D. The 264 ESTs that detected orthologous loci on all three homoeologs using one restriction enzyme were used to construct a consensus physical map. The physical distribution of ESTs was uneven on chromosomes with a tendency toward higher densities in the distal halves of chromosome arms. About 43% of the wheat group 6 ESTs identified rice homologs upon comparisons of genome sequences. Fifty-eight percent of these ESTs were present on rice chromosome 2 and the remaining were on other rice chromosomes. Even within the group 6 bins, rice chromosomal blocks identified by 1-6 wheat ESTs were homologous to up to 11 rice chromosomes. These rice-block contigs were used to resolve the order of wheat ESTs within each bin.

Construction and evaluation of cDNA libraries for large-scale expressed sequence tag sequencing in wheat (Triticum aestivum L.)

A total of 37 original cDNA libraries and 9 derivative libraries enriched for rare sequences were produced from Chinese Spring wheat (Triticum aestivum L.), five other hexaploid wheat genotypes (Cheyenne, Brevor, TAM W101, BH1146, Butte 86), tetraploid durum wheat (T. turgidum L.), diploid wheat (T. monococcum L.), and two other diploid members of the grass tribe Triticeae (Aegilops speltoides Tausch and Secale cereale L.). The emphasis in the choice of plant materials for library construction was reproductive development subjected to environmental factors that ultimately affect grain quality and yield, but roots and other tissues were also included. Partial cDNA expressed sequence tags (ESTs) were examined by various measures to assess the quality of these libraries. All ESTs were processed to remove cloning system sequences and contaminants and then assembled using CAP3. Following these processing steps, this assembly yielded 101,107 sequences derived from 89,043 clones, which defined 16,740 contigs and 33,213 singletons, a total of 49,953 "unigenes." Analysis of the distribution of these unigenes among the libraries led to the conclusion that the enrichment methods were effective in reducing the most abundant unigenes and to the observation that the most diverse libraries were from tissues exposed to environmental stresses including heat, drought, salinity, or low temperature.