Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Leaf shape is considered to be one of the most significant agronomic traits in crop breeding. However, the molecular basis underlying leaf morphogenesis in cotton is still largely unknown. In this study, through genetic mapping and molecular investigation using a natural cotton mutant cu with leaves curling upward, the causal gene GHCU is successfully identified as the key regulator of leaf flattening. Knockout of GHCU or its homolog in cotton and tobacco using CRISPR results in abnormal leaf shape. It is further discovered that GHCU facilitates the transport of the HD protein KNOTTED1-like (KNGH1) from the adaxial to the abaxial domain. Loss of GHCU function restricts KNGH1 to the adaxial epidermal region, leading to lower auxin response levels in the adaxial boundary compared to the abaxial. This spatial asymmetry in auxin distribution produces the upward-curled leaf phenotype of the cu mutant. By analysis of single-cell RNA sequencing and spatiotemporal transcriptomic data, auxin biosynthesis genes are confirmed to be expressed asymmetrically in the adaxial-abaxial epidermal cells. Overall, these findings suggest that GHCU plays a crucial role in the regulation of leaf flattening through facilitating cell-to-cell trafficking of KNGH1 and hence influencing the auxin response level.

The reniform nematode, Rotylenchulus reniformis, is a major yield-limiting pest of upland cotton (Gossypium hirsutum) in the United States that has been steadily increasing in incidence in many states. Reniform nematode-resistant cotton cultivars have recently become commercially available for cotton producers; however, few field trials have evaluated their efficacy as a nematode management tool. The aim of this study was to evaluate reniform nematode population development, plant growth, and seed cotton yield of reniform nematode-resistant cotton cultivars in two nematode-infested fields in Louisiana. Replicated small-plot field trials were conducted in St. Joseph, LA (NERS field) and Winnsboro, LA (MRRS field) during the 2022 and 2023 growing seasons. In 2022, cultivars evaluated included: (1) DP 1646 B2XF (susceptible/tolerant), (2) DP 2141NR B3XF (resistant), (3) PHY 332 W3FE (resistant), (4) PHY 411 W3FE (resistant), and (5) PHY 443 W3FE (resistant). In 2023, an additional susceptible cotton cultivar, PHY 340 W3FE, was also included. All nematode-resistant cotton cultivars evaluated provided suppression of reniform nematode population development relative to that of the susceptible cotton cultivars, with suppression of nematode soil population densities at harvest ranging from 49 - 81% relative to DP 1646 B2XF. The resistant cultivar PHY 411 W3FE provided the most consistent suppression of reniform nematode population development, reducing reniform nematode soil population densities at harvest in both field locations and both trial years. In contrast, DP 2141NR B3XF only reduced soil population densities at harvest in the NERS field in 2023. Despite relatively consistent nematode suppression and improvements in plant vigor ratings and canopy coverage associated with the resistant cotton cultivars, a yield increase was only observed with PHY 332 W3FE and PHY 411 W3FE planted at the NERS field in 2023. Despite strong resistance to reniform nematode in the evaluated cotton cultivars, nematode soil population densities still increased during the growing season in plots planted with resistant cotton cultivars, emphasizing the need for additional management tactics to use alongside host resistance. This study indicates that new reniform nematode-resistant cotton cultivars show promising potential to reduce nematode population development during the growing season in Louisiana.

Naturally-colored brown cotton (NBC) fiber is an environmentally friendly raw source of fiber for textile applications. The fiber of some NBC cultivars exhibits flame-retardant properties, which can be used in textiles that require flame resistance. Proanthocyanidins or their derivatives are responsible for the brown pigment in NBC; however, how flame retardancy is related to pigmentation in NBC is poorly understood. To gain insight into brown pigment biosynthesis, we conducted comparative transcripts and metabolites profiling analysis of developing cotton fibers between the brown (MC-BL) and white (MC-WL) cotton near-isogenic lines (NILs), genetically different only in the Lc1 locus. In this study, mass spectrometry was used to detect metabolites in BL and WL developing fibers at 8, 12, 16, 20, 24, 36, and 40 days post anthesis (DPA) and mature fibers. Transcripts analysis was performed at two critical fiber developmental points, 8 DPA (fiber elongation) and 20 DPA (secondary cell wall deposition). We found 5836 (ESI MS positive mode) and 4541 (ESI MS negative mode) metabolites significantly different accumulated between BL and WL. Among them, 142 were known non-redundant metabolites, including organic acids, amino acids, and derivatives of the phenylpropanoid pathway. Transcript analysis determined 1691 (8 DPA) and 5073 (20 DPA) differentially expressed genes (DEGs) between BL and WL, with the majority of DEGs down-regulated at 20 DPA. Organic acids of the citric acid cycle were induced, while most of the detected amino acids were reduced in the MC-BL line. Both cis- and trans-stereoisomers of flavan-3-ols were detected in developing MC-WL and MC-BL fibers; however, the gallocatechin and catechin accumulated multiple times higher. Gas chromatography-mass spectrometry (GC-MS) analysis of fatty acids determined that palmitic acid long-chain alcohols were the main constituents of waxes of mature fibers. Energy-dispersive X-ray spectrometry (EDS) analysis of mature fibers revealed that potassium accumulated three times greater in MC-BL than in MC-WL mature fibers. This study provides novel insights into the biosynthesis of pigments and its association with flame retardancy in NBC fibers.

Cotton is one of the most important natural fibers used in the textile industry worldwide. It is important to identify the key factors involved in cotton fiber development. In this study, zinc finger protein8 (GhZFP8) encoding a C2H2 transcription factor (TF) was cloned from cotton. qPCR showed that the transcripts of GhZFP8 in cotton were detected in the leaves and fibers at 3, 6, and 30 days post-anthesis (DPA), but not in the roots, stems, or flowers. The overexpression of GhZFP8 increased the trichome number on the siliques, leaves, and inflorescence, but inhibited the growth. The expression of trichome development and cell-elongation-related genes decreased obviously in GhZFP8 overexpressor Arabidopsis. Indole-3-acetic acid (IAA) and 1-Aminocyclopropanecarboxylic acid (ACC) contents were much higher in GhZFP8 overexpressors than that found in the wild type, but the gibberellin (GA) content was lower. The interference of GhZFP8 in cotton caused smaller bolls and shorter fibers than that of the control. The results of DNA affinity purification (DAP)-seq showed that GhZFP8 could bind to the promoter, exon, intron, and intergenic region of the target genes, which are involved in photosynthesis, signal transduction, synthesis of biomass, etc. Our findings implied that GhZFP8 processed multiple biological functions and regulated the development of cotton fiber.

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.

BACKGROUND: Universal stress proteins (USPs) are a class of stress-induced proteins that play a crucial role in biotic and abiotic stress responses. These proteins have previously been reported to participate directly in responses to various stress and protect plants against unfavorable environmental conditions. However, there is limited research on USPs in cotton, and systematic characterization of USPs in Gossypium species is lacking.
RESULTS: In the present study, the USP genes in Gossypium hirsutum were systematically identified and classified into six distinct subfamilies. The expansion of USPs in Gossypium species is mainly caused by dispersed duplication and whole genome duplication. Notably, the USPs that have expanded through allotetraploidization events are highly conserved in the allotetraploid species. The promoter regions of GhUSPs contain a diverse range of cis-acting elements associated with stress response. The RNA-Seq analysis and RT-qPCR assays revealed a significant induction of numerous GhUSPs expressions in response to various abiotic stresses. The co-expression network of GhUSPs revealed their involvement in stress response.
CONCLUSIONS: This study systematically analyzed the biological characteristics of GhUSPs and their response to abiotic stress. These findings serve as a theoretical basis for facilitating the breeding of cotton varieties in future research.

Substitution lines of the cotton Gossypium hirsutum L. involving chromosomes of the tetraploid species G. barbadense L., G. tomentosum Nutt. ex Seem., and G. mustelinum Miers ex Watt. are a valuable source for breeding, increasing the genetic diversity of G. hirsutum. The substitution of certain G. hirsutum L. chromosomes with G. barbadense chromosomes affect fibre elongation, fibre yield, fibre strength, and micronaire. To increase the efficiency of creating lines, it is necessary to study the nature of the introgression of alien chromosomes into the G. hirsutum L. genome. As a result of molecular genetic analysis of BC2F1 hybrids obtained from crossing monosomic lines of the cotton G. hirsutum from the cytogenetic collection of Uzbekistan with monosomic backcross hybrids BC1F1 G. hirsutum × G. barbadense on the same chromosomes, genetic differences between the hybrids in the profile of chromosome-specific microsatellite SSR markers were found. The predominant introgression of chromosomes 4, 6 and 12 of the At-subgenome and 22 of the Dt-subgenome of G. barbadense was revealed, while chromosomes 2 and 7 of the At-subgenome and 18 of the Dt- subgenome of G. barbadense were characterized by elimination. Among them, chromosomes 7 of the At- subgenome and 18 of the Dt-subgenome of G. barbadense were eliminated in the first backcross generation. In this work, two lines, CS- B06 and CS-B07, from the American cytogenetic collection with a putative substitution involving chromosomes 6 and 7 of the At-subgenome were analysed. The presence of only polymorphic alleles from the species G. hirsutum and the absence of polymorphic alleles from the species G. barbadense were revealed, which showed the absence of substitution involving these chromosomes. BC2F1 hybrids with monosomy for both G. barbadense and G. hirsutum chromosomes were characterized by regular pairing of chromosomes and high meiotic indexes. However, many hybrids were characterized by a decrease in pollen fertility. Two hybrids with monosomy for chromosome 7 of the At-subgenome of G. hirsutum and chromosome 6 of the At-subgenome of G. barbadense had the greatest reduction in pollen viability (70.09 ± 1.57 and 75.00 ± 1.66 %, respectively). Thus, this work shows a specific feature in the introgression of individual chromosomes of the cotton species G. barbadense into the cotton G. hirsutum genome.
Линии хлопчатника Gossypium hirsutum L. с чужеродным замещением хромосом тетраплоидных ви- дов G. barbadense L., G. tomentosum Nutt. ex Seem., G. mustelinum Miers ex Watt. являются ценным источником для селекции, увеличивающим генетическое разнообразие G. hirsutum. Замещение определенных хромосом хлоп- чатника вида G. hirsutum L. хромосомами вида G. barbadense оказывает влияние на удлинение, выход и прочность волокна, микронейр. Для повышения эффективности процесса создания линий необходимо изучение характера интрогрессии чужеродных хромосом в геном G. hirsutum L. В результате молекулярно-генетического анализа гибридов BC2F1, полученных от скрещиваний моносомных линий хлопчатника G. hirsutum цитогенетической коллекции Узбекистана с моносомными беккроссными гибридами BC1F1 G. hirsutum × G. barbadense по одинако- вым хромосомам, обнаружены генетические различия по профилю хромосом-специфичных микросателлитных SSR-маркеров между гибридами. Выявлена преимущественная интрогрессия хромосом 4, 6, 12 At-субгенома и 22 Dt-субгенома G. barbadense, тогда как хромосомы 2, 7 At-субгенома и 18 Dt-субгенома G. barbadense характе- ризовались элиминацией, среди них хромосомы 7 At-субгенома и 18 Dt-субгенома G. barbadense элиминировали уже в первом беккроссном поколении. В настоящей работе проанализированы две линии, CS-B06 и CS-B07, аме- риканской цитогенетической коллекции с предполагаемым замещением по хромосомам 6 и 7 Аt-субгенома. Об- наружены присутствие только полиморфных аллелей вида G. hirsutum и отсутствие полиморфных аллелей вида G. barbadense, что показало отсутствие замещения по этим хромосомам. Гибриды BC2F1 с моносомией как по хромосомам G. barbadense, так и по хромосомам G. hirsutum характеризовались регулярной конъюгацией хромо- сом и высоким мейотическим индексом. Однако многие гибриды отличались снижением фертильности пыльцы. Два гибрида с моносомией по хромосоме 7 At-субгенома G. hirsutum и хромосоме 6 At-субгенома G. barbadense имели наибольшую редукцию в жизнеспособности пыльцы (70.09 ± 1.57 и 75.00 ± 1.66 % соответственно). Таким образом, в этой работе показана особенность в интрогрессии индивидуальных хромосом хлопчатника вида G. barbadense в геном хлопчатника G. hirsutum.

A-galactosidases (AGALs), the oligosaccharide (RFO) catabolic genes of the raffinose family, play crucial roles in plant growth and development and in adversity stress. They can break down the non-reducing terminal galactose residues of glycolipids and sugar chains. In this study, the whole genome of AGALs was analyzed. Bioinformatics analysis was conducted to analyze members of the AGAL family in Gossypium hirsutum, Gossypium arboreum, Gossypium barbadense, and Gossypium raimondii. Meanwhile, RT-qPCR was carried out to analyze the expression patterns of AGAL family members in different tissues of terrestrial cotton. It was found that a series of environmental factors stimulated the expression of the GhAGAL3 gene. The function of GhAGAL3 was verified through virus-induced gene silencing (VIGS). As a result, GhAGAL3 gene silencing resulted in milder wilting of seedlings than the controls, and a significant increase in the raffinose content in cotton, indicating that GhAGAL3 responded to NaCl stress. The increase in raffinose content improved the tolerance of cotton. Findings in this study lay an important foundation for further research on the role of the GhAGAL3 gene family in the molecular mechanism of abiotic stress resistance in cotton.

The cotton fleahopper (Pseudatomoscelis seriatus Reuter) is considered a highly economically damaging pest of cotton (Gossypium hirsutum L.) in Texas and Oklahoma. Current control methods rely heavily on the use of foliar-applied chemical insecticides, but considering the cost of insecticides and the critical timeliness of applications, chemical control methods are often not optimized to reduce potential yield losses from this pest. The Bacillus thuringiensis (Bt) Mpp51Aa2 (formerly Cry51Aa2.834_16) protein has proven effective against thrips and plant bugs with piercing and sucking feeding behaviors, but the impact of this toxin on cotton fleahoppers has not been investigated. To evaluate the Mpp51Aa2 trait effectiveness towards the cotton fleahopper, field trials were conducted in 2019, 2020, and 2021, comparing a cotton cultivar containing the Mpp51Aa2 trait to a non-traited isoline cultivar under insecticide-treated and untreated conditions. Populations of cotton fleahopper nymphs and adults were estimated weekly by visually inspecting cotton terminals. Square retention was also assessed during the first week of bloom to provide some insight on how the Bt trait may influence yield. While cotton fleahopper population differences between the traited and non-traited plants were not consistently noted during the pre-bloom squaring period, there was a consistent increase in square retention in cotton expressing Mpp51Aa2 relative to non-traited cotton. Additionally, cotton expressing Mpp51Aa2 offered similar square protection relative to non-traited cotton treated with insecticides for the cotton fleahopper. These findings indicate that the Mpp51Aa2 protein should provide benefits of delayed nymphal growth, population suppression, and increased square retention.