Genome-wide identification of barley MCs (metacaspases) and their possible roles in boron-induced programmed cell death

Effects of foliar iron oxide nanoparticles (Fe3O4) application on photosynthetic parameters, distribution of mineral elements, magnetic behaviour, and photosynthetic genes in tomato (Solanum lycopersicum var. cerasiforme) plants

This study aims to examine the effect of foliar magnetic iron oxide (Fe3O4) nanoparticles (IONP) application on the physiology, photosynthetic parameters, magnetic character, and mineral element distribution of cherry tomatoes (Solanum lycopersicum var. cerasiforme). The IONP suspension (500 mg L-1) was sprayed once (S1), twice (S2), thrice (S3), and four times (S4) a week on seedlings. Upon 21 days of the treatments, photosynthetic parameters (chlorophyll, carotenoids, photosynthetic yield, electron transport rate) were elucidated. Inductively-coupled plasma-optical emission spectrometer (ICP-OES) and vibrating sample magnetometer (VSM) were used to determine the mineral elements and abundance of magnetic power in the seedlings. In addition, the RT-qPCR method was performed to quantify the expressions of photosystem-related (PsaC, PsbP6, and PsbQ) and ferritin-coding (Fer-1 and Fer-2) genes. Results revealed that the physiological and photosynthetic indices were improved upon S1 treatment. The optimal dosage of IONP spraying enhances chlorophyll, carotenoid, electron transport rate (ETR), and effective photochemical quantum yield of photosystem II (Y(II)) but substantially diminishes non-photochemical quenching (NPQ). However, frequent IONP applications (S2, S3, and S4) caused growth retardation and suppressed the photosynthetic parameters, suggesting a toxic effect of IONP in recurrent treatments. Fer-1 and Fer-2 expressions were strikingly increased by IONP applications, suggesting an attempt to neutralize the excess amount of Fe ions by ferritin. Nevertheless, frequent IONP treatment fluctuated the mineral distribution and caused growth inhibition. Although low-repeat foliar applications of IONP (S1 in this study) may help improve plant growth, consecutive applications (S2, S3, and S4) should be avoided.

Uptake and bioaccumulation of iron oxide nanoparticles (Fe3O4) in barley (Hordeum vulgare L.): effect of particle-size

Root-to-shoot translocation of nanoparticles (NPs) is a matter of interest due to their possible unprecedented effects on biota. Properties of NPs, such as structure, surface charge or coating, and size, determine their uptake by cells. This study investigates the size effect of iron oxide (Fe3O4) NPs on plant uptake, translocation, and physiology. For this purpose, Fe3O4 NPs having about 10 and 100 nm in average sizes (namely NP10 and NP100) were hydroponically subjected to barley (Hordeum vulgare L.) in different doses (50, 100, and 200 mg/L) at germination (5 days) and seedling (3 weeks) stages. Results revealed that particle size does not significantly influence the seedlings' growth but improves germination. The iron content in root and leaf tissues gradually increased with increasing NP10 and NP100 concentrations, revealing their root-to-shoot translocation. This result was confirmed by vibrating sample magnetometry analysis, where the magnetic signals increased with increasing NP doses. The translocation of NPs enhanced chlorophyll and carotenoid contents, suggesting their contribution to plant pigmentation. On the other hand, catalase activity and H2O2 production were higher in NP10-treated roots compared to NP100-treated ones. Besides, confocal microscopy revealed that NP10 leads to cell membrane damages. These findings showed that Fe3O4 NPs were efficiently taken up by the roots and transported to the leaves regardless of the size factor. However, small-sized Fe3O4 NPs may be more reactive due to their size properties and may cause cell stress and membrane damage. This study may help us better understand the size effect of NPs in nanoparticle-plant interaction.

Impact of magnetic field on the translocation of iron oxide nanoparticles (Fe3O4) in barley seedlings (Hordeum vulgare L.)

The effect and contribution of an external magnetic field (MF) on the uptake and translocation of nanoparticles (NPs) in plants have been investigated in this study. Barley was treated with iron oxide NPs (Fe3O4, 500 mg/L, 50-100 nm) and grown under various MF strengths (20, 42, 125, and 250 mT). The root-to-shoot translocation of NPs was assessed using a vibrating sample magnetometer (VSM) and inductively coupled plasma optical emission spectrometry (ICP-OES). Additionally, plant phenological parameters, such as germination, protein and chlorophyll content, and photosynthetic and nutritional status, were examined. The results demonstrated that the external MF significantly enhances the uptake of NPs through the roots. The uptake was higher at lower MF strengths (20 and 42 mT) than at higher MF strengths (125 and 250 mT). The root and shoot iron (Fe) contents were approximately 2.5-3-fold higher in the 250 mT application compared to the control. Furthermore, the MF treatments significantly increased micro-elements such as Mn, Zn, Cu, Mo, and B (P < 0.005). This effect could be attributed to the disruption of cell membranes at the root tip cells caused by both the MF and NPs. Moreover, the MF treatments improved germination rates by 28%, total protein content, and photosynthetic parameters. These findings show that magnetic field application helps the effective transport of magnetic NPs, which could be essential for NPs-mediated drug delivery, plant nutrition, and genetic transformation applications.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03727-4.

Regulation of Phenolic Compound Production by Light Varying in Spectral Quality and Total Irradiance

Photosynthetically active radiation (PAR) is an important environmental cue inducing the production of many secondary metabolites involved in plant oxidative stress avoidance and tolerance. To examine the complex role of PAR irradiance and specific spectral components on the accumulation of phenolic compounds (PheCs), we acclimated spring barley (Hordeum vulgare) to different spectral qualities (white, blue, green, red) at three irradiances (100, 200, 400 µmol m⁻² s⁻¹). We confirmed that blue light irradiance is essential for the accumulation of PheCs in secondary barley leaves (in UV-lacking conditions), which underpins the importance of photoreceptor signals (especially cryptochrome). Increasing blue light irradiance most effectively induced the accumulation of B-dihydroxylated flavonoids, probably due to the significantly enhanced expression of the F3′H gene. These changes in PheC metabolism led to a steeper increase in antioxidant activity than epidermal UV-A shielding in leaf extracts containing PheCs. In addition, we examined the possible role of miRNAs in the complex regulation of gene expression related to PheC biosynthesis.

Near-Isogenic Barley Lines Show Enhanced Susceptibility to Powdery Mildew Infection Following High-Temperature Stress

Barley cultivation is adversely affected by high-temperature stress, which may modulate plant defense responses to pathogens such as barley powdery mildew (Blumeria graminis f. sp. hordei, Bgh). Earlier research focused mainly on the influence of short-term heat stress (heat shock) of barley on Bgh infection. In this study, our aim was to investigate the effects of both short- and long-term heat stress (35 °C from 30 s to 5 days) on Bgh infection in the barley cultivar Ingrid and its near-isogenic lines containing different powdery mildew resistance genes (Mla12, Mlg, and mlo5) by analyzing symptom severity and Bgh biomass with RT-qPCR. The expression of selected barley defense genes (BAX inhibitor-1, Pathogenesis- related protein-1b, Respiratory burst oxidase homologue F2, and Heat shock protein 90-1) was also monitored in plants previously exposed to heat stress followed by inoculation with Bgh. We demonstrated that pre-exposure to short- and long-term heat stress negatively affects the resistance of all resistant lines manifested by the appearance of powdery mildew symptoms and increased Bgh biomass. Furthermore, prolonged heat stress (48 and 120 h) enhanced both Bgh symptoms and biomass in susceptible wild-type Ingrid. Heat stress suppressed and delayed early defense gene activation in resistant lines, which is a possible reason why resistant barley became partially susceptible to Bgh.

Magnetic field effects on the magnetic properties, germination, chlorophyll fluorescence, and nutrient content of barley (Hordeum vulgare L.)

The magnetic field (MF) interacts with biological systems and has the potential to increase germination, plant growth and productivity. Although it is known as a low cost and promising approach, the mechanism that increases growth is not fully understood yet. In this study, the effect of different MF strengths (20, 42, 125, and 250 mT) was investigated on barley (Hordeum vulgare L.). In addition to phenological parameters, possible cell damage, electron transport rate, chlorophyll fluorescence, magnetic character and elemental status of tissues were determined. Results showed that lower strengths (≤125 mT) of MF treatment improve germination. Confocal microscopy analyzes revealed MF-induced cell membrane damage in roots that could alter the elemental content of tissues. Elemental analyzes found that the content of macroelements (Ca, Mg, P, and K) are gradually reduced with increasing MF forces; in opposite the microelement contents (Fe, B, Cu, Mn, Zn, and Mo) are increased in roots. Diamagnetism is the dominant magnetic character in all root and leaf samples. However, the roots became surprisingly superparamagnetic in 250 mT application. It seems that MF treatment at higher strength (250 mT in this study) could influence the orientation of magnetic moments. These findings suggest that MF application: i) can alter the magnetic character of plants, ii) enhances the germination, photosynthetic machinery, and growth, and iii) affects the nutrient uptake and abundance in tissues, depending on the MF strength. This comprehensive study can help in understanding the interaction of magnetic field with plants.

Programmed Cell Death in Developing Brachypodium distachyon Grain

The normal developmental sequence in a grass grain entails the death of several maternal and filial tissues in a genetically regulated process termed programmed cell death (PCD). The progression and molecular aspects of PCD in developing grains have been reported for domesticated species such as barley, rice, maize and wheat. Here, we report a detailed investigation of PCD in the developing grain of the wild model species Brachypodium distachyon. We detected PCD in developing Brachypodium grains using molecular and histological approaches. We also identified in Brachypodium the orthologs of protease genes known to contribute to grain PCD and surveyed their expression. We found that, similar to cereals, PCD in the Brachypodium nucellus occurs in a centrifugal pattern following anthesis. However, compared to cereals, the rate of post-mortem clearance in the Brachypodium nucellus is slower. However, compared to wheat and barley, mesocarp PCD in Brachypodium proceeds more rapidly in lateral cells. Remarkably, Brachypodium mesocarp PCD is not coordinated with endosperm development. Phylogenetic analysis suggests that barley and wheat possess more vacuolar processing enzymes that drive nucellar PCD compared to Brachypodium and rice. Our expression analysis highlighted putative grain-specific PCD proteases in Brachypodium. Combined with existing knowledge on grain PCD, our study suggests that the rate of nucellar PCD moderates grain size and that the pattern of mesocarp PCD influences grain shape.

Evolutionary Understanding of Metacaspase Genes in Cultivated and Wild Oryza Species and Its Role in Disease Resistance Mechanism in Rice

Metacaspases (MCs), a class of cysteine-dependent proteases found in plants, fungi, and protozoa, are predominately involved in programmed cell death processes. In this study, we identified metacaspase genes in cultivated and wild rice species. Characterization of metacaspase genes identified both in cultivated subspecies of Oryza sativa, japonica, and indica and in nine wild rice species was performed. Extensive computational analysis was conducted to understand gene structures, phylogenetic relationships, cis-regulatory elements, expression patterns, and haplotypic variations. Further, the haplotyping study of metacaspase genes was conducted using the whole-genome resequencing data publicly available for 4726 diverse genotype and in-house resequencing data generated for north-east Indian rice lines. Sequence variations observed among wild and cultivated rice species for metacaspase genes were used to understand the duplication and neofunctionalization events. The expression profiles of metacaspase genes were analyzed using RNA-seq transcriptome profiling in rice during different developmental stages and stress conditions. Real-time quantitative PCR analysis of candidate metacaspase genes in rice cultivars Pusa Basmati-1 in response to Magnaporthe oryzae infection indicated a significant role in the disease resistance mechanism. The information provided here will help to understand the evolution of metacaspases and their role under stress conditions in rice.

Pumpkin rootstock improves the growth and development of watermelon by enhancing uptake and transport of boron and regulating the gene expression

Boron (B) is an essential trace element that plays a vital role in metabolic and physiological functions of higher plants. The adequate supply of B is important for plant growth and development. Grafting is a technique used to improve the ion uptake and plant growth. In this study, a commercial watermelon cultivar "Zaojia 8424" [Citrullus lanatus (Thunb.) Matsum. and Nakai.] was grafted onto pumpkin (Cucurbita maxima × Cucurbita moschata) rootstock cv. "Qingyan Zhenmu No.1" with an aim to investigate the response of grafted plants to different levels of B supply (0.25 μM, 25 μM and 75 μM B) in the nutrient solution. Self-grafted watermelon plants were used as control. Pumpkin rootstock improved the plant growth, chlorophyll and carotenoid contents, photosynthetic assimilation, stomatal conductance, transpiration rate, B accumulation and up-regulated the expression of NIP5;1, NIP6;1 and B transporter (BOR2, BOR4) genes in the roots and leaves at 25 μM B compared with self-grafted watermelon plants. Moreover, pumpkin rootstock reduced the oxidative stress and cell damage by reducing H₂O₂ and MDA contents, and down-regulating the expression of PDCD2-1, PDCD2-2 genes. Moreover, it enhanced the antioxidant activity of watermelon by up-regulating the expression of SOD1, SOD2, CAT2-1, and CAT2-2 genes. Based on these observations, we concluded that pumpkin rootstock has ability to improve the plant growth of watermelon by enhancing the B uptake. This study may help adjust the B concentration in the nutrient medium for watermelon production where pumpkin grafted plants are utilized.

Engineered magnetic nanoparticles enhance chlorophyll content and growth of barley through the induction of photosystem genes

This study investigates the impact of an engineered magnetic nanoparticle (MNP) on a crop plant. For this purpose, a sonochemical synthetic approach was utilized in order to dope magnetic elements (Co and Nd) into technologically important iron oxide NPs. After being characterized by using TEM, SEM, and XRD instruments, the MNPs were hydroponically applied to barley plants with varying doses (from 125 to 1000 mg/L) both in germination (4 days) and early growing stages (3 weeks). Physiological responses, as well as expression of photosystem marker genes, were assessed. Compared to the untreated control, MNP treatment enhanced germination rate (~ 31%), tissue growth (8% in roots, 16% in shoots), biomass (~ 21%), and chlorophyll (a, b) (~ 20%), and carotenoids (~ 22%) pigments. In general, plants showed the highest growth enhancement at 125 or 250 mg/L treatment. However, higher doses diminished the growth indices. Compared to the control, the catalase activity was significantly reduced in the leaves (~ 33%, p < 0.005) but stimulated in the roots (~ 46%, p < 0.005). All tested photosystem marker genes (BCA, psbA, and psaA) were overexpressed in MNP-treated leaves than non-treated control. Moreover, the gene expressions were found to be proportionally increased with increasing MNP doses, indicating a positive correlation between MNPs and the photosynthetic machinery, which could contribute to the enhancement of plant growth.

Plausible association between drought stress tolerance of barley (Hordeum vulgare L.) and programmed cell death via MC1 and TSN1 genes

Studying the drought-responsive transcriptome is of high interest as it can serve as a blueprint for stress adaptation strategies. Despite extensive studies in this area, there are still many details to be uncovered, such as the importance of each gene involved in the stress response as well as the relationship between these genes and the physiochemical processes governing stress tolerance. This study was designed to address such important details and to gain insights into molecular responses of barley (Hordeum vulgare L.) to drought stress. To that, we combined RNA-seq data analysis with field and greenhouse drought experiments in a systems biology approach. RNA-sequence analysis identified a total of 665 differentially expressed genes (DEGs) belonging to diverse functional categories. A gene network was derived from the DEGs, which comprised of a total of 131 nodes and 257 edges. Gene network topology analysis highlighted two programmed cell death (PCD) modulating genes, MC1 (metacaspase 1) and TSN1 (Tudor-SN 1), as important (hub) genes in the predicted network. Based on the field trial, a drought-tolerant and a drought-susceptible barley genotype was identified from eight tested cultivars. Identified genotypes exhibited different physiochemical characteristics, including proline content, chlorophyll concentration, percentage of electrolyte leakage and malondialdehyde content as well as expression profiles of MC1 and TSN1 genes. Machine learning and correspondence analysis revealed a significant relationship between drought tolerance and measured characteristics in the context of PCD. Our study provides new insights which bridge barley drought tolerance to PCD through MC1 and TSN1 pathway.

Identification of key genes and transcription factors in ageing Arabidopsis papilla cells by transcriptome analysis

Programmed cell death (PCD) play essential roles in plant growth and development. Stigmatic papilla cells form an indispensable organ for plant reproduction. The lifetime of papilla cells is tightly controlled, and the developmental PCD (dPCD) process is involved in papilla cell death. Hence, papilla cell death is a good model for studying on PCD process. In this study, the dPCD signal was visualized in dying papilla cells by detecting the GUS signal of the PCD-related reporter gene BIFUNCTIONAL NUCLEASE 1 (BFN1). We found that the GUS was not expressed at young stage, but strongly expressed in papilla cells at the ageing stage, indicating the PCD process was triggered to terminate the papilla cell fate. Given this, the RNA-Seq data set, which covered the information of the whole lifespan of papilla cells, was analyzed aiming to understand which genes and pathways were involved in papilla cell death. 37 differential expressed genes (DEGs) were isolated. Moreover, the pathways related to energy production and transportation, autophagy, and plant hormone signal transduction were considered as the key pathways involved in the papilla cell death. 9 types, total of 104 transcriptional factors (TFs) were identified as well. Finally, a putative working model of papilla cell death was integrated. The findings herein will enrich the knowledge of the dPCD-mediated pathway in regulating plant organ/tissue growth, development, senescence, and death. Our study will provide some referential gene resources for studying on the dPCD in other plant organs or tissues.

Barley long non-coding RNAs (lncRNA) responsive to excess boron

Long non-coding RNA (lncRNA) has a misleading name, since although they do not encode proteins, they may encode small peptides. Such transcripts are emerging as regulatory molecules. With the advent of next-generation sequencing technologies and novel bioinformatics tools, a tremendous amount of lncRNAs have been identified in several plant species. Recent reports demonstrated roles of plant lncRNAs such as development and environmental response. Here, we reported a genome-wide discovery of ~8000 barley lncRNAs and measured their expression pattern upon excessive boron (B) treatment. According to the tissue-based comparison, leaves have a greater number of B-responsive differentially expressed lncRNAs than the root. Functional annotation of the coding transcripts, which were co-expressed with lncRNAs, revealed that molecular function of the ion transport, establishment of localization, and response to stimulus significantly enriched only in the leaf. On the other hand, 32 barley endogenous target mimics (eTM) as lncRNAs, which potentially decoy the transcriptional suppression activity of 18 miRNAs, were obtained. Also, six lncRNAs, differentially expressed upon B-treatment, were selected and quantitatively analyzed in both B-sensitive and B-tolerant cultivars treated by excess B-level. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis confirmed the B-responsive expressional changes obtained by RNA sequencing. Notably, some lncRNAs (i.e., TCONS_00045190 and TCONS_00056415) over-expressed only in B-tolerant cultivar upon excess B treatment. Presented data including identification, expression measurement, and functional characterization of barley lncRNAs suggest that B-stress response might also be regulated by lncRNA expression, via cooperative interaction of miRNA-eTM-coding target transcript modules.

Uptake and translocation of magnetite (Fe3O4) nanoparticles and its impact on photosynthetic genes in barley (Hordeum vulgare L.)

This study investigates the fate and impact of iron oxide or magnetite (Fe₃O₄, ∼13 nm in size) nanoparticles (NPs) in barley (Hordeum vulgare L.), a common crop cultivated around the world. Barley seedlings were grown in hydroponic culture for three weeks to include NPs (125, 250, 500, and 1000 mg/L). Transmission electron microscopy (TEM) and vibrating sample magnetometer (VSM) techniques were used to assess their uptake and translocation. Photosynthesis marker genes were quantified by RT-qPCR. Results revealed that increasing doses of Fe₃O₄ NPs were gradually enhanced the plant growth up to 500 mg/L, which promoted the fresh weight (FW) respectively ∼19% and ∼88% for leaf and root tissues than the ones for control. No phytotoxic effect was recorded even at high NPs doses. NPs inclusion increased some phenological parameters such as chlorophyll, total soluble protein, number of chloroplasts, and dry weight. High NPs doses dramatically reduced the catalase activity and hydrogen peroxide content, suggesting a possible function of NPs as nanozyme in vivo. TEM observations showed that Fe₃O₄ NPs penetrated and internalized in the root cells. In leaves, they were mostly existed at the surrounding cell wall, suggesting their translocation from root to shoot without cellular penetration. Further analysis by using VSM confirmed the existence of Fe₃O₄ NPs in leaves which result in dramatic alterations of the photosystem genes (PetA, psaA, BCA and psbA). In conclusion, barley plants uptake and translocate Fe₃O₄ NPs, which promoted the plant growth probably due to the promoted gene expression and efficient photosynthetic activity.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

Metacaspase gene family in Rosaceae genomes: Comparative genomic analysis and their expression during pear pollen tube and fruit development

Metacaspase (MC), which is discovered gene family with distant caspase homologs in plants, fungi, and protozoa, may be involved in programmed cell death (PCD) processes during plant development and respond abiotic and biotic stresses. To reveal the evolutionary relationship of MC gene family in Rosaceae genomes, we identified 8, 7, 8, 12, 12, and 23 MC genes in the genomes of Fragaria vesca, Prunus mume, Prunus persica, Pyrus communis, Pyrus bretschneideri and Malus domestica, respectively. Phylogenetic analysis suggested that the MC genes could be grouped into three clades: Type I*, Type I and Type II, which was supported by gene structure and conserved motif analysis. Microsynteny analysis revealed that MC genes present in the corresponding syntenic blocks of P. communis, P. bretschneideri and M. domestica, and further suggested that large-scale duplication events play an important role in the expansion of MC gene family members in these three genomes than other Rosaceae plants (F. vesca, P. mume and P. persica). RNA-seq data showed the specific expression patterns of PbMC genes in response to drought stress. The expression analysis of MC genes demonstrated that PbMC01 and PbMC03 were able to be detected in all four pear pollen tubes and seven fruit development stages. The current study highlighted the evolutionary relationship and duplication of the MC gene family in these six Rosaceae genomes and provided appropriate candidate genes for further studies in P. bretschneideri.

Comparative Genomic and Transcriptomic Analysis Suggests the Evolutionary Dynamic of GH3 Genes in Gramineae Crops

Glycoside hydrolase 3 (GH3) gene family belongs to auxin-responsive gene families and is tightly linked with hormone homeostasis and signaling pathways. However, our knowledge about the evolutionary dynamic of GH3 genes in Gramineae crops is limited. In this study, a comparative genomic and transcriptomic analysis was conducted to study evolutionary patterns and the driving selective forces of GH3 gene family in six representative Gramineae crops, namely, Setaria italica (Si), Zea mays (Zm), Sorghum bicolor (Sb), Hordeum vulgare (Hv), Brachypodium distachyon (Bd), and Oryza sativa ssp. japonica (Os). A total of 17, 13, 11, 9, 8, and 11 GH3 proteins (GH3s) were identified in Si, Zm, Sb, Hv, Bd, and Os, respectively. Phylogenetic, conserved motif, and gene structural analyses could divide all GH3s into two groups (I and II), and all GH3s consisted of seven orthogroups (Ors) on the basis of Or identification result. We further found that genes in the same Or showed similar sequence and structural features, whereas genes in the same groups exhibited intrinsic differences in exon numbers and intron lengths. These results revealed GH3 genes in the same groups have been differentiated. Obvious differences in total numbers of GH3 genes, Ors, and duplication events among these six tested Gramineae crops reflected lineage-specific expansions and homologous gene loss/gain of GH3 gene family during the evolutionary process. In addition, selective force and expression analyses indicated that all GH3 genes were constrained by strong purifying selection, and GH3 genes in conserved Ors showed higher expression levels than that in unconserved Ors. The current study highlighted different evolutionary patterns of GH3 genes in Gramineae crops resulted from different evolutionary rates and duplication events and provided a vital insight into the functional divergence of GH3 genes.