Genome-wide identification and expression analysis of GRAS family transcription factors in tea plant (Camellia sinensis)

Identification of candidate single-nucleotide polymorphisms (SNPs) and genes associated with sugarcane leaf scald disease

Leaf scald, caused by Xanthomonas albilineans, is a severe disease affecting sugarcane worldwide. One of the most practical ways to control it is by developing resistant sugarcane cultivars. It is essential to identify genes associated with the response to leaf scald. A panel of 170 sugarcane genotypes was evaluated for resistance to leaf scald in field conditions for 2 years, followed by a 1-year greenhouse experiment. The phenotypic evaluation data showed a wide continuous distribution, with heritability values ranging from 0.58 to 0.84. Thirteen single nucleotide polymorphisms (SNPs) were identified, significantly associated with leaf scald resistance. Among these, eight were stable across multiple environments and association models. The candidate genes identified and validated based on RNA-seq and qRT-PCR included two genes that encode NB-ARC leucine-rich repeat (LRR)-containing domain disease-resistance protein. These findings provide a basis for developing marker-assisted selection strategies in sugarcane breeding programs.

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.

CsPAT1, a GRAS transcription factor, promotes lignin accumulation by antagonistic interacting with CsWRKY13 in tea plants

Lignin is an important component of plant cell walls and plays crucial roles in the essential agronomic traits of tea quality and tenderness. However, the molecular mechanisms underlying the regulation of lignin biosynthesis in tea plants remain unclear. CsWRKY13 acts as a negative regulator of lignin biosynthesis in tea plants. In this study, we identified a GRAS transcription factor, phytochrome A signal transduction 1 (CsPAT1), that interacts with CsWRKY13. Silencing CsPAT1 expression in tea plants and heterologous overexpression in Arabidopsis demonstrated that CsPAT1 positively regulates lignin accumulation. Further investigation revealed that CsWRKY13 directly binds to the promoters of CsPAL and CsC4H and suppresses transcription of CsPAL and CsC4H. CsPAT1 indirectly affects the promoter activities of CsPAL and CsC4H by interacting with CsWRKY13, thereby facilitating lignin biosynthesis in tea plants. Compared with the expression of CsWRKY13 alone, the co-expression of CsPAT1 and CsWRKY13 in Oryza sativa significantly increased lignin biosynthesis. Conversely, compared with the expression of CsPAT1 alone, the co-expression of CsPAT1 and CsWRKY13 in O. sativa significantly reduced lignin accumulation. These results demonstrated the antagonistic regulation of the lignin biosynthesis pathway by CsPAT1 and CsWRKY13. These findings improve our understanding of lignin biosynthesis mechanisms in tea plants and provide insights into the role of the GRAS transcription factor family in lignin accumulation.

Comprehensive Analysis of the GRAS Gene Family in Paulownia fortunei and the Response of DELLA Proteins to Paulownia Witches' Broom

The GRAS (GAI\RGA\SCL) gene family encodes plant-specific transcription factors that play crucial roles in plant growth and development, stress tolerance, and hormone network regulation. Plant dwarfing symptom is mainly regulated by DELLA proteins of the GRAS gene subfamily. In this study, the association between the GRAS gene family and Paulownia witches' broom (PaWB) was investigated. A total of 79 PfGRAS genes were identified using bioinformatics methods and categorized into 11 groups based on amino acid sequences. Tandem duplication and fragment duplication were found to be the main modes of amplification of the PfGRAS gene family. Gene structure analysis showed that more than 72.1% of the PfGRASs had no introns. The genes PfGRAS12/18/58 also contained unique DELLA structural domains; only PfGRAS12, which showed significant response to PaWB phytoplasma infection in stems, showed significant tissue specificity and responded to gibberellin (GA3) in PaWB-infected plants. We found that the internodes were significantly elongated under 100 µmol·L-1 GA3 treatment for 30 days. The subcellular localization analysis indicated that PfGRAS12 is located in the nucleus and cell membrane. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays confirmed that PfGRAS12 interacted with PfJAZ3 in the nucleus. Our results will lay a foundation for further research on the functions of the PfGRAS gene family and for genetic improvement and breeding of PaWB-resistant trees.

Genome-wide identification and characterization of GRAS gene family in pigeonpea (Cajanus cajan (L.) Millspaugh)

The GRAS proteins are plant-specific transcription factors (TFs) that play a crucial role in various plant physiological processes, including tissue development and stress responses. To date, GRAS family has been comprehensively characterized in Arabidopsis, soybean, rice, chickpea and other plant species. To understand the structural and functional aspects of pigeonpea (C. cajan), we identified 60 putative GRAS (CcGRAS) genes from pigeonpea genome and further analysed their physicochemical properties, subcellular locations, evolutionary classification, exon-intron structures, conserved domains, gene duplication events and cis-promoter regions. Based on the sequence similarity, CcGRAS family was clustered into 9 subfamilies and the genes with a similar structure and motif distribution were clustered in the same group. The gene duplication studies revealed that these genes were derived from tandem and dispersed duplication events. The cis-promoter regulatory analysis of CcGRAS genes indicated the presence of three types of cis-acting elements including light-responsive, hormone-responsive and plant growth and development related. The expression profiling of CcGRAS genes revealed their tissue-specific functions and differential nature. Collectively, this study highlights relevant functional and regulatory elements of GRAS family in pigeonpea creating a significant resource for future functional studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03782-x.

Overexpressing CsPSY1 Gene of Tea Plant, Encoding a Phytoene Synthase, Improves α-Carotene and β-Carotene Contents in Carrot

Tea plants (Camellia sinensis (L.) O. Kuntze) belong to Theaceae family, in the section Thea. Tea plants are widely distributed in subtropical and tropical regions in the word. α-carotene and β-carotene in the tea leaves belong to carotenoids, which are associated with the aroma and color of the tea. Phytoene synthase (PSY) is a rate-limiting enzyme in carotenoids biosynthesis. We identified three CsPSY genes in 'Shuchazao', named CsPSY1, CsPSY2, and CsPSY3. Structural analysis of three CsPSY genes showed that CsPSY1 had a longer intro structure. The cis-acting elements of CsPSYs promoter were mainly associated with light-responsiveness, abiotic stress-responsiveness, and hormone-responsiveness. CsPSY1 exhibited expression in all tissues of the tea plants, whereas CsPSY2 and CsPSY3 were trace expression levels in all tissues. The positive expression of CsPSY1 under hormonal and abiotic stresses suggested its role in plant development and defense responses. The amino acid sequence of CsPSY1 was highly conserved in eight tea cultivars. The recombinant vector pCAMBIA1301-CsPSY1 was constructed to stabilize the overexpression of CsPSY1 in carrot. The contents of α-carotene and β-carotene in transgenic carrot callus were significantly increased. This study provides a foundational basis for further research on the function of CsPSYs and carotenoids accumulation in tea plants.

Genome-wide analysis of the GRAS gene family in Liriodendron chinense reveals the putative function in abiotic stress and plant development

Introduction

GRAS genes encode plant-specific transcription factors that play essential roles in plant growth and development. However, the members and the function of the GRAS gene family have not been reported in Liriodendron chinense. L. chinense, a tree species in the Magnolia family that produces excellent timber for daily life and industry. In addition, it is a good relict species for plant evolution research.

Methods

Therefore, we conducted a genome-wide study of the LcGRAS gene family and identified 49 LcGRAS genes in L. chinense.

Results

We found that LcGRAS could be divided into 13 sub-groups, among which there is a unique branch named HAM-t. We carried out RNA sequencing analysis of the somatic embryos from L. chinense and found that LcGRAS genes are mainly expressed after heart-stage embryo development, suggesting that LcGRAS may have a function during somatic embryogenesis. We also investigated whether GRAS genes are responsive to stress by carrying out RNA sequencing (RNA-seq) analysis, and we found that the genes in the PAT subfamily were activated upon stress treatment, suggesting that these genes may help plants survive stressful environments. We found that PIF was downregulated and COR was upregulated after the transient overexpression of PATs, suggesting that PAT may be upstream regulators of cold stress.

Discussion

Collectively, LcGRAS genes are conserved and play essential roles in plant development and adaptation to abiotic stress.

Genome-wide identification and expression analysis of the GRAS family under low-temperature stress in bananas

Introduction

GRAS, named after GAI, RGA, and SCR, is a class of plant-specific transcription factors family that plays a crucial role in growth and development, signal transduction, and various stress responses.

Methods

To understand the biological functions of the banana GRAS gene family, a genome-wide identification and bioinformatics analysis of the banana GRAS gene family was performed based on information from the M. acuminata, M. balbisiana, and M. itinerans genomic databases.

Result

In the present study, we identified 73 MaGRAS, 59 MbGRAS, and 58 MiGRAS genes in bananas at the whole-genome scale, and 56 homologous genes were identified in the three banana genomes. Banana GRASs can be classified into 10 subfamilies, and their gene structures revealed that most banana GRAS gDNAs lack introns. The promoter sequences of GRASs had a large number of cis-acting elements related to plant growth and development, phytohormone, and adversity stress responsiveness. The expression pattern of seven key members of MaGRAS response to low-temperature stress and different tissues was also examined by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The microRNAs-MaGRASs target prediction showed perfect complementarity of seven GRAS genes with the five mac-miRNAs. The expression of all seven genes was lowest in roots, and the expression of five genes was highest in leaves during low-temperature stress. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 was significantly lower under low-temperature stress compared to the control, except for MaSCL27-2, which was slightly higher than the 28°C control at 4 h. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 dropped to the lowest levels at 24 h, 12 h, and 4 h, respectively. The MaSCL27-4 and MaSCL6-2 expression was intermittently upregulated, rising to the highest expression at 24h, while the expression of MaSCL22 was less variable, remaining at the control level with small changes.

Discussion

In summary, it is tentatively hypothesized that the GRAS family has an important function in low-temperature stress in bananas. This study provides a theoretical basis for further analyzing the function of the banana GRAS gene and the resistance of bananas to cold temperatures.

Light Intensity Regulates Low-Temperature Adaptability of Tea Plant through ROS Stress and Developmental Programs

Low-temperature stress limits global tea planting areas and production efficiency. Light is another essential ecological factor that acts in conjunction with temperature in the plant life cycle. However, it is unclear whether the differential light environment affects the low temperature adaptability of tea plant (Camellia sect. Thea). In this study, tea plant materials in three groups of light intensity treatments showed differentiated characteristics for low-temperature adaptability. Strong light (ST, 240 μmol·m^-2·s^-1) caused the degradation of chlorophyll and a decrease in peroxidase (POD), superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), and polyphenol oxidase (PPO) activities, as well as an increase in soluble sugar, soluble protein, malondialdehyde (MDA), and relative conductivity in tea leaves. In contrast, antioxidant enzyme activities, chlorophyll content, and relative conductivity were highest in weak light (WT, 15 μmol·m^-2·s^-1). Damage was observed in both ST and WT materials relative to moderate light intensity (MT, 160 μmol·m^-2·s^-1) in a frost resistance test. Chlorophyll degradation in strong light was a behavior that prevented photodamage, and the maximum photosynthetic quantum yield of PS II (Fv/Fm) decreased with increasing light intensity. This suggests that the browning that occurs on the leaf surface of ST materials through frost may have been stressed by the previous increase in reactive oxygen species (ROS). Frost intolerance of WT materials is mainly related to delayed tissue development and tenderness holding. Interestingly, transcriptome sequencing revealed that stronger light favors starch biosynthesis, while cellulose biosynthesis is enhanced in weaker light. It showed that light intensity mediated the form of carbon fixation in tea plant, and this was associated with low-temperature adaptability.

Genome-wide characterization of the U-box gene in Camellia sinensis and functional analysis in transgenic tobacco under abiotic stresses

Plants U-box genes are crucial for plant survival, and they extensively regulate plant growth, reproduction and development as well as coping with stress and other processes. In this study, we identified 92 CsU-box genes through genome-wide analysis in the tea plant (Camellia sinensis), all of them contained the conserved U-box domain and were divided into 5 groups, which supported by the further genes structure analysis. The expression profiles in eight tea plant tissues and under abiotic and hormone stresses were analyzed using the TPIA database. 7 CsU-box genes (CsU-box27/28/39/46/63/70/91) were selected to verify and analyze expression patterns under PEG-induced drought and heat stress in tea plant respectively, the qRT-PCR results showed consistent with transcriptome datasets; and the CsU-box39 were further heterologous expressed in tobacco to perform gene function analysis. Phenotypic analyses of overexpression transgenic tobacco seedlings and physiological experiments revealed that CsU-box39 positively regulated the plant response to drought stress. These results lay a solid foundation for studying the biological function of CsU-box, and will provide breeding strategy basis for tea plant breeders.

Molecular regulation of immunity in tea plants

Tea, which is mainly produced using the young leaves and buds of tea plants (Camellia sinensis (L.) O. Kuntze), is one of the most common non-alcoholic beverages consumed in the world. The standard of tea mostly depends on the variety and quality of tea plants, which generally grow in subtropical areas, where the warm and humid conditions are also conducive to the occurrence of diseases. In fighting against pathogens, plants rely on their sophisticated innate immune systems which has been extensively studied in model plants. Many components involved in pathogen associated molecular patterns (PAMPs) triggered immunity (PTI) and effector triggered immunity (ETI) have been found. Nevertheless, the molecular regulating network against pathogens (e.g., Pseudopestalotiopsis sp., Colletotrichum sp. and Exobasidium vexans) causing widespread disease (such as grey blight disease, anthracnose, and blister blight) in tea plants is still unclear. With the recent release of the genome data of tea plants, numerous genes involved in tea plant immunity have been identified, and the molecular mechanisms behind tea plant immunity is being studied. Therefore, the recent achievements in identifying and cloning functional genes/gene families, in finding crucial components of tea immunity signaling pathways, and in understanding the role of secondary metabolites have been summarized and the opportunities and challenges in the future studies of tea immunity are highlighted in this review.

Genome-wide analysis and identification of microRNAs in Medicago truncatula under aluminum stress

Numerous studies have shown that plant microRNAs (miRNAs) play key roles in plant growth and development, as well as in response to biotic and abiotic stresses; however, the role of miRNA in legumes under aluminum (Al) stress have rarely been reported. Therefore, here, we aimed to investigate the role of miRNAs in and their mechanism of Al tolerance in legumes. To this end, we sequenced a 12-strand-specific library of Medicago truncatula under Al stress. A total of 195.80 M clean reads were obtained, and 876 miRNAs were identified, of which, 673 were known miRNAs and 203 were unknown. A total of 55 miRNAs and their corresponding 2,502 target genes were differentially expressed at various time points during Al stress. Further analysis revealed that mtr-miR156g-3p was the only miRNA that was significantly upregulated at all time points under Al stress and could directly regulate the expression of genes associated with root cell growth. Three miRNAs, novel_miR_135, novel_miR_182, and novel_miR_36, simultaneously regulated the expression of four Al-tolerant transcription factors, GRAS, MYB, WRKY, and bHLH, at an early stage of Al stress, indicating a response to Al stress. In addition, legume-specific miR2119 and miR5213 were involved in the tolerance mechanism to Al stress by regulating F-box proteins that have protective effects against stress. Our results contribute to an improved understanding of the role of miRNAs in Al stress in legumes and provide a basis for studying the molecular mechanisms of Al stress regulation.

Topics and trends in fresh tea (Camellia sinensis) leaf research: A comprehensive bibliometric study

Tea plant (Camellia sinensis) is a widely cultivated cash crop and tea is a favorite functional food in the world. Fresh tea leaves (FTLs) play a critical role in bridging the two fields closely related to tea cultivation and tea processing, those are, tea plant biology and tea biochemistry. To provide a comprehensive overview of the development stages, authorship collaboration, research topics, and hotspots and their temporal evolution trends in the field of FTLs research, we conducted a bibliometric analysis, based on 971 publications on FTLs-related research published during 2001-2021 from Web of Science Core Collection. CiteSpace, R package Bibliometrix, and VOSviewer were employed in this research. The results revealed that the development history can be roughly divided into three stages, namely initial stage, slow development stage and rapid development stage. Journal of Agricultural & Food Chemistry published most articles in this field, while Frontiers in Plant Science held the highest total citations and h-index. The most influential country, institution, and author in this field was identified as China, the Chinese Academy of Agricultural Sciences, and Xiaochun Wan, respectively. FTLs-related research can be categorized into three main topics: the regulation mechanism of key genes, the metabolism and features of essential compounds, and tea plants' growth and stress responses. The most concerning hotspots are the application of advanced technologies, essential metabolites, leaf color variants, and effective cultivation treatments. There has been a shift from basic biochemical and enzymatic studies to studies of molecular mechanisms that depend on multi-omics technologies. We also discussed the future development in this field. This study provides a comprehensive summary of the research field, making it easier for researchers to be informed about its development history, status, and trends.

Genome-wide identification and expression analysis of the GRAS gene family in Dendrobium chrysotoxum

The GRAS gene family encodes transcription factors that participate in plant growth and development phases. They are crucial in regulating light signal transduction, plant hormone (e.g. gibberellin) signaling, meristem growth, root radial development, response to abiotic stress, etc. However, little is known about the features and functions of GRAS genes in Orchidaceae, the largest and most diverse angiosperm lineage. In this study, genome-wide analysis of the GRAS gene family was conducted in Dendrobium chrysotoxum (Epidendroideae, Orchidaceae) to investigate its physicochemical properties, phylogenetic relationships, gene structure, and expression patterns under abiotic stress in orchids. Forty-six DchGRAS genes were identified from the D. chrysotoxum genome and divided into ten subfamilies according to their phylogenetic relationships. Sequence analysis showed that most DchGRAS proteins contained conserved VHIID and SAW domains. Gene structure analysis showed that intronless genes accounted for approximately 70% of the DchGRAS genes, the gene structures of the same subfamily were the same, and the conserved motifs were also similar. The Ka/Ks ratios of 12 pairs of DchGRAS genes were all less than 1, indicating that DchGRAS genes underwent negative selection. The results of cis-acting element analysis showed that the 46 DchGRAS genes contained a large number of hormone-regulated and light-responsive elements as well as environmental stress-related elements. In addition, the real-time reverse transcription quantitative PCR (RT-qPCR) experimental results showed significant differences in the expression levels of 12 genes under high temperature, drought and salt treatment, among which two members of the LISCL subfamily (DchGRAS13 and DchGRAS15) were most sensitive to stress. Taken together, this paper provides insights into the regulatory roles of the GRAS gene family in orchids.

Genome-wide identification, expression and salt stress tolerance analysis of the GRAS transcription factor family in Betula platyphylla

The GRAS gene family is a plant-specific family of transcription factors and play a vital role in many plant growth processes and abiotic stress responses. Nevertheless, the functions of the GRAS gene family in woody plants, especially in Betula platyphylla (birch), are hardly known. In this study, we performed a genome-wide analysis of 40 BpGRAS genes (BpGRASs) and identified typical GRAS domains of most BpGRASs. The BpGRASs were unevenly distributed on 14 chromosomes of birch and the phylogenetic analysis of six species facilitated the clustering of 265 GRAS proteins into 17 subfamilies. We observed that closely related GRAS homologs had similar conserved motifs according to motif analysis. Besides, an analysis of the expression patterns of 26 BpGRASs showed that most BpGRASs were highly expressed in the leaves and responded to salt stress. Six BpGRASs were selected for cis-acting element analysis because of their significant upregulation under salt treatment, indicating that many elements were involved in the response to abiotic stress. This result further confirmed that these BpGRASs might participate in response to abiotic stress. Transiently transfected birch plants with transiently overexpressed 6 BpGRASs and RNAi-silenced 6 BpGRASs were generated for gain- and loss-of-function analysis, respectively. In addition, overexpression of BpGRAS34 showed phenotype resistant to salt stress, decreased the cell death and enhanced the reactive oxygen species (ROS) scavenging capabilities and proline content under salt treatment, consistent with the results in transiently transformed birch plants. This study is a systematic analysis of the GRAS gene family in birch plants, and the results provide insight into the molecular mechanism of the GRAS gene family responding to abiotic stress in birch plants.

Genome-Wide Identification and Expression Analysis of Isopentenyl transferase Family Genes during Development and Resistance to Abiotic Stresses in Tea Plant (Camellia sinensis)

The tea plant is an important economic crop and is widely cultivated. Isopentenyl transferase (IPT) is the first and rate-limiting enzyme of cytokinin (CK) signaling, which plays key roles in plant development and abiotic stress. However, the IPT gene family in tea plants has not been systematically investigated until now. The phylogenetic analyses, gene structures, and conserved domains were predicted here. The results showed that a total of 13 CsIPT members were identified from a tea plant genome database and phylogenetically classified into four groups. Furthermore, 10 CsIPT members belonged to plant ADP/ATP-IPT genes, and 3 CsIPTs were tRNA-IPT genes. There is a conserved putative ATP/GTP-binding site (P-loop motif) in all the CsIPT sequences. Based on publicly available transcriptome data as well as through RNA-seq and qRT-PCR analysis, the CsIPT genes which play key roles in the development of different tissues were identified, respectively. Furthermore, CsIPT6.2 may be involved in the response to different light treatments. CsIPT6.4 may play a key role during the dormancy and flush of the lateral buds. CsIPT5.1 may play important regulatory roles during the development of the lateral bud, leaf, and flower. CsIPT5.2 and CsIPT6.2 may both play key roles for increased resistance to cold-stress, whereas CsIPT3.2 may play a key role in improving resistance to high-temperature stress as well as drought-stress and rewatering. This study could provide a reference for further studies of CsIPT family’s functions and could contribute to tea molecular breeding.

Isolation and Characterization of an LBD Transcription Factor CsLBD39 from Tea Plant (Camellia sinensis) and Its Roles in Modulating Nitrate Content by Regulating Nitrate-Metabolism-Related Genes

Nitrate nitrogen is an important nitrogen source for tea plants’ growth and development. LBD transcription factors play important roles in response to the presence of nitrate in plants. The functional study of LBD transcription factors in tea plants remains limited. In this study, the LBD family gene CsLBD39 was isolated and characterized from tea plants. Sequence analysis indicated that CsLBD39 contained a highly conserved CX₂CX₆CX₃CX domain. The phylogenetic tree assay showed that CsLBD39 belonged to class II subfamily of the LBD family. CsLBD39 was highly expressed in flowers and root; we determined that its expression could be induced by nitrate treatment. The CsLBD39 protein was located in the nucleus and has transcriptional activation activity in yeast. Compared with the wild type, overexpression of CsLBD39 gene in Arabidopsis resulted in smaller rosettes, shorter main roots, reduced lateral roots and lower plant weights. The nitrate content and the expression levels of genes related to nitrate transport and regulation were decreased in transgenic Arabidopsis hosting CsLBD39 gene. Compared with the wild type, CsLBD39 overexpression in transgenic Arabidopsis had smaller cell structure of leaves, shorter diameter of stem cross section, and slender and compact cell of stem longitudinal section. Under KNO₃ treatment, the contents of nitrate, anthocyanins, and chlorophyll in leaves, and the content of nitrate in roots of Arabidopsis overexpressing CsLBD39 were reduced, the expression levels of nitrate transport and regulation related genes were decreased. The results revealed that CsLBD39 may be involved in nitrate signal transduction in tea plants as a negative regulator and laid the groundwork for future studies into the mechanism of nitrate response.

Genome-wide study of the GRAS gene family in Hibiscus hamabo Sieb. et Zucc and analysis of HhGRAS14-induced drought and salt stress tolerance in Arabidopsis

GRAS proteins are widely distributed plant-specific transcription factors. In this study, we identified 59 GRAS proteins (HhGRASs) from the genomic and transcriptomic datasets of Hibiscus hamabo Sieb. et Zucc. These proteins were phylogenetically divided into nine subfamilies. RNA-seq analysis revealed that most HhGRASs were expressed in response to abiotic stresses. Results from quantitative real-time PCR analysis of nine selected HhGRASs suggested that HhGRAS14 was significantly upregulated under multiple abiotic stresses; therefore, this gene was selected for further study. Silencing HhGRAS14 in H. hamabo reduced the tolerance to drought and salt stress, while overexpression in Arabidopsis thaliana significantly increased the tolerance to drought and salt and reduced the sensitivity to abscisic acid (ABA). In summary, we analyzed the GRAS family of proteins in semi-mangrove plants for the first time and identified a gene that responds to drought and salt stress, which provided the basis for a comprehensive analysis of GRAS genes and insight into the abiotic stress response mechanism in H. hamabo.

Genome-wide survey and expression analysis of GRAS transcription factor family in sweetpotato provides insights into their potential roles in stress response

BACKGROUND

The plant-specific GRAS transcription factors play pivotal roles in various adverse environmental conditions. Numerous GRAS genes have been explored and characterized in different plants, however, comprehensive survey on GRASs in sweetpotato is lagging.

RESULTS

In this study, 72 putative sweetpotato IbGRAS genes with uneven distribution were isolated on 15 chromosomes and classified into 12 subfamilies supported by gene structures and motif compositions. Moreover, both tandem duplication and segmental duplication events played critical roles in the expansion of sweetpotato GRAS genes, and the collinearity between IbGRAS genes and the related orthologs from nine other plants further depicted evolutionary insights into GRAS gene family. RNA-seq analysis under salt stress and qRT-PCR detection of 12 selected IbGRAS genes demonstrated their significant and varying inductions under multiple abiotic stresses (salt, drought, heat and cold) and hormone treatments (ABA, ACC and JA). Consistently, the promoter regions of IbGRAS genes harbored a series of stress- and hormone-associated cis-acting elements. Among them, IbGRAS71, the potential candidate for breeding tolerant plants, was characterized as having transactivation activity in yeasts, while IbGRAS-2/-4/-9 did not. Moreover, a complex interaction relationship between IbGRASs was observed through the interaction network analysis and yeast two-hybrid assays.

CONCLUSIONS

Our results laid a foundation for further functional identifications of IbGRAS genes, and multiple members may serve as potential regulators for molecular breeding of tolerant sweetpotato.

Genome-Wide Analyses of Tea Plant Stress-Associated Proteins (SAPs) Reveal the Role of CsSAP12 in Increased Drought Tolerance in Transgenic Tomatoes

Plant stress-associated proteins (SAPs) contain A20/AN1 zinc finger domains and are involved in plant response to abiotic stresses. In this study, we aimed to explore the biological function of tea plant CsSAPs. A total of 14 CsSAP genes were identified in the tea plant genome using a reference genome database (Camellia sinensis var. sinensis). The CsSAPs were divided into the following two groups: Group I, containing one AN1 domain and/or one A20 domain; and Group II, containing two AN1 domains and/or two C2H2 domains. The sequence alignments and conserved domains analysis indicated that the CsSAPs were highly structurally conserved in terms of amino acid sequence and protein structure. The CsSAPs showed different transcript levels in spatio-temporal expression and in response to cold and drought stress in tea plants. Furthermore, the expression of CsSAP12 was considerably upregulated under drought stress. The overexpression of CsSAP12 in transgenic tomatoes showed increased tolerance to drought stress compared with the wild type. Altogether, the results showed that CsSAP12 might be involved in drought stress. Thus, CsSAP12 might be a target gene in genetic engineering to improve drought tolerance in tea plants.

Transcription factor CsWRKY40 regulates L-theanine hydrolysis by activating the CsPDX2.1 promoter in tea leaves during withering

L-Theanine is a crucial secondary metabolite in tea and positively determines the potential quality and health benefits of tea products. Previous work found the content of L-theanine decreased during withering process, while the specific mechanism is still unknown. Here, weighted gene co-expression network analysis (WGCNA) was performed based on the transcriptome data obtained previously. The key hydrolysis gene CsPDX2.1 in L-theanine metabolism and seven candidate transcription factors were screened out. Among those transcription factors, CsWRKY40 presented the strongest activation on the CsPDX2.1 promoter (373.18-fold) by binding to W box element based on the dual luciferase assay and EMSA results. Meanwhile, CsWRKY40 protein was located in the nucleoplasm, while CsPDX2.1 was found in both the nucleoplasm and cytoplasm. Furthermore, it was confirmed that the water loss of tea leaves was the critical factor affecting the contents of ABA and L-theanine by activating the expression of CsPDX2.1 and CsPDX2.1 based on the analysis of the withering model, water-retention model and water-loss model. Our results provide a new insight into revealing the regulation mechanism of L-theanine hydrolysis metabolism.

Genome-wide identification studies - A primer to explore new genes in plant species

Genome data have accumulated rapidly in recent years, doubling roughly after every 6 months due to the influx of next-generation sequencing technologies. A plethora of plant genomes are available in comprehensive public databases. This easy access to data provides an opportunity to explore genome datasets and recruit new genes in various plant species not possible a decade ago. In the past few years, many gene families have been published using these public datasets. These genome-wide studies identify and characterize gene members, gene structures, evolutionary relationships, expression patterns, protein interactions and gene ontologies, and predict putative gene functions using various computational tools. Such studies provide meaningful information and an initial framework for further functional elucidation. This review provides a concise layout of approaches used in these gene family studies and demonstrates an outline for employing various plant genome datasets in future studies.

Exploring the GRAS gene family in common bean (Phaseolus vulgaris L.): characterization, evolutionary relationships, and expression analyses in response to abiotic stresses

Genome-wide identification reveals 55 PvuGRAS genes belonging to 16 subfamilies and their gene structures and evolutionary relationships were characterized. Expression analyses highlight their prominence in plant growth, development and abiotic stress responses. GRAS proteins comprise a plant-specific transcription factor family involved in multiple growth regulatory pathways and environmental cues including abiotic/biotic stresses. Despite its crucial importance, characterization of this gene family is still elusive in common bean. A systematic genome-wide scan identified 55 PvuGRAS genes unevenly anchored to the 11 common bean chromosomes. Segmental duplication appeared to be the key driving force behind expansion of this gene family that underwent purifying selection during evolution. Computational investigation unraveled their intronless organization and identified similar motif composition within the same subfamily. Phylogenetic analyses clustered the PvuGRAS proteins into 16 phylogenetic clades and established extensive orthologous relationships with Arabidopsis and rice. Analysis of the upstream promoter region uncovered cis-elements responsive to growth, development, and abiotic stresses that may account for their differential expression. The identified SSRs could serve as putative molecular markers facilitating future breeding programs. 37 PvuGRAS transcripts were post-transcriptionally regulated by different miRNA families, miR171 being the major player preferentially targeting members of the HAM subfamily. Global expression profile based on RNA-seq data indicates a clade specific expression pattern in various tissues and developmental stages. Additionally, nine PvuGRAS genes were chosen for further qPCR analyses under drought, salt, and cold stress suggesting their involvement in acclimation to environmental stimuli. Combined, the present results significantly contribute to the current understanding of the complexity and biological function of the PvuGRAS gene family. The resources generated will provide a solid foundation in future endeavors for genetic improvement in common bean.

The GRAS gene family and its roles in seed development in litchi (Litchi chinensis Sonn)

BACKGROUND

The GRAS gene family plays crucial roles in multiple biological processes of plant growth, including seed development, which is related to seedless traits of litchi (Litchi chinensis Sonn.). However, it hasn't been fully identified and analyzed in litchi, an economic fruit tree cultivated in subtropical regions.

RESULTS

In this study, 48 LcGRAS proteins were identified and termed according to their chromosomal location. LcGRAS proteins can be categorized into 14 subfamilies through phylogenetic analysis. Gene structure and conserved domain analysis revealed that different subfamilies harbored various motif patterns, suggesting their functional diversity. Synteny analysis revealed that the expansion of the GRAS family in litchi may be driven by their tandem and segmental duplication. After comprehensively analysing degradome data, we found that four LcGRAS genes belong to HAM subfamily were regulated via miR171-mediated degradation. The various expression patterns of LcGRAS genes in different tissues uncovered they were involved in different biological processes. Moreover, the different temporal expression profiles of LcGRAS genes between abortive and bold seed indicated some of them were involved in maintaining the normal development of the seed.

CONCLUSION

Our study provides comprehensive analyses on GRAS family members in litchi, insight into a better understanding of the roles of GRAS in litchi development, and lays the foundation for further investigations on litchi seed development.

Comprehensive identification and expression analysis of GRAS gene family under abiotic stress and phytohormone treatments in Pearl millet

Pearl millet is an important C4 cereal plant that possesses enormous capacity to survive under extreme climatic conditions. It serves as a major food source for people in arid and semiarid regions of south-east Asia and Africa. GRAS is an important transcription factor gene family of plant that play a critical role in regulating developmental processes, stress responses and phytohormonal signalling. In the present study, we have identified a total number of 57 GRAS members in pearl millet. Phylogenetic analysis clustered all the PgGRAS genes into eight groups (GroupI-GroupVIII). Motif analysis has shown that all the PgGRAS proteins had conserved GRAS domains and gene structure analysis revealed a high structural diversity among PgGRAS genes. Expression patterns of PgGRAS genes in different tissues (leaf, stem and root) and under various abiotic stress (drought, heat and salinity) were determined. Further, expression analysis was also carried out in response to various hormones (SA, MeJA, GA and ABA). The results provide a clear understanding of GRAS transcription factor family in pearl millet, and lay a good foundation for the functional characterisation of GRAS genes in pearl millet.

Methane-rich water induces bulblet formation of scale cuttings in Lilium davidii var. unicolor by regulating the signal transduction of phytohormones and their levels

Previous studies have shown that methane (CH₄ ) has promoting roles in the adventitious root (AR) and lateral root (LR) formation in plants. However, whether CH₄ could trigger the bulblet formation in scale cutting of Lilium davidii var. unicolor has not been elucidated. To gain insight into the effect of CH₄ on the bulblet formation, different concentrations (1, 10, 50, and 100%) of methane-rich water (MRW) and distilled water were applied to treat the scale cuttings of Lilium. We observed that treatment with 100% MRW obviously induced the bulblet formation in scale cuttings. To explore the mechanism of CH₄ -induced bulblet formation, the transcriptome of scales was analyzed. A total of 2078 differentially expressed genes (DEGs) were identified. The DEGs were classified into different metabolic pathways, especially phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant signal transduction. Of these, approximately 38 candidate DEGs involved in the plant signal transduction were further studied. In addition, the expression of AP2-ERF/ERF, WRKY, GRAS, ARF, and NAC transcription factors (TFs) was changed by MRW treatment, suggesting their potential involvement in bulblet formation. As for hormones, exogenous IAA, GA, and ABA could induce the bulblet formation. Additional experiments suggested that MRW could increase the endogenous IAA, GA, and JA levels, but decrease the levels of ABA during bulblet formation, which showed that higher IAA, GA, JA levels and lower ABA content might facilitate bulblet formation. In addition, the levels of endogenous hormones were consistent with the expression level of genes involved in phytohormone signal transduction. Overall, this study has revealed that CH₄ might improve the bulblet formation of cutting scales in Lilium by regulating the expression of genes related to phytohormone signal transduction and TFs, as well as by changing the endogenous hormone levels.

Genome wide identification and expression pattern analysis of the GRAS family in quinoa

GRAS, a key transcription factor in plant growth and development, has not yet been reported in quinoa. Therefore, this study used the latest quinoa genomic data to identify and analyse GRAS genes in quinoa: 52 GRAS genes were identified in quinoa, these being unevenly distributed on 19 chromosomes. Fragment duplication and tandem duplication events were the main reasons for the expansion of the GRAS gene family in quinoa. Protein sequence analysis showed that there were some differences in amino acid numbers and isoelectric points amongst different subfamilies, and the main secondary structures were α-helix and random coil. The CqGRAS gene was divided into 14 subfamilies based on results from phylogenetic analysis. The genes located in the same subfamily had similar gene structures, conserved motifs, and three-level models. Promoter region analysis showed that the GRAS family genes contained multiple homeostasis elements that responded to hormones and adversity. GO enrichment indicated that CqGRAS genes were involved in biological processes, cell components, and molecular functions. By analysing the expression of CqGRAS genes in different tissues and different treatments, it was found that GRAS genes had obvious differential expression in different tissues and stress, which indicates that GRAS genes had tissue or organ expression specificity and thus might play an important role in response to stress. These results laid a foundation for further functional research on the GRAS gene family in quinoa.

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.

Genome Identification and Expression Analysis of GRAS Family Related to Development, Hormone and Pathogen Stress in Brachypodium distachyon

GRAS transcription factors are widely present in the plant kingdom and play important roles in regulating multiple plant physiological processes. Brachypodium distachyon is a model for grasses for researching plant-pathogen interactions. However, little is known about the BdGRAS family genes involved in plant response to biotic stress. In this study, we identified 63 genes of the GRAS family in B. distachyon. The phylogenetic analysis showed that BdGRAS genes were divided into ten subfamilies and unevenly distributed on five chromosomes. qRT-PCR results showed that the BdGRAS family genes were involved in the growth and development of B. distachyon. Moreover, the expression of the HAM subfamily genes of BdGRAS changed during the interaction between B. distachyon and Magnaporthe oryzae. Interestingly, BdGRAS31 in the HAM subfamily was regulated by miR171 after inoculation with M. oryzae. These results provide insight into the potential functions of the BdGRAS family in disease resistance.

Conservation and diversification of HAIRY MERISTEM gene family in land plants

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

MAPK cascade gene family in Camellia sinensis: In-silico identification, expression profiles and regulatory network analysis

BACKGROUND

Mitogen Activated Protein Kinase (MAPK) cascade is a fundamental pathway in organisms for signal transduction. Though it is well characterized in various plants, there is no systematic study of this cascade in tea.

RESULT

In this study, 5 genes of Mitogen Activated Protein Kinase Kinase (MKK) and 16 genes of Mitogen Activated Protein Kinase (MPK) in Camellia sinensis were found through a genome-wide search taking Arabidopsis thaliana as the reference genome. Also, phylogenetic relationships along with structural analysis which includes gene structure, location as well as protein conserved motifs and domains, were systematically examined and further, predictions were validated by the results. The plant species taken for comparative study clearly displayed segmental duplication, which was a significant candidate for MAPK cascade expansion. Also, functional interaction was carried out in C. sinensis based on the orthologous genes in Arabidopsis. The expression profiles linked to various stress treatments revealed wide involvement of MAPK and MAPKK genes from Tea in response to various abiotic factors. In addition, the expression of these genes was analysed in various tissues.

CONCLUSION

This study provides the targets for further comprehensive identification, functional study, and also contributed for a better understanding of the MAPK cascade regulatory network in C. sinensis.

Genome-wide identification and characterization of GRAS genes in soybean (Glycine max)

BACKGROUND

GRAS proteins are crucial transcription factors, which are plant-specific and participate in various plant biological processes. Thanks to the rapid progress of the whole genome sequencing technologies, the GRAS gene families in different plants have been broadly explored and studied. However, comprehensive research on the soybean (Glycine max) GRAS gene family is relatively lagging.

RESULTS

In this study, 117 Glycine max GRAS genes (GmGRAS) were identified. Further phylogenetic analyses showed that the GmGRAS genes could be categorized into nine gene subfamilies: DELLA, HAM, LAS, LISCL, PAT1, SCL3, SCL4/7, SCR and SHR. Gene structure analyses turned out that the GmGRAS genes lacked introns and were relatively conserved. Conserved domains and motif patterns of the GmGRAS members in the same subfamily or clade exhibited similarities. Notably, the expansion of the GmGRAS gene family was driven both by gene tandem and segmental duplication events. Whereas, segmental duplications took the major role in generating new GmGRAS genes. Moreover, the synteny and evolutionary constraints analyses of the GRAS proteins among soybean and distinct species (two monocots and four dicots) provided more detailed evidence for GmGRAS gene evolution. Cis-element analyses indicated that the GmGRAS genes may be responsive to diverse environmental stresses and regulate distinct biological processes. Besides, the expression patterns of the GmGRAS genes were varied in various tissues, during saline and dehydration stresses and during seed germination processes.

CONCLUSIONS

We conducted a systematic investigation of the GRAS genes in soybean, which may be valuable in paving the way for future GmGRAS gene studies and soybean breeding.

The Apple microR171i-SCARECROW-LIKE PROTEINS26.1 Module Enhances Drought Stress Tolerance by Integrating Ascorbic Acid Metabolism

Drought stress severely restricts crop yield and quality. Small noncoding RNAs play critical roles in plant growth, development, and stress responses by regulating target gene expression, but their roles in drought stress tolerance in apple (Malus domestica) are poorly understood. Here, we identified various small noncoding RNAs and their targets from the wild apple species Malus sieversii via high-throughput sequencing and degradome analysis. Forty known microRNAs (miRNAs) and eight new small noncoding RNAs were differentially expressed in response to 2 or 4 h of drought stress treatment. We experimentally verified the expression patterns of five selected miRNAs and their targets. We established that one miRNA, mdm-miR171i, specifically targeted and degraded SCARECROW-LIKE PROTEINS26 1 (MsSCL26 1) transcripts. Both knockout of mdm-miR171i and overexpression of MsSCL26 1 improved drought stress tolerance in the cultivated apple line 'GL-3' by regulating the expression of antioxidant enzyme genes, especially that of MONODEHYDROASCORBATE REDUCTASE, which functions in metabolism under drought stress. Transient expression analysis demonstrated that MsSCL26.1 activates MsMDHAR transcription by positively regulating the activity of the P1 region in its promoter. Therefore, the miR171i-SCL26 1 module enhances drought stress tolerance in apple by regulating antioxidant gene expression and ascorbic acid metabolism.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance

BACKGROUND

Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.).

RESULTS

In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca²⁺. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response.

CONCLUSIONS

A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

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

Cucumber (Cucumis sativus L.) is one of the most important cucurbit vegetables but is often subjected to stress during cultivation. GRAS (gibberellic acid insensitive, repressor of GAI, and scarecrow) genes encode a family of transcriptional factors that regulate plant growth and development. In the model plant Arabidopsis thaliana, GRAS family genes function in formation of axillary meristem and root radial structure, phytohormone (gibberellin) signal transduction, light signal transduction and abiotic/biological stress. In this study, a gene family was comprehensively analyzed from the aspects of evolutionary tree, gene structure, chromosome location, evolutionary and expression pattern by means of bioinformatics; 37 GRAS gene family members have been screened from cucumber. We reconstructed an evolutionary tree based on multiple sequence alignment of the typical GRAS domain and conserved motif sequences with those of other species (A. thaliana and Solanum lycopersicum). Cucumber GRAS family was divided into 10 groups according to the classification of Arabidopsis and tomato genes. We conclude that tandem and segmental duplication have played important roles in the expansion and evolution of the cucumber GRAS (CsaGRAS) family. Expression patterns of CsaGRAS genes in different tissues and under cold treatment, combined with gene ontology annotation and interaction network analysis, revealed potentially different functions for CsaGRAS genes in response to cold tolerance, with members of the SHR, SCR and DELLA subfamilies likely playing important roles. In conclusion, this study provides valuable information and candidate genes for improving cucumber tolerance to cold stress.

Genome-Wide Analysis of the GRAS Gene Family in Barley (Hordeum vulgare L.)

The GRAS (named after first three identified proteins within this family, GAI, RGA, and SCR) family contains plant-specific genes encoding transcriptional regulators that play a key role in gibberellin (GA) signaling, which regulates plant growth and development. Even though GRAS genes have been characterized in some plant species, little research is known about the GRAS genes in barley (Hordeum vulgare L.). In this study, we observed 62 GRAS members from barley genome, which were grouped into 12 subgroups by using phylogenomic analysis together with the GRAS genes from Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and rice (Oryza sativa). Chromosome localization and gene structure analysis suggested that duplication events and abundant presence of intronless genes might account for the massive expansion of GRAS gene family in barley. The analysis of RNA-seq data indicates the expression pattern of GRAS genes in various tissues at different stages in barley. Noteworthy, our qRT-PCR analysis showed the expression of 18 candidate GRAS genes abundantly in the developing inflorescence, indicating their potential roles in the barley inflorescence development and reproduction. Collectively, our evolutionary and expression analysis of GRAS family are useful for future functional characterization of GA signaling in barley and agricultural improvement.

Genome-wide investigation of WRKY transcription factors in sweet osmanthus and their potential regulation of aroma synthesis

WRKY transcription factors, one of the largest transcription factor families, play important roles in regulating the synthesis of secondary metabolites. In sweet osmanthus (Osmanthus fragrans), the monoterpenes have been demonstrated as the most important volatile compounds, and the W-box, which is the cognate binding site of WRKY transcription factors, could be identified in most of the terpene-synthesis-related genes' promoters. However, the role of the WRKY family in terpene synthesis in sweet osmanthus has rarely been examined. In this study, 154 WRKY genes with conserved WRKY domain were identified and classified into three groups. The group II was further divided into five subgroups, and almost all members of IId contained a plant zinc cluster domain. Eight OfWRKYs (OfWRKY7/19/36/38/42/84/95/139) were screened from 20 OfWRKYs for their flower-specific expression patterns in different tissues. Simultaneously, the expression patterns of OfWRKYs and emission patterns of volatile compounds during the flowering process were determined and gas chromatography-mass spectrometry results showed that monoterpenes, such as linalool and ocimene, accounted for the highest proportion, contributing to the floral scent of sweet osmanthus in two cultivars. In addition, correlation analysis revealed the expression patterns of OfWRKYs (OfWRKY7/19/36/139) were each correlated with distinct monoterpenes (linalool, linalool derivatives, ocimene and ocimene derivatives). Subcellular localization analysis showed that p35S::GFP-OfWRKY7/38/95/139 were localized in the nucleus and OfWRKY139 had very strong transactivation activity. Collectively, the results indicated potential roles of OfWRKY139 and OfWRKYs with plant zinc cluster domain in regulating synthesis of aromatic compounds in sweet osmanthus, laying the foundation for use of OfWRKYs to improve the aroma of ornamental plants.

Genome-wide identification and expression of GRAS gene family members in cassava

BACKGROUND

Cassava is highly tolerant to stressful conditions, especially drought stress conditions; however, the mechanisms underlying this tolerance are poorly understood. The GRAS gene family is a large family of transcription factors that are involved in regulating the growth, development, and stress responses of plants. Currently, GRAS transcription factors have not been systematically studied in cassava, which is the sixth most important crop in the world.

RESULTS

Seventy-seven MeGRAS genes were identified from the cassava genome database. Phylogenetic analysis revealed that the MeGRAS proteins could be divided into 14 subfamilies. The gene structure and motif compositions of the proteins were considerably conserved within the same subfamily. Duplication events, particularly segmental duplication, were identified as the main driving force for GRAS gene expansion in cassava. Global expression analysis revealed that MeGRAS genes exhibited similar or distinct expression profiles within different tissues among different varieties. Moreover, qRT-PCR analysis revealed the expression patterns of MeGRAS genes in response to abiotic stress (drought, salt, cold, and H₂O₂), and the results suggest that these genes may have multiple functions.

CONCLUSION

This study is the first to provide comprehensive information on GRAS gene family members in cassava. The data will increase our understanding of both the molecular basis and the effects of GRAS genes. In addition, the results will contribute further to identifying the responses to various environmental conditions and provide insights into the potential functions of GRAS genes.

The Effects of Ultraviolet A/B Treatments on Anthocyanin Accumulation and Gene Expression in Dark-Purple Tea Cultivar 'Ziyan' (Camellia sinensis)

‘Ziyan’ is a novel anthocyanin-rich tea cultivar with dark purple young shoots. However, how its anthocyanin accumulation is affected by environmental factors, such as ultraviolet (UV), remains unclear. In this study, we observed that UV light treatments stimulated anthocyanin accumulation in ‘Ziyan’ leaves, and we further analyzed the underlying mechanisms at gene expression and enzyme activity levels. In addition, the catechins and chlorophyll contents of young shoots under different light treatments were also changed. The results showed that the contents of total anthocyanins and three major anthocyanin molecules, i.e., delphinidin, cyanidin, and pelargonidin, were significantly higher in leaves under UV-A, UV-B, and UV-AB treatments than those under white light treatment alone. However, the total catechins and chlorophyll contents in these purple tea plant leaves displayed the opposite trends. The anthocyanin content was the highest under UV-A treatment, which was higher by about 66% than control. Compared with the white light treatment alone, the enzyme activities of chalcone synthase (CHS), flavonoid 3′,5′-hydroxylase (F3′5′H), and anthocyanidin synthase (ANS) under UV treatments increased significantly, whereas the leucoanthocyanidin reductase (LAR) and anthocyanidin reductase (ANR) activities reduced. There was no significant difference in dihydroflavonol 4-reductase (DFR) activity under all treatments. Comparative transcriptome analyses unveiled that there were 565 differentially expressed genes (DEGs) of 29,648 genes in three pair-wise comparisons (white light versus UV-A, W vs. UV-A; white light versus UV-B, W vs. UV-A; white light versus UV-AB, W vs. UV-AB). The structural genes in anthocyanin pathway such as flavanone 3-hydroxylase (F3H), F3′5′H, DFR, and ANS, and regulatory gene TT8 were upregulated under UV-A treatment; F3′5′H, DFR, ANS, and UFGT and regulatory genes EGL1 and TT2 were upregulated under UV-AB treatment. However, most structural genes involved in phenylpropanoid and flavonoid pathways were downregulated under UV-B treatment compared with control. The expression of LAR and ANR were repressed in all UV treatments. Our results indicated that UV-A and UV-B radiations can induce anthocyanin accumulation in tea plant ‘Ziyan’ by upregulating the structural and regulatory genes involved in anthocyanin biosynthesis. In addition, UV radiation repressed the expression levels of LAR, ANR, and FLS, resulting in reduced ANR activity and a metabolic flux shift toward anthocyanin biosynthesis.

Differentially expressed bZIP transcription factors confer multi-tolerances in Gossypium hirsutum L

Basic leucine zipper (bZIP) transcription factor plays an important role in various biological processes, such as response to biotic and abiotic stresses. In this study we performed a systematic investigation and analysis of bZIP gene family in Gossypium hirsutum to predict their functions in response to different abiotic stresses. A total of 207 bZIP genes were identified from Gossypium hirsutum genome and classified into 13 subfamilies through phylogenetic analysis, which was testified by the analysis of conserved motifs and exon-intron structures. Annotation of GHbZIPs was performed based on well-studied Arabidopsis bZIPs to speculate the gene function. RNA-seq analysis was conducted to identify the co-expressed and differentially expressed bZIPs under cold, heat, salt and PEG treatments. Promoter analysis and interaction network of GHbZIP proteins demonstrated that ABA-activated signaling pathway was pivotal in the regulation of GHbZIPs, and GHbZIPs involved in ER stress were supposed to function through interaction with other GHbZIPs and ABA pathway. Cis-elements in the upstream and downstream of GHbZIPs interaction network were also discussed. These findings provided us with clues about functions of bZIP in Gossypium hirsutum.

Selection of appropriate reference genes for quantitative real-time reverse transcription PCR in Betula platyphylla under salt and osmotic stress conditions

Selecting appropriate reference genes is vital to normalize gene expression analysis in birch (Betula platyphylla) under different abiotic stress conditions using quantitative real-time reverse transcription PCR (qRT-PCR). In this study, 11 candidate birch reference genes (ACT, TUA, TUB, TEF, 18S rRNA, EF1α, GAPDH, UBC, YLS8, SAND, and CDPK) were selected to evaluate the stability of their expression in different tissues and under different abiotic stress conditions. Three statistical algorithms (GeNorm, NormFinder, and BestKeeper) were used to analyze the stability of the 11 candidate reference genes to identify the most appropriate one. The results indicated that EF-1α was the most stable reference gene in different birch tissues, ACT was the most stable reference gene for normal conditions, ACT and TEF were the most stable reference genes for salt stress treatment, TUB was the most stable reference gene for osmotic stress treatment, and ACT was the most appropriate choice in all samples of birch. In conclusion, the most appropriate reference genes varied among different experimental conditions. However, in this study, ACT was the optimum reference gene in all experimental groups, except in the different tissues group. GAPDH was the least stable candidate reference gene in all experimental conditions. In addition, three stress-induced genes (BpGRAS1, BpGRAS16, and BpGRAS19) were chosen to verify the stability of the selected reference genes in different tissues and under salt stress. This study laid the foundation for the selection of appropriate reference gene(s) for future gene expression pattern studies in birch.

Identification of CBF Transcription Factors in Tea Plants and a Survey of Potential CBF Target Genes under Low Temperature

C-repeat binding factors (CBFs) are key signaling genes that can be rapidly induced by cold and bind to the C-repeat/dehydration-responsive motif (CRT/DRE) in the promoter region of the downstream cold-responsive (COR) genes, which play a vital role in the plant response to low temperature. However, the CBF family in tea plants has not yet been elucidated, and the possible target genes regulated by this family under low temperature are still unclear. In this study, we identified five CsCBF family genes in the tea plant genome and analyzed their phylogenetic tree, conserved domains and motifs, and cis-elements. These results indicate that CsCBF3 may be unique in the CsCBF family. This is further supported by our findings from the low-temperature treatment: all the CsCBF genes except CsCBF3 were significantly induced after treatment at 4 °C. The expression profiles of eight tea plant tissues showed that CsCBFs were mainly expressed in winter mature leaves, roots and fruits. Furthermore, 685 potential target genes were identified by transcriptome data and CRT/DRE element information. These target genes play a functional role under the low temperatures of winter through multiple pathways, including carbohydrate metabolism, lipid metabolism, cell wall modification, circadian rhythm, calcium signaling, transcriptional cascade, and hormone signaling pathways. Our findings will further the understanding of the stress regulatory network of CsCBFs in tea plants.

Genome-wide survey and expression analyses of the GRAS gene family in Brassica napus reveals their roles in root development and stress response

Genome-wide identification, classification, expression analyses, and functional characterization of GRAS genes in oil crop, Brassica napus, indicate their importance in root development and stress response. GRAS proteins are a plant-specific transcription factor gene family involved in tissues development and stress response. We classified 87 putative GRAS genes in the Brassica napus genome (BnGRASs) into 13 subfamilies by phylogenetic analysis. The C-terminal GRAS domains of Brassica napus (B. napus) proteins were less conserved among subfamilies, but were conserved within each subfamily. A series of analyses revealed that 89.7% of the BnGRASs did not have intron insertions, and 24 specific-motifs were found at the N-terminal. A highly conserved microRNA 171 (miRNA171) target was observed specifically in the HAM subfamily across land plants. A total of 868 pairs of interaction proteins were predicted, the primary of which were transcription factors involved in transcriptional regulation and signal transduction. Integrated comparative analysis of GRAS genes across 26 species of algae, mosses, ferns, gymnosperms, and angiosperms revealed that this gene family originated in early mosses and was classified into 19 subfamilies, 14 of which may have originated prior to bryophyte evolution. RNA-Seq analysis demonstrated that most BnGRASs were widely expressed in different tissues/organs at different stages in B. napus, and 24 BnGRASs were highly/specifically expressed in roots. Results from a qRT-PCR analysis suggested that two BnGRASs belonging to SCR and LISCL subfamilies potentially have important roles in the stress response of roots.

Genome-wide identification and expression analyses of the LEA protein gene family in tea plant reveal their involvement in seed development and abiotic stress responses

Late embryogenesis abundant (LEA) proteins are widely known to be present in higher plants and are believed to play important functional roles in embryonic development and abiotic stress responses. However, there is a current lack of systematic analyses on the LEA protein gene family in tea plant. In this study, a total of 48 LEA genes were identified using Hidden Markov Model profiles in C. sinensis, and were classified into seven distinct groups based on their conserved domains and phylogenetic relationships. Genes in the CsLEA_2 group were found to be the most abundant. Gene expression analyses revealed that all the identified CsLEA genes were expressed in at least one tissue, and most had higher expression levels in the root or seed relative to other tested tissues. Nearly all the CsLEA genes were found to be involved in seed development, and thirty-nine might play an important role in tea seed maturation concurrent with dehydration. However, only sixteen CsLEA genes were involved in seed desiccation, and furthermore, most were suppressed. Additionally, forty-six CsLEA genes could be induced by at least one of the tested stress treatments, and they were especially sensitive to high temperature stress. Furthermore, it was found that eleven CsLEA genes were involved in tea plant in response to all tested abiotic stresses. Overall, this study provides new insights into the formation of CsLEA gene family members and improves our understanding on the potential roles of these genes in normal development processes and abiotic stress responses in tea plant, particularly during seed development and desiccation. These results are beneficial for future functional studies of CsLEA genes that will help preserve the recalcitrant tea seeds for a long time and genetically improve tea plant.

Identification and Expression Analysis of GRAS Transcription Factors to Elucidate Candidate Genes Related to Stolons, Fruit Ripening and Abiotic Stresses in Woodland Strawberry (Fragaria vesca)

The cultivated strawberry (Fragaria × ananassa), an allo-octoploid with non-climacteric fleshy fruits, is a popular Rosaceae horticultural crop worldwide that is mainly propagated via stolons during cultivation. Woodland strawberry (Fragaria vesca), one of the four diploid progenitor species of cultivated strawberry, is widely used as a model plant in the study of Rosaceae fruit trees, non-climacteric fruits and stolons. One GRAS transcription factor has been shown to regulate stolon formation; the other GRAS proteins in woodland strawberry remain unknown. In this study, we identified 54 FveGRAS proteins in woodland strawberry, and divided them into 14 subfamilies. Conserved motif analysis revealed that the motif composition of FveGRAS proteins was conserved within each subfamily, but diverged widely among subfamilies. We found 56 orthologous pairs of GRAS proteins between woodland strawberry and Arabidopsis thaliana, 47 orthologous pairs between woodland strawberry and rice and 92 paralogous pairs within woodland strawberry. The expression patterns of FveGRAS genes in various organs and tissues, and changes therein under cold, heat and GA₃ treatments, were characterized using transcriptomic analysis. The results showed that 34 FveGRAS genes were expressed with different degrees in at least four organs, including stolons; only a few genes displayed organ-specific expression. The expression levels of 16 genes decreased, while that of four genes increased during fruit ripening; FveGRAS54 showed the largest increase in expression. Under cold, heat and GA₃ treatments, around half of the FveGRAS genes displayed increased or decreased expression to some extent, suggesting differing functions of these FveGRAS genes in the responses to cold, heat and GAs. This study provides insight into the potential functions of FveGRAS genes in woodland strawberry. A few FveGRAS genes were identified as candidate genes for further study, in terms of their functions in stolon formation, fruit ripening and abiotic stresses.

Genome-wide Identification, Classification, Expression and Duplication Analysis of GRAS Family Genes in Juglans regia L

Fifty-two GRAS genes are identified in walnut genome. Based on the evolutionary relationship and motif analysis, the walnut GRAS gene family was divided into eight subfamilies, and the sequence features analysis of JrGRAS proteins showed that the JrGRAS protein sequences were both conserved and altered during the evolutionary process. Gene duplication analysis indicated that seven GRAS genes in walnut have orthologous genes in other species, and five of them occurred duplicated events in walnut genome. Expression pattern analysis of the GRAS family genes in walnut showed that two JrGRAS genes (JrCIGRa-b and JrSCL28a) were differentially expressed between flower bud and leaf bud (p < 0.01), and two JrGRAS genes (JrCIGRa-b and JrSCL13b-d) were differentially expressed between the different development stages of flower buds transition (p < 0.01), besides, three hub genes (JrGAIa, JrSCL3f and JrSHRc) were identified by co-expression analysis, which suggested these GRAS genes may play an important role in regulating the development of apical meristem in walnut. This study laid a foundation for further understanding of the function of GRAS family genes in walnut.

Identification of the GRAS gene family in the Brassica juncea genome provides insight into its role in stem swelling in stem mustard

GRAS transcription factors are known to play important roles in plant signal transduction and development. A comprehensive study was conducted to explore the GRAS family in the Brassica juncea genome. A total of 88 GRAS genes were identified which were categorized into nine groups according to the phylogenetic analysis. Gene structure analysis showed a high group-specificity, which corroborated the gene grouping results. The chromosome distribution and sequence analysis suggested that gene duplication events are vital for the expansion of GRAS genes in the B. juncea genome. The changes in evolution rates and amino acid properties among groups might be responsible for their functional divergence. Interaction networks and cis-regulatory elements were analyzed including DELLA and eight interaction proteins (including four GID1, two SLY1, and two PIF3 proteins) that are primarily involved in light and hormone signaling. To understand their regulatory role in growth and development, the expression profiles of BjuGRASs and interaction genes were examined based on transcriptome data and qRT-PCR, and selected genes (BjuGRAS3, 5, 7, 8, 10, BjuB006276, BjuB037910, and BjuA021658) had distinct temporal expression patterns during stem swelling, indicating that they possessed diverse regulatory functions during the developmental process. These results contribute to our understanding on the GRAS gene family and provide the basis for further investigations on the evolution and functional characterization of GRAS genes.

Genome-wide identification, characterization, interaction network and expression profile of GRAS gene family in sweet orange (Citrus sinensis)

GRAS genes are suggested to be grouped into plant-specific transcriptional regulatory families that have been reported to participate in multiple processes, including plant development, phytohormone signaling, the formation of symbiotic relationships, and response to environmental signals. GRAS genes have been characterized in a number of plant species, but little is known about this gene family in Citrus sinensis. In this study, we identified a total of 50 GRAS genes and characterized the gene structures, conserved motifs, genome localizations and cis-elements within their promoter regions. According to their structural and phylogenetic features, the identified sweet orange GRAS members were divided into 11 subgroups, of which subfamily CsGRAS34 was sweet orange-specific. Based on publicly available RNA-seq data generated from callus, flower, leaf and fruit in sweet orange, we found that some sweet orange GRAS genes exhibited tissue-specific expression patterning. Three of the six members of subfamily AtSHR, particularly CsGRAS9, and two of the six members of subfamily AtPAT1 were preferentially expressed in leaf. Moreover, protein-protein interactions with CsGRAS were predicted. Gene expression analysis was performed under conditions of phosphate deficiency, and GA3 and NaCl treatment to identify the potential functions of GRAS members in regulating stress and hormone responses. This study provides the first comprehensive understanding of the GRAS gene family in the sweet orange genome. As such, the study generates valuable information for further gene function analysis and identifying candidate genes to improve abiotic stress tolerance in citrus plants.

Genomic and transcriptomic analyses of HD-Zip family transcription factors and their responses to abiotic stress in tea plant (Camellia sinensis)

Tea plant (Camellia sinensis (L.) O. Kuntze) is a perennial evergreen woody plant, and its leaves contain various beneficial ingredients and have healthy efficacy. HD-Zip (homeodomain-leucine zipper) transcription factors (TFs) are widely distributed in plants and play an important role in plant growth and environmental response. To date, knowledge on HD-Zip gene family in tea plant is still limited. In this study, 33 HD-Zip TFs were selected based on the genomic and transcriptomic databases of tea plant. The conserved domains and common motifs of these TFs were predicted and analyzed. These 33 Cshdz TFs were divided into four groups (HD-Zip I, HD-Zip II, HD-Zip III, and HD-Zip IV). The interaction network of the HD-Zip proteins of tea plant was established based on the data of Arabidopsis. In addition, the expression levels of these Cshdz genes in tea plant cv. 'Longjing43' were detected and analyzed under five abiotic stress treatments. Results showed that the different expression profiles of Cshdz genes were associated with different abiotic stress treatments. Our findings suggested a potential relationship between the resistance of tea plant and its Cshdz genes.