Transcriptome analysis of newly classified bZIP transcription factors of Brassica rapa in cold stress response

Comparative Metabolome and Transcriptome Analysis of Rapeseed (Brassica napus L.) Cotyledons in Response to Cold Stress

Cold stress affects the seed germination and early growth of winter rapeseed, leading to yield losses. We employed transmission electron microscopy, physiological analyses, metabolome profiling, and transcriptome sequencing to understand the effect of cold stress (0 °C, LW) on the cotyledons of cold-tolerant (GX74) and -sensitive (XY15) rapeseeds. The mesophyll cells in cold-treated XY15 were severely damaged compared to slightly damaged cells in GX74. The fructose, glucose, malondialdehyde, and proline contents increased after cold stress in both genotypes; however, GX74 had significantly higher content than XY15. The pyruvic acid content increased after cold stress in GX74, but decreased in XY15. Metabolome analysis detected 590 compounds, of which 32 and 74 were differentially accumulated in GX74 (CK vs. cold stress) and XY15 (CK vs. cold stressed). Arachidonic acid and magnoflorine were the most up-accumulated metabolites in GX74 subjected to cold stress compared to CK. There were 461 and 1481 differentially expressed genes (DEGs) specific to XY15 and GX74 rapeseeds, respectively. Generally, the commonly expressed genes had higher expressions in GX74 compared to XY15 in CK and cold stress conditions. The expression changes in DEGs related to photosynthesis-antenna proteins, chlorophyll biosynthesis, and sugar biosynthesis-related pathways were consistent with the fructose and glucose levels in cotyledons. Compared to XY15, GX74 showed upregulation of a higher number of genes/transcripts related to arachidonic acid, pyruvic acid, arginine and proline biosynthesis, cell wall changes, reactive oxygen species scavenging, cold-responsive pathways, and phytohormone-related pathways. Taken together, our results provide a detailed overview of the cold stress responses in rapeseed cotyledons.

In-Depth Characterization of bZIP Genes in the Context of Endoplasmic Reticulum (ER) Stress in Brassica campestris ssp. chinensis

Numerous studies have been conducted to investigate the genomic characterization of bZIP genes and their involvement in the cellular response to endoplasmic reticulum (ER) stress. These studies have provided valuable insights into the coordinated cellular response to ER stress, which is mediated by bZIP transcription factors (TFs). However, a comprehensive and systematic investigations regarding the role of bZIP genes and their involvement in ER stress response in pak choi is currently lacking in the existing literature. To address this knowledge gap, the current study was initiated to elucidate the genomic characteristics of bZIP genes, gain insight into their expression patterns during ER stress in pak choi, and investigate the protein-to-protein interaction of bZIP genes with the ER chaperone BiP. In total, 112 members of the BcbZIP genes were identified through a comprehensive genome-wide analysis. Based on an analysis of sequence similarity, gene structure, conserved domains, and responsive motifs, the identified BcbZIP genes were categorized into 10 distinct subfamilies through phylogenetic analysis. Chromosomal location and duplication events provided insight into their genomic context and evolutionary history. Divergence analysis estimated their evolutionary history with a predicted divergence time ranging from 0.73 to 80.71 million years ago (MYA). Promoter regions of the BcbZIP genes were discovered to exhibit a wide variety of cis-elements, including light, hormone, and stress-responsive elements. GO enrichment analysis further confirmed their roles in the ER unfolded protein response (UPR), while co-expression network analysis showed a strong relationship of BcbZIP genes with ER-stress-responsive genes. Moreover, gene expression profiles and protein-protein interaction with ER chaperone BiP further confirmed their roles and capacity to respond to ER stress in pak choi.

Genome-Wide Identification of bZIP Transcription Factors in Cymbidium ensifolium and Analysis of Their Expression under Low-Temperature Stress

The basic leucine zipper (bZIP) transcription factors constitute the most widely distributed and conserved eukaryotic family. They play crucial roles in plant growth, development, and responses to both biotic and abiotic stresses, exerting strong regulatory control over the expression of downstream genes. In this study, a genome-wide characterization of the CebZIP transcription factor family was conducted using bioinformatic analysis. Various aspects, including physicochemical properties, phylogenetics, conserved structural domains, gene structures, chromosomal distribution, gene covariance relationships, promoter cis-acting elements, and gene expression patterns, were thoroughly analyzed. A total of 70 CebZIP genes were identified from the C. ensifolium genome, and they were randomly distributed across 18 chromosomes. The phylogenetic tree clustered them into 11 subfamilies, each exhibiting complex gene structures and conserved motifs arranged in a specific order. Nineteen pairs of duplicated genes were identified among the 70 CebZIP genes, with sixteen pairs affected by purifying selection. Cis-acting elements analysis revealed a plethora of regulatory elements associated with stress response, plant hormones, and plant growth and development. Transcriptome and qRT-PCR results demonstrated that the expression of CebZIP genes was universally up-regulated under low temperature conditions. However, the expression patterns varied among different members. This study provides theoretical references for identifying key bZIP genes in C. ensifolium that confer resistance to low-temperature stress, and lays the groundwork for further research into their broader biological functions.

Characterization of the basic leucine zipper transcription factor family of Neoporphyra haitanensis and its role in acclimation to dehydration stress

BACKGROUND

Neoporphyra haitanensis, a major marine crop native to southern China, grows in the harsh intertidal habitats of rocky coasts. The thallus can tolerate fluctuating and extreme environmental stresses, for example, repeated desiccation/rehydration due to the turning tides. It is also a typical model system for investigating stress tolerance mechanisms in intertidal seaweed. The basic leucine zipper (bZIP) transcription factors play important roles in the regulation of plants' responses to environmental stress stimuli. However, little information is available regarding the bZIP family in the marine crop Nh. haitanensis.

RESULTS

We identified 19 bZIP genes in the Nh. haitanensis genome and described their conserved domains. Based on phylogenetic analysis, these 19 NhhbZIP genes, distributed unevenly on the 11 superscaffolds, were divided into four groups. In each group, there were analogous exon/intron numbers and motif compositions, along with diverse exon lengths. Cross-species collinearity analysis indicated that 17 and 9 NhhbZIP genes were orthologous to bZIP genes in Neopyropia yezoensis and Porphyra umbilicalis, respectively. Evidence from RNA sequencing (RNA-seq) data showed that the majority of NhhbZIP genes (73.68%) exhibited transcript abundance in all treatments. Furthermore, genes NN 2, 4 and 5 showed significantly altered expression in response to moderate dehydration, severe dehydration, and rehydration, respectively. Gene co-expression network analysis of the representative genes was carried out, followed by gene set enrichment analysis. Two NhhbZIP genes collectively responding to dehydration and rehydration and their co-expressing genes mainly participated in DNA repair, DNA metabolic process, and regulation of helicase activity. Two specific NhhbZIP genes responding to severe dehydration and their corresponding network genes were mainly involved in macromolecule modification, cellular catabolic process, and transmembrane transport. Three specific NhhbZIP genes responding to rehydration and their co-expression gene networks were mainly involved in the regulation of the cell cycle process and defense response.

CONCLUSIONS

This study provides new insights into the structural composition, evolution, and function of the NhhbZIP gene family. Our results will help us to further study the functions of bZIP genes in response to dehydration and rehydration in Nh. haitanensis and improve Nh. haitanensis in southern China.

Role of bZIP transcription factors in the regulation of plant secondary metabolism

MAIN CONCLUSION

This study provides an overview of the structure, classification, regulatory mechanisms, and biological functions of the basic (region) leucine zipper transcription factors and their molecular mechanisms in flavonoid, terpenoid, alkaloid, phenolic acid, and lignin biosynthesis. Basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors (TFs) in eukaryotic organisms. The bZIP TFs are widely distributed in plants and play important roles in plant growth and development, photomorphogenesis, signal transduction, resistance to pathogenic microbes, biotic and abiotic stress, and secondary metabolism. Moreover, the expression of bZIP TFs not only promotes or inhibits the accumulation of secondary metabolites in medicinal plants, but also affects the stress response of plants to the external adverse environment. This paper describes the structure, classification, biological function, and regulatory mechanisms of bZIP TFs. In addition, the molecular mechanism of bZIP TFs regulating the biosynthesis of flavonoids, terpenoids, alkaloids, phenolic acids, and lignin are also elaborated. This review provides a summary for in-depth study of the molecular mechanism of bZIP TFs regulating the synthesis pathway of secondary metabolites and plant molecular breeding, which is of significance for the generation of beneficial secondary metabolites and the improvement of plant varieties.

Kiwifruit bZIP transcription factor AcePosF21 elicits ascorbic acid biosynthesis during cold stress

Cold stress seriously affects plant development, resulting in heavy agricultural losses. L-ascorbic acid (AsA, vitamin C) is an antioxidant implicated in abiotic stress tolerance and metabolism of reactive oxygen species (ROS). Understanding whether and how cold stress elicits AsA biosynthesis to reduce oxidative damage is important for developing cold-resistant plants. Here, we show that the accumulation of AsA in response to cold stress is a common mechanism conserved across the plant kingdom, from single-cell algae to angiosperms. We identified a basic leucine zipper domain (bZIP) transcription factor (TF) of kiwifruit (Actinidia eriantha Benth.), AcePosF21, which was triggered by cold and is involved in the regulation of kiwifruit AsA biosynthesis and defense responses against cold stress. AcePosF21 interacted with the R2R3-MYB TF AceMYB102 and directly bound to the promoter of the gene encoding GDP-L-galactose phosphorylase 3 (AceGGP3), a key conduit for regulating AsA biosynthesis, to up-regulate AceGGP3 expression and produce more AsA, which neutralized the excess ROS induced by cold stress. On the contrary, VIGS or CRISPR-Cas9-mediated editing of AcePosF21 decreased AsA content and increased the generation of ROS in kiwifruit under cold stress. Taken together, we illustrated a model for the regulatory mechanism of AcePosF21-mediated regulation of AceGGP3 expression and AsA biosynthesis to reduce oxidative damage by cold stress, which provides valuable clues for manipulating the cold resistance of kiwifruit.

Genome-wide identification and expression profiling of the bZIP gene family in Betula platyphylla and the functional characterization of BpChr04G00610 under low-temperature stress

The basic leucine zipper (bZIP) gene, which plays a significant role in the regulation of tolerance to biotic/abiotic stresses, has been characterized in many plant species. Betula platyphylla is a significant afforestation species. To elucidate the stress resistance mechanism of birch, previous studies identified some stress resistance genes. However, the genome-wide identification and characterization of bZIP gene family in the birch have not been reported. Here, the 56 BpbZIP genes were identified and classified into 13 groups in birch. Cis-element analysis showed that the promoters of 56 family genes contained 108 elements, of which 16 were shared by 13 groups. There were 8 pairs of fragment repeats and 1 pair of tandem repeats, indicating that duplication may be the major reason for the amplification of the BpbZIP gene family. Tissue-specific of BpbZIP genes showed 18 genes with the highest expression in roots, 15 in flowers, 11 in xylem and 9 in leaves. In addition, five differentially expressed bZIP genes were identified from the RNA-seq data of birch under low-temperature stress, and the co-expressed differentially expressed genes were further screened. The analysis of gene ontology (GO) enrichment of each co-expression regulatory network showed that they were related to membrane lipids and cell walls. Furthermore, the transient overexpression of BpChr04G00610 decreased the ROS scavenging ability of birch under low-temperature stress, suggesting that it may be more sensitive to low-temperature. In conclusion, this study provides a basis for the study of the function of BpbZIP genes.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa)

Zinc-finger protein (ZFP) transcription factors are among the largest families of transcription factors in plants. They participate in various biological processes such as apoptosis, autophagy, and stemness maintenance and play important roles in regulating plant growth and development and the response to stress. To elucidate the functions of ZFP genes in the low-temperature response of winter (Brassica rapa L.) B. rapa, this study identified 141 members of the C2H2 ZFP gene family from B. rapa, which are heterogeneously distributed on 10 chromosomes and have multiple cis-acting elements related to hormone regulation and abiotic stress of adversity. Most of the genes in this family contain only one CDS, and genes distributed in the same evolutionary branch share mostly the same motifs and are highly conserved in the evolution of cruciferous species. The genes were significantly upregulated in the roots and growth cones of ‘Longyou-7’, indicating that they play a role in the stress-response process of winter B. rapa. The expression level of the Bra002528 gene was higher in the strongly cold-resistant varieties than in the weakly cold-resistant varieties after low-temperature stress. The survival rate and BrZAT12 gene expression of trans-BrZAT12 Arabidopsis thaliana (Arabidopsis) were significantly higher than those of the wild-type plants at low temperature, and the enzyme activities in vivo were higher than those of the wild-type plants, indicating that the BrZAT12 gene could improve the cold resistance of winter B. rapa. BrZAT12 expression and superoxide dismutase and ascorbate peroxidase enzyme activities were upregulated in winter B. rapa after exogenous ABA treatment. BrZAT12 expression and enzyme activities decreased after the PD98059 treatment, and BrZAT12 expression and enzyme activities were higher than in the PD98059 treatment but lower than in the control after both treatments together. It is speculated that BrZAT12 plays a role in the ABA signaling process in which MAPKK is involved. This study provides a theoretical basis for the resolution of cold-resistance mechanisms in strong winter B. rapa.

Genome-wide identification and expression characterization of the DoG gene family of moso bamboo (Phyllostachys edulis)

BACKGROUND

The DoG (Delay of Germination1) family plays a key regulatory role in seed dormancy and germination. However, to date, there is no complete genomic overview of the DoG gene family of any economically valuable crop, including moso bamboo (Phyllostachys edulis), and no studies have been conducted to characterize its expression profile. To identify the DoG gene members of moso bamboo (PeDoG) and to investigate their family structural features and tissue expression profile characteristics, a study was conducted. Based on the whole genome and differential transcriptome data, in this investigation, we have scrutinized the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the PeDoG1 family of moso bamboo.

RESULTS

The DoG family genes of moso bamboo were found distributed across 16 chromosomal scaffolds with 24 members. All members were found to carry DoG1 structural domains, while 23 members additionally possessed basic leucine zipper (bZIP) structural domains. We could divide the PeDoG genes into three subfamilies based on phylogenetic relationships. Covariance analysis revealed that tandem duplication was the main driver of amplification of the PeDoG genes. The upstream promoter of these genes containing several cis-acting elements indicates a plausible role in abiotic stress and hormone induction. Gene expression pattern according to transcriptome data revealed participation of the PeDoG genes in tissue and organ development. Analysis using Short Time-series Expression Miner (STEM) tool revealed that the PeDoG gene family is also associated with rapid early shoot growth. Gene ontology (GO) and KEGG analyses showed a dual role of the PeDoG genes. We found that PeDoGs has a possible role as bZIP transcription factors by regulating Polar like1 (PL1) gene expression, and thereby playing a disease response role in moso bamboo. Quantitative gene expression of the PeDoG genes revealed that they were abundantly expressed in roots and leaves, and could be induced in response to gibberellin (GA).

CONCLUSION

In this study, we found that the PeDoG genes are involved in a wide range of activities such as growth and development, stress response and transcription. This forms the first report of PeDoG genes and their potential roles in moso bamboo.

Genome-wide identification and expression profiling of basic leucine zipper transcription factors following abiotic stresses in potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important food crop that is grown and consumed worldwide. The growth and productivity of this crop are severely affected by various abiotic stresses. Basic leucine zipper (bZIP) transcription factors (TFs) in plants are well known for their function during growth and development. However, systematic and in-depth identification and functional characterization of the bZIP gene family of potato is lacking. In the current study, we identified a total of 90 bZIPs (StbZIP) distributed on 12 linkage groups of potato. Based on the previous functional annotation and classification of bZIPs in Arabidopsis, wheat, and rice, a phylogenetic tree of potato bZIPs was constructed and genes were categorized into various functional groups (A to I, S, and U) as previously annotated in Arabidopsis thaliana. Analyses of the transcript sequence (RNA-seq) data led to identifying a total of 18 candidate StbZIPs [four in roots, eight in the tuber, six in mesocarp and endocarp] that were expressed in a tissue-specific manner. Differential expression analysis under the various abiotic conditions (salt, mannitol, water, and heat stress) and treatment with phytohormones (ABA, GA, IAA, and BAP) led to the identification of forty-two [thirteen under salt stress, two under mannitol stress, ten under water stress, and eighteen under heat stress], and eleven [eight and three StbZIPs upon treatment with ABA, and IAA, respectively] candidate StbZIPs, respectively. Using sequence information of candidate StbZIPs, a total of 22 SSR markers were also identified in this study. In conclusion, the genome-wide identification analysis coupled with RNA-Seq expression data led to identifying candidate StbZIPs, which are dysregulated, and may play a pivotal role under various abiotic stress conditions. This study will pave the way for future functional studies using forward and reverse genetics to improve abiotic stress tolerance in potato.

Genome-wide identification and expression profiling of basic leucine zipper transcription factors following abiotic stresses in potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important food crop that is grown and consumed worldwide. The growth and productivity of this crop are severely affected by various abiotic stresses. Basic leucine zipper (bZIP) transcription factors (TFs) in plants are well known for their function during growth and development. However, systematic and in-depth identification and functional characterization of the bZIP gene family of potato is lacking. In the current study, we identified a total of 90 bZIPs (StbZIP) distributed on 12 linkage groups of potato. Based on the previous functional annotation and classification of bZIPs in Arabidopsis, wheat, and rice, a phylogenetic tree of potato bZIPs was constructed and genes were categorized into various functional groups (A to I, S, and U) as previously annotated in Arabidopsis thaliana. Analyses of the transcript sequence (RNA-seq) data led to identifying a total of 18 candidate StbZIPs [four in roots, eight in the tuber, six in mesocarp and endocarp] that were expressed in a tissue-specific manner. Differential expression analysis under the various abiotic conditions (salt, mannitol, water, and heat stress) and treatment with phytohormones (ABA, GA, IAA, and BAP) led to the identification of forty-two [thirteen under salt stress, two under mannitol stress, ten under water stress, and eighteen under heat stress], and eleven [eight and three StbZIPs upon treatment with ABA, and IAA, respectively] candidate StbZIPs, respectively. Using sequence information of candidate StbZIPs, a total of 22 SSR markers were also identified in this study. In conclusion, the genome-wide identification analysis coupled with RNA-Seq expression data led to identifying candidate StbZIPs, which are dysregulated, and may play a pivotal role under various abiotic stress conditions. This study will pave the way for future functional studies using forward and reverse genetics to improve abiotic stress tolerance in potato.

The bZIP gene family in watermelon: genome-wide identification and expression analysis under cold stress and root-knot nematode infection

The basic leucine zipper (bZIP) family transcription factors play crucial roles in regulating plant development and stress response. In this study, we identified 62 ClabZIP genes from watermelon genome, which were unevenly distributed across the 11 chromosomes. These ClabZIP proteins could be classified into 13 groups based on the phylogenetic relationships, and members in the same group showed similar compositions of conserved motifs and gene structures. Transcriptome analysis revealed that a number of ClabZIP genes have important roles in the melatonin (MT) induction of cold tolerance. In addition, some ClabZIP genes were induced or repressed under red light (RL) or root-knot nematode infection according to the transcriptome data, and the expression patterns of several ClabZIP genes were further verified by quantitative real-time PCR, revealing their possible roles in RL induction of watermelon defense against nematode infection. Our results provide new insights into the functions of different ClabZIP genes in watermelon and their roles in response to cold stress and nematode infection.

Quantitative Proteomics Analysis of Lettuce (Lactuca sativa L.) Reveals Molecular Basis-Associated Auxin and Photosynthesis with Bolting Induced by High Temperature

Bolting is a key process in the growth and development of lettuce (Lactuca sativa L.). A high temperature can induce early bolting, which decreases both the quality and production of lettuce. However, knowledge of underlying lettuce bolting is still lacking. To better understand the molecular basis of bolting, a comparative proteomics analysis was conducted on lettuce stems, during the bolting period induced by a high temperature (33 °C) and a control temperature (20 °C) using iTRAQ-based proteomics, phenotypic measures, and biological verifications using qRT-PCR and Western blot. The high temperature induced lettuce bolting, while the control temperature did not. Of the 5454 identified proteins, 619 proteins presented differential abundance induced by high-temperature relative to the control group, of which 345 had an increased abundance and 274 had a decreased abundance. Proteins with an abundance level change were mainly enriched in pathways associated with photosynthesis and tryptophan metabolism involved in auxin (IAA) biosynthesis. Moreover, among the proteins with differential abundance, proteins associated with photosynthesis and tryptophan metabolism were increased. These findings indicate that a high temperature enhances the function of photosynthesis and IAA biosynthesis to promote the process of bolting, which is in line with the physiology and transcription level of IAA metabolism. Our data provide a first comprehensive dataset for gaining novel understanding of the molecular basis underlying lettuce bolting induced by high temperature. It is potentially important for further functional analysis and genetic manipulation for molecular breeding to breed new cultivars of lettuce to restrain early bolting, which is vital for improving vegetable quality.

Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L

BACKGROUND

Brassica napus (B. napus) encompasses diverse transcription factors (TFs), but thorough identification and characterization of TF families, as well as their transcriptional responsiveness to multifarious stresses are still not clear.

RESULTS

Totally 2167 TFs belonging to five families were genome-widely identified in B. napus, including 518 BnAP2/EREBPs, 252 BnbZIPs, 721 BnMYBs, 398 BnNACs and 278 BnWRKYs, which contained some novel members in comparison with existing results. Sub-genome distributions of BnAP2/EREBPs and BnMYBs indicated that the two families might have suffered from duplication and divergence during evolution. Synteny analysis revealed strong co-linearity between B. napus and its two ancestors, although chromosomal rearrangements have occurred and 85 TFs were lost. About 7.6% and 9.4% TFs of the five families in B. napus were novel genes and conserved genes, which both showed preference on the C sub-genome. RNA-Seq revealed that more than 80% TFs were abiotic stress inducible and 315 crucial differentially expressed genes (DEGs) were screened out. Network analysis revealed that the 315 DEGs are highly co-expressed. The homologous gene network in A. thaliana revealed that a considerable amount of TFs could trigger the differential expression of targeted genes, resulting in a complex clustered network with clusters of genes responsible for targeted stress responsiveness.

CONCLUSIONS

We identified and characterized five TF families in B. napus. Some crucial members and regulatory networks involved in different abiotic stresses have been explored. The investigations deepen our understanding of TFs for stress tolerance in B. napus.

Transcriptome Analysis in Chinese Cabbage (Brassica rapa ssp. pekinensis) Provides the Role of Glucosinolate Metabolism in Response to Drought Stress

Although drought stress is one of the most limiting factors in growth and production of Chinese cabbage (Brassica rapa L. ssp. pekinensis), the underlying biochemical and molecular causes are poorly understood. In the present study, to address the mechanisms underlying the drought responses, we analyzed the transcriptome profile of Chinese cabbage grown under drought conditions. Drought stress transcriptionally activated several transcription factor genes, including AP2/ERFs, bHLHs, NACs and bZIPs, and was found to possibly result in transcriptional variation in genes involved in organic substance metabolic processes. In addition, comparative expression analysis of selected BrbZIPs under different stress conditions suggested that drought-induced BrbZIPs are important for improving drought tolerance. Further, drought stress in Chinese cabbage caused differential acclimation responses in glucosinolate metabolism in leaves and roots. Analysis of stomatal aperture indicated that drought-induced accumulation of glucosinolates in leaves directly or indirectly controlled stomatal closure to prevent water loss, suggesting that organ-specific responses are essential for plant survival under drought stress condition. Taken together, our results provide information important for further studies on molecular mechanisms of drought tolerance in Chinese cabbage.

Evolutionary and expression analyses of soybean basic Leucine zipper transcription factor family

BACKGROUND

Soybean, a major legume crop native to East Asia, presents a wealth of resources for utilization. The basic leucine zipper (bZIP) transcription factors play important roles in various biological processes including developmental regulation and responses to environmental stress stimuli. Currently, little information is available regarding the bZIP family in the legume crop soybean.

RESULTS

Using a genome-wide domain analysis, we identified 160 GmbZIP genes in soybean genome, named from GmbZIP1 to GmbZIP160. These 160GmbZIP genes, distributed unevenly across 20 chromosomes, were grouped into 12 subfamilies based on phylogenetic analysis. Gene structure and conserved motif analyses showed that GmbZIP within the same subfamily shared similar intron-exon organizations and motif composition. Syntenic and phylogenetic analyses identified 40 Arabidopsis bZIP genes and 83 soybean bZIP genes as orthologs. By investigating the expression profiling of GmbZIP in different tissues and under drought and flooding stresses, we showed that a majority of GmbZIP (83.44%) exhibited transcript abundance in all examined tissues and 75.6% displayed transcript changes after drought and flooding treatment, suggesting that GmbZIP may play a broad role in soybean development and response to water stress.

CONCLUSIONS

One hundred sixty GmbZIP genes were identified in soybean genome. Our results provide insights for the evolutionary history of bZIP family in soybean and shed light on future studies on the function of bZIP genes in response to water stress in soybean.

Two Cassava Basic Leucine Zipper (bZIP) Transcription Factors (MebZIP3 and MebZIP5) Confer Disease Resistance against Cassava Bacterial Blight

Basic domain-leucine zipper (bZIP) transcription factor, one type of conserved gene family, plays an important role in plant development and stress responses. Although 77 MebZIPs have been genome-wide identified in cassava, their in vivo roles remain unknown. In this study, we analyzed the expression pattern and the function of two MebZIPs (MebZIP3 and MebZIP5) in response to pathogen infection. Gene expression analysis indicated that MebZIP3 and MebZIP5 were commonly regulated by flg22, Xanthomonas axonopodis pv. manihotis (Xam), salicylic acid (SA), and hydrogen peroxide (H₂O₂). Subcellular localization analysis showed that MebZIP3 and MebZIP5 are specifically located in cell nucleus. Through overexpression in tobacco, we found that MebZIP3 and MebZIP5 conferred improved disease resistance against cassava bacterial blight, with more callose depositions. On the contrary, MebZIP3- and MebZIP5-silenced plants by virus-induced gene silencing (VIGS) showed disease sensitive phenotype, lower transcript levels of defense-related genes and less callose depositions. Taken together, this study highlights the positive role of MebZIP3 and MebZIP5 in disease resistance against cassava bacterial blight for further utilization in genetic improvement of cassava disease resistance.

De Novo Assembly and Transcriptome Analysis of Bulb Onion (Allium cepa L.) during Cold Acclimation Using Contrasting Genotypes

Bulb onion (Allium cepa) is the second most widely cultivated and consumed vegetable crop in the world. During winter, cold injury can limit the production of bulb onion. Genomic resources available for bulb onion are still very limited. To date, no studies on heritably durable cold and freezing tolerance have been carried out in bulb onion genotypes. We applied high-throughput sequencing technology to cold (2°C), freezing (-5 and -15°C), and control (25°C)-treated samples of cold tolerant (CT) and cold susceptible (CS) genotypes of A. cepa lines. A total of 452 million paired-end reads were de novo assembled into 54,047 genes with an average length of 1,331 bp. Based on similarity searches, these genes were aligned with entries in the public non-redundant (nr) database, as well as KEGG and COG database. Differentially expressed genes (DEGs) were identified using log₁₀ values with the FPKM method. Among 5,167DEGs, 491 genes were differentially expressed at freezing temperature compared to the control temperature in both CT and CS libraries. The DEG results were validated with qRT-PCR. We performed GO and KEGG pathway enrichment analyses of all DEGs and iPath interactive analysis found 31 pathways including those related to metabolism of carbohydrate, nucleotide, energy, cofactors and vitamins, other amino acids and xenobiotics biodegradation. Furthermore, a large number of molecular markers were identified from the assembled genes, including simple sequence repeats (SSRs) 4,437 and SNP substitutions of transition and transversion types of CT and CS. Our study is the first to provide a transcriptome sequence resource for Allium spp. with regard to cold and freezing stress. We identified a large set of genes and determined their DEG profiles under cold and freezing conditions using two different genotypes. These data represent a valuable resource for genetic and genomic studies of Allium spp.

Glutathione Transferases Superfamily: Cold-Inducible Expression of Distinct GST Genes in Brassica oleracea

Plants, as sessile organisms, can suffer serious growth and developmental consequences under cold stress conditions. Glutathione transferases (GSTs, EC 2.5.1.18) are ubiquitous and multifunctional conjugating proteins, which play a major role in stress responses by preventing oxidative damage by reactive oxygen species (ROS). Currently, understanding of their function(s) during different biochemical and signaling pathways under cold stress condition remain unclear. In this study, using combined computational strategy, we identified 65 Brassica oleracea glutathione transferases (BoGST) and characterized them based on evolutionary analysis into 11 classes. Inter-species and intra-species duplication was evident between BoGSTs and Arabidopsis GSTs. Based on localization analyses, we propose possible pathways in which GST genes are involved during cold stress. Further, expression analysis of the predicted putative functions for GST genes were investigated in two cold contrasting genotypes (cold tolerance and susceptible) under cold condition, most of these genes were highly expressed at 6 h and 1 h in the cold tolerant (CT) and cold susceptible (CS) lines, respectively. Overall, BoGSTU19, BoGSTU24, BoGSTF10 are candidate genes highly expressed in B. oleracea. Further investigation of GST superfamily in B. oleracea will aid in understanding complex mechanism underlying cold tolerance in plants.

The phylogeny of C/S1 bZIP transcription factors reveals a shared algal ancestry and the pre-angiosperm translational regulation of S1 transcripts

Basic leucine zippers (bZIPs) form a large plant transcription factor family. C and S1 bZIP groups can heterodimerize, fulfilling crucial roles in seed development and stress response. S1 sequences also harbor a unique regulatory mechanism, termed Sucrose-Induced Repression of Translation (SIRT). The conservation of both C/S1 bZIP interactions and SIRT remains poorly characterized in non-model species, leaving their evolutionary origin uncertain and limiting crop research. In this work, we explored recently published plant sequencing data to establish a detailed phylogeny of C and S1 bZIPs, investigating their intertwined role in plant evolution, and the origin of SIRT. Our analyses clarified C and S1 bZIP orthology relationships in angiosperms, and identified S1 sequences in gymnosperms. We experimentally showed that the gymnosperm orthologs are regulated by SIRT, tracing back the origin of this unique regulatory mechanism to the ancestor of seed plants. Additionally, we discovered an earlier S ortholog in the charophyte algae Klebsormidium flaccidum, together with a C ortholog. This suggests that C and S groups originated by duplication from a single algal proto-C/S ancestor. Based on our observations, we propose a model wherein the C/S1 bZIP dimer network evolved in seed plants from pre-existing C/S bZIP interactions.

Genome-Wide Analysis of the bZIP Gene Family Identifies Two ABI5-Like bZIP Transcription Factors, BrABI5a and BrABI5b, as Positive Modulators of ABA Signalling in Chinese Cabbage

bZIP (basic leucine zipper) transcription factors coordinate plant growth and development and control responses to environmental stimuli. The genome of Chinese cabbage (Brassica rapa) encodes 136 putative bZIP transcription factors. The bZIP transcription factors in Brassica rapa (BrbZIP) are classified into 10 subfamilies. Phylogenetic relationship analysis reveals that subfamily A consists of 23 BrbZIPs. Two BrbZIPs within subfamily A, Bra005287 and Bra017251, display high similarity to ABI5 (ABA Insensitive 5). Expression of subfamily A BrbZIPs, like BrABI5a (Bra005287/BrbZIP14) and BrABI5b (Bra017251/BrbZIP13), are significantly induced by the plant hormone ABA. Subcellular localization assay reveal that both BrABI5a and BrABI5b have a nuclear localization. BrABI5a and BrABI5b could directly stimulate ABA Responsive Element-driven HIS (a HIS3 reporter gene, which confers His prototrophy) or LUC (LUCIFERASE) expression in yeast and Arabidopsis protoplast. Deletion of the bZIP motif abolished BrABI5a and BrABI5b transcriptional activity. The ABA insensitive phenotype of Arabidopsis abi5-1 is completely suppressed in transgenic lines expressing BrABI5a or BrABI5b. Overall, these results suggest that ABI5 orthologs, BrABI5a and BrABI5b, have key roles in ABA signalling in Chinese cabbage.

Analysis of basic leucine zipper genes and their expression during bud dormancy in peach (Prunus persica)

Dormancy is a biological characteristic developed to resist the cold conditions in winter. The bZIP transcription factors are present exclusively in eukaryotes and have been identified and classified in many species. bZIP proteins are known to regulate numerous biological processes, however, the role of bZIP in bud dodormancy has not been studied extensively. In total, 50 PpbZIP transcription factor-encoding genes were identified and categorized them into 10 groups (A-I and S). Similar intron/exon structures, additional conserved motifs, and DNA-binding site specificity supported our classification scheme. Additionally, chromosomal distribution and collinearity analyses suggested that expansion of the PpbZIP transcription factor family was due to segment/chromosomal duplications. We also predicted the dimerization properties based on characteristic features of the leucine zipper and classified PpbZIP proteins into 23 subfamilies. Furthermore, qRT-PCR results indicated that PpbZIPs genes may be involved in regulating dormancy. The same gene of different species might participate in different regulating networks through interactions with specific partners. Our expression profiling results complemented the microarray data, suggesting that co-expression patterns of bZIP transcription factors during dormancy differed among deciduous fruit trees. Our findings further clarify the molecular characteristics of the PpbZIP transcription factor family, including potential gene functions during dormancy. This information may facilitate further research on the evolutionary history and biological functions of bZIP proteins in peach and other rosaceae plants.

Genome-Wide Identification and Characterization of bZIP Transcription Factors in Brassica oleracea under Cold Stress

Cabbages (Brassica oleracea L.) are an important vegetable crop around world, and cold temperature is among the most significant abiotic stresses causing agricultural losses, especially in cabbage crops. Plant bZIP transcription factors play diverse roles in biotic/abiotic stress responses. In this study, 119 putative BolbZIP transcription factors were identified using amino acid sequences from several bZIP domain consensus sequences. The BolbZIP members were classified into 63 categories based on amino acid sequence similarity and were also compared with BrbZIP and AtbZIP transcription factors. Based on this BolbZIP identification and classification, cold stress-responsive BolbZIP genes were screened in inbred lines, BN106 and BN107, using RNA sequencing data and qRT-PCR. The expression level of the 3 genes, Bol008071, Bol033132, and Bol042729, was significantly increased in BN107 under cold conditions and was unchanged in BN106. The upregulation of these genes in BN107, a cold-susceptible inbred line, suggests that they might be significant components in the cold response. Among three identified genes, Bol033132 has 97% sequence similarity to Bra020735, which was identified in a screen for cold-related genes in B. rapa and a protein containing N-rich regions in LCRs. The results obtained in this study provide valuable information for understanding the potential function of BolbZIP transcription factors in cold stress responses.

Genome-wide characterization and analysis of bZIP transcription factor gene family related to abiotic stress in cassava

The basic leucine zipper (bZIP) transcription factor family plays crucial roles in various aspects of biological processes. Currently, no information is available regarding the bZIP family in the important tropical crop cassava. Herein, 77 bZIP genes were identified from cassava. Evolutionary analysis indicated that MebZIPs could be divided into 10 subfamilies, which was further supported by conserved motif and gene structure analyses. Global expression analysis suggested that MebZIPs showed similar or distinct expression patterns in different tissues between cultivated variety and wild subspecies. Transcriptome analysis of three cassava genotypes revealed that many MebZIP genes were activated by drought in the root of W14 subspecies, indicating the involvement of these genes in the strong resistance of cassava to drought. Expression analysis of selected MebZIP genes in response to osmotic, salt, cold, ABA, and H₂O₂ suggested that they might participate in distinct signaling pathways. Our systematic analysis of MebZIPs reveals constitutive, tissue-specific and abiotic stress-responsive candidate MebZIP genes for further functional characterization in planta, yields new insights into transcriptional regulation of MebZIP genes, and lays a foundation for understanding of bZIP-mediated abiotic stress response.

Global Gene-Expression Analysis to Identify Differentially Expressed Genes Critical for the Heat Stress Response in Brassica rapa

Genome-wide dissection of the heat stress response (HSR) is necessary to overcome problems in crop production caused by global warming. To identify HSR genes, we profiled gene expression in two Chinese cabbage inbred lines with different thermotolerances, Chiifu and Kenshin. Many genes exhibited >2-fold changes in expression upon exposure to 0.5– 4 h at 45°C (high temperature, HT): 5.2% (2,142 genes) in Chiifu and 3.7% (1,535 genes) in Kenshin. The most enriched GO (Gene Ontology) items included ‘response to heat’, ‘response to reactive oxygen species (ROS)’, ‘response to temperature stimulus’, ‘response to abiotic stimulus’, and ‘MAPKKK cascade’. In both lines, the genes most highly induced by HT encoded small heat shock proteins (Hsps) and heat shock factor (Hsf)-like proteins such as HsfB2A (Bra029292), whereas high-molecular weight Hsps were constitutively expressed. Other upstream HSR components were also up-regulated: ROS-scavenging genes like glutathione peroxidase 2 (BrGPX2, Bra022853), protein kinases, and phosphatases. Among heat stress (HS) marker genes in Arabidopsis, only exportin 1A (XPO1A) (Bra008580, Bra006382) can be applied to B. rapa for basal thermotolerance (BT) and short-term acquired thermotolerance (SAT) gene. CYP707A3 (Bra025083, Bra021965), which is involved in the dehydration response in Arabidopsis, was associated with membrane leakage in both lines following HS. Although many transcription factors (TF) genes, including DREB2A (Bra005852), were involved in HS tolerance in both lines, Bra024224 (MYB41) and Bra021735 (a bZIP/AIR1 [Anthocyanin-Impaired-Response-1]) were specific to Kenshin. Several candidate TFs involved in thermotolerance were confirmed as HSR genes by real-time PCR, and these assignments were further supported by promoter analysis. Although some of our findings are similar to those obtained using other plant species, clear differences in Brassica rapa reveal a distinct HSR in this species. Our data could also provide a springboard for developing molecular markers of HS and for engineering HS tolerant B. rapa.

Phosphoproteomic analysis of the response of maize leaves to drought, heat and their combination stress

Drought and heat stress, especially their combination, greatly affect crop production. Many studies have described transcriptome, proteome and phosphoproteome changes in response of plants to drought or heat stress. However, the study about the phosphoproteomic changes in response of crops to the combination stress is scare. To understand the mechanism of maize responses to the drought and heat combination stress, phosphoproteomic analysis was performed on maize leaves by using multiplex iTRAQ-based quantitative proteomic and LC-MS/MS methods. Five-leaf-stage maize was subjected to drought, heat or their combination, and the leaves were collected. Globally, heat, drought and the combined stress significantly changed the phosphorylation levels of 172, 149, and 144 phosphopeptides, respectively. These phosphopeptides corresponded to 282 proteins. Among them, 23 only responded to the combined stress and could not be predicted from their responses to single stressors; 30 and 75 only responded to drought and heat, respectively. Notably, 19 proteins were phosphorylated on different sites in response to the single and combination stresses. Of the seven significantly enriched phosphorylation motifs identified, two were common for all stresses, two were common for heat and the combined stress, and one was specific to the combined stress. The signaling pathways in which the phosphoproteins were involved clearly differed among the three stresses. Functional characterization of the phosphoproteins and the pathways identified here could lead to new targets for the enhancement of crop stress tolerance, which will be particularly important in the face of climate change and the increasing prevalence of abiotic stressors.

Phosphoproteomic analysis of the response of maize leaves to drought, heat and their combination stress

Drought and heat stress, especially their combination, greatly affect crop production. Many studies have described transcriptome, proteome and phosphoproteome changes in response of plants to drought or heat stress. However, the study about the phosphoproteomic changes in response of crops to the combination stress is scare. To understand the mechanism of maize responses to the drought and heat combination stress, phosphoproteomic analysis was performed on maize leaves by using multiplex iTRAQ-based quantitative proteomic and LC-MS/MS methods. Five-leaf-stage maize was subjected to drought, heat or their combination, and the leaves were collected. Globally, heat, drought and the combined stress significantly changed the phosphorylation levels of 172, 149, and 144 phosphopeptides, respectively. These phosphopeptides corresponded to 282 proteins. Among them, 23 only responded to the combined stress and could not be predicted from their responses to single stressors; 30 and 75 only responded to drought and heat, respectively. Notably, 19 proteins were phosphorylated on different sites in response to the single and combination stresses. Of the seven significantly enriched phosphorylation motifs identified, two were common for all stresses, two were common for heat and the combined stress, and one was specific to the combined stress. The signaling pathways in which the phosphoproteins were involved clearly differed among the three stresses. Functional characterization of the phosphoproteins and the pathways identified here could lead to new targets for the enhancement of crop stress tolerance, which will be particularly important in the face of climate change and the increasing prevalence of abiotic stressors.

Transcriptome analysis of wheat inoculated with Fusarium graminearum

Plants are frequently exposed to microorganisms like fungi, bacteria, and viruses that cause biotic stresses. Fusarium head blight (FHB) is an economically risky wheat disease, which occurs upon Fusarium graminearum (Fg) infection. Moderately susceptible (cv. "Mizrak 98") and susceptible (cv. "Gun 91") winter type bread wheat cultivars were subjected to transcriptional profiling after exposure to Fg infection. To examine the early response to the pathogen in wheat, we measured gene expression alterations in mock and pathogen inoculated root crown of moderately susceptible (MS) and susceptible cultivars at 12 hours after inoculation (hai) using 12X135K microarray chip. The transcriptome analyses revealed that out of 39,179 transcripts, 3668 genes in microarray were significantly regulated at least in one time comparison. The majority of differentially regulated transcripts were associated with disease response and the gene expression mechanism. When the cultivars were compared, a number of transcripts and expression alterations varied within the cultivars. Especially membrane related transcripts were detected as differentially expressed. Moreover, diverse transcription factors showed significant fold change values among the cultivars. This study presented new insights to understand the early response of selected cultivars to the Fg at 12 hai. Through the KEGG analysis, we observed that the most altered transcripts were associated with starch and sucrose metabolism and gluconeogenesis pathways.