Comparative transcriptomics in the Triticeae

Convergently selected NPF2.12 coordinates root growth and nitrogen use efficiency in wheat and barley

Understanding the genetic and molecular function of nitrate sensing and acquisition across crop species will accelerate breeding of cultivars with improved nitrogen use efficiency (NUE). Here, we performed a genome-wide scan using wheat and barley accessions characterized under low and high N inputs that uncovered the NPF2.12 gene, encoding a homolog of the Arabidopsis nitrate transceptor NRT1.6 and other low-affinity nitrate transporters that belong to the MAJOR FACILITATOR SUPERFAMILY. Next, it is shown that variations in the NPF2.12 promoter correlated with altered NPF2.12 transcript levels where decreased gene expression was measured under low nitrate availability. Multiple field trials revealed a significantly enhanced N content in leaves and grains and NUE in the presence of the elite allele TaNPF2.12^TT grown under low N conditions. Furthermore, the nitrate reductase encoding gene NIA1 was up-regulated in npf2.12 mutant upon low nitrate concentrations, thereby resulting in elevated levels of nitric oxide (NO) production. This increase in NO correlated with the higher root growth, nitrate uptake, and N translocation observed in the mutant when compared to wild-type. The presented data indicate that the elite haplotype alleles of NPF2.12 are convergently selected in wheat and barley that by inactivation indirectly contribute to root growth and NUE by activating NO signaling under low nitrate conditions.

Characterization of the heterotrimeric G protein gene families in Triticum aestivum and related species

This study characterizes the heterotrimeric G protein gene families in Triticum aestivum, their tissue-specific expression patterns during development and in response to biotic and abiotic stress conditions. There are three Gα genes, three Gβ and 12 Gγ genes, totaling 18 genes encoding heterotrimeric G proteins in the hexaploid wheat genome. Each haploid genome of the hexaploid T. aestivum has a single gene encoding the α subunit of the heterotrimeric G protein complex, GA1, a single Gβ and four Gγ genes. Each gene has three homeologous copies in the A, B and D genomes. The physical interaction between the Gβ (Gpb) and two Gγ subunits, Gpg1 and Gpg2, was shown through bimolecular fluorescence complementation (BiFC). The gene expression in response to biotic and abiotic stresses showed both up-regulation and down-regulation of members of the gene families. Gγ2-B and Gγ2-D are significantly upregulated during heat stress, GA1-D is upregulated by cold stress and Gγ1-A and Gγ1-D were upregulated by Fusarium graminearum inoculation in a F. graminearum resistant cultivar. This suggests that these members may play roles in biotic and abiotic signaling pathways and the roles of these genes within these pathways need further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03156-9.

Annotation and Molecular Characterisation of the TaIRO3 and TaHRZ Iron Homeostasis Genes in Bread Wheat (Triticum aestivum L.)

Effective maintenance of plant iron (Fe) homoeostasis relies on a network of transcription factors (TFs) that respond to environmental conditions and regulate Fe uptake, translocation, and storage. The iron-related transcription factor 3 (IRO3), as well as haemerythrin motif-containing really interesting new gene (RING) protein and zinc finger protein (HRZ), are major regulators of Fe homeostasis in diploid species like Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa L.), but remain uncharacterised in hexaploid bread wheat (Triticum aestivum L.). In this study, we have identified, annotated, and characterised three TaIRO3 homoeologs and six TaHRZ1 and TaHRZ2 homoeologs in the bread wheat genome. Protein analysis revealed that TaIRO3 and TaHRZ proteins contain functionally conserved domains for DNA-binding, dimerisation, Fe binding, or polyubiquitination, and phylogenetic analysis revealed clustering of TaIRO3 and TaHRZ proteins with other monocot IRO3 and HRZ proteins, respectively. Quantitative reverse-transcription PCR analysis revealed that all TaIRO3 and TaHRZ homoeologs have unique tissue expression profiles and are upregulated in shoot tissues in response to Fe deficiency. After 24 h of Fe deficiency, the expression of TaHRZ homoeologs was upregulated, while the expression of TaIRO3 homoeologs was unchanged, suggesting that TaHRZ functions upstream of TaIRO3 in the wheat Fe homeostasis TF network.

Genome-wide identification and characterization of the GDP-L-galactose phosphorylase gene family in bread wheat

BACKGROUND

Ascorbate is a powerful antioxidant in plants and an essential micronutrient for humans. The GDP-L-galactose phosphorylase (GGP) gene encodes the rate-limiting enzyme of the L-galactose pathway-the dominant ascorbate biosynthetic pathway in plants-and is a promising gene candidate for increasing ascorbate in crops. In addition to transcriptional regulation, GGP production is regulated at the translational level through an upstream open reading frame (uORF) in the long 5'-untranslated region (5'UTR). The GGP genes have yet to be identified in bread wheat (Triticum aestivum L.), one of the most important food grain sources for humans.

RESULTS

Bread wheat chromosomal groups 4 and 5 were found to each contain three homoeologous TaGGP genes on the A, B, and D subgenomes (TaGGP2-A/B/D and TaGGP1-A/B/D, respectively) and a highly conserved uORF was present in the long 5'UTR of all six genes. Phylogenetic analyses demonstrated that the TaGGP genes separate into two distinct groups and identified a duplication event of the GGP gene in the ancestor of the Brachypodium/Triticeae lineage. A microsynteny analysis revealed that the TaGGP1 and TaGGP2 subchromosomal regions have no shared synteny suggesting that TaGGP2 may have been duplicated via a transposable element. The two groups of TaGGP genes have distinct expression patterns with the TaGGP1 homoeologs broadly expressed across different tissues and developmental stages and the TaGGP2 homoeologs highly expressed in anthers. Transient transformation of the TaGGP coding sequences in Nicotiana benthamiana leaf tissue increased ascorbate concentrations more than five-fold, confirming their functional role in ascorbate biosynthesis in planta.

CONCLUSIONS

We have identified six TaGGP genes in the bread wheat genome, each with a highly conserved uORF. Phylogenetic and microsynteny analyses highlight that a transposable element may have been responsible for the duplication and specialized expression of GGP2 in anthers in the Brachypodium/Triticeae lineage. Transient transformation of the TaGGP coding sequences in N. benthamiana demonstrated their activity in planta. The six TaGGP genes and uORFs identified in this study provide a valuable genetic resource for increasing ascorbate concentrations in bread wheat.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

Genome-wide and evolutionary analysis of the class III peroxidase gene family in wheat and Aegilops tauschii reveals that some members are involved in stress responses

Background

The class III peroxidase (PRX) gene family is a plant-specific member of the PRX superfamily that is closely related to various physiological processes, such as cell wall loosening, lignification, and abiotic and biotic stress responses. However, its classification, evolutionary history and gene expression patterns are unclear in wheat and Aegilops tauschii.

Results

Here, we identified 374, 159 and 169 PRXs in Triticum aestivum, Triticum urartu and Ae. tauschii, respectively. Together with PRXs detected from eight other plants, they were classified into 18 subfamilies. Among subfamilies V to XVIII, a conserved exon-intron structure within the “001” exon phases was detected in the PRX domain. Based on the analysis, we proposed a phylogenetic model to infer the evolutionary history of the exon-intron structures of PRX subfamilies. A comparative genomics analysis showed that subfamily VII could be the ancient subfamily that originated from green algae (Chlamydomonas reinhardtii). Further integrated analysis of chromosome locations and collinearity events of PRX genes suggested that both whole genome duplication (WGD) and tandem duplication (TD) events contributed to the expansion of T. aestivum PRXs (TaePRXs) during wheat evolution. To validate functions of these genes in the regulation of various physiological processes, the expression patterns of PRXs in different tissues and under various stresses were studied using public microarray datasets. The results suggested that there were distinct expression patterns among different tissues and PRXs could be involved in biotic and abiotic responses in wheat. qRT-PCR was performed on samples exposed to drought, phytohormone treatments and Fusarium graminearum infection to validate the microarray predictions. The predicted subcellular localizations of some TaePRXs were consistent with the confocal microscopy results. We predicted that some TaePRXs had hormone-responsive cis-elements in their promoter regions and validated these predicted cis-acting elements by sequencing promoters.

Conclusion

In this study, identification, classification, evolution, and expression patterns of PRXs in wheat and relative plants were performed. Our results will provide information for further studies on the evolution and molecular mechanisms of wheat PRXs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6006-5) contains supplementary material, which is available to authorized users.

Metabolic engineering of bread wheat improves grain iron concentration and bioavailability

Bread wheat (Triticum aestivum L.) is cultivated on more land than any other crop and produces a fifth of the calories consumed by humans. Wheat endosperm is rich in starch yet contains low concentrations of dietary iron (Fe) and zinc (Zn). Biofortification is a micronutrient intervention aimed at increasing the density and bioavailability of essential vitamins and minerals in staple crops; Fe biofortification of wheat has proved challenging. In this study we employed constitutive expression (CE) of the rice (Oryza sativa L.) nicotianamine synthase 2 (OsNAS2) gene in bread wheat to up-regulate biosynthesis of two low molecular weight metal chelators - nicotianamine (NA) and 2'-deoxymugineic acid (DMA) - that play key roles in metal transport and nutrition. The CE-OsNAS2 plants accumulated higher concentrations of grain Fe, Zn, NA and DMA and synchrotron X-ray fluorescence microscopy (XFM) revealed enhanced localization of Fe and Zn in endosperm and crease tissues, respectively. Iron bioavailability was increased in white flour milled from field-grown CE-OsNAS2 grain and positively correlated with NA and DMA concentrations.

Subtilase activity and gene expression during germination and seedling growth in barley

Proteases play a main role in the mobilization of storage proteins during seed germination. Until today, there is little information about the involvement of serine proteases, particularly subtilases, in the germination of barley grains. The aims of the present work were to study the contribution of serine proteases to the total proteolytic activity induced during germination of barley grains and evaluate the specific involvement of subtilases in this process. Proteolytic activity assayed against azocasein in the presence of specific inhibitors, showed that serine proteases contributed between 10 and 20% of total activity along germination. Subtilase activity increased from day 1 after imbibition with a peak between days 4-5. Moreover, in vivo determination of subtilase activity in germinating grains revealed increasing activity along germination mainly localized in the seed endosperm and developing rootlets. Finally, the expression of 19 barley genes encoding subtilases was measured by real time PCR during germination. Three of the analyzed genes increased their expression along germination, five showed a transient induction, one was down-regulated, nine remained unchanged and one was not expressed. The present work demonstrates the involvement of subtilases in germination of barley grains and describes the positive association of eight subtilase genes to this process.

Characterization of the Esi3/RCI2/PMP3 gene family in the Triticeae

Background

Members of the Early Salt Induced 3 (Esi3/RCI2/PMP3) gene family in plants have been shown to be induced in response to both biotic and abiotic stresses and to enhance stress tolerance in both transgenic plants and Saccharomyces cerevisiae. Esi3 was first identified as a salt stress induced gene in the salt tolerant wild wheat grass, Lophopyrum elongatum, and subsequently homologous genes in many other species were found to be members of the gene family. These include Arabidopsis thaliana and Oryza sativa where they are referred to as Rare Cold Inducible 2 (RCI2), and Zea mays where they are referred to as Plasma Membrane Protein 3 (PMP3). This study characterizes the Esi3 family members in Triticum aestivum and explores the tissue specific expression patterns of the gene family members as well as their response to a variety of environmental stresses.

Results

The Esi3 gene family was found to have a total of 29 family members comprised of ten paralogous groups in the hexaploid T. aestivum. Each paralogous group contains three homeologous copies, one in each of the A, B and D genomes with the exception of Esi3–2 which is missing the B copy. The genes of the Esi3 gene family were also identified in four other monocot species, Aegilops tauschii, Hordeum vulgare, Secale cereale and Sorghum bicolor, and were confirmed or corrected for Brachypodium distachyon, Oryza sativa and Zea mays, as well as the dicot Arabidopsis thaliana. Gene expression of the Esi3s was analyzed using tissue-specific, abiotic and biotic stress RNA-Seq 454 sequence libraries and Affymetrix microarray data for T. aestivum.

Conclusions

Members of nearly all paralogous groups of the Esi3 genes in T. aestivum have altered gene expression in response to abiotic or biotic stress conditions. In addition, there are modest differences in gene expression among homeologous members of the gene family. This suggests that the Esi3 gene family plays an important role in the plants response to the stresses presented in this study. The Esi3–9 in T. aestivum has a unique N terminal extension placing it into Group III, a new group for the Esi3/RCI2/PMP3 gene family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5311-8) contains supplementary material, which is available to authorized users.

Identification of nicotianamine synthase genes in Triticum monococcum and their expression under different Fe and Zn concentrations

In graminaceous plants, nicotianamine (NA) is an important component of metal acquisition. NA is synthesized from S-adenosyl-l-methionine (SAM) catalyzed by nicotianamine synthase (NAS). Here, eight Triticum monococcum NAS (TmNAS) genes were cloned and characterized. Amino acid sequence analysis showed that TmNAS genes had high sequence identity with those from Triticum aestivum, Zea mays, Oryza sativa and Hordeum vulgare. Phylogenetic analysis showed that NAS genes were classified into two distinct groups, e.g. group I and group II. Expression analysis demonstrated that two of the TmNAS genes in group II were highly expressed in shoot tissues, and the other six TmNAS genes in group I were expressed in root tissues. Further analysis indicated that root-specific TmNAS genes were up-regulated under conditions of Fe- or Zn-deficiency growth, while shoot-specific TmNAS genes were up-regulated under conditions of Fe- or Zn-sufficiency. These results help us understand the NAS genes in T. monococcum and provide novel genetic resources for improving Fe and Zn concentrations in common wheat.

Asymmetrical changes of gene expression, small RNAs and chromatin in two resynthesized wheat allotetraploids

Polyploidy occurs in some animals and all flowering plants, including important crops such as wheat. The consequences of polyploidy in crops remain elusive, partly because their progenitors are unknown. Using two resynthesized wheat allotetraploids S^l S^l AA and AADD with known diploid progenitors, we analyzed mRNA and small RNA transcriptomes in the endosperm, compared transcriptomes between endosperm and root in AADD, and examined chromatin changes in the allotetraploids. In the endosperm, there were more non-additively expressed genes in S^l S^l AA than in AADD. In AADD, non-additively expressed genes were developmentally regulated, and the majority (62-70%) were repressed. The repressed genes in AADD included a group of histone methyltransferase gene homologs, which correlated with reduced histone H3K9me2 levels and activation of various transposable elements in AADD. In S^l S^l AA, there was a tendency for expression dominance of S^l over A homoeologs, but the histone methyltransferase gene homologs were additively expressed, correlating with insignificant changes in histone H3K9me2 levels. Moreover, more 24-nucleotide small inferring RNAs (siRNAs) in the A subgenome were disrupted in AADD than in S^l S^l AA, which were associated with expression changes of siRNA-associated genes. Our results indicate that asymmetrical changes in siRNAs, chromatin modifications, transposons and gene expression coincide with unstable AADD genomes and stable S^l S^l AA genomes, which could help explain the evolutionary trajectories of wheat allotetraploids formed by different progenitors.

Genomic Imprinting Was Evolutionarily Conserved during Wheat Polyploidization

Genomic imprinting is an epigenetic phenomenon that causes genes to be differentially expressed depending on their parent of origin. To evaluate the evolutionary conservation of genomic imprinting and the effects of ploidy on this process, we investigated parent-of-origin-specific gene expression patterns in the endosperm of diploid (Aegilops spp), tetraploid, and hexaploid wheat (Triticum spp) at various stages of development via high-throughput transcriptome sequencing. We identified 91, 135, and 146 maternally or paternally expressed genes (MEGs or PEGs, respectively) in diploid, tetraploid, and hexaploid wheat, respectively, 52.7% of which exhibited dynamic expression patterns at different developmental stages. Gene Ontology enrichment analysis suggested that MEGs and PEGs were involved in metabolic processes and DNA-dependent transcription, respectively. Nearly half of the imprinted genes exhibited conserved expression patterns during wheat hexaploidization. In addition, 40% of the homoeolog pairs originating from whole-genome duplication were consistently maternally or paternally biased in the different subgenomes of hexaploid wheat. Furthermore, imprinted expression was found for 41.2% and 50.0% of homolog pairs that evolved by tandem duplication after genome duplication in tetraploid and hexaploid wheat, respectively. These results suggest that genomic imprinting was evolutionarily conserved between closely related Triticum and Aegilops species and in the face of polyploid hybridization between species in these genera.

A comparative analysis of nonhost resistance across the two Triticeae crop species wheat and barley

BACKGROUND

Nonhost resistance (NHR) protects plants against a vast number of non-adapted pathogens which implicates a potential exploitation as source for novel disease resistance strategies. Aiming at a fundamental understanding of NHR a global analysis of transcriptome reprogramming in the economically important Triticeae cereals wheat and barley, comparing host and nonhost interactions in three major fungal pathosystems responsible for powdery mildew (Blumeria graminis ff. ssp.), cereal blast (Magnaporthe sp.) and leaf rust (Puccinia sp.) diseases, was performed.

RESULTS

In each pathosystem a significant transcriptome reprogramming by adapted- or non-adapted pathogen isolates was observed, with considerable overlap between Blumeria, Magnaporthe and Puccinia. Small subsets of these general pathogen-regulated genes were identified as differentially regulated between host and corresponding nonhost interactions, indicating a fine-tuning of the general pathogen response during the course of co-evolution. Additionally, the host- or nonhost-related responses were rather specific for each pair of adapted and non-adapted isolates, indicating that the nonhost resistance-related responses were to a great extent pathosystem-specific. This pathosystem-specific reprogramming may reflect different resistance mechanisms operating against non-adapted pathogens with different lifestyles, or equally, different co-option of the hosts by the adapted isolates to create an optimal environment for infection. To compare the transcriptional reprogramming between wheat and barley, putative orthologues were identified. Within the wheat and barley general pathogen-regulated genes, temporal expression profiles of orthologues looked similar, indicating conserved general responses in Triticeae against fungal attack. However, the comparison of orthologues differentially expressed between host and nonhost interactions revealed fewer commonalities between wheat and barley, but rather suggested different host or nonhost responses in the two cereal species.

CONCLUSIONS

Taken together, our results suggest independent co-evolutionary forces acting on host pathosystems mirrored by barley- or wheat-specific nonhost responses. As a result of evolutionary processes, at least for the pathosystems investigated, NHR appears to rely on rather specific plant responses.

Genome-wide identification, classification, evolutionary analysis and gene expression patterns of the protein kinase gene family in wheat and Aegilops tauschii

In this study we systematically identified and classified PKs in Triticum aestivum, Triticum urartu and Aegilops tauschii. Domain distribution and exon-intron structure analyses of PKs were performed, and we found conserved exon-intron structures within the exon phases in the kinase domain. Collinearity events were determined, and we identified various T. aestivum PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed that some PKs might participate in the signaling pathways of stress response and developmental processes. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of Fusarium graminearum to validate the prediction of microarray. The protein kinase (PK) gene superfamily is one of the largest families in plants and participates in various plant processes, including growth, development, and stress response. To better understand wheat PKs, we conducted genome-wide identification, classification, evolutionary analysis and expression profiles of wheat and Ae. tauschii PKs. We identified 3269, 1213 and 1448 typical PK genes in T. aestivum, T. urartu and Ae. tauschii, respectively, and classified them into major groups and subfamilies. Domain distributions and gene structures were analyzed and visualized. Some conserved intron-exon structures within the conserved kinase domain were found in T. aestivum, T. urartu and Ae. tauschii, as well as the primitive land plants Selaginella moellendorffii and Physcomitrella patens, revealing the important roles and conserved evolutionary history of these PKs. We analyzed the collinearity events of T. aestivum PKs and identified PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed tissue-specific and stress-specific expression profiles, hinting that some wheat PKs may regulate abiotic and biotic stress response signaling pathways. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of F. graminearum to validate the prediction of microarray. Our results will provide the foundational information for further studies on the molecular functions of wheat PKs.

Characterisation of the nicotianamine aminotransferase and deoxymugineic acid synthase genes essential to Strategy II iron uptake in bread wheat (Triticum aestivum L.)

Iron (Fe) uptake in graminaceous plant species occurs via the release and uptake of Fe-chelating compounds known as mugineic acid family phytosiderophores (MAs). In the MAs biosynthetic pathway, nicotianamine aminotransferase (NAAT) and deoxymugineic acid synthase (DMAS) enzymes catalyse the formation of 2'-deoxymugineic acid (DMA) from nicotianamine (NA). Here we describe the identification and characterisation of six TaNAAT and three TaDMAS1 genes in bread wheat (Triticum aestivum L.). The coding sequences of all six TaNAAT homeologs consist of seven exons with ≥88.0% nucleotide sequence identity and most sequence variation present in the first exon. The coding sequences of the three TaDMAS1 homeologs consist of three exons with ≥97.8% nucleotide sequence identity. Phylogenetic analysis revealed that the TaNAAT and TaDMAS1 proteins are most closely related to the HvNAAT and HvDMAS1 proteins of barley and that there are two distinct groups of TaNAAT proteins-TaNAAT1 and TaNAAT2 -that correspond to the HvNAATA and HvNAATB proteins, respectively. Quantitative reverse transcription-PCR analysis revealed that the TaNAAT2 genes are expressed at highest levels in anther tissues whilst the TaNAAT1 and TaDMAS1 genes are expressed at highest levels in root tissues of bread wheat. Furthermore, the TaNAAT1, TaNAAT2 and TaDMAS1 genes were differentially regulated by plant Fe status and their expression was significantly upregulated in root tissues from day five onwards during a seven-day Fe deficiency treatment. The identification and characterization of the TaNAAT1, TaNAAT2 and TaDMAS1 genes provides a valuable genetic resource for improving bread wheat growth on Fe deficient soils and enhancing grain Fe nutrition.

Identification and molecular characterization of the nicotianamine synthase gene family in bread wheat

Nicotianamine (NA) is a non-protein amino acid involved in fundamental aspects of metal uptake, transport and homeostasis in all plants and constitutes the biosynthetic precursor of mugineic acid family phytosiderophores (MAs) in graminaceous plant species. Nicotianamine synthase (NAS) genes, which encode enzymes that synthesize NA from S-adenosyl-L-methionine (SAM), are differentially regulated by iron (Fe) status in most plant species and plant genomes have been found to contain anywhere from 1 to 9 NAS genes. This study describes the identification of 21 NAS genes in the hexaploid bread wheat (Triticum aestivum L.) genome and their phylogenetic classification into two distinct clades. The TaNAS genes are highly expressed during germination, seedling growth and reproductive development. Fourteen of the clade I NAS genes were up-regulated in root tissues under conditions of Fe deficiency. Protein sequence analyses revealed the presence of endocytosis motifs in all of the wheat NAS proteins as well as chloroplast, mitochondrial and secretory transit peptide signals in four proteins. These results greatly expand our knowledge of NAS gene families in graminaceous plant species as well as the genetics underlying Fe nutrition in bread wheat.

Fast-forward genetics by radiation hybrids to saturate the locus regulating nuclear-cytoplasmic compatibility in Triticum

The nuclear-encoded species cytoplasm specific (scs) genes control nuclear-cytoplasmic compatibility in wheat (genus Triticum). Alloplasmic cells, which have nucleus and cytoplasm derived from different species, produce vigorous and vital organisms only when the correct version of scs is present in their nucleus. In this study, bulks of in vivo radiation hybrids segregating for the scs phenotype have been genotyped by sequencing with over 1.9 million markers. The high marker saturation obtained for a critical region of chromosome 1D allowed identification of 3318 reads that mapped in close proximity of the scs. A novel in silico approach was deployed to extend these short reads to sequences of up to 70 Kb in length and identify candidate open reading frames (ORFs). Markers were developed to anchor the short contigs containing ORFs to a radiation hybrid map of 650 individuals with resolution of 288 Kb. The region containing the scs locus was narrowed to a single Bacterial Artificial Chromosome (BAC) contig of Aegilops tauschii. Its sequencing and assembly by nano-mapping allowed rapid identification of a rhomboid gene as the only ORF existing within the refined scs locus. Resequencing of this gene from multiple germplasm sources identified a single nucleotide mutation, which gives rise to a functional amino acid change. Gene expression characterization revealed that an active copy of this rhomboid exists on all homoeologous chromosomes of wheat, and depending on the specific cytoplasm each copy is preferentially expressed. Therefore, a new methodology was applied to unique genetic stocks to rapidly identify a strong candidate gene for the control of nuclear-cytoplasmic compatibility in wheat.

Comprehensive analyses of the annexin gene family in wheat

Background

Annexins are an evolutionarily conserved multigene family of calcium-dependent phospholipid binding proteins that play important roles in stress resistance and plant development. They have been relatively well characterized in model plants Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), but nothing has been reported in hexaploid bread wheat (Triticum aestivum) and barely (Hordeum vulgare), which are the two most economically important plants.

Results

Based on available genomic and transcriptomic data, 25 and 11 putative annexin genes were found through in silico analysis in wheat and barley, respectively. Additionally, eight and 11 annexin genes were identified from the draft genome sequences of Triticum urartu and Aegilops tauschii, progenitor for the A and D genome of wheat, respectively. By phylogenetic analysis, annexins in these four species together with other monocots and eudicots were classified into six different orthologous groups. Pi values of each of Ann1–12 genes among T. aestivum, T. urartu, A. tauschii and H. vulgare species was very low, with the exception of Ann2 and Ann5 genes. Ann2 gene has been under positive selection, but Ann6 and Ann7 have been under purifying selection among the four species in their evolutionary histories. The nucleotide diversities of Ann1–12 genes in the four species were 0.52065, 0.59239, 0.60691 and 0.53421, respectively. No selective pressure was operated on annexin genes in the same species. Gene expression patterns obtained by real-time PCR and re-analyzing the public microarray data revealed differential temporal and spatial regulation of annexin genes in wheat under different abiotic stress conditions such as salinity, drought, cold and abscisic acid. Among those genes, TaAnn10 is specifically expressed in the anther but fails to be induced by low temperature in thermosensitive genic male sterile lines, suggesting that specific down-regulation of TaAnn10 is associated with conditional male sterility in wheat.

Conclusions

This study analyzed the size and composition of the annexin gene family in wheat and barley, and investigated differential tissue-specific and stress responsive expression profiles of the gene family in wheat. These results provided significant information for understanding the diverse roles of plant annexins and opened a new avenue for functional studies of cold induced male sterility in wheat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2750-y) contains supplementary material, which is available to authorized users.

Expansion Mechanisms and Evolutionary History on Genes Encoding DNA Glycosylases and Their Involvement in Stress and Hormone Signaling

DNA glycosylases catalyze the release of methylated bases. They play vital roles in the base excision repair pathway and might also function in DNA demethylation. At least three families of DNA glycosylases have been identified, which included 3′-methyladenine DNA glycosylase (MDG) I, MDG II, and HhH-GPD (Helix–hairpin–Helix and Glycine/Proline/aspartate (D)). However, little is known on their genome-wide identification, expansion, and evolutionary history as well as their expression profiling and biological functions. In this study, we have genome-widely identified and evolutionarily characterized these family members. Generally, a genome encodes only one MDG II gene in most of organisms. No MDG I or MDG II gene was detected in green algae. However, HhH-GPD genes were detectable in all available organisms. The ancestor species contain small size of MDG I and HhH-GPD families. These two families were mainly expanded through the whole-genome duplication and segmental duplication. They were evolutionarily conserved and were generally under purifying selection. However, we have detected recent positive selection among the Oryza genus, which might play roles in species divergence. Further investigation showed that expression divergence played important roles in gene survival after expansion. All of these family genes were expressed in most of developmental stages and tissues in rice plants. High ratios of family genes were downregulated by drought and fungus pathogen as well as abscisic acid (ABA) and jasmonic acid (JA) treatments, suggesting a negative regulation in response to drought stress and pathogen infection through ABA- and/or JA-dependent hormone signaling pathway.

Genetic Diversity and Genome Wide Association Study of β-Glucan Content in Tetraploid Wheat Grains

Non-starch polysaccharides (NSPs) have many health benefits, including immunomodulatory activity, lowering serum cholesterol, a faecal bulking effect, enhanced absorption of certain minerals, prebiotic effects and the amelioration of type II diabetes. The principal components of the NSP in cereal grains are (1,3;1,4)-β-glucans and arabinoxylans. Although (1,3;1,4)-β-glucan (hereafter called β-glucan) is not the most representative component of wheat cell walls, it is one of the most important types of soluble fibre in terms of its proven beneficial effects on human health. In the present work we explored the genetic variability of β-glucan content in grains from a tetraploid wheat collection that had been genotyped with a 90k-iSelect array, and combined this data to carry out an association analysis. The β-glucan content, expressed as a percentage w/w of grain dry weight, ranged from 0.18% to 0.89% across the collection. Our analysis identified seven genomic regions associated with β-glucan, located on chromosomes 1A, 2A (two), 2B, 5B and 7A (two), confirming the quantitative nature of this trait. Analysis of marker trait associations (MTAs) in syntenic regions of several grass species revealed putative candidate genes that might influence β-glucan levels in the endosperm, possibly via their participation in carbon partitioning. These include the glycosyl hydrolases endo-β-(1,4)-glucanase (cellulase), β-amylase, (1,4)-β-xylan endohydrolase, xylanase inhibitor protein I, isoamylase and the glycosyl transferase starch synthase II.

Maize maintains growth in response to decreased nitrate supply through a highly dynamic and developmental stage-specific transcriptional response

Elucidation of the gene networks underlying the response to N supply and demand will facilitate the improvement of the N uptake efficiency of plants. We undertook a transcriptomic analysis of maize to identify genes responding to both a non-growth-limiting decrease in NO3- provision and to development-based N demand changes at seven representative points across the life cycle. Gene co-expression networks were derived by cluster analysis of the transcript profiles. The majority of NO3--responsive transcription occurred at 11 (D11), 18 (D18) and 29 (D29) days after emergence, with differential expression predominating in the root at D11 and D29 and in the leaf at D18. A cluster of 98 probe sets was identified, the expression pattern of which is similar to that of the high-affinity NO3- transporter (NRT2) genes across the life cycle. The cluster is enriched with genes encoding enzymes and proteins of lipid metabolism and transport, respectively. These are candidate genes for the response of maize to N supply and demand. Only a few patterns of differential gene expression were observed over the entire life cycle; however, the composition of the classes of the genes differentially regulated at individual time points was unique, suggesting tightly controlled regulation of NO3--responsive gene expression.

Genome-Wide Comparative Analysis of Flowering-Related Genes in Arabidopsis, Wheat, and Barley

Early flowering is an important trait influencing grain yield and quality in wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) in short-season cropping regions. However, due to large and complex genomes of these species, direct identification of flowering genes and their molecular characterization remain challenging. Here, we used a bioinformatic approach to predict flowering-related genes in wheat and barley from 190 known Arabidopsis (Arabidopsis thaliana (L.) Heynh.) flowering genes. We identified 900 and 275 putative orthologs in wheat and barley, respectively. The annotated flowering-related genes were clustered into 144 orthologous groups with one-to-one, one-to-many, many-to-one, and many-to-many orthology relationships. Our approach was further validated by domain and phylogenetic analyses of flowering-related proteins and comparative analysis of publicly available microarray data sets for in silico expression profiling of flowering-related genes in 13 different developmental stages of wheat and barley. These further analyses showed that orthologous gene pairs in three critical flowering gene families (PEBP, MADS, and BBX) exhibited similar expression patterns among 13 developmental stages in wheat and barley, suggesting similar functions among the orthologous genes with sequence and expression similarities. The predicted candidate flowering genes can be confirmed and incorporated into molecular breeding for early flowering wheat and barley in short-season cropping regions.

Evolutionary, Comparative and Functional Analyses of the Brassinosteroid Receptor Gene, BRI1, in Wheat and Its Relation to Other Plant Genomes

Brassinosteroids (BRs) are plant hormones, fundamental for the growth and development of plants. A trans-membrane protein receptor kinase, Brassinosteroid-Insensitive 1 (BRI1), is known to interact with BRs and be directly involved in plant development. This study investigates the structural organization of BRI1 orthologs in several taxa, with a specific interest in Triticum aestivum. True orthologs of Arabidopsis thaliana BRI1 (AtBRI1) from seven-plant species showed sequence identity ranging from 54% to 95% at the protein level. All gene sequences lacked introns, leading to speculation that post-transcriptional processing in TaBRI1 is similar to AtBRI1. Based on in silico analysis, a single copy of BRI1 was present in each of the three wheat genomes on the long arm of chromosome 3. Domain structure of BRI1 orthologs among different taxa showed multiple leucine rich repeats (LRRs), an island domain (ID), a juxtamembrane/transmembrane domain (JTMD), a catalytic kinase domain (KD), C and N-Terminal domains. The KD showed the highest level of conservation while the LRRs and JTMD were most variable. Phosphorylation of residues in the juxtamembrane domain, known to be involved in the activation of the KD, is conserved in TaBRI1. While TaBRI1 has well-defined differences in the ID and LRR domains, many residues involved in ligand binding are conserved. The activation loop present in the KD showed 100% conservation in all taxa. Despite residue differences, hydrophobicity was conserved in the BR binding pocket across taxa, suggesting that function may not differ as drastically as residue identity may suggest. Predicted 3D structure of AtBRI1 and TaBRI1 showed a conserved super helical assembly, a feature essential in protein-protein interactions. An unrooted phylogram showed TaBRI1 in the monocot clade to be distinct from that of dicots. New insight in the structure and functions of BRI1 may help in targeting BR pathway for crop improvement.

A wheat ABC transporter contributes to both grain formation and mycotoxin tolerance

The mycotoxin deoxynivalenol (DON) acts as a disease virulence factor for Fusarium fungi, and tolerance of DON enhances wheat resistance to Fusarium head blight (FHB) disease. Two variants of an ATP-binding cassette (ABC) family C transporter gene were cloned from DON-treated wheat mRNA, namely TaABCC3.1 and TaABCC3.2. These represent two of three putative genes identified on chromosomes 3A, 3B, and 3D of the wheat genome sequence. Variant TaABCC3.1 represents the DON-responsive transcript previously associated with DON resistance in wheat. PCR-based mapping and in silico sequence analyses located TaABCC3.1 to the short arm of wheat chromosome 3B (not within the FHB resistance quantitative trait locus Fhb1). In silico analyses of microarray data indicated that TaABCC3 genes are expressed in reproductive tissue and roots, and in response to the DON producer Fusarium graminearum. Gene expression studies showed that TaABCC3.1 is activated as part of the early host response to DON and in response to the FHB defence hormone jasmonic acid. Virus-induced gene silencing (VIGS) confirmed that TaABCC3 genes contributed to DON tolerance. VIGS was performed using two independent viral construct applications: one specifically targeted TaABCC3.1 for silencing, while the other targeted this gene and the chromosome 3A homeologue. In both instances, VIGS resulted in more toxin-induced discoloration of spikelets, compared with the DON effects in non-silenced spikelets at 14 d after toxin treatment (≥2.2-fold increase, P<0.05). Silencing by both VIGS constructs enhanced head ripening, and especially so in DON-treated heads. VIGS of TaABCC3 genes also reduced the grain number by more than 28% (P<0.05), both with and without DON treatment, and the effects were greater for the construct that targeted the two homeologues. Hence, DON-responsive TaABCC3 genes warrant further study to determine their potential as disease resistance breeding targets and their function in grain formation and ripening.

Molecular basis of adaptation to high soil boron in wheat landraces and elite cultivars

Environmental constraints severely restrict crop yields in most production environments, and expanding the use of variation will underpin future progress in breeding. In semi-arid environments boron toxicity constrains productivity, and genetic improvement is the only effective strategy for addressing the problem. Wheat breeders have sought and used available genetic diversity from landraces to maintain yield in these environments; however, the identity of the genes at the major tolerance loci was unknown. Here we describe the identification of near-identical, root-specific boron transporter genes underlying the two major-effect quantitative trait loci for boron tolerance in wheat, Bo1 and Bo4 (ref. 2). We show that tolerance to a high concentration of boron is associated with multiple genomic changes including tetraploid introgression, dispersed gene duplication, and variation in gene structure and transcript level. An allelic series was identified from a panel of bread and durum wheat cultivars and landraces originating from diverse agronomic zones. Our results demonstrate that, during selection, breeders have matched functionally different boron tolerance alleles to specific environments. The characterization of boron tolerance in wheat illustrates the power of the new wheat genomic resources to define key adaptive processes that have underpinned crop improvement.

Evolution of the BBAA component of bread wheat during its history at the allohexaploid level

Subgenome integrity in bread wheat (Triticum aestivum; BBAADD) makes possible the extraction of its BBAA component to restitute a novel plant type. The availability of such a ploidy-reversed wheat (extracted tetraploid wheat [ETW]) provides a unique opportunity to address whether and to what extent the BBAA component of bread wheat has been modified in phenotype, karyotype, and gene expression during its evolutionary history at the allohexaploid level. We report here that ETW was anomalous in multiple phenotypic traits but maintained a stable karyotype. Microarray-based transcriptome profiling identified a large number of differentially expressed genes between ETW and natural tetraploid wheat (Triticum turgidum), and the ETW-downregulated genes were enriched for distinct Gene Ontology categories. Quantitative RT-PCR analysis showed that gene expression differences between ETW and a set of diverse durum wheat (T. turgidum subsp durum) cultivars were distinct from those characterizing tetraploid cultivars per se. Pyrosequencing revealed that the expression alterations may occur to either only one or both of the B and A homoeolog transcripts in ETW. A majority of the genes showed additive expression in a resynthesized allohexaploid wheat. Analysis of a synthetic allohexaploid wheat and diverse bread wheat cultivars revealed the rapid occurrence of expression changes to the BBAA subgenomes subsequent to allohexaploidization and their evolutionary persistence.

Characterization of the caleosin gene family in the Triticeae

Background

The caleosin genes encode proteins with a single conserved EF hand calcium-binding domain and comprise small gene families found in a wide range of plant species. Some members of the gene family have been shown to be upregulated by environmental stresses including low water availability and high salinity. Caleosin 3 from wheat has been shown to interact with the α-subunit of the heterotrimeric G proteins, and to act as a GTPase activating protein (GAP). This study characterizes the size and diversity of the gene family in wheat and related species and characterizes the differential tissue-specific expression of members of the gene family.

Results

A total of 34 gene family members that belong to eleven paralogous groups of caleosins were identified in the hexaploid bread wheat, T. aestivum. Each group was represented by three homeologous copies of the gene located on corresponding homeologous chromosomes, except the caleosin 10, which has four gene copies. Ten gene family members were identified in diploid barley, Hordeum vulgare, and in rye, Secale cereale, seven in Brachypodium distachyon, and six in rice, Oryza sativa. The analysis of gene expression was assayed in triticale and rye by RNA-Seq analysis of 454 sequence sets and members of the gene family were found to have diverse patterns of gene expression in the different tissues that were sampled in rye and in triticale, the hybrid hexaploid species derived from wheat and rye. Expression of the gene family in wheat and barley was also previously determined by microarray analysis, and changes in expression during development and in response to environmental stresses are presented.

Conclusions

The caleosin gene family had a greater degree of expansion in the Triticeae than in the other monocot species, Brachypodium and rice. The prior implication of one member of the gene family in the stress response and heterotrimeric G protein signaling, points to the potential importance of the caleosin gene family. The complexity of the family and differential expression in various tissues and under conditions of abiotic stress suggests the possibility that caleosin family members may play diverse roles in signaling and development that warrants further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-239) contains supplementary material, which is available to authorized users.

A detailed gene expression study of the Miscanthus genus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes

BACKGROUND

The Miscanthus genus of perennial C4 grasses contains promising biofuel crops for temperate climates. However, few genomic resources exist for Miscanthus, which limits understanding of its interesting biology and future genetic improvement. A comprehensive catalog of expressed sequences were generated from a variety of Miscanthus species and tissue types, with an emphasis on characterizing gene expression changes in spring compared to fall rhizomes.

RESULTS

Illumina short read sequencing technology was used to produce transcriptome sequences from different tissues and organs during distinct developmental stages for multiple Miscanthus species, including Miscanthus sinensis, Miscanthus sacchariflorus, and their interspecific hybrid Miscanthus × giganteus. More than fifty billion base-pairs of Miscanthus transcript sequence were produced. Overall, 26,230 Sorghum gene models (i.e., ~ 96% of predicted Sorghum genes) had at least five Miscanthus reads mapped to them, suggesting that a large portion of the Miscanthus transcriptome is represented in this dataset. The Miscanthus × giganteus data was used to identify genes preferentially expressed in a single tissue, such as the spring rhizome, using Sorghum bicolor as a reference. Quantitative real-time PCR was used to verify examples of preferential expression predicted via RNA-Seq. Contiguous consensus transcript sequences were assembled for each species and annotated using InterProScan. Sequences from the assembled transcriptome were used to amplify genomic segments from a doubled haploid Miscanthus sinensis and from Miscanthus × giganteus to further disentangle the allelic and paralogous variations in genes.

CONCLUSIONS

This large expressed sequence tag collection creates a valuable resource for the study of Miscanthus biology by providing detailed gene sequence information and tissue preferred expression patterns. We have successfully generated a database of transcriptome assemblies and demonstrated its use in the study of genes of interest. Analysis of gene expression profiles revealed biological pathways that exhibit altered regulation in spring compared to fall rhizomes, which are consistent with their different physiological functions. The expression profiles of the subterranean rhizome provides a better understanding of the biological activities of the underground stem structures that are essentials for perenniality and the storage or remobilization of carbon and nutrient resources.

Developmental transcriptional profiling reveals key insights into Triticeae reproductive development

Despite their importance, there remains a paucity of large-scale gene expression-based studies of reproductive development in species belonging to the Triticeae. As a first step to address this deficiency, a gene expression atlas of triticale reproductive development was generated using the 55K Affymetrix GeneChip(®) wheat genome array. The global transcriptional profiles of the anther/pollen, ovary and stigma were analyzed at concurrent developmental stages, and co-expressed as well as preferentially expressed genes were identified. Data analysis revealed both novel and conserved regulatory factors underlying Triticeae floral development and function. This comprehensive resource rests upon detailed gene annotations, and the expression profiles are readily accessible via a web browser.

Prospects for advancing defense to cereal rusts through genetical genomics

Rusts are one of the most severe threats to cereal crops because new pathogen races emerge regularly, resulting in infestations that lead to large yield losses. In 1999, a new race of stem rust, Puccinia graminis f. sp. tritici (Pgt TTKSK or Ug99), was discovered in Uganda. Most of the wheat and barley cultivars grown currently worldwide are susceptible to this new race. Pgt TTKSK has already spread northward into Iran and will likely spread eastward throughout the Indian subcontinent in the near future. This scenario is not unique to stem rust; new races of leaf rust (Puccinia triticina) and stripe rust (Puccinia striiformis) have also emerged recently. One strategy for countering the persistent adaptability of these pathogens is to stack complete- and partial-resistance genes, which requires significant breeding efforts in order to reduce deleterious effects of linkage drag. These varied resistance combinations are typically more difficult for the pathogen to defeat, since they would be predicted to apply lower selection pressure. Genetical genomics or expression Quantitative Trait Locus (eQTL) analysis enables the identification of regulatory loci that control the expression of many to hundreds of genes. Integrated deployment of these technologies coupled with efficient phenotyping offers significant potential to elucidate the regulatory nodes in genetic networks that orchestrate host defense responses. The focus of this review will be to present advances in genetical genomic experimental designs and analysis, particularly as they apply to the prospects for discovering partial disease resistance alleles in cereals.

Conserved expression of natural antisense transcripts in mammals

Background

Recent studies had found thousands of natural antisense transcripts originating from the same genomic loci of protein coding genes but from the opposite strand. It is unclear whether the majority of antisense transcripts are functional or merely transcriptional noise.

Results

Using the Affymetrix Exon array with a modified cDNA synthesis protocol that enables genome-wide detection of antisense transcription, we conducted large-scale expression analysis of antisense transcripts in nine corresponding tissues from human, mouse and rat. We detected thousands of antisense transcripts, some of which show tissue-specific expression that could be subjected to further study for their potential function in the corresponding tissues/organs. The expression patterns of many antisense transcripts are conserved across species, suggesting selective pressure on these transcripts. When compared to protein-coding genes, antisense transcripts show a lesser degree of expression conservation. We also found a positive correlation between the sense and antisense expression across tissues.

Conclusion

Our results suggest that natural antisense transcripts are subjected to selective pressure but to a lesser degree compared to sense transcripts in mammals.

Tissue-specific transcriptomic profiling of Sorghum propinquum using a rice genome array

Sorghum (Sorghum bicolor) is one of the world's most important cereal crops. S. propinquum is a perennial wild relative of S. bicolor with well-developed rhizomes. Functional genomics analysis of S. propinquum, especially with respect to molecular mechanisms related to rhizome growth and development, can contribute to the development of more sustainable grain, forage, and bioenergy cropping systems. In this study, we used a whole rice genome oligonucleotide microarray to obtain tissue-specific gene expression profiles of S. propinquum with special emphasis on rhizome development. A total of 548 tissue-enriched genes were detected, including 31 and 114 unique genes that were expressed predominantly in the rhizome tips (RT) and internodes (RI), respectively. Further GO analysis indicated that the functions of these tissue-enriched genes corresponded to their characteristic biological processes. A few distinct cis-elements, including ABA-responsive RY repeat CATGCA, sugar-repressive TTATCC, and GA-responsive TAACAA, were found to be prevalent in RT-enriched genes, implying an important role in rhizome growth and development. Comprehensive comparative analysis of these rhizome-enriched genes and rhizome-specific genes previously identified in Oryza longistaminata and S. propinquum indicated that phytohormones, including ABA, GA, and SA, are key regulators of gene expression during rhizome development. Co-localization of rhizome-enriched genes with rhizome-related QTLs in rice and sorghum generated functional candidates for future cloning of genes associated with rhizome growth and development.

Clusters of genes encoding fructan biosynthesizing enzymes in wheat and barley

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

TaWIR1 contributes to post-penetration resistance to Magnaporthe oryzae, but not Blumeria graminis f. sp. tritici, in wheat

Members of the Wheat-Induced Resistance 1 (TaWIR1) gene family are highly induced in response to a wide range of pathogens. Homologues have been identified in barley, but not in Brachypodium, whereas, in rice, only distant WIR1 candidates are known. Phylogenetic analysis placed TaWIR1a and TaWIR1b within a distinct clade of wheat transcripts, whereas TaWIR1c clustered with HvWIR1 genes. Transcripts of all three TaWIR1 genes were strongly induced by a wheat-adapted isolate of Magnaporthe oryzae. Virus-induced gene silencing of the TaWIR1 gene family had no effect on the initial penetration of epidermal cells by M. oryzae. However, following the establishment of an infection site, the fungus was able to grow more extensively within the leaf tissue, relative to control leaves, indicating a role for the TaWIR1 gene family in the cell-to-cell movement of M. oryzae. In contrast, the silencing of TaWIR1 transcripts had no effect on epidermal cell penetration by a wheat-adapted isolate of Blumeria graminis, or on the subsequent growth of hyphae. Differential transcription of TaWIR1 genes was also seen in epidermal peels, relative to the remaining leaf tissue, following inoculation with M. oryzae.

HvHMA2, a P(1B)-ATPase from barley, is highly conserved among cereals and functions in Zn and Cd transport

Manipulation of crops to improve their nutritional value (biofortification) and optimisation of plants for removal of toxic metals from contaminated soils (phytoremediation) are major goals. Identification of membrane transporters with roles in zinc and cadmium transport would be useful for both aspects. The P_1B-ATPases play important roles in heavy metal allocation and detoxification in Arabidopsis and it is now important to elucidate their roles in monocots. We identified nine P_1B-ATPases in barley and this study focuses on the functional characterization of HvHMA2, providing evidence for its role in heavy metal transport. HvHMA2 was cloned using information from EST analysis and 5′ RACE. It possesses the conserved aspartate that is phosphorylated during the reaction cycle of P-type pumps and has motifs and key residues characteristic of P_1B-ATPases, falling into the P_1B-2 subclass. Homologous sequences occur in three major sub-families of the Poaceae (Gramineae). Heterologous expression in Saccharomyces cerevisiae demonstrates that HvHMA2 functions as a Zn and Cd pump. Mutagenesis studies show that proposed cation coordination sites of the P_1B-2 pumps are crucial for the metal responses conferred by HvHMA2 in yeast. HvHMA2 expression suppresses the Zn-deficient phenotype of the Arabidopsis hma2hma4 mutant indicating that HvHMA2 functions as a Zn pump in planta and could play a role in root to shoot Zn transport. When expressed in Arabidopsis, HvHMA2 localises predominantly to the plasma membrane.

Global transgenerational gene expression dynamics in two newly synthesized allohexaploid wheat (Triticum aestivum) lines

Background

Alteration in gene expression resulting from allopolyploidization is a prominent feature in plants, but its spectrum and extent are not fully known. Common wheat (Triticum aestivum) was formed via allohexaploidization about 10,000 years ago, and became the most important crop plant. To gain further insights into the genome-wide transcriptional dynamics associated with the onset of common wheat formation, we conducted microarray-based genome-wide gene expression analysis on two newly synthesized allohexaploid wheat lines with chromosomal stability and a genome constitution analogous to that of the present-day common wheat.

Results

Multi-color GISH (genomic in situ hybridization) was used to identify individual plants from two nascent allohexaploid wheat lines between Triticum turgidum (2n = 4x = 28; genome BBAA) and Aegilops tauschii (2n = 2x = 14; genome DD), which had a stable chromosomal constitution analogous to that of common wheat (2n = 6x = 42; genome BBAADD). Genome-wide analysis of gene expression was performed for these allohexaploid lines along with their parental plants from T. turgidum and Ae. tauschii, using the Affymetrix Gene Chip Wheat Genome-Array. Comparison with the parental plants coupled with inclusion of empirical mid-parent values (MPVs) revealed that whereas the great majority of genes showed the expected parental additivity, two major patterns of alteration in gene expression in the allohexaploid lines were identified: parental dominance expression and non-additive expression. Genes involved in each of the two altered expression patterns could be classified into three distinct groups, stochastic, heritable and persistent, based on their transgenerational heritability and inter-line conservation. Strikingly, whereas both altered patterns of gene expression showed a propensity of inheritance, identity of the involved genes was highly stochastic, consistent with the involvement of diverse Gene Ontology (GO) terms. Nonetheless, those genes showing non-additive expression exhibited a significant enrichment for vesicle-function.

Conclusions

Our results show that two patterns of global alteration in gene expression are conditioned by allohexaploidization in wheat, that is, parental dominance expression and non-additive expression. Both altered patterns of gene expression but not the identity of the genes involved are likely to play functional roles in stabilization and establishment of the newly formed allohexaploid plants, and hence, relevant to speciation and evolution of T. aestivum.

Transcriptome analysis of H2O2-treated wheat seedlings reveals a H2O2-responsive fatty acid desaturase gene participating in powdery mildew resistance

Hydrogen peroxide (H(2)O(2)) plays important roles in plant biotic and abiotic stress responses. However, the effect of H(2)O(2) stress on the bread wheat transcriptome is still lacking. To investigate the cellular and metabolic responses triggered by H(2)O(2), we performed an mRNA tag analysis of wheat seedlings under 10 mM H(2)O(2) treatment for 6 hour in one powdery mildew (PM) resistant (PmA) and two susceptible (Cha and Han) lines. In total, 6,156, 6,875 and 3,276 transcripts were found to be differentially expressed in PmA, Han and Cha respectively. Among them, 260 genes exhibited consistent expression patterns in all three wheat lines and may represent a subset of basal H(2)O(2) responsive genes that were associated with cell defense, signal transduction, photosynthesis, carbohydrate metabolism, lipid metabolism, redox homeostasis, and transport. Among genes specific to PmA, 'transport' activity was significantly enriched in Gene Ontology analysis. MapMan classification showed that, while both up- and down- regulations were observed for auxin, abscisic acid, and brassinolides signaling genes, the jasmonic acid and ethylene signaling pathway genes were all up-regulated, suggesting H(2)O(2)-enhanced JA/Et functions in PmA. To further study whether any of these genes were involved in wheat PM response, 19 H(2)O(2)-responsive putative defense related genes were assayed in wheat seedlings infected with Blumeria graminis f. sp. tritici (Bgt). Eight of these genes were found to be co-regulated by H(2)O(2) and Bgt, among which a fatty acid desaturase gene TaFAD was then confirmed by virus induced gene silencing (VIGS) to be required for the PM resistance. Together, our data presents the first global picture of the wheat transcriptome under H(2)O(2) stress and uncovers potential links between H(2)O(2) and Bgt responses, hence providing important candidate genes for the PM resistance in wheat.

Screening and identification of seed-specific genes using digital differential display tools combined with microarray data from common wheat

Background

Wheat is one of the most important cereal crops for human beings, with seeds being the tissue of highly economic value. Various morphogenetic and metabolic processes are exclusively associated with seed maturation. The goal of this study was to screen and identify genes specifically expressed in the developing seed of wheat with an integrative utilization of digital differential display (DDD) and available online microarray databases.

Results

A total of 201 unigenes were identified as the results of DDD screening and microarray database searching. The expressions of 6 of these were shown to be seed-specific by qRT-PCR analysis. Further GO enrichment analysis indicated that seed-specific genes were mainly associated with defense response, response to stress, multi-organism process, pathogenesis, extracellular region, nutrient reservoir activity, enzyme inhibitor activity, antioxidant activity and oxidoreductase activity. A comparison of this set of genes with the rice (Oryza sativa) genome was also performed and approximately three-fifths of them have rice counterparts. Between the counterparts, around 63% showed similar expression patterns according to the microarray data.

Conclusions

In conclusion, the DDD screening combined with microarray data analysis is an effective strategy for the identification of seed-specific expressed genes in wheat. These seed-specific genes screened during this study will provide valuable information for further studies about the functions of these genes in wheat.

Modeling and phylogeny analysis of bread wheat MnSOD

Superoxide dismutase (SOD) acts as first line of defense against oxidative and genetic stress. Manganese superoxide dismutase (MnSOD), found in mitochondria or peroxisomes, contains Mn(III) at the active site. Therefore, it is of interest to study MnSOD from bread wheat (a grain crop). However, a structure model is not yet solved for bread wheat MnSOD. Hence, we describe the structure model of bread wheat MnSOD developed using homology model. The model provides molecular insight to metal binding molecular function towards the understanding of oxidative stress resistance in plants. The distinction of bread wheat (a monocot) MnSOD from dicots is also shown using phylogenetic analysis.

A guideline to family-wide comparative state-of-the-art quantitative RT-PCR analysis exemplified with a Brassicaceae cross-species seed germination case study

Comparative biology includes the comparison of transcriptome and quantitative real-time RT-PCR (qRT-PCR) data sets in a range of species to detect evolutionarily conserved and divergent processes. Transcript abundance analysis of target genes by qRT-PCR requires a highly accurate and robust workflow. This includes reference genes with high expression stability (i.e., low intersample transcript abundance variation) for correct target gene normalization. Cross-species qRT-PCR for proper comparative transcript quantification requires reference genes suitable for different species. We addressed this issue using tissue-specific transcriptome data sets of germinating Lepidium sativum seeds to identify new candidate reference genes. We investigated their expression stability in germinating seeds of L. sativum and Arabidopsis thaliana by qRT-PCR, combined with in silico analysis of Arabidopsis and Brassica napus microarray data sets. This revealed that reference gene expression stability is higher for a given developmental process between distinct species than for distinct developmental processes within a given single species. The identified superior cross-species reference genes may be used for family-wide comparative qRT-PCR analysis of Brassicaceae seed germination. Furthermore, using germinating seeds, we exemplify optimization of the qRT-PCR workflow for challenging tissues regarding RNA quality, transcript stability, and tissue abundance. Our work therefore can serve as a guideline for moving beyond Arabidopsis by establishing high-quality cross-species qRT-PCR.

Fusarium graminearum from expression analysis to functional assays

Fusarium graminearum, the causal agent of head blight of wheat, was the third filamentous fungus to have a completed genome sequence. Since the release of the genome sequence in 2003, F. graminearum has become a model for studies of genomics and transcriptomics, mycotoxins, fungal population genetics, gene function, and sexual development. Herein we present the techniques we have used in our laboratory to perform expression analyses on life cycle stages of F. graminearum and techniques to functionally characterize those genes identified as potentially interesting.

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.

Vitamin E biosynthesis: functional characterization of the monocot homogentisate geranylgeranyl transferase

The biosynthesis of the tocotrienol and tocopherol forms of vitamin E is initiated by prenylation of homogentisate. Geranylgeranyl diphosphate (GGDP) is the prenyl donor for tocotrienol synthesis, whereas phytyl diphosphate (PDP) is the prenyl donor for tocopherol synthesis. We have previously shown that tocotrienol synthesis is initiated in monocot seeds by homogentisate geranylgeranyl transferase (HGGT). This enzyme is related to homogentisate phytyltransferase (HPT), which catalyzes the prenylation step in tocopherol synthesis. Here we show that monocot HGGT is localized in the plastid and expressed primarily in seed endosperm. Despite the close structural relationship of monocot HGGT and HPT, these enzymes were found to have distinct substrate specificities. Barley (Hordeum vulgare cv. Morex) HGGT expressed in insect cells was six times more active with GGDP than with PDP, whereas the Arabidopsis HPT was nine times more active with PDP than with GGDP. However, only small differences were detected in the apparent Km values of barley HGGT for GGDP and PDP. Consistent with its in vitro substrate properties, barley HGGT generated a mixture of tocotrienols and tocopherols when expressed in the vitamin E-null vte2-1 mutant lacking a functional HPT. Relative levels of tocotrienols and tocopherols produced in vte2-1 differed between organs and growth stages, reflective of the composition of plastidic pools of GGDP and PDP. In addition, HGGT was able to functionally substitute for HPT to rescue vte2-1-associated phenotypes, including reduced seed viability and increased fatty acid oxidation of seed lipids. Overall, we show that monocot HGGT is biochemically distinct from HPT, but can replace HPT in important vitamin E-related physiological processes.

Lowly expressed genes in Arabidopsis thaliana bear the signature of possible pseudogenization by promoter degradation

Pseudogenes are defined as nonfunctional DNA sequences with homology to functional protein-coding genes, and they typically contain nonfunctional mutations within the presumptive coding region. In theory, pseudogenes can also be caused by mutations in upstream regulatory regions, appearing as open reading frames with attenuated expression. In this study, we identified 1,939 annotated protein-coding genes with little evidence of expression in Arabidopsis thaliana and characterized their molecular evolutionary characteristics. On average, this set of genes was shorter than expressed genes and evolved with a 2-fold higher rate of nonsynonymous substitutions. The divergence of upstream sequences, based on ortholog comparisons to A. lyrata, was also higher than expressed genes, suggesting that these lowly expressed genes could be examples of pseudogenization by promoter disablement, often due to transposable element insertion. We complemented our empirical study by extending the models of Force et al. (Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J. 1999. Preservation of duplicate genes by complementary, degenerative mutations. Genetics 151:1531-1545.) to derive the probability of promoter disablements after gene duplication.

A computational-based update on microRNAs and their targets in barley (Hordeum vulgare L.)

Background

Many plant species have been investigated in the last years for the identification and characterization of the corresponding miRNAs, nevertheless extensive studies are not yet available on barley (at the time of this writing). To extend and to update information on miRNAs and their targets in barley and to identify candidate polymorphisms at miRNA target sites, the features of previously known plant miRNAs have been used to systematically search for barley miRNA homologues and targets in the publicly available ESTs database. Matching sequences have then been related to Unigene clusters on which most of this study was based.

Results

One hundred-fifty-six microRNA mature sequences belonging to 50 miRNA families have been found to significantly match at least one EST sequence in barley. As expected on the basis of phylogenetic relations, miRNAs putatively orthologous to those of Triticum are significantly over-represented inside the set of identified barley microRNA mature sequences. Many previously known and several putatively new miRNA/target pairs have been identified. When the predicted microRNA targets were grouped into functional categories, biological processes previously known to be regulated by miRNAs, such as development and response to biotic and abiotic stress, have been highlighted and most of the target molecular functions were related to transcription regulation. Candidate microRNA coding genes have been reported and genetic variation (SNPs/indels) both in functional regions of putative miRNAs (mature sequence) and at miRNA target sites has been found.

Conclusions

This study has provided an update of the information on barley miRNAs and their targets representing a foundation for future studies. Many of previously known plant microRNAs have homologues in barley with expected important roles during development, nutrient deprivation, biotic and abiotic stress response and other important physiological processes. Putative polymorphisms at miRNA target sites have been identified and they can represent an interesting source for the identification of functional genetic variability.

Low red/far-red ratios delay spike and stem growth in wheat

The responses to low red light/far-red light (R/FR) ratios simulating dense stands were evaluated in wheat (Triticum aestivum L) cultivars released at different times in the 20th century and consequently resulting from an increasingly prolonged breeding and selection history. While tillering responses to the R/FR ratio were unaffected by the cultivars, low R/FR ratios reduced grain yield per plant (primarily grain number and secondarily grain weight per plant) particularly in modern cultivars. Low R/FR ratios delayed spike growth and development, reduced the expression of spike marker genes, accelerated the development of florets already initiated, and reduced the number of fertile florets at anthesis. It is noteworthy that low R/FR ratios did not promote stem or leaf sheath growth and therefore the observed reduction of yield cannot be accounted for as a consequence of divergence of resources towards increased plant stature. It is proposed that the regulation of yield components by the R/FR ratio could help plants to adjust to the limited availability of resources under crop conditions.

Multilevel regulation and signalling processes associated with adaptation to terminal drought in wild emmer wheat

Low water availability is the major environmental factor limiting crop productivity. Transcriptome analysis was used to study terminal drought response in wild emmer wheat, Triticum dicoccoides, genotypes contrasting in their productivity and yield stability under drought stress. A total of 5,892 differentially regulated transcripts were identified between drought and well-watered control and/or between drought resistant (R) and drought susceptible (S) genotypes. Functional enrichment analyses revealed that multilevel regulatory and signalling processes were significantly enriched among the drought-induced transcripts, in particular in the R genotype. Therefore, further analyses were focused on selected 221 uniquely expressed or highly abundant transcripts in the R genotype, as potential candidates for drought resistance genes. Annotation of the 221 genes revealed that 26% of them are involved in multilevel regulation, including: transcriptional regulation, RNA binding, kinase activity and calcium and abscisic acid signalling implicated in stomatal closure. Differential expression patterns were also identified in genes known to be involved in drought adaptation pathways, such as: cell wall adjustment, cuticular wax deposition, lignification, osmoregulation, redox homeostasis, dehydration protection and drought-induced senescence. These results demonstrate the potential of wild emmer wheat as a source for candidate genes for improving drought resistance.

wDBTF: an integrated database resource for studying wheat transcription factor families

BACKGROUND

Transcription factors (TFs) regulate gene expression by interacting with promoters of their target genes and are classified into families based on their DNA-binding domains. Genes coding for TFs have been identified in the sequences of model plant genomes. The rice (Oryza sativa spp. japonica) genome contains 2,384 TF gene models, which represent the mRNA transcript of a locus, classed into 63 families.

RESULTS

We have created an extensive list of wheat (Triticum aestivum L) TF sequences based on sequence homology with rice TFs identified and classified in the Database of Rice Transcription Factors (DRTF). We have identified 7,112 wheat sequences (contigs and singletons) from a dataset of 1,033,960 expressed sequence tag and mRNA (ET) sequences available. This number is about three times the number of TFs in rice so proportionally is very similar if allowance is made for the hexaploidy of wheat. Of these sequences 3,820 encode gene products with a DNA-binding domain and thus were confirmed as potential regulators. These 3,820 sequences were classified into 40 families and 84 subfamilies and some members defined orphan families. The results were compiled in the Database of Wheat Transcription Factor (wDBTF), an inventory available on the web http://wwwappli.nantes.inra.fr:8180/wDBFT/. For each accession, a link to its library source and its Affymetrix identification number is provided. The positions of Pfam (protein family database) motifs were given when known.

CONCLUSIONS

wDBTF collates 3,820 wheat TF sequences validated by the presence of a DNA-binding domain out of 7,112 potential TF sequences identified from publicly available gene expression data. We also incorporated in silico expression data on these TFs into the database. Thus this database provides a major resource for systematic studies of TF families and their expression in wheat as illustrated here in a study of DOF family members expressed during seed development.