Comprehensive expression analysis suggests overlapping and specific roles of rice glutathione S-transferase genes during development and stress responses

A novel OsGST gene encoding 9glutathione reductase negatively regulates cadmium accumulation in rice

Cadmium (Cd) accumulates in rice and then moves up the food chain, causing serious health problems for humans. Glutathione S-transferase (GST) binds exogenous hazardous compounds to glutathione (GSH), which performs a variety of roles in plant responses to Cd stress. Here, Cd stimulated the transcripts of a novel OsGST gene, and the OsGST protein, which was localized in the nucleus and cytoplasm, was also induced by Cd. In OsGST deletion mutant lines generated by CRISPR/Cas9, more Cd was accumulated, and Cd hypersensitive phenotypes were observed, while transgenic lines overexpressing OsGST exhibited enhanced Cd tolerance and less Cd accumulation. Further analysis indicated that the osgst mutants exhibited considerably greater reactive oxygen species (ROS) and higher GSH level, and the antioxidant activity associated genes' expression were down-regulated, imply that OsGST controlled rice Cd accumulation and resistance through preserving the equilibrium of the GSH and redox in rice.

Expression pattern of the poplar GSTU family members in response to Alternaria alternate and PdbGSTU10 confers A. alternate resistance to Populus davidiana × P. bolleana

Plant tau glutathione S-transferase (GSTU) is a kind of multiple functions enzyme, but its specific roles in poplar disease resistance remain uncertain. In this study, 27 PdbGSTU-encoding genes from Populus davidiana × P. bollena were cloned and their protein architectures and phylogenetic relationships were subsequently analyzed. Expression analysis revealed that PdbGSTUs were differentially expressed under Alternaria alternate infection. Then, the PdbGSTU10 was further induced by phytohormones and H2O2, especially salicylic acid (SA), indicating its potential role in the pathogen defense of poplar. Subsequently, gain- and loss-of-function assays showed that overexpressed PdbGSTU10 activated antioxidant enzymes and significantly decreased reactive oxygen species (ROS) content, ultimately improving the resistance to A. alternate in poplar. Conversely, silencing PdbGSTU10 had the opposite effect. Moreover, overexpressed PdbGSTU10 also increased the content of SA and induced the expression of SA signal-related genes. These results showed that PdbGSTU10 may enhance disease resistance in poplar by scavenging ROS and affecting the SA signaling pathway. Our findings contribute to the understanding of the functions of GSTU in woody plants, particularly in disease resistance.

Genome-Wide Identification of Glutathione S-Transferase Family from Dendrobium officinale and the Functional Characterization of DoGST5 in Cadmium Tolerance

Glutathione S-transferases (GSTs) are members of a protein superfamily with diverse physiological functions, including cellular detoxification and protection against oxidative damage. However, there is limited research on GSTs responding to cadmium (Cd) stress. This study classified 46 GST genes in Dendrobium officinale (D. officinale) into nine groups using model construction and domain annotation. Evolutionary analysis revealed nine subfamilies with diverse physical and chemical properties. Prediction of subcellular localization revealed that half of the GST members were located in the cytoplasm. According to the expression analysis of GST family genes responding to Cd stress, DoGST5 responded significantly to Cd stress. Transient expression of DoGST5-GFP in tobacco leaves revealed that DoGST5 was localized in the cytoplasm. DoGST5 overexpression in Arabidopsis enhanced Cd tolerance by reducing Cd-induced H2O2 and O2- levels. These findings demonstrate that DoGST5 plays a critical role in enhancing Cd tolerance by balancing reactive oxygen species (ROS) levels, offering potential applications for improving plant adaptability to heavy metal stress.

Cr(VI) behaves differently than Cr(III) in the uptake, translocation and detoxification in rice roots

Excessive accumulation of chromium (Cr) causes severe damage to both physiological and biochemical processes and consequently growth repression in plants. Hexavalent chromium [Cr(VI)]-elicited alterations in plants have been widely elucidated at either physiological or molecular level, whereas little is known about trivalent chromium [Cr(III)]. Here, we found that both Cr(III) and Cr(VI) significantly inhibited root growth in rice plants. However, rice plants under Cr(VI) showed significantly less inhibition in root growth than those under Cr(III) at low levels, which might be attributed to the different hormetic effects of Cr(III) and Cr(VI) on rice plants. It was unexpected that Cr(III) could be actively taken up by rice roots similarly to Cr(VI); whereas they exhibited different kinetic uptake patterns. Furthermore, root-to-shoot Cr translocation under Cr(VI) was much lower than that under Cr(III). These results indicate that the uptake, translocation, and toxicity of Cr(III) differed greatly from those of Cr(VI). Transcriptome profiling of rice roots revealed that a series of gene families involved in detoxification, including ATP-binding cassette (ABC) transporters, multidrug and toxic compound extrusion proteins (MATEs), and Tau class glutathione S-transferases (GSTUs), were significantly associated with Cr accumulation and detoxification in rice roots. In addition, much more members of these gene families were upregulated by Cr(VI) compared to Cr(III), suggesting their vital roles in Cr uptake, translocation, and detoxification, especially under Cr(VI) stress. Further comparison of gstu9 and gstu10/50 mutants with their wild type confirmed that GSTUs play complex roles in the intracellular Cr transport and redox homeostasis during Cr(III) or Cr(VI) stress. Taken together, our findings provides new insights into the differential behaviors of Cr(III) and Cr(VI) in rice roots, as well as new candidate genes such as OsABCs and OsGSTUs, to further elucidate the mechanisms of the uptake, translocation, and detoxification of Cr(III) and Cr(VI).

Genome-wide study of glutathione transferases and their regulation in flufenacet susceptible and resistant black-grass (Alopecurus myosuroides Huds.)

BACKGROUND

Glutathione transferases (GSTs) are enzymes with a wide range of functions, including herbicide detoxification. Up-regulation of GSTs and their detoxification activity enables the grass weed black-grass (Alopecurus myosuroides Huds.) to metabolize the very-long-chain fatty acid synthesis inhibitor flufenacet and other herbicides leading to multiple herbicide resistance. However, the genomic organization and regulation of GSTs genes is still poorly understood.

RESULTS

In this genome-wide study the location and expression of 115 GSTs were investigated using a recently published black-grass genome. Particularly, the most abundant GSTs of class tau and phi were typically clustered and often followed similar expression patterns but possessed divergent upstream regulatory regions. Similarities were found in the promoters of the most up-regulated GSTs, which are located next to each other in a cluster. The binding motif of the E2F/DP transcription factor complex in the promoter of an up-regulated GST was identical in susceptible and resistant plants, however, adjacent sequences differed. This led to a stronger binding of proteins to the motif of the susceptible plant, indicating repressor activity.

CONCLUSIONS

This study constitutes the first analysis dealing with the genomic investigation of GST genes found in black-grass and their transcriptional regulation. It highlights the complexity of the evolution of GSTs in black-grass, their duplication and divergence over time. The large number of GSTs allows weeds to detoxify a broad spectrum of herbicides. Ultimately, more research is needed to fully elucidate the regulatory mechanisms of GST expression. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Glutathione S-transferase activity facilitates rice tolerance to the barnyard grass root exudate DIMBOA

BACKGROUND

In paddy fields, the noxious weed barnyard grass secretes 2,4-dihydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one (DIMBOA) to interfere with rice growth. Rice is unable to synthesize DIMBOA. Rice cultivars with high or low levels of allelopathy may respond differently to DIMBOA.

RESULTS

In this study, we found that low concentrations of DIMBOA (≤ 0.06 mM) promoted seedling growth in allelopathic rice PI312777, while DIMBOA (≤ 0.08 mM) had no significant influence on the nonallelopathic rice Lemont. DIMBOA treatment caused changes in the expression of a large number of glutathione S-transferase (GST) proteins, which resulting in enrichment of the glutathione metabolic pathway. This pathway facilitates plant detoxification of heterologous substances. The basal levels of GST activity in Lemont were significantly higher than those in PI312777, while GST activity in PI312777 was slightly induced by increasing DIMBOA concentrations. Overexpression of GST genes (Os09g0367700 and Os01g0949800) in these two cultivars enhanced rice resistance to DIMBOA.

CONCLUSIONS

Taken together, our results indicated that different rice accessions with different levels of allelopathy have variable tolerance to DIMBOA. Lemont had higher GST activity, which helped it tolerate DIMBOA, while PI312777 had lower GST activity that was more inducible. The enhancement of GST expression facilitates rice tolerance to DIMBOA toxins from barnyard grass root exudates.

DcGST1, encoding a glutathione S-transferase activated by DcMYB7, is the main contributor to anthocyanin pigmentation in purple carrot

The color of purple carrot taproots mainly depends on the anthocyanins sequestered in the vacuoles. Glutathione S-transferases (GSTs) are key enzymes involved in anthocyanin transport. However, the precise mechanism of anthocyanin transport from the cytosolic surface of the endoplasmic reticulum (ER) to the vacuoles in carrots remains unclear. In this study, we conducted a comprehensive analysis of the carrot genome, leading to the identification of a total of 41 DcGST genes. Among these, DcGST1 emerged as a prominent candidate, displaying a strong positive correlation with anthocyanin pigmentation in carrot taproots. It was highly expressed in the purple taproot tissues of purple carrot cultivars, while it was virtually inactive in the non-purple taproot tissues of purple and non-purple carrot cultivars. DcGST1, a homolog of Arabidopsis thaliana TRANSPARENT TESTA 19 (TT19), belongs to the GSTF clade and plays a crucial role in anthocyanin transport. Using the CRISPR/Cas9 system, we successfully knocked out DcGST1 in the solid purple carrot cultivar 'Deep Purple' ('DPP'), resulting in carrots with orange taproots. Additionally, DcMYB7, an anthocyanin activator, binds to the DcGST1 promoter, activating its expression. Compared with the expression DcMYB7 alone, co-expression of DcGST1 and DcMYB7 significantly increased anthocyanin accumulation in carrot calli. However, overexpression of DcGST1 in the two purple carrot cultivars did not change the anthocyanin accumulation pattern or significantly increase the anthocyanin content. These findings improve our understanding of anthocyanin transport mechanisms in plants, providing a molecular foundation for improving and enhancing carrot germplasm.

Translocation, enzymatic reduction and toxicity of dimethylarsenate in rice

Dimethylarsenate [DMAs(V)] can be produced by some soil microorganisms through methylation of inorganic arsenic (As), especially in anoxic paddy soils. DMAs(V) is more phytotoxic than inorganic As and can cause the physiological disorder straighthead disease in rice. Rice cultivars vary widely in the resistance to DMAs(V), but the mechanism remains elusive. Here, we investigated the differences in DMAs(V) uptake, translocation, and reduction to dimethylarsenite [DMAs(III)], as well as the effects on the metabolome, between two rice cultivars Mars and Zhe733. We found that Mars was 11-times more resistant to DMAs(V) than Zhe733. Mars accumulated more DMAs(V) in the roots, whereas Zhe733 translocated more DMAs(V) to the shoots and reduced more DMAs(V) to DMAs(III). DMAs(III) was more toxic than DMAs(V). Using heterologous expression and in vitro enzyme assays, we showed that the glutathione-S-transferases OsGSTU17 and OsGSTU50 were able to reduce DMAs(V) to DMAs(III). The expression levels of OsGSTU17 and OsGSTU50 were higher in the shoot of Zhe733 compared to Mars. Metabolomic analysis in rice shoots showed that glutathione (GSH) metabolism was perturbed by DMAs(V) toxicity in Zhe733. Application of exogenous GSH significantly alleviated the toxicity of DMAs(V) in Zhe733. Taken together, the results suggest that Mars is more resistant to DMAs(V) than Zhe733 because of a lower root-to-shoot translocation and a smaller capacity to reduce DMAs(V) to DMAs(III).

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Functional significance of serine 13 in the active site of glutathione S-transferase F3 from Oryza sativa

Plant glutathione S-transferase (GST, EC 2.5.1.18) is an enzyme that detoxifies various electrophilic compounds including herbicides and organic pollutants by catalyzing the formation of conjugates with reduced glutathione (GSH). Although the structure and function of the GST subunits in rice, an important food in Asia, are not well understood, they are crucial for herbicide development. To investigate the role of active site residues in rice Phi-class GSTF3 (OsGSTF3), evolutionarily conserved serine residues were replaced with alanine using site-directed mutagenesis to obtain the mutants S13A, S38A, S69A, and S169A. These four mutants were expressed in Escherichia coli and purified to electrophoretic homogeneity using immobilized GSH affinity chromatography. Mutation of Ser13 to Ala resulted in substantial reductions in specific activities and kcat/Km values for the GSH-[1-chloro-2,4-dinitrobenzene (CDNB)] conjugation reaction. In contrast, mutations of Ser38, Ser69, and Ser169 to Ala had little effect on the activities and kinetic parameters. Additionally, the mutation of Ser13 to Ala significantly affected the KmGSH and I50 values of S-hexylglutathione and S-(2,4-dinitrophenyl)glutathione, which compete with GSH and the product of GSH-CDNB conjugation, respectively. A pH-log (kcat/KmCDNB) plot was used to estimate the pKa value of GSH in the enzyme-GSH complex of the wild-type enzyme, which was approximately 6.9. However, the pKa value of GSH in the enzyme-GSH complex of the S13A mutant was approximately 8.7, which was about 1.8 pK units higher than that of the wild-type enzyme. OsGSTF3 was also crystallized for crystallographic study, and the structure analyses revealed that Ser13 is located in the active site and that its side chain is in close proximity to the thiol group of glutathione bound in the enzyme. Based on these substitution effects on kinetic parameters, the dependence of kinetic parameters on the pH and 3-dimensional structure, it was suggested that Ser13 in rice OsGSTF3 is the residue responsible for catalytic activity by lowering the pKa of GSH in the enzyme-GSH complex and enhancing the nucleophilicity of the GSH thiol in the active site.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

Genome-wide identification and characterization of Glutathione S-Transferases (GSTs) and their expression profile under abiotic stresses in tobacco (Nicotiana tabacum L.)

BACKGROUND

Glutathione S-transferases (GSTs) are large and multifunctional proteases that play an important role in detoxification, protection against biotic and abiotic stresses, and secondary metabolite transportation which is essential for plant growth and development. However, there is limited research on the identification and function of NtGSTs.

RESULTS

This study uses K326 and other six tobacco varieties (Hongda, HG, GDH11, Va116, VG, and GDH88) as materials to conduct comprehensive genome-wide identification and functional characterization of the GST gene in tobacco. A total of 59 NtGSTs were identified and classified into seven subfamilies via the whole-genome sequence analysis, with the Tau type serving as the major subfamily. The NtGSTs in the same branch of the evolutionary tree had similar exon/intron structure and motif constitution. There were more than 42 collinear blocks between tobacco and pepper, tomato, and potato, indicating high homology conservation between them. Twelve segmental duplicated gene pairs and one tandem duplication may have had a substantial impact on the evolution and expansion of the tobacco GST gene family. The RT-qPCR results showed that the expression patterns of NtGSTs varied significantly among tissues, varieties, and multiple abiotic stresses, suggesting that NtGST genes may widely respond to various abiotic stresses and hormones in tobacco, including NtGSTF4, NtGSTL1, NtGSTZ1, and NtGSTU40.

CONCLUSIONS

This study provides a comprehensive analysis of the NtGST gene family, including structures and functions. Many NtGSTs play a critical regulatory role in tobacco growth and development, and responses to abiotic stresses. These findings offer novel and valuable insights for understanding the biological function of NtGSTs and the reference materials for cultivating highly resistant varieties and enhancing the yield and quality of crops.

iTRAQ based proteomic analysis of rice lines having single or stacked blast resistance genes: Pi54/Pi54rh during incompatible interaction with Magnaporthe oryzae

Deployment of single or multiple blast resistance (R) genes in rice plant is considered to be the most promising approach to enhance resistance against blast disease caused by fungus Magnaporthe oryzae. At the proteome level, relatively little information about R gene mediated defence mechanisms for single and stacking resistance characteristics is available. The overall objective of this study is to look at the proteomics of rice plants that have R genes; Pi54, Pi54rh and stacked Pi54 + Pi54rh in response to rice blast infection. In this study 'isobaric tag for relative and absolute quantification' (iTRAQ)-based proteomics analysis was performed in rice plants at 72-h post inoculation with Magnaporthe oryzae and various differentially expressed proteins were identified in these three transgenic lines in comparison to wild type during resistance response to blast pathogen. Through STRING analysis, the observed proteins were further examined to anticipate their linked partners, and it was shown that several defense-related proteins were co-expressed. These proteins can be employed as targets in future rice resistance breeding against Magnaporthe oryzae. The current study is the first to report a proteomics investigation of rice lines that express single blast R gene Pi54, Pi54rh and stacked (Pi54 + Pi54rh) during incompatible interaction with Magnaporthe oryzae. The differentially expressed proteins indicated that secondary metabolites, reactive oxygen species-related proteins, phenylpropanoid, phytohormones and pathogenesis-related proteins have a substantial relationship with the defense response against Magnaporthe oryzae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01327-3.

Understanding the response in Pugionium cornutum (L.) Gaertn. seedling leaves under drought stress using transcriptome and proteome integrated analysis

Background

Drought is one of the crucial constraints limiting horticultural plant’s production and development around the world. Pugionium cornutum is an annual or biennial xerophyte with strong environmental adaptability and drought resistance; however, the mechanisms with respect to response to drought stress remain largely unclear.

Methods

After seedling emergence, the gravimetric method was used to control soil relative water content (SRWC). Drought stress was applied to the six-leaf stage P. cornutum seedlings. The soil water content of different drought stress levels (L) was controlled by gravimetric method as follows: control (L1): 70–75% SRWC; moderate drought level (L2): 40–45% SRWC; severe drought level (L3): 30–35% SRWC, and the water was added to different drought stress levels at about 18:00 p.m. every day. The experiment ended when the leaves of P. cornutum showed severe wilting (10-leaf stage). Samples were harvested and stored at −80 °C for physiological determination, and transcriptomic and proteomic sequencing.

Results

Compared with L1, the leaves of P. cornutum seedlings were increasingly wilted after drought treatment; the SRWC of the drought-stress leaves decreased notably while the leaf water potential was rose; the proline, malondialdehyde (MDA) content increased with the continuous drought treatment but peroxidase (POD) activity decreased. Besides, 3,027 differential genes (DGs) and 196 differential proteins (DPs), along with 1,943 DGs and 489 DPs were identified in L2-L1 and L3-L1, respectively. The transcriptome and proteome integrated analysis manifested that only 30 and 70 were commonly regulated both in L2-L1 and L3-L1, respectively. Of which, 24 and 61 DGs or DPs showed the same trend including sHSPs, APX2, GSTU4, CML42, and POD, etc. However, most of DGs or DPs were regulated only at the transcriptome or proteome level mainly including genes encoding signal pathway (PYR1, PYLs, SnRK2J, PLC2, CDPK9/16/29, CML9, MAPKs), transcription factors (WRKYs, DREB2A, NAC055, NAC072, MYB and, HB7) and ion channel transporters (ALMT4, NHX1, NHX2 and TPK2). These genes or proteins were involved in multiple signaling pathways and some important metabolism processes, which offers valuable information on drought-responsive genes and proteins for further study in P. cornutum.

PgMYB1 Positively Regulates Anthocyanin Accumulation by Activating PgGSTF6 in Pomegranate

The peel color of pomegranates is an important exterior quality that determines market value. Anthocyanins are biosynthesized in the cytosol and then transported to the vacuole for storage. However, the molecular mechanism that determines the color variation between red and white pomegranates remains unclear. In this study, we identified an R2R3-MYB protein (PgMYB1) that interacts with the PgGSTF6 promoter and regulates its transcriptional expression, thus promoting the accumulation of anthocyanins in pomegranate. The expression of PgMYB1 and PgGSTF6 was positively correlated with the anthocyanin content in red and white pomegranates. Further investigation showed that the knockdown of PgMYB1 in red pomegranate ‘Taishanhong’ (TSH), by the virus-induced gene-silencing system, inhibited anthocyanin accumulation. Together, our results indicate that PgMYB1 controls the transport of anthocyanin via PgGSTF6 and thus promotes anthocyanin accumulation in red pomegranates. Our results have a certain reference value for further clarifying the regulation of anthocyanin synthesis and transport in pomegranate fruits.

Genome-Wide Identification and Expression Profiling of Glutathione S-Transferase Gene Family in Foxtail Millet (Setaria italica L.)

Glutathione S-transferases (GSTs) are a critical superfamily of multifunctional enzymes in plants. As a ligand or binding protein, GSTs regulate plant growth and development and detoxification. Foxtail millet (Setaria italica (L.) P. Beauv) could respond to abiotic stresses through a highly complex multi-gene regulatory network in which the GST family is also involved. However, GST genes have been scarcely studied in foxtail millet. Genome-wide identification and expression characteristics analysis of the foxtail millet GST gene family were conducted by biological information technology. The results showed that 73 GST genes (SiGSTs) were identified in the foxtail millet genome and were divided into seven classes. The chromosome localization results showed uneven distribution of GSTs on the seven chromosomes. There were 30 tandem duplication gene pairs belonging to 11 clusters. Only one pair of SiGSTU1 and SiGSTU23 were identified as fragment duplication genes. A total of ten conserved motifs were identified in the GST family of foxtail millet. The gene structure of SiGSTs is relatively conservative, but the number and length of exons of each gene are still different. The cis-acting elements in the promoter region of 73 SiGST genes showed that 94.5% of SiGST genes possessed defense and stress-responsive elements. The expression profiles of 37 SiGST genes covering 21 tissues suggested that most SiGST genes were expressed in multiple organs and were highly expressed in roots and leaves. By qPCR analysis, we found that 21 SiGST genes were responsive to abiotic stresses and abscisic acid (ABA). Taken together, this study provides a theoretical basis for identifying foxtail millet GST family information and improving their responses to different stresses.

Genes encoding cytochrome P450 monooxygenases and glutathione S-transferases associated with herbicide resistance evolved before the origin of land plants

Cytochrome P450 (CYP) monooxygenases and glutathione S-transferases (GST) are enzymes that catalyse chemical modifications of a range of organic compounds. Herbicide resistance has been associated with higher levels of CYP and GST gene expression in some herbicide-resistant weed populations compared to sensitive populations of the same species. By comparing the protein sequences of 9 representative species of the Archaeplastida-the lineage which includes red algae, glaucophyte algae, chlorophyte algae, and streptophytes-and generating phylogenetic trees, we identified the CYP and GST proteins that existed in the common ancestor of the Archaeplastida. All CYP clans and all but one land plant GST classes present in land plants evolved before the divergence of streptophyte algae and land plants from their last common ancestor. We also demonstrate that there are more genes encoding CYP and GST proteins in land plants than in algae. The larger numbers of genes among land plants largely results from gene duplications in CYP clans 71, 72, and 85 and in the GST phi and tau classes [1,2]. Enzymes that either metabolise herbicides or confer herbicide resistance belong to CYP clans 71 and 72 and the GST phi and tau classes. Most CYP proteins that have been shown to confer herbicide resistance are members of the CYP81 family from clan 71. These results demonstrate that the clan and class diversity in extant plant CYP and GST proteins had evolved before the divergence of land plants and streptophyte algae from a last common ancestor estimated to be between 515 and 474 million years ago. Then, early in embryophyte evolution during the Palaeozoic, gene duplication in four of the twelve CYP clans, and in two of the fourteen GST classes, led to the large numbers of CYP and GST proteins found in extant land plants. It is among the genes of CYP clans 71 and 72 and GST classes phi and tau that alleles conferring herbicide resistance evolved in the last fifty years.

Genome-wide in silico identification of glutathione S-transferase (GST) gene family members in fig (Ficus carica L.) and expression characteristics during fruit color development

Glutathione S-transferase (GSTs), a large and diverse group of multi-functional enzymes (EC 2.5.1.18), are associated with cellular detoxification, various biotic and abiotic stress responses, as well as secondary metabolites transportation. Here, 53 members of the FcGST gene family were screened from the genome database of fig (Ficus carica), which were further classified into five subfamilies, and the tau and phi were the major subfamilies. These genes were unevenly distributed over all the 13 chromosomes, and 12 tandem and one segmental duplication may contribute to this family expansion. Syntenic analysis revealed that FcGST shared closer genetic evolutionary origin relationship with species from the Ficus genus of the Moraceae family, such as F. microcarpa and F. hispida. The FcGST members of the same subfamily shared similar gene structure and motif distribution. The α helices were the chief structure element in predicted secondary and tertiary structure of FcGSTs proteins. GO and KEGG indicated that FcGSTs play multiple roles in glutathione metabolism and stress reactions as well as flavonoid metabolism. Predictive promoter analysis indicated that FcGSTs gene may be responsive to light, hormone, stress stimulation, development signaling, and regulated by MYB or WRKY. RNA-seq analysis showed that several FcGSTs that mainly expressed in the female flower tissue and peel during 'Purple-Peel' fig fruit development. Compared with 'Green Peel', FcGSTF1, and FcGSTU5/6/7 exhibited high expression abundance in the mature fruit purple peel. Additionally, results of phylogenetic sequences analysis, multiple sequences alignment, and anthocyanin content together showed that the expression changes of FcGSTF1, and FcGSTU5/6/7 may play crucial roles in fruit peel color alteration during fruit ripening. Our study provides a comprehensive overview of the GST gene family in fig, thus facilitating the further clarification of the molecular function and breeding utilization.

Computational insights into diverse aspects of glutathione S-transferase gene family in Papaver somniferum

Plant glutathione S-transferases are an ancient protein superfamily having antioxidant activity. These proteins are primarily involved in diverse plant functions such as plant growth and development, secondary metabolism, signaling pathways and defense against biotic and abiotic stresses. The current study aimed to comprehensively identify and characterize the GST gene family in the medicinally important crop Papaver somniferum. A total of 93 GST proteins were identified belonging to eight GST classes and found to be majorly localized in the cytoplasm. All GST genes were found on eleven opium chromosomes. Gene duplication analysis showed segmental duplication as a key factor for opium GST gene family expansion under strong purifying selection. Phylogenetic analysis with gymnosperm, angiosperm and bryophyte revealed the evolution of GSTs earlier than their division into separate groups and also prior to the divergence of monocot and dicot. The secondary structure prediction showed the dominance of α-helices indicative of PsomGSTs as structurally stable and elastic proteins. Gene architecture showed the conservation of number of exons across the classes. MEME analysis revealed only a few class specific and many across class conserved motifs. Ser was found to be the active site residue of tau, phi, theta and zeta class and Cys was catalytic residue of DHAR, lambda and GHR class. Promoter analyses identified many cis-acting regulatory elements related to hormonal, cellular, stress and light response functions. Ser was the key phosphorylation site. Only three glycosylation sites were found across the 93 PsomGSTs. 3D structure prediction was also performed and was validated. Interactome analyses revealed the correlation of PsomGSTs with glutathione metabolizing proteins. Gene enrichment analysis and KEGG pathway analyzed the involvement of PsomGSTs in three major pathways i.e. glutathione metabolism, tyrosine metabolism and ascorbate metabolism. The outcome revealed high model quality of PsomGSTs. The results of the current study will be of potential significance to understand the functional and structural importance of the GST gene family in opium, a medicinally important crop.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

Insight into the toxic effects, bioconcentration and oxidative stress of acetamiprid on Rana nigromaculata tadpoles

Pesticide pollution has been identified as a factor in the amphibian population decrease. Acetamiprid is a common neonicotinoid pesticide that poses a risk to amphibians due to its high water solubility and inability to be digested. However, there is little research on acetamiprid's toxicity in amphibians, particularly on its biochemical toxic effects. In this study, we investigated the acute toxicity, bioenrichment-elimination, biochemical parameters and metabolism of acetamiprid in Rana nigromaculata tadpoles. The results indicated that acetamiprid is harmful to Rana nigromaculata tadpoles, with an LC₅₀ = 18.49 mg L^-1 of 96 h for acute toxicity. Acetamiprid showed rapid accumulation and low bioconcentration levels in tadpoles, with bioconcentration factors (BCFs) < 1. In the elimination process, the concentration of acetamiprid decreased rapidly, with the elimination half-life t_1/2 values < 1 d. Additionally, oxidative stress was observed in tadpoles, with biochemical parameters such as superoxide dismutase (SOD), catalase (CAT) and malondialdehyde (MDA) being significantly altered. Nontargeted metabolomics revealed significant changes in biomolecules such as lipids, organic acids and nucleotides in tadpoles, and these metabolites influence pathways including serine and threonine metabolism, histidine metabolism, linoleic acid metabolism and sphingolipid metabolism. These results indicate that acetamiprid caused toxic effects on Rana nigromaculata tadpoles. Our study provides a better understanding of the fate and risk of acetamiprid in amphibians, as well as guidelines for its rational use.

Effects of modified clay on the formation of Phaeocystis globosa colony revealed by physiological and transcriptomic analyses

The harmful algal bloom (HAB) species Phaeocystis globosa is commonly observed in global temperate and tropical oceans, and colonies of P. globosa exhibit a dominant morphotype during blooms. The use of polyaluminium chloride modified clay (PAC-MC) is an effective mitigation strategy for P. globosa blooms. Although previous studies have found that PAC-MC can stimulate P. globosa colony formation at low concentrations and inhibit it at higher concentrations, the underlying mechanisms of these effects are poorly understood. Here, we comprehensively compared the physiochemical indices and transcriptomic response of residual P. globosa cells after treatment with two concentrations of PAC-MC. The results showed that PAC-MC induced oxidative stress, photosynthetic inhibition, and DNA damage in residual cells. Moreover, it could activate antioxidant responses and enhance the repair of photosynthetic structure and DNA damage in cells. The biosynthesis of polysaccharides was enhanced and genes associated with cell motility were down-regulated after treatment with PAC-MC, resulting in the accumulation of colonial matrixes. After treatment with a low concentration of PAC-MC (0.1 g/L), the residual cells were slightly stressed, including physical damage, oxidative stress and other damage, and polysaccharide synthesis was enhanced to promote colony formation to alleviate environmental stress. Moreover, the damage to residual cells was slight; thus, normal cell function provided abundant energy and matter for colony formation. After treatment with a high concentration of PAC-MC (0.5 g/L), the residual cells suffered severe damage, which disrupted normal physiological processes and inhibited cell proliferation and colony formation. The present study elucidated the concentration-dependent mechanism of PAC-MC affecting the formation of P. globosa colonies and provided a reference for the application of PAC-MC to control P. globosa blooms.

Biochemical and Structural Insights on the Poplar Tau Glutathione Transferase GSTU19 and 20 Paralogs Binding Flavonoids

Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in xenobiotic detoxification and/or in specialized metabolism. Populus trichocarpa genome (V4.1 assembly, Phytozome 13) consists of 74 genes coding for full-length GSTs and ten likely pseudogenes. These GSTs are divided into 11 classes, in which the tau class (GSTU) is the most abundant with 54 isoforms. PtGSTU19 and 20, two paralogs sharing more than 91% sequence identity (95% of sequence similarity), would have diverged from a common ancestor of P. trichocarpa and P. yatungensis species. These enzymes display the distinctive glutathione (GSH)-conjugation and peroxidase activities against model substrates. The resolution of the crystal structures of these proteins revealed significant structural differences despite their high sequence identity. PtGSTU20 has a well-defined deep pocket in the active site whereas the bottom of this pocket is disordered in PtGSTU19. In a screen of potential ligands, we were able to identify an interaction with flavonoids. Some of them, previously identified in poplar (chrysin, galangin, and pinocembrin), inhibited GSH-conjugation activity of both enzymes with a more pronounced effect on PtGSTU20. The crystal structures of PtGSTU20 complexed with these molecules provide evidence for their potential involvement in flavonoid transport in P. trichocarpa.

Genome-wide identification and characterization of glutathione S-transferase gene family in Musa acuminata L. AAA group and gaining an insight to their role in banana fruit development

Glutathione S-transferases are a multifunctional protein superfamily that is involved in diverse plant functions such as defense mechanisms, signaling, stress response, secondary metabolism, and plant growth and development. Although the banana whole-genome sequence is available, the distribution of GST genes on banana chromosomes, their subcellular localization, gene structure, their evolutionary relation with each other, conserved motifs, and their roles in banana are still unknown. A total of 62 full-length GST genes with the canonical thioredoxin fold have been identified belonging to nine GST classes, namely tau, phi, theta, zeta, lambda, DHAR, EF1G, GHR, and TCHQD. The 62 GST genes were distributed into 11 banana chromosomes. All the MaGSTs were majorly localized in the cytoplasm. Gene architecture showed the conservation of exon numbers in individual GST classes. Multiple Em for Motif Elicitation analyses revealed few class-specific motifs and many motifs were found in all the GST classes. Multiple sequence alignment of banana GST amino acid sequences with rice, Arabidopsis, and soybean sequences revealed the Ser and Cys as conserved catalytic residues. Gene duplication analyses showed the tandem duplication as a driving force for GST gene family expansion in banana. Cis-regulatory element analysis showed the dominance of light-responsive element followed by stress- and hormone-responsive elements. Expression profiling analyses were also done by RNA-seq data. It was observed that MaGSTs are involved in various stages of fruit development. MaGSTU1 was highly upregulated. The comprehensive and organized studies of MaGST gene family provide groundwork for further functional analysis of MaGST genes in banana at molecular level and further for plant breeding approaches.

An Anthocyanin-Related Glutathione S-Transferase, MrGST1, Plays an Essential Role in Fruit Coloration in Chinese Bayberry (Morella rubra)

Chinese bayberry (Morella rubra) is a fruit tree economically important in China and accumulates abundant amounts of anthocyanins in fruit as it ripens. Owing to the fact that all anthocyanin containing fruit tissues in Chinese bayberry are edible and anthocyanins can provide various health benefits in human body, the mechanisms underpinning anthocyanin accumulation in this fruit are worthy of investigation. It has been known that in plants anthocyanins are synthesized in the cytoplasmic surface of the endoplasmic reticulum and subsequently transported into the vacuole for storage, and glutathione S-transferases (GSTs) have been verified to be involved in this process. But the characterization and functionalization of the GST counterpart in Chinese bayberry is not available. The GST anthocyanin transporter MrGST1 was discovered to be related with anthocyanin accumulation in fruit from distinct developmental stages of "Biqi," a staple cultivar that accumulates over 1 mg/g anthocyanins in ripe fruit. The expression of MrGST1 was well associated with anthocyanin accumulation either in fruit collected at six developmental stages or in ripe fruit from 12 cultivars. MrGST1 was found to be responsible for the transport of anthocyanins but not proanthocyanidins when the Arabidopsis tt19 mutant was functionally complemented. Transient ectopic expression of MrGST1 in combination with MrMYB1.1 and MrbHLH1 dramatically boosted pigmentation in Nicotiana tabacum leaves in contrast to MrMYB1.1 and MrbHLH1. The promoter of MrGST1 comprised eight MYB binding sites (MBSs) according to cis-element analysis. Data from yeast one-hybrid assay and dual-luciferase tests demonstrated that MrMYB1.1 exerted considerable transactivation effect on the MrGST1 promoter by recognizing the MBS4, the fourth MBS from the ATG start site. Our results together provided molecular evidence for the contribution of MrGST1 in regulating anthocyanin accumulation in Chinese bayberry fruit.

LSD3 mediates the oxidative stress response through fine-tuning APX2 activity and the NF-YC15-GSTs module in cassava

Reactive oxygen species (ROS) overproduction leads to oxidative damage under almost all stress conditions. Lesion-Simulating Disease (LSD), a zinc finger protein, is an important negative regulator of ROS accumulation and cell death in plants. However, the in vivo role of LSD in cassava (Manihot esculenta) and the underlying molecular mechanisms remain elusive. Here, we found that MeLSD3 is essential for the oxidative stress response in cassava. MeLSD3 physically interacted with ascorbate peroxidase 2 (MeAPX2), thereby promoting its enzymatic activity. In addition, MeLSD3 also interacted with the nuclear factor YC15 (MeNF-YC15), which also interacted with nuclear factor YA2/4 (MeNF-YA2/4) and nuclear factor YB18 (MeNF-YB18) to form an MeNF-YC15-MeNF-YA2/4-MeNF-YB18 complex. Notably, MeLSD3 positively modulated the transcriptional activation of the MeNF-YC15-MeNF-YA2/4-MeNF-YB18 complex by interacting with the CCAAT boxes of the promoters of glutathione S-transferases U37/U39 (MeGST-U37/U39), activating their transcription. When one or both of MeLSD3 and the MeNF-YC15-MeNF-YA2/4-MeNF-YB18 complex were co-silenced, cassava showed decreased oxidative stress resistance, while overexpression of MeGST-U37/U39 alleviated the oxidative stress-sensitive phenotype of these silenced plants. This study illustrates the dual roles of MeLSD3 in promoting MeAPX2 activity and MeNF-YC15-MeGST-U37/U39 regulation, which underlie the oxidative stress response in cassava.

A tau class GST, OsGSTU5, interacts with VirE2 and modulates the Agrobacterium-mediated transformation in rice

OsGSTU5 interacts and glutathionylates the VirE2 protein of Agrobacterium and its (OsGSTU5) overexpression and downregulation showed a low and high AMT efficiency in rice, respectively. During Agrobacterium-mediated transformation (AMT), T-DNA along with several virulence proteins such as VirD2, VirE2, VirE3, VirD5, and VirF enter the plant cytoplasm. VirE2 serves as a single-stranded DNA binding (SSB) protein that assists the cytoplasmic trafficking of T-DNA inside the host cell. Though the regulatory roles of VirE2 have been established, the cellular reaction of their host, especially in monocots, has not been characterized in detail. This study identified a cellular interactor of VirE2 from the cDNA library of rice. The identified plant protein encoded by the gene cloned from rice was designated OsGSTU5, it interacted specifically with VirE2 in the host cytoplasm. OsGSTU5 was upregulated during Agrobacterium infection and involved in the post-translational glutathionylation of VirE2 (gVirE2). Interestingly, the in silico analysis showed that the 'gVirE2 + ssDNA' complex was structurally less stable than the 'VirE2 + ssDNA' complex. The gel shift assay also confirmed the attenuated SSB property of gVirE2 over VirE2. Moreover, knock-down and overexpression of OsGSTU5 in rice showed increased and decreased T-DNA expression, respectively after Agrobacterium infection. The present finding establishes the role of OsGSTU5 as an important target for modulation of AMT efficiency in rice.

Assessing the toxic effects of bisphenol A in consumed crayfish Astacus leptodactylus using multi biochemical markers

Bisphenol A (BPA), an endocrine-disrupting chemical (EDC), has strong potential for daily exposure to humans and animals due to its persistence and widespread in the environment, so its effects directly concern public health. Although invertebrates represent important components of aquatic ecosystems and are at significant risk of exposure, there is little information about the biological effects of EDCs in these organisms. Astacus leptodactylus used in this study is one of the most consumed and exported freshwater species in Europe. In this study, the 96-h effect of BPA on A. leptodactylus was examined using various biomarkers. The LC₅₀ value of BPA was determined as 96.45 mg L^-1. After 96 h of exposure to BPA, there were increases in superoxide dismutase (SOD), glutathione peroxidase (GPx), glutathione S-transferase (GST), alanine aminotransferase (ALT), aspartate aminotransferase (AST) and lactate dehydrogenase (LDH) activities and levels of malondialdehyde (MDA), and total oxidant status context (TOSC), and there were decreases in the activity of glutathione reductase (GR), carboxylesterase (CaE), acetylcholinesterase (AChE), Na⁺/K⁺ ATPase, Mg²⁺ ATPase, Ca²⁺ ATPase, and total ATPase and the total antioxidant context (TAC). From the results of this study, it can be concluded that BPA has significant toxic effects on A. leptodactylus based on the selected biochemical parameters of antioxidant, cholinergic, detoxification, and metabolic systems in crayfish even at low doses. Thus, it can be said that BPA can seriously threaten the aquatic ecosystem and public health.

A Sweet Cherry Glutathione S-Transferase Gene, PavGST1, Plays a Central Role in Fruit Skin Coloration

Sweet cherry, an economically important horticultural crop, has strong antioxidant activity. The fruits contain compounds potentially beneficial to human health—particularly anthocyanins, which are synthesized in cytosol and predominantly accumulated in vacuoles. Although anthocyanin levels differ among dark-red, blush, and yellow sweet cherry cultivars, the regulatory mechanism of anthocyanin transport and accumulation is not well understood in this species. In this study, we identified 53 glutathione S-transferase genes (PavGSTs) from sweet cherry and found that PavGST1 expression was well correlated with anthocyanin accumulation in cultivars with different fruit skin colors. TRV-mediated virus-induced silencing of PavGST1 decreased anthocyanin accumulation in sweet cherry fruits and downregulated the expressions of anthocyanin biosynthetic and regulatory genes. In addition, transient overexpression of PavGST1 promoted anthocyanin accumulation. Furthermore, yeast one-hybrid and dual-luciferase assays revealed that PavMYB10.1 and PavMYB75 directly bind to different MYB binding sites of the PavGST1 promoter (MBS-1 and MBS-3) to activate PavGST1 transcription. According to our results, PavGST1 plays a central role in sweet cherry fruit anthocyanin accumulation. Our findings provide novel insights into the coordinative regulatory mechanisms of PavGST1 and PavMYBs in anthocyanin accumulation in sweet cherry.

A tau class glutathione-S-transferase (OsGSTU5) confers tolerance against arsenic toxicity in rice by accumulating more arsenic in root

Arsenic (As) considered as one of the hazardous metalloid that hampers various physiological activities in rice. To study the mechanism of As tolerance in rice, one differentially expressed tau class glutathione-S-transferase (OsGSTU5) has been selected and transgenic rice plants with knockdown (KD) and overexpressing (OE) OsGSTU5 were generated. Our results suggested that KD lines became less tolerant to As stress than WT plants, while OE lines showed enhanced tolerance to As. Under As toxicity, OE and KD lines showed enhanced and reduced antioxidant activities such as, SOD, PRX and catalase, respectively indicating its role in ROS homeostasis. In addition, higher malondialdehyde content, poor photosynthetic parameters and higher reactive oxygen species (ROS) in KD plant, suggests that knockdown of OsGSTU5 renders KD plants more susceptible to oxidative damage. Also, the relative expression profile of various transporters such as OsABCC1 (As sequestration), Lsi2 and Lsi6 (As translocaters) and GSH dependent activity of GSTU5 suggests that GSTU5 might help in chelation of As with GSH and sequester it into the root vacuole using OsABCC1 transporter and thus limits the upward translocation of As towards shoot. This study suggests the importance of GSTU5 as a good target to improve the As tolerance in rice.

Genome-Wide Identification, Evolution and Expression Analysis of the Glutathione S-Transferase Supergene Family in Euphorbiaceae

Euphorbiaceae, a family of plants mainly grown in the tropics and subtropics, is also widely distributed all over the world and is well known for being rich in rubber, oil, medicinal materials, starch, wood and other economically important plant products. Glutathione S-transferases (GSTs) constitute a family of proteins encoded by a large supergene family and are widely expressed in animals, bacteria, fungi and plants, but with few reports of them in Euphorbiaceae plants. These proteins participate in and regulate the detoxification and oxidative stress response of heterogeneous organisms, resistance to stress, growth and development, signal transduction and other related processes. In this study, we identified and analyzed the whole genomes of four species of Euphorbiaceae, namely Ricinus communis, Jatropha curcas, Hevea brasiliensis, and Manihot esculenta, which have high economic and practical value. A total of 244 GST genes were identified. Based on their sequence characteristics and conserved domain types, the GST supergene family in Euphorbiaceae was classified into 10 subfamilies. The GST supergene families of Euphorbiaceae and Arabidopsis have been found to be highly conserved in evolution, and tandem repeats and translocations in these genes have made the greatest contributions to gene amplification here and have experienced strong purification selection. An evolutionary analysis showed that Euphorbiaceae GST genes have also evolved into new subtribes (GSTO, EF1BG, MAPEG), which may play a specific role in Euphorbiaceae. An analysis of expression patterns of the GST supergene family in Euphorbiaceae revealed the functions of these GSTs in different tissues, including resistance to stress and participation in herbicide detoxification. In addition, an interaction analysis was performed to determine the GST gene regulatory mechanism. The results of this study have laid a foundation for further analysis of the functions of the GST supergene family in Euphorbiaceae, especially in stress and herbicide detoxification. The results have also provided new ideas for the study of the regulatory mechanism of the GST supergene family, and have provided a reference for follow-up genetics and breeding work.

Genome-wide analysis of the Glutathione S-Transferase family in wild Medicago ruthenica and drought-tolerant breeding application of MruGSTU39 gene in cultivated alfalfa

Transformation of MruGSTU39 in M. ruthenica and alfalfa enhanced growth and survival of transgenic plants by up-regulating GST and glutathione peroxidase activity to detoxify ROS under drought stress. Glutathione S-transferases (GSTs) are ubiquitous supergene family which play crucial roles in detoxification of reactive oxygen species (ROS). Despite studies on GSTs, few studies have focused on them in perennial, wild plant species with high tolerance to environmental stress. Here, we identified 66 MruGST genes from the genome of Medicago ruthenica, a perennial legume species native to temperate grasslands with high tolerance to environmental stress. These genes were divided into eight classes based on their conserved domains, phylogenetic tree and gene structure, with the tau class being the most numerous. Duplication analysis revealed that GST family in M. ruthenica was expanded by segmental and tandem duplication. Several drought-responsive MruGSTs were identified by transcriptomic analyses. Of them, expression of MruGSTU39 was up-regulated much more in a tolerant accession by drought stress. Transformation of MruGSTU39 in M. ruthenica and alfalfa (Medicago sativa) enhanced growth and survival of transgenic seedlings than their wild-type counterparts under drought. We demonstrated that MruGSTU39 can detoxify ROS to reduce its damage to membrane by up-regulating activities of GST and glutathione peroxidase. Our findings provide full-scale knowledge on GST family in the wild legume M. ruthenica with high tolerance to drought, and highlight improvement tolerance of legume forages to drought using genomic information of M. ruthenica.

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

CsGSTU8, a Glutathione S-Transferase From Camellia sinensis, Is Regulated by CsWRKY48 and Plays a Positive Role in Drought Tolerance

Glutathione S-transferases (GSTs) constitute a large family of enzymes with a wide range of cellular functions. Recently, plant GSTs have gained a great deal of attention due to their involvement in the detoxification of electrophilic xenobiotics and peroxides under adverse environmental conditions, such as salt, cold, UV-B and drought stress. A previous study reported that a GST gene (CsGSTU8) in tea plant was distinctly induced in response to drought, suggesting this gene plays a critical role in the drought stress response. In this study, by using quantitative real-time PCR (qRT-PCR) and β-glucuronidase (GUS) reporter lines, we further demonstrated that CsGSTU8 was upregulated in response to drought stress and exogenous abscisic acid (ABA) treatments. Overexpression of CsGSTU8 in Arabidopsis resulted in enhanced drought tolerance as indicated by the improved scavenging of excess amounts of reactive oxygen species (ROS) under drought conditions. Furthermore, we found that CsWRKY48 acts as a transcriptional activator and that its expression is induced in response to drought stress and ABA treatment. Electrophoretic mobility shift assays (EMSAs), dual-luciferase (LUC) assays and transient expression assays in tea plant leaves revealed that CsWRKY48 directly binds to the W-box elements in the promoter of CsGSTU8 and activates its expression. Taken together, our results provide additional knowledge of drought stress responses in tea plant.

Acclimation to cold stress reduces injury from low temperature and bispyribac-sodium on rice

BACKGROUND

In subtropical areas, early planting exposes rice seedlings to cold stress, impairing seedling growth and making them more vulnerable to other stresses including herbicide injury. The objectives of this work were: to evaluate the effect of cold stress on bispyribac-sodium selectivity in rice; to determine the mechanisms of cold tolerance in sensitive ('Epagri 109') and tolerant ('IRGA 424') rice cultivars; and to ascertain that cold acclimatization influences bispyribac-sodium selectivity in rice.

RESULTS

Prolonged cold stress caused high lipid peroxidation, increased rice injury, and stunted growth. Short-term acclimation with cold stress reduced rice injury with bispyribac-sodium. Total phenols were upregulated in rice exposed to cold stress. Prolonged cold stress increased the superoxide dismutase and catalase activity in IRGA 424. Antioxidant activity was higher in the cold-tolerant than in the cold-sensitive cultivar. Only catalase activity was responsive to bispyribac-sodium. OsRAN2, OsGSTL2, and CYP72A21 were upregulated by cold and herbicide stress in both cultivars. OsGSTL2 was upregulated more in IRGA 424 than in Epagri 109. OsFAD8 was upregulated in cold-sensitive rice exposed to short-duration cold stress but was not responsive to bispyribac-sodium.

CONCLUSION

Cold stress reduces bispyribac-sodium selectivity in rice. Short-term acclimation to cold stress reduces the effect of cold stress and enhances bispyribac-sodium selectivity. The tolerance of rice (IRGA 424) to cold stress is due to differential induction of protection genes CYP72A21 and OsGSTL2 associated with herbicide metabolism, together with the accumulation of total phenols and higher activity of antioxidant enzymes.

Identification and functional characterization of RsGST1, an anthocyanin-related glutathione S-transferase gene in radish

Anthocyanins not only affect the quality of horticultural crops but are also vital for human health. Glutathione transferase family members (GSTs) are enzymes that help to control plant development and stress responses, and are also involved in anthocyanin accumulation. In this study, we targeted a phi (F) class glutathione S-transferase gene RsGST1 (RSG01330.t1) as a crucial gene in the accumulation of anthocyanins in radish. RsGST1 expression was closely associated with anthocyanin content in the skin and flesh of taproot from different color type radish cultivars. Furthermore, RsGST1 was able to restore anthocyanin accumulation in Arabidopsis tt19 mutants, indicating that RsGST1 has a similar function as AtTT19, a gene responsible for the transport of anthocyanins in Arabidopsis. Transient overexpression of RsGST1 together with the key anthocyanin biosynthesis regulator RsMYB1a in radish leaves significantly enhanced anthocyanin biosynthesis compared with in plants that overexpressed RsMYB1a alone. Dual-luciferase and yeast one-hybrid assays revealed that RsMYB1a binds to promotor and activates the expression of RsGST1. Altogether, these results provide molecular evidence indicating that RsGST1 and RsMYB1a coordinate radish anthocyanin accumulation.

Integrative Transcriptomic and Proteomic Analysis Reveals an Alternative Molecular Network of Glutamine Synthetase 2 Corresponding to Nitrogen Deficiency in Rice (Oryza sativa L.)

Nitrogen (N) is an essential nutrient for plant growth and development. The root system architecture is a highly regulated morphological system, which is sensitive to the availability of nutrients, such as N. Phenotypic characterization of roots from LY9348 (a rice variety with high nitrogen use efficiency (NUE)) treated with 0.725 mM NH₄NO₃ (1/4N) was remarkable, especially primary root (PR) elongation, which was the highest. A comprehensive analysis was performed for transcriptome and proteome profiling of LY9348 roots between 1/4N and 2.9 mM NH₄NO₃ (1N) treatments. The results indicated 3908 differential expression genes (DEGs; 2569 upregulated and 1339 downregulated) and 411 differential abundance proteins (DAPs; 192 upregulated and 219 downregulated). Among all DAPs in the proteome, glutamine synthetase (GS2), a chloroplastic ammonium assimilation protein, was the most upregulated protein identified. The unexpected concentration of GS2 from the shoot to the root in the 1/4N treatment indicated that the presence of an alternative pathway of N assimilation regulated by GS2 in LY9348 corresponded to the low N signal, which was supported by GS enzyme activity and glutamine/glutamate (Gln/Glu) contents analysis. In addition, N transporters (NRT2.1, NRT2.2, NRT2.3, NRT2.4, NAR2.1, AMT1.3, AMT1.2, and putative AMT3.3) and N assimilators (NR2, GS1;1, GS1;2, GS1;3, NADH-GOGAT2, and AS2) were significantly induced during the long-term N-deficiency response at the transcription level (14 days). Moreover, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that phenylpropanoid biosynthesis and glutathione metabolism were significantly modulated by N deficiency. Notably, many transcription factors and plant hormones were found to participate in root morphological adaptation. In conclusion, our study provides valuable information to further understand the response of rice roots to N-deficiency stress.

Identification of Glutathione S-Transferase Genes in Hami Melon (Cucumis melo var. saccharinus) and Their Expression Analysis Under Cold Stress

As a group of multifunctional enzymes, glutathione S-transferases (GSTs) participate in oxidative stress resistance and cellular detoxification. Here, we identified 39 CmGST genes with typical binding sites from the Hami melon genome, and they can be classified into seven subfamilies. Their molecular information, chromosomal locations, phylogenetic relationships, synteny relationships, gene structures, protein-protein interactions, structure of 3-D models, and expression levels under cold stress were analyzed. Expression analysis indicates that cold-tolerant Jia Shi-310 (JS) had higher GST enzyme activities and expression levels of 28 stress-related genes under cold stress. Some CmGSTs belonging to Tau, Phi, and DHAR classes play significant roles under cold stress, and they could be regarded as candidate genes for further studies. The present study systematically investigated the characterization of the Hami melon GST gene family, extending our understanding of Hami melon GST mediated stress-response mechanisms in this worldwide fruit.

Fumarylacetoacetate hydrolase is required for fertility in rice

The rice OsFAH gene functions identically to that of Arabidopsis SSCD1 encoding FAH. Loss of OsFAH causes rice sterility. Fumarylacetoacetate hydrolase (FAH) is the last enzyme in the tyrosine (Tyr) degradation pathway that is crucial for animals. By genetic analysis of the mutant of Short-day Sensitive Cell Death 1 gene encoding Arabidopsis FAH, we first found the pathway also plays a critical role in plants (Han et al., Plant Physiol 162:1956-1964, 2013). To further understand the role of the Tyr degradation pathway in plants, we investigated a biological function of the rice FAH. Firstly, the cDNA of rice FAH gene (OsFAH) was cloned and confirmed to be able to rescue the Arabidopsis Short-day Sensitive Cell Death 1 mutant defective in the FAH. Then, we identified the OsFAH T-DNA insertion mutant and generated the OsFAH RNA interference lines, and found that loss of OsFAH results in rice sterility. Furthermore, we analyzed expression of the OsFAH gene in roots, stems, leaves and young panicles at booting stage of rice and found that its transcript level was highest in young panicles and lowest in roots. In addition, the expression analysis of β-glucuronidase driven by OsFAH promoter in transgenic Arabidopsis showed that the OsFAH promoter was highly active in aerial tissues in vegetative stage, and sepals, filaments and stigma in reproductive stage. These results suggested that FAH plays an important role in rice fertility.

A comprehensive transcriptome analysis of contrasting rice cultivars highlights the role of auxin and ABA responsive genes in heat stress response

Sensing a change in ambient temperature is key to survival among all living organisms. Temperature fluctuations due to climate change are a matter of grave concern since it adversely affects growth and eventually the yield of crop plants, including two of the major cereals, i.e., rice and wheat. Thus, to understand the response of rice seedlings to elevated temperatures, we performed microarray-based transcriptome analysis of two contrasting rice cultivars, Annapurna (heat tolerant) and IR64 (heat susceptible), by subjecting their seedlings to 37 °C and 42 °C, sequentially. The transcriptome analyses revealed a set of uniquely regulated genes and related pathways in red rice cultivar Annapurna, particularly associated with auxin and ABA as a part of heat stress response in rice. The changes in expression of few auxin and ABA associated genes, such as OsIAA13, OsIAA20, ILL8, OsbZIP12, OsPP2C51, OsDi19-1 and OsHOX24, among others, were validated under high-temperature conditions using RT-qPCR. In particular, the expression of auxin-inducible SAUR genes was enhanced considerably at both elevated temperatures. Further, using genes that expressed inversely under heat vs. cold temperature conditions, we built a regulatory network between transcription factors (TF) such as HSFs, NAC, WRKYs, bHLHs or bZIPs and their target gene pairs and determined regulatory coordination in their expression under varying temperature conditions. Our work thus provides useful insights into temperature-responsive genes, particularly under elevated temperature conditions, and could serve as a resource of candidate genes associated with thermotolerance or downstream components of temperature sensors in rice.

Genetic variation and gene expression of anthocyanin synthesis and transport related enzymes in Oryza sativa against thiocyanate

Plants exposed to environmental contaminants often synthesize anthocyanins (ATHs) as an approach to safeguard themselves from adverse impact. However, the overload of ATHs in plant cells can threaten their growth and development through proteins oxidization and intercalating with DNAs inside cells. In the present study, a microcosm hydroponic experiment was conducted using rice seedlings to investigate the molecular signaling pathways involved in regulating and controlling ATHs synthesis and transport exposed to thiocyanate (SCN^-). Our results indicated that SCN^- exposure significantly (p < 0.05) increased the expression of ATHs synthesis related genes (i.e., PAL, CHS, ANS, UFGT genes) in rice tissues, altered the activities of these ATHs synthesis related enzymes, and consequently elevated the ATHs content. However, SCN^- exposure significantly decreased the expression of ATHs transport related genes (i.e., GST, ABC, MATE genes) in rice seedlings, suggesting that SCN^- exposure have restrained ATHs transport from cytosol to vacuole in cells, eventually posing a significant adverse effect on cells survival. Our findings highlight on one of the plant aspects in managing the toxicity triggered by secondary metabolites under stress conditions.

Genome-wide identification and expression profiling of glutathione S-transferase family under multiple abiotic and biotic stresses in Medicago truncatula L

Glutathione transferases (GSTs) constitute an ancient, ubiquitous, multi-functional antioxidant enzyme superfamily that has great importance on cellular detoxification against abiotic and biotic stresses as well as plant development and growth. The present study aimed to a comprehensive genome-wide identification and functional characterization of GST family in one of the economically important legume plants-Medicago truncatula. Here, we have identified a total of ninety-two putative MtGST genes that code for 120 proteins. All these members were classified into twelve classes based on their phylogenetic relationship and the presence of structural conserved domain/motif. Among them, 7 MtGST gene pairs were identified to have segmental duplication. Expression profiling of MtGST transcripts revealed their high level of organ/tissue-specific expression in most of the developmental stages and anatomical tissues. The transcripts of MtGSTU5, MtGSTU8, MtGSTU17, MtGSTU46, and MtGSTU47 showed significant up-regulation in response to various abiotic and biotic stresses. Moreover, transcripts of MtGSTU8, MtGSTU14, MtGSTU28, MtGSTU30, MtGSTU34, MtGSTU46 and MtGSTF8 were found to be highly upregulated in response to drought treatment for 24h and 48h. Among the highly stress-responsive MtGST members, MtGSTU17 showed strong affinity towards its conventional substrates reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB) with the lowest binding energy of-5.7 kcal/mol and -6.5 kcal/mol, respectively. Furthermore, the substrate-binding site residues of MtGSTU17 were found to be highly conserved. These findings will facilitate the further functional and evolutionary characterization of GST genes in Medicago.

Comparative Analysis of the Glutathione S-Transferase Gene Family of Four Triticeae Species and Transcriptome Analysis of GST Genes in Common Wheat Responding to Salt Stress

Glutathione S-transferases (GSTs) are ancient proteins encoded by a large gene family in plants, which play multiple roles in plant growth and development. However, there has been little study on the GST genes of common wheat (Triticum aestivum) and its relatives (Triticum durum, Triticum urartu, and Aegilops tauschii), which are four important species of Triticeae. Here, a genome-wide comprehensive analysis of this gene family was performed on the genomes of common wheat and its relatives. A total of 346 GST genes in T. aestivum, 226 in T. durum, 104 in T. urartu, and 105 in Ae. tauschii were identified, and all members were divided into ten classes. Transcriptome analysis was used to identify GST genes that respond to salt stress in common wheat, which revealed that the reaction of GST genes is not sensitive to low and moderate salt concentrations but is sensitive to severe concentrations of the stressor, and the GST genes related to salt stress mainly come from the Tau and Phi classes. Six GST genes which respond to different salt concentrations were selected and validated by a qRT-PCR assay. These findings will not only provide helpful information about the function of GST genes in Triticeae species but also offer insights for the future application of salt stress resistance breeding in common wheat.

LhGST is an anthocyanin-related glutathione S-transferase gene in Asiatic hybrid lilies (Lilium spp.)

LhGST, an anthocyanin-related GST gene, was identified from Asiatic hybrid lilies. Expression and functional analyses demonstrated that LhGST might be involved in anthocyanin sequestration in lily tepals. Anthocyanins are responsible for the pink, red and purple pigmentation of flowers in Asiatic hybrid lilies, synthesized at the cytoplasmic surface of the endoplasmic reticulum (ER) and then transported to the vacuole. To date, the mechanism involved in the intracellular transport of anthocyanins in lilies has not been well elucidated. Here, full-length glutathione S-transferase gene (LhGST) was identified from lilies. Expression analysis revealed that LhGST was positively correlated with anthocyanin accumulation. Phylogenetic tree analysis showed that LhGST clustered with other anthocyanin-related GSTs in the same phi clade. Moreover, functional complementation of an Arabidopsis tt19 mutant demonstrated that LhGST might be involved in anthocyanin accumulation in lily tepals. Additionally, according to phenotype analysis, LhGST was found to be correlated with the transport of anthocyanin in lilies by virus-induced gene silencing (VIGS) assay. In addition, cis-element analysis of the LhGST promoter showed the presence of ABA-, auxin-, MeJA-, gibberellin-, light-, and stress-responsive elements and an MYB recognition site (MRS, CCGTTG). Yeast one-hybrid and dual-luciferase report assays revealed that the promoter of LhGST was activated by LhMYB12-lat, which is a key R2R3-MYB transcription factor that regulates anthocyanin biosynthesis in lilies. In conclusion, our results revealed that LhGST plays a key role in anthocyanin transport and accumulation in the tepals of lilies.

Integrated analysis of co-expression, conserved genes and gene families reveal core regulatory network of heat stress response in Cleistogenes songorica, a xerophyte perennial desert plant

BACKGROUND

As global warming continues, heat stress (HS) is becoming an increasingly significant factor limiting plant growth and reproduction, especially for cool-season grass species. The objective of this study was to determine the transcriptional regulatory network of Cleistogenes songorica under HS via transcriptome profiling, identify of gene families and comparative analysis across major Poaceae species.

RESULTS

Physiological analysis revealed significantly decreased leaf relative water content (RWC) but increased proline (Pro) content in C. songorica under 24 h of HS. Transcriptome profiling indicated that 16,028 and 14,645 genes were differentially expressed in the shoots and roots of C. songorica under HS, respectively. Two subgenomes of C. songorica provide equal contribution under HS on the basis of the distribution and expression of differentially expressed genes (DEGs). Furthermore, 216 DEGs were identified as key evolutionarily conserved genes involved in the response to HS in C. songorica via comparative analysis with genes of four Poaceae species; these genes were involved in the 'response to heat' and 'heat acclimation'. Notably, most of the conserved DEGs belonged to the heat-shock protein (HSP) superfamily. Similar results were also obtained from co-expression analysis. Interestingly, hub-genes of co-expression analysis were found to overlap with conserved genes, especially heat-shock protein (HSP). In C. songorica, 84 HSP and 32 heat-shock transcription factor (HSF) genes were identified in the allotetraploid C. songorica genome, and might have undergone purifying selection during evolutionary history based on syntenic and phylogenetic analysis. By analysing the expression patterns of the CsHSPs and CsHSFs, we found that the transcript abundance of 72.7% of the CsHSP genes and of 62.5% of the CsHSF genes changed under heat stress in both the shoots and roots. Finally, a core regulatory network of HS was constructed on the basis of the CsHSP, CsHSF and other responsive genes in C. songorica.

CONCLUSIONS

Regulatory network and key genes were comprehensively analysed and identified in C. songorica under HS. This study improves our knowledge of thermotolerance mechanisms in native grasses, and also provides candidate genes for potential applications in the genetic improvement of grasses.

Genome-wide identification, characterization, and expression analysis related to autotoxicity of the GST gene family in Cucumis melo L

Glutathione S-transferase (GST) plays an important role in plant resistance to biotic and abiotic stresses. In this paper, the characteristics of melon GST gene family members were analyzed from a genome-wide perspective. Forty-nine GSTs were identified in melon genome, belonging to eight classes. Through the phylogenetic analysis of GST proteins in melon and other plants, it was found that members from the same subfamily in different species clustered together, indicating that the subfamilies of GST have diversified before the divergence within these species. The results of chromosome mapping showed that GSTs were present in all chromosomes except for chromosome 5. Gene replication events played an important role in the expansion and evolution of melon GST gene family. Ten GSTs with significant differential expression were screened in the transcriptome database related to melon autotoxicity stress. The differential expression of these 10 GSTs was detected in roots and leaves of melon seedlings treated with cinnamic acid. The relative expression level of CmGSTU7, CmGSTU10, CmGSTU18, CmGSTF2 and CmGSTL1 in roots of melon seedlings was significantly higher than that in control group. It suggested that the five GSTs might play an important role in cinnamic acid mediated autotoxicity stress in melon. The results of this paper were helpful to reveal the evolution and functional succession of GST family and further understand the response of GST to autotoxicity stress in melon.

Genome-wide identification and comparative analysis of GST gene family in apple (Malus domestica) and their expressions under ALA treatment

Anthocyanins, a subclass of flavonoids, are synthesized at the cytoplasmic surface of the endoplasmic reticulum (ER), which then accumulate in vacuoles. Plant glutathione S-transferase (GST) genes are involved in anthocyanin transportation. Here, a total of 52, 42, 50, and 29 GST genes were identified from apple, pear, peach, and strawberry, respectively, through a comprehensive genome-wide survey. Based on phylogenetic analyses, the GST proteins of the four crops could be divided into the classes Phi, Tau, DHAR, TCHQD, and Lambda. The structure and chromosomal distribution of apple GST genes were further analyzed. The GST gene family expansion in apple likely occurred through tandem duplications, and purifying selection played a pivotal role in the evolution of GST genes. Synteny analysis showed strong microsynteny between apple and Arabidopsis/strawberry, but no microsynteny was detected between apple/strawberry/Arabidopsis and rice. Aminolevulinic acid (ALA), a key precursor of tetrapyrrole compounds, can significantly improve anthocyanin accumulation in fruits, Using RNA-seq and qRT-PCR analysis, we found that ALA treatment led to the differential expression of GST genes in apples. MdGSTF12 was strongly induced by ALA, suggesting that MdGSTF12 may play a role in ALA-induced anthocyanin accumulation. These results provide a detailed overview of GST genes in four Rosaceae species and indicate that GSTs are involved in ALA-induced anthocyanin accumulation.