1. Introduction
Cotton (
Gossypium spp.) is a globally cultivated and economically significant fiber crop, playing a pivotal role in the textile industry and making substantial contributions to the global economy. The quality of cotton fibers directly influences the excellence of cotton-based textiles. Derived from single-celled epidermal seed trichomes produced by four domesticated species within the
Gossypium genus, cotton stands as the most widely utilized plant textile globally. The genus, encompassing over 50 species spread across the tropics and subtropics, exhibits rich species diversity and morphological variation [
1]. This diversity, especially in the context of the domestication experiment involving four species and two ploidy states within
Gossypium, serves as a natural and powerful system for studying genetic, genomic, and genotype-to-phenotype transitions. Notably,
Gossypium hirsutum cotton exhibits exceptional adaptability to various environmental conditions and is renowned for its capacity to yield substantial quantities of high-quality fiber, constituting approximately 95% of all planted cotton [
2], while
Gossypium barbadense produces luxury textiles with fine, lengthy, and robust fibers [
3]. Additionally,
Gossypium hirsutum requires fewer resources and displays enhanced resistance to pests and diseases compared to
Gossypium barbadense. Despite domesticated diploids producing inferior fibers, they persist in local cultivation due to their adaptation to specific regional conditions [
3].
Each cotton fiber is a single, elongated cotton seed coat epidermal cell. The differentiation and development of cotton fiber cells are highly complex, divided into five overlapping stages: fiber initiation [days post-anthesis (DPA), -1~3 DPA], fiber cell elongation (1~16 DPA), transition (16~20 DPA), secondary cell wall (SCW) synthesis (20~40 DPA), and fiber cell maturation (40~50 DAP) [
4,
5,
6] (
Figure 1). Transition is a distinct developmental stage that occurs between primary and secondary wall synthesis. Over several days, a stage-specific transcriptome underpins the distinctive cellular and biochemical status of the fiber cell. The galacturonosyltransferases (GAUTs) gene family, which is a critical participant in the pectin synthesis pathway, plays an important role in elongation, transition, and cell wall synthesis [
7]. The often mentioned ‘maturation’ phase of fiber development is not well characterized. After SCW cellulose synthesis stops, the fiber continues to dehydrate within the closed boll and the lateral packing of cellulose increases (
Figure 1A). Morphological and embryological studies, coupled with advanced microscopy techniques, have significantly contributed to our understanding of cotton fiber development. In
Gossypium hirsutum, long lint fibers primarily initiate before or on the day of anthesis, while ovule epidermal cells starting at or after 3 DPA produce shorter fibers known as linters or fuzz (
Figure 1B). Fiber mutants with inhibited fiber initial growth have demonstrated that the first detectable sign of fiber initiation occurs on the day of anthesis (0 DPA), with approximately 25% of epidermal cells contributing to fiber initials [
8]. Vibrational sum frequency generation, attenuated total refection infrared (ATR-IR), Fourier transform Raman (FT-Raman) spectroscopy and X-ray diffraction (XRD) were all used to study the cellulose component in mature and air-dry fibers from two species,
Gossypium hirsutum and
Gossypium barbadense [
9].
Understanding the intricate processes of initiation and elongation in cotton fibers involves the coordination of multiple genes and pathways. The integration of next generation sequencing technology, including mRNA sequencing and genome sequencing with innovative phenotype evaluation methods, has emerged as a transformative tool in comprehending the biological mechanisms governing cotton fiber development. For example, the 102
TBL (
Trichome Birefringence Like) gene family members were evaluated for their involvement in fiber development from
Gossypium hirsutum with two
GhTBL genes (
GhTBL7 and
GhTBL58) showing differential expression at 10 DPA fiber. Silencing these genes resulted in a significant reduction in the fiber length at 10 DPA, suggesting a potential role for these genes in fiber elongation [
10]. Recently, a total of 143
GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in
Gossypium hirsutum. A gene expression profile provided essential clues for the function of
GhHMA genes in cotton fiber development and response to various abiotic stresses.
GhHMA26 was predominantly expressed in 10 DPA fiber cells and the relative expression was higher than any other
GhHMAs which indicate that
GhHMA26 may positively regulate fiber elongation [
11]. Genome-wide analysis of the serine carboxypeptidase-like protein family reveals
Ga09G1039 is involved in fiber elongation in cotton with overexpression of
Ga09G1039 significantly increasing the length of stem trichomes [
12]. RNA sequencing (RNA-seq) data and N6-methyladenosine sequencing (m6A-seq) data showed that methylation of m6A affected the mRNA stability of these fiber elongation-related genes including the transcription factor
GhMYB44. Overexpression of GhMYB44 reduced fiber elongation, whereas the silencing of
GhMYB44 produced longer fibers [
13]. Comparative phosphor-proteomic analysis between two cotton varieties,
J7-1 and
J14-1, revealed that phosphorylation of sucrose synthase
GhSUS2 by Ca
2+ dependent protein kinases
GhCPK84/93 affects cotton fiber development. Moreover, ABA could promote the transcription and translation of
GhCPK84 and
GhCPK93, thereby enhancing the phosphorylation of GhSUS2 to impede fiber elongation [
14]. A genome-wide association study (GWAS) identified a ~6.2 kb insertion, larINDELFZ, positioned at the end of chromosome 8 in fuzzless
Gossypium arboreum. Comprised of a ~5.0 kb repetitive sequence and a ~1.2 kb fragment translocated from chromosome 12, this remote insertion was predicted to function as an enhancer located ~18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). This finding unveiled a novel regulator of fiber/trichome development, shedding light on the significance of noncoding sequences in cotton. The large-fragment insertion activates the GaFZ gene and is associated with fuzz and trichome reduction in
Gossypium arboreum [
15]. Genome-wide exploration identified a total of 125, 73 and 71 full-length
Catharanthus roseus receptor-like kinase 1-like (
CrRLK1L) family genes in
Gossypium hirsutum,
Gossypium arboreum and
Gossypium raimondii, respectively, with some of
GhCrRLK1Ls preferentially expressed in fibers at the different stages.
GhCrRLK1L104 was highly expressed in fibers at 30 DPA and overexpression of
GhCrRLK1L104 in
Arabidopsis increased the trichomes length which indicate its function in cell elongation [
16].
Recently, fiber osmoregulation has emerged as a crucial player in regulating cotton fiber initiation and elongation. A two-year field experiment was conducted to determine whether potassium ameliorates
Gossypium hirsutum fiber length by regulating osmotic and K
+/Na
+ homeostasis under salt stress [
17]. Another study showed that low soil available phosphorus (AP) contents (P
0: 3 ±0.5; P
1: 6 ± 0.5 mg kg
-1) inhibited the fiber cell elongation leading to reduce the maximum velocity of fiber elongation (V
Lmax) and fiber length, mainly due to lower malate content and V-H
+-PPase activities [
18]. A cell wall–localized β-1,3-glucanase, GhGLU18 was found to promote fiber elongation and cell wall thickening by degrading callose and enhancing polysaccharide metabolism [
19]. Calcium is also involved in fiber development. In vitro ovule culture demonstrated Ca
2+ rescued the shorter-fiber phenotype of
GhIQD10 overexpression lines.
GhIQD10 was expressed mainly in the transition period of cotton fiber development.
GhIQD10 interacted with
GhCaM7 and the interaction was inhibited by Ca
2+ [
20].
Phytohormones, including gibberellin acid (GA), auxin, cytokinin, brassinosteroid (BR), abscisic acid (ABA), ethylene, jasmonic acid (JA), cytokinin, salicylic acid, and strigolactone (SL), are small endogenous signaling molecules in plants [
21]. Many of these hormones directly participate in fiber initiation and elongation. For example, strigolactones are a class of carotenoid-derived plant hormones that modulate cotton fiber elongation and secondary cell wall thickening. The endogenous SLs were significantly higher in fibers 20 DPA. Exogenous SLs significantly increased fiber length and cell wall thickness [
22]. Fiber-specific expression of
GhOR1Del, a positive regulator of carotenoid accumulation, was found to upregulate the carotenoid level in cotton fiber simultaneously increasing the contents of carotenoids, ABA, and ethylene in elongating fibers [
23].
This manuscript provides a comprehensive review of the essential roles played by various transcription factors, such as MYB, WRKY, HD-ZIP and bHLH transcription factors, and phytohormones including auxin, GA, BR, JA, ethylene, ABA and cytokinin during fiber initiation and elongation. The review offers an in-depth analysis of their contributions to cotton fiber development, emphasizing their interplay and regulatory mechanisms. Understanding these factors holds immense potential for advancing our knowledge and optimizing cotton cultivation for improved fiber quality and yield.
2. Fiberless Mutants Identification
Cotton fibers originate from single cells within embryo epidermal cells, with only 25~30% of them ultimately developing into fibers [
5]. The initiation and final development of fibers play a crucial role in determining cotton yield. Stewart (1975) observed a fiber density of approximately 3,300 fibers per mm
2, with the ratio of fiber initials to total epidermal cells being 1:3.7 at anthesis [
24]. Two types of fibers exist: lint fibers and fuzz fibers, with lint fibers possessing higher economic value. Interestingly, there is no observable phenotype difference when both types of fibers initiate on the epidermal surface [
25]. Lint fibers initiate on or a day before the day of anthesis (-1 ~ 0 DPA) and elongate to 2~3.5 cm, while fuzz fibers initiate 3~5 DPA and only reach around 5~10 mm in length [
24].
Several fiberless or fuzzless mutants, including four dominant (
Li1,
Li2,
N1, and
Fbl) and three recessive (
n2,
sma-4 (ha), and
sma-4 (fz)) mutants, have been studied to identify genes and understand their interactions in the molecular mechanism of fiber initiation and elongation [
26]. By crossing fuzzless and/or lintless mutants, the
N1,
n2,
Li3, and
Li4 loci were identified to control the presence or absence of lint or fuzz. For instance,
N1N1 confers the presence of fuzz,
n2n2 inhibits fuzz initiation and development, and duplicate gene pairs
Li3Li3 and
Li4Li4 determine the presence of lint. Homozygosity for
li3li3 and
li4li4 may also inhibit fuzz development [
27]. Comparative scanning electron microscopy studies of fiber development in a normal TM-1 genotype and the near-isogenic
Li1 mutant at 1 DAP and 3 DAP revealed minimal differences during early stages, suggesting that
Li1 gene expression occurs later, probably during the elongation phase [
28]. Cross-pollination of
N1,
n2, and
n3 in upland cotton lines produced fiberless seeds, such as
MD17 and
SL1-7-1 [
29]. The mutant exhibited lower short fiber content and better yarn quality than the wild type cultivar [
30]. The
Li2 short fiber mutation is located within a terminal deletion of chromosome 18 in cotton [
31]. siRNA-induced silencing of a family of
RanBP1s inhibits the elongation of cotton fiber cells in the
Li2 mutant [
32].
XZ142FLM, a natural fiberless mutant with well-studied genes responsible for the fiberless phenotype, stands in contrast to
GhVIN1-RNAi (
GhVIN1i), one of the few fiberless cotton lines associated with sugar metabolism and signaling generated through reverse genetics. Comparative transcriptome analysis between the natural fiberless mutant
XZ142FLM and the transgenic fiberless line
GhVIN1i, obtained by RNAi silencing of
GhVIN1, identified common differentially expressed genes (DEGs) in ovules during fiber initiation. The respective DEGs were enriched in several identical pathways related to fiber initiation, revealing shared molecular regulatory networks controlling fiber initiation [
33]. Several genes related to fiber initiation or elongation, such as
GhMML3,
GhVIN1,
GhMYB25,
GhHD-1, and
GhHOX3, exhibited similar expression patterns in mutant
XZ142FLM and the RNAi mutant
GhVIN1i during fiber early development. This suggests the operation of similar mechanisms for fiber initiation in these two fiberless lines. The study sheds light on the regulatory networks mediated by
GhMML3 and
GhVIN1 in controlling fiber initiation in cotton [
33].
In a study investigating natural antisense transcripts and siRNA control over fiber development, researchers discovered that small RNA derived from the
GhMML3_A12 locus can induce self-cleavage of
GhMML3_A12 mRNA, leading to the production of naked seeds and subsequent inhibition of lint fiber in
N1 plants [
34]. Employing a map-based cloning strategy for the first time in tetraploid cotton, they successfully cloned the naked seed mutant gene (
N1) encoding a
MYBMIXTA-like transcription factor 3 (
MML3)/
GhMYB25-like on chromosome A12, known as
GhMML3_A12, associated with fuzz fiber development [
34]. Phenotypic and genotypic analysis of MYB25-like alleles in cottons exhibiting various fiber phenotypes and their crossed progeny revealed that both
MYB25-like_At and
MYB25-like_Dt are linked to lint development. Fuzz development, on the other hand, is primarily determined by the expression level of
MYB25-like_Dt at approximately 3 DPA, making
MYB25-like_Dt a strong candidate for
N2. Recently, cotton microtubule-associated protein
GhMAP20L5 was reported to mediate fiber elongation through the interaction with the tubulin
GhTUB13. In the RNA silencing plants,
GhMAP20L5 expression was repressed at least 28.1% in comparison with its null plants at different development stages (0 DPA, 6 DPA, 12 DPA, and 18 DPA) which resulted in reducing fiber elongation rate, fiber length and lint percentage [
35].