Genome-wide characterization and phylogenetic analysis of GSK gene family in three species of cotton: evidence for a role of some GSKs in fiber development and responses to stress

Genome-Wide and Expression Pattern Analysis of the DVL Gene Family Reveals GhM_A05G1032 Is Involved in Fuzz Development in G. hirsutum

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

HbBIN2 Functions in Plant Cold Stress Resistance through Modulation of HbICE1 Transcriptional Activity and ROS Homeostasis in Hevea brasiliensis

Cold stress poses significant limitations on the growth, latex yield, and ecological distribution of rubber trees (Hevea brasiliensis). The GSK3-like kinase plays a significant role in helping plants adapt to different biotic and abiotic stresses. However, the functions of GSK3-like kinase BR-INSENSITIVE 2 (BIN2) in Hevea brasiliensis remain elusive. Here, we identified HbBIN2s of Hevea brasiliensis and deciphered their roles in cold stress resistance. The transcript levels of HbBIN2s are upregulated by cold stress. In addition, HbBIN2s are present in both the nucleus and cytoplasm and have the ability to interact with the INDUCER OF CBF EXPRESSION1(HbICE1) transcription factor, a central component in cold signaling. HbBIN2 overexpression in Arabidopsis displays decreased tolerance to chilling stress with a lower survival rate and proline content but a higher level of electrolyte leakage (EL) and malondialdehyde (MDA) than wild type under cold stress. Meanwhile, HbBIN2 transgenic Arabidopsis treated with cold stress exhibits a significant increase in the accumulation of reactive oxygen species (ROS) and a decrease in the activity of antioxidant enzymes. Further investigation reveals that HbBIN2 inhibits the transcriptional activity of HbICE1, thereby attenuating the expression of C-REPEAT BINDING FACTOR (HbCBF1). Consistent with this, overexpression of HbBIN2 represses the expression of CBF pathway cold-regulated genes under cold stress. In conclusion, our findings indicate that HbBIN2 functions as a suppressor of cold stress resistance by modulating HbICE1 transcriptional activity and ROS homeostasis.

Auxin signaling module OsSK41-OsIAA10-OsARF regulates grain yield traits in rice

Auxin is an important phytohormone in plants, and auxin signaling pathways in rice play key roles in regulating its growth, development, and productivity. To investigate how rice grain yield traits are regulated by auxin signaling pathways and to facilitate their application in rice improvement, we validated the functional relationships among regulatory genes such as OsIAA10, OsSK41, and OsARF21 that are involved in one of the auxin (OsIAA10) signaling pathways. We assessed the phenotypic effects of these genes on several grain yield traits across two environments using knockout and/or overexpression transgenic lines. Based on the results, we constructed a model that showed how grain yield traits were regulated by OsIAA10 and OsTIR1, OsAFB2, and OsSK41 and OsmiR393 in the OsSK41-OsIAA10-OsARF module and by OsARF21 in the transcriptional regulation of downstream auxin response genes in the OsSK41-OsIAA10-OsARF module. The population genomic analyses revealed rich genetic diversity and the presence of major functional alleles at most of these loci in rice populations. The strong differentiation of many major alleles between Xian/indica and Geng/japonica subspecies and/or among modern varieties and landraces suggested that they contributed to improved productivity during evolution and breeding. We identified several important aspects associated with the genetic and molecular bases of rice grain and yield traits that were regulated by auxin signaling pathways. We also suggested rice auxin response factor (OsARF) activators as candidate target genes for improving specific target traits by overexpression and/or editing subspecies-specific alleles and by searching and pyramiding the 'best' gene allelic combinations at multiple regulatory genes in auxin signaling pathways in rice breeding programs.

Regulatory network of GSK3-like kinases and their role in plant stress response

Glycogen synthase kinase 3 (GSK3) family members are evolutionally conserved Ser/Thr protein kinases in mammals and plants. In plants, the GSK3s function as signaling hubs to integrate the perception and transduction of diverse signals required for plant development. Despite their role in the regulation of plant growth and development, emerging research has shed light on their multilayer function in plant stress responses. Here we review recent advances in the regulatory network of GSK3s and the involvement of GSK3s in plant adaptation to various abiotic and biotic stresses. We also discuss the molecular mechanisms underlying how plants cope with environmental stresses through GSK3s-hormones crosstalk, a pivotal biochemical pathway in plant stress responses. We believe that our overview of the versatile physiological functions of GSK3s and underlined molecular mechanism of GSK3s in plant stress response will not only opens further research on this important topic but also provide opportunities for developing stress-resilient crops through the use of genetic engineering technology.

Genomic identification of cotton SAC genes branded ovule and stress-related key genes in Gossypium hirsutum

SAC genes have been identified to play a variety of biological functions and responses to various stresses. Previously, SAC genes have been recognized in animals and Arabidopsis. For the very first time, we identified 157 SAC genes in eight cotton species including three diploids and five tetraploids with 23 SAC members in G. hirsutum. Evolutionary analysis classified all cotton SAC gene family members into five distinct groups. Cotton SAC genes showed conserved sequence logos and WGD or segmental duplication. Multiple synteny and collinearity analyses revealed gene family expansion and purifying selection pressure during evolution. G. hirsutum SAC genes showed uneven chromosomal distribution, multiple exons/introns, conserved protein motifs, and various growth and stress-related cis-elements. Expression pattern analysis revealed three GhSAC genes (GhSAC3, GhSAC14, and GhSAC20) preferentially expressed in flower, five genes (GhSAC1, GhSAC6, GhSAC9, GhSAC13, and GhSAC18) preferentially expressed in ovule and one gene (GhSAC5) preferentially expressed in fiber. Similarly, abiotic stress treatment verified that GhSAC5 was downregulated under all stresses, GhSAC6 and GhSAC9 were upregulated under NaCl treatment, and GhSAC9 and GhSAC18 were upregulated under PEG and heat treatment respectively. Overall, this study identified key genes related to flower, ovule, and fiber development and important genetic material for breeding cotton under abiotic stress conditions.

GhBES1 mediates brassinosteroid regulation of leaf size by activating expression of GhEXO2 in cotton (Gossypium hirsutum)

We proposed a working model of BR to promote leaf size through cell expansion. In the BR signaling pathway, GhBES1 affects cotton leaf size by binding to and activating the expression of the E-box element in the GhEXO2 promoter region. Brassinosteroid (BR) is an essential phytohormone that controls plant growth. However, the mechanisms of BR regulation of leaf size remain to be determined. Here, we found that the BR deficient cotton mutant pagoda1 (pag1) had a smaller leaf size than wild-type CRI24. The expression of EXORDIUM (GhEXO2) gene, was significantly downregulated in pag1. Silencing of BRI1-EMS-SUPPRESSOR 1 (GhBES1), inhibited leaf cell expansion and reduced leaf size. Overexpression of GhBES1.4 promoted leaf cell expansion and enlarged leaf size. Expression analysis showed GhEXO2 expression positively correlated with GhBES1 expression. In plants, altered expression of GhEXO2 promoted leaf cell expansion affecting leaf size. Furthermore, GhBES1.4 specifically binds to the E-box elements in the GhEXO2 promoter, inducing its expression. RNA-seq data revealed many down-regulated genes related to cell expansion in GhEXO2 silenced plants. In summary, we discovered a novel mechanism of BR regulation of leaf size through GhBES1 directly activating the expression of GhEXO2.

Genome-Wide Characterization and Phylogenetic Analysis of GSK Genes in Maize and Elucidation of Their General Role in Interaction with BZR1

Glycogen synthase kinase-3 (GSK-3) is a nonreceptor serine/threonine protein kinase that is involved in diverse processes, including cell development, photomorphogenesis, biotic and abiotic stress responses, and hormone signaling. In contrast with the deeply researched GSK family in Arabidopsis and rice, maize GSKs’ common bioinformatic features and protein functions are poorly understood. In this study, we identified 11 GSK genes in the maize (Zea mays L.) genome via homologous alignment, which we named Zeama;GSKs (ZmGSKs). The results of ZmGSK protein sequences, conserved motifs, and gene structures showed high similarities with each other. The phylogenetic analyses showed that a total of 11 genes from maize were divided into four clades. Furthermore, semi-quantitative RT-PCR analysis of the GSKs genes showed that ZmGSK1, ZmGSK2, ZmGSK4, ZmGSK5, ZmGSK8, ZmGSK9, ZmGSK10, and ZmGSK11 were expressed in all tissues; ZmGSK3, ZmGSK6, and ZmGSK7 were expressed in a specific organization. In addition, GSK expression profiles under hormone treatments demonstrated that the ZmGSK genes were induced under BR conditions, except for ZmGSK2 and ZmGSK5. ZmGSK genes were regulated under ABA conditions, except for ZmGSK1 and ZmGSK8. Finally, using the yeast two-hybrid and BiFC assay, we determined that clads II (ZmGSK1, ZmGSK4, ZmGSK7, ZmGSK8, and ZmGSK11) could interact with ZmBZR1. The results suggest that clade II of ZmGSKs is important for BR signaling and that ZmGSK1 may play a dominant role in BR signaling as the counterpart to BIN2. This study provides a foundation for the further study of GSK3 functions and could be helpful in devising strategies for improving maize.

Genome-wide identification and expression analysis of the GSK gene family in wheat (Triticum aestivum L.)

BACKGROUND

Plant glycogen synthase kinase 3/shaggy kinase (GSK3) proteins contain the conserved kinase domain and play a pivotal role in the regulation of plant growth and abiotic stress responses. Nonetheless, genome-wide analysis of the GSK gene family in wheat (Triticum aestivum L.) has not been reported.

METHODS AND RESULTS

Using high-quality wheat genome sequences, a comprehensive genome-wide characterization of the GSK gene family in wheat was conducted. Their phylogenetics, chromosome location, gene structure, conserved domains, promoter cis-elements, gene duplications, and network interactions were systematically analyzed. In this study, we identified 22 GSK genes in wheat genome that were unevenly distributed on nine wheat chromosomes. Based on phylogenetic analysis, the GSK genes from Arabidopsis, rice, barley, and wheat were clustered into four subfamilies. Gene structure and conserved protein motif analysis revealed that GSK proteins in the same subfamily share similar motif structures and exon/intron organization. Results from gene duplication analysis indicate that four segmental duplications events contribute to the expansion of the wheat GSK gene family. Promoter analysis indicated the participation of TaSK genes in response to the hormone, light and abiotic stress, and plant growth and development. Furthermore, gene network analysis found that five TaSKs were involved in the regulatory network and 130 gene pairs of network interactions were identified. The heat map generated from the available transcriptomic data revealed that the TaSKs exhibited preferential expression in specific tissues and different expression patterns under abiotic stress conditions. Moreover, results from qRT-PCR analysis revealed that the randomly selected TaSK genes were abundantly expressed in spikes and grains at one specific developmental stage, as well as in responding to drought and salt stress.

CONCLUSIONS

These findings clearly depicted the evolutionary processes and the characteristics, and expression profiles of the GSK gene family in wheat, revealed their role in wheat development and response to abiotic stress responses.

Glycogen synthase kinases in model and crop plants - From negative regulators of brassinosteroid signaling to multifaceted hubs of various signaling pathways and modulators of plant reproduction and yield

Glycogen synthase kinases, also known as SHAGGY-like Kinases (GSKs/SKs), are highly conserved serine/threonine protein kinases present both in animals and plants. Plant genomes contain multiple homologs of the GSK3 genes which participate in various biological processes. Plant GSKs/SKs, and their best known representative in Arabidopsis thaliana - Brassinosteroid Insentisive2 (BIN2/SK21) in particular, were first identified as components of the brassinosteroid (BR) signaling pathway. As phytohormones, BRs regulate a wide range of physiological processes in plants - from germination, cell division, elongation and differentiation to leaf senescence, and response to environmental stresses. The GSKs/SKs proteins belong to a group of several highly conserved components of the BR signaling which evolved early during evolution of this molecular relay. However, recent reports indicated that the GSKs/SKs proteins are also implicated in signaling pathways of other phytohormones and stress-response processes. As a consequence, the GSKs/SKs proteins became hubs of various signaling pathways and modulators of plant development and reproduction. Thus, it is very important to understand molecular mechanisms regulating activity of the GSKs/SKs proteins, but also to get insights into role of the GSKs/SKs proteins in modulation of stability and activity of various substrate proteins which participate in the numerous signaling pathways. Although elucidation of these aspects is still in progress, this review presents a comprehensive and detailed description of these processes and their implications for regulation of development, stress response, and reproduction of model and crop species. The GSKs/SKs proteins and their activity are modulated through phosphorylation and de-phosphorylation reactions which are regulated by various proteins. Importantly, both phosphorylations and de-phosphorylations may have positive and negative effects on the activity of the GSKs/SKs proteins. Additionally, the activity of the GSKs/SKs proteins is positively regulated by reactive oxygen species, whereas it is negatively regulated through ubiquitylation, deacetylation, and nitric oxide-mediated nitrosylation. On the other hand, the GSKs/SKs proteins interact with proteins representing various signaling pathways, and on the basis of the complicated network of interactions the GSKs/SKs proteins differentially regulate various physiological, developmental, stress response, and yield-related processes.

Genome-wide characterization and expression analysis of cystathionine β-synthase genes in plant development and abiotic stresses of cotton (Gossypium spp.)

Cystathionine β-synthase (CBS) domains containing proteins (CDCPs) form a large family and play roles in development via regulation of the thioredoxin system as well as abiotic and biotic stress responses of plant. However, the comprehensive study of CBS genes remained elusive in cotton. Here, we identified 237 CBS genes in 11 plant species and the phylogenetic analysis categorized CBS genes into four groups. Whole-genome or segmental with dispersed duplication events contributed to GhCBS gene family expansion. Moreover, orthologous/paralogous genes among three cotton species (G. hirsutum, G. arboreum, and G. raimondii) were detected from the syntenic map among eight plant species. Strong purifying selection for dicotyledonous and monocotyledonous CBS genes, and cis-elements related to plant growth and development, abiotic and hormonal response were observed. Transcriptomic data and qRT-PCR validation of 12 GhCBS genes indicated their critical role in ovule development as most of the genes showed high enrichment. Further, some of GhCBS (GhCBS5, GhCBS16, GhCBS17, GhCBS24, GhCBS25, GhCBS26, and GhCBS52) genes were regulated under various abiotic and hormonal treatments for different time points and involve in ovule and fiber development which provided key genes for future cotton breeding programs. In addition, transgenic tobacco plants overexpressing GhCBS4 transiently exhibited higher water and chlorophyll content indicating improved tolerance toward drought stress. Overall, this study provides the characterization of GhCBS genes for plant growth, abiotic and hormonal stresses, thereby, intimating their significance in cotton molecular breeding for resistant cultivars.

GSK3s: nodes of multilayer regulation of plant development and stress responses

Glycogen synthase kinase 3 (GSK3) family members are highly conserved serine/threonine protein kinases in eukaryotes. Unlike animals, plants have evolved with multiple homologs of GSK3s involved in a diverse array of biological processes. Emerging evidence suggests that GSK3s act as signaling hubs for integrating perception and transduction of diverse signals required for plant development and responses to abiotic and biotic cues. Here we review recent advances in understanding the molecular interactions between GSK3s and an expanding spectrum of their upstream regulators and downstream substrates in plants. We further discuss how GSK3s act as key signaling nodes of multilayer regulation of plant development and stress response through either being regulated at the post-translational level or regulating their substrates via phosphorylation.

Genome-wide analysis of ZAT gene family revealed GhZAT6 regulates salt stress tolerance in G. hirsutum

High salt environments can induce stress in different plants. The genes containing the ZAT domain constitute a family that belongs to a branch of the C₂H₂ family, which plays a vital role in responding to abiotic stresses. In this study, we identified 169 ZAT genes from seven plant species, including 44 ZAT genes from G. hirsutum. Phylogenetic tree analysis divided ZAT genes in six groups with conserved gene structure, protein motifs. Two C₂H₂ domains and an EAR domain and even chromosomal distribution on At and Dt sub-genome chromosomes of G. hirsutum was observed. GhZAT6 was primarily expressed in the root tissue and responded to NaCl and ABA treatments. Subcellular localization found that GhZAT6 was located in the nucleus and demonstrated transactivation activity during a transactivation activity assay. Arabidopsis transgenic lines overexpressing the GhZAT6 gene showed salt tolerance and grew more vigorously than WT on MS medium supplemented with 100 mmol NaCl. Additionally, the silencing of the GhZAT6 gene in cotton plants showed more obvious leaf wilting than the control plants, which were subjected to 400 mmol NaCl treatment. Next, the expressions of GhAPX1, GhFSD1, GhFSD2, and GhSOS3 were significantly lower in the GhZAT6-silenced plants treated with NaCl than the control. Based on these findings, GhZAT6 may be involved in the ABA pathway and mediate salt stress tolerance by regulating ROS-related gene expression.

BIN2 negatively regulates plant defence against Verticillium dahliae in Arabidopsis and cotton

Verticillium wilt is caused by the soil-borne vascular pathogen Verticillium dahliae, and affects a wide range of economically important crops, including upland cotton (Gossypium hirsutum). Previous studies showed that expression levels of BIN2 were significantly down-regulated during infestation with V. dahliae. However, the underlying molecular mechanism of BIN2 in plant regulation against V. dahliae remains enigmatic. Here, we characterized a protein kinase GhBIN2 from Gossypium hirsutum, and identified GhBIN2 as a negative regulator of resistance to V. dahliae. The Verticillium wilt resistance of Arabidopsis and cotton were significantly enhanced when BIN2 was knocked down. Constitutive expression of BIN2 attenuated plant resistance to V. dahliae. We found that BIN2 regulated plant endogenous JA content and influenced the expression of JA-responsive marker genes. Further analysis revealed that BIN2 interacted with and phosphorylated JAZ family proteins, key repressors of the JA signalling pathway in both Arabidopsis and cotton. Spectrometric analysis and site-directed mutagenesis showed that BIN2 phosphorylated AtJAZ1 at T196, resulting in the degradation of JAZ proteins. Collectively, these results show that BIN2 interacts with JAZ proteins and plays a negative role in plant resistance to V. dahliae. Thus, BIN2 may be a potential target gene for genetic engineering against Verticillium wilt in crops.

Identification and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

Identification of GhLOG gene family revealed that GhLOG3 is involved in regulating salinity tolerance in cotton (Gossypium hirsutum L.)

Cytokinin (CK) is an important plant hormone that promotes plant cell division and differentiation, and participates in salt response under osmotic stress. LOGs (LONELY GUY) are CK-activating enzymes involved in CK synthesis. The LOG gene family has not been comprehensively characterized in cotton. In this study we identified 151 LOG genes from nine plant species, including 28 LOG genes in Gossypium hirsutum. Phylogenetic analysis divided LOG genes into three groups. Exon/intron structures and protein motifs of GhLOG genes were highly conserved. Synteny analysis revealed that several gene loci were highly conserved between the A and D sub-genomes of G. hirsutum with purifying selection pressure during evolution. Expression profiles showed that most LOG genes were constitutively expressed in eight different tissues. Furthermore, LOG genes can be regulated by abiotic stresses and phytohormone treatments. Moreover, subcellular localization revealed that GhLOG3_At resides inside the cell membrane. Overexpression of GhLOG3 enhanced salt tolerance in Arabidopsis. Virus-induced gene silencing (VIGS) of GhLOG3_At in cotton enhanced sensitivity of plants to salt stress with increased H₂O₂ contents and decreased chlorophyll and proline (PRO) activity. Our results suggested that GhLOG3_At induces salt stress tolerance in cotton, and provides a basis for the use of CK synthesis genes to regulate cotton growth and stress resistance.

Identification and Characterization of the ERF Subfamily B3 Group Revealed GhERF13.12 Improves Salt Tolerance in Upland Cotton

The APETALA2 (AP2)/ethylene response factor plays vital functions in response to environmental stimulus. The ethylene response factor (ERF) subfamily B3 group belongs to the AP2/ERF superfamily and contains a single AP2/ERF domain. Phylogenetic analysis of the ERF subfamily B3 group genes from Arabdiposis thaliana, Gossypium arboreum, Gossypium hirsutum, and Gossypium raimondii made it possible to divide them into three groups and showed that the ERF subfamily B3 group genes are conserved in cotton. Collinearity analysis identified172 orthologous/paralogous gene pairs between G. arboreum and G. hirsutum; 178 between G. hirsutum and G. raimondii; and 1,392 in G. hirsutum. The GhERF subfamily B3 group gene family experienced massive gene family expansion through either segmental or whole genome duplication events, with most genes showing signature compatible with the action of purifying selection during evolution. Most G. hirsutum ERF subfamily B3 group genes are responsive to salt stress. GhERF13.12 transgenic Arabidopsis showed enhanced salt stress tolerance and exhibited regulation of related biochemical parameters and enhanced expression of genes participating in ABA signaling, proline biosynthesis, and ROS scavenging. In addition, the silencing of the GhERF13.12 gene leads to increased sensitivity to salt stress in cotton. These results indicate that the ERF subfamily B3 group had remained conserved during evolution and that GhERF13.12 induces salt stress tolerance in Arabidopsis and cotton.

Comparative analysis of the SPL gene family in five Rosaceae species: Fragaria vesca, Malus domestica, Prunus persica, Rubus occidentalis, and Pyrus pyrifolia

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are very important for the plant growth and development. Here 15 RoSPLs were identified in Rubus occidentalis. The conserved domains and motifs, phylogenetic relationships, posttranscriptional regulation, and physiological function of the 92 SPL family genes in Fragaria vesca, Malus domestica, Prunus persica, R. occidentalis, and Pyrus pyrifolia were analyzed. Sequence alignment and phylogenetic analysis showed the SPL proteins had sequence conservation, some FvSPLs could be lost or developed, and there was a closer relationship between M. domestica and P. pyrifolia, F. vesca and R. occidentalis, respectively. Genes with similar motifs clustering together in the same group had their functional redundancy. Based on the function of SPLs in Arabidopsis thaliana, these SPLs could be involved in vegetative transition from juvenile to adult, morphological change in the reproductive phase, anthocyanin biosynthesis, and defense stress. Forty-eight SPLs had complementary sequences of miR156, of which nine PrpSPLs in P. persica and eight RoSPLs in R. occidentalis as the potential targets of miR156 were reported for the first time, suggesting the conservative regulatory effects of miR156 and indicating the roles of miR156-SPL modules in plant growth, development, and defense response. It provides a basic understanding of SPLs in Rosaceae plants.

Genome-wide identification and expression analysis of the GSK gene family in Solanum tuberosum L. under abiotic stress and phytohormone treatments and functional characterization of StSK21 involvement in salt stress

Plant Glycogen Synthase Kinase 3 (GSK3)/SHAGGY-like kinase (GSK) proteins play important roles in modulating growth, development, and stress responses in several plant species. However, little is known about the members of the potato GSK (StGSK) family. Here, nine StGSK genes were identified and phylogenetically grouped into four clades. Gene duplication analysis revealed that segmental duplication contributed to the expansion of the StGSK family. Gene structure and motif pattern analyses indicated that similar exon/intron and motif organizations were found in StGSKs from the same clade. Conserved motif and kinase activity analyses indicated that the StGSKs encode active protein kinases, and they were shown to be distributed throughout whole cells. Cis-acting regulatory element analysis revealed the presence of many growth-, hormone-, and stress-responsive elements within the promoter regions of the StGSKs, which is consistent with their expression in different organs, and their altered expression in response to hormone and stress treatments. Association network analysis indicated that various proteins, including two confirmed BES1 family transcription factors, potentially interact with StGSKs. Overexpression of StSK21 provides enhanced sensitivity to salt stress in Arabidopsis thaliana plants. Overall, these results reveal that StGSK proteins are active protein kinases with purported functions in regulating growth, development, and stress responses.

Silencing of HvGSK1.1-A GSK3/SHAGGY-Like Kinase-Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions

Glycogen synthase kinase 3 (GSK3) is a highly conserved kinase present in all eukaryotes and functions as a key regulator of a wide range of physiological and developmental processes. The kinase, known in land plants as GSK3/SHAGGY-like kinase (GSK), is a key player in the brassinosteroid (BR) signaling pathway. The GSK genes, through the BRs, affect diverse developmental processes and modulate responses to environmental factors. In this work, we describe functional analysis of HvGSK1.1, which is one of the GSK3/SHAGGY-like orthologs in barley. The RNAi-mediated silencing of the target HvGSK1.1 gene was associated with modified expression of its paralogs HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1 in plants grown in normal and in salt stress conditions. Low nucleotide similarity between the silencing fragment and barley GSK genes and the presence of BR-dependent transcription factors’ binding sites in promoter regions of barley and rice GSK genes imply an innate mechanism responsible for co-regulation of the genes. The results of the leaf inclination assay indicated that silencing of HvGSK1.1 and the changes of GSK paralogs enhanced the BR-dependent signaling in the plants. The strongest phenotype of transgenic lines with downregulated HvGSK1.1 and GSK paralogs had greater biomass of the seedlings grown in normal conditions and salt stress as well as elevated kernel weight of plants grown in normal conditions. Both traits showed a strong negative correlation with the transcript level of the target gene and the paralogs. The characteristics of barley lines with silenced expression of HvGSK1.1 are compatible with the expected phenotypes of plants with enhanced BR signaling. The results show that manipulation of the GSK-encoding genes provides data to explore their biological functions and confirm it as a feasible strategy to generate plants with improved agricultural traits.

Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Brassinosteroid (BR), a steroid phytohormone, whose signaling transduction pathways include a series of phosphorylation and dephosphorylation events, and GSK3s are the main negative regulator kinases. BRs have been shown to play vital roles in cotton fiber elongation. However, the underlying mechanism is still elusive. In this study, fibers of a BR-defective mutant Pagoda 1 (pag1), and its corresponding wild-type (ZM24) were selected for a comparative global phosphoproteome analysis at critical developmental time points: fast-growing stage (10 days after pollination (DPA)) and secondary cell wall synthesis stage (20 DPA). Based on the substrate characteristics of GSK3, 900 potential substrates were identified. Their GO and KEGG annotation results suggest that BR functions in fiber development by regulating GhSKs (GSK3s of Gossypium hirsutum L.) involved microtubule cytoskeleton organization, and pathways of glucose, sucrose and lipid metabolism. Further experimental results revealed that among the GhSK members identified, GhSK13 not only plays a role in BR signaling pathway, but also functions in developing fiber by respectively interacting with an AP2-like ethylene-responsive factor GhAP2L, a nuclear transcription factor Gh_DNF_YB19, and a homeodomain zipper member GhHDZ5. Overall, our phosphoproteomic research advances the understanding of fiber development controlled by BR signal pathways especially through GhSKs, and also offers numbers of target proteins for improving cotton fiber quality.

Casparian strip membrane domain proteins in Gossypium arboreum: genome-wide identification and negative regulation of lateral root growth

Background

Root systems are critical for plant growth and development. The Casparian strip in root systems is involved in stress resistance and maintaining homeostasis. Casparian strip membrane domain proteins (CASPs) are responsible for the formation of Casparian strips.

Results

To investigate the function of CASPs in cotton, we identified and characterized 48, 54, 91 and 94 CASPs from Gossypium arboreum, Gossypium raimondii, Gossypium barbadense and Gossypium hirsutum, respectively, at the genome-wide level. However, only 29 common homologous CASP genes were detected in the four Gossypium species. A collinearity analysis revealed that whole genome duplication (WGD) was the primary reason for the expansion of the genes of the CASP family in the four cotton species. However, dispersed duplication could also contribute to the expansion of the GaCASPs gene family in the ancestors of G. arboreum. Phylogenetic analysis was used to cluster a total of 85 CASP genes from G. arboreum and Arabidopsis into six distinct groups, while the genetic structure and motifs of CASPs were conserved in the same group. Most GaCASPs were expressed in diverse tissues, with the exception of that five GaCASPs (Ga08G0113, Ga08G0114, Ga08G0116, Ga08G0117 and Ga08G0118) that were highly expressed in root tissues. Analyses of the tissue and subcellular localization suggested that GaCASP27 genes (Ga08G0117) are membrane protein genes located in the root. In the GaCASP27 silenced plants and the Arabidopsis mutants, the lateral root number significantly increased. Furthermore, GaMYB36, which is related to root development was found to regulate lateral root growth by targeting GaCASP27.

Conclusions

This study provides a fundamental understanding of the CASP gene family in cotton and demonstrates the regulatory role of GaCASP27 on lateral root growth and development.

Genome-wide identification and expression analysis of PUB genes in cotton

Background

The U-box gene encodes a ubiquitin ligase that contain U-box domain. The plant U-box gene (PUB) plays an important role in the response to stresses, but few reports about PUBs in cotton were available. Therefore research on PUBs is of great importance and a necessity when studying the mechanisms of stress- tolerance in cotton.

Results

In this study, we identified 93, 96, 185 and 208 PUBs from four sequenced cotton species G. raimondii (D₅), G. arboreum (A₂), G. hirsutum (AD₁) and G. barbadense (AD2), respectively. Prediction analysis of subcellular localization showed that the PUBs in cotton were widely localized in cells, but primarily in the nucleus. The PUBs in cotton were classified into six subfamilies (A-F) on the basis of phylogenetic analysis, which was testified by the analysis of conserved motifs and exon-intron structures. Chromosomal localization analysis showed that cotton PUBs were unevenly anchored on all chromosomes, varying from 1 to 14 per chromosome. Through multiple sequence alignment analysis, 3 tandem duplications and 28 segmental duplications in cotton genome D₅, 2 tandem duplications and 25 segmental duplications in A₂, and 143 homologous gene pairs in A₂ and D₅ were found; however no tandem duplications in A₂ or D₅ were found. Additionally, 105, 14 and 17 homologous gene pairs were found in the intra-subgenome of At and Dt, At sub-genome and Dt sub-genome of G. hirsutum, respectively. Functional analysis of GhPUB85A and GhPUB45D showed that these genes positively responded to abiotic stresses, but the expression patterns were different. In addition, although the expression levels of these two homologous genes were similar, their contributions were different when responding to stresses, specifically showing different responses to abiotic stresses and functional differences between the two subgenomes of G. hirsutum.

Conclusions

This study reported the genome-wide identification, structure, evolution and expression analysis of PUBs in cotton, and the results showed that the PUBs were highly conserved throughout the evolutionary history of cotton. All PUB genes were involved in the response to abiotic stresses (including salt, drought, hot and cold) to varying degrees.

MiR-21 upregulation increases IL-8 expression and tumorigenesis program in airway epithelial cells exposed to cigarette smoke

BACKGROUND

Cigarette smoke exposure, increasing Toll-like receptor 4 (TLR4) and reactive oxygen species (ROS), promotes inflammatory responses in airway epithelial cells. Chronic inflammation, microRNA (miRNA), and oxidative stress are associated with cancer development.

AIMS

The present study was aimed to explore whether cigarette smoke exposure, altering miR-21 expression, promoted inflammatory responses and tumorigenesis processes in airway epithelial cells.

METHODS

Airway normal and cancer epithelial cells (16HBE and A549) were exposed to cigarette smoke extracts (CSE) or with/without agomiR-21, and then it was assessed: a) miR-21 expression; b) signal transducer and activator of transcription 3 (STAT3) nuclear protein expression and ERK1/2 activation; c) IL-8 gene expression and protein release. An antagonist of TLR4 (CLI-095) and the antioxidant flavonoid, apigenin, were also included to evaluate miR-21 expression in CSE exposed cells.

RESULTS

It was demonstrated that: a) A549 cells constitutively expressed higher levels of miR-21 and IL-8; b) CSE increased STAT3 nuclear expression in 16HBE; c) in both cell lines, CSE and agomiR-21 increased: miR-21 expression; ERK1/2 activation and IL-8 gene expression and protein release; d) TLR4 inhibition counteracted the effects of CSE on miR-21 in A549; e) apigenin reduced miR-21 and IL-8 gene expression in both cell lines.

CONCLUSIONS

Data herein provided identified for the first time new mechanisms supporting the crucial role of cigarette smoke-induced miR-21 expression in the amplification of inflammatory responses and in tumorigenesis processes within the airways.