Pungency related gene network in Allium sativum L., response to sulfur treatments

Pungency of garlic (Allium sativum L.) is generated from breakdown of the alk(en)yl cysteine sulphoxide (CSO), alliin and its subsequent breakdown to allicin under the activity of alliinase (All). Based on recent evidence, two other important genes including Sulfite reductase (SiR) and Superoxide dismutase (SOD) are thought to be related to sulfur metabolism. These three gene functions are in sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. In this research, the order and priority of simultaneous expression of three genes including All, SiR and SOD were measured on some garlic ecotypes of Iran, collected from Zanjan, Hamedan and Gilan, provinces under sulfur concentrations (0, 6, 12, 24 and 60 g/ per experimental unit: pot) using real-time quantitative PCR (RT-qPCR) analysis. For understanding the network interactions between studied genes and other related genes, in silico gene network analysis was constructed to investigate various mechanisms underlying stimulation of A. sativum L. to cope with imposed sulfur. Complicated network including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-networks to have a deeper insight. Analysis of q-RT-PCR data revealed the highest expression in All and SiR genes respectively. To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. Complex sub-network exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed their involvement in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic. The findings of this research can be used to extend further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency for health benefits and industry demands.