Genome-wide analysis and expression profiling of the phospholipase D gene family in Gossypium arboreum

Genome-wide characterization of phospholipase D family genes in allotetraploid peanut and its diploid progenitors revealed their crucial roles in growth and abiotic stress responses

Abiotic stresses such as cold, drought and salinity are the key environmental factors that limit the yield and quality of oil crop peanut. Phospholipase Ds (PLDs) are crucial hydrolyzing enzymes involved in lipid mediated signaling and have valuable functions in plant growth, development and stress tolerance. Here, 22, 22 and 46 PLD genes were identified in Arachis duranensis, Arachis ipaensis and Arachis hypogaea, respectively, and divided into α, β, γ, δ, ε, ζ and φ isoforms. Phylogenetic relationships, structural domains and molecular evolution proved the conservation of PLDs between allotetraploid peanut and its diploid progenitors. Almost each A. hypogaea PLD except for AhPLDα6B had a corresponding homolog in A. duranensis and A. ipaensis genomes. The expansion of Arachis PLD gene families were mainly attributed to segmental and tandem duplications under strong purifying selection. Functionally, the most proteins interacting with AhPLDs were crucial components of lipid metabolic pathways, in which ahy-miR3510, ahy-miR3513-3p and ahy-miR3516 might be hub regulators. Furthermore, plenty of cis-regulatory elements involved in plant growth and development, hormones and stress responses were identified. The tissue-specific transcription profiling revealed the broad and unique expression patterns of AhPLDs in various developmental stages. The qRT-PCR analysis indicated that most AhPLDs could be induced by specific or multiple abiotic stresses. Especially, AhPLDα3A, AhPLDα5A, AhPLDβ1A, AhPLDβ2A and AhPLDδ4A were highly up-regulated under all three abiotic stresses, whereas AhPLDα9A was neither expressed in 22 peanut tissues nor induced by any abiotic stresses. This genome-wide study provides a systematic analysis of the Arachis PLD gene families and valuable information for further functional study of candidate AhPLDs in peanut growth and abiotic stress responses.

Characterization and evolutionary diversification of the phospholipase D gene family in mosses

Plant phospholipase D (PLD) exerts important roles in various biological processes, such as intracellular signaling and morphological development. Our knowledge about early land plant PLDs is still underdeveloped. In this study, we identified 84 PLD genes in six mosses, i.e., Physcomitrella patens, Ceratodon purpureus, Fontinalis antipyretica, Pleurozium schreberi, Sphagnum magellanicum, and Sphagnum fallax. These PLDs were classified into four clades (I–IV). We showed that PLD underwent rapid expansion in mosses. A total of six conserved domains and two core HKD motifs were detected. Structure analysis uncovered that the moss PLDs from within a clade generally exhibited similar exon-intron organization. Cis-elements prediction and expression analyses indicated that P. patens PLDs had key roles in stress responsiveness and plant development. Particularly, about half of the P. patens PLDs (e.g., PpPLD1, PpPLD2, and PpPLD5) were differentially expressed under biotic and abiotic stresses. We also determined the expression pattern of P. patens PLD genes in various tissues and at different stages of development. Although the moss, clubmoss, liverwort, and fern PLDs evolved largely under functional constraints, we found episodic positive selection in the moss PLDs, e.g., C. purpureus PLD2 and P. patens PLD11. We infer that the evolutionary force acting on the PLDs may have facilitated moss colonization of land. Our work provides valuable insights into the diversification of moss PLD genes, and can be used for future studies of their functions.

Genome-wide in silico identification of phospholipase D (PLD) gene family from Corchorus capsularis and Corchorus olitorius: reveals their responses to plant stress

Background

Plant grows in nature facing various types of abiotic stresses for their normal growth and development. During abiotic stress, plants evolve different types of mechanisms to survive in a hostile environment. Phospholipase D (PLD) plays important role in the regulation of diverse cellular processes including stress responses in plants. Member of PLD genes are well studied in different model plants; however, their functions in the jute are not clear yet.

Result

In the present study, a total of 12 and 11 PLD genes were identified in the genome of C. capsularis and C. olitorius, respectively. The presence of the two conserved HKD motifs in PLD genes except for CoPLDδ-2 in jute suggests their strong lipase activity. Twenty different motifs were found in the identified PLD genes, and PLD-β1, PLD-γ1, and all members of PLD-δ1 of both jute species contained the highest number of motifs. Phylogenetic analysis showed the close evolutionary relationship among the five groups of jute PLD proteins along with the PLD proteins from Arabidopsis. Tissue-specific expression pattern of PLDα1-2, PLD-α2, PLDβ1, PLDγ1, and PLDδ1 of two jute species suggested their involvement in plant growth and development. However, the expression pattern of PLDα1-2, PLDα1-3, PLD-α4, PLDδ1, and PLDδ3 indicated their association during waterlogging stress. In addition, PLD-α2, PLDβ1, and PLDδ2 seemed to be involved in drought stress as well as salinity stress.

Conclusion

This genome-wide identification of jute PLD genes from C. capsularis and C. olitorius will help to further functional characterization of the PLD genes for developing stress-tolerant jute variety.

Genome-Wide Analysis of the Phospholipase D Family in Five Cotton Species, and Potential Role of GhPLD2 in Fiber Development and Anther Dehiscence

Phospholipase D (PLD) and its hydrolysis product phosphatidic acid play an important role in the regulation of several cellular processes, including root growth, pollen tube elongation, and microtubule reorganization. Here, we systematically identified and analyzed the membership, characterization, and evolutionary relationship of PLDs in five species of cotton. The results of the transcriptomic analysis suggested that the evaluated PLD genes showed high expression levels in anther tissue and during the fiber initiation and elongation periods. Quantitative real-time polymerase chain reaction showed differential expression of GhPLD genes in the anthers of photoperiod sensitive male sterility mutant 5 (psm5). Previous research on multiple stable quantitative trait loci also suggests the role of PLD genes in the fiber development. Further analyses showed that GhPLD2 protein is localized to the plasma membrane. The virus-induced gene silencing of GhPLD2 in cotton seedlings repressed its expression by 40-70%, which led to a reduction in reactive oxygen species (ROS) levels, 22% anther indehiscence, and disrupted fiber initiation and elongation. Thus, we inferred that GhPLD2 may promote ROS production, which, in turn, may regulate anther dehiscence and fiber development.

Cotton GhSSI2 isoforms from the stearoyl acyl carrier protein fatty acid desaturase family regulate Verticillium wilt resistance

Lipids are major and essential constituents of plant cells and provide energy for various metabolic processes. However, the function of the lipid signal in defence against Verticillium dahliae, a hemibiotrophic pathogen, remains unknown. Here, we characterized 19 conserved stearoyl-ACP desaturase family proteins from upland cotton (Gossypium hirsutum). We further confirmed that GhSSI2 isoforms, including GhSSI2-A, GhSSI2-B, and GhSSI2-C located on chromosomes A10, D10, and A12, respectively, played a dominant role to the cotton 18:1 (oleic acid) pool. Suppressing the expression of GhSSI2s reduced the 18:1 level, which autoactivated the hypersensitive response (HR) and enhanced cotton Verticillium wilt and Fusarium wilt resistance. We found that low 18:1 levels induced phenylalanine ammonia-lyase-mediated salicylic acid (SA) accumulation and activated a SA-independent defence response in GhSSI2s-silenced cotton, whereas suppressing expression of GhSSI2s affected PDF1.2-dependent jasmonic acid (JA) perception but not the biosynthesis and signalling cascade of JA. Further investigation showed that structurally divergent resistance-related genes and nitric oxide (NO) signal were activated in GhSSI2s-silenced cotton. Taken together, these results indicate that SA-independent defence response, multiple resistance-related proteins, and elevated NO level play an important role in GhSSI2s-regulated Verticillium wilt resistance. These findings broaden our knowledge regarding the lipid signal in disease resistance and provide novel insights into the molecular mechanism of cotton fungal disease resistance.

Genome-Wide Analysis and Expression Profiling of the Phospholipase D Gene Family in Solanum tuberosum

Phospholipase D (PLD) is the most important phospholipid hydrolase in plants, which can hydrolyze phospholipids into phosphatidic acid (PA) and choline. When plants encounter low temperature, drought and high salt stress, phospholipase D and its products play an important role in regulating plant growth and development and coping with stress. In this study, 16 members of StPLD gene family were identified in potato genome, which were distributed in α, β, δ, and ζ subfamilies, and their expression patterns under salt, high temperature, drought, and ABA stress were detected by qRT-PCR method. Gene expression analysis showed that the expression of StPLD genes in potato was upregulated and downregulated to varying degrees under the four stresses, indicating that the PLD gene family is involved in the interaction of potato plant hormones and abiotic stress signals. Chromosome distribution showed that StPLD gene was unevenly distributed on 8 chromosomes, and only one pair of tandem repeat genes was found. All StPLD promoters contain hormone and stress-related cis-regulatory elements to respond to different stresses. Structural analysis showed that StPLD genes in the same subgroup had a similar exon-intron structure. Our study provides a valuable reference for further research of the function and structure of PLD gene.

Identification, evolution, expression analysis of phospholipase D (PLD) gene family in tea (Camellia sinensis)

Phospholipase D (PLD) (EC 3.1.4.4) plays important roles in plants growth, development, and response to environmental stresses. Tea plant (Camellia sinensis) is the most important non-alcoholic beverage in the world with health benefits, but tea production decreases in response to environmental stresses such as cold and drought. Therefore, a genome-wide analysis of the C. sinensis PLD gene family (CsPLDs) was carried out. In the current study, identification, evolutionary relationship, duplication, selection pressure, gene structure, promoter analysis, transcript-targeted miRNA, and simple sequence repeat markers prediction, RNA-seq data analysis, and three-dimensional structure of the CsPLDs have been investigated using bioinformatics tools. 15 PLDs were identified from the tea genome which belongs to five groups, including CsPLDα, CsPLDβ, CsPLDδ, CsPLDε, and CsPLDζ. Both segmental and tandem duplications have occurred in the CsPLD gene family. Ka/Ks ratio for the most duplicated pair genes was less than 1 which implies negative selection to conserve their function during the tea evolution. 68 cis-elements have been found in CsPLDs indicating the contribution of these genes in response to environmental stresses. Likewise, 72 SSR loci and 96 miRNA molecules in 14 and 15 CsPLDs have been detected. According to RNA-seq data, the highest expression in all tissues under various abiotic stresses was related to CsPLDα1. Besides, a three-dimensional structure of the CsPLDα1 was evaluated to better understand its biological activity. This research provides comprehensive information that could be useful in future studies to develop stress-tolerant tea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01007-0.

Characterization of a PLDζ2 Homology Gene from Developing Castor Bean Endosperm

Castor oil contains approximately 90% ricinoleic acid (RA) which is stored mainly in the form of tri-ricinoleic acid containing triacylglycerols (TAG). Ricinoleate is synthesized from oleate (18:1n-9) esterified to the sn-2 position of phosphatidylcholine (PtdCho) catalyzed by oleoyl-12-hydroxylase. PtdCho-derived diacylglycerol (DAG) is an important substrate pool for TAG synthesis, and the interconversion between PtdCho and DAG has been shown to play a critical role in channeling hydroxy fatty acids (HFA) to TAG. Although phospholipase D (PLD) has been reported to catalyze the hydrolysis of PtdCho to produce phosphatidic acid which can then be converted to DAG, its potential functions in the channeling of RA from PtdCho to DAG and the assembly of RA on TAG is largely unknown. In the present study, 11 PLD genes were identified from the Castor Bean Genome Database. Gene expression analysis indicated that RcPLD9 is expressed at relatively high levels in developing seeds compared to other plant tissues. Sequence and phylogenetic analyses revealed that RcPLD9 is a homolog of Arabidopsis PLDζ2. Overexpression of RcPLD9 in the Arabidopsis CL7 line producing C18-HFA resulted in RA content reductions in the polar lipid fraction (mainly PtdCho) and mono-HFA-TAG, but increased RA content in di-HFA-TAG. Since part of RA in di-HFA-TAG is derived from HFA-DAG, the results indicated that RcPLD9 facilitates the channeling of RA from PtdCho to DAG for its assembly on TAG in developing seeds.

Genome-Wide Analysis of Phospholipase D Gene Family and Profiling of Phospholipids under Abiotic Stresses in Brassica napus

Oil crop Brassica napus is subjected to environmental stresses such as drought, cold and salt. Phospholipase Ds (PLDs) have vital roles in regulation of plant growth, development and stress tolerance. In this study, 32 BnaPLD genes were identified and classified into six subgroups depending on the conserved protein structures. High similarity in gene and protein structures exists between BnaPLDs and AtPLDs. Gene expression analysis showed that BnaPLDα1s and BnaPLDδs had higher expression than other PLDs. BnaPLDα1 and BnaPLDδ were significantly induced by abiotic stresses including dehydration, NaCl, abscisic acid (ABA) and 4�C. Lipidomic analysis showed that the content of main membrane phospholipids decreased gradually under stresses, except phosphatidylglycerol increased under the treatment of ABA and phosphatidylethanolamine increased under 4�C. Correspondingly, their product of phosphatidic acid increased often with a transient peak at 8 h. The plant height of mutants of PLDα1 was significantly reduced. Agronomic traits such as yield, seed number, silique number and branches were significantly impaired in PLDα1 mutants. These results indicate that there is a large family of PLD genes in B. napus, especially BnaPLDα1s and BnaPLDδs may play important roles in membrane lipids remodeling and maintaining of the growth and stress tolerance of B. napus.

Phospholipase D and phosphatidic acid in plant immunity

Phospholipase D (PLD) hydrolyzes membrane phospholipids to generate phosphatidic acid (PA). Both PLD and its lipid product PA are involved in various physiological processes, including plant response to pathogens. The PLD family is comprised of multiple members in higher plants, and PLDs have been reported to play positive and/or negative roles in plant immunity, depending on the types of pathogens and specific PLDs involved. Individual PLDs have distinguishable biochemical properties, such as Ca²⁺ and phosphatidylinositide requirements. In addition, PLDs and PA are found to interact with various proteins in hormone and stress signaling. The different biochemical and regulatory properties of PLDs and PA shed light on the mechanisms for the functional diversity of PLDs in plant defense signaling and response.

Genome-wide analysis of BES1 genes in Gossypium revealed their evolutionary conserved roles in brassinosteroid signaling

Brassinosteroids (BRs), which are essential phytohormones for plant growth and development, are important for cotton fiber development. Additionally, BES1 transcription factors are critical for BR signal transduction. However, cotton BES1 family genes have not been comprehensively characterized. In this study, we identified 11 BES1 genes in G. arboreum, 11 in G. raimondii, 16 in G. barbadense, and 22 in G. hirsutum. The BES1 sequences were significantly conserved in the Arabidopsis thaliana, rice, and upland cotton genomes. A total of 94 BES1 genes from 10 different plant species were divided into three clades according to the neighbor-joining and minimum-evolution methods. Moreover, the exon/intron patterns and motif distributions were highly conserved among the A. thaliana and cotton BES1 genes. The collinearity among the orthologs from the At and Dt subgenomes was estimated. Segmental duplications in the At and Dt subgenomes were primarily responsible for the expansion of the cotton BES1 gene family. Of the GhBES1 genes, GhBES1.4_At/Dt exhibited BL-induced expression and was predominantly expressed in fibers. Furthermore, Col-0/mGhBES1.4_At plants produced curled leaves with long and bent petioles. These transgenic plants also exhibited decreased hypocotyl sensitivity to brassinazole and constitutive BR induced/repressed gene expression patterns. The constitutive BR responses of the plants overexpressing mGhBES1.4_At were similar to those of the bes1-D mutant.

Cotton functional genomics reveals global insight into genome evolution and fiber development

Due to the economic value of natural textile fiber, cotton has attracted much research attention, which has led to the publication of two diploid genomes and two tetraploid genomes. These big data facilitate functional genomic study in cotton, and allow researchers to investigate cotton genome structure, gene expression, and protein function on the global scale using high-throughput methods. In this review, we summarized recent studies of cotton genomes. Population genomic analyses revealed the domestication history of cultivated upland cotton and the roles of transposable elements in cotton genome evolution. Alternative splicing of cotton transcriptomes was evaluated genome-widely. Several important gene families like MYC, NAC, Sus and GhPLDα1 were systematically identified and classified based on genetic structure and biological function. High-throughput proteomics also unraveled the key functional proteins correlated with fiber development. Functional genomic studies have provided unprecedented insights into global-scale methods for cotton research.

'Bending' models of halotropism: incorporating protein phosphatase 2A, ABCB transporters, and auxin metabolism

Salt stress causes worldwide reductions in agricultural yields, a problem that is exacerbated by the depletion of global freshwater reserves and the use of contaminated or recycled water (i.e. effluent water). Additionally, salt stress can occur as cultivated areas are subjected to frequent rounds of irrigation followed by periods of moderate to severe evapotranspiration, which can result in the heterogeneous aggregation of salts in agricultural soils. Our understanding of the later stages of salt stress and the mechanisms by which salt is transported out of cells and roots has greatly improved over the last decade. The precise mechanisms by which plant roots perceive salt stress and translate this perception into adaptive, directional growth away from increased salt concentrations (i.e. halotropism), however, are not well understood. Here, we provide a review of the current knowledge surrounding the early responses to salt stress and the initiation of halotropism, including lipid signaling, protein phosphorylation cascades, and changes in auxin metabolism and/or transport. Current models of halotropism have focused on the role of PIN2- and PIN1-mediated auxin efflux in initiating and controlling halotropism. Recent studies, however, suggest that additional factors such as ABCB transporters, protein phosphatase 2A activity, and auxin metabolism should be included in the model of halotropic growth.

Genome-wide identification and characterization of phospholipase C gene family in cotton (Gossypium spp.)

Phospholipase C (PLC) are important regulatory enzymes involved in several lipid and Ca²⁺-dependent signaling pathways. Previous studies have elucidated the versatile roles of PLC genes in growth, development and stress responses of many plants, however, the systematic analyses of PLC genes in the important fiber-producing plant, cotton, are still deficient. In this study, through genome-wide survey, we identified twelve phosphatidylinositol-specific PLC (PI-PLC) and nine non-specific PLC (NPC) genes in the allotetraploid upland cotton Gossypium hirsutum and nine PI-PLC and six NPC genes in two diploid cotton G. arboretum and G.raimondii, respectively. The PI-PLC and NPC genes of G. hirsutum showed close phylogenetic relationship with their homologous genes in the diploid cottons and Arabidopsis. Segmental and tandem duplication contributed greatly to the formation of the gene family. Expression profiling indicated that few of the PLC genes are constitutely expressed, whereas most of the PLC genes are preferentially expressed in specific tissues and abiotic stress conditions. Promoter analyses further implied that the expression of these PLC genes might be regulated by MYB transcription factors and different phytohormones. These results not only suggest an important role of phospholipase C members in cotton plant development and abiotic stress response but also provide good candidate targets for future molecular breeding of superior cotton cultivars.

Molecular characterization of GhPLDα1 and its relationship with secondary cell wall thickening in cotton fibers

Phospholipase D (PLD) hydrolyzes phospholipids to generate a free polar head group (e.g., choline) and a second messenger phosphatidic acid and plays diverse roles in plant growth and development, including seed germination, leaf senescence, root hair growth, and hypocotyl elongation. However, the function of PLD in cotton remains largely unexplored. Here, the comprehensive molecular characterization of GhPLDα1 was explored with its role in upland cotton (Gossypium hirsutum) fiber development. The GhPLDα1 gene was cloned successfully, and a sequence alignment showed that GhPLDα1 contains one C2 domain and two HKD (HxKxxxxD) domains. Quantitative reverse transcriptase-polymerase chain reaction measured the expression of GhPLDα1 in various cotton tissues with the highest level in fibers at 20 days post anthesis (d.p.a.). Fluorescent microscopy and immunoblotting in tobacco epidermis showed the GhPLDα1 distribution in both cell membranes and the cytoplasm. An activity assay indicated changes in PLDα enzyme activity in developing fiber cells with a peak level at 20 d.p.a., coinciding with the onset of cellulose accumulation and the increased H₂O₂ content during fiber development. Furthermore, the inhibition of PLDα activity obviously decreased the cellulose and H₂O₂ contents of in vitro-cultured cotton fibers. These results provide important evidence explaining the relationship of GhPLDα1 with secondary cell wall thickening in cotton fibers in that GhPLDα1 may correlate with the increased H₂O₂ content at the onset of secondary cell wall thickening, ultimately promoting cellulose biosynthesis.

Phospholipases Dα and δ are involved in local and systemic wound responses of cotton (G. hirsutum)

Phospholipases D (PLDs) catabolize structural phospholipids to produce phosphatidic acid (PtdOH), a lipid playing central role in signalling pathways in animal, yeast and plant cells. In animal cells two PLD genes have been studied while in model plant Arabidopsis twelve genes exist, classified in six classes (α-ζ). This underlines the role of these enzymes in plant responses to environmental stresses. However, information concerning the PLD involvement in the widely cultivated and economically important cotton plant responses is very limited. The aim of this report was to study the activity of conventional cotton PLD and its participation in plant responses to mechanical wounding, which resembles both biotic and abiotic stresses. PLDα activity was identified and further characterized by transphosphatidylation reaction. Upon wounding, cotton leaf responses consist of an acute in vitro increase of PLDα activity in both wounded and systemic tissue. However, determination of the in vivo PtdOH levels under the same wounding conditions revealed a rapid PtdOH formation only in wounded leaves and a late response of a PtdOH increase in both tissues. Εxpression analysis of PLDα and PLDδ isoforms showed mRNA accumulation of both isoforms in the wounded tissue, but only PLDδ exerts a high and sustainable expression in systemic leaves, indicating that this isoform is mainly responsible for the systemic wound-induced PtdOH production. Therefore, our data suggest that PLDα and PLDδ isoforms are involved in different steps in cotton wound signalling.

•

PLDα activity and PtdOH levels rapidly increase in wounded cotton leaves.

•

PLDα is also activated rapidly in systemic tissue.

•

Doubling of PtdOH levels occurs as a late response in both wounded and systemic tissue.

•

PLDδ (but not PLDα) exerts a high and sustainable expression in systemic leaves.

•

PLDα and PLDδ are involved in different steps in cotton wound signalling.

Genome-Wide Comparative Analysis of the Phospholipase D Gene Families among Allotetraploid Cotton and Its Diploid Progenitors

In this study, 40 phospholipase D (PLD) genes were identified from allotetraploid cotton Gossypium hirsutum, and 20 PLD genes were examined in diploid cotton Gossypium raimondii. Combining with 19 previously identified Gossypium arboreum PLD genes, a comparative analysis was performed among the PLD gene families among allotetraploid and two diploid cottons. Based on the orthologous relationships, we found that almost each G. hirsutum PLD had a corresponding homolog in the G. arboreum and G. raimondii genomes, except for GhPLDβ3A, whose homolog GaPLDβ3 may have been lost during the evolution of G. arboreum after the interspecific hybridization. Phylogenetic analysis showed that all of the cotton PLDs were unevenly classified into six numbered subgroups: α, β/γ, δ, ε, ζ and φ. An N-terminal C2 domain was found in the α, β/γ, δ and ε subgroups, while phox homology (PX) and pleckstrin homology (PH) domains were identified in the ζ subgroup. The subgroup φ possessed a single peptide instead of a functional domain. In each phylogenetic subgroup, the PLDs showed high conservation in gene structure and amino acid sequences in functional domains. The expansion of GhPLD and GrPLD gene families were mainly attributed to segmental duplication and partly attributed to tandem duplication. Furthermore, purifying selection played a critical role in the evolution of PLD genes in cotton. Quantitative RT-PCR documented that allotetraploid cotton PLD genes were broadly expressed and each had a unique spatial and developmental expression pattern, indicating their functional diversification in cotton growth and development. Further analysis of cis-regulatory elements elucidated transcriptional regulations and potential functions. Our comparative analysis provided valuable information for understanding the putative functions of the PLD genes in cotton fiber.

Identification and Profiling of microRNAs Expressed in Elongating Cotton Fibers Using Small RNA Deep Sequencing

Plant microRNAs (miRNAs) have been shown to play essential roles in the regulation of gene expression. In this study, small RNA deep sequencing was applied to explore novel miRNAs expressed in elongating cotton fibers. A total of 46 novel and 96 known miRNAs, primarily derived from the corresponding specific loci in genome of Gossypium arboreum, were identified. 64 miRNAs were shown to be differentially expressed during the fiber elongation process; 16 were predicted to be novel miRNAs while the remaining 48 belong to known miRNA families. Furthermore, RLM-5' RACE (RNA ligase-mediated rapid amplification of 5'-cDNA ends) experiments identified the targets of eight important miRNAs, and the expression levels of these target genes were confirmed to be negatively correlated with the expression patterns of their corresponding miRNAs. We propose a potential functional network mediated through these eight miRNAs to illustrate their important functions in fiber elongation. Our study provides novel insights into the dynamic profiles of these miRNAs and a basis for investigating the regulatory mechanisms involved in the elongation of cotton fibers.