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

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Phospholipase C: Diverse functions in plant biotic stress resistance and fungal pathogenicity

Phospholipase C (PLC) generates various second messenger molecules and mediates phospholipid hydrolysis. In recent years, the important roles of plant and fungal PLC in disease resistance and pathogenicity, respectively, have been determined. However, the roles of PLC in plants and fungi are unintegrated and relevant literature is disorganized. This makes it difficult for researchers to implement PLC-based strategies to improve disease resistance in plants. In this comprehensive review, we summarize the structure, classification, and phylogeny of the PLCs involved in plant biotic stress resistance and fungal pathogenicity. PLCs can be divided into two groups, nonspecific PLC (NPC) and phosphatidylinositol-specific PLC (PI-PLC), which present marked differences in phylogenetic evolution. The products of PLC genes in fungi play significant roles in physiological activity and pathogenesis, whereas those encoded by plant PLC genes mediate the immune response to fungi. This review provides a perspective for the future control of plant fungal diseases.

Drought and heat stress mediated activation of lipid signaling in plants: a critical review

Lipids are a principal component of plasma membrane, acting as a protective barrier between the cell and its surroundings. Abiotic stresses such as drought and temperature induce various lipid-dependent signaling responses, and the membrane lipids respond differently to environmental challenges. Recent studies have revealed that lipids serve as signal mediators forreducing stress responses in plant cells and activating defense systems. Signaling lipids, such as phosphatidic acid, phosphoinositides, sphingolipids, lysophospholipids, oxylipins, and N-acylethanolamines, are generated in response to stress. Membrane lipids are essential for maintaining the lamellar stack of chloroplasts and stabilizing chloroplast membranes under stress. However, the effects of lipid signaling targets in plants are not fully understood. This review focuses on the synthesis of various signaling lipids and their roles in abiotic stress tolerance responses, providing an essential perspective for further investigation into the interactions between plant lipids and abiotic stress.

Grapevine-Associated Lipid Signalling Is Specifically Activated in an Rpv3 Background in Response to an Aggressive P. viticola Pathovar

Vitis vinifera L. is highly susceptible to the biotrophic pathogen Plasmopara viticola. To control the downy mildew disease, several phytochemicals are applied every season. Recent European Union requirements to reduce the use of chemicals in viticulture have made it crucial to use alternative and more sustainable approaches to control this disease. Our previous studies pinpoint the role of fatty acids and lipid signalling in the establishment of an incompatible interaction between grapevine and P. viticola. To further understand the mechanisms behind lipid involvement in an effective defence response we have analysed the expression of several genes related to lipid metabolism in three grapevine genotypes: Chardonnay (susceptible); Regent (tolerant), harbouring an Rpv3-1 resistance loci; and Sauvignac (resistant) that harbours a pyramid of Rpv12 and Rpv3-1 resistance loci. A highly aggressive P. viticola isolate was used (NW-10/16). Moreover, we have characterised the grapevine phospholipases C and D gene families and monitored fatty acid modulation during infection. Our results indicate that both susceptible and resistant grapevine hosts did not present wide fatty acid or gene expression modulation. The modulation of genes associated with lipid signalling and fatty acids seems to be specific to Regent, which raises the hypothesis of being specifically linked to the Rpv3 loci. In Sauvignac, the Rpv12 may be dominant concerning the defence response, and, thus, this genotype may present the activation of other pathways rather than lipid signalling.

Identification and expression analysis of phospholipase C family genes between different male fertility accessions in pepper

Phospholipase C (PLC) is one of the major lipid-hydrolyzing enzymes, involved in lipid-mediating signal pathway. PLCs have been found to play a significant role in the growth and development of plants. In this study, the genome-wide identification and characteristic analysis of CaPLC family genes in pepper were conducted and the expression of two CaPLC genes were investigated. The results showed that a total of 11 CaPLC family genes were systematically identified, which were distributed on five chromosomes and divided into two groups based on their evolutionary relevance. Some cis-elements responding to different hormones and stresses were screened in the promoters of CaPLC genes. Quantitative real-time PCR indicated that the expression of CaPIPLC1 and CaPIPLC5 in flowers were dozens of times higher than in other tissues. In addition, with the development of flower buds, the relative expressions of CaPIPLC1 and CaPIPLC5 gradually increased in fertile materials R1 and F₁. However, no expression of CaPIPLC1 and CaPIPLC5 were detected at all developmental stages of cytoplasmic male sterile lines (CMS) compared with fertile accessions. The study revealed the number and characteristics of the CaPLC family genes, which supplied a basic and systematic understanding of CaPLC family. In addition, these findings provided new insights into the role of CaPLC genes in pollen development and fertility restoration in pepper.

32Pi Labeled Transgenic Wheat Shows the Accumulation of Phosphatidylinositol 4,5-bisphosphate and Phosphatidic Acid Under Heat and Osmotic Stress

The ensuing heat stress drastically affects wheat plant growth and development, consequently compromising its grain yield. There are many thermoregulatory processes/mechanisms mediated by ion channels, lipids, and lipid-modifying enzymes that occur in the plasma membrane and the chloroplast. With the onset of abiotic or biotic stresses, phosphoinositide-specific phospholipase C (PI-PLC), as a signaling enzyme, hydrolyzes phosphatidylinositol 4,5-bisphosphate (PIP₂) to generate inositol 1,4,5-trisphosphate (IP₃) and diacylglycerol (DAG) which is further phosphorylated into phosphatidic acid (PA) as a secondary messenger and is involved in multiple processes. In the current study, a phospholipase C (PLC) signaling pathway was investigated in spring wheat (Triticum aestivum L.) and evaluated its four AtPLC5 overexpressed (OE)/transgenic lines under heat and osmotic stresses through ³²P_i radioactive labeling. Naturally, the wheat harbors only a small amount of PIP₂. However, with the sudden increase in temperature (40°C), PIP₂ levels start to rise within 7.5 min in a time-dependent manner in wild-type (Wt) wheat. While the Phosphatidic acid (PA) level also elevated up to 1.6-fold upon exposing wild-type wheat to heat stress (40°C). However, at the anthesis stage, a significant increase of ∼4.5-folds in PIP₂ level was observed within 30 min at 40°C in AtPLC5 over-expressed wheat lines. Significant differences in PIP₂ level were observed in Wt and AtPLC5-OE lines when treated with 1200 mM sorbitol solution. It is assumed that the phenomenon might be a result of the activation of PLC/DGK pathways. Together, these results indicate that heat stress and osmotic stress activate several lipid responses in wild-type and transgenic wheat and can explain heat and osmotic stress tolerance in the wheat plant.

Emerging role of phospholipase C mediated lipid signaling in abiotic stress tolerance and development in plants

Environmental stimuli are primarily perceived at the plasma membrane. Stimuli perception leads to membrane disintegration and generation of molecules which trigger lipid signaling. In plants, lipid signaling regulates important biological functions however, the molecular mechanism involved is unclear. Phospholipases C (PLCs) are important lipid-modifying enzymes in eukaryotes. In animals, PLCs by hydrolyzing phospholipids, such as phosphatidylinositol-4,5-bisphosphate [PI(4,5)P₂] generate diacylglycerol (DAG) and inositol- 1,4,5-trisphosphate (IP₃). However, in plants their phosphorylated variants i.e., phosphatidic acid (PA) and inositol hexakisphosphate (IP₆) are proposed to mediate lipid signaling. Specific substrate preferences divide PLCs into phosphatidylinositol-PLC (PI-PLC) and non-specific PLCs (NPC). PLC activity is regulated by various cellular factors including, calcium (Ca²⁺) concentration, phospholipid substrate, and post-translational modifications. Both PI-PLCs and NPCs are implicated in plants' response to stresses and development. Emerging evidences show that PLCs regulate structural and developmental features, like stomata movement, microtubule organization, membrane remodelling and root development under abiotic stresses. Thus, crucial insights are provided into PLC mediated regulatory mechanism of abiotic stress responses in plants. In this review, we describe the structure and regulation of plant PLCs. In addition, cellular and physiological roles of PLCs in abiotic stresses, phosphorus deficiency, aluminium toxicity, pollen tube growth, and root development are discussed.

Genome-wide characterization and expression profiling of the Phospholipase C (PLC) gene family in three orchids of economic importance

Background

Phospholipases hydrolyze glycerophospholipids and generate diverse lipid-derived molecules with secondary messenger activity. Out of these, phospholipase C (PLC) specifically cleaves the phospholipids at ester linkages and yields diacylglycerol (DAG) and phosphorylated head groups. PLCs are classified further as phosphatidylinositol-specific PLCs (PI-PLCs) and non-specific PLCs with biased specificity for phosphatidylcholine (NPC/PC-PLC).

Results

In the present report, we identified and characterized PLC genes in the genomes of three orchids, Phalaenopsis equestris (seven PePLCs), Dendrobium catenatum (eight DcPLCs), and Apostasia shenzhenica (seven AsPLCs). Multiple sequence alignment analysis confirmed the presence of conserved X and Y catalytic domains, calcium/lipid-binding domain (C2 domain) at the C terminal region, and EF-hand at the N-terminal region in PI-PLC proteins and esterase domain in PC-PLC. Systematic phylogenetic analysis established the relationship of the PLC protein sequences and clustered them into two groups (PI-PLC and PC-PLC) along with those of Arabidopsis thaliana and Oryza sativa. Gene architecture studies showed the presence of nine exons in all PI-PLC genes while the number varied from one to five in PC-PLCs. RNA-seq-based spatio-temporal expression profile for PLC genes was generated, which showed that PePC-PLC1, PePC-PLC2A, DcPC-PLC1A, DcPC-PLC1B, DcPC-PLC2, DcPC-PLC1B, and AsPC-PLC1 had significant expression in all reproductive and vegetative tissues. The expression profile is matched to their upstream cis-regulatory promoter elements, which indicates that PLC genes have a role in various growth and development processes and during stress responses.

Conclusions

The present study unwrapped the opportunity for functional characterization of selected PLC genes in planta for plant improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00217-z.

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.

GhPIPLC2D promotes cotton fiber elongation by enhancing ethylene biosynthesis

Inositol-1,4,5-trisphosphate (IP₃) is an important second messenger and one of the products of phosphoinositide-specific phospholipase C (PIPLC)-mediated phosphatidylinositol (4,5) bisphosphate (PIP₂) hydrolysis. However, the function of IP₃ in cotton is unknown. Here, we characterized the function of GhPIPLC2D in cotton fiber elongation. GhPIPLC2D was preferentially expressed in elongating fibers. Suppression of GhPIPLC2D transcripts resulted in shorter fibers and decreased IP₃ accumulation and ethylene biosynthesis. Exogenous application of linolenic acid (C18:3) and phosphatidylinositol (PI), the precursor of IP₃, improved IP₃ and myo-inositol-1,2,3,4,5,6-hexakisphosphate (IP₆) accumulation, as well as ethylene biosynthesis. Moreover, fiber length in GhPIPLC2D-silenced plant was reduced after exogenous application of IP₆ and ethylene. These results indicate that GhPIPLC2D positively regulates fiber elongation and IP₃ promotes fiber elongation by enhancing ethylene biosynthesis. Our study broadens our understanding of the function of IP₃ in cotton fiber elongation and highlights the possibility of cultivating better cotton varieties by manipulating GhPIPLC2D in the future.

•

GhPIPLC2D positively regulates cotton fiber elongation

•

GhPIPLC2D cleaves PIP₂ into IP₃, which could be phosphorylated to IP₆

•

IP₆ enhances fiber elongation via improving ethylene biosynthesis

Genome-Wide Investigation of the Phospholipase C Gene Family in Zea mays

Phospholipase C (PLC) is one of the main hydrolytic enzymes in the metabolism of phosphoinositide and plays an important role in a variety of signal transduction processes responding to plant growth, development, and stress. Although the characteristics of many plant PLCs have been studied, PLC genes of maize have not been comprehensively identified. According to the study, five phosphatidylinositol-specific PLC (PI-PLC) and six non-specific PLC (NPC) genes were identified in maize. The PI-PLC and NPC genes of maize are conserved compared with homologous genes in other plants, especially in evolutionary relationship, protein sequences, conserved motifs, and gene structures. Transient expression of ZmPLC-GFP fusion protein in Arabidopsis protoplast cells showed that ZmPLCs are multi-localization. Analyses of transcription levels showed that ZmPLCs were significantly different under various different tissues and abiotic stresses. Association analysis shown that some ZmPLCs significantly associated with agronomic traits in 508 maize inbred lines. These results contribute to study the function of ZmPLCs and to provide good candidate targets for the yield and quality of superior maize cultivars.

Two Aquaporin Genes, GhPIP2;7 and GhTIP2;1, Positively Regulate the Tolerance of Upland Cotton to Salt and Osmotic Stresses

Aquaporins (AQPs) facilitate the transport of water and small molecules across intrinsic membranes and play a critical role in abiotic stresses. In this study, 111, 54, and 56 candidate AQP genes were identified in Gossypium hirsutum (AD₁), Gossypium arboreum (A₂), and Gossypium raimondii (D₅), respectively, and were further classified into five subfamilies, namely, plasma intrinsic protein (PIP), tonoplast intrinsic protein (TIP), nodulin 26-like intrinsic protein (NIP), small basic intrinsic protein (SIP), and uncategorized X intrinsic protein (XIP). Transcriptome analysis and quantitative real-time PCR (qRT-PCR) revealed some high-expression GhPIPs and GhTIPs (PIP and TIP genes in G. hirsutum, respectively) in drought and salt stresses. GhPIP2;7-silenced plants decreased in the chlorophyll content, superoxide dismutase (SOD) activity, and peroxidase (POD) activity comparing the mock control (empty-vector) under 400 mM NaCl treatment, which indicated a positive regulatory role of GhPIP2;7 in salt tolerance of cotton. The GhTIP2;1-silenced cotton plants were more sensitive to osmotic stress. GhTIP2;1-overexpressed plants exhibited less accumulation of H₂O₂ and malondialdehyde but higher proline content under osmotic stress. In summary, our study elucidates the positive regulatory roles of two GhAQPs (GhPIP2;7 and GhTIP2;1) in salt and osmotic stress responses, respectively, and provides a new gene resource for future research.

Expression and evolution of the phospholipase C gene family in Brachypodium distachyon

BACKGROUND

Phospholipase C (PLC) is an enzyme that hydrolyzes phospholipids and plays an important role in plant growth and development. The Brachypodium distachyon is a model plant of Gramineae, but the research on PLC gene family of Brachypodium has not been reported.

OBJECTIVE

This study was performed to identify the PLC family gene in Brachypodium and to determine the expression profiles of PLCs under the abiotic stress.

METHODS

Complete genome sequences and transcriptomes of Brachypodium were downloaded from the PLAZA. The hidden Markov model-based profile of the conserved PLC domain was submitted as a query to identify all potential PLC domain sequences with HMMER software. Expression profiles of BdPLCs were obtained based on the qRT-PCR analysis.

RESULTS

There were 8 PLC genes in Brachypodium (BdPI-PLCs 1-4 and BdNPCs 1-4). All members of BdPI-PLC had three conserved domains of X, Y, and C2, and no EF-hand was found. All BdNPCs contained a phosphatase domain. BdPI-PLC genes were distributed on Chr1, Chr2 and Chr4, with different types and numbers of cis-regulatory elements in their respective gene promoters. Phylogenetic analysis showed that the genetic relationship between Brachypodium and rice was closer than Arabidopsis. The expression patterns of BdPI-PLC gene under abiotic stresses (drought, low temperature, high temperature and salt stress) were up-regulated, indicated their important roles in response to low temperature, high temperature, drought and salt stresses.

CONCLUSIONS

This study provides comprehensive information for the study of Brachypodium PLC gene family and lays a foundation for further research on the molecular mechanism of Brachypodium stress adaptation.

Genome-Wide Identification and Expression Profile Analysis of the Phospholipase C Gene Family in Wheat (Triticum aestivum L.)

Phospholipid-hydrolyzing enzymes include members of the phospholipase C (PLC) family that play important roles in regulating plant growth and responding to stress. In the present study, a systematic in silico analysis of the wheat PLC gene family revealed a total of 26 wheat PLC genes (TaPLCs). Phylogenetic and sequence alignment analyses divided the wheat PLC genes into 2 subfamilies, TaPI-PLC (containing the typical X, Y, and C2 domains) and TaNPC (containing a phosphatase domain). TaPLC expression patterns differed among tissues, organs, and under abiotic stress conditions. The transcript levels of 8 TaPLC genes were validated through qPCR analyses. Most of the TaPLC genes were sensitive to salt stress and were up-regulated rapidly, and some were sensitive to low temperatures and drought. Overexpression of TaPI-PLC1-2B significantly improved resistance to salt and drought stress in Arabidopsis, and the primary root of P1-OE was significantly longer than that of the wild type under stress conditions. Our results not only provide comprehensive information for understanding the PLC gene family in wheat, but can also provide a solid foundation for functional characterization of the wheat PLC gene family.

Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Brassinosteroid (BR), a steroid phytohormone, whose signaling transduction pathways include a series of phosphorylation and dephosphorylation events, and GSK3s are the main negative regulator kinases. BRs have been shown to play vital roles in cotton fiber elongation. However, the underlying mechanism is still elusive. In this study, fibers of a BR-defective mutant Pagoda 1 (pag1), and its corresponding wild-type (ZM24) were selected for a comparative global phosphoproteome analysis at critical developmental time points: fast-growing stage (10 days after pollination (DPA)) and secondary cell wall synthesis stage (20 DPA). Based on the substrate characteristics of GSK3, 900 potential substrates were identified. Their GO and KEGG annotation results suggest that BR functions in fiber development by regulating GhSKs (GSK3s of Gossypium hirsutum L.) involved microtubule cytoskeleton organization, and pathways of glucose, sucrose and lipid metabolism. Further experimental results revealed that among the GhSK members identified, GhSK13 not only plays a role in BR signaling pathway, but also functions in developing fiber by respectively interacting with an AP2-like ethylene-responsive factor GhAP2L, a nuclear transcription factor Gh_DNF_YB19, and a homeodomain zipper member GhHDZ5. Overall, our phosphoproteomic research advances the understanding of fiber development controlled by BR signal pathways especially through GhSKs, and also offers numbers of target proteins for improving cotton fiber quality.

Non-specific phospholipase C (NPC): an emerging class of phospholipase C in plant growth and development

Non-specific phospholipase C (NPC) is a novel class of phospholipase C found only in bacteria and higher plants. NPC hydrolyzes major phospholipid classes such as phosphatidylcholine (PC) and phosphatidylethanolamine (PE) to produce diacylglycerol (DAG) and a corresponding phosphate-containing polar head group. Originally known as a toxin in certain bacteria to invade the host cell, this class of phospholipase has been well-investigated in bacteriology. Since the first discovery of eukaryotic NPC in Arabidopsis in 2005, this emerging class of phospholipase has received greater attention in plant biology in elucidating the biochemical characteristics and physiological function in the context of plant growth regulation and stress response. Particularly in the last few years, there has been significant progress made in understanding the fundamental character of 6 NPC isoforms in Arabidopsis, as well as novel function in other plant models. Now that research with plant NPC is entering into a new phase, this review aims to summarize recent progress in plant NPC along with some future perspectives.

The Arabidopsis thaliana non-specific phospholipase C2 is involved in the response to Pseudomonas syringae attack

Background and Aims

The non-specific phospholipase C (NPC) is a new member of the plant phospholipase family that reacts to abiotic environmental stresses, such as phosphate deficiency, high salinity, heat and aluminium toxicity, and is involved in root development, silicon distribution and brassinolide signalling. Six NPC genes (NPC1-NPC6) are found in the Arabidopsis genome. The NPC2 isoform has not been experimentally characterized so far.

Methods

The Arabidopsis NPC2 isoform was cloned and heterologously expressed in Escherichia coli. NPC2 enzyme activity was determined using fluorescent phosphatidylcholine as a substrate. Tissue expression and subcellular localization were analysed using GUS- and GFP-tagged NPC2. The expression patterns of NPC2 were analysed via quantitative real-time PCR. Independent homozygous transgenic plant lines overexpressing NPC2 under the control of a 35S promoter were generated, and reactive oxygen species were measured using a luminol-based assay.

Key Results

The heterologously expressed protein possessed phospholipase C activity, being able to hydrolyse phosphatidylcholine to diacylglycerol. NPC2 tagged with GFP was predominantly localized to the Golgi apparatus in Arabidopsis roots. The level of NPC2 transcript is rapidly altered during plant immune responses and correlates with the activation of multiple layers of the plant defence system. Transcription of NPC2 decreased substantially after plant infiltration with Pseudomonas syringae, flagellin peptide flg22 and salicylic acid treatments and expression of the effector molecule AvrRpm1. The decrease in NPC2 transcript levels correlated with a decrease in NPC2 enzyme activity. NPC2-overexpressing mutants showed higher reactive oxygen species production triggered by flg22.

Conclusions

This first experimental characterization of NPC2 provides new insights into the role of the non-specific phospholipase C protein family. The results suggest that NPC2 is involved in the response of Arabidopsis to P. syringae attack.

Structural, expression and evolutionary analysis of the non-specific phospholipase C gene family in Gossypium hirsutum

BACKGROUND

Nonspecific phospholipase C (NPC), which belongs to a phospholipase C subtype, is a class of phospholipases that hydrolyzes the primary membrane phospholipids, such as phosphatidylcholine, to yield sn-1, 2-diacylglycerol and a phosphorylated head-group. NPC plays multiple physiological roles in lipid metabolism and signaling in plants. To fully understand the putative roles of NPC genes in upland cotton, we cloned NPC genes from Gossypium hirsutum and carried out structural, expression and evolutionary analysis.

RESULTS

Eleven NPC genes were cloned from G. hirsutum, which were found on chromosomes scaffold269.1, D03, A07, D07, A08, D11, and scaffold3511_A13. All GhNPCs had typical phosphoesterase domains and have hydrolase activity that acts on ester bonds. GhNPCs were annotated as phospholipase C, which was involved in glycerophospholipid metabolism, ether lipid metabolism, and biosynthesis of secondary metabolites. These GhNPCs showed differential expression patterns in distinct plant tissues and in response to various types of stress (low-phosphate, salt, drought, and abscisic acid). They also had different types and numbers of cis-element. GhNPCs could be classified into four subfamilies. Four pairs of GhNPCs were generated by whole-genome duplication and they underwent purifying selection.

CONCLUSIONS

Our results suggested that GhNPCs are involved in regulating key abiotic stress responses and ABA signaling transduction, and they may have various functional roles for different members under complex abiotic stress conditions. Functional divergence may be the evolutionary driving force for the retention of four pairs of duplicate NPCs. Our analysis provides a solid foundation for the further functional characterization of the GhNPC gene family, and leads to potential applications in the genetic improvement of cotton cultivars.