A Nck-associated protein 1-like protein affects drought sensitivity by its involvement in leaf epidermal development and stomatal closure in rice

MAPKKK28 functions upstream of the MKK1-MPK1 cascade to regulate abscisic acid responses in rice

The mitogen-activated protein kinase (MAPK) cascade (MAPKKK-MAPKK-MAPK) plays a critical role in biotic and abiotic stress responses and abscisic acid (ABA) signalling. A previous study has shown that the ABA-activated MKK1-MPK1 cascade is essential in regulating ABA response and stress tolerance in rice. However, the specific MAPKKK upstream of the MKK1-MPK1 cascade in ABA signalling remains unknown. Here, we identified that MAPKKK28, a previously uncharacterized member of the rice MEKK family, is involved in regulating ABA responses, including seed germination, root growth, stomatal closure, and the tolerance to oxidative stress and osmotic stress. We found that MAPKKK28 directly interacts with and phosphorylates MKK1. Further analysis indicated that the activation of both MKK1 and MPK1 depends on MAPKKK28 in ABA signalling. Genetic analysis revealed that MAPKKK28 functions upstream of the MKK1-MPK1 cascade to positively regulate ABA responses and enhance tolerance to oxidative and osmotic stress. These results not only reveal a new complete MAPK cascade in plants but also uncover its importance in ABA signalling.

Engineering quantitative stomatal trait variation and local adaptation potential by cis-regulatory editing

Cis-regulatory element editing can generate quantitative trait variation that mitigates extreme phenotypes and harmful pleiotropy associated with coding sequence mutations. Here, we applied a multiplexed CRISPR/Cas9 approach, informed by bioinformatic datasets, to generate genotypic variation in the promoter of OsSTOMAGEN, a positive regulator of rice stomatal density. Engineered genotypic variation corresponded to broad and continuous variation in stomatal density, ranging from 70% to 120% of wild-type stomatal density. This panel of stomatal variants was leveraged in physiological assays to establish discrete relationships between stomatal morphological variation and stomatal conductance, carbon assimilation and intrinsic water use efficiency in steady-state and fluctuating light conditions. Additionally, promoter alleles were subjected to vegetative drought regimes to assay the effects of the edited alleles on developmental response to drought. Notably, the capacity for drought-responsive stomatal density reprogramming in stomagen and two cis-regulatory edited alleles was reduced. Collectively our data demonstrate that cis-regulatory element editing can generate near-isogenic trait variation that can be leveraged for establishing relationships between anatomy and physiology, providing a basis for optimizing traits across diverse environments.

NAP1L1 Promotes Endometrial Cancer Progression via EP300-Mediated DDX5 Promoter Acetylation

Endometrial cancer is one of the predominant tumors of the female reproductive system. In this current study, we investigated the functions and related mechanisms of nucleosome assembly protein 1 like 1 (NAP1L1)/ DEAD-box helicase 5 (DDX5) in endometrial cancer. This retrospective study analyzed the medical records of patients with endometrial cancer, collected tissue samples for NAP1L1 and DDX5 staining, and conducted survival analysis using the Kaplan-Meier method. To evaluate the impact of NAP1L1 and/or DDX5 on cellular processes in endometrial cancer cells, several techniques were employed. These included Cell Counting Kit-8 assay, wound healing assay, Transwell assay, as well as overexpression or knockdown of target gene expression. Additionally, chromatin immunoprecipitation, dual luciferase reporter gene, and coimmunoprecipitation (Co-IP) assay were utilized to confirm the interaction between NAP1L1, E1A-binding protein p300 (EP300), and DDX5. Furthermore, qRT-PCR, Western blot, and Co-IP assay were performed to analyze the modulation of NAP1L1/DDX5 in Wnt/β-catenin. NAP1L1 and DDX5 expression were upregulated in endometrial cancer tissues, and correlated with poor prognosis. NAP1L1/DDX5 promoted endometrial cancer cell proliferation, migration, and invasion. NAP1L1 promotes acetylation and transcription by recruiting EP300 to the DDX5 promoter. DDX5 could activate Wnt/β-catenin signal by binding to β-catenin. In animal models, knockdown of NAP1L1 inhibits endometrial cancer tumor growth and lung metastasis. To sum up, our study demonstrated that NAP1L1 promoted the malignant phenotypes of endometrial cancer cells via recruiting EP300 to promote DDX5 acetylation, thus activating the Wnt/β-catenin signaling pathway. Implications: Our research findings indicate that targeting the NAP1L1/EP300/DX5 axis might be a new potential treatment option for endometrial cancer.

The F-Box Protein TaFBA1 Positively Regulates Drought Resistance and Yield Traits in Wheat

Environmental stresses, including drought stress, seriously threaten food security. Previous studies reported that wheat F-box protein, TaFBA1, responds to abiotic stresses in tobacco. Here, we generated transgenic wheat with enhanced (overexpression, OE) or suppressed (RNA interference, RNAi) expression of TaFBA1. The TaFBA1-OE seedlings showed enhanced drought tolerance, as measured by survival rate and fresh weight under severe drought stress, whereas the RNAi plants showed the opposite phenotype. Furthermore, the OE plants had stronger antioxidant capacity compared to WT and RNAi plants and maintained stomatal opening, which resulted in higher water loss under drought stress. However, stronger water absorption capacity in OE roots contributed to higher relative water contents in leaves under drought stress. Moreover, the postponed stomatal closure in OE lines helped to maintain photosynthesis machinery to produce more photoassimilate and ultimately larger seed size. Transcriptomic analyses conducted on WT and OE plants showed that genes involved in antioxidant, fatty acid and lipid metabolism and cellulose synthesis were significantly induced by drought stress in the leaves of OE lines. Together, our studies determined that the F-box protein TaFBA1 modulated drought tolerance and affected yield in wheat and the TaFBA1 gene could provide a desirable target for further breeding of wheat.

CmoPIP1-4 confers drought tolerance in pumpkin by altering hydrogen sulfide signaling

Drought is a major limiting factor for the growth and development of pumpkins. Plasma membrane intrinsic proteins (PIPs) are major water channels that play a crucial role in the regulation of cellular water status and solute trafficking during drought conditions. CmoPIP1-4 is a plasma membrane-localized protein that is significantly upregulated in roots and leaves under drought-stress conditions. In this study, the overexpression of CmoPIP1-4 enhances drought resistance in yeast. In contrast, CRISPR-mediated CmoPIP1-4 knockout in pumpkin roots increased drought sensitivity. This increased drought sensitivity of CmoPIP1-4 knockout plants is associated with a decline in the levels of hydrogen sulfide (H2S) and abscisic acid (ABA), accompanied by an increase in water loss caused by greater levels of transpiration and stomatal conductance. In addition, the sensitivity of CmoPIP1-4 CRISPR plants is further aggravated by reduced antioxidative enzyme activity, decreased proline and sugar contents, and extensive root damage. Furthermore, expression profiles of genes such as CmoHSP70s, CmoNCED3, CmoNCED4, and others involved in metabolic activities were markedly reduced in CmoPIP1-4 CRISPR plants. Moreover, we also discovered an interaction between the drought-responsive gene CmoDCD and CmoPIP1-4, indicating their potential role in activating H2S-mediated signaling in pumpkin, which could confer drought tolerance. The findings of our study collectively demonstrate CmoPIP1-4 plays a crucial role in the regulation of H2S-mediated signaling, influencing stomatal density and aperture in pumpkin plants, and thereby enhancing their drought tolerance.

Molecular Mechanisms and Regulatory Pathways Underlying Drought Stress Response in Rice

Rice is a staple food for 350 million people globally. Its yield thus affects global food security. Drought is a serious environmental factor affecting rice growth. Alleviating the inhibition of drought stress is thus an urgent challenge that should be solved to enhance rice growth and yield. This review details the effects of drought on rice morphology, physiology, biochemistry, and the genes associated with drought stress response, their biological functions, and molecular regulatory pathways. The review further highlights the main future research directions to collectively provide theoretical support and reference for improving drought stress adaptation mechanisms and breeding new drought-resistant rice varieties.

Paralog editing tunes rice stomatal density to maintain photosynthesis and improve drought tolerance

Rice (Oryza sativa) is of paramount importance for global nutrition, supplying at least 20% of global calories. However, water scarcity and increased drought severity are anticipated to reduce rice yields globally. We explored stomatal developmental genetics as a mechanism for improving drought resilience in rice while maintaining yield under climate stress. CRISPR/Cas9-mediated knockouts of the positive regulator of stomatal development STOMAGEN and its paralog EPIDERMAL PATTERNING FACTOR-LIKE10 (EPFL10) yielded lines with ∼25% and 80% of wild-type stomatal density, respectively. epfl10 lines with moderate reductions in stomatal density were able to conserve water to similar extents as stomagen lines but did not suffer from the concomitant reductions in stomatal conductance, carbon assimilation, or thermoregulation observed in stomagen knockouts. Moderate reductions in stomatal density achieved by editing EPFL10 present a climate-adaptive approach for safeguarding yield in rice. Editing the paralog of STOMAGEN in other species may provide a means for tuning stomatal density in agriculturally important crops beyond rice.

Phosphorylation of OsRbohB by the protein kinase OsDMI3 promotes H2O2 production to potentiate ABA responses in rice

In rice, the Ca²⁺/calmodulin-dependent protein kinase OsDMI3 is an important positive regulator of abscisic acid (ABA) signaling. In ABA signaling, H₂O₂ is required for ABA-induced activation of OsDMI3, which in turn increase H₂O₂ production. However, how OsDMI3 regulates H₂O₂ production in ABA signaling remains unknown. Here we show that OsRbohB is the main NADPH oxidase involved in ABA-induced H₂O₂ production and ABA-mediated physiological responses. OsDMI3 directly interacts with and phosphorylates OsRbohB at Ser-191, which is OsDMI3-mediated site-specific phosphorylation in ABA signaling. Further analyses revealed that OsDMI3-mediated OsRbohB Ser-191 phosphorylation positively regulates the activity of NADPH oxidase and the production of H₂O₂ in ABA signaling, thereby enhancing the sensitivity of seed germination and root growth to ABA and plant tolerance to water stress and oxidative stress. Moreover, we discovered that the OsDMI3-mediated OsRbohB phosphorylation and H₂O₂ production is dependent on the sucrose non-fermenting 1-related protein kinases SAPK8/9/10, which phosphorylate OsRbohB at Ser-140 in ABA signaling. Taken together, these results not only reveal an important regulatory mechanism that directly activates Rboh for ABA-induced H₂O₂ production but also uncover the importance of this regulatory mechanism in ABA signaling.

PtrVCS2 Regulates Drought Resistance by Changing Vessel Morphology and Stomatal Closure in Populus trichocarpa

Drought has severe effects on plant growth, forest productivity, and survival throughout the world. Understanding the molecular regulation of drought resistance in forest trees can enable effective strategic engineering of novel drought-resistant genotypes of tree species. In this study, we identified a gene, PtrVCS2, encoding a zinc finger (ZF) protein of the ZF-homeodomain transcription factor in Populus trichocarpa (Black Cottonwood) Torr. & A. Gray. ex Hook. Overexpression of PtrVCS2 (OE-PtrVCS2) in P. trichocarpa resulted in reduced growth, a higher proportion of smaller stem vessels, and strong drought-resistance phenotypes. Stomatal movement experiments revealed that the OE-PtrVCS2 transgenics showed lower stomata apertures than wild-type plants under drought conditions. RNA-seq analysis of the OE-PtrVCS2 transgenics showed that PtrVCS2 regulates the expression of multiple genes involved in regulation of stomatal opening and closing, particularly the PtrSULTR3;1-1 gene, and several genes related to cell wall biosynthesis, such as PtrFLA11-12 and PtrPR3-3. Moreover, we found that the water use efficiency of the OE-PtrVCS2 transgenic plants was consistently higher than that of wild type plants when subjected to chronic drought stress. Taken together, our results suggest that PtrVCS2 plays a positive role in improving drought adaptability and resistance in P. trichocarpa.

Saussurea involucrata PIP2;4 improves growth and drought tolerance in Nicotiana tabacum by increasing stomatal density and sensitivity

Aquaporins, the major facilitators of water transport across membranes, are involved in growth and development and adaptation to drought stress in plants. In this study, a plasma membrane intrinsic protein (SiPIP2;4) was cloned from Saussurea involucrata, a cold-tolerant hardy herb. The expression of SiPIP2;4 increased the stomatal density and sensitivity of tobacco (Nicotiana tabacum), thus, affecting the plant's growth and resistance to the diverse water environment. The higher stomatal density under well-watered conditions effectively promoted the photosynthetic rate, which led to the rapid growth of transgenic lines. The stomata in the transgenic lines responded more sensitively to the vapor pressure deficit than the wild-type under different levels of ambient humidity. Their stomatal apertures positively correlated with the ambient humidity. Under drought conditions, the overexpression of SiPIP2;4 promoted rapid stomatal closure, reduced water dissipation, and enhanced drought tolerance. These results indicate that SiPIP2;4 regulates the density and sensitivity of plant stomata, thus, playing an important role in balancing plant growth and stress tolerance. This suggests that SiPIP2;4 has the potential to serve as a genetic resource for crop improvement.

Multi-model genome-wide association studies of leaf anatomical traits and vein architecture in rice

Introduction

The anatomy of rice leaves is closely related to photosynthesis and grain yield. Therefore, exploring insight into the quantitative trait loci (QTLs) and alleles related to rice flag leaf anatomical and vein traits is vital for rice improvement.

Methods

Here, we aimed to explore the genetic architecture of eight flag leaf traits using one single-locus model; mixed-linear model (MLM), and two multi-locus models; fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage disequilibrium iteratively nested keyway (BLINK). We performed multi-model GWAS using 329 rice accessions of RDP1 with 700K single-nucleotide polymorphisms (SNPs) markers.

Results

The phenotypic correlation results indicated that rice flag leaf thickness was strongly correlated with leaf mesophyll cells layer (ML) and thickness of both major and minor veins. All three models were able to identify several significant loci associated with the traits. MLM identified three non-synonymous SNPs near NARROW LEAF 1 (NAL1) in association with ML and the distance between minor veins (IVD) traits.

Discussion

Several numbers of significant SNPs associated with known gene function in leaf development and yield traits were detected by multi-model GWAS performed in this study. Our findings indicate that flag leaf traits could be improved via molecular breeding and can be one of the targets in high-yield rice development.

Tomato ARPC1 regulates trichome morphology and density and terpene biosynthesis

Based on transcriptomic analysis of wild-type and mutant tomato plants, ARPC1 was found to be important for trichome formation and development and it plays a key role in terpene synthesis. Trichomes are protruding epidermal cells in plant species. They function as the first defense layer against biotic and abiotic stresses. Despite the essential role of tomato trichomes in defense against herbivores, the understanding of their development is still incomplete. Therefore, the aim of this study was to identify genes involved in trichome formation and morphology and terpene synthesis, using transcriptomic techniques. To achieve this, we examined leaf morphology and compared the expression levels of some putative genes involved in trichome formation between wild-type (WT) and hairless-3 (hl-3) tomato mutant. The hl-3 plants displayed swollen and distorted trichomes and reduced trichome density (type I and IV) and terpene synthesis compared with that of the WT plants. Gene expression analysis showed that Actin-Related Protein Component1 (ARPC1) was expressed more highly in the WT than in the hl-3 mutant, indicating its critical role in trichome morphology and density. Additionally, the expression of MYC1 and several terpene synthase genes (TPS9, 12, 20), which are involved in type VI trichome initiation and terpene synthesis, was lower in the hl-3 mutant than in the WT plants. Moreover, transformation of the hl-3 mutant with WT ARPC1 restored normal trichome structure and density, and terpene synthesis. Structural and amino acid sequence analysis showed that there was a missplicing mutation in the hl-3 mutant, which was responsible for the abnormal trichome structure and density, and impaired terpene synthesis. Overall, the findings of this study demonstrated that ARPC1 is involved in regulating trichome structure and terpene synthesis in tomato.

Identification of the Regulators of Epidermis Development under Drought- and Salt-Stressed Conditions by Single-Cell RNA-Seq

As sessile organisms, plants constantly face challenges from the external environment. In order to meet these challenges and survive, plants have evolved a set of sophisticated adaptation strategies, including changes in leaf morphology and epidermal cell development. These developmental patterns are regulated by both light and hormonal signaling pathways. However, our mechanistic understanding of the role of these signaling pathways in regulating plant response to environmental stress is still very limited. By applying single-cell RNA-Seq, we determined the expression pattern of PHYTOCHROME INTERACTING FACTOR (PIF) 1, PIF3, PIF4, and PIF5 genes in leaf epidermal pavement cells (PCs) and guard cells (GCs). PCs and GCs are very sensitive to environmental stress, and our previous research suggests that these PIFs may be involved in regulating the development of PCs, GCs, and leaf morphology under environmental stress. Growth analysis showed that pif1/3/4/5 quadruple mutant maintained tolerance to drought and salt stress, and the length to width ratio of leaves and petiole length under normal growth conditions were similar to those of wild-type (WT) plants under drought and salt treatment. Analysis of the developmental patterns of PCs and GCs, and whole leaf morphology, further confirmed that these PIFs may be involved in mediating the development of epidermal cells under drought and salt stress, likely by regulating the expression of MUTE and TOO MANY MOUTHS (TMM) genes. These results provide new insights into the molecular mechanism of plant adaptation to adverse growth environments.

CmNF-YB8 affects drought resistance in chrysanthemum by altering stomatal status and leaf cuticle thickness

Drought is a major abiotic stress that limits plant growth and development. Adaptive mechanisms have evolved to mitigate drought stress, including the capacity to adjust water loss rate and to modify the morphology and structure of the epidermis. Here, we show that the expression of CmNF-YB8, encoding a nuclear factor Y (NF-Y) B-type subunit, is lower under drought conditions in chrysanthemum (Chrysanthemum morifolium). Transgenic chrysanthemum lines in which transcript levels of CmNF-YB8 were reduced by RNA interference (CmNF-YB8-RNAi) exhibited enhanced drought resistance relative to control lines, whereas lines overexpressing CmNF-YB8 (CmNF-YB8-OX) were less tolerant to drought. Compared to wild type (WT), CmNF-YB8-RNAi plants showed reduced stomatal opening and a thicker epidermal cuticle that correlated with their water loss rate. We also identified genes involved in stomatal adjustment (CBL-interacting protein kinase 6, CmCIPK6) and cuticle biosynthesis (CmSHN3) that are more highly expressed in CmNF-YB8-RNAi lines than in WT, CmCIPK6 being a direct downstream target of CmNF-YB8. Virus-induced gene silencing of CmCIPK6 or CmSHN3 in the CmNF-YB8-RNAi background abolished the effects of CmNF-YB8-RNAi on stomatal closure and cuticle deposition, respectively. CmNF-YB8 thus regulates CmCIPK6 and CmSHN3 expression to alter stomatal movement and cuticle thickness in the leaf epidermis, thereby affecting drought resistance.

GhVLN4 is involved in multiple stress responses and required for resistance to Verticillium wilt

As structural and signaling platform in plant cell, the actin cytoskeleton is regulated by diverse actin binding proteins (ABPs). Villins are one type of major ABPs responsible for microfilament bundling, which have proved to play important roles in plant growth and development. However, the function of villins in stress tolerance is poorly understood. Here, we report the function of cotton GhVLN4 in Verticillium wilt resistance and abiotic stress tolerance. The expression of GhVLN4 was up-regulated by gibberellin, ethylene, ABA, salicylic acid, jasmonate, NaCl, PEG, and Verticillium dahliae treatment, suggesting the involvement of GhVLN4 in multiple stress and hormone responses and signaling. Virus-induced gene silencing GhVLN4 made cotton more susceptible to V. dahliae characterized by the preferential colonization and rapid growth of the fungus in both phloem and xylem of the infected stems. Arabidopsis overexpressing GhVLN4 exhibited higher resistance to V. dahliae, salt and drought than the wild-type plants. The enhanced resistance to V. dahliae is likely related to the upregulated components in SA signaling pathway; the improved tolerance to salt and drought is characterized by upregulation of the components both in ABA- related and ABA-independent signal pathways, along with altered stomatal aperture under drought. Our findings demonstrate that GhVLN4 may play important roles in regulating plant tolerance to both biotic and abiotic stresses.

Weighted gene coexpression network analysis-based identification of key modules and hub genes associated with drought sensitivity in rice

BACKGROUND

Drought stress is an adverse factor with deleterious effects on several aspects of rice growth. However, the mechanism underlying drought resistance in rice remains unclear. To understand the molecular mechanism of the drought response in rice, drought-sensitive CSSL (Chromosome Single-substitution Segment Line) PY6 was used to map QTLs of sensitive phenotypes and to reveal the impact of the QTLs on transcriptional profiling.

RESULTS

The QTL dss-1 was mapped onto the short arm of chromosome 1 of rice. According to transcriptomic analysis, the identified differentially expressed genes (DEGs) exhibited a downregulated pattern and were mainly enriched in photosynthesis-related GO terms, indicating that photosynthesis was greatly inhibited under drought. Further, according to weighted gene coexpression network analysis (WGCNA), specific gene modules (designating a group of genes with a similar expression pattern) were strongly correlated with H₂O₂ (4 modules) and MDA (3 modules), respectively. Likewise, GO analysis revealed that the photosynthesis-related GO terms were consistently overrepresented in H₂O₂-correlated modules. Functional annotation of the differentially expressed hub genes (DEHGs) in the H₂O₂ and MDA-correlated modules revealed cross-talk between abiotic and biotic stress responses for these genes, which were annotated as encoding WRKYs and PR family proteins, were notably differentially expressed between PY6 and PR403.

CONCLUSIONS

We speculated that drought-induced photosynthetic inhibition leads to H₂O₂ and MDA accumulation, which can then trigger the reprogramming of the rice transcriptome, including the hub genes involved in ROS scavenging, to prevent oxidative stress damage. Our results shed light on and provide deep insight into the drought resistance mechanism in rice.

GmNAP1 is essential for trichome and leaf epidermal cell development in soybean

KEY MESSAGE

Map-based cloning revealed that two novel soybean distorted trichome mutants were due to loss function of GmNAP1 gene, which affected the trichome morphology and pavement cell ploidy by regulating actin filament assembly. Trichomes increase both biotic and abiotic stress resistance in soybean. In this study, Gmdtm1-1 and Gmdtm1-2 mutants with shorter trichomes and bigger epidermal pavement cells were isolated from an ethyl methylsulfonate mutagenized population. Both of them had reduced plant height and smaller seeds. Map-based cloning and bulked segregant analysis identified that a G-A transition at the 3' boundary of the sixth intron of Glyma.20G019300 in the Gmdtm1-1 mutant and another G-A transition mutation at the 5' boundary of the fourteenth intron of Glyma.20G019300 in Gmdtm1-2; these mutations disrupted spliceosome recognition sites creating truncated proteins. Glyma.20G019300 encodes a Glycine max NCK-associated protein 1 homolog (GmNAP1) in soybean. Further analysis revealed that the GmNAP1 involved in actin filament assembling and genetic information processing pathways during trichome and pavement cell development. This study shows that GmNAP1 plays an important role in soybean growth and development and agronomic traits.

Genome-wide identification and functional characterization of cotton (Gossypium hirsutum) MAPKKK gene family in response to drought stress

BACKGROUND

Mitogen-activated protein kinase kinase kinases (MAPKKKs) are significant components in the MAPK signal pathway and play essential roles in regulating plants against drought stress. To explore MAPKKK gene family functioning in cotton response and resistance to drought stress, we conducted a systematic analysis of GhMAPKKKs.

RESULTS

In this study, 157 nonredundant GhMAPKKKs (including 87 RAFs, 46 MEKKs and 24 ZIKs) were identified in cotton (Gossypium hirsutum). These GhMAPKKK genes are unevenly distributed on 26 chromosomes, and segmental duplication is the major way for the enlargement of MAPKKK family. Furthermore, members within the same subfamily share a similar gene structure and motif composition. A lot of cis-elements relevant to plant growth and response to stresses are distributed in promoter regions of GhMAPKKKs. Additionally, these GhMAPKKKs show differential expression patterns in cotton tissues. The transcription levels of most genes were markedly altered in cotton under heat, cold and PEG treatments, while the expressions of some GhMAPKKKs were induced in cotton under drought stress. Among these drought-induced genes, we selected GhRAF4 and GhMEKK12 for further functional characterization by virus-induced gene silencing (VIGS) method. The experimental results indicated that the gene-silenced cotton displayed decreased tolerance to drought stress. Malondialdehyde (MDA) content was higher, but proline accumulation, relative leaf water content and activities of superoxide dismutase (SOD) and peroxidase (POD) were lower in the gene-silenced cotton, compared with those in the controls, under drought stress.

CONCLUSION

Collectively, a systematic survey of gene structure, chromosomal location, motif composition and evolutionary relationship of MAPKKKs were performed in upland cotton (Gossypium hirsutum). The following expression and functional study showed that some of them take important parts in cotton drought tolerance. Thus, the data presented here may provide a foundation for further investigating the roles of GhMAPKKKs in cotton response and resistance to drought stress.