Three genomes differentially contribute to the biosynthesis of benzoxazinones in hexaploid wheat

The transcription factor TaMYB31 regulates the benzoxazinoid biosynthetic pathway in wheat

Benzoxazinoids are specialized metabolites that are highly abundant in staple crops, such as maize and wheat. Although their biosynthesis has been studied for several decades, the regulatory mechanisms of the benzoxazinoid pathway remain unknown. Here, we report that the wheat transcription factor MYB31 functions as a regulator of benzoxazinoid biosynthesis genes. A transcriptomic analysis of tetraploid wheat (Triticum turgidum) tissue revealed the up-regulation of two TtMYB31 homoeologous genes upon aphid and caterpillar feeding. TaMYB31 gene silencing in the hexaploid wheat Triticum aestivum significantly reduced benzoxazinoid metabolite levels and led to susceptibility to herbivores. Thus, aphid progeny production, caterpillar body weight gain, and spider mite oviposition significantly increased in TaMYB31-silenced plants. A comprehensive transcriptomic analysis of hexaploid wheat revealed that the TaMYB31 gene is co-expressed with the target benzoxazinoid-encoded Bx genes under several biotic and environmental conditions. Therefore, we analyzed the effect of abiotic stresses on benzoxazinoid levels and discovered a strong accumulation of these compounds in the leaves. The results of a dual fluorescence assay indicated that TaMYB31 binds to the Bx1 and Bx4 gene promoters, thereby activating the transcription of genes involved in the benzoxazinoid pathway. Our finding is the first report of the transcriptional regulation mechanism of the benzoxazinoid pathway in wheat.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Genetic Resources of Cereal Crops for Aphid Resistance

The genetic resources of cereal crops in terms of resistance to aphids are reviewed. Phytosanitary destabilization led to a significant increase in the harmfulness of this group of insects. The breeding of resistant plant genotypes is a radical, the cheapest, and environmentally safe way of pest control. The genetic homogeneity of crops hastens the adaptive microevolution of harmful organisms. Both major and minor aphid resistance genes of cereal plants interact with insects differentially. Therefore, rational breeding envisages the expansion of the genetic diversity of cultivated varieties. The possibilities of replenishing the stock of effective resistance genes by studying the collection of cultivated cereals, introgression, and creating mutant forms are considered. The interaction of insects with plants is subject to the gene-for-gene relationship. Plant resistance genes are characterized by close linkage and multiple allelism. The realizing plant genotype depends on the phytophage biotype. Information about the mechanisms of constitutional and induced plant resistance is discussed. Resistance genes differ in terms of stability of expression. The duration of the period when varieties remain resistant is not related either to its phenotypic manifestation or to the number of resistance genes. One explanation for the phenomenon of durable resistance is the association of the virulence mutation with pest viability.

Polyploidization affects the allelic variation of jasmonate-regulated protein Ta-JA1 belonging to the monocot chimeric jacalin (MCJ) family in wild emmer wheat

The jasmonate-regulated protein Ta-JA1 belongs to the monocot chimeric jacalin (MCJ) family and plays a vital role in stress resistance in wheat. However, the impact of wheat polyploidization on Ta-JA1 remains unclear. In this study, 149 members of the MCJ family were identified among members of Triticeae using a genome-wide approach. The genes were resolved into three clades; MCJ genes in each clade were derived from different donor genes during evolution. Segmental duplication may have been the primary driver, compared with tandem duplication, of expansion in the MCJ family of wheat. Gene loss and acquisition occurred during tetraploidization, and the core expansion of the family occurred after tetraploidization. Sequencing data for 2104 accessions of T. aestivum and 99 accessions of T. dicoccoides showed that Ta-JA1-2A and Ta-JA1 were highly conserved in common wheat, and four alleles (TdJA1-Ax2, TdJA1-Ay2, TdJA1-Ax3, and TdJA1-Ay3) were detected in T. dicoccoides. Using gene-specific markers, one AsJA1-B allele was detected in 11 Ae. speltoides accessions and one TuJA1-Ax1 allele was detected in 70 T. urartu accessions. Six alleles were detected on chromosome 2A: TdJA1-Ax1 (13 accessions), TdJA1-Ay1 (57 accessions), TdJA1-Ax2 (23 accessions), TdJA1-Ay2 (42 accessions), TdJA1-Ax3 (29 accessions), and TdJA1-Ay3 (251 accessions). Only one allele (TdJA1-B) on chromosome 2B was detected in 415 T. dicoccoides accessions. A geographical distribution analysis revealed that Israel hosted higher allelic variation than other regions. Quantitative reverse transcription PCR analysis indicated that divergence in expression has occurred among Ta-JA1 alleles and, notably, TdJA1-Ax1 and TdJA1-Ay1 showed significantly higher expression levels than the other four allelic types in T. dicoccoides. The present results contribute to an improved understanding of the effects of polyploidization on the MCJ gene family and the functions of Ta-JA1, and may be useful to enrich common wheat germplasm resources.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds—namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

Physio-molecular traits of contrasting bread wheat genotypes associated with 15N influx exhibiting homeolog expression bias in nitrate transporter genes under different external nitrate concentrations

The high affinity nitrate transport system is a potential target for improving nitrogen use efficiency of bread wheat growing either under optimal or limiting nitrate concentration. Nitrate uptake is one of the most important traits to take into account to improve nitrogen use efficiency in wheat (Triticum aestivum L.). In this study, we aimed to gain an insight into the regulation of NO₃^- -uptake and translocation systems in two contrasting wheat genotypes [K9107(K9) vs. Choti Lerma (CL)]. Different conditions, such as NO₃^--uptake rates, soil-types, N-free solid external media, and external NO₃^- levels at the seedling stage, were considered. We also studied the contribution of homeolog expression of five genes encoding two nitrate transporters in the root tissue, along with their overall transcript expression levels relative to specific external nitrate availability. We observed that K9107 had a higher ¹⁵N influx than Choti Lerma under both limiting as well as optimum external N conditions in vermiculite-perlite (i.e., N-free solid) medium, with the improved translocation efficiency in Choti Lerma. However, in different soil types, different levels of ¹⁵N-enrichment in both the genotypes were found. Our results also demonstrated that the partitioning of dry matter in root and shoot was different under these growing conditions. Moreover, K9107 showed significantly higher relative expression of TaNRT2.1 at the lowest and TaNPF6.1 and TaNPF6.2 at the highest external nitrate concentrations. We also observed genotype-specific and nitrate starvation-dependent homeolog expression bias in all five nitrate transporter genes. Our data suggest that K9107 had a higher NO₃^- influx capacity, involving different nitrate transporters, than Choti Lerma at the seedling stage.

The wheat dioxygenase BX6 is involved in the formation of benzoxazinoids in planta and contributes to plant defense against insect herbivores

Benzoxazinoids are plant specialized metabolites with defense properties, highly abundant in wheat (Triticum), one of the world's most important crops. The goal of our study was to characterize dioxygenase BX6 genes in tetraploid and hexaploid wheat genotypes and to elucidate their effects on defense against herbivores. Phylogenetic analysis revealed four BX6 genes in the hexaploid wheat T. aestivum, but only one ortholog was found in the tetraploid (T. turgidum) wild emmer wheat and the cultivated durum wheat. Transcriptome sequencing of durum wheat plants, damaged by either aphids or caterpillars, revealed that several BX genes, including TtBX6, were upregulated upon caterpillar feeding, relative to the undamaged control plants. A virus-induced gene silencing approach was used to reduce the expression of BX6 in T. aestivum plants, which exhibited both reduced transcript levels and reduced accumulation of different benzoxazinoids. To elucidate the effect of BX6 on plant defense, bioassays with different herbivores feeding on BX6-silenced leaves were conducted. The results showed that plants with silenced BX6 were more susceptible to aphids and the two-spotted spider mite than the control. Overall, our study indicates that wheat BX6 is involved in benzoxazinoid formation in planta and contributes to plant resistance against insect herbivores.

Engineering of benzoxazinoid biosynthesis in Arabidopsis thaliana: Metabolic and physiological challenges

Plant specialised metabolites constitute a layer of chemical defence. Classes of the defence compounds are often restricted to a certain taxon of plants, e.g. benzoxazinoids (BX) are characteristically detected in grasses. BXs confer wide-range defence by controlling herbivores and microbial pathogens and are allelopathic compounds. In the crops maize, wheat and rye high concentrations of BXs are synthesised at an early developmental stage. By transfer of six Bx-genes (Bx1 to Bx5 and Bx8) it was possible to establish the biosynthesis of 2,4-dihydroxy-1,4-benzoxazin-3-one glucoside (GDIBOA) in a concentration of up to 143 nmol/g dry weight in Arabidopsis thaliana. Our results indicate that inefficient channeling of substrates along the pathway and metabolisation of intermediates in host plants might be a general drawback for transgenic establishment of specialised metabolite biosynthesis pathways. As a consequence, BX levels required for defence are not obtained in Arabidopsis. We could show that indolin-2-one (ION), the first specific intermediate, is phytotoxic and is metabolised by hydroxylation and glycosylation by a wide spectrum of plants. In Arabidopsis, metabolic stress due to the enrichment of ION leads to elevated levels of salicylic acid (SA) and in addition to its intrinsic phytotoxicity, ION affects plant morphology indirectly via SA. We could show that Bx3 has a crucial role in the evolution of the pathway, first based on its impact on flux into the pathway and, second by C3-hydroxylation of the phytotoxic ION. Thereby BX3 interferes with a supposedly generic detoxification system towards the non-specific intermediate.

Activation of Cryptic Secondary Metabolite Biosynthesis in Bamboo Suspension Cells by a Histone Deacetylase Inhibitor

Plants have evolved a diverse array of secondary metabolite biosynthetic pathways. Undifferentiated plant cells, however, tend to biosynthesize secondary metabolites to a lesser extent and sometimes not at all. This phenomenon in cultured cells is associated with the transcriptional suppression of biosynthetic genes due to epigenetic alterations, such as low histone acetylation levels and/or high DNA methylation levels. Here, using cultured cells of bamboo (Bambusa multiplex; Bm) as a model system, we investigated the effect of histone deacetylase (HDAC) inhibitors on the activation of cryptic secondary metabolite biosynthesis. The Bm suspension cells cultured in the presence of an HDAC inhibitor, suberoyl bis-hydroxamic acid (SBHA), exhibited strong biosynthesis of some compounds that are inherently present at very low levels in Bm cells. Two major compounds induced by SBHA were isolated and were identified as 3-O-p-coumaroylquinic acid (1) and 3-O-feruloylquinic acid (2). Their productivities depended on the type of basal culture medium, initial cell density, and culture period, as well as the SBHA concentration. The biosynthesis of these two compounds was also induced by another HDAC inhibitor, trichostatin A. These results demonstrate the usefulness of HDAC inhibitors to activate cryptic secondary metabolite biosynthesis in cultured plant cells.

The Effectiveness of Physical and Chemical Defense Responses of Wild Emmer Wheat Against Aphids Depends on Leaf Position and Genotype

The bird cherry-oat aphid (Rhopalosiphum padi) is one of the most destructive insect pests in wheat production. To reduce aphid damage, wheat plants have evolved various chemical and physical defense mechanisms. Although these mechanisms have been frequently reported, much less is known about their effectiveness. The tetraploid wild emmer wheat (WEW; Triticum turgidum ssp. dicoccoides), one of the progenitors of domesticated wheat, possesses untapped resources from its numerous desirable traits, including insect resistance. The goal of this research was to determine the effectiveness of trichomes (physical defense) and benzoxazinoids (BXDs; chemical defense) in aphid resistance by exploiting the natural diversity of WEW. We integrated a large dataset composed of trichome density and BXD abundance across wheat genotypes, different leaf positions, conditions (constitutive and aphid-induced), and tissues (whole leaf and phloem sap). First, we evaluated aphid reproduction on 203 wheat accessions and found large variation in this trait. Then, we chose eight WEW genotypes and one domesticated durum wheat cultivar for detailed quantification of the defense mechanisms across three leaves. We discovered that these defense mechanisms are influenced by both leaf position and genotype, where aphid reproduction was the highest on leaf-1 (the oldest), and trichome density was the lowest. We compared the changes in trichome density and BXD levels upon aphid infestation and found only minor changes relative to untreated plants. This suggests that the defense mechanisms in the whole leaf are primarily anticipatory and unlikely to contribute to aphid-induced defense. Next, we quantified BXD levels in the phloem sap and detected a significant induction of two compounds upon aphid infestation. Moreover, evaluating aphid feeding patterns showed that aphids prefer to feed on the oldest leaf. These findings revealed the dynamic response at the whole leaf and phloem levels that altered aphid feeding and reproduction. Overall, they suggested that trichomes and the BXD 2,4-dihydroxy-7- methoxy-1,4-benzoxazin-3-one (DIMBOA) levels are the main factors determining aphid resistance, while trichomes are more effective than BXDs. Accessions from the WEW germplasm, rich with trichomes and BXDs, can be used as new genetic sources to improve the resistance of elite wheat cultivars.

Increased benzoxazinoid (Bx) levels in wheat seedlings via jasmonic acid treatment and etiolation and their effects on Bx genes including Bx6

Wheat accumulates benzoxazinoid (Bx) as a defensive compound. While Bx occurs at high concentrations, particularly in the early growth stages, its mechanism of regulation remains unclear. In the present study, we first examined the effects of several plant hormones on Bx concentrations in wheat seedlings. Among the compounds tested, jasmonate (JA) elevated the concentrations of DIMBOA-Glc (2-β-D-glucoside of 2,4-dihydroxy-7-methoy-1,4-benzoxazin-3-one), the primary Bx species in intact wheat seedlings, without a significant increase in HDMBOA-Glc (4-O-methyl-DIMBOA-Glc), which is known to be upregulated by stresses. In addition, growing the plants in the dark increased DIMBOA-Glc levels. Quantification of the Bx-biosynthetic genes showed that TaBx8 (UDP-Glc:Bx glucosyltrasferase) was influenced by neither JA nor etiolation, indicating that TaBx8 is under the regulation mechanism distinct from the mechanisms influencing the others. In addition, none of the other gene expression patterns exhibited considerable correlation with DIMBOA-Glc accumulation. Since there was no correlation between transcript levels of the genes involved in Bx biosynthesis and Bx accumulation, other factors may control the levels of Bx in wheat. In the course of gene analyses, we isolated TaBx6, one of the last two genes that had not been identified in wheat in the DIMBOA-Glc biosynthetic pathway. All the four TaBx6 genes cloned in the present study were expressed in Escherichia coli and characterized their activity.

The Roots of Rye (Secale cereale L.) Are Capable of Synthesizing Benzoxazinoids

According to current opinion, the first step of benzoxazinoids (BXs) synthesis, that is, the conversion of indole-3-glycerol phosphate to indole, occurs exclusively in the photosynthesising parts of plants. However, the results of our previous work and some other studies suggest that this process may also occur in the roots. In this study, we provide evidence that the first step of BXs synthesis does indeed occur in the roots of rye seedlings. We detected ScBx1 transcripts, BX1 enzyme, and six BXs (2-hydroxy-1,4-benzoxazin-3-one, 2,4-dihydroxy-1,4-benzoxazin-3-one, (2R)-2-O-β-d-glucopyranosyl-4-hydroxy-(2H)-1,4-benzoxazin-3(4H)-one glucoside, 2,4-dihydroxy- 7-methoxy-1,4-benzoxazin-3-one, 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one glucoside, and 6-methoxy-2-benzoxazolinone) in the roots developed from seeds deprived of the coleoptile at 2 days after sowing (i.e., roots without contact with aerial parts). In roots regenerated in vitro, both ScBx1 transcripts and BX1 enzyme were detected at a low but still measurable levels. Thus, BXs are able to be synthesised in both the roots and above-ground parts of rye plants.

Analysis of the chromosomal clustering of Fusarium-responsive wheat genes uncovers new players in the defence against head blight disease

There is increasing evidence that some functionally related, co-expressed genes cluster within eukaryotic genomes. We present a novel pipeline that delineates such eukaryotic gene clusters. Using this tool for bread wheat, we uncovered 44 clusters of genes that are responsive to the fungal pathogen Fusarium graminearum. As expected, these Fusarium-responsive gene clusters (FRGCs) included metabolic gene clusters, many of which are associated with disease resistance, but hitherto not described for wheat. However, the majority of the FRGCs are non-metabolic, many of which contain clusters of paralogues, including those implicated in plant disease responses, such as glutathione transferases, MAP kinases, and germin-like proteins. 20 of the FRGCs encode nonhomologous, non-metabolic genes (including defence-related genes). One of these clusters includes the characterised Fusarium resistance orphan gene, TaFROG. Eight of the FRGCs map within 6 FHB resistance loci. One small QTL on chromosome 7D (4.7 Mb) encodes eight Fusarium-responsive genes, five of which are within a FRGC. This study provides a new tool to identify genomic regions enriched in genes responsive to specific traits of interest and applied herein it highlighted gene families, genetic loci and biological pathways of importance in the response of wheat to disease.

Transcriptome Analysis Reveals a Potential Role of Benzoxazinoid in Regulating Stem Elongation in the Wheat Mutant qd

The stems of cereal crops provide both mechanical support for lodging resistance and a nutrient supply for reproductive organs. Elongation, which is considered a critical phase for yield determination in winter wheat (Triticum aestivum L.), begins from the first node detectable to anthesis. Previously, we characterized a heavy ion beam triggered wheat mutant qd, which exhibited an altered stem elongation pattern without affecting mature plant height. In this study, we further analyzed mutant stem developmental characteristics by using transcriptome data. More than 40.87 Mb of clean reads including at least 36.61 Mb of unique mapped reads were obtained for each biological sample in this project. We utilized our transcriptome data to identify 124,971 genes. Among these genes, 4,340 differentially expressed genes (DEG) were identified between the qd and wild-type (WT) plants. Compared to their WT counterparts, qd plants expressed 2,462 DEGs with downregulated expression levels and 1878 DEGs with upregulated expression levels. Using DEXSeq, we identified 2,391 counting bins corresponding to 1,148 genes, and 289 of them were also found in the DEG analysis, demonstrating differences between qd and WT. The 5,199 differentially expressed genes between qd and WT were employed for GO and KEGG analyses. Biological processes, including protein-DNA complex subunit organization, protein-DNA complex assembly, nucleosome organization, nucleosome assembly, and chromatin assembly, were significantly enriched by GO analysis. However, only benzoxazinoid biosynthesis pathway-associated genes were enriched by KEGG analysis. Genes encoding the benzoxazinoid biosynthesis enzymes Bx1, Bx3, Bx4, Bx5, and Bx8_9 were confirmed to be differentially expressed between qd and WT. Our results suggest that benzoxazinoids could play critical roles in regulating the stem elongation phenotype of qd.

Cytochrome P450-mediated herbicide metabolism in plants: current understanding and prospects

Cytochrome P450s (P450s) have been at the center of herbicide metabolism research as a result of their ability to endow selectivity in crops and resistance in weeds. In the last 20 years, ≈30 P450s from diverse plant species have been revealed to possess herbicide-metabolizing function, some of which were demonstrated to play a key role in plant herbicide sensitivity. Recent research even demonstrated that some P450s from crops and weeds metabolize numerous herbicides from various chemical backbones, which highlights the importance of P450s in the current agricultural systems. However, due to the enormous number of plant P450s and the complexity of their function, expression and regulation, it remains a challenge to fully explore the potential of P450-mediated herbicide metabolism in crop improvement and herbicide resistance mitigation. Differences in the substrate specificity of each herbicide-metabolizing P450 are now evident. Comparisons of the substrate specificity and protein structures of P450s will be beneficial for the discovery of selective herbicides and may lead to the development of crops with higher herbicide tolerance by transgenics or genome-editing technologies. Furthermore, the knowledge will help design sound management strategies for weed resistance including the prediction of cross-resistance patterns. Overcoming the ambiguity of P450 function in plant xenobiotic pathways will unlock the full potential of this enzyme family in advancing global agriculture and food security. © 2020 Society of Chemical Industry.

Benzoxazinoids Biosynthesis in Rye (Secale cereale L.) Is Affected by Low Temperature

Benzoxazinoids (BXs) are specialized metabolites with protective properties that are synthesized predominantly by Poaceae species, including rye (Secale cereale). Among factors known to influence BXs production, prolonged low temperature has not been studied previously. In this study, the influence of cultivation at 4 °C, which is essential for vernalization, on the concentration of BXs (HBOA, DIBOA, GDIBOA, DIMBOA, GDIMBOA, and MBOA) and the expression level of genes involved in the BX biosynthesis pathway (ScBx1–ScBx5 and ScIgl) in three rye inbred lines was investigated. After cultivation for seven weeks at 4 °C, the expression level of all analyzed genes and BX concentrations had decreased compared with those at the initiation of treatment (21 days after germination) in control and cold-treated plants. At this time point, the decrease in BX concentrations and gene expression was lower in cold-treated plants than in untreated plants. In contrast, at 77 days after germination, the gene expression levels and BX concentrations in untreated plants had generally increased. Investigation of the vernalization impact on rye BXs accumulation, as well as on Bx gene expression, may aid with determination of the most suitable winter lines and cultivars of rye for cultivation and breeding purposes.

Transcriptome Dynamics of the Inflorescence in Reciprocally Formed Allopolyploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary mechanism and is prevalent among land plants. Most polyploid species examined have multiple origins, which provide genetic diversity and may enhance the success of polyploids. In some polyploids, recurrent origins can result from reciprocal crosses between the same diploid progenitors. Although great progress has been made in understanding the genetic consequences of polyploidy, the genetic implications of reciprocal polyploidization remain poorly understood, especially in natural polyploids. Tragopogon (Asteraceae) has become an evolutionary model system for studies of recent and recurrent polyploidy. Allotetraploid T. miscellus has formed reciprocally in nature with resultant distinctive floral and inflorescence morphologies (i.e., short- vs. long-liguled forms). In this study, we performed comparative inflorescence transcriptome analyses of reciprocally formed T. miscellus and its diploid parents, T. dubius and T. pratensis. In both forms of T. miscellus, homeolog expression of ∼70% of the loci showed vertical transmission of the parental expression patterns (i.e., parental legacy), and ∼20% of the loci showed biased homeolog expression, which was unbalanced toward T. pratensis. However, 17.9% of orthologous pairs showed different homeolog expression patterns between the two forms of T. miscellus. No clear effect of cytonuclear interaction on biased expression of the maternal homeolog was found. In terms of the total expression level of the homeologs studied, 22.6% and 16.2% of the loci displayed non-additive expression in short- and long-liguled T. miscellus, respectively. Unbalanced expression level dominance toward T. pratensis was observed in both forms of T. miscellus. Significantly, genes annotated as being involved in pectin catabolic processes were highly expressed in long-liguled T. miscellus relative to the short-liguled form, and the majority of these differentially expressed genes were transgressively down-regulated in short-liguled T. miscellus. Given the known role of these genes in cell expansion, they may play a role in the differing floral and inflorescence morphologies of the two forms. In summary, the overall inflorescence transcriptome profiles are highly similar between reciprocal origins of T. miscellus. However, the dynamic homeolog-specific expression and non-additive expression patterns observed in T. miscellus emphasize the importance of reciprocal origins in promoting the genetic diversity of polyploids.

Comparative functional genomics analysis of cytochrome P450 gene superfamily in wheat and maize

BACKGROUND

The cytochrome P450s (CYP450s) as the largest enzyme family of plant metabolism participate in various physiological processes, whereas no study has demonstrated interest in comprehensive comparison of the genes in wheat and maize. Genome-wide survey, characterization and comparison of wheat and maize CYP450 gene superfamily are useful for genetic manipulation of the Gramineae crops.

RESULTS

In total, 1285 and 263 full-length CYP450s were identified in wheat and maize, respectively. According to standard nomenclature, wheat CYP450s (TaCYP450s) were categorized into 45 families, while maize CYP450s (ZmCYP450s) into 43 families. A comprehensive analysis of wheat and maize CYP450s, involved in functional domains, conserved motifs, phylogeny, gene structures, chromosome locations and duplicated events was performed. The result showed that each family/subfamily in both species exhibited characteristic features, suggesting their phylogenetic relationship and the potential divergence in their functions. Functional divergence analysis at the amino acid level of representative clans CYP51, CYP74 and CYP97 in wheat, maize and rice identified some critical amino acid sites that are responsible for functional divergence of a gene family. Expression profiles of Ta-, ZmCYP450s were investigated using RNA-seq data, which contribute to infer the potential functions of the genes during development and stress responses. We found in both species CYP450s had preferential expression in specific tissues, and many tissue-specific genes were identified. Under water-deficit condition, 82 and 39 significantly differentially expressed CYP450s were respectively detected in wheat and maize. These genes may have some roles in protecting plants against drought damage. Thereinto, fourteen CYP450s were selected to validate their expression level through qRT-PCR. To further elucidating molecular mechanisms of CYP450 action, gene co-expression network was constructed. In total, 477 TaCYP450s were distributed in 22 co-expression modules, and some co-expressed genes that likely take part in the same biochemical pathway were identified. For instance, the expression of TaCYP74A98_4D was highly correlated with TaLOX9, TaLOX36, TaLOX39, TaLOX44 and TaOPR8, and all of them may be involved in jasmonate (JA) biosynthesis. TaCYP73A201_3A showed coexpression with TaPAL1.25, TaCCoAOMT1.2, TaCOMT.1, TaCCR1.6 and TaLAC5, which probably act in the wheat stem and/or root lignin synthesis pathway.

CONCLUSION

Our study first established systematic information about evolutionary relationship, expression pattern and function characterization of CYP450s in wheat and maize.

TaCYP81D5, one member in a wheat cytochrome P450 gene cluster, confers salinity tolerance via reactive oxygen species scavenging

As one of the largest gene families in plants, the cytochrome P450 monooxygenase genes (CYPs) are involved in diverse biological processes including biotic and abiotic stress response. Moreover, P450 genes are prone to expanding due to gene tandem duplication during evolution, resulting in generations of novel alleles with the neo-function or enhanced function. Here, the bread wheat (Triticum aestivum) gene TaCYP81D5 was found to lie within a cluster of five tandemly arranged CYP81D genes, although only a single such gene (BdCYP81D1) was present in the equivalent genomic region of the wheat relative Brachypodium distachyon. The imposition of salinity stress could up-regulate TaCYP81D5, but the effect was abolished in plants treated with an inhibitor of reactive oxygen species synthesis. In SR3, a wheat cultivar with an elevated ROS content, the higher expression and the rapider response to salinity of TaCYP81D5 were related to the chromatin modification. Constitutively expressing TaCYP81D5 enhanced the salinity tolerance both at seedling and reproductive stages of wheat via accelerating ROS scavenging. Moreover, an important component of ROS signal transduction, Zat12, was proven crucial in this process. Though knockout of solely TaCYP81D5 showed no effect on salinity tolerance, knockdown of BdCYP81D1 or all TaCYP81D members in the cluster caused the sensitivity to salt stress. Our results provide the direct evidence that TaCYP81D5 confers salinity tolerance in bread wheat and this gene is prospective for crop improvement.

Comparative transcriptomic and metabolic analysis of wild and domesticated wheat genotypes reveals differences in chemical and physical defense responses against aphids

BACKGROUND

Young wheat plants are continuously exposed to herbivorous insect attack. To reduce insect damage and maintain their growth, plants evolved different defense mechanisms, including the biosynthesis of deterrent compounds named benzoxazinoids, and/or trichome formation that provides physical barriers. It is unclear whether both of these mechanisms are equally critical in providing an efficient defense for wheat seedlings against aphids-an economically costly pest in cereal production.

RESULTS

In this study, we compared the transcriptome, metabolome, benzoxazinoids, and trichome density of three selected wheat genotypes, with a focus on differences related to defense mechanisms. We chose diverse wheat genotypes: two tetraploid wheat genotypes, domesticated durum 'Svevo' and wild emmer 'Zavitan,' and one hexaploid bread wheat, 'Chinese Spring.' The full transcriptomic analysis revealed a major difference between the three genotypes, while the clustering of significantly different genes suggested a higher similarity between the two domesticated wheats than between either and the wild wheat. A pathway enrichment analysis indicated that the genes associated with primary metabolism, as well as the pathways associated with defense such as phytohormones and specialized metabolites, were different between the three genotypes. Measurement of benzoxazinoid levels at the three time points (11, 15, and 18 days after germination) revealed high levels in the two domesticated genotypes, while in wild emmer wheat, they were below detection level. In contrast to the benzoxazinoid levels, the trichome density was dramatically higher in the wild emmer than in the domesticated wheat. Lastly, we tested the bird cherry-oat aphid's (Rhopalosiphum padi) performance and found that Chinese Spring is more resistant than the tetraploid genotypes.

CONCLUSIONS

Our results show that benzoxazinoids play a more significant defensive role than trichomes. Differences between the abundance of defense mechanisms in the wild and domesticated plants were observed in which wild emmer possesses high physical defenses while the domesticated wheat genotypes have high chemical defenses. These findings provide new insights into the defense adaptations of wheat plants against aphids.

Molecular Characterization of Auxin Efflux Carrier- ABCB1 in hexaploid wheat

Auxin is an important phytohormone that regulates response, differentiation, and development of plant cell, tissue, and organs. Along with its local production, long-distance transport coordinated by the efflux/influx membrane transporters is instrumental in plant development and architecture. In the present study, we cloned and characterized a wheat (Triticum aestivum) auxin efflux carrier ABCB1. The TaABCB1 was physically localized to the proximal 15% of the short arm of wheat homoeologous group 7 chromosomes. Size of the Chinese spring (CS) homoeologs genomic copies ranged from 5.3–6.2 kb with the 7A copy being the largest due to novel insertions in its third intron. The three homoeologous copies share 95–97% sequence similarity at the nucleotide, 98–99% amino acid, and overall Q-score of 0.98 at 3-D structure level. Though detected in all analyzed tissues, TaABCB1 predominantly expressed in the meristematic tissues likely due to the presence of meristem-specific activation regulatory element identified in the promoter region. RNAi plants of TaABCB1 gene resulted in reduced plant height and increased seed width. Promoter analysis revealed several responsive elements detected in the promoter region including that for different hormones as auxin, gibberellic acid, jasmonic acid and abscisic acid, light, and circadian regulated elements.

Deletion of the benzoxazinoid detoxification gene NAT1 in Fusarium graminearum reduces deoxynivalenol in spring wheat

Benzoxazinoid (Bx) metabolites produced by wheat and other members of the Poaceae have activity against Fusarium sp. that cause cereal diseases including Fusarium head blight (FHB) on wheat and barley. Certain Bx metabolites can be detoxified by Fusarium sp. with the arylamine N-acetyltransferase NAT1. Investigation of this pathway may reveal strategies for increasing FHB resistance, such as selection for higher levels of Bx metabolites within existing germplasm and/or engineering fungal susceptibility via host induced silencing of NAT1. We assessed the reactions of fifteen wheat cultivars or breeding lines adapted to the Northwestern United States to infection with F. graminearum Δnat1 mutants that should be sensitive to Bx metabolites. Significant differences were noted in disease severity and deoxynivalenol (DON) among the cultivars 21 d after inoculation with either mutant or wildtype (PH1) strains. Mutant vs. wildtype strains did not result in significant variation for infection severity (as measured by % infected florets), but inoculation with Δnat1 mutants vs. wildtype resulted in significantly lower DON concentrations in mature kernels (p < 0.0001). Of the cultivars tested, HRS3419 was the most resistant cultivar to PH1 (severity = 62%, DON = 45 ppm) and Δnat1 mutants (severity = 61%, DON = 30 ppm). The cultivar most susceptible to infection was Kelse with PH1 (severity = 100%, DON = 292 ppm) and Δnat1 mutants (severity = 100%, DON = 158 ppm). We hypothesized that sub-lethal Bx metabolite levels may suppress DON production in F. graminearum Δnat1 mutants. In vitro assays of Bx metabolites BOA, MBOA, and DIMBOA at 30 μM did not affect growth, but did reduce DON production by Δnat1 and PH1. Although the levels of Bx metabolites are likely too low in the wheat cultivars we tested to suppress FHB, higher levels of Bx metabolites may contribute towards reductions in DON and FHB.

Cereal aphids differently affect benzoxazinoid levels in durum wheat

Aphids are major pests in cereal crops that cause direct and indirect damage leading to yield reduction. Despite the fact that wheat provides 20% of the world’s caloric and protein diet, its metabolic responses to aphid attack, in general, and specifically its production of benzoxazinoid defense compounds are poorly understood. The objective of this study was to compare the metabolic diversity of durum wheat seedlings (Triticum turgidum ssp. durum) under attack by three different cereal aphids: i) the English grain aphid (Sitobion avenae Fabricius), ii) the bird cherry-oat aphid (Rhopalosiphum padi L.), and iii) the greenbug aphid (Schizaphis graminum Rondani), which are some of the most destructive aphid species to wheat. Insect progeny bioassays and metabolic analyses using chromatography/Q-Exactive/mass spectrometry non-targeted metabolomics and a targeted benzoxazinoid profile were performed on infested leaves. The insect bioassays revealed that the plants were susceptible to S. graminum, resistant to S. avenae, and mildly resistant to R. padi. The metabolic analyses of benzoxazinoids suggested that the predominant metabolites DIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin- 3-one) and its glycosylated form DIMBOA-glucoside (Glc) were significantly induced upon both S. avenae, and R. padi aphid feeding. However, the levels of the benzoxazinoid metabolite HDMBOA-Glc (2-hydroxy-4,7-dimethoxy-1,4-benzoxazin-3-one glucoside) were enhanced due to the feeding of S. avenae and S. graminum aphids, to which Svevo was the most resistant and the most susceptible, respectively. The results showed a partial correlation between the induction of benzoxazinoids and aphid reproduction. Overall, our observations revealed diverse metabolic responses of wheat seedlings to cereal aphid feeding.

Three genomes differentially contribute to the seedling lateral root number in allohexaploid wheat: evidence from phenotype evolution and gene expression

Common wheat is an allohexaploid (BBAADD) that originated from the hybridization and polyploidization of the diploid Aegilops tauschii (DD) with the allotetraploid Triticum turgidum (BBAA). Phenotypic changes often arise with the formation and evolution of allopolyploid wheat, but little is known about the evolution of root traits in different wheat species with varying ploidy levels. Here, we reported that the lateral root number on the primary root (LRNPR) of synthetic and natural allohexaploid wheats (BBAADD) is significantly higher than that of their allotetraploid (BBAA) and diploid (AA and SS) progenitors, but is much lower than that of their diploid (DD) progenitors. The expression of the wheat gene TaLBD16, an ortholog of the Arabidopsis LATERAL ORGAN BOUNDARIES-DOMAIN16/ASYMMETRIC LEAVES2-LIKE18 (LBD16), which is involved in lateral root development in Arabidopsis, was positively correlated with the LRNPR in diploid and allopolyploid wheats. In natural and synthetic allohexaploid wheats, the transcript of the TaLBD16 from the D genome (TaLBD16-D) was relatively more abundant compared with TaLBD16-A and TaLBD16-B. Consistent with the observed variation in LRNPR, the divergence in the expression of TaLBD16 homoeologous genes occurred before the formation of polyploidy wheat. Collectively, our observations indicate that the D genome played a crucial role in the increased lateral root number of allohexaploid wheats compared with their allotetraploid progenitors, and that TaLBD16-D was one of the key genes involved in the formation of lateral root number during wheat evolution.

Plant Protection by Benzoxazinoids—Recent Insights into Biosynthesis and Function

Benzoxazinoids (BXs) are secondary metabolites present in many Poaceae including the major crops maize, wheat, and rye. In contrast to other potentially toxic secondary metabolites, BXs have not been targets of counter selection during breeding and the effect of BXs on insects, microbes, and neighbouring plants has been recognised. A broad knowledge about the mode of action and metabolisation in target organisms including herbivorous insects, aphids, and plants has been gathered in the last decades. BX biosynthesis has been elucidated on a molecular level in crop cereals. Recent advances, mainly made by investigations in maize, uncovered a significant diversity in the composition of BXs within one species. The pattern can be specific for single plant lines and dynamic changes triggered by biotic and abiotic stresses were observed. Single BXs might be toxic, repelling, attractive, and even growth-promoting for insects, depending on the particular species. BXs delivered into the soil influence plant and microbial communities. Furthermore, BXs can possibly be used as signalling molecules within the plant. In this review we intend to give an overview of the current data on the biosynthesis, structure, and function of BXs, beyond their characterisation as mere phytotoxins.

Transcriptome analysis of Brachypodium during fungal pathogen infection reveals both shared and distinct defense responses with wheat

Fusarium crown rot (FCR) of wheat and barley, predominantly caused by the fungal pathogen Fusarium pseudograminearum, is a disease of economic significance. The quantitative nature of FCR resistance within cultivated wheat germplasm has significantly limited breeding efforts to enhanced FCR resistance in wheat. In this study, we characterized the molecular responses of Brachypodium distachyon (Brachypodium hereafter) to F. pseudograminearum infection using RNA-seq to determine whether Brachypodium can be exploited as a model system towards better understanding of F. pseudograminearum-wheat interaction. The transcriptional response to infection in Brachypodium was strikingly similar to that previously reported in wheat, both in shared expression patterns of wheat homologs of Brachypodium genes and functional overlap revealed through comparative gene ontology analysis in both species. Metabolites produced by various biosynthetic pathways induced in both wheat and Brachypodium were quantified, revealing a high degree of overlap between these two species in metabolic response to infection but also showed Brachypodium does not produce certain defence-related metabolites found in wheat. Functional analyses of candidate genes identified in this study will improve our understanding of resistance mechanisms and may lead to the development of new strategies to protect cereal crops from pathogen infection.

Altered expression of the TaRSL2 gene contributed to variation in root hair length during allopolyploid wheat evolution

Altered expression of the TaRSL2 gene was positively correlated with variation in root hair length during allopolyploid wheat evolution, and overexpression of TaRSL2 in Arabidopsis increases root hair length. Root hairs aid nutrient and water uptake and anchor the plant in the soil. Allopolyploid wheats display significant growth vigor in terms of root hair length compared to their diploid progenitors, but little is known about the molecular basis of variation in root hair length during wheat allopolyploidization. Here, we isolated three orthologs of the Arabidopsis root hair gene ROOT HAIR DEFECTIVE SIX-LIKE 2 (AtRSL2) in allohexaploid wheat, designated TaRSL2-4A, TaRSL2-4B and TaRSL2-4D. The deduced polypeptides of these three TaRSL2 homoeologous genes shared high similarity, and a conserved basic helix-loop-helix (bHLH) domain was present in their C-terminal regions. Notably, the expression of TaRSL2 was positively correlated with root hair length of wheat accessions with different ploidy levels. Moreover, ectopic overexpression of TaRSL2-4D in Arabidopsis could increase root hair length. We found that the transcript levels of TaRSL2 homoeologous genes dynamically changed during allopolyploid wheat evolution, implicating the complexity of the underlying molecular mechanism. Collectively, we propose that altered expression of the TaRSL2 gene contributed to variation in root hair length in allopolyploid wheats.

Biosynthesis and chemical transformation of benzoxazinoids in rye during seed germination and the identification of a rye Bx6-like gene

Benzoxazinoids are secondary metabolites with plant defense properties and possible health-promoting effects in humans. In this study, the transcriptional activity of ScBx genes (ScBx1-ScBx5; ScBx6-like), involved in benzoxazinoid biosynthesis, was analyzed during germination and early seedling development in rye. Our results showed that ScBx genes had highest levels of expression at 24-30 h after germination, followed by a decrease at later stages. For ScBx1-ScBx5 genes expression was higher in shoots compared with root tissues and vice versa for ScBx6-like gene transcripts. Moreover, methylated forms of benzoxazinoids accumulated in roots rather than in shoots during seedling development, in particular reaching high levels of HMBOA-glc in roots. Chemical profiles of benzoxazinoid accumulation in the developing seedling reflected the combined effects of de novo biosynthesis of the compounds as well as the turnover of compounds either pre-stored in the embryo or de novo biosynthesized. Bioinformatic analysis, together with the differential distribution of ScBx6-like transcripts in root and shoot tissues, suggested the presence of a ZmBx6 homolog encoding a 2-oxoglutarate dependent dehydrogenase in rye. The ScBx6-like cDNA was expressed in E. coli for functional characterization in vitro. LC-MS/MS analysis showed that the purified enzyme was responsible for the oxidation of DIBOA-glc into TRIBOA-glc, strongly suggesting the ScBX6-like enzyme in rye to be a functional ortholog of maize ZmBX6.

Evolution of Gene Expression Balance Among Homeologs of Natural Polyploids

Polyploidy is a major evolutionary process in eukaryotes, yet the expression balance of homeologs in natural polyploids is largely unknown. To study this expression balance, the expression patterns of 2180 structurally well-characterized genes of wheat were studied, of which 813 had the expected three copies and 375 had less than three. Copy numbers of the remaining 992 ranged from 4 to 14, including homeologs, orthologs, and paralogs. Of the genes with three structural copies corresponding to homeologs, 55% expressed from all three, 38% from two, and the remaining 7% expressed from only one of the three copies. Homeologs of 76-87% of the genes showed differential expression patterns in different tissues, thus have evolved different gene expression controls, possibly resulting in novel functions. Homeologs of 55% of the genes showed tissue-specific expression, with the largest percentage (14%) in the anthers and the smallest (7%) in the pistils. The highest number (1.72/3) of homeologs/gene expression was in the roots and the lowest (1.03/3) in the anthers. As the expression of homeologs changed with changes in structural copy number, about 30% of the genes showed dosage dependence. Chromosomal location also impacted expression pattern as a significantly higher proportion of genes in the proximal regions showed expression from all three copies compared to that present in the distal regions.

The defence-associated transcriptome of hexaploid wheat displays homoeolog expression and induction bias

Bread wheat (Triticum aestivum L.) is an allopolyploid species containing three ancestral genomes. Therefore, three homoeologous copies exist for the majority of genes in the wheat genome. Whether different homoeologs are differentially expressed (homoeolog expression bias) in response to biotic and abiotic stresses is poorly understood. In this study, we applied a RNA-seq approach to analyse homoeolog-specific global gene expression patterns in wheat during infection by the fungal pathogen Fusarium pseudograminearum, which causes crown rot disease in cereals. To ensure specific detection of homoeologs, we first optimized read alignment methods and validated the results experimentally on genes with known patterns of subgenome-specific expression. Our global analysis identified widespread patterns of differential expression among homoeologs, indicating homoeolog expression bias underpins a large proportion of the wheat transcriptome. In particular, genes differentially expressed in response to Fusarium infection were found to be disproportionately contributed from B and D subgenomes. In addition, we found differences in the degree of responsiveness to pathogen infection among homoeologous genes with B and D homoeologs exhibiting stronger responses to pathogen infection than A genome copies. We call this latter phenomenon as 'homoeolog induction bias'. Understanding how homoeolog expression and induction biases operate may assist the improvement of biotic stress tolerance in wheat and other polyploid crop species.

The Fusarium crown rot pathogen Fusarium pseudograminearum triggers a suite of transcriptional and metabolic changes in bread wheat (Triticum aestivum L.)

Background and Aims

Fusarium crown rot caused by the fungal pathogen Fusarium pseudograminearum is a disease of wheat and barley, bearing significant economic cost. Efforts to develop effective resistance to this disease have been hampered by the quantitative nature of resistance and a lack of understanding of the factors associated with resistance and susceptibility. Here, we aimed to dissect transcriptional responses triggered in wheat by F. pseudograminearum infection.

Methods

We used an RNA-seq approach to analyse host responses during a compatible interaction and identified >2700 wheat genes differentially regulated after inoculation with F. pseudograminearum . The production of a few key metabolites and plant hormones in the host during the interaction was also analysed.

Key Results

Analysis of gene ontology enrichment showed that a disproportionate number of genes involved in primary and secondary metabolism, signalling and transport were differentially expressed in infected seedlings. A number of genes encoding pathogen-responsive uridine-diphosphate glycosyltransferases (UGTs) potentially involved in detoxification of the Fusarium mycotoxin deoxynivalenol (DON) were differentially expressed. Using a F. pseudograminearum DON-non-producing mutant, DON was shown to play an important role in virulence during Fusarium crown rot. An over-representation of genes involved in the phenylalanine, tryptophan and tyrosine biosynthesis pathways was observed. This was confirmed through metabolite analyses that demonstrated tryptamine and serotonin levels are induced after F. pseudograminearum inoculation.

Conclusions

Overall, the observed host response in bread wheat to F. pseudograminearum during early infection exhibited enrichment of processes related to pathogen perception, defence signalling, transport and metabolism and deployment of chemical and enzymatic defences. Additional functional analyses of candidate genes should reveal their roles in disease resistance or susceptibility. Better understanding of host responses contributing to resistance and/or susceptibility will aid the development of future disease improvement strategies against this important plant pathogen.

Identification and VIGS-based characterization of Bx1 ortholog in rye (Secale cereale L.)

The first step of the benzoxazinoid (BX) synthesis pathway is catalyzed by an enzyme with indole-3-glycerol phosphate lyase activity encoded by 3 genes, Bx1, TSA and Igl. A gene highly homologous to maize and wheat Bx1 has been identified in rye. The goal of the study was to analyze the gene and to experimentally verify its role in the rye BX biosynthesis pathway as a rye ortholog of the Bx1 gene. Expression of the gene showed peak values 3 days after imbibition (dai) and at 21 dai it was undetectable. Changes of the BX content in leaves were highly correlated with the expression pattern until 21 dai. In plants older than 21 dai despite the undetectable expression of the analyzed gene there was still low accumulation of BXs. Function of the gene was verified by correlating its native expression and virus-induced silencing with BX accumulation. Barley stripe mosaic virus (BSMV)-based vectors were used to induce transcriptional (TGS) and posttranscriptional (PTGS) silencing of the analyzed gene. Both strategies (PTGS and TGS) significantly reduced the transcript level of the analyzed gene, and this was highly correlated with lowered BX content. Inoculation with virus-based vectors specifically induced expression of the analyzed gene, indicating up-regulation by biotic stressors. This is the first report of using the BSMV-based system for functional analysis of rye gene. The findings prove that the analyzed gene is a rye ortholog of the Bx1 gene. Its expression is developmentally regulated and is strongly induced by biotic stress. Stable accumulation of BXs in plants older than 21 dai associated with undetectable expression of ScBx1 indicates that the function of the ScBx1 in the BX biosynthesis is redundant with another gene. We anticipate that the unknown gene is a putative ortholog of the Igl, which still remains to be identified in rye.

Function and application of a non-ester-hydrolyzing carboxylesterase discovered in tulip

Plants have evolved secondary metabolite biosynthetic pathways of immense rich diversity. The genes encoding enzymes for secondary metabolite biosynthesis have evolved through gene duplication followed by neofunctionalization, thereby generating functional diversity. Emerging evidence demonstrates that some of those enzymes catalyze reactions entirely different from those usually catalyzed by other members of the same family; e.g. transacylation catalyzed by an enzyme similar to a hydrolytic enzyme. Tuliposide-converting enzyme (TCE), which we recently discovered from tulip, catalyzes the conversion of major defensive secondary metabolites, tuliposides, to antimicrobial tulipalins. The TCEs belong to the carboxylesterase family in the α/β-hydrolase fold superfamily, and specifically catalyze intramolecular transesterification, but not hydrolysis. This non-ester-hydrolyzing carboxylesterase is an example of an enzyme showing catalytic properties that are unpredictable from its primary structure. This review describes the biochemical and physiological aspects of tulipalin biogenesis, and the diverse functions of plant carboxylesterases in the α/β-hydrolase fold superfamily.

Distribution of the tryptophan pathway-derived defensive secondary metabolites gramine and benzoxazinones in Poaceae

The Poaceae is a large taxonomic group consisting of approximately 12,000 species and is classified into 12 subfamilies. Gramine and benzoxazinones (Bxs), which are biosynthesized from the tryptophan pathway, are well-known defensive secondary metabolites in the Poaceae. We analyzed the presence or absence of garamine and Bxs in 64 species in the Poaceae by LC-MS/MS. We found that Hordeum brachyantherum and Hakonechloa macra accumulated gramine, but the presence of gramine was limited to small groups of species. We also detected Bxs in four species in the Pooideae and six species in the Panicoideae. In particular, four species in the Paniceae tribe in Panicoideae accumulaed Bxs, indicating that this tribe is a center of the Bx distribution. Bxs were absent in the subfamilies other than Pooideae and Panicoideae. These findings provide an overview of biased distribution of gramine and Bxs in Poaceae species.

The natural plant stress elicitor cis-jasmone causes cultivar-dependent reduction in growth of the stink bug, Euschistus heros and associated changes in flavonoid concentrations in soybean, Glycine max

To test the hypothesis that the plant stress related elicitor cis-jasmone (cJ) provides protection in soybean pods against the seed-sucking stink bug pest, Euschistus heros, the growth of E. heros on cJ-treated pods was investigated using three soybean cultivars differing in insect susceptibility, i.e. BRS 134 (susceptible), IAC 100 (resistant) and Dowling (resistant). E. heros showed reduced weight gain when fed cJ-treated Dowling, whereas no effect on weight gain was observed when fed other treated cultivars. Using analysis of variance, a three factor (cultivar x treatment x time) interaction was observed with concentrations of the flavonoid glycosides daidzin and genistin, and their corresponding aglycones, daidzein and genistein. There were increases in genistein and genistin concentrations in cJ-treated Dowling at 144 and 120 h post treatment, respectively. Higher concentrations of malonyldaidzin and malonylgenistin in Dowling, compared to BRS 134 and IAC 100, were observed independently of time, the highest concentrations being observed in cJ-treated seeds. Levels of glycitin and malonylglycitin were higher in BRS 134 and IAC 100 compared to Dowling. Canonical variate analysis indicated daidzein (in the first two canonical variates) and genistein (in the first only) as important discriminatory variables. These results suggest that cJ treatment leads to an increase in the levels of potentially defensive isoflavonoids in immature soybean seeds, but the negative effect upon E. heros performance is cultivar-dependent.

Distinct expression and function of carotenoid metabolic genes and homoeologs in developing wheat grains

BACKGROUND

β-carotene, the most active provitamin A molecule produced by plants, plays important roles in human nutrition and health. β-carotene does not usually accumulate in the endosperm (i.e. flour) of mature wheat grains, which is a major food source of calories for humans. Therefore, enriching β-carotene accumulation in wheat grain endosperm will enable a sustainable dietary supplementation of provitamin A. Several metabolic genes affecting β-carotene accumulation have already been isolated from wheat, including phytoene synthase 1 (PSY1), lycopene ε-cyclase (LCYe) and carotenoid β-ring hydroxylase1/2 (HYD1/2).

RESULTS

In this work, we cloned and biochemically characterized two carotenoid cleavage dioxygenases (CCDs), CCD1 and CCD4, from wheat. While CCD1 homoeologs cleaved β-apo-8'-carotenal, β-carotene, lutein and zeaxanthin into apocarotenoid products, CCD4 homoeologs were inactive towards these substrates in in vitro assays. When analyzed by real-time qPCR, PSY1, LCYe, HYD1/2 and CCD1/4 homoeologs showed distinct expression patterns in vegetative tissues and sections of developing tetraploid and hexaploid wheat grains, suggesting that carotenoid metabolic genes and homoeologs are differentially regulated at the transcriptional level in wheat.

CONCLUSIONS

The CCD1/4 enzyme activity and the spatial-temporal gene expression data provide critical insights into the specific carotenoid metabolic gene homoeologs that control β-carotene accumulation in wheat grain endosperm, thus establishing the knowledge base for generation of wheat varieties with enhanced β-carotene in the endosperm through breeding and genome editing approaches.

Altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat

Polyploidy is a major driving force in plant evolution and speciation. Phenotypic changes often arise with the formation, natural selection and domestication of polyploid plants. However, little is known about the consequence of hybridization and polyploidization on root hair development. Here, we report that root hair length of synthetic and natural allopolyploid wheats is significantly longer than those of their diploid progenitors, whereas no difference is observed between allohexaploid and allotetraploid wheats. The expression of wheat gene TaRSL4, an orthologue of AtRSL4 controlling the root hair development in Arabidopsis, was positively correlated with the root hair length in diploid and allotetraploid wheats. Moreover, transcript abundance of TaRSL4 homoeologue from A genome (TaRSL4-A) was much higher than those of other genomes in natural allopolyploid wheat. Notably, increased root hair length by overexpression of the TaRSL4-A in wheat led to enhanced shoot fresh biomass under nutrient-poor conditions. Our observations indicate that increased root hair length in allohexaploid wheat originated in the allotetraploid progenitors and altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat.

Genome wide identification of C1-2i zinc finger proteins and their response to abiotic stress in hexaploid wheat

The C1-2i wheat Q-type C2H2 zinc finger protein (ZFP) transcription factor subclass has been reported to play important roles in plant stress responses. This subclass of ZFPs has not been studied in hexaploid wheat (Triticum aestivum) and we aimed to identify all members of this subclass and evaluate their responses to different abiotic stresses causing oxidative stress. Exploiting the recently published wheat draft genome sequence, we identified 53 members (including homoeologs from A, B and D genomes) of the C1-2i wheat Q-type C2H2 ZFPs (TaZFPs) representing 21 genes. Evolution analysis revealed that 9 TaZFPs members are directly inherited from the parents Triticum urartu and Aegilops tauschii, while 15 diverged through neoploidization events. This TaZFP subclass is responsive to the oxidative stress generator H2O2 and to high light, drought stress and flooding. Most TaZFPs are responsive to H2O2 (37/53), high light (44/53), flooding (31/53) or drought (37/53); 32 TaZFPs were up-regulated by at least 3 stresses and 16 were responsive to all stresses tested. A large number of these TaZFPs were physically mapped on different wheat draft genome sequences with known markers useful for QTL mapping. Our results show that the C1-2i subclass of TaZFPs is associated with responses to different abiotic stresses and that most TaZFPs (30/53 or 57 %) are located on group 5 chromosomes known to be involved in environment adaptation. Detailed characterization of these novel wheat TaZFPs and their association to QTL or eQTL may help to design wheat cultivars with improved tolerance to abiotic stress.

Structural characteristics of ScBx genes controlling the biosynthesis of hydroxamic acids in rye (Secale cereale L.)

Benzoxazinoids (BX) are major secondary metabolites of gramineous plants that play an important role in disease resistance and allelopathy. They also have many other unique properties including anti-bacterial and anti-fungal activity, and the ability to reduce alfa–amylase activity. The biosynthesis and modification of BX are controlled by the genes Bx1 ÷ Bx10, GT and glu, and the majority of these Bx genes have been mapped in maize, wheat and rye. However, the genetic basis of BX biosynthesis remains largely uncharacterized apart from some data from maize and wheat. The aim of this study was to isolate, sequence and characterize five genes (ScBx1, ScBx2, ScBx3, ScBx4 and ScBx5) encoding enzymes involved in the synthesis of DIBOA, an important defense compound of rye. Using a modified 3D procedure of BAC library screening, seven BAC clones containing all of the ScBx genes were isolated and sequenced. Bioinformatic analyses of the resulting contigs were used to examine the structure and other features of these genes, including their promoters, introns and 3’UTRs. Comparative analysis showed that the ScBx genes are similar to those of other Poaceae species, especially to the TaBx genes. The polymorphisms present both in the coding sequences and non-coding regions of ScBx in relation to other Bx genes are predicted to have an impact on the expression, structure and properties of the encoded proteins.

Transcript suppression of TaGW2 increased grain width and weight in bread wheat

Bread wheat (Triticum aestivum L.) is a major staple crop in the world. Grain weight is a major factor of grain yield in wheat, and the identification of candidate genes associated with grain weight is very important for high-yield breeding of wheat. TaGW2 is an orthologous gene of rice OsGW2 that negatively regulates the grain width and weight in rice. There are three TaGW2 homoeologs in bread wheat, TaGW2A, TaGW2B, and TaGW2D. In this study, a specific TaGW2-RNA interference (RNAi) cassette was constructed and transformed into a Chinese bread wheat variety 'Shi 4185' with small grain. The transcript levels of TaGW2A, TaGW2B, and TaGW2D were simultaneously downregulated in TaGW2-RNAi transgenic wheat lines. Compared with the controls, TaGW2-underexpressing transgenic lines displayed significantly increases in the grain width and weight, suggesting that TaGW2 negatively regulated the grain width and weight in bread wheat. Further transcript analysis showed that in different bread wheat accessions, the transcript abundance of TaGW2A was negatively associated with the grain width.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Functional genomics of seed dormancy in wheat: advances and prospects

Seed dormancy is a mechanism underlying the inability of viable seeds to germinate under optimal environmental conditions. To achieve rapid and uniform germination, wheat and other cereal crops have been selected against dormancy. As a result, most of the modern commercial cultivars have low level of seed dormancy and are susceptible to preharvest sprouting when wet and moist conditions occur prior to harvest. As it causes substantial loss in grain yield and quality, preharvest sprouting is an ever-present major constraint to the production of wheat. The significance of the problem emphasizes the need to incorporate an intermediate level of dormancy into elite wheat cultivars, and this requires detailed dissection of the mechanisms underlying the regulation of seed dormancy and preharvest sprouting. Seed dormancy research in wheat often involves after-ripening, a period of dry storage during which seeds lose dormancy, or comparative analysis of seeds derived from dormant and non-dormant cultivars. The increasing development in wheat genomic resources along with the application of transcriptomics, proteomics, and metabolomics approaches in studying wheat seed dormancy have extended our knowledge of the mechanisms acting at transcriptional and post-transcriptional levels. Recent progresses indicate that some of the molecular mechanisms are associated with hormonal pathways, epigenetic regulations, targeted oxidative modifications of seed mRNAs and proteins, redox regulation of seed protein thiols, and modulation of translational activities. Given that preharvest sprouting is closely associated with seed dormancy, these findings will significantly contribute to the designing of efficient strategies for breeding preharvest sprouting tolerant wheat.

Expression analysis of individual homoeologous wheat genome- and rye genome-specific transcripts in a 2BS.2RL wheat-rye translocation

Wheat-rye translocations are widely used in wheat breeding to confer resistance against abiotic and biotic stress. Studying gene expression in wheat-rye translocations is complicated due to the presence of homoeologous genes in hexaploid wheat and high levels of synteny between wheat and rye chromatin. To distinguish transcripts expressed from each of the three wheat genomes and those from rye chromatin, genomic probes generated from diploid progenitors of wheat and rye were synthesized on a custom array. A total of 407 transcripts showed homoeologous genome ('A', 'B' or 'D' genome)- or rye genome ('R')-specific differential expression, based on unequal values of probe hybridization. In a 2BS.2RL wheat-rye translocation, thirteen of the 407 transcripts showed preferential expressions from rye chromatin. As well as quantifying variation in homoeologous transcript in wheat-rye translocations, this study also provides a potential aid to examine the contribution of the subgenomes to complex allohexapolyploids.

Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures

Polyploidy has occurred throughout the evolutionary history of plants and led to diversification and plant ecological adaptation. Functional plasticity of duplicate genes is believed to play a major role in the environmental adaptation of polyploids. In this context, we characterized genome-wide homoeologous gene expression in Coffea arabica, a recent allopolyploid combining two subgenomes that derive from two closely related diploid species, and investigated its variation in response to changing environment. The transcriptome of leaves of C. arabica cultivated at different growing temperatures suitable for one or the other parental species was examined using RNA-sequencing. The relative contribution of homoeologs to gene expression was estimated for 9959 and 10,628 genes in warm and cold conditions, respectively. Whatever the growing conditions, 65% of the genes showed equivalent levels of homoeologous gene expression. In 92% of the genes, relative homoeologous gene expression varied < 10% between growing temperatures. The subgenome contributions to the transcriptome appeared to be only marginally altered by the different conditions (involving intertwined regulations of homeologs) suggesting that C. arabica's ability to tolerate a broader range of growing temperatures than its diploid parents does not result from differential use of homoeologs.

Overexpression of three TaEXPA1 homoeologous genes with distinct expression divergence in hexaploid wheat exhibit functional retention in Arabidopsis

Common wheat is a hexaploid species with most of the genes present as triplicate homoeologs. Expression divergences of homoeologs are frequently observed in wheat as well as in other polyploid plants. However, little is known about functional variances among homologous genes arising from polyploidy. Expansins play diverse roles in plant developmental processes related to the action of cell wall loosening. Expression of the three TaEXPA1 homoeologs varied dynamically at different stages and organs, and epigenetic modifications contribute to the expression divergence of three TaEXPA1 homoeologs during wheat development. Nevertheless, their functions remain to be clarified. We found that over expression of TaEXPA1-A, -B and -D produced similar morphological changes in transgenic Arabidopsis plants, including increased germination and growth rate during seedling and adult stages, indicating that the proteins encoded by these three wheat TaEXPA1 homoeologs have similar (or conserved) functions in Arabidopsis. Collectively, our present study provided an example of a set of homoeologous genes expression divergence in different developmental stages and organs in hexaploid wheat but functional retention in transgenic Arabidopsis plants.

Epigenetic modification contributes to the expression divergence of three TaEXPA1 homoeologs in hexaploid wheat (Triticum aestivum)

Common wheat is a hexaploid species with most of the genes present as triplicate homoeologs. Expression divergences of homoeologs are frequently observed in wheat, as well as in other polyploid plants. However, the mechanisms underlying this phenomenon are poorly understood. Expansin genes play important roles in the regulation of cell size, as well as organ size. We found that all three TaEXPA1 homoeologs were silenced in seedling roots. In seedling leaves, TaEXPA1-A and TaEXPA1-D were expressed, but TaEXPA1-B was silenced. Further analysis revealed that silencing of TaEXPA1-B in leaves occurred after the formation of the hexaploid. Chromatin immunoprecipitation assays revealed that the transcriptional silencing of three TaEXPA1 homoeologs in roots was correlated with an increased level of H3K9 dimethylation and decreased levels of H3K4 trimethylation and H3K9 acetylation. Reactivation of TaEXPA1-A and TaEXPA1-D expression in leaves was correlated with increased levels of H3K4 trimethylation and H3K9 acetylation, and decreased levels of H3K9 dimethylation in their promoters, respectively. Moreover, a higher level of cytosine methylation was detected in the promoter region of TaEXPA1-B, which may contribute to its silencing in leaves. We demonstrated that epigenetic modifications contribute to the expression divergence of three TaEXPA1 homoeologs during wheat development.

Benzoxazinoids in rye allelopathy - from discovery to application in sustainable weed control and organic farming

The allelopathic potency of rye (Secale cereale L.) is due mainly to the presence of phytotoxic benzoxazinones-compounds whose biosynthesis is developmentally regulated, with the highest accumulation in young tissue and a dependency on cultivar and environmental influences. Benzoxazinones can be released from residues of greenhouse-grown rye at levels between 12 and 20 kg/ha, with lower amounts exuded by living plants. In soil, benzoxazinones are subject to a cascade of transformation reactions, and levels in the range 0.5-5 kg/ha have been reported. Starting with the accumulation of less toxic benzoxazolinones, the transformation reactions in soil primarily lead to the production of phenoxazinones, acetamides, and malonamic acids. These reactions are associated with microbial activity in the soil. In addition to benzoxazinones, benzoxazolin-2(3H)-one (BOA) has been investigated for phytotoxic effects in weeds and crops. Exposure to BOA affects transcriptome, proteome, and metabolome patterns of the seedlings, inhibits germination and growth, and can induce death of sensitive species. Differences in the sensitivity of cultivars and ecotypes are due to different species-dependent strategies that have evolved to cope with BOA. These strategies include the rapid activation of detoxification reactions and extrusion of detoxified compounds. In contrast to sensitive ecotypes, tolerant ecotypes are less affected by exposure to BOA. Like the original compounds BOA and MBOA, all exuded detoxification products are converted to phenoxazinones, which can be degraded by several specialized fungi via the Fenton reaction. Because of their selectivity, specific activity, and presumably limited persistence in the soil, benzoxazinoids or rye residues are suitable means for weed control. In fact, rye is one of the best cool season cover crops and widely used because of its excellent weed suppressive potential. Breeding of benzoxazinoid resistant crops and of rye with high benzoxazinoid contents, as well as a better understanding of the soil persistence of phenoxazinones, of the weed resistance against benzoxazinoids, and of how allelopathic interactions are influenced by cultural practices, would provide the means to include allelopathic rye varieties in organic cropping systems for weed control.