APEX2-Mediated Proximity Protein Labeling in Dictyostelium.

Largely due to its simplicity, while being more like human cells compared to other experimental models, Dictyostelium continues to be of great use to discover basic molecular mechanisms and signaling pathways underlying evolutionarily conserved biological processes. However, the identification of new protein interactions implicated in signaling pathways can be particularly challenging in Dictyostelium due to its extremely fast signaling kinetics coupled with the dynamic nature of signaling protein interactions. Recently, the proximity labeling method using engineered ascorbic acid peroxidase 2 (APEX2) in mammalian cells was shown to allow the detection of weak and/or transient protein interactions and also to obtain spatial and temporal resolution. Here, we describe a protocol for successfully using the APEX2-proximity labeling method in Dictyostelium. Coupled with the identification of the labeled proteins by mass spectrometry, this method expands Dictyostelium's proteomics toolbox and should be widely useful for identifying interacting partners involved in a variety of biological processes in Dictyostelium.

14-3-3 proteins as a major hub for plant immunity.

14-3-3 proteins, ubiquitously present in eukaryotic cells, are regulatory proteins involved in a plethora of cellular processes. In plants, they have been studied in the context of metabolism, development, and stress responses. Recent studies have highlighted the pivotal role of 14-3-3 proteins in regulating plant immunity. The ability of 14-3-3 proteins to modulate immune responses is primarily attributed to their function as interaction hubs, mediating protein-protein interactions and thereby regulating the activity and overall function of their binding partners. Here, we shed light on how 14-3-3 proteins contribute to plant defense mechanisms, the implications of their interactions with components of plant immunity cascades, and the potential for leveraging this knowledge for crop improvement strategies.

Metabolomic Approach to Identify the Potential Metabolites from Alpinia malaccensis for Treating SARS-CoV-2 Infection.

The advent of the new coronavirus, leading to the SARS-CoV-2 pandemic, has presented a substantial worldwide health hazard since its inception in the latter part of 2019. The severity of the current pandemic is exacerbated by the occurrence of re-infection or co-infection with SARS-CoV-2. Hence, comprehending the molecular process underlying the pathophysiology of sepsis and discerning possible molecular targets for therapeutic intervention holds significant importance. For the first time, 31 metabolites were tentatively identified by GC-MS analysis from Alpinia malaccensis. On the other hand, five phenolic compounds were identified and quantified from the plant in HPLC-DAD analysis, including (-) epicatechin, rutin hydrate, rosmarinic acid, quercetin, and kaempferol. Nine GC-MS and five HPLC-identified metabolites had shown interactions with 45 and 30 COVID-19-associated human proteins, respectively. Among the proteins, PARP1, FN1, PRKCA, EGFR, ALDH2, AKR1C3, AHR, and IKBKB have been found as potential therapeutic targets to mitigate SARS-CoV-2 infection. KEGG pathway analysis also showed a strong association of FN1, EGFR, and IKBKB genes with SARS-CoV-2 viral replication and cytokine overexpression due to viral infection. Protein-protein interaction (PPI) analysis also showed that TP53, MMP9, FN1, EGFR, and NOS2 proteins are highly related to the genes involved in COVID-19 comorbidity. These proteins showed interaction with the plant phytoconstituents as well. As the study offers a robust network-based procedure for identifying biomolecules relevant to COVID-19 disease, A. malaccensis could be a good source of effective therapeutic agents against COVID-19 and related viral diseases.

Identification of Differentially Expressed Genes and Protein-Protein Interaction in Patients With COVID-19 and Diabetes Peripheral Neuropathy: A Bioinformatics and System Biology Approach.

The coronavirus disease 2019 (COVID-19) pandemic has had a significant impact globally, resulting in a higher death toll and persistent health issues for survivors, particularly those with pre-existing medical conditions. Numerous studies have demonstrated a strong correlation between catastrophic COVID-19 results and diabetes. To gain deeper insights, we analysed the transcriptome dataset from COVID-19 and diabetic peripheral neuropathic patients. Using the R programming language, differentially expressed genes (DEGs) were identified and classified based on up and down regulations. The overlaps of DEGs were then explored between these groups. Functional annotation of those common DEGs was performed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Bio-Planet, Reactome, and Wiki pathways. A protein-protein interaction (PPI) network was created with bioinformatics tools to understand molecular interactions. Through topological analysis of the PPI network, we determined hub gene modules and explored gene regulatory networks (GRN). Furthermore, the study extended to suggesting potential drug molecules for the identified mutual DEG based on the comprehensive analysis. These approaches may contribute to understanding the molecular intricacies of COVID-19 in diabetic peripheral neuropathy patients through insights into potential therapeutic interventions.

Regulation of RHOV signaling by interaction with SH3 domain-containing adaptor proteins and phosphorylation by PKA.

RHOV and RHOU are considered atypical Rho-family small GTPases because of the existence of N- and C-terminal extension regions, abnormal GDP/GTP cycling, and post-translational modification. Particularly, RHOV and RHOU both have a proline-rich (PR) motif in the N-terminal region. It has been reported that the PR motif of RHOU interacts with GRB2, a SH3 domain-containing adaptor protein, and regulates its activity through EGF receptor signaling. However, it is unknown whether RHOV, like RHOU, interacts with SH3 domain-containing adaptor proteins. In this study, we investigated the interactions between RHOV and SH3 domain-containing adaptor proteins, including GRB2 and NCK2. The RHOV-induced serum response factor (SRF)-dependent gene transcriptional activity was attenuated in cells co-expressing either GRB2 or NCK2 compared to cells expressing RHOV alone. From the results of experiments using various gene mutants of RHOV and GRB2, it appears that the PR motif of the N-terminal region of RHOV is the crucial binding site for the SH3 domain-containing proteins. Furthermore, we found that Ser25 in the N-terminal region of RHOV is phosphorylated by PKA and that its phosphorylation is suppressed by interaction with NCK2 but not GRB2. We have found a novel regulatory mechanism for the phosphorylation of RHOV and its interaction with SH3 domain-containing adaptor proteins.

Characterization of the UDP-glycosyltransferase UGT33D1 in silkworm Bombyx mori.

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are important metabolizing enzymes functioning by adding a sugar moiety to a small lipophilic substrate molecule and play critical roles in drug/toxin metabolism for all realms of life. In this study, the silkworm Bombyx mori UGT33D1 gene was characterized in detail. UGT33D1 was found localized in the endoplasmic reticulum (ER) compartment just like other animal UGTs and was mainly expressed in the silkworm midgut. We first reported that UGT33D1 was important to BmNPV infection, as silencing UGT33D1 inhibited the BmNPV infection in silkworm BmN cells, while overexpressing the gene promoted viral infection. The molecular pathways regulated by UGT33D1 were analysed via transcriptome sequencing upon UGT33D1 knockdown, highlighting the important role of the gene in maintaining a balanced oxidoreductive state of the organism. In addition, proteins that physically interact with UGT33D1 were identified through immunoprecipitation and mass spectrometry analysis, which includes tubulin, elongation factor, certain ribosomal proteins, histone proteins and zinc finger proteins that had been previously reported for human UGT-interacting proteins. This study provided preliminary but important functional information on UGT33D1 and is hoped to trigger deeper investigations into silkworm UGTs and their functional mechanisms.

Interface-guided phenotyping of coding variants in the transcription factor RUNX1.

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

Molecular mechanism of contactin 2 homophilic interaction.

Contactin 2 (CNTN2) is a cell adhesion molecule involved in axon guidance, neuronal migration, and fasciculation. The ectodomains of CNTN1-CNTN6 are composed of six Ig domains (Ig1-Ig6) and four FN domains. Here, we show that CNTN2 forms transient homophilic interactions (KD ∼200 nM). Cryo-EM structures of full-length CNTN2 and CNTN2_Ig1-Ig6 reveal a T-shaped homodimer formed by intertwined, parallel monomers. Unexpectedly, the horseshoe-shaped Ig1-Ig4 headpieces extend their Ig2-Ig3 tips outwards on either side of the homodimer, while Ig4, Ig5, Ig6, and the FN domains form a central stalk. Cross-linking mass spectrometry and cell-based binding assays confirm the 3D assembly of the CNTN2 homodimer. The interface mediating homodimer formation differs between CNTNs, as do the homophilic versus heterophilic interaction mechanisms. The CNTN family thus encodes a versatile molecular platform that supports a very diverse portfolio of protein interactions and that can be leveraged to strategically guide neural circuit development.

[Applications of native mass spectrometry and ultraviolet photodissociation in protein structure and interaction analysis].

Dynamic changes in the structures and interactions of proteins are closely correlated with their biological functions. However, the precise detection and analysis of these molecules are challenging. Native mass spectrometry (nMS) introduces proteins or protein complexes into the gas phase by electrospray ionization, and then performs MS analysis under near-physiological conditions that preserve the folded state of proteins and their complexes in solution. nMS can provide information on stoichiometry, assembly, and dissociation constants by directly determining the relative molecular masses of protein complexes through high-resolution MS. It can also integrate various MS dissociation technologies, such as collision-induced dissociation (CID), surface-induced dissociation (SID), and ultraviolet photodissociation (UVPD), to analyze the conformational changes, binding interfaces, and active sites of protein complexes, thereby revealing the relationship between their interactions and biological functions. UVPD, especially 193 nm excimer laser UVPD, is a rapidly evolving MS dissociation method that can directly dissociate the covalent bonds of protein backbones with a single pulse. It can generate different types of fragment ions, while preserving noncovalent interactions such as hydrogen bonds within these ions, thereby enabling the MS analysis of protein structures with single-amino-acid-site resolution. This review outlines the applications and recent progress of nMS and UVPD in protein dynamic structure and interaction analyses. It covers the nMS techniques used to analyze protein-small-molecule ligand interactions, the structures of membrane proteins and their complexes, and protein-protein interactions. The discussion on UVPD includes the analysis of gas-phase protein structures and interactions, as well as alterations in protein dynamic structures, and interactions resulting from mutations and ligand binding. Finally, this review describes the future development prospects for protein analysis by nMS and new-generation advanced extreme UV light sources with higher brightness and shorter pulses.

MediatorWeb: a protein-protein interaction network database for the RNA polymerase II Mediator complex.

The protein-protein interaction (PPI) network of the Mediator complex is very tightly regulated and depends on different developmental and environmental cues. Here, we present an interactive platform for comparative analysis of the Mediator subunits from humans, baker's yeast Saccharomyces cerevisiae, and model plant Arabidopsis thaliana in a user-friendly web-interface database called MediatorWeb. MediatorWeb provides an interface to visualize and analyze the PPI network of Mediator subunits. The database facilitates downloading the untargeted and unweighted network of Mediator complex, its submodules, and individual Mediator subunits to better visualize the importance of individual Mediator subunits or their submodules. Further, MediatorWeb offers network visualization of the Mediator complex and interacting proteins that are functionally annotated. This feature provides clues to understand functions of Mediator subunits in different processes. In an additional tab, MediatorWeb provides quick access to secondary and tertiary structures, as well as residue-level contact information for Mediator subunits in each of the three model organisms. Another useful feature of MediatorWeb is detection of interologs based on orthologous analyses, which can provide clues to understand the functions of Mediator complex in less explored kingdoms. Thus, MediatorWeb and its features can help the user to understand the role of Mediator complex and its subunits in the transcription regulation of gene expression.

A high-efficiency transient expression system mediated by Agrobacterium tumefaciens in Spinacia oleracea leaves.

BACKGROUND: Optimization of a highly efficient transient expression system is critical for the study of gene function, particularly in those plants in which stable transformation methods are not widely available. Agrobacterium tumefaciens­mediated transient transformation is a simple and low-cost method that has been developed and applied to a wide variety of plant species. However, the transient expression in spinach (Spinacia oleracea L.) is still not reported. RESULTS: We developed a transient expression system in spinach leaves of the Sp75 and Sp73 varieties. Several factors influencing the transformation efficiency were optimized such as Agrobacterium strain, spinach seedling stage, leaf position, and the expression time after injection. Agrobacterium strain GV3101 (pSoup-p19) was more efficient than AGL1 in expressing recombinant protein in spinach leaves. In general, Sp75 leaves were more suitable than Sp73 leaves, regardless of grow stage. At four-leaf stage, higher intensity and efficiency of transient expression were observed in group 1 (G1) of Sp75 at 53 h after injection (HAI) and in G1 of Sp73 at 64 HAI. At six-leaf stage of Sp75, group 3 (G3) at 72 HAI were the most effective condition for transient expression. Using the optimized expression system, we detected the subcellular localization of a transcriptional co-activator SoMBF1c and a NADPH oxidase SoRbohF. We also detected the interaction of the protein kinase SoCRK10 and the NADPH oxidase SoRbohB. CONCLUSION: This study established a method of highly efficient transient expression mediated by Agrobacterium in spinach leaves. The transient expression system will facilitate the analysis of gene function and lay a solid foundation for molecular design breeding of spinach.

The identification, evolutionary analysis, and immune roles of Rab family members in red swamp crayfish, Procambarus clarkii.

The Rab GTPase constitutes the largest family of small GTPases that regulate intracellular trafficking. Different eukaryotes possess varying numbers of Rab paralogs. However, limited knowledge exists regarding the evolutionary pattern of Rab family in most major eukaryotic supergroups. This study cloned 24 Rab genes from transcriptome data of Procambarus clarkii haemocytes. The multiple sequence alignment and phylogenetic tree analysis revealed a relatively high degree of conservation for PcRab. Furthermore, PcRab exhibited similarities in motif composition with all members showing presence of G, PM, RabF, and RabSF motifs. The tertiary structure indicated that PcRab proteins mainly consisted of α-helices and ß-strands, and most PcRab proteins shared similar tertiary structures, and it was indicated that they have similar protein characteristics. Protein-protein interaction prediction identified a total of 20 interacting proteins involved in vesicle trafficking, phagocytosis, and signal transduction with 193 interactions. Expression analysis showed wide expression patterns for PcRab in P. clarkii organs. Upon infection by white spot syndrome virus and Aeromonas veronii, significant induction was observed for PcRab gene expression levels, indicating their involvement in pathogen response mechanisms. The present study represents the pioneering effort in comprehensively identifying and cloning the Rab family genes in crustacean, followed by a systematic investigation into their evolutionary patterns and immune response upon pathogen infection. The results provided valuable insights for further investigation into the molecular mechanism underlying the response of P. clarkii to pathogen infection.

A proximity proteomics pipeline with improved reproducibility and throughput.

Proximity labeling (PL) via biotinylation coupled with mass spectrometry (MS) captures spatial proteomes in cells. Large-scale processing requires a workflow minimizing hands-on time and enhancing quantitative reproducibility. We introduced a scalable PL pipeline integrating automated enrichment of biotinylated proteins in a 96-well plate format. Combining this with optimized quantitative MS based on data-independent acquisition (DIA), we increased sample throughput and improved protein identification and quantification reproducibility. We applied this pipeline to delineate subcellular proteomes across various compartments. Using the 5HT2A serotonin receptor as a model, we studied temporal changes of proximal interaction networks induced by receptor activation. In addition, we modified the pipeline for reduced sample input to accommodate CRISPR-based gene knockout, assessing dynamics of the 5HT2A network in response to perturbation of selected interactors. This PL approach is universally applicable to PL proteomics using biotinylation-based PL enzymes, enhancing throughput and reproducibility of standard protocols.

Simulated complexes formed from a set of postsynaptic proteins suggest a localised effect of a hypomorphic Shank mutation.

BACKGROUND: The postsynaptic density is an elaborate protein network beneath the postsynaptic membrane involved in the molecular processes underlying learning and memory. The postsynaptic density is built up from the same major proteins but its exact composition and organization differs between synapses. Mutations perturbing protein: protein interactions generally occurring in this network might lead to effects specific for cell types or processes, the understanding of which can be especially challenging. RESULTS: In this work we use systems biology-based modeling of protein complex distributions in a simplified set of major postsynaptic proteins to investigate the effect of a hypomorphic Shank mutation perturbing a single well-defined interaction. We use data sets with widely variable abundances of the constituent proteins. Our results suggest that the effect of the mutation is heavily dependent on the overall availability of all the protein components of the whole network and no trivial correspondence between the expression level of the directly affected proteins and overall complex distribution can be observed. CONCLUSIONS: Our results stress the importance of context-dependent interpretation of mutations. Even the weakening of a generally occurring protein: protein interaction might have well-defined effects, and these can not easily be predicted based only on the abundance of the proteins directly affected. Our results provide insight on how cell-specific effects can be exerted by a mutation perturbing a generally occurring interaction even when the wider interaction network is largely similar.

Multilayered proteomics reveals that JAM-A promotes breast cancer progression via regulation of amino acid transporter LAT1.

Recent studies have shown that transmembrane-type tight junction proteins are upregulated in various cancers compared with their levels in normal tissues and are involved in cancer progression, suggesting that they are potential therapeutic targets. Here, we demonstrated the expression profile and a novel role of junctional adhesion molecule-A (JAM-A) in breast cancer. Immunohistochemistry of surgical specimens showed that JAM-A was highly expressed from carcinoma in situ lesions, as in other adenocarcinomas, with higher expression in invasive carcinomas. High expression of JAM-A contributed to malignant aspects such as lymph node metastasis and lymphatic involvement positivity. In breast cancer cells, JAM-A expression status affects malignant potentials including proliferation and migration. Multilayered proteomics revealed that JAM-A interacts with the amino acid transporter LAT1 in breast cancer cells. JAM-A regulates the expression of LAT1 and interacts with it on the whole cell membrane, leading to enhanced amino acid uptake to promote tumor growth. Double high expression of JAM-A and LAT1 predicts poor prognosis in patients with breast cancer. Of note, an antibody against an extracellular domain of JAM-A suppressed the proliferation of breast cancer cells. Our findings indicate the possibility of JAM-A-targeted therapy ideally combined with LAT1-targeted therapy as a new therapeutic strategy against breast cancer.

What Can Be Learned from the Partitioning Behavior of Proteins in Aqueous Two-Phase Systems?

This review covers the analytical applications of protein partitioning in aqueous two-phase systems (ATPSs). We review the advancements in the analytical application of protein partitioning in ATPSs that have been achieved over the last two decades. Multiple examples of different applications, such as the quality control of recombinant proteins, analysis of protein misfolding, characterization of structural changes as small as a single-point mutation, conformational changes upon binding of different ligands, detection of protein-protein interactions, and analysis of structurally different isoforms of a protein are presented. The new approach to discovering new drugs for a known target (e.g., a receptor) is described when one or more previous drugs are already available with well-characterized biological efficacy profiles.

Mechanism of oxymatrine in the treatment of cryptosporidiosis through TNF/NF-κB signaling pathway based on network pharmacology and experimental validation.

Cryptosporidiosis is a worldwide zoonotic disease. Oxymatrine, an alkaloid extracted and isolated from the plant bitter ginseng, has been reported to have therapeutic effects on cryptosporidiosis. However, the underlying mechanism of its action remains unclear. In this study, we utilized network pharmacology and experimental validation to investigate the mechanism of oxymatrine in the treatment of cryptosporidiosis. First, the potential targets of drugs and diseases were predicted by TCMSP, Gene Cards, and other databases. Following the intersection of drug-disease targets, the DAVID database was used to implement the enrichment analysis of GO functions and KEGG pathways, and then the network diagram of "intersected target-KEGG" relationship was constructed. Autodock 4.2.6 software was used to carry out the molecular docking of core targets to drug components. Based on the establishment of a mouse model of cryptosporidiosis, the validity of the targets in the TNF/NF-κB signaling pathway was confirmed using Western blot analysis and Quantitative Rea-ltime-PCR. A total of 41 intersectional targets of oxymatrine and Cryptosporidium were generated from the results, and five core targets were screened out by network analysis, including RELA, AKT1, ESR1, TNF, and CASP3. The enrichment analysis showed that oxymatrine could regulate multiple gene targets, mediate TNF, Apoptpsis, IL-17, NF-κB and other signaling pathways. Molecular docking experiments revealed that oxymatrine was tightly bound to core targets with stable conformation. Furthermore, we found through animal experiments that oxymatrine could regulate the mRNA and protein expression of IL-6, NF-κB, and TNF-α in the intestinal tissues of post-infected mice through the TNF/NF-κB signaling pathway. Therefore, it can be concluded that oxymatrine can regulate the inflammatory factors TNF-α, NF-κB, and IL-6 through the TNF/NF-κB signaling pathway for the treatment of cryptosporidiosis. This prediction has also been validated by network pharmacology and animal experiments.

N-terminal truncated phospholipase A1 accessory protein PlaS from Serratia marcescens alleviates inhibitory on host cell growth and enhances PlaA1 enzymatic activity.

Phospholipase A1 (PLA1) is a kind of specific phospholipid hydrolase widely used in food, medical, textile. However, limitations in its expression and enzymatic activity have prompted the investigation of the phospholipase-assisting protein PlaS. In this study, we elucidate the role of PlaS in enhancing the expression and activity of PlaA1 through N-terminal truncation. Our research demonstrates that truncating the N-terminal region of PlaS effectively overcomes its inhibitory effect on host cells, resulting in improved cell growth and increased protein solubility of the protein. The yeast two-hybrid assay confirms the interaction between PlaA1 and N-terminal truncated PlaS (∆N27 PlaS), highlighting their binding capabilities. Furthermore, in vitro studies using Biacore analysis reveal a concentration-dependent and specific binding between PlaA1 and ∆N27 PlaS, exhibiting high affinity. Molecular docking analysis provides insights into the hydrogen bond interactions between ∆N27 PlaS and PlaA1, identifying key amino acid residues crucial for their binding. Finally, the enzyme activity of PLA1 was boost to 8.4 U/mL by orthogonal test. Study significantly contributes to the understanding of the interaction mechanism between PlaS and PlaA1, offering potential strategies for enhancing PlaA1 activity through protein engineering approaches.

Galectin-3-Binding Protein Inhibits Extracellular Heparan 6-O-Endosulfatase Sulf-2.

Human extracellular 6-O-endosulfatases Sulf-1 and Sulf-2 are the only enzymes that post-synthetically alter the 6-O sulfation of heparan sulfate proteoglycans (HSPG), which regulates interactions of HSPG with many proteins. Oncogenicity of Sulf-2 in different cancers has been documented, and we have shown that Sulf-2 is associated with poor survival outcomes in head and neck squamous cell carcinoma (HNSCC). Despite its importance, limited information is available on direct protein-protein interactions of the Sulf-2 protein in the tumor microenvironment. In this study, we used monoclonal antibody (mAb) affinity purification and mass spectrometry to identify galectin-3-binding protein (LG3BP) as a highly specific binding partner of Sulf-2 in the conditioned media of HNSCC cell lines. We validated their direct interaction in vitro using recombinant proteins and have shown that the chondroitin sulfate (CS) covalently bound to the Sulf-2 influences the binding to LG3BP. We confirmed the importance of the CS chain for the interaction by generating a mutant Sulf-2 protein that lacks the CS. Importantly, we have shown that the LG3BP inhibits Sulf-2 activity in vitro in a concentration-dependent manner. As a consequence, the addition of LG3BP to a spheroid cell culture inhibited the invasion of the HNSCC cells into Matrigel. Thus, Sulf-2 interaction with LG3BP may regulate the physiological activity of the Sulf-2 enzyme as well as its activity in the tumor microenvironment.

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

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.