Cotton is an important economic crop affected by different abiotic stresses at different developmental stages. Salinity limits the growth and productivity of crops worldwide. Na⁺/H⁺ antiporters play a key role during the plant development and in its tolerance to salt stress. The aim of the present study was a genome-wide characterization and expression pattern analysis under the salinity stress of the sodium-proton antiporter (NHX) of Gossypium barbadense in comparison with Gossypium hirsutum. In G. barbadense, 25 NHX genes were identified on the basis of the Na⁺_H⁺ exchanger domain. All except one of the G. barbadense NHX transporters have an Amiloride motif that is a known inhibitor of Na⁺ ions in plants. A phylogenetic analysis inferred three classes of GbNHX genes—viz., Vac (GbNHX1, 2 and 4), Endo (GbNHX6), and PM (GbNHX7). A high number of the stress-related cis-acting elements observed in promoters show their role in tolerance against abiotic stresses. The Ka/Ks values show that the majority of GbNHX genes are subjected to strong purifying selection under the course of evolution. To study the functional divergence of G. barbadense NHX transporters, the real-time gene expression was analyzed under salt stress in the root, stem, and leaf tissues. In G. barbadense, the expression was higher in the stem, while in G. hirsutum the leaf and root showed a high expression. Moreover, our results revealed that NHX2 homologues in both species have a high expression under salinity stress at higher time intervals, followed by NHX7. The protein-protein prediction study revealed that GbNHX7 is involved in the CBL-CIPK protein interaction pathway. Our study also provided valuable information explaining the molecular mechanism of Na⁺ transport for the further functional study of Gossypium NHX genes.

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Crop molecular breeding primarily focuses on increasing the trait of plant yield. An elongator-associated protein, KTI12, is closely associated with plant biomass and yield. KTI12 is involved in developmental processes of most organs, including the leaf, root, flower, and seed, through regulating cell division and differentiation. Previous work has shown that in upland cotton (Gossypium hirsutum), GhKTI12 regulates plant height, flowering, and tolerance to salt and drought stress. However, little is known about the molecular regulation mechanism of GhKTI12 in plant developmental processes. In this study, we identified the main GhKTI12 (Gh_D02G144400) gene and transformed it into tobacco (Nicotonia tabacum cv NC89). From seven transgenic lines, we obtained three (OE5, OE6 and OE8) with high expression of GhKTI12; compared with wild type plants, these three lines exhibited larger plant size, later flowering, and higher seed yield. Microscopic observation revealed that the number of leaf epidermal cells and stem parenchyma cells was increased by ~55%. Biochemical analysis showed that chlorophyll content and starch accumulation were significantly increased in younger leaves at the top canopy of transgenic plants, which may contribute to improved photosynthetic rate and, in turn, increased seed yield. To understand the molecular mechanism of GhKTI12 in transgenic plants development, two lines (OE6 and OE8) with higher expression levels of GhKTI12 were used as representative plants to conduct RNA-seq analysis. Through transcriptome analysis of the plant’s shoot apical meristematic tissue of these two lines, we identified 518 upregulated genes and 406 downregulated genes common to both overexpression lines. A large number of cellular component genes associated with cell division and differentiation, such as RD21, TET8, KTN80, AOX1, AOX2, CP1, and KIC, were found to be upregulated, and genes showing the most downregulation included MADS-box genes related to flowering time, such as MADS6, AP1, AP3, AGL8, AGL6, SEP1, and SEP2. Downregulation of these genes caused delayed flowering time and longer vegetative stage during development. Combined with the upregulation of the yield-related gene RD21, the GhKTI12 transgenic plants could produce a higher seed yield. We here show that the overexpression of GhKTI12 could positively improve key agronomic traits in tobacco by regulating cell proliferation, photosynthesis, and organ development, and suggest that homologs of GhKTI12 may also be important in the genetic improvement of other crop plants.

Although genetic transformation has opened up a new era for cotton molecular breeding, it still suffers from the limitation problem of long transformation periods, which slows down the generation of new cotton germplasms. In this study, LT gene (SV40 large T antigen), which promotes the transformation efficiency of animal cells, was codon-optimized. Its overexpression vector was transformed into cotton. It was observed that EC (embryogenic callus) formation period was 33% shorter and transformation efficiency was slightly higher in the LT T0 generation than that of control. RNA-seq data of NEC (non-embryonic callus) and EC from LT and control revealed that more DEGs (differential expression genes) in NEC were identified than that of EC, indicating LT mainly functioned in NEC. Further KEGG, GO, and transcription factor analyses showed that DEGs were significantly enriched in brassinosteroid biosynthesis pathways and that bHLH, MYB, and AP2/ERF were the top three gene families, which are involved in EC formation. In addition, the key genes related to the auxin pathway were differentially expressed only in LT overexpression NEC, which caused early response, biosynthesis, and transportation of the hormone, resulting in EC earlier formation. In summary, the results demonstrated that LT can promote somatic embryogenesis in cotton, which provides a new strategy for improving cotton transformation and shortening EC formation time.

Cold stress adversely affects cotton growth and development often causing morphological abnormalities during germination and seedling establishment. Nevertheless, the molecular and metabolic mechanisms conferring cold tolerance in cotton remain largely unclear. The study aimed to examine the morphological and physiological responses of two semi-wild Gossypium hirsutum accessions, the cold resistant Latifolium (CR-307) and cold susceptible Punctatum (CS-131) under cold stress. The experiment was conducted at two temperature regimes: 28-30 °C (control) and 4 °C for 96 h (cold stress). In comparison to the control, primary root length (PRL) decreased by 63% in CR-307 and by 73% in CS-131. Furthermore, biochemical analysis confirmed that CR-307 exhibited a stronger antioxidant defense system, showing higher SOD, POD, and CAT activities and reduced H₂O₂ accumulation under cold stress compared with CS-131. Transcriptome analysis identified 3135 differentially expressed genes (DEGs) associated with cold response regulation, mainly belonging to MYB, AP2/ERF, bHLH, WRKY, bZIP, and HSF transcription factor families. Cold stress also altered pathways related to phenylpropanoid, proline, flavonoid, and cofactor biosynthesis, reflecting adaptive molecular reprogramming. Moreover, 328 lipids, 105 organic acids, 87 phenylpropanoids and 39 flavonoids were detected as differentially expressed metabolites (DEMs) that contribute to cell wall reinforcement, antioxidant defense, transcriptional regulation and membrane stabilization under cold stress. Overall, the study provides comprehensive insights into the molecular and metabolic mechanisms underlying cold tolerance in cotton, offering potential targets for developing cold-resilient cultivars.