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

Genome-wide expression analysis of LACS gene family implies GhLACS25 functional responding to salt stress in cotton

BACKGROUND

Long-chain acyl-coenzyme A synthetase (LACS) is a type of acylating enzyme with AMP-binding, playing an important role in the growth, development, and stress response processes of plants.

RESULTS

The research team identified different numbers of LACS in four cotton species (Gossypium hirsutum, Gossypium barbadense, Gossypium raimondii, and Gossypium arboreum). By analyzing the structure and evolutionary characteristics of the LACS, the GhLACS were divided into six subgroups, and a chromosome distribution map of the family members was drawn, providing a basis for further research classification and positioning. Promoter cis-acting element analysis showed that most GhLACS contain plant hormones (GA, MeJA) or non-biological stress-related cis-elements. The expression patterns of GhLACS under salt stress treatment were analyzed, and the results showed that GhLACS may significantly participate in salt stress response through different mechanisms. The research team selected 12 GhLACSs responsive to salt stress for tissue expression analysis and found that these genes are expressed in different tissues.

CONCLUSIONS

There is a certain diversity of LACS among different cotton species. Analysis of promoter cis-acting elements suggests that GhLACS may be involved in regulating plant growth, development and stress response processes. GhLACS25 was selected for in-depth study, which confirmed its significant role in salt stress response through virus-induced gene silencing (VIGS) and induced expression in yeast cells.

Genome-Wide Identification of Calmodulin-Binding Protein 60 Gene Family and the Function of GhCBP60B in Cotton Growth and Development and Abiotic Stress Response

The calmodulin-binding protein 60 (CBP60) family is a gene family unique to plants, and its members play a crucial role in plant defense responses to pathogens and growth and development. Considering that cotton is the primary source of natural cotton textile fiber, the functional study of its CBP60 gene family members is critical. In this research, we successfully identified 162 CBP60 members from the genomes of 21 species. Of these, 72 members were found in four cotton species, divided into four clades. To understand the function of GhCBP60B in cotton in depth, we conducted a detailed analysis of its sequence, structure, cis-acting elements, and expression patterns. Research results show that GhCBP60B is located in the nucleus and plays a crucial role in cotton growth and development and response to salt and drought stress. After using VIGS (virus-induced gene silencing) technology to conduct gene silencing experiments, we found that the plants silenced by GhCBP60B showed dwarf plants and shortened stem nodes, and the expression of related immune genes also changed. In further abiotic stress treatment experiments, we found that GhCBP60B-silenced plants were more sensitive to drought and salt stress, and their POD (peroxidase) activity was also significantly reduced. These results imply the vital role of GhCBP60B in cotton, especially in regulating plant responses to drought and salt stress. This study systematically analyzed CBP60 gene family members through bioinformatics methods and explored in depth the biological function of GhCBP60B in cotton. These research results lay a solid foundation for the future use of the GhCBP60B gene to improve cotton plant type and its drought and salt resistance.

Role of bZIP Transcription Factors in Response to NaCl Stress in Tamarix ramosissima under Exogenous Potassium (K+)

Salt stress is a significant environmental factor affecting plant growth and development, with NaCl stress being one of the most common types of salt stress. The halophyte, Tamarix ramosissima Ledeb (T. ramosissima), is frequently utilized for the afforestation of saline-alkali soils. Indeed, there has been limited research and reports by experts and scholars on the regulatory mechanisms of basic leucine zipper (bZIP) genes in T. ramosissima when treated with exogenous potassium (K+) to alleviate the effects of NaCl stress. This study focused on the bZIP genes in T. ramosissima roots under NaCl stress with additional KCl applied. We identified key candidate genes and metabolic pathways related to bZIP and validated them through quantitative real-time PCR (qRT-PCR). The results revealed that under NaCl stress with additional KCl applied treatments at 0 h, 48 h, and 168 h, based on Pfam protein domain prediction and physicochemical property analysis, we identified 20 related bZIP genes. Notably, four bZIP genes (bZIP_2, bZIP_6, bZIP_16, and bZIP_18) were labeled with the plant hormone signal transduction pathway, showing a predominant up-regulation in expression levels. The results suggest that these genes may mediate multiple physiological pathways under NaCl stress with additional KCl applied at 48 h and 168 h, enhancing signal transduction, reducing the accumulation of ROS, and decreasing oxidative damage, thereby enhancing the tolerance of T. ramosissima to NaCl stress. This study provides gene resources and a theoretical basis for further breeding of salt-tolerant Tamarix species and the involvement of bZIP transcription factors in mitigating NaCl toxicity.

Genome-wide identification of ETHYLENE INSENSITIVE 2 in Triticeae species reveals that TaEIN2-4D.1 regulates cadmium tolerance in Triticum aestivum

ETHYLENE INSENSITIVE 2 (EIN2), as the core component of the ethylene signaling pathway, can widely regulate plant growth, development, and stress responses. However, the comprehensive study and function of EIN2 in wheat Cadmium (Cd) stress remain largely unexplored. Here, we identified 33 EIN2 genes and designated as TaEIN2-2B to TaEIN2-Un.3 in Triticum aestivum. The analysis of cis-regulatory elements in promoter regions and RNA-Seq showed that TaEIN2s were functionally related to plant growth and development, as well as the response to biotic and abiotic stress. qRT-PCR analysis of TaEIN2s indicated their sensitivity to Cd stress. Compared with WT plants, TaEIN2-4D.1-RNAi transgenic wheat lines showed enhanced shoot and root elongation, dry weight and chlorophyll accumulation, together with a reduced accumulation of Cd in wheat grain. In addition, TaEIN2-4D.1-RNAi transgenic wheat lines showed enhanced Reactive Oxygen Species (ROS) scavenging capacity compared with WT plants. In conclusion, our research indicates that TaEIN2 plays a key role in response to cadmium stress in wheat, which provides valuable information for crop improvement.

Differential Gene Expression Responses to Salt and Drought Stress in Tall Fescue (Festuca arundinacea Schreb.)

Understanding gene expression kinetics and the underlying physiological mechanisms in stress combinations is a challenge for the purpose of stress resistance breeding. The novelty of this study is correlating the physiological mechanisms with the expression of key target genes in tall fescue under a combination of various salinity and osmotic stress treatments. Four drought- and salt-responsive genes belonging to different crucial pathways evaluated included one transcription factor FabZIP69, one for the cytosolic polyamine synthetase FaADC1, one for ABA signaling FaCYP707A1, and another one for the specific Na+/H+ plasma membrane antiporter FaSOS1 involve in osmotic homeostasis. FaSOS1, FaCYP707A1, and FabZIP69 were induced early at 6 h after NaCl treatment, while FaSOS1 and FaCYP707A1 were transcribed gradually after exposure to PEG. However, stress interactions showed a significantly increased expression in all genes. Expression of these genes was positively correlated to Pro, SSs, IL, DPPH, and antioxidant enzyme activity and negatively correlated with RWC, total Chl, and MSI. Chemical analyses showed that tall fescue plants exposed to the combination of stresses exhibited increased quantity of reactive oxygen species (H2O2), EL and DPPH, and higher levels of antioxidant enzyme activities (CAT, and SOD), Pro, and SSs content, compared with control seedlings. Under dual-stress conditions, the expression of FabZIP69 was effective in controlling the expression of FaSOS1 and FaADC1 genes differently.

GhCKX14 responding to drought stress by modulating antioxi-dative enzyme activity in Gossypium hirsutum compared to CKX family genes

BACKGROUND

Cytokinin oxidase/dehydrogenase (CKX) plays a vital role in response to abiotic stress through modulating the antioxidant enzyme activities. Nevertheless, the biological function of the CKX gene family has yet to be reported in cotton.

RESULT

In this study, a total of 27 GhCKXs were identified by the genome-wide investigation and distributed across 18 chromosomes. Phylogenetic tree analysis revealed that CKX genes were clustered into four clades, and most gene expansions originated from segmental duplications. The CKXs gene structure and motif analysis displayed remarkably well conserved among the four groups. Moreover, the cis-acting elements related to the abiotic stress, hormones, and light response were identified within the promoter regions of GhCKXs. Transcriptome data and RT-qPCR showed that GhCKX genes demonstrated higher expression levels in various tissues and were involved in cotton's abiotic stress and phytohormone response. The protein-protein interaction network indicates that the CKX family probably participated in redox regulation, including oxidoreduction or ATP levels, to mediate plant growth and development. Functionally identified via virus-induced gene silencing (VIGS) found that the GhCKX14 gene improved drought resistance by modulating the antioxidant-related activitie.

CONCLUSIONS

In this study, the CKX gene family members were analyzed by bioinformatics, and validates the response of GhCKX gene to various phytohormone treatment and abiotic stresses. Our findings established the foundation of GhCKXs in responding to abiotic stress and GhCKX14 in regulating drought resistance in cotton.

Genome-Wide Identification and Functional Analysis of RF2 Gene Family and the Critical Role of GhRF2-32 in Response to Drought Stress in Cotton

Cotton is an important natural fiber crop. The RF2 gene family is a member of the bZIP transcription factor superfamily, which plays an important role in plant resistance to environmental stresses. In this paper, the RF2 gene family of four cotton species was analyzed genome-wide, and the key gene RF2-32 was cloned for functional verification. A total of 113 RF2 genes were identified in the four cotton species, and the RF2 family was relatively conserved during the evolution of cotton. Chromosome mapping and collinear analysis indicated that fragment replication was the main expansion mode of RF2 gene family during evolution. Cis-element analysis showed that there were many elements related to light response, hormone response and abiotic stress response in the promoters of RF2 genes. The transcriptome and qRT-PCR analysis of RF2 family genes in upland cotton showed that RF2 family genes responded to salt stress and drought stress. GhRF2-32 protein was localized in the cell nucleus. Silencing the GhRF2-32 gene showed less leaf wilting and increased total antioxidant capacity under drought and salt stress, decreased malondialdehyde content and increased drought and salt tolerance. This study revealed the evolutionary and functional diversity of the RF2 gene family, which laid a foundation for the further study of stress-resistant genes in cotton.

Identification of salt stress-tolerant candidate genes in the BC2F2 population at the seedling stages of G. hirsutum and G. darwinii using NGS-based bulked segregant analysis

Salinity is a major threat to the yield and productivity of cotton seedlings. In the present study, we developed a BC₂F₂ population of cotton plants from Gossypium darwinii (5-7) and Gossypium hirsutum (CCRI 12-4) salt-susceptible parents to identify salt-resistant candidate genes. The Illumina HiSeq™ strategy was used with bulked segregant analysis. Salt-resistant and salt-susceptible DNA bulks were pooled by using 30 plants from a BC₂F₂ population. Next-generation sequencing (NGS) technology was used for the sequencing of parents and both bulks. Four significant genomic regions were identified: the first genomic region was located on chromosome 18 (1.86 Mb), the second and third genomic regions were on chromosome 25 (1.06 Mb and 1.94 Mb, respectively), and the fourth was on chromosome 8 (1.41 Mb). The reads of bulk1 and bulk2 were aligned to the G. darwinii and G. hirsutum genomes, respectively, leading to the identification of 20,664,007 single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels). After the screening, 6,573 polymorphic markers were obtained after filtration of the candidate regions. The SNP indices in resistant and susceptible bulks and Δ(SNP-index) values of resistant and susceptible bulks were measured. Based on the higher Δ(SNP-index) value, six effective polymorphic SNPs were selected in a different chromosome. Six effective SNPs were linked to five candidate genes in four genomic regions. Further validation of these five candidate genes was carried out using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), resulting in an expression profile that showed two highly upregulated genes in the salt-tolerant species G. darwinii, i.e., Gohir.D05G367800 and Gohir.D12G239100; however, the opposite was shown in G. hirsutum, for which all genes, except one, showed partial expression. The results indicated that Gohir.D05G367800 and Gohir.D12G239100 may be salt-tolerant genes. We are confident that this study could be helpful for the cloning, transformation, and development of salt-resistant cotton varieties.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

Comprehensive genomic identification of cotton starch synthase genes reveals that GhSS9 regulates drought tolerance

Introduction

Starch metabolism is involved in the stress response. Starch synthase (SS) is the key enzyme in plant starch synthesis, which plays an indispensable role in the conversion of pyrophosphoric acid to starch. However, the SS gene family in cotton has not been comprehensively identified and systematically analyzed.

Result

In our study, a total of 76 SS genes were identified from four cotton genomes and divided into five subfamilies through phylogenetic analysis. Genetic structure analysis proved that SS genes from the same subfamily had similar genetic structure and conserved sequences. A cis-element analysis of the SS gene promoter showed that it mainly contains light response elements, plant hormone response elements, and abiotic stress elements, which indicated that the SS gene played key roles not only in starch synthesis but also in abiotic stress response. Furthermore, we also conducted a gene interaction network for SS proteins. Silencing GhSS9 expression decreased the resistance of cotton to drought stress. These findings suggested that SS genes could be related to drought stress in cotton, which provided theoretical support for further research on the regulation mechanism of SS genes on abiotic starch synthesis and sugar levels.

GhAAO2 was observed responding to NaHCO3 stress in cotton compared to AAO family genes

BACKGROUND

Abscisic acid (ABA) is an important stress hormone, the changes of abscisic acid content can alter plant tolerance to stress, abscisic acid is crucial for studying plant responses to abiotic stress. The abscisic acid aldehyde oxidase (AAO) plays a vital role in the final step in the synthesis of abscisic acid, therefore, understanding the function of AAO gene family is of great significance for plants to response to abiotic stresses.

RESULT

In this study, 6, 8, 4 and 4 AAO genes were identified in four cotton species. According to the structural characteristics of genes and the traits of phylogenetic tree, we divided the AAO gene family into 4 clades. Gene structure analysis showed that the AAO gene family was relatively conservative. The analysis of cis-elements showed that most AAO genes contained cis-elements related to light response and plant hormones. Tissue specificity analysis under NaHCO₃ stress showed that GhAAO2 gene was differentially expressed in both roots and leaves. After GhAAO2 gene silencing, the degree of wilting of seedlings was lighter than that of the control group, indicating that GhAAO2 could respond to NaHCO₃ stress.

CONCLUSIONS

In this study, the AAO gene family was analyzed by bioinformatics, the response of GhAAO gene to various abiotic stresses was preliminarily verified, and the function of the specifically expressed gene GhAAO2 was further verified. These findings provide valuable information for the study of potential candidate genes related to plant growth and stress.

Evolutionary Relationships and Divergence of Filamin Gene Family Involved in Development and Stress in Cotton (Gossypium hirsutum L.)

Filamin protein is characterized by an N-terminal actin-binding domain that is followed by 24 Ig (immunoglobulin)-like repeats, which act as hubs for interactions with a variety of proteins. In humans, this family has been found to be involved in cancer cell invasion and metastasis and can be involved in a variety of growth signal transduction processes, but it is less studied in plants. Therefore, in this study, 54 Filamin gene family members from 23 plant species were investigated and divided into two subfamilies: FLMN and GEX2. Subcellular localization showed that most of the Filamin gene family members were located in the cell membrane. A total of 47 Filamin gene pairs were identified, most of which were whole-genome copies. Through the analyses of cis-acting elements, expression patterns and quantitative fluorescence, it was found that GH_ A02G0519 and GH_ D02G0539 are mainly expressed in the reproductive organs of upland cotton, and their interacting proteins are also related to the fertilization process, whereas GH_A02G0216 and GH_D02G0235 were related to stress. Thus, it is speculated that two genes of the GEX2 subfamily (GH_A02G0519 and GH_D02G0539) may be involved in the reproductive development of cotton and may affect the fertilization process of cotton. This study provides a theoretical basis for the further study of the cotton Filamin gene family.

Identification and Functional Analysis of bZIP Genes in Cotton Response to Drought Stress

The basic leucine zipper (bZIP) transcription factors, which harbor a conserved bZIP domain composed of two regions, a DNA-binding basic region and a Leu Zipper region, operate as important switches of transcription networks in eukaryotes. However, this gene family has not been systematically characterized in cotton (Gossypium hirsutum). Here, we identified 197 bZIP family members in cotton. The chromosome distribution pattern indicates that the GhbZIP genes have undergone 53 genome-wide segmental and 7 tandem duplication events which contribute to the expansion of the cotton bZIP family. Phylogenetic analysis showed that cotton GhbZIP proteins cluster into 13 subfamilies, and homologous protein pairs showed similar characteristics. Inspection of the DNA-binding basic region and leucine repeat heptads within the bZIP domains indicated different DNA-binding site specificities as well as dimerization properties among different groups. Comprehensive expression analysis indicated the most highly and differentially expressed genes in root and leaf that might play significant roles in cotton response to drought stress. GhABF3D was identified as a highly and differentially expressed bZIP family gene in cotton leaf and root under drought stress treatments that likely controls drought stress responses in cotton. These data provide useful information for further functional analysis of the GhbZIP gene family and its potential application in crop improvement.

Genome-wide identification of GAD family genes suggests GhGAD6 functionally respond to Cd2+ stress in cotton

Glutamate decarboxylase (GAD) mainly regulated the biosynthesis of γ-aminobutyric acid (GABA) and played an important role in plant growth and stress resistance. To explore the potential function of GAD in cotton growth, the genome-wide identification, structure, and expression analysis of GAD genes were performed in this study. There were 10, 9, 5, and 5 GAD genes identified in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. GAD was divided into four clades according to the protein motif composition, gene structure, and phylogenetic relationship. The segmental duplication was the main way of the GAD gene family evolution. Most GhGADs respond to abiotic stress. Clade Ⅲ GAD was induced by Cd²⁺ stress, especially GhGAD6, and silencing GhGAD6 would lead to more serious Cd²⁺ poisoning in cotton. The oxidative damage caused by Cd²⁺ stress was relieved by increasing the GABA content. It was speculated that the decreased expression of GhGAD6 reduced the content of GABA in vivo and caused the accumulation of ROS. This study will further expand our understanding of the relationship between the evolution and function of the GhGAD gene family and provide new genetic resources for cotton breeding under environmental stress and phytoremediation.

Genome-wide identification, evolution and function analysis of UGTs superfamily in cotton

Glycosyltransferases mainly catalyse the glycosylation reaction in living organisms and widely exists in plants. UGTs have been identified from G. raimondii, G. arboreum and G. hirsutum. However, Genome-wide systematic analysis of UGTs superfamily have not been studied in G. barbadense. 752 UGTs were identified from four cotton species and grouped into 18 clades, of which R was newly discovered clades. Most UGTs were clustered at both ends of the chromosome and showed a heterogeneous distribution. UGT proteins were widely distributed in cells, with the highest distribution in chloroplasts. UGTs of the same clade shared similar intron/exon structural features. During evolution, the gene family has undergone strong selection for purification. UGTs were significantly enriched in “transcriptional activity (GO:0016758)” and “metabolic processes (GO:0008152)”. Genes from the same clade differed in function under various abiotic stresses. The analysis of cis-acting element and qRT–PCR may indicate that GHUGTs play important roles in plant growth, development and abiotic stress. We further found that GHUGT74-2 plays an important role under submergence. The study broadens the understanding of UGTs in terms of gene characteristics, evolutionary processes, and gene function in cotton and provides a new way to systematically and globally understand the structure–function relationship of multigene families in the evolutionary process.

Genome-Wide Expression Profiling and Networking Reveals an Imperative Role of IMF-Associated Novel CircRNAs as ceRNA in Pigs

Intramuscular fat (IMF) deposition is a biological process that has a strong impact on the nutritional and sensorial properties of meat, with relevant consequences on human health. Pork loins determine the effects of marbling on the sensory attributes and meat quality properties, which differ among various pig breeds. This study explores the crosstalk of non-coding RNAs with mRNAs and analyzes the potential pathogenic role of IMF-associated competing endogenous RNA (ceRNA) in IMF tissues, which offer a framework for the functional validation of key/potential genes. A high-throughput whole-genome transcriptome analysis of IMF tissues from longissimus dorsi muscles of Large White (D_JN) and Laiwu (L_JN) pigs resulted in the identification of 283 differentially expressed circRNAs (DECs), including two key circRNAs (circRNA-23437, circRNA-08840) with potential binding sites for multiple miRNAs regulating the whole network. The potential ceRNA mechanism identified the DEC target miRNAs-mRNAs involved in lipid metabolism, fat deposition, meat quality, and metabolic syndrome via the circRNA-miRNA-mRNA network, concluding that ssc-mir-370 is the most important target miRNA shared by both key circRNAs. TGM2, SLC5A6, ECI1, FASN, PER1, SLC25A34, SOD1, and COL5A3 were identified as hub genes through an intensive protein-protein interaction (PPI) network analysis of target genes acquired from the ceRNA regulatory network. Functional enrichments, pathway examinations, and qRT-PCR analyses infer their implications in fat/cholesterol metabolism, insulin secretion, and fatty acid biosynthesis. Here, circRNAs and miRNA sequencing accompanied by computational techniques were performed to analyze their expressions in IMF tissues from the longissimus dorsi muscles of two pig breeds. Their target gene evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations, and structural advances with high-throughput protein modeling, following genomic organizations, will provide new insights into the underlying molecular mechanisms of adipocyte differentiation and IMF deposition and a much-needed qualitative framework for future research to improve meat quality and its role as a biomarker to treat lipid metabolic syndromes.

Structure and character analysis of cotton response regulator genes family reveals that GhRR7 responses to draught stress

Background

Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS).

Method

From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed.

Result

A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H₂O₂ and other reactive oxygen removal systems.

Conclusion

This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00394-2.

Identification and Expression Analysis of bZIP Members under Abiotic Stress in Mung Bean (Vigna radiata)

The main aim of this study was to identify the bZIP family members in mung bean and explore their expression patterns under several abiotic stresses, with the overarching goal of elucidating their biological functions. Results identified 75 bZIP members in mung bean, which were unevenly distributed in the chromosomes (1-11), and all had a highly conserved bZIP domain. Phylogenetic analysis divided the members into 10 subgroups, with members in the same subgroup having similar structure and motif. The cis-acting elements in the promoter region revealed that most of the bZIP members might have the connection with abscisic acid, ethylene, and stress responsive elements. The transcriptome data demonstrated that bZIP members could respond to salt stress at different degrees in leaves, but the expression patterns could vary at different time points under stress. Differentially expressed genes (DEGs), such as VrbZIP12, VrbZIP37, and VrZIP45, were annotated into the plant hormone signal transduction pathway, which might be regulated the expression of abiotic stress-related gene (ABF). Quantitative real-time polymerase chain reaction (qRT-PCR) was applied to determine the expression of bZIP members in roots and leaves under drought, alkali, and low-temperature stress. Results showed that bZIP members respond differently to diverse stresses, and their expression was tissue-specific, which suggests that they may have different regulatory mechanism in different tissues. Overall, this study will provide a reference for further research on the functions of bZIP members in mung bean.

Genome-wide expression analysis of carboxylesterase (CXE) gene family implies GBCXE49 functional responding to alkaline stress in cotton

BACKGROUND

Carboxylesterase (CXE) is a type of hydrolase with α/β sheet hydrolase activity widely found in animals, plants and microorganisms, which plays an important role in plant growth, development and resistance to stress.

RESULTS

A total of 72, 74, 39, 38 CXE genes were identified in Gossypium barbadense, Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively. The gene structure and expression pattern were analyzed. The GBCXE genes were divided into 6 subgroups, and the chromosome distribution of members of the family were mapped. Analysis of promoter cis-acting elements showed that most GBCXE genes contain cis-elements related to plant hormones (GA, IAA) or abiotic stress. These 6 genes we screened out were expressed in the root, stem and leaf tissues. Combined with the heat map, GBCXE49 gene was selected for subcellular locate and confirmed that the protein was expressed in the cytoplasm.

CONCLUSIONS

The collinearity analysis of the CXE genes of the four cotton species in this family indicated that tandem replication played an indispensable role in the evolution of the CXE gene family. The expression patterns of GBCXE gene under different stress treatments indicated that GBCXE gene may significantly participate in the response to salt and alkaline stress through different mechanisms. Through the virus-induced gene silencing technology (VIGS), it was speculated that GBCXE49 gene was involved in the response to alkaline stress in G. barbadense.

Molecular insights into sensing, regulation and improving of heat tolerance in plants

Climate-change-mediated increase in temperature extremes has become a threat to plant productivity. Heat stress-induced changes in growth pattern, sensitivity to pests, plant phonologies, flowering, shrinkage of maturity period, grain filling, and increased senescence result in significant yield losses. Heat stress triggers multitude of cellular, physiological and molecular responses in plants beginning from the early sensing followed by signal transduction, osmolyte synthesis, antioxidant defense, and heat stress-associated gene expression. Several genes and metabolites involved in heat perception and in the adaptation response have been isolated and characterized in plants. Heat stress responses are also regulated by the heat stress transcription factors (HSFs), miRNAs and transcriptional factors which together form another layer of regulatory circuit. With the availability of functionally validated candidate genes, transgenic approaches have been applied for developing heat-tolerant transgenic maize, tobacco and sweet potato. In this review, we present an account of molecular mechanisms of heat tolerance and discuss the current developments in genetic manipulation for heat tolerant crops for future sustainable agriculture.

Genome-wide survey of the dehydrin genes in bread wheat (Triticum aestivum L.) and its relatives: identification, evolution and expression profiling under various abiotic stresses

BACKGROUND

Bread wheat (Triticum aestivum) is an important staple cereal grain worldwide. The ever-increasing environmental stress makes it very important to mine stress-resistant genes for wheat breeding programs. Therefore, dehydrin (DHN) genes can be considered primary candidates for such programs, since they respond to multiple stressors.

RESULTS

In this study, we performed a genome-wide analysis of the DHN gene family in the genomes of wheat and its three relatives. We found 55 DHN genes in T. aestivum, 31 in T. dicoccoides, 15 in T. urartu, and 16 in Aegilops tauschii. The phylogenetic, synteny, and sequence analyses showed we can divide the DHN genes into five groups. Genes in the same group shared similar conserved motifs and potential function. The tandem TaDHN genes responded strongly to drought, cold, and high salinity stresses, while the non-tandem genes respond poorly to all stress conditions. According to the interaction network analysis, the cooperation of multiple DHN proteins was vital for plants in combating abiotic stress.

CONCLUSIONS

Conserved, duplicated DHN genes may be important for wheat being adaptable to a different stress conditions, thus contributing to its worldwide distribution as a staple food. This study not only highlights the role of DHN genes help the Triticeae species against abiotic stresses, but also provides vital information for the future functional studies in these crops.

The Evolution and Expression Profiles of EC1 Gene Family during Development in Cotton

Fertilization is essential to sexual reproduction of flowering plants. EC1 (EGG CELL 1) proteins have a conserved cysteine spacer characteristic and play a crucial role in double fertilization process in many plant species. However, to date, the role of EC1 gene family in cotton is fully unknown. Hence, detailed bioinformatics analysis was explored to elucidate the biological mechanisms of EC1 gene family in cotton. In this study, we identified 66 genes in 10 plant species in which a total of 39 EC1 genes were detected from cotton genome. Phylogenetic analysis clustered the identified EC1 genes into three families (I-III) and all of them contain Prolamin-like domains. A good collinearity was observed in the synteny analysis of the orthologs from cotton genomes. Whole-genome duplication was determined to be one of the major impetuses for the expansion of the EC1 gene family during the process of evolution. qRT-PCR analysis showed that EC1 genes were highly expressed in reproductive tissues under multiple stresses, signifying their potential role in enhancing stress tolerance or responses. Additionally, gene interaction networks showed that EC1 genes may be involved in cell stress and response transcriptional regulator in the synergid cells and activate the expression of genes required for pollen tube guidance. Our results provide novel functional insights into the evolution and functional elucidation of EC1 gene family in cotton.

Genome-wide analysis of the CalS gene family in cotton reveals their potential roles in fiber development and responses to stress

Callose deposition occurs during plant growth and development, as well as when plants are under biotic and abiotic stress. Callose synthase is a key enzyme for the synthesis of callose. In this study, 27, 28, 16, and 15 callose synthase family members were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium raimondii, and Gossypium arboreum using the sequence of Arabidopsis callose synthase. The CalSs were divided into five groups by phylogenetic, gene structure, and conservative motif analysis. The conserved motifs and gene structures of CalSs in each group were highly similar. Based on the analysis of cis-acting elements, it is inferred that GhCalSs were regulated by abiotic stress. WGD/Segmental duplication promoted the amplification of the CalS gene in cotton, and purification selection had an important function in the CalS family. The transcriptome data and qRT-PCR under cold, heat, salt, and PEG treatments showed that GhCalSs were involved in abiotic stress. The expression patterns of GhCalSs were different in various tissues. We predicted that GhCalS4, which was highly expressed in fibers, had an important effect on fiber elongation. Hence, these results help us understand the role of GhCalSs in fiber development and stress response.

Components and Functional Diversification of Florigen Activation Complexes in Cotton

In shoot apex cells of rice, a hexameric florigen activation complex (FAC), comprising flowering locus T (FT), 14-3-3 and the basic leucine zipper transcription factor FD, activates downstream target genes and regulates several developmental transitions, including flowering. The allotetraploid cotton (Gossypium hirsutum L.) contains only one FT locus in both of the A- and D-subgenomes. However, there is limited information regarding cotton FACs. Here, we identified a 14-3-3 protein that interacts strongly with GhFT in the cytoplasm and the nuclei, and five FD homoeologous gene pairs were characterized. In vivo, all five GhFD proteins interacted with Gh14-3-3 and GhFT in the nucleus. GhFT, 14-3-3 and all the GhFDs interacted in the nucleus as well, suggesting that they formed a ternary complex. Virus-induced silencing of GhFD1, -2 and -4 in cotton delayed flowering and inhibited the expression of floral meristem identity genes. Silencing GhFD3 strongly decreased lateral root formation, suggesting a function in lateral root development. GhFD overexpression in Arabidopsis and transcriptional activation assays suggested that FACs containing GhFD1 and GhFD2 function mainly in promoting flowering with partial functional redundancy. Moreover, GhFD3 was specifically expressed in lateral root meristems and dominantly activated the transcription of auxin response factor genes, such as ARF19. Thus, the diverse functions of FACs may depend on the recruited GhFD. Creating targeted genetic mutations in the florigen system using Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins (Cas) genome editing may fine-tune flowering and improve plant architecture.

Genome-wide expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton

Cotton is the most important crop for the production of natural fibres used in the textile industry. High salinity, drought, cold and high temperature represent serious abiotic stresses, which seriously threaten cotton production. Phospholipase A_S has an irreplaceable role in lipid signal transmission, growth and development and stress events. Phospholipase A can be divided into three families: PLA1, PLA2 and pPLA. Among them, the PLA1 family is rarely studied in plants. In order to study the potential functions of the PLA1 family in cotton, the bioinformatics analysis of the PLA1 family was correlated with cotton adversity, and tissue-specific analysis was performed. Explore the structure-function relationship of PLA1 members. It is found that the expression of GbPLA1-32 gene is affected by a variety of environmental stimuli, indicating that it plays a very important role in stress and hormone response, and closely associates the cotton adversity with this family. Through further functional verification, we found that virus-induced GbPLA1-32 gene silencing (VIGS) caused Gossypium barbadense to be sensitive to salt stress. This research provides an important basis for further research on the molecular mechanism of cotton resistance to abiotic stress.

Genome-wide analysis of bZIP, BBR, and BZR transcription factors in Triticum aestivum

Transcription factors are regulatory proteins known to modulate gene expression. These are the critical component of signaling pathways and help in mitigating various developmental and stress responses. Among them, bZIP, BBR, and BZR transcription factor families are well known to play a crucial role in regulating growth, development, and defense responses. However, limited data is available on these transcription factors in Triticum aestivum. In this study, bZIP, BBR, and BZR sequences from Brachypodium distachyon, Oryza sativa, Oryza barthii, Oryza brachyantha, T. aestivum, Triticum urartu, Sorghum bicolor, Zea mays were retrieved, and dendrograms were constructed to analyze the evolutionary relatedness among them. The sequences clustered into one group indicated a degree of evolutionary correlation highlighting the common lineage of cereal grains. This analysis also exhibited that these genes were highly conserved among studied monocots emphasizing their common ancestry. Furthermore, these transcription factor genes were evaluated for envisaging conserved motifs, gene structure, and subcellular localization in T. aestivum. This comprehensive computational analysis has provided an insight into transcription factor evolution that can also be useful in developing approaches for future functional characterization of these genes in T. aestivum. Furthermore, the data generated can be beneficial in future for genetic manipulation of economically important plants.

Divergent Response Strategies of CsABF Facing Abiotic Stress in Tea Plant: Perspectives From Drought-Tolerance Studies

In plants, the bZIP family plays vital roles in various biological processes, including seed maturation, flower development, light signal transduction, pathogen defense, and various stress responses. Tea, as a popular beverage, is widely cultivated and has withstood a degree of environmental adversity. Currently, knowledge of the bZIP gene family in tea plants remains very limited. In this study, a total of 76 CsbZIP genes in tea plant were identified for the whole genome. Phylogenetic analysis with Arabidopsis counterparts revealed that CsbZIP proteins clustered into 13 subgroups, among which 13 ABFs related to the ABA signaling transduction pathway were further identified by conserved motif alignment and named CsABF1-13, these belonged to the A and S subgroups of CsbZIP and had close evolutionary relationships, possessing uniform or similar motif compositions. Transcriptome analysis revealed the expression profiles of CsABF genes in different tissues (bud, young leaf, mature leaf, old leaf, stem, root, flower, and fruit) and under diverse environmental stresses (drought, salt, chilling, and MeJA). Several CsABF genes with relatively low tissue expression, including CsABF1, CsABF5, CsABF9, and CsABF10, showed strong expression induction in stress response. Thirteen CsABF genes, were examined by qRT-PCR in two tea plant cultivars, drought-tolerant "Taicha 12" and drought-sensitive "Fuyun 6", under exogenous ABA and drought stress. Furthermore, CsABF2, CsABF8, and CsABF11, were screened out as key transcription factors regulating drought tolerance of tea cultivars. Subsequently, some potential target genes regulated by CsABFs were screened by co-expression network and enrichment analysis. This study update CsbZIP gene family and provides a global survey of the ABF gene family in tea plant. The resolution of the molecular mechanism of drought resistance in different varieties could be helpful for improving stress resistance in tea plant via genetic engineering.

Genome-Wide Characterization and Expression Analysis of bZIP Gene Family Under Abiotic Stress in Glycyrrhiza uralensis

bZIP gene family is one of the largest transcription factor families. It plays an important role in plant growth, metabolic, and environmental response. However, complete genome-wide investigation of bZIP gene family in Glycyrrhiza uralensis remains unexplained. In this study, 66 putative bZIP genes in the genome of G. uralensis were identified. And their evolutionary classification, physicochemical properties, conserved domain, functional differentiation, and the expression level under different stress conditions were further analyzed. All the members were clustered into 13 subfamilies (A–K, M, and S). A total of 10 conserved motifs were found in GubZIP proteins. Members from the same subfamily shared highly similar gene structures and conserved domains. Tandem duplication events acted as a major driving force for the evolution of bZIP gene family in G. uralensis. Cis-acting elements and protein–protein interaction networks showed that GubZIPs in one subfamily are involved in multiple functions, while some GubZIPs from different subfamilies may share the same functional category. The miRNA network targeting GubZIPs showed that the regulation at the transcriptional level may affect protein–protein interaction networks. We suspected that domain-mediated interactions may categorize a protein family into subfamilies in G. uralensis. Furthermore, the tissue-specific gene expression patterns of GubZIPs were analyzed using the public RNA-seq data. Moreover, gene expression level of 66 bZIP family members under abiotic stress treatments was quantified by using qRT-PCR. The results of this study may serve as potential candidates for functional characterization in the future.

Genome-wide identification and function analysis of HMAD gene family in cotton (Gossypium spp.)

BACKGROUND

The abiotic stress such as soil salinization and heavy metal toxicity has posed a major threat to sustainable crop production worldwide. Previous studies revealed that halophytes were supposed to tolerate other stress including heavy metal toxicity. Though HMAD (heavy-metal-associated domain) was reported to play various important functions in Arabidopsis, little is known in Gossypium.

RESULTS

A total of 169 G. hirsutum genes were identified belonging to the HMAD gene family with the number of amino acids ranged from 56 to 1011. Additionally, 84, 76 and 159 HMAD genes were identified in each G. arboreum, G. raimondii and G. barbadense, respectively. The phylogenetic tree analysis showed that the HMAD gene family were divided into five classes, and 87 orthologs of HMAD genes were identified in four Gossypium species, such as genes Gh_D08G1950 and Gh_A08G2387 of G. hirsutum are orthologs of the Gorai.004G210800.1 and Cotton_A_25987 gene in G. raimondii and G. arboreum, respectively. In addition, 15 genes were lost during evolution. Furthermore, conserved sequence analysis found the conserved catalytic center containing an anion binding (CXXC) box. The HMAD gene family showed a differential expression levels among different tissues and developmental stages in G. hirsutum with the different cis-elements for abiotic stress.

CONCLUSIONS

Current study provided important information about HMAD family genes under salt-stress in Gossypium genome, which would be useful to understand its putative functions in different species of cotton.

Genome-wide identification of CK gene family suggests functional expression pattern against Cd2+ stress in Gossypium hirsutum L

Choline kinase (CK) gene plays an important role in plants growth, development and resistance to stress. It mainly regulates the biosynthesis of phosphatidylcholine. This study aims to explore the structure-function relationship, and to provide a framework for functional validation and biochemical characterization of various CK genes. Our analysis showed that 87 CK genes were identified in cotton and 7 diploid plants, of which 43 genes encode CK proteins in 4 cotton species, and 13 genes were identified in Gossypium hirsutum. Most of GhCK genes are affected by the abiotic stress conditions, indicating the importance of CK proteins for plant development and response to abiotic stress. RT-qPCR analysis showed the tissue specificity of GhCK genes in response to Cd²⁺ and other abiotic stresses. Under Cd²⁺ stress, the expression level of GhCK gene family members has undergone different changes. The expression level of GhCK5 was enhanced, indicating that Cd²⁺ stress caused the increase of phosphatidylcholine content, which in turn reacted on the plant cell membrane, finally reached the absorption of Cd²⁺ into plant cells to repair Cd²⁺ the purpose of contaminated soil. This study will further broaden our understanding of the association between evolution and function of the GhCK gene family.

The Roles of CCCH Zinc-Finger Proteins in Plant Abiotic Stress Tolerance

Zinc-finger proteins, a superfamily of proteins with a typical structural domain that coordinates a zinc ion and binds nucleic acids, participate in the regulation of growth, development, and stress adaptation in plants. Most zinc fingers are C2H2-type or CCCC-type, named after the configuration of cysteine (C) and histidine (H); the less-common CCCH zinc-finger proteins are important in the regulation of plant stress responses. In this review, we introduce the domain structures, classification, and subcellular localization of CCCH zinc-finger proteins in plants and discuss their functions in transcriptional and post-transcriptional regulation via interactions with DNA, RNA, and other proteins. We describe the functions of CCCH zinc-finger proteins in plant development and tolerance to abiotic stresses such as salt, drought, flooding, cold temperatures and oxidative stress. Finally, we summarize the signal transduction pathways and regulatory networks of CCCH zinc-finger proteins in their responses to abiotic stress. CCCH zinc-finger proteins regulate the adaptation of plants to abiotic stress in various ways, but the specific molecular mechanisms need to be further explored, along with other mechanisms such as cytoplasm-to-nucleus shuttling and post-transcriptional regulation. Unraveling the molecular mechanisms by which CCCH zinc-finger proteins improve stress tolerance will facilitate the breeding and genetic engineering of crops with improved traits.

Genomic insight into the divergence and adaptive potential of a forgotten landrace G. hirsutum L. purpurascens

Wild progenitors are an excellent source for strengthening the genetic basis and accumulation of desirable variation lost because of directional selection and adaptation in modern cultivars. Here, we re-evaluate a landrace of Gossypium hirsutum, formerly known as Gossypium purpurascens. Our study seeks to understand the genomic structure, variation, and breeding potential of this landrace, providing potential insights into the biogeographic history and genomic changes likely associated with domestication. A core set of accessions, including current varieties, obsolete accessions, G. purpurascens, and other geographical landraces, are subjected to genotyping along with multilocation phenotyping. Population fixation statistics suggests a marked differentiation between G. purpurascens and three other groups, emphasizing the divergent genomic behavior of G. purpurascens. Phylogenetic analysis establishes the primitive nature of G. purpurascens, identifying it as a vital source of functional variation, the inclusion of which in the upland cotton (cultivated G. hirsutum) gene pool may broaden the genetic basis of modern cultivars. Genome-wide association results indicate multiple loci associated with domestication regions corresponding to flowering and fiber quality. Moreover, the conserved nature of G. purpurascens can also provide insights into the evolutionary process of G. hirsutum.

Cotton DMP gene family: characterization, evolution, and expression profiles during development and stress

Members of DOMAIN OF UNKNOWN FUNCTION 679 membrane protein (DMP) gene family, a type of plant-specific membrane proteins, have been proposed to function in various physiological processes such as reproductive development and senescence in plants. Here, a total of 174 DMP genes were identified and analyzed in 16 plant species (including 58 DMPs in four cotton species). Phylogenetic analysis showed that these DMPs could be clustered into five subfamilies (I-V). 137 duplicated cotton gene pairs were identified and most duplicate events were formed by whole-genome duplication (WGD)/segmental duplications. Expression analysis revealed that most of cotton DMPs were mainly expressed in the reproductive organs (the sepal, petal, pistil and anther) and the fiber of secondary cell wall stage. GhDMPs promoter regions containing the different cis-elements also showed different responses to abiotic stress. In addition, gene interaction networks showed that DMPs, as an endomembrane system, were involved in plant senescence process and flower reproductive development. We speculated GhDMP8-A/-D, GbDMP8-A/-D could be used as some candidate gene for inducing cotton haploid. This genome-wide study provides a systematic analysis of the cotton DMP gene family, and further insights towards understanding the potential functions of candidate genes.

Cotton transcriptome analysis reveals novel biological pathways that eliminate reactive oxygen species (ROS) under sodium bicarbonate (NaHCO3) alkaline stress

Alkaline stress is one of the abiotic stresses limiting cotton production. Though RNA-Seq analyses, have been conducted to investigate genome-wide gene expression in response to alkaline stress in plants, the response of sodium bicarbonate (NaHCO₃) stress-related genes in cotton has not been reported. To explore the mechanisms of cotton response to this alkaline stress, we used next-generation sequencing (NGS) technology to study transcriptional changes of cotton under NaHCO₃ alkaline stress. A total of 18,230 and 11,177 differentially expressed genes (DEGs) were identified in cotton roots and leaves, respectively. Gene ontology (GO) analysis indicated the enrichment of DEGs involved in various stimuli or stress responses. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs associated with plant hormone signal transduction, amino acid biosynthesis, and biosynthesis of secondary metabolites were regulated in response to the NaHCO₃ stress. We further analyzed genes enriched in secondary metabolic pathways and found that secondary metabolites were regulated to eliminate the reactive oxygen species (ROS) and improve the cotton tolerance to the NaHCO₃ stress. In this study, we learned that the toxic effect of NaHCO₃ was more profound than that of NaOH at the same pH. Thus, Na⁺, HCO₃^- and pH had a great impact on the growth of cotton plant. The novel biological pathways and candidate genes for the cotton tolerance to NaHCO₃ stress identified from the study would be useful in the genetic improvement of the alkaline tolerance in cotton.

Regulatory Network of Cotton Genes in Response to Salt, Drought and Wilt Diseases (Verticillium and Fusarium): Progress and Perspective

In environmental conditions, crop plants are extremely affected by multiple abiotic stresses including salinity, drought, heat, and cold, as well as several biotic stresses such as pests and pathogens. However, salinity, drought, and wilt diseases (e.g., Fusarium and Verticillium) are considered the most destructive environmental stresses to cotton plants. These cause severe growth interruption and yield loss of cotton. Since cotton crops are central contributors to total worldwide fiber production, and also important for oilseed crops, it is essential to improve stress tolerant cultivars to secure future sustainable crop production under adverse environments. Plants have evolved complex mechanisms to respond and acclimate to adverse stress conditions at both physiological and molecular levels. Recent progresses in molecular genetics have delivered new insights into the regulatory network system of plant genes, which generally includes defense of cell membranes and proteins, signaling cascades and transcriptional control, and ion uptake and transport and their relevant biochemical pathways and signal factors. In this review, we mainly summarize recent progress concerning several resistance-related genes of cotton plants in response to abiotic (salt and drought) and biotic (Fusarium and Verticillium wilt) stresses and classify them according to their molecular functions to better understand the genetic network. Moreover, this review proposes that studies of stress related genes will advance the security of cotton yield and production under a changing climate and that these genes should be incorporated in the development of cotton tolerant to salt, drought, and fungal wilt diseases (Verticillium and Fusarium).

Identification and expression profiling of MYB transcription factors related to l-theanine biosynthesis in Camellia sinensis

The MYB proteins belong to a large family of transcription factors in plant genomes and play significant roles in primary and secondary metabolism. Although several CsMYB genes have been identified in Camellia sinensis, few CsMYBs involved in l-theanine biosynthesis have been analyzed. In this study, we screened and identified 20 CsMYBs related to l-theanine biosynthesis. Transcriptomic analysis revealed that the expression profiles of the CsMYBs were positively or negatively related to dynamic changes in the l-theanine content. Validation of selected l-theanine biosynthetic and CsMYB genes was conducted by qRT-PCR. The results illustrated that most of the structural and CsMYB genes were downregulated with a decrease in the l-theanine levels. Protein-protein interaction networks of CsMYB5, CsMYB12 and CsMYB94 proteins demonstrated that they might form complexes with bHLH and WD 40 proteins. Multiple DNA-binding sites of the R2R3-MYB protein were observed in promoter regions of structural genes, indicating CsMYB family proteins might be involved in l-theanine metabolism via the attachment of AC elements. Moreover, CsMYB73 demonstrated binding specificity to the promoter region of CsGDH2 (CsGDH2-pro). These findings provide fundamental understanding of specific members of the CsMYBs related to the l-theanine biosynthesis pathway.