Genome wide and evolutionary analysis of heat shock protein 70 proteins in tomato and their role in response to heat and drought stress.

Heat shock protein 70 (HSP70) proteins play a crucial role in mitigating the detrimental effects of abiotic stresses in plants. In the present study, 21 full length non-redundant SlHSP70 genes were detected and characterized in tomato (Solanum lycopersicum L.). The SlHSP70 genes were classified into four groups based on phylogenetic analysis. Similarities were observed in gene features and motif structures of SlHSP70s belonging to the same group. SlHSP70 genes were unevenly and unequally mapped on 11 chromosomes. Segmental and tandem duplication are the main events that have contributed to the expansion of the SlHSP70 genes. A large number of groups and sub-groups were generated during comparative analysis of HSP70 genes in multiple plant species including tomato. These findings indicated a common ancestor which created diverse sub-groups prior to a mono-dicot split. The selection pressure on specific codons was identified through a maximum-likelihood approach and we found some important coding sites in the coding region of all groups. Diversifying positive selection was indirectly associated with evolutionary changes in SlHSP70 proteins and suggests that gene evolution modulated the tomato domestication event. In addition, expression analysis using RNA-seq revealed that 21 SlHSP70 genes were differentially expressed in response to drought and heat stress. SlHSP70-5 was down-regulated by heat treatment and up-regulated by drought stress. Furthermore, the expression of some of the duplicate genes was partially redundant, while others showed functional diversity. Our results indicate the diverse role of HSP70 gene family in S. lycopersicum under drought and heat stress conditions and open the gate for further investigation of HSP70 gene family functions, especially under drought and heat stress.

GmMAX2-D14 and -KAI interaction-mediated SL and KAR signaling play essential roles in soybean root nodulation.

Despite of important functions of strigolactones (SLs) and karrikins (KARs) in plant development, plant-parasite and plant-fungi interactions, their roles in soybean-rhizobia interaction remain elusive. SL/KAR signaling genes GmMAX2a, GmD14s, and GmKAIs are activated by rhizobia infection. GmMAX2a restored atmax2 root hair defects and soybean root hairs were changed in GmMAX2a overexpression (GmMAX2a-OE) or knockdown (GmMAX2a-KD) mutants. GmMAX2a-KD gave fewer, whereas GmMAX2a-OE produced more nodules than GUS hairy roots. Mutation of GmMAX2a in its KD or OE transgenic hairy roots affected the rhizobia infection-induced increases in early nodulation gene expression. Both mutant hairy roots also displayed the altered auxin, jasmonate and abscisic acid levels, as further verified by transcriptomic analyses of their synthetic genes. Overexpression of an auxin synthetic gene GmYUC2a also affected SL and KAR signaling genes. GmMAX2a physically interacted with SL/KAR receptors GmD14s, GmKAIs, and GmD14Ls with different binding affinities, depending on variations in the critical amino acids, forming active D14/KAI-SCFMAX2 complexes. The knockdown mutant roots of the nodule-specifically expressing GmKAIs and GmD14Ls gave fewer nodules, with altered expression of several early nodulation genes. The expression levels of GmKAIs, and GmD14Ls were markedly changed in GmMAX2a mutant roots, so did their target repressor genes GmD53s and GmSMAX1s. Thus, SL and KAR signaling were involved in soybean-rhizobia interaction and nodulation partly through interactions with hormones, and this may explain the different effects of MXA2 orthologs on legume determinate and indeterminate nodulation. The study provides fresh insights into the roles of GmMAX2-mediated SL/KAR signaling in soybean root hair and nodule formation.

Isoflavone malonyl-CoA acyltransferase GmMaT2 is involved in nodulation of soybean by modifying synthesis and secretion of isoflavones.

Malonyl-CoA:flavonoid acyltransferases (MaTs) modify isoflavones, but only a few have been characterized for activity and assigned to specific physiological processes. Legume roots exude isoflavone malonates into the rhizosphere, where they are hydrolyzed into isoflavone aglycones. Soybean GmMaT2 was highly expressed in seeds, root hairs, and nodules. GmMaT2 and GmMaT4 recombinant enzymes used isoflavone 7-O-glucosides as acceptors and malonyl-CoA as an acyl donor to generate isoflavone glucoside malonates. GmMaT2 had higher activity towards isoflavone glucosides than GmMaT4. Overexpression in hairy roots of GmMaT2 and GmMaT4 produced more malonyldaidzin, malonylgenistin, and malonylglycitin, and resulted in more nodules than control. However, only GmMaT2 knockdown (KD) hairy roots showed reduced levels of malonyldaidzin, malonylgenistin, and malonylglycitin, and, likewise, reduced nodule numbers. These were consistent with the up-regulation of only GmMaT2 by rhizobial infection, and higher expression levels of early nodulation genes in GmMaT2- and GmMaT4-overexpressing roots, but lower only in GmMaT2-KD roots compared with control roots. Higher malonyl isoflavonoid levels in transgenic hairy roots were associated with higher levels of isoflavones in root exudates and more nodules, and vice versa. We suggest that GmMaT2 participates in soybean nodulation by catalyzing isoflavone malonylation and affecting malonyl isoflavone secretion for activation of Nod factor and nodulation.

Comparative analysis of the mitochondrial proteins reveals complex structural and functional relationships in Fasciola species.

Mitochondria is a cellular source of energy, appears to play an essential role in dealing with cellular stress induced by environmental stimuli. The genetic diversity of mitochondrial genes involved in oxidative phosphorylation affecting the production of cellular energy and regional adaptation to various ecological (climatic) pressures affecting amino acid sequences (variants of protein). However, little is known about the combined effect of protein changes on cell-level metabolic alterations in simultaneous exposure to various environmental conditions, including mitochondrial dysfunction and oxidative stress induction. The present study was designed to address this issue by analyzing the mitochondrial proteins in Fasciola species including Cytochrome oxidase (COX1, COX2, COX3, and CYTB) and NADH dehydrogenase (ND1, ND2, ND3, ND4, ND5, and ND6). Mitochondrial proteins were used for detailed computational investigation, using available standard bioinformatics tools to exploit structural and functional relationships. These proteins in Fasciola hepatica, Fasciola gigentica, and Fasciola jacksoni were functionally annotated using public databases. The results showed that the protein of COX1 of F. hepatica, F. gigantica, and F. jacksoni consist of 510, 513, and 517 amino acids, respectively. The alignment of proteins showed that these proteins are conserved in the same regions at ten positions in COX and CYTB proteins while at twelve locations in NADH. Three-dimensional structure of COX, CYTB, and NADH proteins were compared and showed differences in additional conserved and binding sites in COX and CYTB proteins as compared to NADH in three species of Fasciola. These results based on the amino acid diversity pattern were used to identify sites in the enzyme and the variations in mitochondrial proteins among Fasciola species. Our study provides valuable information for future experimental studies, including identification of therapeutic, diagnostic, and immunoprophylactic interests with novel mitochondrial proteins.

Adaptation to host-specific bacterial pathogens drive rapid evolution of novel PhoP/PhoQ regulation pathway modulating the virulence.

The presence of the PhoP-PhoQ system is usually different in various bacterial groups, suggesting that PhoP can control the expression of different genes in species. However, little is known about the evolution of the PhoP-PhoQ system among bacterial pathogens. Here, we study the evolution of PhoP and PhoQ regulation in 15 species of Enterobacteriaceae family. We have determined that the regulatory objectives adopted by PhoP and PhoQ are mainly different, due to the result of horizontal gene transfer events and even the change in the genetic content between closely related species. We have compared many possibilities tests (M1 vs. M2 and M7 with M8) to determine the positive selection. Estimating parameters at M1 and M2, with positive selection in M2 of the two proteins. The proportions of positive selection sites significant with ω = 4.53076 for PhoP and ω = 4.21041 PhQ. M8 was significant for PhoP and PhQ proteins. To further confirm the positive selection results, we used the Selecton server to confer positive selection on individual sites using the Mechanistic-Empirical Combination model, and we noticed that several sites had been identified under selection pressure during the evolution. There was a strong indication for the positive selection in bacterial genes of PhoP and PhoQ showed the results. By the use of REL and IFEL, the positive selection for PhoP was detected 14 and 11 sites respectively at different codon positions. The positively selected sites of amino acids such as Arginine, Alanine, Lysine, and Leucine are more important for the production of signals. Our results suggest that the positive selection of PhoP-PhoQ genes in host adaptation during evolution raises an intriguing possibility causes subtle variations in actions of PhoP-PhoQ and also increases the opportunities that cause modification in protein structure for the evolution of increasing pathogenicity in bacterial pathogens.

A genome-wide approach to the comprehensive analysis of GASA gene family in Glycine max.

A vital role of short amino acid gene family, gibberellic acid stimulated arabidopsis (GASA), has been reported in plant growth and development. Although, little information is available about these cysteine rich short proteins in different plant species and this is the first comprehensive approach to exploit available genomic data and to analyze the GASA family in G. max. The phylogenetic and sequence composition analysis distributed the 37 identified GmGASA genes into three groups. Further investigation of the tissue expression pattern, phylogenetic analysis, motif, gene structure, chromosome distributions, duplication patterns, positive-selection pressure and cis-element analysis of 37 GmGASA genes. A conserved GASA domain was found in all identified GmGASA genes and exhibited similar characteristics. The online gene expression profile based analysis of GmGASA genes reveled that these genes were highly expressed in almost all soybean parts and some have high expression in flower which indicates that GmGASA genes displayed special or distinct expression pattern among different tissues. The segmental duplication was found in five pairs from 37 GmGASA genes and was distributed on 15 different chromosomes. The Ka/Ks ratio of 5 pairs of segmentally duplicated gene indicated that after the occurrence of duplication events, the duplicated gene pairs were purified and selected after restrictive functional differentiation. This investigated study of GmGASA gene will useful to support the statement about GASA genes role during flower induction in flowering plants.

Strigolactones promote rhizobia interaction and increase nodulation in soybean (Glycine max).

Strigolactones (SLs) play an important role in controlling root growth, shoot branching, and plant-symbionts interaction. Despite the importance, the components of SL biosynthesis and signaling have not been unequivocally explored in soybean. Here we identified the putative components of SL synthesis enzymes GmMAX1a and GmMAX4a with tissue expression patterns and were apparently regulated by rhizobia infection and changed during nodule development. GmMAX1a and GmMAX4a were further characterized in soybean nodulation with knockdown transgenic hairy roots. GmMAX1a and GmMAX4a knockdown lines exhibit decreased nodule number and expression levels of several nodulation genes required for nodule development. Hormone analysis showed that GmMAX1a and GmMAX4a knockdown hairy roots had increased physiological level of ABA and JA but significantly decreased auxin content. This study not only revealed the conservation of SL biosynthesis but also showed close interactions between SL and other hormone signaling in controlling plant development and legume-rhizobia interaction.

Functional characterization of soybean strigolactone biosynthesis and signaling genes in Arabidopsis MAX mutants and GmMAX3 in soybean nodulation.

BACKGROUND: Strigolactones (SLs) play important roles in controlling root growth, shoot branching, and plant-symbionts interaction. Despite the importance, the components of SL biosynthesis and signaling have not been unequivocally explored in soybean. RESULTS: Here we identified the putative components of SL synthetic enzymes and signaling proteins in soybean genome. Soybean genome contains conserved MORE AXILLARY BRANCHING (MAX) orthologs, GmMAX1s, GmMAX2s, GmMAX3s, and GmMAX4s. The tissue expression patterns are coincident with SL synthesis in roots and signaling in other tissues under normal conditions. GmMAX1a, GmMAX2a, GmMAX3b, and GmMAX4a expression in their Arabidopsis orthologs' mutants not only restored most characteristic phenotypes, such as shoot branching and shoot height, leaf shape, primary root length, and root hair growth, but also restored the significantly changed hormone contents, such as reduced JA and ABA contents in all mutant leaves, but increased auxin levels in atmax1, atmax3 and atmax4 mutants. Overexpression of these GmMAXs also altered the hormone contents in wild-type Arabidopsis. GmMAX3b was further characterized in soybean nodulation with overexpression and knockdown transgenic hairy roots. GmMAX3b overexpression (GmMAX3b-OE) lines exhibited increased nodule number while GmMAX3b knockdown (GmMAX3b-KD) decreased the nodule number in transgenic hairy roots. The expression levels of several key nodulation genes were also altered in GmMAX3b transgenic hairy roots. GmMAX3b overexpression hairy roots had reduced ABA, but increased JA levels, with no significantly changed auxin content, while the contrast changes were observed in GmMAX3b-KD lines. Global gene expression in GmMAX3b-OE or GmMAX3b-KD hairy roots also revealed that altered expression of GmMAX3b in soybean hairy roots changed several subsets of genes involved in hormone biosynthesis and signaling and transcriptional regulation of nodulation processes. CONCLUSIONS: This study not only revealed the conservation of SL biosynthesis and signaling in soybean, but also showed possible interactions between SL and other hormone synthesis and signaling during controlling plant development and soybean nodulation. GmMAX3b-mediated SL biosynthesis and signaling may be involved in soybean nodulation by affecting both root hair formation and its interaction with rhizobia.

Identification, Classification and Characterization Analysis of FBXL Gene in Cotton.

F-box/LR (FBXL), Leucine-rich repeats in F-box proteins, belongs to the Skp1-Cullin1-F-box protein (SCF) E3 ligase family. FBXL genes play important roles in plant growth, such as plant hormones, responses to environmental stress, and floral organ development. Here, a total of 518 FBXL genes were identified and analyzed in six plant species. Phylogenetic analysis showed that AtFBXLs, VvFBXLs, and GrFBXLs were clustered into three subfamilies (Ⅰ-Ⅲ). Based on the composition of the F-box domain and carboxyl-terminal amino acid sequence, FBXL proteins were classified into three types (Type-A/-B/-C). Whole-genome duplication (WGD) along with tandem duplications and segmental contributed to the expansion of this gene family. The result indicates that four cotton species are also divided into three subfamilies. FBXLs in cotton were classified into three clades by phylogenetic and structural analyses. Furthermore, expression analyses indicated that the expression patterns of GhFBXLs in different cotton tissues were different. The highly expressed of GH_A07G2363 in 5-8 mm anthers, indicates that this gene might play a role in the reproductive process, providing candidate genes for future studies on cotton fertility materials. This study provides an original functional opinion and a useful interpretation of the FBXL protein family in cotton.

Genome-wide analysis of sucrose synthase family in soybean and their expression in response to abiotic stress and seed development.

The sucrose synthase (SS) is an important enzyme family which play a vital role in sugar metabolism to improve the fruit quality of the plants. In many plant species, the members of SS family have been investigated but the detailed information is not available in legumes particularly and Glycine max specifically. In the present study, we found thirteen SS members (GmSS1-GmSS13) in G. max genome. High conserved regions were present in the GmSS sequences that may due to the selection pressure during evolutionary events. The segmental duplication was the major factor to increase the number of GmSS family members. The identified thirteen GmSS genes were divided into Class I, Class II and Class III with variable numbers of genes in each class. The protein interaction of GmSS gave the co-expression of sucrose synthase with glucose-1-phosphate adenylyltransferase while SLAC and REL test found number of positive sites in the coding sequences of SS family members. All the GmSS family members except GmSS7 and few of class III members, were highly expressed in all the soybean tissues. The expression of the class I members decreased during seed development, whireas, the class II members expression increased during the seed developing, may involve in sugar metabolism during seed development. Solexa sequencing libraries of acidic condition (pH 4.2) stress samples showed that the expression of class I GmSS genes increased 1- to 2-folds in treated samples than control. The differential expression pattern was observed between the members of a paralogous. This study provides detailed genome-wide analysis of GmSS family in soybean that will provide new insights for future evolutionary and soybean breeding to improve the plant growth and development.

Molecular Characterization of Infectious Bronchitis Virus Strain HH06 Isolated in a Poultry Farm in Northeastern China.

Spike (S) glycoprotein is an important virulent factor for coronaviruses (CoVs), and variants of CoVs have been characterized based on S gene analysis. We present phylogenetic relationship of an isolated infectious bronchitis virus (IBV) strain with reference to the available genome and protein sequences based on network, multiple sequence, selection pressure, and evolutionary fingerprinting analysis in People's Republic of China. One hundred and elven strains of CoVs i.e., Alphacoronaviruses (Alpha-CoVs; n = 12), Betacoronaviruses (Beta-CoVs; n = 37), Gammacoronaviruses (Gamma-CoVs; n = 46), and Deltacoronaviruses (Delta-CoVs; n = 16) were selected for this purpose. Phylogenetically, SARS-CoV-2 and SARS-CoVs clustered together with Bat-CoVs and MERS-CoV of Beta-CoVs (C). The IBV HH06 of Avian-CoVs was closely related to Duck-CoV and partridge S14, LDT3 (teal and chicken host). Beluga whale-CoV (SW1) and Bottlenose dolphin-CoVs of mammalian origin branched distantly from other animal origin viruses, however, making group with Avian-CoVs altogether into Gamma-CoVs. The motif analysis indicated well-conserved domains on S protein, which were similar within the same phylogenetic class and but variable at different domains of different origins. Recombination network tree indicated SARS-CoV-2, SARS-CoV, and Bat-CoVs, although branched differently, shared common clades. The MERS-CoVs of camel and human origin spread branched into a different clade, however, was closely associated closely with SARS-CoV-2, SARS-CoV, and Bat-CoVs. Whereas, HCoV-OC43 has human origin and branched together with bovine CoVs with but significant distant from other CoVs like SARS CoV-2 and SARS-CoV of human origin. These findings explain that CoVs' constant genetic recombination and evolutionary process that might maintain them as a potential veterinary and human epidemic threat.

CsTCPs regulate shoot tip development and catechin biosynthesis in tea plant (Camellia sinensis).

The growth of leaves and biosynthesis of characteristic secondary metabolites are critically important for tea production and quality control. However, little is known about the coordinated regulation of leaf development and catechin biosynthesis in tea plants. Here, we reported that TCP TFs are involved in both catechin biosynthesis and leaf development. An integrated analysis of catechin profiling and CsTCP expression in different tissues of plants under various environmental conditions at different developmental stages indicated significant correlations between the transcript levels of CIN-type TCPs and catechin production. CIN-type CsTCP3 and CsTCP4 and PCF-type CsTCP14 interacted with the MYB-bHLH-WD40 repeat (MBW) complex by forming a CsTCP3-CsTT8 heterodimer and modulating the transactivation activity of the promoters of anthocyanin synthase (CsANS1) and anthocyanidin reductase (CsANR1). Four types of microRNA/target modules, miR319b/CsTCP3-4, miR164b/CsCUC, miR396/CsGRF-GIF, and miR165b/HD-ZIPIII ones, were also identified and characterized for their functions in the regulation of the development of tea plant shoot tips and leaf shape. The results of these modules were reflected by their different expression patterns in developing buds and leaves that had distinctly different morphologies in three different tea plant varieties. Their roles in the regulation of catechin biosynthesis were also further verified by manipulation of microRNA319b (miR319b), which targets the transcripts of CsTCP3 and CsTCP4. Thus, CsTCPs represent at least one of these important groups of TFs that can integrate tea plant leaf development together with secondary metabolite biosynthesis. Our study provides new insight into shoot tip development and catechin production in tea plants and lays a foundation for further mechanistic understanding of the regulation of tea plant leaf development and secondary metabolism.

Genome-Wide Analysis of Serine Carboxypeptidase-Like Acyltransferase Gene Family for Evolution and Characterization of Enzymes Involved in the Biosynthesis of Galloylated Catechins in the Tea Plant (Camellia sinensis).

Tea (Camellia sinensis L.) leaves synthesize and concentrate a vast array of galloylated catechins (e.g., EGCG and ECG) and non-galloylated catechins (e.g., EGC, catechin, and epicatechin), together constituting 8%-24% of the dry leaf mass. Galloylated catechins account for a major portion of soluble catechins in tea leaves (up to 75%) and make a major contribution to the astringency and bitter taste of the green tea, and their pharmacological activity for human health. However, the catechin galloylation mechanism in tea plants is largely unknown at molecular levels. Previous studies indicated that glucosyltransferases and serine carboxypeptidase-like acyltransferases (SCPL) might be involved in the process. However, details about the roles of SCPLs in the biosynthesis of galloylated catechins remain to be elucidated. Here, we performed the genome-wide identification of SCPL genes in the tea plant genome. Several SCPLs were grouped into clade IA, which encompasses previously characterized SCPL-IA enzymes with an acylation function. Twenty-eight tea genes in this clade were differentially expressed in young leaves and vegetative buds. We characterized three SCPL-IA enzymes (CsSCPL11-IA, CsSCPL13-IA, CsSCPL14-IA) with galloylation activity toward epicatechins using recombinant enzymes. Not only the expression levels of these SCPLIA genes coincide with the accumulation of galloylated catechins in tea plants, but their recombinant enzymes also displayed ß-glucogallin:catechin galloyl acyltransferase activity. These findings provide the first insights into the identities of genes encoding glucogallin:catechin galloyl acyltransferases with an active role in the biosynthesis of galloylated catechins in tea plants.

Adaptive Molecular Evolution of AKT3 Gene for Positive Diversifying Selection in Mammals.

The V-Akt Murine Thymoma Viral Oncogene Homolog 3 (AKT3) gene is of the serine/threonine-protein kinase family and influences the production of milk fats and cholesterol by acting on the sterol administrative area restricting protein (SREBP). The AKT3 gene is highly preserved in animals, and during lactation in cattle, its expression increases. The AKT3 gene is expressed in the digestive system, mammary gland, and immune cells. A phylogenetic investigation was performed to clarify the evolutionary role of AKT3, by maximum probability. The AKT3 gene sequence data of various mammalian species was evident even with animals undergoing breeding selection. From 39 mammalian species studied, there was a signal of positive diversifying selection with Hominidae at 13Q, 16G, 23R, 24P, 121P, 294K, 327V, 376L, 397K, 445T, and 471F among other codon sites of the AKT3 gene. These sites were codes for amino acids such as arginine, proline, lysine, and leucine indicating major roles for the function of immunological proteins, and in particular, the study highlighted the importance of changes in gene expression of AKT3 on immunity.

Hydrogen sulfide enhances rice tolerance to nickel through the prevention of chloroplast damage and the improvement of nitrogen metabolism under excessive nickel.

Hydrogen sulfide (H2S) modulates plant tolerance to abiotic stresses, but its regulatory effects on nitrogen metabolism and chloroplast protection under nickel (Ni) stress in crop plants remain elusive. Taking this into account, we investigated the potential roles of sodium hydrosulfide (NaHS), a H2S generator, in the improvement of growth performance of rice plants under Ni stress. Results showed that NaHS successfully reversed the adverse effects of Ni, as reflected in plant growth and biomass, and photosynthesis attributes including photosynthetic rates, stomatal conductance, transpiration rate, internal CO2 concentration and photosynthetic pigment contents. NaHS generated H2S plays a crucial role in controlling the photosynthetic machinery of rice as evidenced by the ultrastructure of chloroplast viewed under transmission electron microscope (TEM). The reduced content of Ni in roots and leaves of NaHS-supplemented Ni-stressed plants has revealed the restricted uptake and accumulation of Ni. A rescue of NaHS to the Ni-induced decline in nitrate (NO3-) content and the activities NO3- biosynthesizing enzymes nitrate reductase, nitrite reductase, glutamate synthase, glutamate oxaloacetate transaminase, glutamine synthetase, and glutamate pyruvate transaminase in leaves indicated a positive role of H2S on NO3- metabolism in rice under Ni stress. NaHS application also reverted Ni-mediated increases in ammonium (NH4+) content and glutamate dehydrogenase activity, implying H2S-induced alleviation of NH4+ toxicity. The regulatory effects of H2S on nitrogen metabolism was further confirmed by increased and decreased transcript abundance of NO3- and NH4+ metabolism associated genes, respectively. Our study suggests a decisive role of H2S in controlling Ni toxicity as elucidated by the novel findings such as enhanced gas exchanged parameters, Ni homeostasis and chloroplast protection. Moreover, this article highlights the significance of H2S in controlling chloroplast biogenesis and nitrogen metabolism in rice crop under Ni stress.

Nitric oxide induces rice tolerance to excessive nickel by regulating nickel uptake, reactive oxygen species detoxification and defense-related gene expression.

Soil contamination with nickel (Ni) is a persistent threat to crop production worldwide. The present study examined the putative roles of nitric oxide (NO) in improving Ni-tolerance in rice. Our findings showed that application of exogenous sodium nitroprusside (SNP), a NO donor, significantly improved the growth performance of rice seedlings when grown under excessive Ni. The enhanced Ni-tolerance of rice prompted by SNP could be ascribed to its ability to regulate Ni uptake, decrease Ni-induced oxidative stress as evidenced by reduced levels of hydrogen peroxide, malondialdehyde, and electrolyte leakage in Ni-stressed plants. The positive roles of NO against Ni-toxicity also reflected through its protective effects on photosynthetic pigments, soluble proteins and proline. SNP also boosted antioxidant capacity in Ni-stressed plants by maintaining increased levels of ascorbate, enhanced activities of ROS-detoxifying enzymes, particularly peroxidase (POD) and catalase (CAT) in both roots and shoots compared with Ni-stressed alone plants. Moreover, SNP treatment also upregulated the transcript levels of CAT, POD, ascorbate peroxidase, glutathione reductase and superoxide dismutase genes in shoots under Ni-stress. Using different sulfide compounds and NO scavenger cPTIO, we also provided evidence that NO, rather than other byproducts of SNP, contributed to the improved performance of rice seedlings under Ni-stress. Collectively, our results conclude that exogenous SNP-mediated modulation of endogenous NO enhanced rice tolerance to Ni-stress by restricting Ni accumulation, maintaining photosynthetic performance and reducing oxidative damage through improved antioxidant system, thereby suggesting NO as an effective stress regulator in mitigating Ni-toxicity in economically important rice, and perhaps in other crop plants.