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

Potential role of the regulatory miR1119-MYC2 module in wheat (Triticum aestivum L.) drought tolerance

MicroRNA (miRNA)-target gene modules are essential components of plants' abiotic stress signalling pathways Little is known about the drought-responsive miRNA-target modules in wheat, but systems biology approaches have enabled the prediction of these regulatory modules and systematic study of their roles in responses to abiotic stresses. Using such an approach, we sought miRNA-target module(s) that may be differentially expressed under drought and non-stressed conditions by mining Expressed Sequence Tag (EST) libraries of wheat roots and identified a strong candidate (miR1119-MYC2). We then assessed molecular and physiochemical differences between two wheat genotypes with contrasting drought tolerance in a controlled drought experiment and assessed possible relationships between their tolerance and evaluated traits. We found that the miR1119-MYC2 module significantly responds to drought stress in wheat roots. It is differentially expressed between the contrasting wheat genotypes and under drought versus non-stressed conditions. We also found significant associations between the module's expression profiles and ABA hormone content, water relations, photosynthetic activities, H₂O₂ levels, plasma membrane damage, and antioxidant enzyme activities in wheat. Collectively, our results suggest that a regulatory module consisting of miR1119 and MYC2 may play an important role in wheat's drought tolerance.

MicroRNA miR1118 contributes to wheat (Triticum aestivum L.) salinity tolerance by regulating the Plasma Membrane Intrinsic Proteins1;5 (PIP1;5) gene

microRNAs (miRNAs) are important regulators of various adaptive stress responses in crops; however, many details about associations among miRNAs, their target genes and physiochemical responses of crops under salinity stress remain poorly understood. We designed this study in a systems biology context and used a collection of computational, experimental and statistical procedures to uncover some regulatory functions of miRNAs in the response of the important crop, wheat, to salinity stress. Accordingly, under salinity conditions, wheat roots' Expressed Sequence Tag (EST) libraries were computationally mined to identify the most reliable differentially expressed miRNA and its related target gene(s). Then, molecular and physiochemical evaluations were carried out in a separate salinity experiment using two contrasting wheat genotypes. Finally, the association between changes in measured characteristics and wheat salinity tolerance was determined. From the results, miR1118 was assigned as a reliable salinity-responsive miRNA in wheat roots. The expression profiles of miR1118 and its predicted target gene, Plasma Membrane Intrinsic Proteins1,5 (PIP1;5), significantly differed between wheat genotypes. Moreover, results revealed that expression profiles of miR1118 and PIP1;5 significantly correlate to Relative Water Content (RWC), root hydraulic conductance (Lp), photosynthetic activities, plasma membrane damages, osmolyte accumulation and ion homeostasis of wheat. Our results suggest a plausible regulatory node through miR1118 adjusting the wheat water status, maintaining ion homeostasis and mitigating membrane damages, mainly through the PIP1;5 gene, under salinity conditions. To our knowledge, this is the first report on the role of miR1118 and PIP1;5 in wheat salinity response.

A systems biology study unveils the association between a melatonin biosynthesis gene, O-methyl transferase 1 (OMT1) and wheat (Triticum aestivum L.) combined drought and salinity stress tolerance

Enhanced levels of endogenous melatonin in the root of wheat, mainly through the OMT1 gene, augment the antioxidant system, reestablish redox homeostasis and are associated with combined stress tolerance. A systems biology approach, including a collection of computational analyses and experimental assays, led us to uncover some aspects of a poorly understood phenomenon, namely wheat (Triticum aestivum L.) combined drought and salinity stress tolerance. Accordingly, a cross-study comparison of stress experiments was performed via a meta-analysis of Expressed Sequence Tags (ESTs) data from wheat roots to uncover the overlapping gene network of drought and salinity stresses. Identified differentially expressed genes were functionally annotated by gene ontology enrichment analysis and gene network analysis. Among those genes, O-methyl transferase 1 (OMT1) was highlighted as a more important (hub) gene in the dual-stress response gene network. Afterwards, the potential roles of OMT1 in mediating physiochemical indicators of stress tolerance were investigated in two wheat genotypes differing in abiotic stress tolerance. Regression analysis and correspondence analysis (CA) confirmed that the expression profiles of the OMT1 gene and variations in melatonin content, antioxidant enzyme activities, proline accumulation, H₂O₂ and malondialdehyde (MDA) contents are significantly associated with combined stress tolerance. These results reveal that the OMT1 gene may contribute to wheat combined drought and salinity stress tolerance through augmenting the antioxidant system and re-establishing redox homeostasis, probably via the regulation of melatonin biosynthesis as a master regulator molecule. Our findings provide new insights into the roles of melatonin in wheat combined drought and salinity stress tolerance and suggest a novel plausible regulatory node through the OMT1 gene to improve multiple-stress tolerant crops.

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.

Two wheat (Triticum aestivum) pathogenesis-related 10 (PR-10) transcripts with distinct patterns of abundance in different organs

PR-10 genes encode small, acidic, intracellular proteins that respond to abiotic and biotic stimuli. Transgenic expression of PR-10 genes has been shown to enhance early seedling growth of dicots in saline environments. To identify candidate PR-10 genes in cereals for increasing stress tolerance, we conducted phylogenetic analyses and real-time polymerase chain reaction of representatives of the two major clades of putative PR-10 genes in wheat. We observed that the abundance of BQ752893 was generally greater than the abundance of CV778999, particularly when measured in roots across four wheat genotypes. However, CV778999 transcripts were more abundant than BQ752893 in flag leaves. These data suggest that the transcripts define two functionally divergent groups of PR-10 type genes in wheat, both of which may be suitable targets for biotechnological manipulation under different circumstances.

The biocontrol bacterium Pseudomonas fluorescens Pf29Arp strain affects the pathogenesis-related gene expression of the take-all fungus Gaeumannomyces graminis var. tritici on wheat roots

The main effects of antagonistic rhizobacteria on plant pathogenic fungi are antibiosis, fungistasis or an indirect constraint through the induction of a plant defence response. To explore different biocontrol mechanisms, an in vitro confrontation assay was conducted with the rhizobacterium Pseudomonas fluorescens Pf29Arp as a biocontrol agent of the fungus Gaeumannomyces graminis var. tritici (Ggt) on wheat roots. In parallel with the assessment of disease extension, together with the bacterial and fungal root colonization rates, the transcript levels of candidate fungal pathogenicity and plant-induced genes were monitored during the 10-day infection process. The bacterial inoculation of wheat roots with the Pf29Arp strain reduced the development of Ggt-induced disease expressed as attack frequency and necrosis length. The growth rates of Ggt and Pf29Arp, monitored through quantitative polymerase chain reaction of DNA amounts with a part of the Ggt 18S rDNA gene and a specific Pf29Arp strain detection probe, respectively, increased throughout the interactions. Bacterial antagonism and colonization had no significant effect on root colonization by Ggt. The expression of fungal and plant genes was quantified in planta by quantitative reverse transcription-polymerase chain reaction during the interactions thanks to the design of specific primers and an innovative universal reference system. During the early stages of the tripartite interaction, several of the fungal genes assayed were down-regulated by Pf29Arp, including two laccases, a β-1,3-exoglucanase and a mitogen-activated protein kinase. The plant host glutathione-S-transferase gene was induced by Ggt alone and up-regulated by Pf29Arp bacteria in interaction with the pathogen. We conclude that Pf29Arp antagonism acts through the alteration of fungal pathogenesis and probably through the activation of host defences.

Analysis of expressed sequence tags from a single wheat cultivar facilitates interpretation of tandem mass spectrometry data and discrimination of gamma gliadin proteins that may play different functional roles in flour

BACKGROUND

The gamma gliadins are a complex group of proteins that together with other gluten proteins determine the functional properties of wheat flour. The proteins have unusually high levels of glutamine and proline and contain large regions of repetitive sequences. While most gamma gliadins are monomeric proteins containing eight conserved cysteine residues, some contain an additional cysteine residue that enables them to be linked with other gluten proteins into large polymers that are critical for flour quality. The ability to differentiate among the gamma gliadins is important for studies of wheat flour quality because proteins with similar sequences can have different effects on functional properties.

RESULTS

The complement of gamma gliadin genes expressed in the wheat cultivar Butte 86 was evaluated by analyzing publicly available expressed sequence tag (EST) data. Eleven contigs were assembled from 153 Butte 86 ESTs. Nine of the contigs encoded full-length proteins and four of the proteins contained nine cysteine residues. Only one of the encoded proteins was a perfect match with a sequence reported in NCBI. Contigs from four different publicly available EST assemblies encoded proteins that were perfect matches with some, but not all, of the Butte 86 gamma gliadins and the complement of identical proteins was different for each assembly. A specialized database that included the sequences of Butte 86 gamma gliadins was constructed for identification of flour proteins by tandem mass spectrometry (MS/MS). In a pilot experiment, proteins corresponding to six Butte 86 gamma gliadin contigs were distinguished by MS/MS, including one containing the extra cysteine residue. Two other proteins were identified as one of two closely related Butte 86 proteins but could not be distinguished unequivocally. Unique peptide tags specific for Butte 86 gamma gliadins are reported.

CONCLUSIONS

Inclusion of cultivar-specific gamma gliadin sequences in databases maximizes the number and quality of peptide identifications and increases sequence coverage of these gamma gliadins by MS/MS. This approach makes it possible to distinguish closely related proteins, to associate individual proteins with sequences of specific genes, and to evaluate proteomic data in a biological context to better address questions about wheat flour quality.

Parallel evolution of nacre building gene sets in molluscs

The capacity to biomineralize is closely linked to the rapid expansion of animal life during the early Cambrian, with many skeletonized phyla first appearing in the fossil record at this time. The appearance of disparate molluscan forms during this period leaves open the possibility that shells evolved independently and in parallel in at least some groups. To test this proposition and gain insight into the evolution of structural genes that contribute to shell fabrication, we compared genes expressed in nacre (mother-of-pearl) forming cells in the mantle of the bivalve Pinctada maxima and the gastropod Haliotis asinina. Despite both species having highly lustrous nacre, we find extensive differences in these expressed gene sets. Following the removal of housekeeping genes, less than 10% of all gene clusters are shared between these molluscs, with some being conserved biomineralization genes that are also found in deuterostomes. These differences extend to secreted proteins that may localize to the organic shell matrix, with less than 15% of this secretome being shared. Despite these differences, H. asinina and P. maxima both secrete proteins with repetitive low-complexity domains (RLCDs). Pinctada maxima RLCD proteins-for example, the shematrins-are predominated by silk/fibroin-like domains, which are absent from the H. asinina data set. Comparisons of shematrin genes across three species of Pinctada indicate that this gene family has undergone extensive divergent evolution within pearl oysters. We also detect fundamental bivalve-gastropod differences in extracellular matrix proteins involved in mollusc-shell formation. Pinctada maxima expresses a chitin synthase at high levels and several chitin deacetylation genes, whereas only one protein involved in chitin interactions is present in the H. asinina data set, suggesting that the organic matrix on which calcification proceeds differs fundamentally between these species. Large-scale differences in genes expressed in nacre-forming cells of Pinctada and Haliotis are compatible with the hypothesis that gastropod and bivalve nacre is the result of convergent evolution. The expression of novel biomineralizing RLCD proteins in each of these two molluscs and, interestingly, sea urchins suggests that the evolution of such structural proteins has occurred independently multiple times in the Metazoa.

Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum

SUMMARY The fungal pathogen Sclerotinia sclerotiorum infects a broad range of dicotyledonous plant species and causes stem rot in Brassica napus. To elucidate the mechanisms underlying the defence response, the patterns of gene expression in the partially resistant B. napus cultivar ZhongYou 821 (ZY821) and the susceptible cultivar Westar were studied using a B. napus oligonucleotide microarray. Although maximum differential gene expression was observed at 48 h post-inoculation (hpi) in both cultivars, increased transcript levels were detected in cv. ZY821 at the earlier stages of infection (6-12 hpi) for many genes, including those encoding defence-associated proteins, such as chitinases, glucanases, osmotins and lectins, as well as genes encoding transcription factors belonging to the zinc finger, WRKY, APETALA2 (AP2) and MYB classes. In both cultivars, genes encoding enzymes involved in jasmonic acid, ethylene and auxin synthesis were induced, as were those for gibberellin degradation. In addition, changes in the expression of genes encoding enzymes involved in carbohydrate and energy metabolism appeared to be directed towards shuttling carbon reserves to the tricarboxylic acid cycle and generating reactive oxygen species. Transcripts from genes encoding enzymes involved in glucosinolate and phenylpropanoid biosynthesis were highly elevated in both cultivars, suggesting that secondary metabolites are also components of the response to S. sclerotiorum in B. napus.

Transcriptome pathways unique to dehydration tolerant relatives of modern wheat

Among abiotic stressors, drought is a major factor responsible for dramatic yield loss in agriculture. In order to reveal differences in global expression profiles of drought tolerant and sensitive wild emmer wheat genotypes, a previously deployed shock-like dehydration process was utilized to compare transcriptomes at two time points in root and leaf tissues using the Affymetrix GeneChip(R) Wheat Genome Array hybridization. The comparison of transcriptomes reveal several unique genes or expression patterns such as differential usage of IP(3)-dependent signal transduction pathways, ethylene- and abscisic acid (ABA)-dependent signaling, and preferential or faster induction of ABA-dependent transcription factors by the tolerant genotype that distinguish contrasting genotypes indicative of distinctive stress response pathways. The data also show that wild emmer wheat is capable of engaging known drought stress responsive mechanisms. The global comparison of transcriptomes in the absence of and after dehydration underlined the gene networks especially in root tissues that may have been lost in the selection processes generating modern bread wheats.

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.

Data mining for miRNAs and their targets in the Triticeae

MicroRNAs (miRNAs) and the mRNA targets of miRNAs were identified by sequence complementarity within a DNA sequence database for species of the Triticeae. Data screening identified 28 miRNA precursor sequences from 15 miRNA families that contained conserved mature miRNA sequences within predicted stem-loop structures. In addition, the identification of 337 target sequences among Triticeae genes provided further evidence of the existence of 26 miRNA families in the cereals. MicroRNA targets included genes that are homologous to known targets in diverse model species as well as novel targets. MicroRNA precursors and targets were identified in 10 related species, though the great majority of them were identified in bread wheat, Triticum aestivum, and barley, Hordeum vulgare, the two species with the largest EST data sets among the Triticeae.

Genome-wide analysis for identification of salt-responsive genes in common wheat

To identify salt-responsive genes in wheat, global expression analysis of transcripts was carried out using oligo-DNA microarrays. Microarrays have been designed from approximately 32,000 unique wheat genes classified from a large number of expressed sequence tags (ESTs). Two-week-old seedlings of wheat were treated with 150 mM NaCl for 1, 6, and 24 h, and their roots and shoots were separately subjected to analyses. Consequently, 5,996 genes showed changes in expression of more than twofold and were classified into 12 groups according to correlations in expression patterns. These salt-responsive genes were assigned functions using the Gene Ontology (GO). Genes assigned to transcription factor, transcription-regulator activity, and DNA-binding functions were preferentially classified into early response groups. On the other hand, those assigned transferase and transporter activity were classified into late response groups. These data suggest that multiple signal transduction pathways in response to salinity exist in wheat. Transcription factors (TFs) which have been reported as participants in salt-tolerant pathway changed their expression levels in response to salt treatment. Among them, only a few TFs show high sequence homologies to genes in rice. These investigations suggest that salt-responsive genes identified by this study are candidates for salt-stress tolerance uniquely in wheat.

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.