Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton.

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

Comprehensive identification and systematical characterization of BRX gene family and the functional of GhBRXL5A in response to salt stress.

BACKGROUND: BRVIS RADIX (BRX) family is a small gene family with the highly conserved plant-specific BRX domains, which plays important roles in plant development and response to abiotic stress. Although BRX protein has been studied in other plants, the biological function of cotton BRX-like (BRXL) gene family is still elusive. RESULT: In this study, a total of 36 BRXL genes were identified in four cotton species. Whole genome or segmental duplications played the main role in the expansion of GhBRXL gene family during evolutionary process in cotton. These BRXL genes were clustered into 2 groups, α and ß, in which structural and functional conservation within same groups but divergence among different groups were found. Promoter analysis indicated that cis-elements were associated with the phytohormone regulatory networks and the response to abiotic stress. Transcriptomic analysis indicated that GhBRXL2A/2D and GhBRXL5A/5D were up/down-regulated in response to the different stress. Silencing of GhBRXL5A gene via virus-induced gene silencing (VIGS) improved salt tolerance in cotton plants. Furthermore, yeast two hybrid analysis suggested homotypic and heterotypic interactions between GhBRXL1A and GhBRXL5D. CONCLUSIONS: Overall, these results provide useful and valuable information for understanding the evolution of cotton GhBRXL genes and their functions in salt stress.

Exploring cotton SFR2's conundrum in response to cold stress.

Cotton is an important agricultural crop to many regions across the globe but is sensitive to low-temperature exposure. The activity of the enzyme SENSITIVE TO FREEZING 2 (SFR2) improves cold tolerance of plants and produces trigalactosylsyldiacylglycerol (TGDG), but its role in cold sensitive plants, such as cotton remains unknown. Recently, it was reported that cotton SFR2 produced very little TGDG under normal and cold conditions. Here, we investigate cotton SFR2 activation and TGDG production. Using multiple approaches in the native system and transformation into Arabidopsis thaliana, as well as heterologous yeast expression, we provide evidence that cotton SFR2 activates differently than previously found among other plant species. We conclude with the hypothesis that SFR2 in cotton is not activated in a similar manner regarding acidification or freezing like Arabidopsis and that other regions of SFR2 protein are critical for activation of the enzyme than previously reported.

GhGASA14 regulates the flowering time of upland cotton in response to GA3.

KEY MESSAGE: The silencing of GhGASA14 and the identification of superior allelic variation in its coding region indicate that GhGASA14 may positively regulate flowering and the response to GA3. Gibberellic acid-stimulated Arabidopsis (GASA), a member of the gibberellin-regulated short amino acid family, has been extensively investigated in several plant species and found to be critical for plant growth and development. However, research on this topic in cotton has been limited. In this study, we identified 38 GhGASAs that were dispersed across 18 chromosomes in upland cotton, and all of these genes had a GASA core domain. Transcriptome expression patterns and qRT-PCR results revealed that GhGASA9 and GhGASA14 exhibited upregulated expression not only in the floral organs but also in the leaves of early-maturing cultivars. The two genes were functionally characterized by virus-induced gene silencing (VIGS), and the budding and flowering times after silencing the target genes were later than those of the control (TRV:00). Compared with that in the water-treated group (MOCK), the flowering period of the different fruiting branches in the GA3-treated group was more concentrated. Interestingly, allelic variation was detected in the coding sequence of GhGASA14 between early-maturing and late-maturing accessions, and the frequency of this favorable allele was greater in high-latitude cotton cultivars than in low-latitude ones. Additionally, a significant linear relationship was observed between the expression level of GhGASA14 and flowering time among the 12 upland cotton accessions. Taken together, these results indicated that GhGASA14 may positively regulate flowering time and respond to GA3. These findings could lead to the use of valuable genetic resources for breeding early-maturing cotton cultivars in the future.

A cotton mitochondrial alternative electron transporter, GhD2HGDH, induces early flowering by modulating GA and photoperiodic pathways.

D-2-hydroxyglutarate dehydrogenase (D2HGDH) is a mitochondrial enzyme containing flavin adenine dinucleotide FAD, existing as a dimer, and it facilitates the specific oxidation of D-2HG to 2-oxoglutarate (2-OG), which is a key intermediate in the tricarboxylic acid (TCA) cycle. A Genome-wide expression analysis (GWEA) has indicated an association between GhD2HGDH and flowering time. To further explore the role of GhD2HGDH, we performed a comprehensive investigation encompassing phenotyping, physiology, metabolomics, and transcriptomics in Arabidopsis thaliana plants overexpressing GhD2HGDH. Transcriptomic and qRT-PCR data exhibited heightened expression of GhD2HGDH in upland cotton flowers. Additionally, early-maturing cotton exhibited higher expression of GhD2HGDH across all tissues than delayed-maturing cotton. Subcellular localization confirmed its presence in the mitochondria. Overexpression of GhD2HGDH in Arabidopsis resulted in early flowering. Using virus-induced gene silencing (VIGS), we investigated the impact of GhD2HGDH on flowering in both early- and delayed-maturing cotton plants. Manipulation of GhD2HGDH expression levels led to changes in photosynthetic pigment and gas exchange attributes. GhD2HGDH responded to gibberellin (GA3) hormone treatment, influencing the expression of GA biosynthesis genes and repressing DELLA genes. Protein interaction studies, including yeast two-hybrid, luciferase complementation (LUC), and GST pull-down assays, confirmed the interaction between GhD2HGDH and GhSOX (Sulfite oxidase). The metabolomics analysis demonstrated GhD2HGDH's modulation of the TCA cycle through alterations in various metabolite levels. Transcriptome data revealed that GhD2HGDH overexpression triggers early flowering by modulating the GA3 and photoperiodic pathways of the flowering core factor genes. Taken together, GhD2HGDH positively regulates the network of genes associated with early flowering pathways.

Investigation of salt tolerance in cotton germplasm by analyzing agro-physiological traits and ERF genes expression.

The development of genotypes that can tolerate high levels of salt is crucial for the efficient use of salt-affected land and for enhancing crop productivity worldwide. Therefore, incorporating salinity tolerance is a critical trait that crops must possess. Salt resistance is a complex character, controlled by multiple genes both physiologically and genetically. To examine the genetic foundation of salt tolerance, we assessed 16 F1 hybrids and their eight parental lines under normal and salt stress (15 dS/m) conditions. Under salt stress conditions significant reduction was observed for plant height (PH), bolls/plant (NBP), boll weight (BW), seed cotton yield (SCY), lint% (LP), fiber length (FL), fiber strength (FS), potassium to sodium ratio (K+/Na+), potassium contents (K+), total soluble proteins (TSP), carotenoids (Car) and chlorophyll contents. Furthermore, the mean values for hydrogen peroxide (H2O2), sodium contents (Na+), catalase (CAT), superoxide dismutase (SOD), peroxidase (POD), and fiber fineness (FF) were increased under salt stress. Moderate to high heritability and genetic advancement was observed for NBP, BW, LP, SCY, K+/Na+, SOD, CAT, POD, Car, TSP, FL, and FS. Mean performance and multivariate analysis of 24 cotton genotypes based on various agro-physiological and biochemical parameters suggested that the genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, along with crosses FBS-FALCON × JSQ-White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 achieved the maximum values for K+/Na+, K+, TSP, POD, Chlb, CAT, Car, LP, FS, FL, PH, NBP, BW, and SCY under salt stress and declared as salt resistant genotypes. The above-mentioned genotypes also showed relatively higher expression levels of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 at 15 dS/m and proved the role of these ERF genes in salt tolerance in cotton. These findings suggest that these genotypes have the potential for the development of salt-tolerant cotton varieties with desirable fiber quality traits.

GhVIM28, a negative regulator identified from VIM family genes, positively responds to salt stress in cotton.

The VIM (belonged to E3 ubiquitin ligase) gene family is crucial for plant growth, development, and stress responses, yet their role in salt stress remains unclear. We analyzed phylogenetic relationships, chromosomal localization, conserved motifs, gene structure, cis-acting elements, and gene expression patterns of the VIM gene family in four cotton varieties. Our findings reveal 29, 29, 17, and 14 members in Gossypium hirsutum (G.hirsutum), Gossypium barbadense (G.barbadense), Gossypium arboreum (G.arboreum), and Gossypium raimondii (G. raimondii), respectively, indicating the maturity and evolution of this gene family. motifs among GhVIMs genes were observed, along with the presence of stress-responsive, hormone-responsive, and growth-related elements in their promoter regions. Gene expression analysis showed varying patterns and tissue specificity of GhVIMs genes under abiotic stress. Silencing GhVIM28 via virus-induced gene silencing revealed its role as a salt-tolerant negative regulator. This work reveals a mechanism by which the VIM gene family in response to salt stress in cotton, identifying a potential negative regulator, GhVIM28, which could be targeted for enhancing salt tolerance in cotton. The objective of this study was to explore the evolutionary relationship of the VIM gene family and its potential function in salt stress tolerance, and provide important genetic resources for salt tolerance breeding of cotton.

GhTKPR1_8 functions to inhibit anther dehiscence and reduce pollen viability in cotton.

Sporopollenin, as the main component of the pollen exine, is a highly resistant polymer that provides structural integrity under unfavourable environmental conditions. Tetraketone α-pyrone reductase 1 (TKPR1) is essential for sporopollenin formation, catalyzing the reduction of tetraketone carbonyl to hydroxylated α-pyrone. The functional role of TKPR1 in male sterility has been reported in flowering plants such as maize, rice, and Arabidopsis. However, the molecular cloning and functional characterization of TKPR1 in cotton remain unaddressed. In this study, we identified 68 TKPR1s from four cotton species, categorized into three clades. Transcriptomics and RT-qPCR demonstrated that GhTKPR1_8 exhibited typical expression patterns in the tetrad stage of the anther. GhTKPR1_8 was localized to the endoplasmic reticulum. Moreover, ABORTED MICROSPORES (GhAMS) transcriptionally activated GhTKPR1_8 as indicated by luciferase complementation tests. GhTKPR1_8-knockdown inhibited anther dehiscence and reduced pollen viability in cotton. Additionally, overexpression of GhTKPR1_8 in the attkpr1 mutant restored its male sterile phenotype. This study offers novel insights into the investigation of TKPR1 in cotton while providing genetic resources for studying male sterility.

Decoding drought resilience: a comprehensive exploration of the cotton Eceriferum (CER) gene family and its role in stress adaptation.

BACKGROUND: The cuticular wax serves as a primary barrier that protects plants from environmental stresses. The Eceriferum (CER) gene family is associated with wax production and stress resistance. RESULTS: In a genome-wide identification study, a total of 52 members of the CER family were discovered in four Gossypium species: G. arboreum, G. barbadense, G. raimondii, and G. hirsutum. There were variations in the physicochemical characteristics of the Gossypium CER (GCER) proteins. Evolutionary analysis classified the identified GCERs into five groups, with purifying selection emerging as the primary evolutionary force. Gene structure analysis revealed that the number of conserved motifs ranged from 1 to 15, and the number of exons varied from 3 to 13. Closely related GCERs exhibited similar conserved motifs and gene structures. Analyses of chromosomal positions, selection pressure, and collinearity revealed numerous fragment duplications in the GCER genes. Additionally, nine putative ghr-miRNAs targeting seven G. hirsutum CER (GhCER) genes were identified. Among them, three miRNAs, including ghr-miR394, ghr-miR414d, and ghr-miR414f, targeted GhCER09A, representing the most targeted gene. The prediction of transcription factors (TFs) and the visualization of the regulatory TF network revealed interactions with GhCER genes involving ERF, MYB, Dof, bHLH, and bZIP. Analysis of cis-regulatory elements suggests potential associations between the CER gene family of cotton and responses to abiotic stress, light, and other biological processes. Enrichment analysis demonstrated a robust correlation between GhCER genes and pathways associated with cutin biosynthesis, fatty acid biosynthesis, wax production, and stress response. Localization analysis showed that most GCER proteins are localized in the plasma membrane. Transcriptome and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) expression assessments demonstrated that several GhCER genes, including GhCER15D, GhCER04A, GhCER06A, and GhCER12D, exhibited elevated expression levels in response to water deficiency stress compared to control conditions. The functional identification through virus-induced gene silencing (VIGS) highlighted the pivotal role of the GhCER04A gene in enhancing drought resistance by promoting increased tissue water retention. CONCLUSIONS: This investigation not only provides valuable evidence but also offers novel insights that contribute to a deeper understanding of the roles of GhCER genes in cotton, their role in adaptation to drought and other abiotic stress and their potential applications for cotton improvement.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton.

BACKGROUND: Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS: In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS: These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

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.

GhiPLATZ17 and GhiPLATZ22, zinc-dependent DNA-binding transcription factors, promote salt tolerance in upland cotton.

KEY MESSAGE: The utilization of transcriptome analysis, functional validation, VIGS, and DAB techniques have provided evidence that GhiPLATZ17 and GhiPLATZ22 play a pivotal role in improving the salt tolerance of upland cotton. PLATZ (Plant AT-rich sequences and zinc-binding proteins) are known to be key regulators in plant growth, development, and response to salt stress. In this study, we comprehensively analyzed the PLATZ family in ten cotton species in response to salinity stress. Gossypium herbaceum boasts 25 distinct PLATZ genes, paralleled by 24 in G. raimondii, 25 in G. arboreum, 46 in G. hirsutum, 48 in G. barbadense, 43 in G. tomentosum, 67 in G. mustelinum, 60 in G. darwinii, 46 in G. ekmanianum, and a total of 53 PLATZ genes attributed to G. stephensii. The PLATZ gene family shed light on the hybridization and allopolyploidy events that occurred during the evolutionary history of allotetraploid cotton. Ka/Ks analysis suggested that the PLATZ gene family underwent intense purifying selection during cotton evolution. Analysis of synteny and gene collinearity revealed a complex pattern of segmental and dispersed duplication events to expand PLATZ genes in cotton. Cis-acting elements and gene expressions revealed that GhiPLATZ exhibited salt stress resistance. Transcriptome analysis, functional validation, virus-induced gene silencing (VIGS), and diaminobenzidine staining (DAB) demonstrated that GhiPLATZ17 and GhiPLATZ22 enhance salt tolerance in upland cotton. The study can potentially advance our understanding of identifying salt-resistant genes in cotton.

Partial root-zone drying irrigation improves intrinsic water-use efficiency and maintains high photosynthesis by uncoupling stomatal and mesophyll conductance in cotton leaves.

Partial root-zone drying irrigation (PRD) can improve water-use efficiency (WUE) without reductions in photosynthesis; however, the mechanism by which this is attained is unclear. To amend that, PRD conditions were simulated by polyethylene glycol 6000 in a root-splitting system and the effects of PRD on cotton growth were studied. Results showed that PRD decreased stomatal conductance (gs) but increased mesophyll conductance (gm). Due to the contrasting effects on gs and gm, net photosynthetic rate (AN) remained unaffected, while the enhanced gm/gs ratio facilitated a larger intrinsic WUE. Further analyses indicated that PRD-induced reduction of gs was related to decreased stomatal size and stomatal pore area in adaxial and abaxial surface which was ascribed to lower pore length and width. PRD-induced variation of gm was ascribed to the reduced liquid-phase resistance, due to increases in chloroplast area facing to intercellular airspaces and the ratio of chloroplast surface area to total mesophyll cell area exposed to intercellular airspaces, as well as to decreases in the distance between cell wall and chloroplast, and between adjacent chloroplasts. The above results demonstrate that PRD, through alterations to stomatal and mesophyll structures, decoupled gs and gm responses, which ultimately increased intrinsic WUE and maintained AN.

Yield-related quantitative trait loci identification and lint percentage hereditary dissection under salt stress in upland cotton.

Salinity is frequently mentioned as a major constraint in worldwide agricultural production. Lint percentage (LP) is a crucial yield-component in cotton lint production. While the genetic factors affect cotton yield in saline soils are still unclear. Here, we employed a recombinant inbred line population in upland cotton (Gossypium hirsutum L.) and investigated the effects of salt stress on five yield and yield component traits, including seed cotton yield per plant, lint yield per plant, boll number per plant, boll weight, and LP. Between three datasets of salt stress (E1), normal growth (E2), and the difference values dataset of salt stress and normal conditions (D-value), 87, 82, and 55 quantitative trait loci (QTL) were detectable, respectively. In total, five QTL (qLY-Chr6-2, qBNP-Chr4-1, qBNP-Chr12-1, qBNP-Chr15-5, qLP-Chr19-2) detected in both in E1 and D-value were salt related QTL, and three stable QTL (qLP-Chr5-3, qLP-Chr13-1, qBW-Chr5-5) were detected both in E1 and E2 across 3 years. Silencing of nine genes within a stable QTL (qLP-Chr5-3) highly expressed in fiber developmental stages increased LP and decreased fiber length (FL), indicating that multiple minor-effect genes clustered on Chromosome 5 regulate LP and FL. Additionally, the difference in LP caused by Gh_A05G3226 is mainly in transcription level rather than in the sequence difference. Moreover, silencing of salt related gene (GhDAAT) within qBNP-Chr4-1 decreased salt tolerance in cotton. Our findings shed light on the regulatory mechanisms underlining cotton salt tolerance and fiber initiation.

Drought decreases cotton fiber strength by altering sucrose flow route.

The potential mechanisms by which drought restricts cotton fiber cell wall synthesis and fiber strength are still not fully understood. Herein, drought experiments were conducted using two cultivars of upland cotton (Gossypium hirsutum), Dexiamian 1 (drought-tolerant) and Yuzaomian 9110 (drought-sensitive). Results showed that drought notably reduced sucrose efflux from cottonseed coats to fibers by down-regulating the expression of GhSWEET10 and GhSWEET15 in outer cottonseed coats, leading to enhanced sucrose accumulation in cottonseed coats but decreased sucrose accumulation in fibers. Within cotton fibers, drought restricted the hydrolysis of sucrose to uridine-5'-diphosphoglucose by suppressing sucrose synthase activity, and drought favored the conversion of uridine-5'-diphosphoglucose to ß-1,3-glucan rather than cellulose by up-regulating GhCALS5. Hence, cellulose content was reduced, which was the main reason for the decreased fiber strength under drought. Moreover, drought promoted lignin synthesis by up-regulating the expression of Gh4CL4, GhPAL9, GhCCR5, GhCAD11, and GhCOMT6, which partly offset the negative influence of reduced cellulose content on fiber strength. Compared with Yuzaomian 9110, the drought-tolerance of Dexiamian 1 was evidenced by the following under drought conditions: (i) greater sucrose flow from seedcoat to fiber, (ii) less ß-1,3-glucan accumulation, and (iii) more lignin biosynthesis. Overall, this study provides new insights into the mechanism of reduced cotton fiber strength induced by drought.

Changes in carbohydrate distribution in cotton photosynthetic organs and increase in boll weight reduce yield loss under high temperature.

Yield of cotton (Gossypium hirsutum) does not always fall with high temperature (HT) even though this induces significant reductions in fruit retention. To investigate the underlying mechanisms, a greenhouse experiment was conducted with two temperature regimes [control treatment, 28 °C; high temperature (HT), 34 °C] for 7 d. Results showed HT did not significantly influence cotton yield, but reduced boll number and increased boll weight. The 13C distribution ratio of the leaf subtending the cotton boll (LSCB) decreased while that of the cotton boll increased under HT. Transcriptomic and proteomic analyses of the LSCB revealed up-regulated genes involved in cytokinin and jasmonic acid synthesis, as well as SWEET15 (GH_D01G0218), which positively regulated photosynthesis and transport photosynthate, ultimately leading to increased boll weight. After 7 d recovery from HT, the 13C distribution ratio of the LSCB increased while that of the cotton boll decreased. However, boll weight still increased, which was related to increased amylase and sucrose phosphate synthase activities and up-regulated sucrose transport genes in the main-stem leaf and capsule wall. Thus, both accelerated sucrose synthesis and transport in the LSCB under HT and increased sucrose supply ability of the main-stem leaf and capsule wall after recovery from HT contributed to an increased boll weight, which finally maintained cotton yield.

Single-cell transcriptome atlas reveals somatic cell embryogenic differentiation features during regeneration.

Understanding somatic cell totipotency remains a challenge facing scientific inquiry today. Plants display remarkable cell totipotency expression, illustrated by single-cell differentiation during somatic embryogenesis (SE) for plant regeneration. Determining cell identity and exploring gene regulation in such complex heterogeneous somatic cell differentiation have been major challenges. Here, we performed high-throughput single-cell sequencing assays to define the precise cellular landscape and revealed the modulation mode of marker genes during embryogenic differentiation in cotton (Gossypium hirsutum L.) as the crop for biotechnology application. We demonstrated that nonembryogenic calli (NEC) and primary embryogenic calli (PEC) tissues were composed of heterogeneous cells that could be partitioned into four broad populations with six distinct cell clusters. Enriched cell clusters and cell states were identified in NEC and PEC samples, respectively. Moreover, a broad repertoire of new cluster-specific genes and associated expression modules were identified. The energy metabolism, signal transduction, environmental adaptation, membrane transport pathways, and a series of transcription factors were preferentially enriched in cell embryogenic totipotency expression. Notably, the SE-ASSOCIATED LIPID TRANSFER PROTEIN (SELTP) gene dose-dependently marked cell types with distinct embryogenic states and exhibited a parabolic curve pattern along the somatic cell embryogenic differentiation trajectory, suggesting that SELTP could serve as a favorable quantitative cellular marker for detecting embryogenic expression at the single-cell level. In addition, RNA velocity and Scissor analysis confirmed the pseudo-temporal model and validated the accuracy of the scRNA-seq data, respectively. This work provides valuable marker-genes resources and defines precise cellular taxonomy and trajectory atlases for somatic cell embryogenic differentiation in plant regeneration.

Screening cotton genotypes for their drought tolerance ability based on the expression level of dehydration-responsive element-binding protein and proline biosynthesis-related genes and morpho-physio-biochemical responses.

Drought stress adversely affects growth, development, productivity, and fiber quality of cotton (Gossypium hirsutum L). Breeding strategies to enhance drought tolerance require an improved knowledge of plant drought responses necessitating proper identification of drought-tolerant genotypes of crops, including cotton. The objective of this study was to classify the selected cotton genotypes for their drought tolerance ability based on morpho-physio-biochemical traits using Hierarchical Ward's cluster analysis. Five genotypes of cotton (Takfa 3, Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5) were selected as plant materials, and were grown under well-watered (WW; 98 ± 2% field capacity) and water-deficit (WD; 50 ± 2% field capacity) conditions for 16 days during the flower initiation stage. Data on morpho-physio-biochemical parameters and gene expression levels for these parameters were collected, and subsequently genotypes were classified either as a drought tolerant or drought susceptible one. Upregulation of GhPRP (proline-rich protein), GhP5CS (Δ1-pyrroline-5-carboxylate synthetase), and GhP5CR (Δ1-pyrroline-5-carboxylate reductase) in relation to free proline enrichment was observed in Takfa 3 genotype under WD condition. An accumulation of free proline, total soluble sugar, and potassium in plants under WD conditions was detected, which played a key role as major osmolytes controlling cellular osmotic potential. Magnesium and calcium concentrations were also enriched in leaves under WD conditions, functioning as essential elements and regulating photosynthetic abilities. Leaf greenness, net photosynthetic rate, stomatal conductance, and transpiration rate were also declined under WD conditions, leading to growth retardation, especially aboveground traits of Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5 genotypes. An increase in leaf temperature (1.1 - 4.0 °C) and crop water stress index (CWSI > 0.75) in relation to stomatal closure and reduced transpiration rate was recorded in cotton genotypes under WD conditions compared with WW conditions. Based on the increase of free proline, soluble sugar, leaf temperature, and CWSI, as well as the decrease of aboveground growth traits and physiological attributes, five genotypes were categorized into two cluster groups: drought tolerant (Takfa 3) and drought susceptible (Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5). The identified drought-tolerant cotton genotype, namely, Takfa 3, may be grown in areas experiencing drought conditions. It is recommended to further validate the yield traits of Takfa 3 under rainfed field conditions in drought-prone environments.

Global gene expression profile and functional analysis reveal the conservation of reproduction-associated gene networks in Gossypium hirsutum.

KEY MESSAGE: Lastly, the bZIP gene family encompasses genes that have been reported to play a role in flower development, such as bZIP14 (FD). Notably, bZIP14 is essential for Flowering Locus T (FT) initiation of floral development in Arabidopsis (Abe et al. 2005). Cotton (Gossypium hirsutum L.) is the world's most extensively cultivated fiber crop. However, its reproductive development is poorly characterized at the molecular level. Thus, this study presents a detailed transcriptomic analysis of G. hirsutum at three different reproductive stages. We provide evidence that more than 64,000 genes are active in G. hirsutum during flower development, among which 94.33% have been assigned to functional terms and specific pathways. Gene set enrichment analysis (GSEA) revealed that the biological process categories of floral organ development, pollen exine formation, and stamen development were enriched among the genes expressed during the floral development of G. hirsutum. Furthermore, we identified putative Arabidopsis homologs involved in the G. hirsutum gene regulatory network (GRN) of pollen and flower development, including transcription factors such as WUSCHEL (WUS), INNER NO OUTER (INO), AGAMOUS-LIKE 66 (AGL66), SPOROCYTELESS/NOZZLE (SPL/NZZ), DYSFUNCTIONAL TAPETUM 1 (DYT1), ABORTED MICROSPORES (AMS), and ASH1-RELATED 3 (ASHR3), which are known crucial genes for plant reproductive success. The cotton MADS-box protein-protein interaction pattern resembles the previously described patterns for AGAMOUS (AG), SEEDSTICK (STK), SHATTERPROOF (SHP), and SEPALLATA3 (SEP3) homolog proteins from Arabidopsis. In addition to serving as a resource for comparative flower development studies, this work highlights the changes in gene expression profiles and molecular networks underlying stages that are valuable for cotton breeding improvement.