Genome-wide characterization, expression profiling, and post-transcriptional study of GASA gene family

GhGASA14 regulates the flowering time of upland cotton in response to GA3

KEY MESSAGE

The silencing of GhGASA14 and the identification of superior allelic variation in its coding region indicate that GhGASA14 may positively regulate flowering and the response to GA3. Gibberellic acid-stimulated Arabidopsis (GASA), a member of the gibberellin-regulated short amino acid family, has been extensively investigated in several plant species and found to be critical for plant growth and development. However, research on this topic in cotton has been limited. In this study, we identified 38 GhGASAs that were dispersed across 18 chromosomes in upland cotton, and all of these genes had a GASA core domain. Transcriptome expression patterns and qRT-PCR results revealed that GhGASA9 and GhGASA14 exhibited upregulated expression not only in the floral organs but also in the leaves of early-maturing cultivars. The two genes were functionally characterized by virus-induced gene silencing (VIGS), and the budding and flowering times after silencing the target genes were later than those of the control (TRV:00). Compared with that in the water-treated group (MOCK), the flowering period of the different fruiting branches in the GA3-treated group was more concentrated. Interestingly, allelic variation was detected in the coding sequence of GhGASA14 between early-maturing and late-maturing accessions, and the frequency of this favorable allele was greater in high-latitude cotton cultivars than in low-latitude ones. Additionally, a significant linear relationship was observed between the expression level of GhGASA14 and flowering time among the 12 upland cotton accessions. Taken together, these results indicated that GhGASA14 may positively regulate flowering time and respond to GA3. These findings could lead to the use of valuable genetic resources for breeding early-maturing cotton cultivars in the future.

Genome-Wide Identification of GAST Family Members and Their Potential Roles in Epicotyl Dormancy in Chinese Cork Oak (Quercus variabilis)

Chinese cork oak (Quercus variabilis Blume) is a widespread tree species with high economic and ecological values. Chinese cork oak exhibits epicotyl dormancy, causing emergence heterogeneity and affecting the quality of seedling cultivation. Gibberellic acid-stimulated transcript (GAST) is a plant-specific protein family that plays a crucial regulatory role in plant growth, development, and seed germination. However, their evolution in Chinese cork oak and roles in epicotyl dormancy are still unclear. Here, a genome-wide identification of the GAST gene family was conducted in Chinese cork oak. Ten QvGAST genes were identified, and nine of them were expressed in seed. The physicochemical properties and promoter cis-acting elements of the selected Chinese cork oak GAST family genes indicated that the cis-acting elements in the GAST promoter are involved in plant development, hormone response, and stress response. Germinated seeds were subjected to gibberellins (GAs), abscisic acid (ABA), and fluridone treatments to show their response during epicotyl dormancy release. Significant changes in the expression of certain QvGAST genes were observed under different hormone treatments. QvGAST1, QvGAST2, QvGAST3, and QvGAST6 exhibited upregulation in response to gibberellin. QvGAST2 was markedly upregulated during the release of epicotyl dormancy in response to GA. These findings suggested that QvGAST2 might play an important role in epicotyl dormancy release. This study provides a basis for further analysis of the mechanisms underlying the alleviation of epicotyl dormancy in Chinese cork oak by QvGASTs genes.

Genome-wide analysis of the Tritipyrum NAC gene family and the response of TtNAC477 in salt tolerance

NAC transcription factors are widely distributed in the plant kingdom and play an important role in the response to various abiotic stresses in plant species. Tritipyrum, an octoploid derived from hybridization of Triticum aestivum (AABBDD) and Thinopyrum elongatum (EE), is an important genetic resource for integrating the desirable traits of Th. elongatum into wheat. In this study, we investigated the tissue distribution and expression of Tritipyrum NAC genes in the whole genomes of T. aestivum and Th. elongatum after obtaining their complete genome sequences. Based on phylogenetic relationships, conserved motifs, gene synthesis, evolutionary analysis, and expression patterns, we identified and characterized 732 Tritipyrum NAC genes. These genes were divided into six main groups (A, B, C, D, E, and G) based on phylogenetic relationships and evolutionary studies, with members of these groups sharing the same motif composition. The 732 TtNAC genes are widely distributed across 28 chromosomes and include 110 duplicated genes. Gene synthesis analysis indicated that the NAC gene family may have a common ancestor. Transcriptome data and quantitative polymerase chain reaction (qPCR) expression profiles showed 68 TtNAC genes to be highly expressed in response to various salt stress and recovery treatments. Tel3E01T644900 (TtNAC477) was particularly sensitive to salt stress and belongs to the same clade as the salt tolerance genes ANAC019 and ANAC055 in Arabidopsis. Pearson correlation analysis identified 751 genes that correlated positively with expression of TtNAC477, and these genes are enriched in metabolic activities, cellular processes, stimulus responses, and biological regulation. TtNAC477 was found to be highly expressed in roots, stems, and leaves in response to salt stress, as confirmed by real-time PCR. These findings suggest that TtNAC477 is associated with salt tolerance in plants and might serve as a valuable exogenous gene for enhancing salt tolerance in wheat.

Genome-Wide Identification and Comprehensive Analysis of the GASA Gene Family in Peanuts (Arachis hypogaea L.) under Abiotic Stress

Peanut (Arachis hypogaea L.) is a globally cultivated crop of significant economic and nutritional importance. The role of gibberellic-acid-stimulated Arabidopsis (GASA) family genes is well established in plant growth, development, and biotic and abiotic stress responses. However, there is a gap in understanding the function of GASA proteins in cultivated peanuts, particularly in response to abiotic stresses such as drought and salinity. Thus, we conducted comprehensive in silico analyses to identify and verify the existence of 40 GASA genes (termed AhGASA) in cultivated peanuts. Subsequently, we conducted biological experiments and performed expression analyses of selected AhGASA genes to elucidate their potential regulatory roles in response to drought and salinity. Phylogenetic analysis revealed that AhGASA genes could be categorized into four distinct subfamilies. Under normal growth conditions, selected AhGASA genes exhibited varying expressions in young peanut seedling leaves, stems, and roots tissues. Notably, our findings indicate that certain AhGASA genes were downregulated under drought stress but upregulated under salt stress. These results suggest that specific AhGASA genes are involved in the regulation of salt or drought stress. Further functional characterization of the upregulated genes under both drought and salt stress will be essential to confirm their regulatory roles in this context. Overall, our findings provide compelling evidence of the involvement of AhGASA genes in the mechanisms of stress tolerance in cultivated peanuts. This study enhances our understanding of the functions of AhGASA genes in response to abiotic stress and lays the groundwork for future investigations into the molecular characterization of AhGASA genes.

Genome-Wide Identification and Characterization of Gibberellic Acid-Stimulated Arabidopsis Gene Family in Pineapple (Ananas comosus)

The gibberellic acid-stimulated Arabidopsis (GASA) gene family plays a crucial role in growth, development, and stress response, and it is specific to plants. This gene family has been extensively studied in various plant species, and its functional role in pineapple has yet to be characterized. In this study, 15 AcGASA genes were identified in pineapple through a genome-wide scan and categorized into three major branches based on a phylogenetic tree. All AcGASA proteins share a common structural domain with 12 cysteine residues, but they exhibit slight variations in their physicochemical properties and motif composition. Predictions regarding subcellular localization suggest that AcGASA proteins are present in the cell membrane, Golgi apparatus, nucleus, and cell wall. An analysis of gene synteny indicated that both tandem and segmental repeats have a significant impact on the expansion of the AcGASA gene family. Our findings demonstrate the differing regulatory effects of these hormones (GA, NAA, IAA, MeJA, and ABA) on the AcGASA genes. We analyzed the expression profiles of GASA genes in different pineapple tissue parts, and the results indicated that AcGASA genes exhibit diverse expression patterns during the development of different plant tissues, particularly in the regulation of floral organ development. This study provides a comprehensive understanding of GASA family genes in pineapple. It serves as a valuable reference for future studies on the functional characterization of GASA genes in other perennial herbaceous plants.

Genome-wide analysis of the AP2/ERF gene family in Tritipyrum and the response of TtERF_B2-50 in salt-tolerance

The AP2/ERF transcription factor is widely distributed across the plant kingdom and plays a crucial role in various abiotic stress responses in plants. Tritipyrum, an octoploid resulting from an intergeneric cross between Triticum aestivum (AABBDD) and Thinopyrum elongatum (EE), is a valuable source of germplasm for incorporating superior traits of Th. elongatum into T. aestivum. With the recent availability of whole -genome sequences for T. aestivum and Th. elongatum, we explored the organization and expression profiling of Tritipyrum AP2/ERF genes across the entire genome. Our investigation identified 543 Tritipyrum AP2/ERF genes, which evolutionary analysis categorized into four major groups (AP2, DREB, ERF, and RAV), whose members share a conserved motif composition. These 543 TtAP2/ERF genes were distributed throughout 28 chromosomes, with 132 duplications. Synteny analysis suggests that the AP2/ERF gene family may have a common ancestor. Transcriptome data and Real-Time PCR expression profiles revealed 43 TtAP2/ERF genes with high expression levels in response to various salt stressors and recovery regimens. Tel2E01T236300 (TtERF_B2-50) was particularly salt stress-sensitive and evolutionarily related to the salt-tolerant gene AtERF7 in A. thaliana. Pearson correlation analysis identified 689 genes positively correlated (R > 0.9) with TtERF_B2-50 expression, enriched in metabolic activities, cellular processes, stimulus response, and biological regulation. Real-time PCR showed that TtERF_B2-50 was highly expressed in roots, stems, and leaves under salt stress. These findings suggest that TtERF_B2-50 may be associated with salt stress tolerance and may serve as a valuable foreign gene for enhancing salt tolerance in wheat.

GASA Proteins: Review of Their Functions in Plant Environmental Stress Tolerance

Gibberellic acid-stimulated Arabidopsis (GASA) gene family is a class of functional cysteine-rich proteins characterized by an N-terminal signal peptide and a C-terminal-conserved GASA domain with 12 invariant cysteine (Cys) residues. GASA proteins are widely distributed among plant species, and the majority of them are involved in the signal transmission of plant hormones, the regulation of plant development and growth, and the responses to different environmental constraints. To date, their action mechanisms are not completely elucidated. This review reports an overview of the diversity, structure, and subcellular localization of GASA proteins, their involvement in hormone crosstalk and redox regulation during development, and plant responses to abiotic and biotic stresses. Knowledge of this complex regulation can be a contribution to promoting multiple abiotic stress tolerance with potential agricultural applications through the engineering of genes encoding GASA proteins and the production of transgenic plants.

Genome-wide identification and expression profiling of snakin/GASA genes under drought stress in barley (Hordeum vulgare L.)

Gibberellic Acid-Stimulated Arabidopsis (GASA) proteins are present in various plants and have a role in plant growth, stress responses, and hormone crosstalk. GASA coding sequences in barley were discovered in this study. We then investigated gene and protein structure, physicochemical characteristics, evolutionary and phylogenetic relationships, promoter region, post-translational modification, and in silico gene expression. Finally, real-time quantitative PCR (RT-qPCR) was used to examine the expression of GASA genes in root and shoot tissues under drought stress. We found 11 GASA genes spread across six of seven chromosomes in the barley genome. A conserved GASA domain and 12-cysteine residues at the C-terminus were included in the proteins. All GASA genes contained secretory signal peptides. The GASA genes in Hordeum vulgare (HvGASA) have been classified into three subfamilies based on evolutionary analysis. According to synteny analyses, segmental duplications are significant in forming the GASA gene family. According to the cis-elements analyses, GASA genes may be induced by a variety of phytohormones and stresses. Tissue-specific expression analysis indicated that GASA genes had varied expression patterns in different tissues. Contrary to common perception, the expression study of GASA genes under biotic and abiotic stresses revealed that GASA genes are more induced by abiotic stresses than biotic stresses. The qPCR confirmed the response of GASA genes to abiotic stresses and showed different expression patterns of these genes under drought stress. Overall, these results can improve our knowledge about the function of GASA genes and provide data for future researches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03545-8.

Genome-wide analysis of the laccase (LAC) gene family in Aeluropus littoralis: A focus on identification, evolution and expression patterns in response to abiotic stresses and ABA treatment

Laccases are plant enzymes with essential functions during growth and development. These monophenoloxidases are involved in lignin polymerization, and their expression respond to environmental stress. However, studies of laccases in some plants and fungi have highlighted that many structural and functional aspects of these genes are still unknown. Here, the laccase gene family in Aeluropus littoralis (AlLAC) is described based on sequence structure and expression patterns under abiotic stresses and ABA treatment. Fifteen non-redundant AlLACs were identified from the A. littoralis genome, which showed differences in physicochemical characteristics and gene structure. Based on phylogenetic analysis, AlLACs and their orthologues were classified into five groups. A close evolutionary relationship was observed between LAC gene family members in rice and A. littoralis. According to the interaction network, AlLACs interact more with proteins involved in biological processes such as iron incorporation into the metallo-sulfur cluster, lignin catabolism, regulation of the symbiotic process and plant-type primary cell wall biogenesis. Gene expression analysis of selected AlLACs using real-time RT (reverse transcription)-PCR revealed that AlLACs are induced in response to abiotic stresses such as cold, salt, and osmotic stress, as well as ABA treatment. Moreover, AlLACs showed differential expression patterns in shoot and root tissues. Our findings indicate that AlLACs are preferentially involved in the late response of A. littoralis to abiotic stress.

Genome-wide identification and characterization of the lettuce GASA family in response to abiotic stresses

BACKGROUND

Lettuce is one of the most extensively farmed vegetables in the world, and it prefers cool growing conditions. High temperatures promote premature bolt formation, reducing quality and yield. The gibberellic acid-stimulated Arabidopsis (GASA) family genes play critical roles in plant growth, development, and stress responses. However, the biological functions of GASA proteins in lettuce have yet to be thoroughly investigated.

RESULTS

Using genome-wide analysis, 20 GASAs were identified in lettuce including, three groups of LsGASA proteins based on the phylogenetic analysis. Except for one, all GASA proteins included a conserved GASA domain with 12 cysteine residues. Cis-element analysis showed that LsGASAs were closely associated with light, phytohormones, and stress resistance. Five segmental and three tandem duplication events were observed in the LsGASA family based on duplication analysis. GASA synteny analysis among lettuce, Arabidopsis, tobacco, and rice revealed that LsGASA5 is highly collinear with all species. Six of the 20 LsGASA showed increased expression patterns at specific time points in the shoot apical meristem when subjected to heat stress. According to gene expression analysis, the majority of GASA were highly expressed in flowers compared to other organs, and six GASA exhibited highly increased expression levels in response to NaCl, abscisic acid, and gibberellin treatment. Furthermore, LsGASA proteins are predominantly found in the plasma membrane and/or the cytosol.

CONCLUSIONS

This study provides a comprehensive characterization of LsGASA genes for their diversity and biological functions. Moreover, our results will be useful for further studies on the function of lettuce GASA in abiotic stress- and heat-induced bolting signaling.

Genome-Wide Characterization of the Sulfate Transporter Gene Family in Oilseed Crops: Camelina sativa and Brassica napus

Sulfate transporters (SULTRs) are responsible for the uptake of sulfate (SO₄^2-) ions in the rhizosphere by roots and their distribution to plant organs. In this study, SULTR family members in the genomes of two oilseed crops (Camelina sativa and Brassica napus) were identified and characterized based on their sequence structures, duplication events, phylogenetic relationships, phosphorylation sites, and expression levels. In total, 36 and 45 putative SULTR genes were recognized in the genomes of C. sativa and B. napus, respectively. SULTR proteins were predicted to be basophilic proteins with low hydrophilicity in both studied species. According to the observed phylogenetic relationships, we divided the SULTRs into five groups, out of which the SULTR 3 group showed the highest variation. Additionally, several duplication events were observed between the SULTRs. The first duplication event occurred approximately five million years ago between three SULTR 3.1 genes in C. sativa. Furthermore, two subunits were identified in the 3D structures of the SULTRs, which demonstrated that the active binding sites differed between C. sativa and B. napus. According to the available RNA-seq data, the SULTRs showed diverse expression levels in tissues and diverse responses to stimuli. SULTR 3 was expressed in all tissues. SULTR 3.1 was more upregulated in response to abiotic stresses in C. sativa, while SULTR 3.3 and SULTR 2.1 were upregulated in B. napus. Furthermore, SULTR 3 and SULTR 4.1 were upregulated in response to biotic stresses in B. napus. Additionally, the qPCR data showed that the SULTRs in C. sativa were involved in the plant's response to salinity. Based on the distribution of cis-regulatory elements in the promoter region, we speculated that SULTRs might be controlled by phytohormones, such as ABA and MeJA. Therefore, it seems likely that SULTR genes in C. sativa have been more heavily influenced by evolutionary processes and have acquired further diversity. The results reveal new insights of the structures and functions of SULTRs in oilseed crops. However, further analyses, related to functional studies, are needed to uncover the role of SULTRs in the plants' development and growth processes, as well as in their response to stimuli.

Genome-Wide Analysis of Aquaporin Gene Family in Triticum turgidum and Its Expression Profile in Response to Salt Stress

During the response of plants to water stresses, aquaporin (AQP) plays a prominent role in membrane water transport based on the received upstream signals. Due to the importance of the AQP gene family, studies have been conducted that investigate the function and regulatory system of these genes. However, many of their molecular aspects are still unknown. This study aims to carry out a genome-wide investigation of the AQP gene family in Triticum turgidum using bioinformatics tools and to investigate the expression patterns of some members in response to salt stress. Our results show that there are 80 TtAQP genes in T. turgidum, which are classified into four main groups based on phylogenetic analysis. Several duplications were observed between the members of the TtAQP gene family, and high diversity in response to post-translational modifications was observed between TtAQP family members. The expression pattern of TtAQP genes disclosed that these genes are primarily upregulated in response to salt stress. Additionally, the qPCR data revealed that TtAQPs are more induced in delayed responses to salinity stress. Overall, our findings illustrate that TtAQP members are diverse in terms of their structure, regulatory systems, and expression levels.

Genome-wide analysis of the Tritipyrum WRKY gene family and the response of TtWRKY256 in salt-tolerance

INTRODUCTION

The transcription factor WRKY is widespread in the plant kingdom and plays a crucial role in diverse abiotic stress responses in plant species. Tritipyrum, an octoploid derived from an intergeneric cross between Triticum aestivum (AABBDD) and Thinopyrum elongatum (EE), is a valuable germplasm resource for introducing superior traits of Th. elongatum into T. aestivum. The recent release of the complete genome sequences of T. aestivum and Th. elongatum enabled us to investigate the organization and expression profiling of Tritipyrum WRKY genes across the entire genome.

RESULTS

In this study, 346 WRKY genes, from TtWRKY1 to TtWRKY346, were identified in Tritipyrum. The phylogenetic analysis grouped these genes into three subfamilies (I-III), and members of the same subfamilies shared a conserved motif composition. The 346 TtWRKY genes were dispersed unevenly across 28 chromosomes, with 218 duplicates. Analysis of synteny suggests that the WRKY gene family may have a common ancestor. Expression profiles derived from transcriptome data and qPCR demonstrated that 54 TtWRKY genes exhibited relatively high levels of expression across various salt stresses and recovery treatments. Tel1E01T143800 (TtWRKY256) is extremely sensitive to salt stress and is on the same evolutionary branch as the salt-tolerant A. thaliana genes AtWRKY25 and AtWRKY33. From 'Y1805', the novel AtWRKY25 was cloned. The Pearson correlation analysis identified 181 genes that were positively correlated (R>0.9) with the expression of TtWRKY256, and these genes were mainly enriched in metabolic processes, cellular processes, response to stimulus, biological regulation, and regulation of biological. Subcellular localization and qRT-PCR analysis revealed that TtWRKY256 was located in the nucleus and was highly expressed in roots, stems, and leaves under salt stress.

DISCUSSION

The above results suggest that TtWRKY256 may be associated with salt stress tolerance in plants and may be a valuable alien gene for improving salt tolerance in wheat.

Genome-wide identification and expression analysis reveals spinach brassinosteroid-signaling kinase (BSK) gene family functions in temperature stress response

Background

Brassinosteroid (BR)- signaling kinase (BSK) is a critical family of receptor-like cytoplasmic kinase for BR signal transduction, which plays important roles in plant development, immunity, and abiotic stress responses. Spinach (Spinacia oleracea) is cold- tolerant but heat- sensitive green leafy vegetable. A study on BSK family members and BSKs- mediated metabolic processes in spinach has not been performed.

Results

We identified and cloned seven SoBSKs in spinach. Phylogenetic and collinearity analyses suggested that SoBSKs had close relationship with dicotyledonous sugar beet (Beta vulgaris) rather than monocotyledons. The analyses of gene structure and conserved protein domain/ motif indicated that most SoBSKs were relative conserved, while SoBSK6 could be a truncated member. The prediction of post-translation modification (PTM) sites in SoBSKs implied their possible roles in signal transduction, redox regulation, and protein turnover of SoBSKs, especially the N-terminal myristoylation site was critical for BSK localization to cell periphery. Cis-acting elements for their responses to light, drought, temperature (heat and cold), and hormone distributed widely in the promoters of SoBSKs, implying the pivotal roles of SoBSKs in response to diverse abiotic stresses and phytohormone stimuli. Most SoBSKs were highly expressed in leaves, except for SoBSK7 in roots. Many SoBSKs were differentially regulated in spinach heat- sensitive variety Sp73 and heat- tolerant variety Sp75 under the treatments of heat, cold, as well as exogenous brassinolide (BL) and abscisic acid (ABA). The bsk134678 mutant Arabidopsis seedlings exhibited more heat tolerance than wild- type and SoBSK1- overexpressed seedlings.

Conclusions

A comprehensive genome- wide analysis of the BSK gene family in spinach presented a global identification and functional prediction of SoBSKs. Seven SoBSKs had relatively- conserved gene structure and protein function domains. Except for SoBSK6, all the other SoBSKs had similar motifs and conserved PTM sites. Most SoBSKs participated in the responses to heat, cold, BR, and ABA. These findings paved the way for further functional analysis on BSK- mediated regulatory mechanisms in spinach development and stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08684-5.

Genome-wide identification and expression profiling analysis of sucrose synthase (SUS) and sucrose phosphate synthase (SPS) genes family in Actinidia chinensis and A. eriantha

Sucrose synthase (SUS) is a common sugar-base transfer enzyme in plants, and sucrose phosphate synthase (SPS) is one of the major enzymes in higher plants that regulates sucrose synthesis. However, information of the SPS and SUS gene families in Actinidia, as well as their evolutionary and functional properties, is limited. According to the SPS and SUS proteins conserved domain of Arabidopsis thaliana, we found 6 SPS genes and 6 SUS genes from A. chinensis (cultivar: 'Hongyang'), and 3 SPS genes and 6 SUS genes from A. eriantha (cultivar: 'White'). The novel CDC50 conserved domains were discovered on AcSUS2, and all members of the gene family contain similar distinctive conserved domains. The majority of SUS and SPS proteins were hydrophilic, lipid-soluble enzymes that were expected to be found in the cytoplasm. The tertiary structure of SPS and SUS protein indicated that there were many tertiary structures in SPS, and there were windmill-type and spider-type tertiary structures in SUS. The phylogenetic tree was created using the neighbor-joining method, and members of the SPS and SUS gene families are grouped into three subgroups. Genes with comparable intron counts, conserved motifs, and phosphorylation sites were clustered together first. SPS and SUS were formed through replication among their own family members. AcSPS1, AcSPS2, AcSPS4, AcSPS5, AcSUS5, AcSUS6, AeSPS3, AeSUS3 and AeSUS4 were the important genes in regulating the synthesis and accumulation of sucrose for Actinidia during the fruit growth stages.

Investigation and Computational Analysis of the Sulfotransferase (SOT) Gene Family in Potato (Solanum tuberosum): Insights into Sulfur Adjustment for Proper Development and Stimuli Responses

Various kinds of primary metabolisms in plants are modulated through sulfate metabolism, and sulfotransferases (SOTs), which are engaged in sulfur metabolism, catalyze sulfonation reactions. In this study, a genome-wide approach was utilized for the recognition and characterization of SOT family genes in the significant nutritional crop potato (Solanum tuberosum L.). Twenty-nine putative StSOT genes were identified in the potato genome and were mapped onto the nine S. tuberosum chromosomes. The protein motifs structure revealed two highly conserved 5'-phosphosulfate-binding (5' PSB) regions and a 3'-phosphate-binding (3' PB) motif that are essential for sulfotransferase activities. The protein-protein interaction networks also revealed an interesting interaction between SOTs and other proteins, such as PRTase, APS-kinase, protein phosphatase, and APRs, involved in sulfur compound biosynthesis and the regulation of flavonoid and brassinosteroid metabolic processes. This suggests the importance of sulfotransferases for proper potato growth and development and stress responses. Notably, homology modeling of StSOT proteins and docking analysis of their ligand-binding sites revealed the presence of proline, glycine, serine, and lysine in their active sites. An expression essay of StSOT genes via potato RNA-Seq data suggested engagement of these gene family members in plants' growth and extension and responses to various hormones and biotic or abiotic stimuli. Our predictions may be informative for the functional characterization of the SOT genes in potato and other nutritional crops.

The BAHD Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome-Wide Identification and Expression Analysis

The benzyl alcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase, and deacetylvindoline 4-O-acetyltransferase (BAHD) enzymes play a critical role in regulating plant metabolites and affecting cell stability. In the present study, members of the BAHD gene family were recognized in the genome of Theobroma cacao and characterized using various bioinformatics tools. We found 27 non-redundant putative tcBAHD genes in cacao for the first time. Our findings indicate that tcBAHD genes are diverse based on sequence structure, physiochemical properties, and function. When analyzed with BAHDs of Gossypium raimondii and Corchorus capsularis clustered into four main groups. According to phylogenetic analysis, BAHD genes probably evolved drastically after their divergence. The divergence time of duplication events with purifying selection pressure was predicted to range from 1.82 to 15.50 MYA. Pocket analysis revealed that serine amino acid is more common in the binding site than other residuals, reflecting its key role in regulating the activity of tcBAHDs. Furthermore, cis-acting elements related to the responsiveness of stress and hormone, particularly ABA and MeJA, were frequently observed in the promoter region of tcBAHD genes. RNA-seq analysis further illustrated that tcBAHD13 and tcBAHD26 are involved in response to Phytophthora megakarya fungi. In conclusion, it is likely that evolutionary processes, such as duplication events, have caused high diversity in the structure and function of tcBAHD genes.

The GASA Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome Wide Identification and Expression Analysis

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants and involved in various physiological and biological processes. These genes also provide resistance to abiotic and biotic stresses, including antimicrobial, antiviral, and antifungal. We are interested in characterizing the GASA gene family and determining its role in various physiological and biological process in Theobroma cacao. Here, we report 17 tcGASA genes distributed on six chromosomes in T. cacao. The gene structure, promoter region, protein structure and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. Among 17 tcGASA genes, nine segmentally duplicating genes were identified which formed four pairs and cluster together in phylogenetic tree. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The differential expression of tcGASAs was recorded between the tolerant and susceptible cultivars of cacao, which indicating their possible role as fungal resistant. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T.cacao and may suggest new target genes for development of fungi-resistant cacao varieties in breeding programs.

Genome-Wide Analysis of Potassium Channel Genes in Rice: Expression of the OsAKT and OsKAT Genes under Salt Stress

Potassium (K+), as a vital element, is involved in regulating important cellular processes such as enzyme activity, cell turgor, and nutrient movement in plant cells, which affects plant growth and production. Potassium channels are involved in the transport and release of potassium in plant cells. In the current study, three OsKAT genes and two OsAKT genes, along with 11 nonredundant putative potassium channel genes in the rice genome, were characterized based on their physiochemical properties, protein structure, evolution, duplication, in silico gene expression, and protein-protein interactions. In addition, the expression patterns of OsAKTs and OsKATs were studied in root and shoot tissues under salt stress using real-time PCR in three rice cultivars. K+ channel genes were found to have diverse functions and structures, and OsKATs showed high genetic divergence from other K+ channel genes. Furthermore, the Ka/Ks ratios of duplicated gene pairs from the K+ channel gene family in rice suggested that these genes underwent purifying selection. Among the studied K+ channel proteins, OsKAT1 and OsAKT1 were identified as proteins with high potential N-glycosylation and phosphorylation sites, and LEU, VAL, SER, PRO, HIS, GLY, LYS, TYR, CYC, and ARG amino acids were predicted as the binding residues in the ligand-binding sites of K+ channel proteins. Regarding the coexpression network and KEGG ontology results, several metabolic pathways, including sugar metabolism, purine metabolism, carbon metabolism, glycerophospholipid metabolism, monoterpenoid biosynthesis, and folate biosynthesis, were recognized in the coexpression network of K+ channel proteins. Based on the available RNA-seq data, the K+ channel genes showed differential expression levels in rice tissues in response to biotic and abiotic stresses. In addition, the real-time PCR results revealed that OsAKTs and OsKATs are induced by salt stress in root and shoot tissues of rice cultivars, and OsKAT1 was identified as a key gene involved in the rice response to salt stress. In the present study, we found that the repression of OsAKTs, OsKAT2, and OsKAT2 in roots was related to salinity tolerance in rice. Our findings provide valuable insights for further structural and functional assays of K+ channel genes in rice.

New Insights Into Structure and Function of TIFY Genes in Zea mays and Solanum lycopersicum: A Genome-Wide Comprehensive Analysis

The TIFY gene family, a key plant-specific transcription factor (TF) family, is involved in diverse biological processes including plant defense and growth regulation. Despite TIFY proteins being reported in some plant species, a genome-wide comparative and comprehensive analysis of TIFY genes in plant species can reveal more details. In the current study, the members of the TIFY gene family were significantly increased by the identification of 18 and six new members using maize and tomato reference genomes, respectively. Thus, a genome-wide comparative analysis of the TIFY gene family between 48 tomato (Solanum lycopersicum, a dicot plant) genes and 26 maize (Zea mays, a monocot plant) genes was performed in terms of sequence structure, phylogenetics, expression, regulatory systems, and protein interaction. The identified TIFYs were clustered into four subfamilies, namely, TIFY-S, JAZ, ZML, and PPD. The PPD subfamily was only detected in tomato. Within the context of the biological process, TIFY family genes in both studied plant species are predicted to be involved in various important processes, such as reproduction, metabolic processes, responses to stresses, and cell signaling. The Ka/Ks ratios of the duplicated paralogous gene pairs indicate that all of the duplicated pairs in the TIFY gene family of tomato have been influenced by an intense purifying selection, whereas in the maize genome, there are three duplicated blocks containing Ka/Ks > 1, which are implicated in evolution with positive selection. The amino acid residues present in the active site pocket of TIFY proteins partially differ in each subfamily, although the Mg or Ca ions exist heterogeneously in the centers of the active sites of all the predicted TIFY protein models. Based on the expression profiles of TIFY genes in both plant species, JAZ subfamily proteins are more associated with the response to abiotic and biotic stresses than other subfamilies. In conclusion, globally scrutinizing and comparing the maize and tomato TIFY genes showed that TIFY genes play a critical role in cell reproduction, plant growth, and responses to stress conditions, and the conserved regulatory mechanisms may control their expression.

Exogenous EBR Ameliorates Endogenous Hormone Contents in Tomato Species under Low-Temperature Stress

Low-temperature stress is a type of abiotic stress that limits plant growth and production in both subtropical and tropical climate conditions. In the current study, the effects of 24-epi-brassinolide (EBR) as analogs of brassinosteroids (BRs) were investigated, in terms of hormone content, antioxidant enzyme activity, and transcription of several cold-responsive genes, under low-temperature stress (9 °C) in two different tomato species (cold-sensitive and cold-tolerant species). Results indicated that the treatment with exogenous EBR increases the content of gibberellic acid (GA3) and indole-3-acetic acid (IAA), whose accumulation is reduced by low temperatures in cold-sensitive species. Furthermore, the combination or contribution of BR and abscisic acid (ABA) as a synergetic interaction was recognized between BR and ABA in response to low temperatures. The content of malondialdehyde (MDA) and proline was significantly increased in both species, in response to low-temperature stress; however, EBR treatment did not affect the MDA and proline content. Moreover, in the present study, the effect of EBR application was different in the tomato species under low-temperature stress, which increased the catalase (CAT) activity in the cold-tolerant species and increased the glutathione peroxidase (GPX) activity in the cold-sensitive species. Furthermore, expression levels of cold-responsive genes were influenced by low-temperature stress and EBR treatment. Overall, our findings revealed that a low temperature causes oxidative stress while EBR treatment may decrease the reactive oxygen species (ROS) damage into increasing antioxidant enzymes, and improve the growth rate of the tomato by affecting auxin and gibberellin content. This study provides insight into the mechanism by which BRs regulate stress-dependent processes in tomatoes, and provides a theoretical basis for promoting cold resistance of the tomato.

Genome-wide comparative analysis of Mg transporter gene family between Triticum turgidum and Camelina sativa

Magnesium (Mg) as a bimetal plays critical roles in biochemical processes, membrane stability, and enzyme activity. Mg transporters (MGTs) are involving in maintaining Mg homeostasis in cells. Although the MGT family members have been identified in different plant species, there is no comprehensive analysis of the other plants' MGT genes. In the current study, 62 and 41 non-redundant putative MGT proteins were recognized into the genome of Camelina sativa, and Triticum turgidum and they were compared based on physicochemical properties, protein structure, expression, and interaction. All identified MGTs were classified into three subgroups, NIPA, CorA, and MRS2/MGT, based on conserved-motifs distribution. The results showed that the secondary structure pattern in NIPA and MRS2 subfamily members in both studied plant species were highly similar. Furthermore, MGTs encompass the conserved structures and the critical sites mainly in the metal ion and Mg²⁺ binding centers as well as the catalytic sites were observed. The highest numbers of protein channels were predicted in CorA proteins in both C. sativa and T. turgidum with 24 and 17 channel numbers, respectively. The Ser, Pro, Gly, Lys, Tyr, and Arg amino acids were predicted as the binding residues in MGTs channel regions. The expression pattern of identified genes demonstrated that MGT genes have diverse tissue-specific expression and stress response expression patterns. Besides, 147 co-expressed genes with MGTs were clustered into the eight co-expression nodes involved in N-glycan biosynthesis, protein processing in the endoplasmic reticulum, carbon metabolism, biosynthesis of amino acids, and endocytosis. In the present study, all interpretations are based on in silico predictions, which can be used in further studies related to functional genomics of MGT genes.

The GASA Gene Family in Theobroma cacao: Genome wide Identification and Expression Analysis.. [DOI: 10.1101/2021.01.27.425041]

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants. The role of the GASA gene family has been reported previously in various physiological and biological processes, such as cell division, root and seed development, stem growth, and fruit ripening. These genes also provide resistance to abiotic and biotic stresses including antimicrobial, antiviral, and antifungal. Here, we report 17 tcGASA genes in Theobroma cacao L. distributed on six chromosomes. The gene structure, promoter-region sequences, protein structure, and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. The nine segmentally duplicating genes form four pairs and cluster together in phylogenetic tree. Purifying selection pressure was recorded on tcGASA, including duplicated genes. Several stress/hormone-responsive cis-regulatory elements were also recognized in the promoter region of tcGASAs. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The black rod disease of genus Phytophthora caused up to 20–25% loss (700,000 metric tons) in world cacao production. The role of tcGASA genes in conferring fungal resistance was also explored based on RNAseq data against Phytophthora megakarya. The differential expression of tcGASA genes was recorded between the tolerant and susceptible cultivars of cacao plants, which were inoculated with the fungus for 24h and 72h. This differential expression indicating possible role of tcGASA genes to fungal resistant in cacao. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T. cacao and provide new target genes for development of fungi-resistant cacao varieties in breeding programs.

The AP2/ERF Gene Family in Triticum durum: Genome-Wide Identification and Expression Analysis under Drought and Salinity Stresses

Members of the AP2/ERF transcription factor family play critical roles in plant development, biosynthesis of key metabolites, and stress response. A detailed study was performed to identify TtAP2s/ERFs in the durum wheat (Triticum turgidum ssp. durum) genome, which resulted in the identification of 271 genes distributed on chromosomes 1A-7B. By carrying 27 genes, chromosome 6A had the highest number of TtAP2s/ERFs. Furthermore, a duplication assay of TtAP2s/ERFs demonstrated that 70 duplicated gene pairs had undergone purifying selection. According to RNA-seq analysis, the highest expression levels in all tissues and in response to stimuli were associated with DRF and ERF subfamily genes. In addition, the results revealed that TtAP2/ERF genes have tissue-specific expression patterns, and most TtAP2/ERF genes were significantly induced in the root tissue. Additionally, 13 TtAP2/ERF genes (six ERFs, three DREBs, two DRFs, one AP2, and one RAV) were selected for further analysis via qRT-PCR of their potential in coping with drought and salinity stresses. The TtAP2/ERF genes belonging to the DREB subfamily were markedly induced under both drought-stress and salinity-stress conditions. Furthermore, docking simulations revealed several residues in the pocket sites of the proteins associated with the stress response, which may be useful in future site-directed mutagenesis studies to increase the stress tolerance of durum wheat. This study could provide valuable insights for further evolutionary and functional assays of this important gene family in durum wheat.

Insights into the genes involved in the ethylene biosynthesis pathway in Arabidopsis thaliana and Oryza sativa

Background

Ethylene is a gaseous plant hormone that acts as a requisite role in many aspects of the plant life cycle, and it is also a regulator of plant responses to abiotic and biotic stresses. In this study, we attempt to provide comprehensive information through analyses of existing data using bioinformatics tools to compare the identified ethylene biosynthesis genes between Arabidopsis (as dicotyledonous) and rice (as monocotyledonous).

Results

The results exposed that the Arabidopsis proteins of the ethylene biosynthesis pathway had more potential glycosylation sites than rice, and 1-aminocyclopropane-1-carboxylate oxidase proteins were less phosphorylated than 1-aminocyclopropane-1-carboxylate synthase and S-adenosylmethionine proteins. According to the gene expression patterns, S-adenosylmethionine genes were more involved in the rice-ripening stage while in Arabidopsis, ACS2, and 1-aminocyclopropane-1-carboxylate oxidase genes were contributed to seed maturity. Furthermore, the result of miRNA targeting the transcript sequences showed that ath-miR843 and osa-miR1858 play a key role to regulate the post-transcription modification of S-adenosylmethionine genes in Arabidopsis and rice, respectively. The discovered cis- motifs in the promoter site of all the ethylene biosynthesis genes of A. thaliana genes were engaged to light-induced response in the cotyledon and root genes, sulfur-responsive element, dehydration, cell cycle phase-independent activation, and salicylic acid. The ACS4 protein prediction demonstrated strong protein-protein interaction in Arabidopsis, as well as, SAM2, Os04T0578000, Os01T0192900, and Os03T0727600 predicted strong protein-protein interactions in rice.

Conclusion

In the current study, the complex between miRNAs with transcript sequences of ethylene biosynthesis genes in A. thaliana and O. sativa were identified, which could be helpful to understand the gene expression regulation after the transcription process. The binding sites of common transcription factors such as MYB, WRKY, and ABRE that control target genes in abiotic and biotic stresses were generally distributed in promoter sites of ethylene biosynthesis genes of A. thaliana. This was the first time to wide explore the ethylene biosynthesis pathway using bioinformatics tools that markedly showed the capability of the in silico study to integrate existing data and knowledge and furnish novel insights into the understanding of underlying ethylene biosynthesis pathway genes that will be helpful for more dissection.