Abstract:Background: Epidermal patterning factor / -like (EPF/EPFL) gene family encodes a class of cysteine-rich secretory peptides, which are widelyfound in terrestrial plants.Multiple studies has indicated that EPF/EPFLs might play significant roles in coordinating plant development and growth, especially as the morphogenesis processes of stoma, awn, stamen, and fruit skin. However, few research on EPF/EPFL gene family was reported in Gossypium. Results: We separately identified 20 G. raimondii, 24 G. arboreum, 44 G. hirsutum, and 44 G. barbadense EPF/EPFL genes in the 4 representative cotton species, which were divided into four clades together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 17 Selaginella moellendorffii ones based on their evolutionary relationships. The similar gene structure and common motifs indicated the high conservation among the EPF/EPFL members, while the uneven distribution in chromosomes implied the variability during the long-term evolutionary process. Hundreds of collinearity relationships were identified from the pairwise comparisons of intraspecifc and interspecific genomes, which illustrated gene duplication might contribute to the expansion of cotton EPF/EPFL gene family. A total of 15 kinds of cis-regulatory elements were predicted in the promoter regions, and divided into three major categories relevant to the biological processes of development and growth, plant hormone response, and abiotic stress response. Having performing the expression pattern analyses with the basic of the published RNA-seq data, we found most of GhEPF/EPFL and GbEPF/EPFL genes presented the relatively low expression levels among the 9 tissues or organs, while showed more dramatically different responses to high/low temperature and salt or drought stresses. Combined with transcriptome data of developing ovules and fibers and quantitative Real-time PCR results (qRT-PCR) of 15 highly expressed GhEPF/EPFL genes, it could be deduced that the cotton EPF/EPFL genes were closely related with fiber development. Additionally, the networks of protein-protein interacting among EPF/EPFLs concentrated on the cores of GhEPF1 and GhEPF7, and thosefunctional enrichment analyses indicated that most of EPF/EPFLs participate in the GO (Gene Ontology) terms of stomatal development and plant epidermis development, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways of DNA or base excision repair. Conclusion: Totally, 132 EPF/EPFL genes were identified for the first time in cotton, whose bioinformatic analyses of cis-regulatory elements and expression patterns combined with qRT-PCR experiments to prove the potential functions in the biological processes of plant growth and responding to abiotic stresses, specifically in the fiber development. These results not only provide comprehensive and valuable information for cotton EPF/EPFL gene family, but also lay solid foundation for screening candidate EPF/EPFL genes in further cotton breeding.

Membrane Interactions of GET1 and GET2 Facilitate Fiber Cell Initiation through the Guided Entry of the TA Protein Pathway in Cotton

EB1C forms dimer and interacts with protein phosphatase 2A (PP2A) to regulate fiber elongation in upland cotton ( Gossypium hirsutum )

Characterization of chromatin accessibility and gene expression reveal the key genes involved in cotton fiber elongation

Metabolomic and Transcriptomic Insights into How Cotton Fiber Transitions to Secondary Wall Synthesis, Represses Lignification, and Prolongs Elongation

Functional Analysis of the Cotton CLE Polypeptide Signaling Gene Family in Plant Growth and Development

Cotton Fiber Elongation Network Revealed by Expression Profiling of Longer Fiber Lines Introgressed with Different Gossypium Barbadense Chromosome Segments

GhWER controls fiber initiation and early elongation by regulating ethylene signaling pathway in cotton (Gossypium hirsutum)

Cotton microtubule-associated protein GhMAP20L5 mediates fiber elongation through the interaction with the tubulin GhTUB13

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

Co-expression Networks Regulating Cotton Fiber Initiation Generated by Comparative Transcriptome Analysis Between Fiberless XZ142FLM and GhVIN1i

Association of GhGeBP genes with fiber quality and early maturity related traits in upland cotton

A Novel Tandem Zinc Finger Protein in Gossypium hirsutum, GhTZF2, Interacts with GhMORF8 to Regulate Cotton Fiber Cell Development

A Genome-Wide Analysis of the CEP Gene Family in Cotton and a Functional Study of GhCEP46-D05 in Plant Development

Genome-wide analyses of member identification, expression pattern, and protein-protein interaction of EPF/EPFL gene family in Gossypium

Genetics And Evolution Of Mixta Genes Regulating Cotton Lint Fiber Development

GhGRF4/GhARF2‐GhGASA24 module regulates fiber cell wall thickness by modulating cellulose biosynthesis in upland cotton (Gossypium hirsutum)

The Mir319-Targeted GhTCP4 Promotes the Transition from Cell Elongation to Wall Thickening in Cotton Fiber.

Transcriptional regulation of phospholipid transport in cotton fiber elongation by GhMYB30D04–GhHD1 interaction complex

Co-Expression Network Analysis and Introgressive Gene Identification for Fiber Length and Strength Reveal Transcriptional Differences in 15 Cotton Chromosome Substitution Segment Lines and Their Upland and Sea Island Parents

R2R3 MYB Transcription Factor GhMYB201 Promotes Cotton Fiber Elongation via Cell Wall Loosening and Very-Long-Chain Fatty Acid Synthesis

The R3-Myb Gene Ghcpc Negatively Regulates Cotton Fiber Elongation