Membrane transporters encoded by NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER (NPF) genes, which play crucial roles in plant growth, development and resistance to various stresses, are involved in the transport of nitrate (NO₃^-) and peptides. In several plant species, NPF genes are involved in the resistance to abiotic stresses; however, whether the whole NPF gene family in cotton contributes to this resistance has not been systematically investigated. Here, 201 genes encoding NPF proteins with a peptide transporter (PTR) domain were confirmed in three different Gossypium species, namely, Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii. The NPF proteins in these three Gossypium species and Arabidopsis thaliana were classified into three different subfamilies via phylogenetic analysis. Among the genes that encode these proteins, most GhNPF genes in the same subfamily contained similar gene structures and conserved domains. Predictions of the promoters of these genes revealed that the cis-acting elements included phytohormone- and light-responsive elements, indicating that some of these genes might be expressed in response to abiotic stress. Furthermore, 52 common potential candidate genes in 98 GhNPFs were predicted to exhibit specific spatiotemporal expression patterns in different tissues based on two RNA sequencing (RNA-seq) datasets. Finally, the gene expression profiles of abiotic stress indicated that 31 GhNPF genes were upregulated in at least one treatment period. Under abiotic stress for 12 and 24 h, the expression of GhNPF8 was upregulated upon cold treatment but downregulated with heat treatment, salt treatment and drought treatment. Furthermore, the expression of genes GhNPF8, GhNPF54 and GhNPF43 peaked at 6 h after heat and salt treatment. These results indicated that these genes exhibit underlying characteristics related to responses to abiotic stress. The verification of NPFs and analysis of their expression profiles in different tissues and in response to different abiotic stresses of cotton provide a basis for further studying the relationship between abiotic stress resistance and nitrogen (N) transport in cotton, as well as identifying candidate genes to facilitate their functional identification.

The boll opening rate (BOR) is an early maturity trait that plays a crucial role in cotton production in China, as BOR has a significant effect on defoliant spraying and picking time of unginned cotton, ultimately determining yield and fiber quality. Therefore, elucidating the genetic basis of BOR and identifying stably associated loci, elite alleles, and potential candidate genes can effectively accelerate the molecular breeding process. In this study, we utilized the mixed linear model (MLM) algorithm to perform a genome-wide association study (GWAS) based on 4,452,629 single-nucleotide polymorphisms (SNPs) obtained through whole-genome resequencing of a natural population of 418 upland cotton accessions and phenotypic BOR data acquired from five environments. A total of 18 SNP loci were identified on chromosome D11 that are stable and significantly associated with BOR in multiple environments. Moreover, a significant SNP peak (23.703–23.826 Mb) was identified, and a GH-D11G2034 gene and favorable allelic variation (GG) related to BOR were found in this genomic region, significantly increasing cotton BOR. Evolutionary studies have shown that GH-D11G2034 may have been subjected to artificial selection throughout the variety selection process. This study provides valuable insights and suggests that the GH-D11G2034 gene and its favorable allelic variation (GG) could be potential targets for molecular breeding to improve BOR in upland cotton. However, further research is needed to validate the function of this gene and explore its potential applications in cotton breeding programs. Overall, this study contributes to the advancement of genetic improvement in early maturity and has important implications for the sustainable development of the cotton industry.

Annexins (ANNs) are a structurally conserved protein family present in almost all plants. In the present study, 27 GhANNs were identified in cotton and were unevenly distributed across 14 chromosomes. Transcriptome data and RT-qPCR results revealed that multiple GhANNs respond to at least two abiotic stresses. Similarly, the expression levels of GhANN4 and GhANN11 were significantly upregulated under heat, cold, and drought stress. Using virus-induced gene silencing (VIGS), functional characterization of GhANN4 and GhANN11 revealed that, compared with those of the controls, the leaf wilting of GhANN4-silenced plants was more obvious, and the activities of catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD) were lower under NaCl and PEG stress. Moreover, the expression of stress marker genes (GhCBL3, GhDREB2A, GhDREB2C, GhPP2C, GhRD20-2, GhCIPK6, GhNHX1, GhRD20-1, GhSOS1, GhSOS2 and GhSnRK2.6) was significantly downregulated in GhANN4-silenced plants after stress. Under cold stress, the growth of the GHANN11-silenced plants was significantly weaker than that of the control plants, and the activities of POD, SOD, and CAT were also lower. However, compared with those of the control, the elasticity and orthostatic activity of the GhANN11-silenced plants were greater; the POD, SOD, and CAT activities were higher; and the GhDREB2C, GhHSP, and GhSOS2 expression levels were greater under heat stress. These results suggest that different GhANN family members respond differently to different types of abiotic stress.