Genome-wide identification and mining elite allele variation of the Monoacylglycerol lipase (MAGL) gene family in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Monoacylglycerol lipase (MAGL) genes belong to the alpha/beta hydrolase superfamily, catalyze the terminal step of triglyceride (TAG) hydrolysis, converting monoacylglycerol (MAG) into free fatty acids and glycerol. RESULTS: In this study, 30 MAGL genes in upland cotton have been identified, which have been classified into eight subgroups. The duplication of GhMAGL genes in upland cotton was predominantly influenced by segmental duplication events, as revealed through synteny analysis. Furthermore, all GhMAGL genes were found to contain light-responsive elements. Through comprehensive association and haplotype analyses using resequencing data from 355 cotton accessions, GhMAGL3 and GhMAGL6 were detected as key genes related to lipid hydrolysis processes, suggesting a negative regulatory effect. CONCLUSIONS: In summary, MAGL has never been studied in upland cotton previously. This study provides the genetic mechanism foundation for the discover of new genes involved in lipid metabolism to improve cottonseed oil content, which will provide a strategic avenue for marker-assisted breeding aimed at incorporating desirable traits into cultivated cotton varieties.

Insights into the Genomic Regions and Candidate Genes of Senescence-Related Traits in Upland Cotton via GWAS.

Senescence is the last stage of plant development and is controlled by both internal and external factors. Premature senescence significantly affects the yield and quality of cotton. However, the genetic architecture underlying cotton senescence remains unclear. In this study, genome-wide association studies (GWAS) were performed based on 3,015,002 high-quality SNP markers from the resequencing data of 355 upland cotton accessions to detect genomic regions for cotton senescence. A total of 977 candidate genes within 55 senescence-related genomic regions (SGRs), SGR1-SGR55, were predicted. Gene ontology (GO) analysis of candidate genes revealed that a set of biological processes was enriched, such as salt stress, ethylene processes, and leaf senescence. Furthermore, in the leaf senescence GO term, one candidate gene was focused on: Gohir.A12G270900 (GhMKK9), located in SGR36, which encodes a protein of the MAP kinase kinase family. Quantitative real-time PCR (qRT-PCR) analysis showed that GhMKK9 was up-regulated in old cotton leaves. Overexpression of GhMKK9 in Arabidopsis accelerated natural leaf senescence. Virus-induced gene silencing (VIGS) of GhMKK9 in cotton increased drought tolerance. These results suggest that GhMKK9 is a positive regulator and might be involved in drought-induced senescence in cotton. The results provide new insights into the genetic basis of cotton senescence and will be useful for improving cotton breeding in the future.

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.).

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

Identification and profiling of microRNAs and differentially expressed genes during anther development between a genetic male-sterile mutant and its wildtype cotton via high-throughput RNA sequencing.

Genetic male sterility (GMS) facilitates hybrid seed production in crops including cotton (Gossypium hirsutum). However, the genetic and molecular mechanisms specifically involved in this developmental process are poorly understood. In this study, small RNA sequencing, degradome sequencing, and transcriptome sequencing were performed to analyze miRNAs and their target genes during anther development in a GMS mutant ('Dong A') and its fertile wildtype (WT). A total of 80 known and 220 novel miRNAs were identified, 71 of which showed differential expressions during anther development. A further degradome sequencing revealed a total of 117 candidate target genes cleaved by 16 known and 36 novel miRNAs. Based on RNA-seq, 24, 11, and 21 predicted target genes showed expression correlations with the corresponding miRNAs at the meiosis, tetrad and uninucleate stages, respectively. In addition, a large number of differentially expressed genes were identified, most of which were involved in sucrose and starch metabolism, carbohydrate metabolism, and plant hormone signal transduction based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The results of our study provide valuable information for further functional investigations of the important miRNAs and target genes involved in genetic male sterility and advance our understanding of miRNA regulatory functions during cotton anther development.

The WRKY transcription factor GhWRKY27 coordinates the senescence regulatory pathway in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Premature senescence can reduce the yield and quality of crops. WRKY transcription factors (TFs) play important roles during leaf senescence, but little is known about their ageing mechanisms in cotton. RESULTS: In this study, a group III WRKY TF, GhWRKY27, was isolated and characterized. The expression of GhWRKY27 was induced by leaf senescence and was higher in an early-ageing cotton variety than in a non-early-ageing cotton variety. Overexpression of GhWRKY27 in Arabidopsis promoted leaf senescence, as determined by reduced chlorophyll content and elevated expression of senescence-associated genes (SAGs). Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that GhWRKY27 interacted with an MYB TF, GhTT2. Putative target genes of GhWRKY27 were identified via chromatin immunoprecipitation followed by sequencing (ChIP-seq). Yeast one-hybrid (Y1H) assay and electrophoretic mobility shift assay (EMSA) revealed that GhWRKY27 binds directly to the promoters of cytochrome P450 94C1 (GhCYP94C1) and ripening-related protein 2 (GhRipen2-2). In addition, the expression patterns of GhTT2, GhCYP94C1 and GhRipen2-2 were identified during leaf senescence. Transient dual-luciferase reporter assay indicated that GhWRKY27 could activate the expression of GhCYP94C1 and GhRipen2-2. CONCLUSIONS: Our work lays the foundation for further study of the functional roles of WRKY genes during leaf senescence in cotton. In addition, our data provide new insights into the senescence-associated mechanisms of WRKY genes in cotton.

Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton.

KEY MESSAGE: Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton. A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

Identification of favorable SNP alleles and candidate genes for traits related to early maturity via GWAS in upland cotton.

BACKGROUND: Early maturity is one of the most important and complex agronomic traits in upland cotton (Gossypium hirsutum L). To dissect the genetic architecture of this agronomically important trait, a population consisting of 355 upland cotton germplasm accessions was genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) approach, of which a subset of 185 lines representative of the diversity among the accessions was phenotypically characterized for six early maturity traits in four environments. A genome-wide association study (GWAS) was conducted using the generalized linear model (GLM) and mixed linear model (MLM). RESULTS: A total of 81,675 SNPs in 355 upland cotton accessions were discovered using SLAF-seq and were subsequently used in GWAS. Thirteen significant associations between eight SNP loci and five early maturity traits were successfully identified using the GLM and MLM; two of the 13 associations were common between the models. By computing phenotypic effect values for the associations detected at each locus, 11 highly favorable SNP alleles were identified for five early maturity traits. Moreover, dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between the numbers of highly favorable alleles and the phenotypic values of the target traits were identified. Most importantly, a major locus (rs13562854) on chromosome Dt3 and a potential candidate gene (CotAD_01947) for early maturity were detected. CONCLUSIONS: This study identified highly favorable SNP alleles and candidate genes associated with early maturity traits in upland cotton. The results demonstrate that GWAS is a powerful tool for dissecting complex traits and identifying candidate genes. The highly favorable SNP alleles and candidate genes for early maturity traits identified in this study should be show high potential for improvement of early maturity in future cotton breeding programs.

Genome-wide analysis of the family 1 glycosyltransferases in cotton.

Family 1 GT, designated as UGT, is the largest and most functionally important multigene family in the plant kingdom. In this study, we carried out a genome-wide identification, analysis, and comparison of 142, 146, and 196 putative UGTs from Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum, respectively. All members present the 44 amino-acid conserved consensus sequence termed the plant secondary product glycosyltransferase motif. According to the phylogenetic relationship among the cotton UGT proteins and those from other species, GrUGTs and GaUGTs could be classified into 16 major phylogenetic groups (A-P), whereas GhUGTs are classified into 15 major phylogenetic groups with a lack of group C. All cotton UGTs are dispersed throughout the chromosomes and are displayed in clusters with the same open reading frame orientation. The expansion of them appears to result from genome duplication and rearrangement. Two conserved introns, A and B, are detected in most of the intron-containing-UGTs in G. raimondii and G. arboreum, whereas only intron A is detected in the intron-containing-UGTs in G. hirsutum. Furthermore, expression patterns of the UGT genes in G. hirsutum wild type and its near isogenic fuzzless-lintless mutant at the stage of fiber initiation were analyzed using the RNA-seq data. Overall, this study not only deepens our understanding of the structure, phylogeny, evolution, and expression of cotton UGT genes, but also provides a solid foundation for further cloning and functional studies of the UGT family genes.

Genome-Wide Analysis of the SAUR Gene Family and Its Expression Profiles in Response to Salt Stress in Santalum album.

The SAUR (small auxin-up RNA) family constitutes a category of genes that promptly respond to the hormone auxin and play a pivotal role in diverse biological processes encompassing plant growth and the response to abiotic stress. Santalum album L., a semi-parasitic evergreen tree, is renowned for its economically valuable essential oils, positioning it among the most prized tree species. In this study, a meticulous identification and comprehensive analysis of 43 SAUR genes was conducted within S. album. Based on phylogenetic relationships, the SaSAUR genes were systematically categorized into five groups. A collinearity analysis revealed intriguing insights, disclosing 14 segmental duplications and 9 tandem duplications within the SaSAUR genes, emphasizing the pivotal role of duplication in the expansion of this gene family. Noteworthy variations in the expression levels of SaSAUR genes were observed by delving into the SaSAUR transcriptome data from various tissues, including leaves, roots, and heartwood, as well as under salt-stress conditions. Notably, SaSAUR08 and SaSAUR13 were significantly upregulated in heartwood compared with roots and leaves, while SaSAUR18 was markedly more expressed in roots compared with heartwood and leaves. Furthermore, SaSAUR27 and SaSAUR28 were found to respond closely to salt stress, hinting at their potential involvement in the salt-stress response mechanism. This research offers a comprehensive investigation of SAUR genes in S. album and establishes a foundation for future exploration of the SAUR gene family, particularly its relation to growth and salt-stress responses.

Genome-wide identification of R-SNARE gene family in upland cotton and function analysis of GhVAMP72l response to drought stress.

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (R-SNAREs) mainly promoted the assembly of the SNARE complex to drive the final membrane fusion step of membrane transport. Previous research on R-SNAREs has mainly focused on development and growth and has rarely been involved in abiotic stress, especially in cotton. Here, we performed a comprehensive analysis of R-SNARE genes in upland cotton. In total, 51 Gh-R-SNARE genes across six phylogenetic groups were unevenly distributed on 21 chromosomes. Cis elements related to plant growth and response to abiotic stress responses were found in the promoter region of Gh-R-SNAREs. Nine Gh-R-SNARE genes were obviously upregulated under drought stress conditions by RNA-seq and qRT-PCR analysis. Among them, GhVAMP72l might be the key candidate gene contributing to drought stress tolerance in cotton by virus-induced gene silencing (VIGS) assay. These results provide valuable insights for the functional analysis of cotton R-SNAREs in response to drought stress and highlight potential beneficial genes for genetic improvement and breeding in cotton.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS.

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

OsSHMT4 Is Required for Synthesis of Rice Storage Protein and Storage Organelle Formation in Endosperm Cells.

Storage proteins are essential for seed germination and seedling growth, as they provide an indispensable nitrogen source and energy. Our previous report highlighted the defective endosperm development in the serine hydroxymethyltransferase 4 (OsSHMT4) gene mutant, floury endosperm20-1 (flo20-1). However, the alterations in storage protein content and distribution within the flo20-1 endosperm remained unclear. Here, the immunocytochemistry analyses revealed a deficiency in storage protein accumulation in flo20-1. Electron microscopic observation uncovered abnormal morphological structures in protein bodies (PBI and PBII) in flo20-1. Immunofluorescence labeling demonstrated that aberrant prolamin composition could lead to the subsequent formation and deposition of atypical structures in protein body I (PBI), and decreased levels of glutelins and globulin resulted in protein body II (PBII) malformation. Further RNA-seq data combined with qRT-PCR results indicated that altered transcription levels of storage protein structural genes were responsible for the abnormal synthesis and accumulation of storage protein, which further led to non-concentric ring structural PBIs and amorphous PBIIs. Collectively, our findings further underscored that OsSHMT4 is required for the synthesis and accumulation of storage proteins and storage organelle formation in endosperm cells.

Characterization and Functional Analysis of GhNAC82, A NAM Domain Gene, Coordinates the Leaf Senescence in Upland Cotton (Gossypium hirsutum L.).

In the process of growth and development, cotton exhibits premature senescence under various abiotic stresses, impairing yield and fiber quality. NAC (NAM, ATAF1,2, and CUC2) protein widely distributed in land plants, play the critical role in responding to abiotic stress and regulating leaf senescence. We have identified and functional analyzed a NAM domain gene GhNAC82 in upland cotton, it was located on the A11 chromosome 4,921,702 to 4,922,748 bp, only containing one exon. The spatio-temporal expression pattern analysis revealed that it was highly expressed in root, torus, ovule and fiber development stage. The results of qRT-PCR validated that GhNAC82 negatively regulated by salt stress, drought stress, H2O2 stress, IAA treatment, and ethylene treatment, positively regulated by the ABA and MeJA treatment. Moreover, heterologous overexpression of GhNAC82 results in leaf premature senescence and delays root system development in Arabidopsis thaliana. The phenotype of delayed-senescence was performed after silencing GhNAC82 by VIGS in premature cotton. Taken together, GhNAC82 was involved in different abiotic stress pathways and play important roles in negatively regulating leaf premature senescence.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton.

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

GhCDPK60 positively regulates drought stress tolerance in both transgenic Arabidopsis and cotton by regulating proline content and ROS level.

Calcium-Dependent Protein Kinases (CDPKs) involved in regulating downstream components of calcium signaling pathways play a role in tolerance to abiotic stresses and seed development in plants. However, functions of only a few cotton CDPKs have been clarified at present. In this study, 80 conserved CDPKs in Gossypium hirsutum L. were identified and characterized, which was divided into four subgroups. Among them, the transcript level of GhCDPK60 was significantly upregulated under drought and several hormone treatments. And we found that the expression levels of several stress-inducible genes down-regulated in GhCDPK60-silence cotton and up-regulated in GhCDPK60-overexpressing Arabidopsis. In addition, physiological analyses demonstrated that GhCDPK60 improved drought stress tolerance by improving the osmotic adjustment ability and reducing the accumulation of reactive oxygen species (ROS) in plants. These findings broaden our understanding of the biological roles of GhCDPK60 and mechanisms underlying drought stress tolerance in cotton.

Uncovering Novel Genomic Regions and Candidate Genes for Senescence-Related Traits by Genome-Wide Association Studies in Upland Cotton (Gossypium hirsutum L.).

Senescence in plants is a complex trait, which is controlled by both genetic and environmental factors and can affect the yield and quality of cotton. However, the genetic basis of cotton senescence remains relatively unknown. In this study, we reported genome-wide association studies (GWAS) based on 185 accessions of upland cotton and 26,999 high-quality single-nucleotide polymorphisms (SNPs) to reveal the genetic basis of cotton senescence. To determine cotton senescence, we evaluated eight traits/indices. Our results revealed a high positive correlation (r>0.5) among SPAD value 20 days after topping (SPAD20d), relative difference of SPAD (RSPAD), nodes above white flower on topping day (NAWF0d), nodes above white flower 7 days after topping (NAWF7d), and number of open bolls on the upper four branches (NB), and genetic analysis revealed that all traits had medium or high heritability ranging from 0.53 to 0.86. Based on a multi-locus method (FASTmrMLM), a total of 63 stable and significant quantitative trait nucleotides (QTNs) were detected, which represented 50 genomic regions (GWAS risk loci) associated with cotton senescence. We observed three reliable loci located on chromosomes A02 (A02_105891088_107196428), D03 (D03_37952328_38393621) and D13 (D13_59408561_60730103) because of their high repeatability. One candidate gene (Ghir_D03G011060) was found in the locus D03_37952328_38393621, and its Arabidopsis thaliana homologous gene (AT5G23040) encodes a cell growth defect factor-like protein (CDF1), which might be involved in chlorophyll synthesis and cell death. Moreover, qRT-PCR showed that the transcript level of Ghir_D03G011060 was down-regulated in old cotton leaves, and virus-induced gene silencing (VIGS) indicated that silencing of Ghir_D03G011060 resulted in leaf chlorosis and promoted leaf senescence. In addition, two candidate genes (Ghir_A02G017660 and Ghir_D13G021720) were identified in loci A02_105891088_107196428 and D13_59408561_60730103, respectively. These results provide new insights into the genetic basis of cotton senescence and will serve as an important reference for the development and implementation of strategies to prevent premature senescence in cotton breeding programs.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses.

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F 2 population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

Transcriptome analysis of nitric oxide-responsive genes in upland cotton (Gossypium hirsutum).

Nitric oxide (NO) is an important signaling molecule with diverse physiological functions in plants. It is therefore important to characterize the downstream genes and signal transduction networks modulated by NO. Here, we identified 1,932 differentially expressed genes (DEGs) responding to NO in upland cotton using high throughput tag sequencing. The results of quantitative real-time polymerase chain reaction (qRT-PCR) analysis of 25 DEGs showed good consistency. Gene Ontology (GO) and KEGG pathway were analyzed to gain a better understanding of these DEGs. We identified 157 DEGs belonging to 36 transcription factor (TF) families and 72 DEGs related to eight plant hormones, among which several TF families and hormones were involved in stress responses. Hydrogen peroxide and malondialdehyde (MDA) contents were increased, as well related genes after treatment with sodium nitroprusside (SNP) (an NO donor), suggesting a role for NO in the plant stress response. Finally, we compared of the current and previous data indicating a massive number of NO-responsive genes at the large-scale transcriptome level. This study evaluated the landscape of NO-responsive genes in cotton and identified the involvement of NO in the stress response. Some of the identified DEGs represent good candidates for further functional analysis in cotton.

Identification of the group IIa WRKY subfamily and the functional analysis of GhWRKY17 in upland cotton (Gossypium hirsutum L.).

WRKY transcription factors play important roles in plant defense, stress response, leaf senescence, and plant growth and development. Previous studies have revealed the important roles of the group IIa GhWRKY genes in cotton. To comprehensively analyze the group IIa GhWRKY genes in upland cotton, we identified 15 candidate group IIa GhWRKY genes in the Gossypium hirsutum genome. The phylogenetic tree, intron-exon structure, motif prediction and Ka/Ks analyses indicated that most group IIa GhWRKY genes shared high similarity and conservation and underwent purifying selection during evolution. In addition, we detected the expression patterns of several group IIa GhWRKY genes in individual tissues as well as during leaf senescence using public RNA sequencing data and real-time quantitative PCR. To better understand the functions of group IIa GhWRKYs in cotton, GhWRKY17 (KF669857) was isolated from upland cotton, and its sequence alignment, promoter cis-acting elements and subcellular localization were characterized. Moreover, the over-expression of GhWRKY17 in Arabidopsis up-regulated the senescence-associated genes AtWRKY53, AtSAG12 and AtSAG13, enhancing the plant's susceptibility to leaf senescence. These findings lay the foundation for further analysis and study of the functions of WRKY genes in cotton.