Evolutionary Analyses of GRAS Transcription Factors in Angiosperms

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

A journey into the world of small RNAs in the arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal (AM) symbiosis is a mutualistic interaction between fungi and most land plants that is underpinned by a bidirectional exchange of nutrients. AM development is a tightly regulated process that encompasses molecular communication for reciprocal recognition, fungal accommodation in root tissues and activation of symbiotic function. As such, a complex network of transcriptional regulation and molecular signaling underlies the cellular and metabolic reprogramming of host cells upon AM fungal colonization. In addition to transcription factors, small RNAs (sRNAs) are emerging as important regulators embedded in the gene network that orchestrates AM development. In addition to controlling cell-autonomous processes, plant sRNAs also function as mobile signals capable of moving to different organs and even to different plants or organisms that interact with plants. AM fungi also produce sRNAs; however, their function in the AM symbiosis remains largely unknown. Here, we discuss the contribution of host sRNAs in the development of AM symbiosis by considering their role in the transcriptional reprogramming of AM fungal colonized cells. We also describe the characteristics of AM fungal-derived sRNAs and emerging evidence for the bidirectional transfer of functional sRNAs between the two partners to mutually modulate gene expression and control the symbiosis.

GRAS gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses

BACKGROUND

The GRAS transcription factor family plays a crucial role in various biological processes in different plants, such as tissue development, fruit maturation, and environmental stress. However, the GRAS family in rye has not been systematically analyzed yet.

RESULTS

In this study, 67 GRAS genes in S. cereale were identified and named based on the chromosomal location. The gene structures, conserved motifs, cis-acting elements, gene replications, and expression patterns were further analyzed. These 67 ScGRAS members are divided into 13 subfamilies. All members include the LHR I, VHIID, LHR II, PFYRE, and SAW domains, and some nonpolar hydrophobic amino acid residues may undergo cross-substitution in the VHIID region. Interested, tandem duplications may have a more important contribution, which distinguishes them from other monocotyledonous plants. To further investigate the evolutionary relationship of the GRAS family, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. The response characteristics of 19 ScGRAS members from different subfamilies to different tissues, grains at filling stages, and different abiotic stresses of rye were systematically analyzed. Paclobutrazol, a triazole-based plant growth regulator, controls plant tissue and grain development by inhibiting gibberellic acid (GA) biosynthesis through the regulation of DELLA proteins. Exogenous spraying of paclobutrazol significantly reduced the plant height but was beneficial for increasing the weight of 1000 grains of rye. Treatment with paclobutrazol, significantly reduced gibberellin levels in grain in the filling period, caused significant alteration in the expression of the DELLA subfamily gene members. Furthermore, our findings with respect to genes, ScGRAS46 and ScGRAS60, suggest that these two family members could be further used for functional characterization studies in basic research and in breeding programmes for crop improvement.

CONCLUSIONS

We identified 67 ScGRAS genes in rye and further analysed the evolution and expression patterns of the encoded proteins. This study will be helpful for further analysing the functional characteristics of ScGRAS genes.

Mapping yield and yield-related traits using diverse common bean germplasm

Common bean (bean) is one of the most important legume crops, and mapping genes for yield and yield-related traits is essential for its improvement. However, yield is a complex trait that is typically controlled by many loci in crop genomes. The objective of this research was to identify regions in the bean genome associated with yield and a number of yield-related traits using a collection of 121 diverse bean genotypes with different yields. The beans were evaluated in replicated trials at two locations, over two years. Significant variation among genotypes was identified for all traits analyzed in the four environments. The collection was genotyped with the BARCBean6K_3 chip (5,398 SNPs), two yield/antiyield gene-based markers, and seven markers previously associated with resistance to common bacterial blight (CBB), including a Niemann-Pick polymorphism (NPP) gene-based marker. Over 90% of the single-nucleotide polymorphisms (SNPs) were polymorphic and separated the panel into two main groups of small-seeded and large-seeded beans, reflecting their Mesoamerican and Andean origins. Thirty-nine significant marker-trait associations (MTAs) were identified between 31 SNPs and 15 analyzed traits on all 11 bean chromosomes. Some of these MTAs confirmed genome regions previously associated with the yield and yield-related traits in bean, but a number of associations were not reported previously, especially those with derived traits. Over 600 candidate genes with different functional annotations were identified for the analyzed traits in the 200-Kb region centered on significant SNPs. Fourteen SNPs were identified within the gene model sequences, and five additional SNPs significantly associated with five different traits were located at less than 0.6 Kb from the candidate genes. The work confirmed associations between two yield/antiyield gene-based markers (AYD1m and AYD2m) on chromosome Pv09 with yield and identified their association with a number of yield-related traits, including seed weight. The results also confirmed the usefulness of the NPP marker in screening for CBB resistance. Since disease resistance and yield measurements are environmentally dependent and labor-intensive, the three gene-based markers (CBB- and two yield-related) and quantitative trait loci (QTL) that were validated in this work may be useful tools for simplifying and accelerating the selection of high-yielding and CBB-resistant bean cultivars.

DWARF AND LOW-TILLERING 2 functions in brassinosteroid signaling and controls plant architecture and grain size in rice

Brassinosteroids (BRs) are a class of steroid phytohormones that control various aspects of plant growth and development. Several transcriptional factors (TFs) have been suggested to play roles in BR signaling. However, their possible relationship remains largely unknown. Here, we identified a rice mutant dwarf and low-tillering 2 (dlt2) with altered plant architecture, increased grain width, and reduced BR sensitivity. DLT2 encodes a GIBBERELLIN INSENSITIVE (GAI)-REPRESSOR OF GAI (RGA)-SCARECROW (GRAS) TF that is mainly localized in the nucleus and has weak transcriptional activity. Our further genetic and biochemical analyses indicate that DLT2 interacts with two BR-signaling-related TFs, DLT and BRASSINAZOLE-RESISTANT 1, and probably modulates their transcriptional activity. These findings imply that DLT2 is implicated in a potentially transcriptional complex that mediates BR signaling and rice development and suggests that DLT2 could be a potential target for improving rice architecture and grain morphology. This work also sheds light on the role of rice GRAS members in regulating numerous developmental processes.

Genome-wide identification, expression analysis, and functional study of the bZIP transcription factor family and its response to hormone treatments in pea (Pisum sativum L.)

BACKGROUND

Basic leucine zipper (bZIP) protein is a plant-specific transcription factor involved in various biological processes, including light signaling, seed maturation, flower development, cell elongation, seed accumulation protein, and abiotic and biological stress responses. However, little is known about the pea bZIP family.

RESULTS

In this study, we identified 87 bZIP genes in pea, named PsbZIP1 ~ PsbZIP87, via homology analysis using Arabidopsis. The genes were divided into 12 subfamilies and distributed unevenly in 7 pea chromosomes. PsbZIPs in the same subfamily contained similar intron/exon organization and motif composition. 1 tandem repeat event and 12 segmental duplication events regulated the expansion of the PsbZIP gene family. To better understand the evolution of the PsbZIP gene family, we conducted collinearity analysis using Arabidopsis thaliana, Oryza sativa Japonica, Fagopyrum tataricum, Solanum lycopersicum, Vitis vinifera, and Brachypodium distachyon as the related species of pea. In addition, interactions between PsbZIP proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of PsbZIP expression was complex. We also evaluated the expression patterns of bZIP genes in different tissues and at different fruit development stages, all while subjecting them to five hormonal treatments.

CONCLUSION

These results provide a deeper understanding of PsbZIP gene family evolution and resources for the molecular breeding of pea. The findings suggested that PsbZIP genes, specifically PSbZIP49, play key roles in the development of peas and their response to various hormones.

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.

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, phylogenetic and expression pattern analysis of HSF family genes in the Rye (Secale cereale L.)

BACKGROUND

Heat shock factor (HSF), a typical class of transcription factors in plants, has played an essential role in plant growth and developmental stages, signal transduction, and response to biotic and abiotic stresses. The HSF genes families has been identified and characterized in many species through leveraging whole genome sequencing (WGS). However, the identification and systematic analysis of HSF family genes in Rye is limited.

RESULTS

In this study, 31 HSF genes were identified in Rye, which were unevenly distributed on seven chromosomes. Based on the homology of A. thaliana, we analyzed the number of conserved domains and gene structures of ScHSF genes that were classified into seven subfamilies. To better understand the developmental mechanisms of ScHSF family during evolution, we selected one monocotyledon (Arabidopsis thaliana) and five (Triticum aestivum L., Hordeum vulgare L., Oryza sativa L., Zea mays L., and Aegilops tauschii Coss.) specific representative dicotyledons associated with Rye for comparative homology mapping. The results showed that fragment replication events modulated the expansion of the ScHSF genes family. In addition, interactions between ScHSF proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of ScHSF expression was complex. A total of 15 representative genes were targeted from seven subfamilies to characterize their gene expression responses in different tissues, fruit developmental stages, three hormones, and six different abiotic stresses.

CONCLUSIONS

This study demonstrated that ScHSF genes, especially ScHSF1 and ScHSF3, played a key role in Rye development and its response to various hormones and abiotic stresses. These results provided new insights into the evolution of HSF genes in Rye, which could help the success of molecular breeding in Rye.

Genome-Wide Identification, Expression and Stress Analysis of the GRAS Gene Family in Phoebe bournei

GRAS genes are important transcriptional regulators in plants that govern plant growth and development through enhancing plant hormones, biosynthesis, and signaling pathways. Drought and other abiotic factors may influence the defenses and growth of Phoebe bournei, which is a superb timber source for the construction industry and building exquisite furniture. Although genome-wide identification of the GRAS gene family has been completed in many species, that of most woody plants, particularly P. bournei, has not yet begun. We performed a genome-wide investigation of 56 PbGRAS genes, which are unequally distributed across 12 chromosomes. They are divided into nine subclades. Furthermore, these 56 PbGRAS genes have a substantial number of components related to abiotic stress responses or phytohormone transmission. Analysis using qRT-PCR showed that the expression of four PbGRAS genes, namely PbGRAS7, PbGRAS10, PbGRAS14 and PbGRAS16, was differentially increased in response to drought, salt and temperature stresses, respectively. We hypothesize that they may help P. bournei to successfully resist harsh environmental disturbances. In this work, we conducted a comprehensive survey of the GRAS gene family in P. bournei plants, and the results provide an extensive and preliminary resource for further clarification of the molecular mechanisms of the GRAS gene family in P. bournei in response to abiotic stresses and forestry improvement.

Functional Study of PgGRAS68-01 Gene Involved in the Regulation of Ginsenoside Biosynthesis in Panax ginseng

Ginseng (Panax ginseng C. A. Meyer) is a perennial herb from the genus Panax in the family Araliaceae. It is famous in China and abroad. The biosynthesis of ginsenosides is controlled by structural genes and regulated by transcription factors. GRAS transcription factors are widely found in plants. They can be used as tools to modify plant metabolic pathways by interacting with promoters or regulatory elements of target genes to regulate the expression of target genes, thereby activating the synergistic interaction of multiple genes in metabolic pathways and effectively improving the accumulation of secondary metabolites. However, there are no reports on the involvement of the GRAS gene family in ginsenoside biosynthesis. In this study, the GRAS gene family was located on chromosome 24 pairs in ginseng. Tandem replication and fragment replication also played a key role in the expansion of the GRAS gene family. The PgGRAS68-01 gene closely related to ginsenoside biosynthesis was screened out, and the sequence and expression pattern of the gene were analyzed. The results showed that the expression of PgGRAS68-01 gene was spatio-temporal specific. The full-length sequence of PgGRAS68-01 gene was cloned, and the overexpression vector pBI121-PgGRAS68-01 was constructed. The ginseng seedlings were transformed by Agrobacterium rhifaciens-mediated method. The saponin content in the single root of positive hair root was detected, and the inhibitory role of PgGRAS68-01 in ginsenoside synthesis is reported.

Genome-wide identification, evolution and transcriptome analysis of GRAS gene family in Chinese chestnut (Castanea mollissima)

GRAS transcription factors play an important role in regulating various biological processes in plant growth and development. However, their characterization and potential function are still vague in Chinese chestnut (Castanea mollissima), an important nut with rich nutrition and high economic value. In this study, 48 CmGRAS genes were identified in Chinese chestnut genome and phylogenetic analysis divided CmGRAS genes into nine subfamilies, and each of them has distinct conserved structure domain and features. Genomic organization revealed that CmGRAS tend to have a representative GRAS domain and fewer introns. Tandem duplication had the greatest contribution for the CmGRAS expansion based on the comparative genome analysis, and CmGRAS genes experienced strong purifying selection pressure based on the Ka/Ks. Gene expression analysis revealed some CmGRAS members with potential functions in bud development and ovule fertility. CmGRAS genes with more homologous relationships with reference species had more cis-acting elements and higher expression levels. Notably, the lack of DELLA domain in members of the DELLA subfamily may cause de functionalization, and the differences between the three-dimensional structures of them were exhibited. This comprehensive study provides theoretical and practical basis for future research on the evolution and function of GRAS gene family.

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.

Genome-Wide Identification, Characterization, and Expression Analysis of GRAS Gene Family in Ginger (Zingiber officinale Roscoe)

GRAS family proteins are one of the most abundant transcription factors in plants; they play crucial roles in plant development, metabolism, and biotic- and abiotic-stress responses. The GRAS family has been identified and functionally characterized in some plant species. However, this family in ginger (Zingiber officinale Roscoe), a medicinal crop and non-prescription drug, remains unknown to date. In the present study, 66 GRAS genes were identified by searching the complete genome sequence of ginger. The GRAS family is divided into nine subfamilies based on the phylogenetic analyses. The GRAS genes are distributed unevenly across 11 chromosomes. By analyzing the gene structure and motif distribution of GRAS members in ginger, we found that the GRAS genes have more than one cis-acting element. Chromosomal location and duplication analysis indicated that whole-genome duplication, tandem duplication, and segmental duplication may be responsible for the expansion of the GRAS family in ginger. The expression levels of GRAS family genes are different in ginger roots and stems, indicating that these genes may have an impact on ginger development. In addition, the GRAS genes in ginger showed extensive expression patterns under different abiotic stresses, suggesting that they may play important roles in the stress response. Our study provides a comprehensive analysis of GRAS members in ginger for the first time, which will help to better explore the function of GRAS genes in the regulation of tissue development and response to stress in ginger.

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 Classification and Phylogenetic Analyses of the GDSL-Type Esterase/Lipase (GELP) Family in Flowering Plants

GDSL-type esterase/lipase (GELP) enzymes have key functions in plants, such as developmental processes, anther and pollen development, and responses to biotic and abiotic stresses. Genes that encode GELP belong to a complex and large gene family, ranging from tens to more than hundreds of members per plant species. To facilitate functional transfer between them, we conducted a genome-wide classification of GELP in 46 plant species. First, we applied an iterative phylogenetic method using a selected set of representative angiosperm genomes (three monocots and five dicots) and identified 10 main clusters, subdivided into 44 orthogroups (OGs). An expert curation for gene structures, orthogroup composition, and functional annotation was made based on a literature review. Then, using the HMM profiles as seeds, we expanded the classification to 46 plant species. Our results revealed the variable evolutionary dynamics between OGs in which some expanded, mostly through tandem duplications, while others were maintained as single copies. Among these, dicot-specific clusters and specific amplifications in monocots and wheat were characterized. This approach, by combining manual curation and automatic identification, was effective in characterizing a large gene family, allowing the establishment of a classification framework for gene function transfer and a better understanding of the evolutionary history of GELP.

GRAS transcription factors emerging regulator in plants growth, development, and multiple stresses

GRAS transcription factors play multifunctional roles in plant growth, development, and resistance to various biotic and abiotic stresses. The structural and functional features of GRAS TFs have been unveiled in the last two decades. A typical GRAS protein contained a C-terminal GRAS domain with a highly variable N-terminal region. Studies on these TFs increase in numbers and are reported to be involved in various important developmental processes such as flowering, root formation, and stress responses. The GRAS TFs and hormone signaling crosstalk can be implicated in plant development and to stress responses. There are relatively few reports about GRAS TFs roles in plants, and no related reviews have been published. In this review, we summarized the features of GRAS TFs, their targets, and the roles these GRAS TFs playing in plant development and multiple stresses.

Genome-wide Identification and Characterization of the GRAS Transcription Factors in Garlic (Allium sativum L.)

GRAS transcription factors play crucial roles in plant growth and development and have been widely explored in many plant species. Garlic (Allium sativum L.) is an important crop owing to its edible and medicinal properties. However, no GRAS transcription factors have been identified in this crop. In this study, 46 garlic GRAS genes were identified and assigned to 16 subfamilies using the GRAS members of Arabidopsis thaliana, Oryza sativa, and Amborella trichopoda as reference queries. Expression analysis revealed that garlic GRAS genes showed distinct differences in various garlic tissues, as well as during different growth stages of the bulbs. Five of these 46 genes were identified as DELLA-like protein-encoding genes and three of which, Asa2G00237.1/Asa2G00240.1 and Asa4G02090.1, responded to exogenous GA3 treatment, and showed a significant association between their transcription abundance and bulb traits in 102 garlic accessions, thereby indicating their role in regulating the growth of garlic bulbs. These results will lay a useful foundation for further investigation of the biological functions of GRAS genes and guiding the genetic breeding of garlic in the future.

Expression and roles of GRAS gene family in plant growth, signal transduction, biotic and abiotic stress resistance and symbiosis formation-a review

The GRAS (derived from GAI, RGA and SCR) gene family consists of plant-specific genes, works as a transcriptional regulator and plays a key part in the regulation of plant growth and development. The past decade has witnessed significant progress in understanding and advances on GRAS transcription factors in various plants. A notable concern is to what extent the mechanisms found in plants, particularly crops, are shared by other species, and what other characteristics are dependent on GRAS transcription factor (TFS)-mediated gene expression. GRAS are involved in many processes that are intimately linked to plant growth regulation. However, GRAS also perform additional roles against environmental stresses, allowing plants to function more efficiently. GRAS increase plant growth and development by improving several physiological processes, such as phytohormone, biosynthetic and signalling pathways. Furthermore, the GRAS gene family plays an important role in response to abiotic stresses, e.g. photooxidative stress. Moreover, evidence shows the involvement of GRAS in arbuscule development during plant-mycorrhiza associations. In this review, the diverse roles of GRAS in plant systems are highlighted that could be useful in enhancing crop productivity through genetic modification, especially of crops. This is the first review to report the role and function of the GRAS gene family in plant systems. Furthermore, a large number of studies are reviewed, and several limitations and research gaps identified that must be addressed in future studies.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

Multifarious and Interactive Roles of GRAS Transcription Factors During Arbuscular Mycorrhiza Development

Arbuscular mycorrhiza (AM) is a mutualistic symbiotic interaction between plant roots and AM fungi (AMF). This interaction is highly beneficial for plant growth, development and fitness, which has made AM symbiosis the focus of basic and applied research aimed at increasing plant productivity through sustainable agricultural practices. The creation of AM requires host root cells to undergo significant structural and functional modifications. Numerous studies of mycorrhizal plants have shown that extensive transcriptional changes are induced in the host during all stages of colonization. Advances have recently been made in identifying several plant transcription factors (TFs) that play a pivotal role in the transcriptional regulation of AM development, particularly those belonging to the GRAS TF family. There is now sufficient experimental evidence to suggest that GRAS TFs are capable to establish intra and interspecific interactions, forming a transcriptional regulatory complex that controls essential processes in the AM symbiosis. In this minireview, we discuss the integrative role of GRAS TFs in the regulation of the complex genetic re-programming determining AM symbiotic interactions. Particularly, research being done shows the relevance of GRAS TFs in the morphological and developmental changes required for the formation and turnover of arbuscules, the fungal structures where the bidirectional nutrient translocation occurs.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Genome-wide identification and expression analysis of the GRAS gene family in response to drought stress in chickpea (Cicer arietinum L.)

The GRAS (gibberellic acid insensitive, repressor of GAI and scarecrow) transcription factors (TFs) regulate diverse biological processes involved in plant growth and development. These TFs are also known to regulate gene expression in response to various abiotic stress factors like cold, drought, etc. In chickpea one of the most devastating abiotic stress factors is terminal drought. The GRAS TF family has not been characterized in chickpea (Cicer arietinum L.) until now. In this study, we report 46 GRAS TF genes (CaGRAS genes) in the chickpea genome. The CaGRAS proteins were categorized into nine subfamilies based on their phylogenetic relationship with known GRAS members of Arabidopsis and soybean. The PAT subfamily was the largest consisting of ten CaGRAS members whereas the LAS subfamily was the smallest with only one member. Gene duplication analysis revealed that segmental duplication was the primary reason for the expansion of this gene family within the chickpea genome. The gene expression levels of CaGRAS genes were analysed using two different chickpea varieties contrasting for drought tolerance trait, i.e., ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive). On exposure to drought stress, the two chickpea genotypes, exhibited differential drought response, which was quantified and estimated in terms of differences in leaf relative water content (RWC). The well-watered or control plants of the drought tolerant variety were able to maintain a higher leaf RWC by the end of the drought stress period, whereas the control plants of the drought sensitive variety continued to show a decline in leaf RWC. The two genotypes also differed in their root morphologies, under well-watered and drought stress conditions. The gene expression analysis revealed a potential role of PAT, SCR, SCL3 and SHR GRAS members in the regulation of differential response to drought, in the root tissues, for both the genotypes. CaGRAS 12 (SCR) was identified as a drought-responsive GRAS TF gene, which could serve as a potential candidate gene for utilization in developing chickpea varieties with improved drought tolerance. This study demonstrates the drought-responsive expression of CaGRAS genes in chickpea and also describes the morpho-physiological response of chickpea plants to drought stress conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03104-z.

A novel SCARECROW-LIKE3 transcription factor LjGRAS36 in Lotus japonicus regulates the development of arbuscular mycorrhizal symbiosis

The symbiosis with arbuscular mycorrhizal (AM) fungi improves plants' nutrient uptake. During this process, transcription factors have been highlighted to play crucial roles. Members of the GRAS transcription factor gene family have been reported involved in AM symbiosis, but little is known about SCARECROW-LIKE3 (SCL3) genes belonging to this family in Lotus japonicus. In this study, 67 LjGRAS genes were identified from the L. japonicus genome, seven of which were clustered in the SCL3 group. Three of the seven LjGRAS genes expression levels were upregulated by AM fungal inoculation, and our biochemical results showed that the expression of LjGRAS36 was specifically induced by AM colonization. Functional loss of LjGRAS36 in mutant ljgras36 plants exhibited a significantly reduced mycorrhizal colonization rate and arbuscular size. Transcriptome analysis showed a deficiency of LjGRAS36 led to the dysregulation of the gibberellic acid signal pathway associated with AM symbiosis. Together, this study provides important insights for understanding the important potential function of SCL3 genes in regulating AM symbiotic development.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01161-z.

Genome-wide and comparative phylogenetic analysis of senescence-associated NAC transcription factors in sunflower (Helianthus annuus)

BACKGROUND

Leaf senescence delay impacts positively in grain yield by maintaining the photosynthetic area during the reproductive stage and during grain filling. Therefore a comprehensive understanding of the gene families associated with leaf senescence is essential. NAC transcription factors (TF) form a large plant-specific gene family involved in regulating development, senescence, and responses to biotic and abiotic stresses. The main goal of this work was to identify sunflower NAC TF (HaNAC) and their association with senescence, studying their orthologous to understand possible functional relationships between genes of different species.

RESULTS

To clarify the orthologous relationships, we used an in-depth comparative study of four divergent taxa, in dicots and monocots, with completely sequenced genomes (Arabidopsis thaliana, Vitis vinifera, Musa acuminata and Oryza sativa). These orthologous groups provide a curated resource for large scale protein sequence annotation of NAC TF. From the 151 HaNAC genes detected in the latest version of the sunflower genome, 50 genes were associated with senescence traits. These genes showed significant differential expression in two contrasting lines according to an RNAseq assay. An assessment of overexpressing the Arabidopsis line for HaNAC001 (a gene of the same orthologous group of Arabidopsis thaliana ORE1) revealed that this line displayed a significantly higher number of senescent leaves and a pronounced change in development rate.

CONCLUSIONS

This finding suggests HaNAC001 as an interesting candidate to explore the molecular regulation of senescence in sunflower.

Genome-Wide Identification, Transcript Profiling and Bioinformatic Analyses of GRAS Transcription Factor Genes in Rice

Our group has previously identified the activation of a GRAS transcription factor (TF) gene in the gain-of-function mutant population developed through activation tagging in rice (in an indica rice variety, BPT 5204) that was screened for water use efficiency. This family of GRAS transcription factors has been well known for their diverse roles in gibberellin signaling, light responses, root development, gametogenesis etc. Recent studies indicated their role in biotic and abiotic responses as well. Although this family of TFs received significant attention, not many genes were identified specifically for their roles in mediating stress tolerance in rice. Only OsGRAS23 (here named as OsGRAS22) was reported to code for a TF that induced drought tolerance in rice. In the present study, we have analyzed the expression patterns of rice GRAS TF genes under abiotic (NaCl and ABA treatments) and biotic (leaf samples infected with pathogens, Xanthomonas oryzae pv. oryzae that causes bacterial leaf blight and Rhizoctonia solani that causes sheath blight) stress conditions. In addition, their expression patterns were also analyzed in 13 different developmental stages. We studied their spatio-temporal regulation and correlated them with the in-silico studies. Fully annotated genomic sequences available in rice database have enabled us to study the protein properties, ligand interactions, domain analysis and presence of cis-regulatory elements through the bioinformatic approach. Most of the genes were induced immediately after the onset of stress particularly in the roots of ABA treated plants. OsGRAS39 was found to be a highly expressive gene under sheath blight infection and both abiotic stress treatments while OsGRAS8, OsSHR1 and OsSLR1 were also responsive. Our earlier activation tagging based functional characterization followed by the genome-wide characterization of the GRAS gene family members in the present study clearly show that they are highly appropriate candidate genes for manipulating stress tolerance in rice and other crop plants.

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.

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 scarecrow-like transcription factor SlSCL3 regulates volatile terpene biosynthesis and glandular trichome size in tomato (Solanum lycopersicum)

Tomato (Solanum lycopersicum L.) type VI glandular trichomes that occur on the surface of leaves, stems, young fruits and flowers produce and store a blend of volatile monoterpenes and sesquiterpenes. These compounds play important roles in the interaction with pathogens and herbivorous insects. Although the function of terpene synthases in the biosynthesis of volatile terpenes in tomato has been comprehensively investigated, the deciphering of their transcriptional regulation is only just emerging. We selected transcription factors that are over-expressed in trichomes based on existing transcriptome data and silenced them individually by virus-induced gene silencing. Of these, SlSCL3, a scarecrow-like (SCL) subfamily transcription factor, led to a significant decrease in volatile terpene content and expression of the corresponding terpene synthase genes when its transcription level was downregulated. Overexpression of SlSCL3 dramatically increased both the volatile terpene content and glandular trichome size, whereas its homozygous mutants showed reduced terpene biosynthesis. However, its heterozygous mutants also showed a significantly elevated volatile terpene content and enlarged glandular trichomes, similar to the overexpression plants. SlSCL3 modulates the expression of terpene biosynthetic pathway genes by transcriptional activation, but neither direct protein-DNA binding nor interaction with known regulators was observed. Moreover, transcript levels of the endogenous copy of SlSCL3 were decreased in the overexpression plants but increased in the heterozygous and homozygous mutants, suggesting feedback repression of its own promoter. Taken together, our results provide new insights into the role of SlSCL3 in the complex regulation of volatile terpene biosynthesis and glandular trichome development in tomato.

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, 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.

Molecular Cloning, Transcriptional Profiling, Subcellular Localization, and miRNA-Binding Site Analysis of Six SCL9 Genes in Poplar

The SCL9 subfamily is a key member of the GRAS family that regulates plant development and stress responses. Nevertheless, the functional role of these genes in the growth and development of poplar still unclear. Here, we reported the six SCL9 genes, which were found to be differentially expressed during poplar adventitious root formation. The full-length sequences of PeSCL9 genes of ‘Nanlin895’ poplar (Populus deltoids × Populus euramericana) were cloned by the RACE technique All PeSCL9 genes lacked introns. RT-qPCR revealed that PeSCL9 genes displayed a dynamic expression pattern in the adventitious root of poplar, according to RT-qPCR data. A series of comprehensive genes characteristics analysis were carried out for six genes by bioinformation. Meanwhile, transient expression analysis of the Populus protoplasts showed that all the PeSCL9 proteins were localized in the nucleus. In addition, the degradome and sRNA of ‘Nanlin895’ poplar in combination were used to predict miRNAs that regulate PeSCL9. It was found that miR396a and miR396c may affect PeSCL9 expression via cleavage, which was further verified by a transient expression experiment in Populus protoplasts. Overall, the development of poplar adventitious root and other tissues was closely related to these six SCL9 genes, and they serve as a starting point for further research into the mechanisms regulating poplar growth and development.

GWAS Mediated Elucidation of Heterosis for Metric Traits in Cotton (Gossypium hirsutum L.) Across Multiple Environments

For about a century, plant breeding has widely exploited the heterosis phenomenon-often considered as hybrid vigor-to increase agricultural productivity. The ensuing F₁ hybrids can substantially outperform their progenitors due to heterozygous combinations that mitigate deleterious mutations occurring in each genome. However, only fragmented knowledge is available concerning the underlying genes and processes that foster heterosis. Although cotton is among the highly valued crops, its improvement programs that involve the exploitation of heterosis are still limited in terms of significant accomplishments to make it broadly applicable in different agro-ecological zones. Here, F₁ hybrids were derived from mating a diverse Upland Cotton germplasm with commercially valuable cultivars in the Line × Tester fashion and evaluated across multiple environments for 10 measurable traits. These traits were dissected into five different heterosis types and specific combining ability (SCA). Subsequent genome-wide predictions along-with association analyses uncovered a set of 298 highly significant key single nucleotide polymorphisms (SNPs)/Quantitative Trait Nucleotides (QTNs) and 271 heterotic Quantitative Trait Nucleotides (hQTNs) related to agronomic and fiber quality traits. The integration of a genome wide association study with RNA-sequence analysis yielded 275 candidate genes in the vicinity of key SNPs/QTNs. Fiber micronaire (MIC) and lint percentage (LP) had the maximum number of associated genes, i.e., each with 45 related to QTNs/hQTNs. A total of 54 putative candidate genes were identified in association with HETEROSIS of quoted traits. The novel players in the heterosis mechanism highlighted in this study may prove to be scientifically and biologically important for cotton biologists, and for those breeders engaged in cotton fiber and yield improvement programs.

A De Novo Transcriptome Assembly of Ceratopteris richardii Provides Insights into the Evolutionary Dynamics of Complex Gene Families in Land Plants

As the closest extant sister group to seed plants, ferns are an important reference point to study the origin and evolution of plant genes and traits. One bottleneck to the use of ferns in phylogenetic and genetic studies is the fact that genome-level sequence information of this group is limited, due to the extreme genome sizes of most ferns. Ceratopteris richardii (hereafter Ceratopteris) has been widely used as a model system for ferns. In this study, we generated a transcriptome of Ceratopteris, through the de novo assembly of the RNA-seq data from 17 sequencing libraries that are derived from two sexual types of gametophytes and five different sporophyte tissues. The Ceratopteris transcriptome, together with 38 genomes and transcriptomes from other species across the Viridiplantae, were used to uncover the evolutionary dynamics of orthogroups (predicted gene families using OrthoFinder) within the euphyllophytes and identify proteins associated with the major shifts in plant morphology and physiology that occurred in the last common ancestors of euphyllophytes, ferns, and seed plants. Furthermore, this resource was used to identify and classify the GRAS domain transcriptional regulators of many developmental processes in plants. Through the phylogenetic analysis within each of the 15 GRAS orthogroups, we uncovered which GRAS family members are conserved or have diversified in ferns and seed plants. Taken together, the transcriptome database and analyses reported here provide an important platform for exploring the evolution of gene families in land plants and for studying gene function in seed-free vascular plants.

A multilayered cross-species analysis of GRAS transcription factors uncovered their functional networks in plant adaptation to the environment

Introduction

Environmental stress is both a major force of natural selection and a prime factor affecting crop qualities and yields. The impact of the GRAS [gibberellic acid-insensitive (GAI), repressor of GA1-3 mutant (RGA), and scarecrow (SCR)] family on plant development and the potential to resist environmental stress needs much emphasis.

Objectives

This study aims to investigate the evolution, expansion, and adaptive mechanisms of GRASs of important representative plants during polyploidization.

Methods

We explored the evolutionary characteristics of GRASs in 15 representative plant species by systematic biological analysis of the genome, transcriptome, metabolite, protein complex map and phenotype.

Results

The GRAS family was systematically identified from 15 representative plant species of scientific and agricultural importance. The detection of gene duplication types of GRASs in all species showed that the widespread expansion of GRASs in these species was mainly contributed by polyploidization events. Evolutionary analysis reveals that most species experience independent genome-wide duplication (WGD) events and that interspecies GRAS functions may be broadly conserved. Polyploidy-related Chenopodium quinoa GRASs (CqGRASs) and Arabidopsis thaliana GRASs (AtGRASs) formed robust networks with flavonoid pathways by crosstalk with auxin and photosynthetic pathways. Furthermore, Arabidopsis thaliana population transcriptomes and the 1000 Plants (OneKP) project confirmed that GRASs are components of flavonoid biosynthesis, which enables plants to adapt to the environment by promoting flavonoid accumulation. More importantly, the GRASs of important species that may potentially improve important agronomic traits were mapped through TAIR and RARGE-II publicly available phenotypic data. Determining protein interactions and target genes contributes to determining GRAS functions.

Conclusion

The results of this study suggest that polyploidy-related GRASs in multiple species may be a target for improving plant growth, development, and environmental adaptation.

Phytohormones in fruit development and maturation

Phytohormones are integral to the regulation of fruit development and maturation. This review expands upon current understanding of the relationship between hormone signaling and fruit development, emphasizing fleshy fruit and highlighting recent work in the model crop tomato (Solanum lycopersicum) and additional species. Fruit development comprises fruit set initiation, growth, and maturation and ripening. Fruit set transpires after fertilization and is associated with auxin and gibberellic acid (GA) signaling. Interaction between auxin and GAs, as well as other phytohormones, is mediated by auxin-responsive Aux/IAA and ARF proteins. Fruit growth consists of cell division and expansion, the former shown to be influenced by auxin signaling. While regulation of cell expansion is less thoroughly understood, evidence indicates synergistic regulation via both auxin and GAs, with input from additional hormones. Fruit maturation, a transitional phase that precipitates ripening, occurs when auxin and GA levels subside with a concurrent rise in abscisic acid (ABA) and ethylene. During fruit ripening, ethylene plays a clear role in climacteric fruits, whereas non-climacteric ripening is generally associated with ABA. Recent evidence indicates varying requirements for both hormones within both ripening physiologies, suggesting rebalancing and specification of roles for common regulators rather than reliance upon one. Numerous recent discoveries pertaining to the molecular basis of hormonal activity and crosstalk are discussed, while we also note that many questions remain such as the molecular basis of additional hormonal activities, the role of epigenome changes, and how prior discoveries translate to the plethora of angiosperm species.

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.

Constitutive Overexpression of RAM1 Leads to an Increase in Arbuscule Density in Brachypodium distachyon

Arbuscular mycorrhizal (AM) symbiosis is a mutually beneficial association of plants and fungi of the subphylum Glomeromycotina. Endosymbiotic AM fungi colonize the inner cortical cells of the roots, where they form branched hyphae called arbuscules that function in nutrient exchange with the plant. To support arbuscule development and subsequent bidirectional nutrient exchange, the root cortical cells undergo substantial transcriptional reprogramming. REDUCED ARBUSCULAR MYCORRHIZA1 (RAM1), previously studied in several dicot plant species, is a major regulator of this cortical cell transcriptional program. Here, we generated ram1 mutants and RAM1 overexpressors in a monocot, Brachypodium distachyon. The AM phenotypes of two ram1 lines revealed that RAM1 is only partly required to enable arbuscule development in B. distachyon Transgenic lines constitutively overexpressing BdRAM1 showed constitutive expression of AM-inducible genes even in the shoots. Following inoculation with AM fungi, BdRAM1-overexpressing plants showed higher arbuscule densities relative to controls, indicating the potential to manipulate the relative proportion of symbiotic interfaces via modulation of RAM1 However, the overexpressors also show altered expression of hormone biosynthesis genes and aberrant growth patterns, including stunted bushy shoots and poor seed set. While these phenotypes possibly provide additional clues about the scope of influence of BdRAM1, they also indicate that directed approaches to increase the density of symbiotic interfaces will require a more focused, potentially cell type specific manipulation of transcription factor gene expression.

DELLA family duplication events lead to different selective constraints in angiosperms

Gibberellic acid (GA) is a major plant hormone involved in several biological processes from the flowering to the symbiosis with microorganisms. Thus, the GA regulation is crucial for plant biology. This regulation occurs via the DELLA proteins that belong to the GRAS transcription factor family. DELLA proteins are characterised by a DELLA N-terminal and a GRAS C-terminal domains. It is well known that DELLA activity appears after the bryophytes divergence and then evolved in the vascular plant lineages. Here we present the phylogeny of DELLA across 75 species belonging to various lineages from algae, liverworts and angiosperms. Our study confirmed two main duplication events, the first occurring before the angiosperms divergence and the other specific to the eudicots lineage. Comparative analysis of DELLA subclades in angiosperms revealed the loss in Poaceae and strong alteration in other species of the DELLA functional domain in the DELLA2 clade. In addition, molecular evolution analysis suggests that each of the clades (named DELLA1.1, DELLA1.2 and DELLA2) evolved differently but copies of each subclade are under strong purifying selection. This also suggests that, although the DELLA functional domain is altered in DELLA2, DELLA2 orthologs are still functional and operate in a different way compared to DELLA1 copies. In angiosperms, additional duplication events occurred and led to duplicate copies in species, genus or family such as in the Fabaceae subfamily Papilionoideae. This duplication led to the formation of additional paralogs in the DELLA1.2 subclade (DELLA1.2.1 and DELLA1.2.2). Interestingly, both copies appeared to be under relaxing selection revealing different evolutionary fate of the DELLA duplicated copies.

Genome-wide identification and analysis of GRAS transcription factors in the bottle gourd genome

GRAS genes belong to the plant-specific transcription factors (TF’s) family that are known to be involved in plant growth and development. In this study, we have identified 37 genes from the bottle gourd genome that encodes for GRAS TF’s. Except for the SCLA, we were able to identify at least one gene from each of the 17 subfamilies. Gene structure and chromosomal analysis showed that maximum seven genes are present on Chr7 followed by six genes on Chr1. The subcellular location analysis revealed that most of the genes were localized in the nucleus, except for a few in chloroplast and mitochondria. Additionally, we have identified one tandem gene duplication event on Chr7 and three major motifs that were present in all the GRAS genes. Furthermore, the protein–protein interaction prediction and gene expression analysis showed five candidate hub-genes interact with various other genes and thus probably control the expression of interacting partners in different plant tissues. Overall, this study provides a comprehensive analysis of GRAS transcription factors in bottle gourd genome which could be further extended to other vegetable crops.

Step by step evolution of Indeterminate Domain (IDD) transcriptional regulators: from algae to angiosperms

INTRODUCTION

The Indeterminate Domain (IDD) proteins are a plant-specific subclass of C2H2 Zinc Finger transcription factors. Some of these transcription factors play roles in diverse aspects of plant metabolism and development, but the function of most of IDD genes is unknown and the molecular evolution of the subfamily has not been explored in detail.

METHODS

In this study, we mined available genome sequences of green plants (Viridiplantae) to reconstruct the phylogeny and then described the motifs/expression patterns of IDD genes.

KEY RESULTS

We identified the complete set of IDD genes of 16 Streptophyta genomes. We found that IDD and its sister clade STOP arose by a duplication at the base of Streptophyta. Once on land, the IDD genes duplicated extensively, giving rise to at least ten lineages. Some of these lineages were lost in extant non-vascular plants and gymnosperms, but all of them were retained in angiosperms, duplicating profoundly in dicots and monocots and acquiring, at the same time, surprising heterogeneity in their C-terminal regions and expression patterns.

CONCLUSIONS

IDDs were present in the last common ancestor of Streptophyta. On land, IDDs duplicated extensively, leading to ten lineages. Later, IDDs were recruited by angiosperms where they diversified greatly in number, C-terminal and expression patterns. Interestingly, such diversification occurred during the evolution of novel traits of the plant body. This study provides a solid framework of the orthology relationships of green land plant IDD transcription factors, thus increasing the accuracy of orthologue identification in model and non-model species and facilitating the identification of agronomically important genes related to plant metabolism and development.

Transcriptome-Wide Identification, Evolutionary Analysis, and GA Stress Response of the GRAS Gene Family in Panax ginseng C. A. Meyer

GRAS transcription factors are a kind of plant-specific transcription factor that have been found in a variety of plants. According to previous studies, GRAS proteins are widely involved in the physiological processes of plant signal transduction, stress, growth and development. The Jilin ginseng (Panax ginseng C.A. Meyer) is a heterogeneous tetraploid perennial herb of the Araliaceae family, ginseng genus. Important information regarding the GRAS transcription factors has not been reported in ginseng. In this study, 59 Panax ginseng GRAS (PgGRAS) genes were obtained from the Jilin ginseng transcriptome data and divided into 13 sub-families according to the classification of Arabidopsis thaliana. Through systematic evolution, structural variation, function and gene expression analysis, we further reveal GRAS’s potential function in plant growth processes and its stress response. The expression of PgGRAS genes responding to gibberellin acids (GAs) suggests that these genes could be activated after application concentration of GA. The qPCR analysis result shows that four PgGRAS genes belonging to the DELLA sub-family potentially have important roles in the GA stress response of ginseng hairy roots. This study provides not only a preliminary exploration of the potential functions of the GRAS genes in ginseng, but also valuable data for further exploration of the candidate PgGRAS genes of GA signaling in Jilin ginseng, especially their roles in ginseng hairy root development and GA stress response.

The hornwort genome and early land plant evolution

Hornworts, liverworts and mosses are three early diverging clades of land plants, and together comprise the bryophytes. Here, we report the draft genome sequence of the hornwort Anthoceros angustus. Phylogenomic inferences confirm the monophyly of bryophytes, with hornworts sister to liverworts and mosses. The simple morphology of hornworts correlates with low genetic redundancy in plant body plan, while the basic transcriptional regulation toolkit for plant development has already been established in this early land plant lineage. Although the Anthoceros genome is small and characterized by minimal redundancy, expansions are observed in gene families related to RNA editing, UV protection and desiccation tolerance. The genome of A. angustus bears the signatures of horizontally transferred genes from bacteria and fungi, in particular of genes operating in stress-response and metabolic pathways. Our study provides insight into the unique features of hornworts and their molecular adaptations to live on land.

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 association study identified novel candidate loci affecting wood formation in Norway spruce

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

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.