Genome-Wide Characterization of GRAS Family and Their Potential Roles in Cold Tolerance of Cucumber (Cucumis sativus L.)

Genome-wide identification and expression analysis of GRAS gene family in Eucalyptus grandis

BACKGROUND

The GRAS gene family is a class of plant-specific transcription factors with important roles in many biological processes, such as signal transduction, disease resistance and stress tolerance, plant growth and development. So far, no information available describes the functions of the GRAS genes in Eucalyptus grandis.

RESULTS

A total of 82 GRAS genes were identified with amino acid lengths ranging from 267 to 817 aa, and most EgrGRAS genes had one exon. Members of the GRAS gene family of Eucalyptus grandis are divided into 9 subfamilies with different protein structures, while members of the same subfamily have similar gene structures and conserved motifs. Moreover, these EgrGRAS genes expanded primarily due to segmental duplication. In addition, cis-acting element analysis showed that this family of genes was involved involved in the signal transduction of various plant hormones, growth and development, and stress response. The qRT-PCR data indicated that 18 EgrGRAS genes significantly responded to hormonal and abiotic stresses. Among them, the expression of EgrGRAS13, EgrGRAS68 and EgrGRAS55 genes was significantly up-regulated during the treatment period, and it was hypothesised that members of the EgrGRAS family play an important role in stress tolerance.

CONCLUSIONS

In this study, the phylogenetic relationship, conserved domains, cis-elements and expression patterns of GRAS gene family of Eucalyptus grandis were analyzed, which filled the gap in the identification of GRAS gene family of Eucalyptus grandis and laid the foundation for analyzing the function of EgrGRAS gene in hormone and stress response.

Genome-wide identification and expression pattern analysis of the kiwifruit GRAS transcription factor family in response to salt stress

BACKGROUND

GRAS is a family of plant-specific transcription factors (TFs) that play a vital role in plant growth and development and response to adversity stress. However, systematic studies of the GRAS TF family in kiwifruit have not been reported.

RESULTS

In this study, we used a bioinformatics approach to identify eighty-six AcGRAS TFs located on twenty-six chromosomes and phylogenetic analysis classified them into ten subfamilies. It was found that the gene structure is relatively conserved for these genes and that fragmental duplication is the prime force for the evolution of AcGRAS genes. However, the promoter region of the AcGRAS genes mainly contains cis-acting elements related to hormones and environmental stresses, similar to the results of GO and KEGG enrichment analysis, suggesting that hormone signaling pathways of the AcGRAS family play a vital role in regulating plant growth and development and adversity stress. Protein interaction network analysis showed that the AcGRAS51 protein is a relational protein linking DELLA, SCR, and SHR subfamily proteins. The results demonstrated that 81 genes were expressed in kiwifruit AcGRAS under salt stress, including 17 differentially expressed genes, 13 upregulated, and four downregulated. This indicates that the upregulated AcGRAS55, AcGRAS69, AcGRAS86 and other GRAS genes can reduce the salt damage caused by kiwifruit plants by positively regulating salt stress, thus improving the salt tolerance of the plants.

CONCLUSIONS

These results provide a theoretical basis for future exploration of the characteristics and functions of more AcGRAS genes. This study provides a basis for further research on kiwifruit breeding for resistance to salt stress. RT-qPCR analysis showed that the expression of 3 AcGRAS genes was elevated under salt stress, indicating that AcGRAS exhibited a specific expression pattern under salt stress conditions.

A genome-wide analysis of the GRAS gene family in upland cotton and a functional study of the role of the GhGRAS55 gene in regulating early maturity in cotton

The members of the GRAS gene family play important roles in regulating plant growth and development, but their functions in regulating early plant maturity traits are still unknown. In this study, we used a series of bioinformatics tools to identify GRAS gene family members and investigate the function of the gene family (GhGRAS55) using a genome-wide database of upland cotton samples. A total of 58 members of the GRAS gene family were identified and screened, which were distributed on 21 chromosomes within the whole cotton genome. The results of the phylogenetic analysis showed that the genes of upland cotton, island cotton, African cotton, Raymond cotton, and Arabidopsis were distributed in subfamilies I-VIII, although subfamily II did not contain any upland cotton or Arabidopsis GRAS family members. The structures and other characteristics of the genes in this family were clarified using bioinformatics technology. The transcriptomic sequencing results for early and late maturing cotton species showed that the expression of most GRAS family genes, such as GhGRAS10, GhGRAS5511, and GhGRAS55, was lower in early maturing species than late maturing species. We also found that cotton plants with GhGRAS55 genes that were silenced by virus-induced gene silencing (VIGS) technology showed early bud emergence phenotypes, so it could be speculated that the GhGRAS55 gene has the function of regulating early maturity in cotton.

The role of shoot-derived RNAs transported to plant root in response to abiotic stresses

A large number of RNA molecules are transported over long-distance between shoots and roots via phloem in higher plants. Mobile RNA signals are important for plants to tackle abiotic stresses. Shoot-derived mobile RNAs can be involved in the response to different developmental or environmental signals in the root. Some environmental conditions such as climate change, water deficit, nutrient deficiency challenge modern agriculture with more expeditious abiotic stress conditions. Root architecture determines the ability of water and nutrient uptake and further abiotic stress tolerance, and shoot tissue also determines the balance between shoot-root relationship in plant growth and adaptations. Thus, it is necessary to understand the roles of shoot-derived RNA signals and their potential function in roots upon abiotic stresses in the model plants (Arabidopsis thaliana and Nicotiana benthamiana) and agricultural crops. In this review, we summarize the so-far discovered shoot-derived mobile RNA transportation to the root under abiotic stress conditions, e.g. drought, cold stress and nutrient deficiencies. Furthermore, we will focus on the biological relevance and the potential roles of these RNAs in root development and stress responses which will be an asset for the future breeding strategies.

Genome-wide identification and characterization of abiotic stress responsive GRAS family genes in oat (Avena sativa)

Background

GRAS transcription factors play a variety of functions in plant growth and development and are named after the first three transcription factors GAI (GIBBERRELLICACIDINSENSITIVE), RGA (REPRESSOROFGAI), and SCR (SCARECROW) found in this family. Oat (Avena sativa) is one of the most important forage grasses in the world. However, there are few reports on the GRAS gene family in oat.

Methods

In order to understand the information and expression pattern of oat GRAS family members, we identified the GRAS members and analyzed their phylogenetic relationship, gene structure, and expression pattern in oat by bioinformatics technology.

Results

The results showed that the oat GRAS family consists of 30 members, and most of the AsGRAS proteins were neutral or acidic proteins. The phylogenetic tree divided the oat GRAS members into four subfamilies, and each subfamily has different conservative domains and functions. Chromosome location analysis suggested that 30 GRAS genes were unevenly distributed on five chromosomes of oat. The results of real-time quantitative reverse transcription-PCR (qRT-PCR) showed that some AsGRAS genes (AsGRAS12, AsGRAS14, AsGRAS21, and AsGRAS24) were all up-regulated with increasing stress treatment time.The results of this study provide a theoretical basis for further research into the corresponding stress of oat. Therefore, further studies concentrating on these AsGRAS genes might reveal the many roles played by GRAS genes in oat.

Evolution and functional analysis of the GRAS family genes in six Rosaceae species

BACKGROUND

GRAS genes formed one of the important transcription factor gene families in plants, had been identified in several plant species. The family genes were involved in plant growth, development, and stress resistance. However, the comparative analysis of GRAS genes in Rosaceae species was insufficient.

RESULTS

In this study, a total of 333 GRAS genes were identified in six Rosaceae species, including 51 in strawberry (Fragaria vesca), 78 in apple (Malus domestica), 41 in black raspberry (Rubus occidentalis), 59 in European pear (Pyrus communis), 56 in Chinese rose (Rosa chinensis), and 48 in peach (Prunus persica). Motif analysis showed the VHIID domain, SAW motif, LR I region, and PFYRE motif were considerably conserved in the six Rosaceae species. All GRAS genes were divided into 10 subgroups according to phylogenetic analysis. A total of 15 species-specific duplicated clades and 3 lineage-specific duplicated clades were identified in six Rosaceae species. Chromosomal localization presented the uneven distribution of GRAS genes in six Rosaceae species. Duplication events contributed to the expression of the GRAS genes, and Ka/Ks analysis suggested the purification selection as a major force during the evolution process in six Rosaceae species. Cis-acting elements and GO analysis revealed that most of the GRAS genes were associated with various environmental stress in six Rosaceae species. Coexpression network analysis showed the mutual regulatory relationship between GRAS and bZIP genes, suggesting the ability of the GRAS gene to regulate abiotic stress in woodland strawberry. The expression pattern elucidated the transcriptional levels of FvGRAS genes in various tissues and the drought and salt stress in woodland strawberry, which were verified by RT-qPCR analysis.

CONCLUSIONS

The evolution and functional analysis of GRAS genes provided insights into the further understanding of GRAS genes on the abiotic stress of Rosaceae species.

Genome-wide identification of RsGRAS gene family reveals positive role of RsSHRc gene in chilling stress response in radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an important worldwide root vegetable crop. Little information of the GRAS gene family was available in radish. Herein, a total of 51 GRAS family members were firstly identified from radish genome, and unevenly located onto nine radish chromosomes. Expression analysis of RsGRAS genes in taproot displayed that RsSCL15a and RsSHRc were highly expressed in the radish cambium, and its expression level was increased with the taproot thickening. Comparative transcriptome analysis revealed that the expression patterns of RsGRAS genes varied upon exposure to different abiotic stresses including heavy metals, salt and heat. The expression level of six RsGRAS genes including RsSHRc was increased under chilling stress in two radish genotypes with different cold tolerance. Further analysis indicated that RsGRAS genes could respond to cold stress rapidly and the expression of RsSHRc was up-regulated at different development stages (cortex splitting and thickening stages) under long-term cold treatment. Transient expression of RsSHRc gene in radish showed that RsSHRc possessed the reliable function of eliminating reactive oxygen species (ROS), inhibiting the formation of malondialdehyde (MDA) and promoting to accumulate proline under cold stress. Together, these findings provided insights into the function of RsGRAS genes in taproot development and chilling stress response in radish.

Integrative Transcriptomic and Phytohormonal Analyses Provide Insights into the Cold Injury Recovery Mechanisms of Tea Leaves

Tea plant is susceptible to low temperature, while the cold injury recovery mechanisms of tea leaves are still unclear. Windbreak has an effective and gradient range of protecting tea plants. Tea plants with increasing cold damage degree have varying recovery status accordingly, which are the ideal objects for investigating the cold injury recovery mechanisms of tea leaves. Here, we investigated the transcriptome and phytohormone profiles of tea leaves with different cold injury degrees in recovery (adjacent to the windbreak), and the levels of chlorophylls, malondialdehyde, major phytohormones as well as the activities of peroxidase (POD) and superoxide dismutase (SOD) were also measured. The results showed the content of total chlorophylls and the activity of POD in mature tea leaves gradually decreased with the distance to windbreak, while SOD showed the opposite. The major phytohormones were highly accumulated in the moderately cold-injured tea leaves. The biosynthesis of abscisic acid (ABA) was enhanced in the moderate cold damaged tea leaves, suggesting that ABA plays an important role in the cold response and resistance of tea plants. The transcriptomic result showed that the samples in different rows were well discriminated, and the pathways of plant-pathogen interaction and flavonoid biosynthesis were enriched based on KEGG analysis. WRKY, GRAS and NAC were the top classes of transcription factors differentially expressed in the different cold-injured tea leaves. Thus, windbreak is effective to protect adjacent tea plants from cold wave, and phytohormones importantly participate in the cold injury recovery of tea leaves.

MdSCL8 as a Negative Regulator Participates in ALA-Induced FLS1 to Promote Flavonol Accumulation in Apples

Apples (Malus domestica) are rich in flavonols, and 5-aminolevulinic acid (ALA) plays an important role in the regulation of plant flavonoid metabolism. To date, the underlying mechanism of ALA promoting flavonol accumulation is unclear. Flavonol synthase (FLS) is a key enzyme in flavonol biosynthesis. In this study, we found that ALA could enhance the promoter activity of MdFLS1 in the ‘Fuji’ apple and improve its expression. With MdFLS1 as bait, we screened a novel transcription factor MdSCL8 by the Yeast One-Hybrid (Y1H) system from the apple cDNA library which we previously constructed. Using luciferase reporter assay and transient GUS activity assay, we verified that MdSCL8 inhibits the activity of MdFLS1 promoter and hinders MdFLS1 expression, thus reducing flavonol accumulation in apple. ALA significantly inhibited MdSCL8 expression. Therefore, ALA promoted the expression of MdFLS1 and the consequent flavonol accumulation probably by down-regulating MdSCL8. We also found that ALA significantly enhanced the gene expression of MdMYB22 and MdHY5, two positive regulators of MdFLS. We further demonstrated that MdMYB22 interacts with MdHY5, but neither of them interacts with MdSCL8. Taken together, our data suggest MdSCL8 as a novel regulator of MdFLS1 and provide important insights into mechanisms of ALA-induced flavonol accumulation in apples.

Genome-wide investigation of the GRAS transcription factor family in foxtail millet (Setaria italica L.)

BACKGROUND

GRAS transcription factors perform indispensable functions in various biological processes, such as plant growth, fruit development, and biotic and abiotic stress responses. The development of whole-genome sequencing has allowed the GRAS gene family to be identified and characterized in many species. However, thorough in-depth identification or systematic analysis of GRAS family genes in foxtail millet has not been conducted.

RESULTS

In this study, 57 GRAS genes of foxtail millet (SiGRASs) were identified and renamed according to the chromosomal distribution of the SiGRAS genes. Based on the number of conserved domains and gene structure, the SiGRAS genes were divided into 13 subfamilies via phylogenetic tree analysis. The GRAS genes were unevenly distributed on nine chromosomes, and members of the same subfamily had similar gene structures and motif compositions. Genetic structure analysis showed that most SiGRAS genes lacked introns. Some SiGRAS genes were derived from gene duplication events, and segmental duplications may have contributed more to GRAS gene family expansion than tandem duplications. Quantitative polymerase chain reaction showed significant differences in the expression of SiGRAS genes in different tissues and stages of fruits development, which indicated the complexity of the physiological functions of SiGRAS. In addition, exogenous paclobutrazol treatment significantly altered the transcription levels of DELLA subfamily members, downregulated the gibberellin content, and decreased the plant height of foxtail millet, while it increased the fruit weight. In addition, SiGRAS13 and SiGRAS25 may have the potential for genetic improvement and functional gene research in foxtail millet.

CONCLUSIONS

Collectively, this study will be helpful for further analysing the biological function of SiGRAS. Our results may contribute to improving the genetic breeding of foxtail millet.

Genome-Wide Identification of GRAS Gene Family and Their Responses to Abiotic Stress in Medicago sativa

Alfalfa (Medicago sativa) is a high-quality legume forage crop worldwide, and alfalfa production is often threatened by abiotic environmental stresses. GRAS proteins are important transcription factors that play a vital role in plant development, as well as in response to environmental stress. In this study, the availability of alfalfa genome "Zhongmu No.1" allowed us to identify 51 GRAS family members, i.e., MsGRAS. MsGRAS proteins could be classified into nine subgroups with distinct conserved domains, and tandem and segmental duplications were observed as an expansion strategy of this gene family. In RNA-Seq analysis, 14 MsGRAS genes were not expressed in the leaf or root, 6 GRAS genes in 3 differentially expressed gene clusters were involved in the salinity stress response in the leaf. Moreover, qRT-PCR results confirmed that MsGRAS51 expression was induced under drought stress and hormone treatments (ABA, GA and IAA) but down-regulated in salinity stress. Collectively, our genome-wide characterization, evolutionary, and expression analysis suggested that the MsGRAS proteins might play crucial roles in response to abiotic stresses and hormonal cues in alfalfa. For the breeding of alfalfa, it provided important information on stress resistance and functional studies on MsGRAS and hormone signaling.

Genome-wide identification, expression analysis, and functional study of the GRAS transcription factor family and its response to abiotic stress in sorghum [Sorghum bicolor (L.) Moench]

Background

GRAS, an important family of transcription factors, have played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. Since the sequencing of the sorghum genome, a plethora of genetic studies were mainly focused on the genomic information. The indepth identification or genome-wide analysis of GRAS family genes, especially in Sorghum bicolor, have rarely been studied.

Results

A total of 81 SbGRAS genes were identified based on the S. bicolor genome. They were named SbGRAS01 to SbGRAS81 and grouped into 13 subfamilies (LISCL, DLT, OS19, SCL4/7, PAT1, SHR, SCL3, HAM-1, SCR, DELLA, HAM-2, LAS and OS4). SbGRAS genes are not evenly distributed on the chromosomes. According to the results of the gene and motif composition, SbGRAS members located in the same group contained analogous intron/exon and motif organizations. We found that the contribution of tandem repeats to the increase in sorghum GRAS members was slightly greater than that of fragment repeats. By quantitative (q) RT-PCR, the expression of 13 SbGRAS members in different plant tissues and in plants exposed to six abiotic stresses at the seedling stage were quantified. We further investigated the relationship between DELLA genes, GAs and grain development in S. bicolor. The paclobutrazol treatment significantly increased grain weight, and affected the expression levels of all DELLA subfamily genes. SbGRAS03 is the most sensitive to paclobutrazol treatment, but also has a high response to abiotic stresses.

Conclusions

Collectively, SbGRAs play an important role in plant development and response to abiotic stress. This systematic analysis lays the foundation for further study of the functional characteristics of GRAS genes of S. bicolor.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07848-z.

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

Genome-Wide Characterization of Cucumber (Cucumis sativus L.) GRAS Genes and Their Response to Various Abiotic Stresses

The GRAS (gibberellic acid insensitive, repressor of GAI, and scarecrow) proteins are a family of plant-specific transcription factors that regulate plant growth, development, and stress response. Currently, the role of GRAS transcription factors in various abiotic stress responses has not been systematically studied in cucumber (Cucumis sativus L.), a popular vegetable crop. Here, we provide a comprehensive bioinformatics analysis of the 35 GRAS genes identified in the cucumber genome. In this study, cucumber genotypes, i.e., “CG104”, which is stress-tolerant, and genotype “CG37”, which is stress-sensitive, were examined to provide insight on potential differences in the GRAS-regulated abiotic stress pathways. Transcriptional analysis by RNA-seq or qRT-PCR of these two genotypes revealed common and divergent functions of CsGRAS genes regulated by low and high temperatures, salinity, and by exposure to the phytohormones gibberellin (GA) and abscisic acid (ABA). Notably, CsGRAS2 (DELLA) and CsGRAS26 (LISCL) were regulated by all abiotic stresses and hormone treatments, suggesting that they may function in the biological cross-talk between multiple signaling pathways. This study provides candidate genes for improving cucumber tolerance to various environmental stresses.