In silico identification and expression analysis of superoxide dismutase (SOD) gene family in Medicago truncatula

Superoxide dismutase gene family in cassava revealed their involvement in environmental stress via genome-wide analysis

Superoxide dismutase (SOD) is a crucial metal-containing enzyme that plays a vital role in catalyzing the dismutation of superoxide anions, converting them into molecular oxygen and hydrogen peroxide, essential for enhancing plant stress tolerance. We identified 8 SOD genes (4 CSODs, 2 FSODs, and 2 MSODs) in cassava. Bioinformatics analyses provided insights into chromosomal location, phylogenetic relationships, gene structure, conserved motifs, and gene ontology annotations. MeSOD genes were classified into two groups through phylogenetic analysis, revealing evolutionary connections. Promoters of these genes harbored stress-related cis-elements. Duplication analysis indicated the functional significance of MeCSOD2/MeCSOD4 and MeMSOD1/MeMSOD2. Through qRT-PCR, MeCSOD2 responded to salt stress, MeMSOD2 to drought, and cassava bacterial blight. Silencing MeMSOD2 increased XpmCHN11 virulence, indicating MeMSOD2 is essential for cassava's defense against XpmCHN11 infection. These findings enhance our understanding of the SOD gene family's role in cassava and contribute to strategies for stress tolerance improvement.

Sesame (Sesamum indicum L.) response to drought stress: susceptible and tolerant genotypes exhibit different physiological, biochemical, and molecular response patterns

Drought is one of the main environmental stresses affecting the quality and quantity of sesame production worldwide. The present study was conducted to investigate the effect of drought stress and subsequent re-watering on physiological, biochemical, and molecular responses of two contrasted sesame genotypes (susceptible vs. tolerant). Results showed that plant growth, photosynthetic rate, stomatal conductance, transpiration rate, and relative water content were negatively affected in both genotypes during water deficit. Both genotypes accumulated more soluble sugars, free amino acids, and proline and exhibited an increased enzyme activity for peroxidase, catalase, superoxide dismutase, and pyruvate dehydrogenase in response to drought damages including increased lipid peroxidation and membrane disruption. However, the tolerant genotype revealed a more extended root system and a more efficient photosynthetic apparatus. It also accumulated more soluble sugars (152%), free amino acids (48%), proline (75%), and antioxidant enzymes while showing lower electrolyte leakage (26%), lipid peroxidation (31%), and starch (35%) content, compared to the susceptible genotype at severe drought. Moreover, drought-related genes such as MnSOD1, MnSOD2, and PDHA-M were more expressed in the tolerant genotype, which encode manganese-dependent superoxide dismutase and the alpha subunit of pyruvate dehydrogenase, respectively. Upon re-watering, tolerant genotype recovered to almost normal levels of photosynthesis, carboxylation efficiency, lipid peroxidation, and electrolyte leakage, while susceptible genotype still suffered critical issues. Overall, these results suggest that a developed root system and an efficient photosynthetic apparatus along with the timely and effective accumulation of protective compounds enabled the tolerant sesame to withstand stress and successfully return to a normal growth state after drought relief. The findings of this study can be used as promising criteria for evaluating genotypes under drought stress in future sesame breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01372-y.

Comparative analysis of physiological variations and genetic architecture for cold stress response in soybean germplasm

Soybean (Glycine max L.) is susceptible to low temperatures. Increasing lines of evidence indicate that abiotic stress-responsive genes are involved in plant low-temperature stress response. However, the involvement of photosynthesis, antioxidants and metabolites genes in low temperature response is largely unexplored in Soybean. In the current study, a genetic panel of diverse soybean varieties was analyzed for photosynthesis, chlorophyll fluorescence and leaf injury parameters under cold stress and control conditions. This helps us to identify cold tolerant (V100) and cold sensitive (V45) varieties. The V100 variety outperformed for antioxidant enzymes activities and relative expression of photosynthesis (Glyma.08G204800.1, Glyma.12G232000.1), GmSOD (GmSOD01, GmSOD08), GmPOD (GmPOD29, GmPOD47), trehalose (GmTPS01, GmTPS13) and cold marker genes (DREB1E, DREB1D, SCOF1) than V45 under cold stress. Upon cold stress, the V100 variety showed reduced accumulation of H₂O₂ and MDA levels and subsequently showed lower leaf injury compared to V45. Together, our results uncovered new avenues for identifying cold tolerant soybean varieties from a large panel. Additionally, we identified the role of antioxidants, osmo-protectants and their posttranscriptional regulators miRNAs such as miR319, miR394, miR397, and miR398 in Soybean cold stress tolerance.

Superoxide dismutase (SOD) family in durum wheat: promising candidates for improving crop resilience

The SOD family has been extensively analyzed at genome wide level in several crops. However, little is known about this family in durum wheat. In this study, a total of 14 TdSOD genes were identified in whole durum wheat genome including 8 TdCu-ZnSODs, 2 TdMnSODs, and 4 TdFeSODs. In silico analysis evinced that TdSOD family members displayed a closer evolutionary relationship, similar gene structure and protein features with their homologs from other plant species. Furthermore, the analysis of their promoter regions revealed the presence of a great number of cis-regulatory elements related to plant development, abiotic and biotic stresses, phytohormones, and several potential binding sites for transcription factors. Interestingly, 3D structure analysis revealed that TdCu-ZnSOD2A-2 and TdCu-ZnSOD2B-2, belonging to the Cu-Zn group, were modeled as copper chaperone for SOD like their homologs from rice and Arabidopsis. The expression profile of eight TdSOD candidate genes was investigated under salt, drought, cold, and ABA treatments. Notably, TdCu-ZnSOD2A-1, TdFeSOD4A-1, and TdFeSOD7A-1 were significantly up-regulated under all stress treatments. On the other hand, TdCu-ZnSOD7B and TdMnSOD2B were strongly expressed in roots and leaves under cold stress and TdCu-ZnSOD2B-2 was particularly up-regulated in leaves under ABA treatment. Ultimately, these findings provide valuable information for the identification of attractive candidate genes to improve wheat resilience.

Genome-wide transcriptome profiling revealed biological macromolecules respond to low temperature stress in Brassica napus L

Brassica napus L. (B. napus) is a vital oilseed crop cultivated worldwide; low temperature (LT) is one of the major stress factors that limit its growth, development, distribution, and production. Even though processes have been developed to characterize LT-responsive genes, only limited studies have exploited the molecular response mechanisms in B. napus. Here the transcriptome data of an elite B. napus variety with LT adaptability was acquired and applied to investigate the gene expression profiles of B. napus in response to LT stress. The bioinformatics study revealed a total of 79,061 unigenes, of which 3,703 genes were differentially expressed genes (DEGs), with 2,129 upregulated and 1,574 downregulated. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis pinpointed that the DEGs were enriched in LT-stress-responsive biological functions and metabolic pathways, which included sugar metabolism, antioxidant defense system, plant hormone signal transduction, and photosynthesis. Moreover, a group of LT-stress-responsive transcription factors with divergent expression patterns under LT was summarized. A combined protein interaction suggested that a complex interconnected regulatory network existed in all detected pathways. RNA-seq data was verified using real-time quantitative polymerase chain reaction analysis. Based on these findings, we presented a hypothesis model illustrating valuable information for understanding the LT response mechanisms in B. napus.

Genome-Wide Characterization of Superoxide Dismutase (SOD) Genes in Daucus carota: Novel Insights Into Structure, Expression, and Binding Interaction With Hydrogen Peroxide (H2O2) Under Abiotic Stress Condition

Superoxide dismutase (SOD) proteins are important antioxidant enzymes that help plants to grow, develop, and respond to a variety of abiotic stressors. SOD gene family has been identified in a number of plant species but not yet in Daucus carota. A total of 9 DcSOD genes, comprising 2 FeSODs, 2 MnSODs, and 5 Cu/ZnSODs, are identified in the complete genome of D. carota, which are dispersed in five out of nine chromosomes. Based on phylogenetic analysis, SOD proteins from D. carota were categorized into two main classes (Cu/ZnSODs and MnFeSODs). It was predicted that members of the same subgroups have the same subcellular location. The phylogenetic analysis was further validated by sequence motifs, exon-intron structure, and 3D protein structures, with each subgroup having a similar gene and protein structure. Cis-regulatory elements responsive to abiotic stresses were identified in the promoter region, which may contribute to their differential expression. Based on RNA-seq data, tissue-specific expression revealed that DcCSD2 had higher expression in both xylem and phloem. Moreover, DcCSD2 was differentially expressed in dark stress. All SOD genes were subjected to qPCR analysis after cold, heat, salt, or drought stress imposition. SODs are antioxidants and play a critical role in removing reactive oxygen species (ROS), including hydrogen peroxide (H2O2). DcSODs were docked with H2O2 to evaluate their binding. The findings of this study will serve as a basis for further functional insights into the DcSOD gene family.

Molecular and Biochemical Analysis of Duplicated Cytosolic CuZn Superoxide Dismutases of Rice and in silico Analysis in Plants

Superoxide dismutases (SODs, EC 1.15.1.1) are ubiquitous antioxidant metalloenzymes important for oxidative stress tolerance and cellular redox environment. Multiple factors have contributed toward the origin and diversity of SOD isoforms among different organisms. In plants, the genome duplication events, responsible for the generation of multiple gene copies/gene families, have also contributed toward the SOD diversity. However, the importance of such molecular events on the characteristics of SODs has not been studied well. This study investigated the effects of divergence on important characteristics of two block-duplicated rice cytosolic CuZn SODs (OsCSD1, OsCSD4), along with in silico assessment of similar events in other plants. The analysis revealed heterogeneity in gene length, regulatory regions, untranslated regions (UTRs), and coding regions of two OsCSDs. An inconsistency in the database-predicted OsCSD1 gene structure was also identified and validated experimentally. Transcript analysis showed differences in the basal levels and stress responsiveness of OsCSD1 and OsCSD4, and indicated the presence of two transcription start sites in the OsCSD1. At the amino acid level, the two OsCSDs showed differences at 18 sites; however, both exist as a homodimer, displaying typical CuZn SOD characteristics, and enhancing the oxidative stress tolerance of Escherichia coli cells. However, OsCSD4 showed higher specific activity as well as stability. The comparison of the two OsCSDs with reported thermostable CSDs from other plants identified regions likely to be associated with stability, while the homology modeling and superposition highlighted structural differences. The two OsCSDs displayed heteromeric interaction capability and forms an enzymatically active heterodimer (OsCSD1:OsCSD4) on co-expression, which may have significance as both are cytosolic. In silico analysis of 74 plant genomes revealed the prevalence of block duplications for multiple CSD copies (mostly cytosolic). The divergence and clustering analysis of CSDs suggested the possibility of an ancestral duplication event in monocots. Conserved SOD features indicating retention of SOD function among CSD duplicates were evident in few monocots and dicots. In most other species, the CSD copies lacked critical features and may not harbor SOD function; however, other feature-associated functions or novel functions might be present. These aspects of divergent CSD copies encoding co-localized CSDs may have implications in plant SOD functions in the cytosol and other organelles.

The Superoxide Dismutase Gene Family in Nicotiana tabacum: Genome-Wide Identification, Characterization, Expression Profiling and Functional Analysis in Response to Heavy Metal Stress

Superoxide dismutases (SODs) play an important role in protecting plants against ROS toxicity induced by biotic and abiotic stress. Recent studies have shown that the SOD gene family is involved in plant growth and development; however, knowledge of the SOD gene family in tobacco is still limited. In the present study, the SOD gene family was systematically characterized in the tobacco genome. Based on the conserved motif and phylogenetic tree, 15 NtSOD genes were identified and classified into three subgroups, including 5 NtCSDs, 7 NtFSDs and 3 NtMSDs. The predicted results of the transport peptide or signal peptide were consistent with their subcellular localization. Most NtSOD genes showed relatively well-maintained exon-intron and motif structures in the same subgroup. An analysis of cis-acting elements in SOD gene promoters showed that NtSOD expression was regulated by plant hormones, defense and stress responses, and light. In addition, multiple transcription factors and miRNAs are predicted to be involved in the regulation of NtSOD gene expression. The qPCR results indicated specific spatial and temporal expression patterns of the NtSOD gene family in different tissues and developmental stages, and this gene family played an important role in protecting against heavy metal stress. The results of functional complementation tests in the yeast mutant suggested that NtCSD1a, NtFSD1e and NtMSD1b scavenge ROS produced by heavy metal stress. This study represents the first genome-wide analysis of the NtSOD gene family, which lays a foundation for a better understanding of the function of the NtSOD gene family and improving the tolerance of plants to heavy metal toxicity.

Genome-Wide Evolution and Comparative Analysis of Superoxide Dismutase Gene Family in Cucurbitaceae and Expression Analysis of Lagenaria siceraria Under Multiple Abiotic Stresses

Superoxide dismutase (SOD) is an important enzyme that serves as the first line of defense in the plant antioxidant system and removes reactive oxygen species (ROS) under adverse conditions. The SOD protein family is widely distributed in the plant kingdom and plays a significant role in plant growth and development. However, the comprehensive analysis of the SOD gene family has not been conducted in Cucurbitaceae. Subsequently, 43 SOD genes were identified from Cucurbitaceae species [Citrullus lanatus (watermelon), Cucurbita pepo (zucchini), Cucumis sativus (cucumber), Lagenaria siceraria (bottle gourd), Cucumis melo (melon)]. According to evolutionary analysis, SOD genes were divided into eight subfamilies (I, II, III, IV, V, VI, VII, VIII). The gene structure analysis exhibited that the SOD gene family had comparatively preserved exon/intron assembly and motif as well. Phylogenetic and structural analysis revealed the functional divergence of Cucurbitaceae SOD gene family. Furthermore, microRNAs 6 miRNAs were predicted targeting 3 LsiSOD genes. Gene ontology annotation outcomes confirm the role of LsiSODs under different stress stimuli, cellular oxidant detoxification processes, metal ion binding activities, SOD activity, and different cellular components. Promoter regions of the SOD family revealed that most cis-elements were involved in plant development, stress response, and plant hormones. Evaluation of the gene expression showed that most SOD genes were expressed in different tissues (root, flower, fruit, stem, and leaf). Finally, the expression profiles of eight LsiSOD genes analyzed by qRT-PCR suggested that these genetic reserves responded to drought, saline, heat, and cold stress. These findings laid the foundation for further study of the role of the SOD gene family in Cucurbitaceae. Also, they provided the potential for its use in the genetic improvement of Cucurbitaceae.

Methyl jasmonate-induced resistance to Delia platura (Diptera: Anthomyiidae) in Lupinus mutabilis

BACKGROUND

Andean lupin (Lupinus mutabilis Sweet) is an important leguminous crop from South America with a high protein content. In Ecuador, lupin yields are severely affected by the infestation of Delia platura larvae on germinating seeds. The application of elicitor molecules with activity against herbivorous insects to control D. platura infestation constitutes an unexplored and promising alternative for chemical insecticides. In this study, methyl jasmonate (MeJA), hexanoic acid, menadione sodium bisulfite, and DL-β-aminobutyric acid were evaluated for their ability to induce resistance against D. platura in three commercial lupin cultivars.

RESULTS

Only seeds pretreated with MeJA significantly impaired insect performance during choice and no-choice assays. Additionally, fitness indicators such as seed germination and growth were not affected by MeJA treatment. To investigate the molecular mechanisms behind the MeJA-mediated resistance, RT-qPCR assays were performed. First, RT-qPCR reference genes were validated, showing that LmUBC was the most stable reference gene. Next, expression analysis over time revealed that MeJA application up-regulated the activity of the jasmonic acid biosynthetic genes LmLOX2 and LmAOS, together with other jasmonate-related defense genes, such as LmTPS1, LmTPS4, LmPI2, LmMBL, LmL/ODC, LmCSD1, and LmPOD.

CONCLUSION

This study indicates that MeJA can be used as an environmentally friendly elicitor molecule to protect Andean lupin from D. platura attack without fitness cost. MeJA application induces plant defense responses to insects in Andean lupin that may be modulated by the onset of terpenoid biosynthesis, proteinase inhibitors, lectins, polyamines, and antioxidative enzymes. © 2021 Society of Chemical Industry.

Genome-wide analysis of the maize superoxide dismutase (SOD) gene family reveals important roles in drought and salt responses

Superoxide dismutase proteins (SODs) are antioxidant enzymes with important roles in abiotic stress responses. The SOD gene family has been systematically analyzed in many plants; however, it is still poorly understood in maize. Here, a bioinformatics analysis of maize SOD gene family was conducted by describing gene structure, conserved motifs, phylogenetic relationships, gene duplications, promoter cis-elements and GO annotations. In total, 13 SOD genes were identified in maize and five members were involved in segmental duplication. Phylogenetic analysis indicated that SODs from maize and other plants comprised two groups, which could be further classified into different subgroups, with most members in the same subgroup having the same subcellular localization. The ZmSOD promoters contained 2-10 stress-responsive cis-elements with different distributions. Heatmap analysis indicated that ZmSODs were expressed in most of the detected tissues and organs. The expression patterns of ZmSODs were investigated under drought and salt treatments by qRT-PCR, and most members were responsive to drought or salt stress, especially some ZmSODs with significant expression changes were identified, such as ZmCSD2 and ZmMSD2, suggesting the important roles of ZmSODs in abiotic stress responses. Our results provide an important basis for further functional study of ZmSODs in future study.

Genome-Wide Analysis and Expression Profile of Superoxide Dismutase (SOD) Gene Family in Rapeseed (Brassica napus L.) under Different Hormones and Abiotic Stress Conditions

Superoxide dismutase (SOD) is an important enzyme that acts as the first line of protection in the plant antioxidant defense system, involved in eliminating reactive oxygen species (ROS) under harsh environmental conditions. Nevertheless, the SOD gene family was yet to be reported in rapeseed (Brassica napus L.). Thus, a genome-wide investigation was carried out to identify the rapeseed SOD genes. The present study recognized 31 BnSOD genes in the rapeseed genome, including 14 BnCSDs, 11 BnFSDs, and six BnMSDs. Phylogenetic analysis revealed that SOD genes from rapeseed and other closely related plant species were clustered into three groups based on the binding domain with high bootstrap values. The systemic analysis exposed that BnSODs experienced segmental duplications. Gene structure and motif analysis specified that most of the BnSOD genes displayed a relatively well-maintained exon-intron and motif configuration within the same group. Moreover, we identified five hormones and four stress- and several light-responsive cis-elements in the promoters of BnSODs. Thirty putative bna-miRNAs from seven families were also predicted, targeting 13 BnSODs. Gene ontology annotation outcomes confirm the BnSODs role under different stress stimuli, cellular oxidant detoxification processes, metal ion binding activities, SOD activity, and different cellular components. Twelve BnSOD genes exhibited higher expression profiles in numerous developmental tissues, i.e., root, leaf, stem, and silique. The qRT-PCR based expression profiling showed that eight genes (BnCSD1, BnCSD3, BnCSD14, BnFSD4, BnFSD5, BnFSD6, BnMSD2, and BnMSD10) were significantly up-regulated under different hormones (ABA, GA, IAA, and KT) and abiotic stress (salinity, cold, waterlogging, and drought) treatments. The predicted 3D structures discovered comparable conserved BnSOD protein structures. In short, our findings deliver a foundation for additional functional investigations on the BnSOD genes in rapeseed breeding programs.

Genome-wide analysis of the superoxide dismutase (SOD) gene family in Zostera marina and expression profile analysis under temperature stress

Superoxide dismutases (SODs) serve as the first line of defense in the plant antioxidant enzyme system, and play a primary role in the removal of reactive oxygen species (ROS). However, our understanding of the functions of the SOD family in Zostera marina is limited. In this study, a systematic analysis was conducted on the characteristics of the SOD genes in Z. marina at the whole-genome level. Five SOD genes were identified, consisting of two Cu/ZnSODs, two FeSODs, and one MnSOD. Phylogenetic analysis showed that ZmSOD proteins could be divided into two major categories (Cu/ZnSODs and Fe-MnSODs). Sequence motifs, gene structure, and the 3D-modeled protein structures further supported the phylogenetic analysis, with each subgroup having similar motifs, exon-intron structures, and protein structures. Additionally, several cis-elements were identified that may respond to biotic and abiotic stresses. Transcriptome analysis revealed expression diversity of ZmSODs in various tissues. Moreover, qRT-PCR analysis showed that the expression level of most ZmSOD genes trended to decreased expression with the increase of temperature, indicating that heat stress inhibits expression of ZmSODs and may result in reduced ability of ZmSODs to scavenge ROS. Our results provide a basis for further functional research on the SOD gene family in Z. marina, which will help to determine the molecular mechanism of ZmSOD genes in response to environmental stress.

Genome-wide identification and characterization of the soybean SOD family during alkaline stress

Background

Superoxide dismutase (SOD) proteins, as one kind of the antioxidant enzymes, play critical roles in plant response to various environment stresses. Even though its functions in the oxidative stress were very well characterized, the roles of SOD family genes in regulating alkaline stress response are not fully reported.

Methods

We identified the potential family members by using Hidden Markov model and soybean genome database. The neighbor-joining phylogenetic tree and exon-intron structures were generated by using software MEGA 5.0 and GSDS online server, respectively. Furthermore, the conserved motifs were analyzed by MEME online server. The syntenic analysis was conducted using Circos-0.69. Additionally, the expression levels of soybean SOD genes under alkaline stress were identified by qRT-PCR.

Results

In this study, we identified 13 potential SOD genes in soybean genome. Phylogenetic analysis suggested that SOD genes could be classified into three subfamilies, including MnSODs (GmMSD1-2), FeSODs (GmFSD1-5) and Cu/ZnSODs (GmCSD1-6). We further investigated the gene structure, chromosomal locations and gene-duplication, conserved domains and promoter cis-elements of the soybean SOD genes. We also explored the expression profiles of soybean SOD genes in different tissues and alkaline, salt and cold stresses, based on the transcriptome data. In addition, we detected their expression patterns in roots and leaves by qRT-PCR under alkaline stress, and found that different SOD subfamily genes may play different roles in response to alkaline stress. These results also confirmed the hypothesis that the great evolutionary divergence may contribute to the potential functional diversity in soybean SOD genes. Taken together, we established a foundation for further functional characterization of soybean SOD genes in response to alkaline stress in the future.

The CYP74 Gene Family in Watermelon: Genome-Wide Identification and Expression Profiling Under Hormonal Stress and Root-Knot Nematode Infection

Allene oxide synthase (AOS) and hydroperoxide lyase (HPL), members of the CYP74 gene family, are branches of the oxylipin pathway and play vital roles in plant responses to a number of stresses. In this study, four HPL genes and one AOS gene were identified in the watermelon genome, which were clustered into three subfamilies (CYP74A, CYP74B and CYP74C). Sequence analysis revealed that most HPL and AOS proteins from various plants contain representative domains, including Helix-I region, Helix-K region (ExxR) and Heme-binding domain. A number of development-, stress-, and hormone-related cis-elements were found in the promoter regions of the ClAOS and ClHPL genes, and the detected ClAOS and ClHPL genes were differentially expressed in different tissues and fruit development stages, as well as in response to various hormones. In addition, red light could enhance the expression of ClAOS in root-knot nematode-infected leaves and roots of watermelon, implying that ClAOS might play a primary role in red light-induced resistance against root-knot nematodes. These findings lay a foundation for understanding the specific function of CYP74 genes in watermelon.

Highlight on the expression and the function of a novel MnSOD from diploid wheat (T. monococcum) in response to abiotic stress and heavy metal toxicity

Superoxide dismutases (SODs) play a pivotal role in improving abiotic stress tolerance in plant cells. A novel manganese superoxide dismutase gene, denoted as TmMnSOD, was identified from Triticum monococcum. The encoded protein displayed high sequence identity with MnSOD family members and was highly homologous to TdMnSOD from durum wheat. Furthermore, the 3D structure analysis revealed that TmMnSOD displayed homotetramer subunit organization, incorporating four Mn²⁺ ions. Notably, TmMnSOD structure contains predominantly alpha helices with three beta sheets. On the other hand, under stress conditions, TmMnSOD transcript level was significantly up-regulated by salt, oxidative and heavy metal stresses. At the functional level, TmMnSOD imparts tolerance of yeast and E. coli cells under diverse stresses. Promoter analysis of TmMnSOD gene showed the presence of a great number of salt and pathogen-responsive cis-regulatory elements, highlighting the interest of this gene in breeding programs towards improved tolerance to salt stress in wheat.

Genome-wide identification, characterization, and expression analysis of superoxide dismutase (SOD) genes in foxtail millet (Setaria italica L.)

Superoxide dismutases (SODs) play important roles in plant growth, development, and response to abiotic stresses. Despite SOD gene families have been identified in various plant species, little is known in foxtail millet (Setaria italica L.). In this study, a systematic analysis of SOD gene family was performed in foxtail millet and the expression pattern of SOD genes in response to abiotic stressors was analyzed at the whole-genomic level. Eight SOD genes were identified in foxtail millet, including 4 Cu/ZnSODs, 3 FeSODs, and 1 MnSOD. These SiSODs are unevenly distributed across 5 of the 9 chromosomes. Phylogenetic analysis showed that SOD proteins could be divided into two major categories (Cu/ZnSODs and Fe-MnSODs), containing seven subgroups, from foxtail millet and other plant species. SOD genes have conserved motif and exon/intron composition in the same subgroup among Setaria italica, Setaria viridis, and Oryza sativa. Additionally, many cis-elements that respond to different stressors were distributed at different densities in the promoters of 8 SiSODs. The expression patterns of SiSODs in different tissues and different abiotic stressors indicated that the SiSODs may play important roles in reactive oxygen species scavenging, caused by various stressors in foxtail millet. This study provides a foundation for the further cloning and functional verification of the SOD gene family response to environmental stimuli in foxtail millet.