Genome-wide identification and expression analysis of the fatty acid desaturase genes in Medicago truncatula

Genome-Wide Characterization of Tomato FAD Gene Family and Expression Analysis under Abiotic Stresses

The fatty acid desaturase (FAD) gene family plays a crucial regulatory role in the resistance process of plant biomembranes. To understand the role of FADs in tomato growth and development, this study identified and analyzed the tomato FAD gene family based on bioinformatics analysis methods. In this study, 26 SlFADs were unevenly distributed on 10 chromosomes. Phylogenetic analysis showed that the SlFAD gene family was divided into six branches, and the exon-intron composition and conserved motifs of SlFADs clustered in the same branch were quite conservative. Several hormone and stress response elements in the SlFAD promoter suggest that the expression of SlFAD members is subject to complex regulation; the construction of a tomato FAD protein interaction network found that SlFAD proteins have apparent synergistic effects with SPA and GPAT proteins. qRT-PCR verification results show that SlFAD participates in the expression of tomato root, stem, and leaf tissues; SlFAD8 is mainly highly expressed in leaves; SlFAD9 plays a vital role in response to salt stress; and SlFAB5 regulates all stages of fruit development under the action of exogenous hormones. In summary, this study provides a basis for a systematic understanding of the SlFAD gene family. It provides a theoretical basis for in-depth research on the functional characteristics of tomato SlFAD genes.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

Comparative transcriptomic and evolutionary analysis of FAD-like genes of Brassica species revealed their role in fatty acid biosynthesis and stress tolerance

BACKGROUND

Fatty acid desaturases (FADs) are involved in regulating plant fatty acid composition by adding double bonds to growing hydrocarbon chain. Apart from regulating fatty acid composition FADs are of great importance, and are involved in stress responsiveness, plant development, and defense mechanisms. FADs have been extensively studied in crop plants, and are broadly classed into soluble and non-soluble fatty acids. However, FADs have not yet been characterized in Brassica carinata and its progenitors.

RESULTS

Here we have performed comparative genome-wide identification of FADs and have identified 131 soluble and 28 non-soluble FADs in allotetraploid B. carinata and its diploid parents. Most soluble FAD proteins are predicted to be resided in endomembrane system, whereas FAB proteins were found to be localized in chloroplast. Phylogenetic analysis classed the soluble and non-soluble FAD proteins into seven and four clusters, respectively. Positive type of selection seemed to be dominant in both FADs suggesting the impact of evolution on these gene families. Upstream regions of both FADs were enriched in stress related cis-regulatory elements and among them ABRE type of elements were in abundance. Comparative transcriptomic data analysis output highlighted that FADs expression reduced gradually in mature seed and embryonic tissues. Moreover, under heat stress during seed and embryo development seven genes remained up-regulated regardless of external stress. Three FADs were only induced under elevated temperature whereas five genes were upregulated under Xanthomonas campestris stress suggesting their involvement in abiotic and biotic stress response.

CONCLUSIONS

The current study provides insights into the evolution of FADs and their role in B. carinata under stress conditions. Moreover, the functional characterization of stress-related genes would exploit their utilization in future breeding programs of B. carinata and its progenitors.

Phytohormones regulate the non-redundant response of ω-3 fatty acid desaturases to low temperatures in Chorispora bungeana

To explore the contributions of ω-3 fatty acid desaturases (FADs) to cold stress response in a special cryophyte, Chorispora bungeana, two plastidial ω-3 desaturase genes (CbFAD7, CbFAD8) were cloned and verified in an Arabidopsis fad7fad8 mutant, before being compared with the microsomal ω-3 desaturase gene (CbFAD3). Though these genes were expressed in all tested tissues of C. bungeana, CbFAD7 and CbFAD8 have the highest expression in leaves, while CbFAD3 was mostly expressed in suspension-cultured cells. Low temperatures resulted in significant increases in trienoic fatty acids (TAs), corresponding to the cooperation of CbFAD3 and CbFAD8 in cultured cells, and the coordination of CbFAD7 and CbFAD8 in leaves. Furthermore, the cold induction of CbFAD8 in the two systems were increased with decreasing temperature and independently contributed to TAs accumulation at subfreezing temperature. A series of experiments revealed that jasmonie acid and brassinosteroids participated in the cold-responsive expression of ω-3 CbFAD genes in both C. bungeana cells and leaves, while the phytohormone regulation in leaves was complex with the participation of abscisic acid and gibberellin. These results point to the hormone-regulated non-redundant contributions of ω-3 CbFADs to maintain appropriate level of TAs under low temperatures, which help C. bungeana survive in cold environments.

Desaturases: Structural and mechanistic insights into the biosynthesis of unsaturated fatty acids

This review highlights the key role of fatty acid desaturases in the synthesis of naturally occurring, more common and not unsaturated fatty acids. The three major classes of fatty acid desaturases, such as acyl-lipid, acyl-acyl carrier protein and acyl-coenzyme A, are described in detail, with particular attention to the cellular localisation, the structure, the substrate and product specificity and the expression and regulation of desaturase genes. The review also gives an insight into the biocatalytic reaction of fatty acid desaturation by covering the general and more class-specific mechanistic studies around the synthesis of unsaturated fatty acids Finally, we conclude the review by looking at the numerous novel applications for desaturases in order to meet the very high demand for polyunsaturated fatty acids, taking into account the opportunity for the development of new, more efficient, easily reproducible, sustainable bioengineering advances in the field.

Genome-Wide Identification of Membrane-Bound Fatty Acid Desaturase Genes in Three Peanut Species and Their Expression in Arachis hypogaea during Drought Stress

As a crop irrigated primarily by rain, the quality and yield of peanuts are significantly limited by drought. To date, many studies have indicated that fatty acid desaturase (FAD) genes enhance plant tolerance to drought stresses. In this study, 16, 15, and 31 FADs were identified in Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. All the FADs were divided into four subfamilies, which had relatively conserved gene structures, motifs, and domains. The synteny relationships and chromosomal position analysis showed that the FADs in subgenome pairs, A. duranensis-A. hypogaea (AA) and A. ipaensis-A. hypogaea (BB), were homologous, and their physical locations were consistent. The Ka/Ks results indicated that nine FAD genes underwent a purifying selection, and Ah|FAD3.2 experienced positive selection during tetraploid peanut speciation. Various cis-acting elements related to hormone signaling and stress responsiveness in promoters and the predicted miRNA targeting Ah|FADs suggested that these genes play crucial roles in drought tolerance. The expression profiles of Ah|FADs in 22 tissues and drought-tolerant and -sensitive cultivars under drought stress suggested that 4 and 6 FADs were putative genes related to oil accumulation and drought, respectively. These findings will help provide insight into the potential functional roles of the FAD genes, which may aid in dealing with plant drought stress.

Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar

Fatty acid desaturases (FADs) modulate carbon–carbon single bonds to form carbon–carbon double bonds in acyl chains, leading to unsaturated fatty acids (UFAs) that have vital roles in plant growth and development and their response to environmental stresses. In this study, a total of 23 Populus trichocarpaFAD (PtFAD) candidates were identified from the poplar genome and clustered into seven clades, including FAB2, FAD2, FAD3/7/8, FAD5, FAD6, DSD, and SLD. The exon–intron compositions and conserved motifs of the PtFADs, clustered into the same clade, were considerably conserved. It was found that segmental duplication events are predominantly attributable to the PtFAD gene family expansion. Several hormone- and stress-responsive elements in the PtFAD promoters implied that the expression of the PtFAD members was complicatedly regulated. A gene expression pattern analysis revealed that some PtFAD mRNA levels were significantly induced by abiotic stress. An interaction proteins and gene ontology (GO) analysis indicated that the PtFADs are closely associated with the UFAs biosynthesis. In addition, the UFA contents in poplars were significantly changed under drought and salt stresses, especially the ratio of linoleic and linolenic acids. The integration of the PtFAD expression patterns and UFA contents showed that the abiotic stress-induced PtFAD3/7/8 members mediating the conversion of linoleic and linolenic acids play vital roles in response to osmotic stress. This study highlights the profiles and functions of the PtFADs and identifies some valuable genes for forest improvements.

Tolerant mechanism of model legume plant Medicago truncatula to drought, salt, and cold stresses

Legume plants produce one-third of the total yield of primary crops and are important food sources for both humans and animals worldwide. Frequent exposure to abiotic stresses, such as drought, salt, and cold, greatly limits the production of legume crops. Several morphological, physiological, and molecular studies have been conducted to characterize the response and adaptation mechanism to abiotic stresses. The tolerant mechanisms of the model legume plant Medicago truncatula to abiotic stresses have been extensively studied. Although many potential genes and integrated networks underlying the M. truncatula in responding to abiotic stresses have been identified and described, a comprehensive summary of the tolerant mechanism is lacking. In this review, we provide a comprehensive summary of the adaptive mechanism by which M. truncatula responds to drought, salt, and cold stress. We also discuss future research that need to be explored to improve the abiotic tolerance of legume plants.

Genome-wide identification of FAD gene family and their contributions to the temperature stresses and mutualistic and parasitic fungi colonization responses in banana

Fatty acid desaturase (FAD) plays important roles in plant growth and development and plant defense processes. In this study, we identified 27 MaFAD genes from the banana genome. According to the amino acid sequence similarities, their encoded proteins could be classified into five subfamilies. This classification is consistently supported by their gene and protein structures, conserved motifs and subcellular localizations. Segmental duplication events were found to play predominant roles in the MaFAD gene family expansion. Thirty miRNAs targeting MaFADs were identified and many hormone- and stress-responsive cis-acting elements and transcription factor binding sites (TFBSs) were identified in their promoters, indicating that the MaFADs expression regulation was very complicated. Gene expression analysis showed that some MaFADs showed significant differential expression in response to high and low temperature. FocTR4 influenced greatly the expression of several MaFADs and greatly induced the fatty acid (FA) accumulations in roots. Although S. indica showed no significant influence on the expression of most MaFADs, it could greatly alleviate the influence of FocTR4 on several MaFADs and FA biosynthesis. Our study revealed that MaFADs contributed greatly to the responses of high and low temperature stresses and mutualistic and parasitic fungi colonization in banana.

Genome-Wide Analysis of TIR-NBS-LRR Gene Family in Potato Identified StTNLC7G2 Inducing Reactive Oxygen Species in Presence of Alternaria solani

Resistance gene analogs (RGAs) comprising NBS-LRR gene family members are considered prominent candidates in the development of disease-resistant genotypes. NBS-LRR gene family comprised a very large number of genes; therefore, members of one subfamily TIR-NBS-LRR (TNL) are identified in the present study from Solanum tuberosum genome, followed by their bioinformatics characterization. The study identified a total of 44 genes encoding 60 TNL transcripts with two prominent clusters at chromosome 1 and chromosome 11. Expression analysis of 14 TNL genes after Alternaria solani infection at 1, 2, 3, 5, and 7 days post inoculation in two disease-tolerant varieties, Kufri Jyoti and Kufri Pukhraj, and one relatively susceptible variety, Kufri Chandramukhi, showed differential expression of many genes including a high expression (>15-fold) of StTNLC6G2T1 and StTNLC11G9T1. Functional characterization of one such gene, StTNLC7G2, reveals involvement in the generation of reactive oxygen species under A. solani attack, implicating its putative role in plant defense via hypersensitive response.

Genome-wide identification and characterization of the heavy metal ATPase (HMA) gene family in Medicago truncatula under copper stress

In nature, the normal growth, development, and quality of plants are significantly affected by many abiotic stresses, such as drought, salinity, low temperature, and heavy metals. Among heavy metals, copper is an essential element for plant growth and development but also has a toxic effect on plants when its concentration is excessive. Therefore, plants have evolved a complex regulatory network to regulate the balance of copper ions in cells. Heavy metal ATPases (HMAs), which transport heavy metals to intracellular compartments or detoxify heavy metals present at excessive concentrations, have been extensively studied in model plant species. However, no comprehensive and systematic surveys of members of the HMA gene family have been conducted in the model legume species Medicago truncatula. Here, nine putative MtHMAs were identified in the M. truncatula genome. These MtHMAs were phylogenetically divided into two distinct groups. The members in each group had a relatively conserved gene structure and motif composition. The number of introns in the MtHMAs varied from 5 to 16, with the majority of these genes containing 8 introns. The expression patterns showed that MtHMAs exhibit preferential or distinct expression patterns among different tissues. Finally, the expression patterns of the members of this gene family were verified in the leaves and roots of plants under Cu stress. Our findings will be valuable for the functional investigation and application of members of this gene family in M. truncatula and other related legume species.

Integrated Analysis of Fatty Acid Metabolism and Transcriptome Involved in Olive Fruit Development to Improve Oil Composition

Olea europaea L. is an important oil crop with excellent nutritional properties. In this study, a full-length transcriptome combined with fatty acid composition was used to investigate the molecular mechanism of fatty acid (FA) metabolism of olive fruits at various stages of development (S1–S5). A total of 34 fatty acids (FAs) were measured using gas chromatography-mass spectrometry (GC-MS). All transcripts of FA metabolism during olive fruit development were studied, including glycolysis, fatty acid synthesis, triacylglycerol synthesis, and FA degradation. A total of 100 transcripts of 11 gene families, 68 transcripts of 12 gene families, 55 transcripts of 7 gene families, and 28 transcripts of 7 gene families were identified as encoding for enzymes involved in FA metabolism. Furthermore, one of the critical reactions in TAG metabolism is the activation of fatty acyl chains to fatty acyl CoA, which is catalyzed by long-chain acyl CoA synthetases (LACS). Phylogenetic analysis showed that 13 putative LACS-encoding genes clustered into five groups, of which two putative transcripts encoding LACS6/7 may participate in FA degradation. The aim of this study was to evaluate the fatty acid from synthesis to degradation pathways during olive fruit development to provide a better understanding of the molecular mechanism of FA metabolism during olive fruit maturation and provide information to improve the synthesis of oil components that are beneficial to human health.

Genome-wide investigation and expression analysis of membrane-bound fatty acid desaturase genes under different biotic and abiotic stresses in sunflower (Helianthus annuus L.)

Membrane-bound fatty acid desaturase (FAD) gene family plays crucial roles in regulation of fatty acid (FA) compositions in plants. Sunflower (Helianthus annuus L.) is an important oilseed crop in the world; however, no comprehensive study on exploring the role of FAD family in relation to stress tolerance in sunflower has been performed yet. In this study, we identified 40 putative FAD genes in H. annuus (HaFAD), which were unevenly distributed across 13 of the total 17 chromosomes. Phylogenetic analysis indicated that HaFAD genes were divided into four subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that tandem duplication events played a crucial role in the expansion of HaFAD gene family. In addition, tissue-specific expression showed that 32 HaFAD genes were widely expressed in various tissues or organs of sunflower. Furthermore, qRT-PCR results revealed significant expression changes of HaFAD genes in response to abiotic (cadmium, drought) and biotic (Orobanche cumana) stresses, suggesting their important functions in response to different stresses. Therefore, our results provide insights into HaFAD gene family in response to different stresses, and some specific up-regulated genes such as HaFAD3.2, HaADS8, HaFAD2.1, and HaADS9 would be the potential candidate genes for the sunflower tolerance breeding.

Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.)

Backgrounds

Fatty acid desaturases (FADs) introduce a double bond into the fatty acids acyl chain resulting in unsaturated fatty acids that have essential roles in plant development and response to biotic and abiotic stresses. Wheat germ oil, one of the important by-products of wheat, can be a good alternative for edible oils with clinical advantages due to the high amount of unsaturated fatty acids. Therefore, we performed a genome-wide analysis of the wheat FAD gene family (TaFADs).

Results

68 FAD genes were identified from the wheat genome. Based on the phylogenetic analysis, wheat FADs clustered into five subfamilies, including FAB2, FAD2/FAD6, FAD4, DES/SLD, and FAD3/FAD7/FAD8. The TaFADs were distributed on chromosomes 2A-7B with 0 to 10 introns. The Ka/Ks ratio was less than one for most of the duplicated pair genes revealed that the function of the genes had been maintained during the evolution. Several cis-acting elements related to hormones and stresses in the TaFADs promoters indicated the role of these genes in plant development and responses to environmental stresses. Likewise, 72 SSRs and 91 miRNAs in 36 and 47 TaFADs have been identified. According to RNA-seq data analysis, the highest expression in all developmental stages and tissues was related to TaFAB2.5, TaFAB2.12, TaFAB2.15, TaFAB2.17, TaFAB2.20, TaFAD2.1, TaFAD2.6, and TaFAD2.8 genes while the highest expression in response to temperature stress was related to TaFAD2.6, TaFAD2.8, TaFAB2.15, TaFAB2.17, and TaFAB2.20. Furthermore, docking simulations revealed several residues in the active site of TaFAD2.6 and TaFAD2.8 in close contact with the docked oleic acid that could be useful in future site-directed mutagenesis studies to increase the catalytic efficiency of them and subsequently improve agronomic quality and tolerance of wheat against environmental stresses.

Conclusions

This study provides comprehensive information that can lead to the detection of candidate genes for wheat genetic modification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07199-1.

Transformation and Characterization of Δ12-Fatty Acid Acetylenase and Δ12-Oleate Desaturase Potentially Involved in the Polyacetylene Biosynthetic Pathway from Bidens pilosa

Bidens pilosa is commonly used as an herbal tea component or traditional medicine for treating several diseases, including diabetes. Polyacetylenes have two or more carbon–carbon triple bonds or alkynyl functional groups and are mainly derived from fatty acid and polyketide precursors. Here, we report the cloning of full-length cDNAs that encode Δ12-fatty acid acetylenase (designated BPFAA) and Δ12-oleate desaturase (designated BPOD) from B. pilosa, which we predicted to play a role in the polyacetylene biosynthetic pathway. Subsequently, expression vectors carrying BPFAA or BPOD were constructed and transformed into B. pilosa via the Agrobacterium-mediated method. Genomic PCR analysis confirmed the presence of transgenes and selection marker genes in the obtained transgenic lines. The copy numbers of transgenes in transgenic lines were determined by Southern blot analysis. Furthermore, 4–5 FAA genes and 2–3 OD genes were detected in wild-type (WT) plants. Quantitative real time-PCR revealed that some transgenic lines had higher expression levels than WT. Western blot analysis revealed OD protein expression in the selected transformants. High-performance liquid chromatography profiling was used to analyze the seven index polyacetylenic compounds, and fluctuation patterns were found.

The expression of several pepper fatty acid desaturase genes is robustly activated in an incompatible pepper-tobamovirus interaction, but only weakly in a compatible interaction

The replication of positive strand RNA viruses in plant cells is markedly influenced by the desaturation status of fatty acid chains in lipids of intracellular plant membranes. At present, little is known about the role of lipid desaturation in the replication of tobamoviruses. Therefore, we investigated the expression of fatty acid desaturase (FAD) genes and the fatty acid composition of pepper leaves inoculated with two different tobamoviruses. Obuda pepper virus (ObPV) inoculation induced a hypersensitive reaction (incompatible interaction) while Pepper mild mottle virus (PMMoV) inoculation caused a systemic infection (compatible interaction). Changes in the expression of 16 FADs were monitored in pepper leaves following ObPV and PMMoV inoculations. ObPV inoculation rapidly and markedly upregulated seven Δ12-FADs that encode enzymes putatively located in the endoplasmic reticulum membrane. In contrast, PMMoV inoculation resulted in a weaker but rapid upregulation of two Δ12-FADs and a Δ15-FAD. The expression of genes encoding plastidial FADs was not influenced neither by ObPV nor by PMMoV. In accordance with gene expression results, a significant accumulation of linoleic acid was observed by gas chromatography-mass spectrometry in ObPV-, but not in PMMoV-inoculated leaves. ObPV inoculation led to a marked accumulation of H₂O₂ in the inoculated leaves. Therefore, the effect of H₂O₂ treatments on the expression of six tobamovirus-inducible FADs was also studied. The expression of these FADs was upregulated to different degrees by H₂O₂ that correlated with ObPV-inducibility of these FADs. These results underline the importance of further studies on the role of pepper FADs in pepper-tobamovirus interactions.

Molecular cloning and functional analysis of a Δ12-fatty acid desaturase from the Antarctic microalga Chlamydomonas sp. ICE-L

Chlamydomonas sp. ICE-L, which can thrive in extreme environments of the Antarctic, could represent a promising alternative for polyunsaturated fatty acid (PUFA) production. A new Δ¹²-fatty acid desaturase (FAD)-encoding gene (Δ ¹² CiFAD), 1269 bp in size, was cloned from Chlamydomonas sp. ICE-L. Bioinformatics analysis showed that Δ ¹² CiFAD-encoded protein was homologous to known FADs with conserved histidine motifs, and localized to the chloroplast. Functional analysis of Δ ¹² CiFAD indicated that recombinant Synechococcus 6803 expressing Δ¹²CiFAD could accumulate C18:2, whereas recombinant Saccharomyces cerevisiae expressing this enzyme could not accumulate C18:2 or any other new fatty acids. These results indicate that Δ¹²CiFAD is a functional enzyme in the chloroplast that can adjust Chlamydomonas sp. ICE-L cell membrane fluidity to adapt to Antarctic extreme low-temperature environments, which give us insights into the frigostable and cold-resistant mechanisms of hypothermic organisms.

Hydrogen Peroxide and Nitric Oxide Crosstalk Mediates Brassinosteroids Induced Cold Stress Tolerance in Medicago truncatula

Brassinosteroids (BRs) play pivotal roles in modulating plant growth, development, and stress responses. In this study, a Medicago truncatula plant pretreated with brassinolide (BL, the most active BR), enhanced cold stress tolerance by regulating the expression of several cold-related genes and antioxidant enzymes activities. Previous studies reported that hydrogen peroxide (H₂O₂) and nitric oxide (NO) are involved during environmental stress conditions. However, how these two signaling molecules interact with each other in BRs-induced abiotic stress tolerance remain largely unclear. BL-pretreatment induced, while brassinazole (BRZ, a specific inhibitor of BRs biosynthesis) reduced H₂O₂ and NO production. Further, application of dimethylthiourea (DMTU, a H₂O₂ and OH- scavenger) blocked BRs-induced NO production, but BRs-induced H₂O₂ generation was not sensitive to 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl-3-oxide (PTIO, a scavenger of NO). Moreover, pretreatment with DMTU and PTIO decreased BL-induced mitochondrial alternative oxidase (AOX) and the photosystem capacity. However, pretreatment with PTIO was found to be more effective than DMTU in reducing BRs-induced increases in Valt, Vt, and MtAOX1 gene expression. Similarly, BRs-induced photosystem II efficiency was found in NO dependent manner than H₂O₂. Finally, we conclude that H₂O₂ was involved in NO generation, whereas NO was found to be crucial in BRs-induced AOX capacity, which further contributed to the protection of the photosystem under cold stress conditions in Medicago truncatula.