Contact toxicity and transcriptomic analysis of terpinen‐4‐ol exposure in Tribolium castaneum

The terpene, terpinen-4-ol (T4ol), exhibits contact toxicity in Tribolium castaneum. However, the molecular mechanisms underlying this toxicity have not been elucidated. This study examined changes in the expression of four classic enzymes after exposure of T. castaneum to T4ol. Acetylcholinesterase and glutathione S-transferase activities were markedly inhibited after exposure to T4ol, while that of the detoxifying enzyme cytochrome oxidase P450 increased markedly. Carboxylesterase activity did not show significant changes. Furthermore, RNA sequencing revealed 260 differentially expressed genes (DEG) between the T4ol-treated and control samples, and qRT-PCR was used to validate the RNA-Seq data. The Gene Ontology analysis classified the DEGs into 36 functional groups, including the immune system processes, response to stimulus, and developmental processes. T4ol altered the response to stimulus and the immune system process of beetles by inducing the expression of the genes Stabilin-1, Attacin 1, and Defensin 1. Furthermore, the DEGs receptor tyrosine kinase Torso-like protein (RTKTsl), Frizzled 4 (Fz4), Protein Wnt-5b, Ecdysone-induced protein 78C (E78), Zinc finger protein GLIS1 (ZFPGLIS1) were classified as participating in beetle development, and Fz4 and Protein Wnt-5b also mapped to the Wnt signaling pathway. This indicated that pathways associated with development are inhibited after exposure to T4ol. T4ol also induced CYP9Z6/GSTs7 overexpression, and RNAi targeting these genes significantly increased larvae mortality on T4ol exposure, supporting the participation of CYP9Z6/GSTs7 in the response to T4ol in T. castaneum. The results of this study will facilitate understanding of the toxic mechanisms of T4ol and provide a basis for controlling the pests of stored products.     

Differential gene expression associated with flower development of mango (Mangifera indica L.) varieties with different shelf-life

Mango (Mangifera indica L.) is one of the most important commercial fruit crop grown in many parts of the world. Major challenges affecting mango trade are short shelf-life, high susceptibility to chilling injury, post-harvest diseases and consumer demand for improved fruit quality. The objective of the present study was to reveal the key regulators present in bud and flower tissues during flower development stage, associated with fruit development and affect the shelf-life of the mango fruit. RNA-sequencing of contrasting genotypes having short and long shelf-life, was carried out. Comparative differential expression pathway studies of long shelf-life (Totapuri) and short shelf-life (Bombay Green) mango genotypes revealed a total of 177 highly differentially expressed genes. Out of 177 total genes, 101 genes from endoplasmic reticulum pathway and very few from gibberellins (3) and jasmonic acid (1) pathway were identified. Genes from endoplasmic reticulum pathway like hsp 90, SRC2, DFRA, CHS, BG3 and ASPG1 mainly up regulated in Bombay Green. Uniprotein B9R8D3 also shows up regulation in Bombay Green. Ethylene insensitive pathway gene EIL1 up regulated in Bombay Green. Gene CAD1 from phenylpropanoid pathway mainly up regulated in Bombay Green. A total of 4 SSRs and 227 SNPs were mined from these pathways specific to the shelf-life. Molecular studies of endoplasmic reticulum, phenylpropanoid, ethylene, polygalacturonase and hormone pathways at the time of bud and flower formation revealed key regulators that determine the shelf-life of mango fruit.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

Potential mechanisms underlying the promoting effects of 3D collagen scaffold culture on stemness and drug resistance of glioma cells

Background3D collagen scaffold culture is a good tool to study glioma metastasis and recurrence in vitro.MethodsThe effect of 3D collagen culture on the colony formation, the sphere formation, and drug sensitivity of glioma cells was observed by soft-agar colony formation assays, sphere formation assays, and CCK-8 assays, respectively. 3D-glioma-drug genes were identified by previous results and online databases. Gene enrichment and PPI analyses were performed by R software and Metacsape. Hub 3D-glioma-drug genes were screened by STRING and Cytoscape. TCGA and CGGA databases and R software were used to analyze the distribution of hub genes in glioma and their effects on the prognosis. Western Blot was used to verify the effect of 3D collagen culture on the expression of hub genes. miRNAs targeting hub genes were predicted by ENCORI.Results3D collagen scaffold culture promoted colony formation, sphere formation, and drug resistance of glioma cells. There were 77 3D-glioma-drug genes screened, and the pathways enriched in the protein interaction network mainly included responses to stressors, DNA damage and repair, and drug metabolism. Hub 3D-glioma-drug genes were AKT1, ATM, CASP3, CCND1, EGFR, PARP1, and TP53. These genes and predicted miRNAs were expressed differentially in glioma samples and partially affected the prognosis of patients with glioma. These findings suggested these hub genes and miRNAs may play a key role in the effects generated by the 3D culture model and become new markers for glioma diagnosis and treatment.     

Differential expression and bioinformatics analysis of exosome circRNAs in pancreatic ductal adenocarcinoma

BackgroundPancreatic ductal adenocarcinoma (PDAC) is a fatal malignant tumor with an unfavorable prognosis. Increasing evidence indicated circRNAs were associated with the pathogenesis and progression of tumors, but data on the expression of serum exosomal circRNAs in PDAC are scarce. This study attempted to explore the prognostic value and function of serum exosomes in PDAC patients.MethodsMicroarray-based circRNA expression was determined in PDAC and paired with normal serum samples, and the intersection of differentially expressed circRNAs (DECs) in serum exosomal samples and GSE79634 tissue samples was conducted. A specific CircRNA database was applied to investigate DECs binding miRNAs. Target genes were predicted using the R package multiMiR. Cox regression analyses were applied for constructing a prognostic model. The immunological characteristics analysis was carried out through the TIMER, QUANTISEQ, XCELL, EPIC, and ssGSEA algorithms.Results15 DECs were finally identified, and a circRNA-miRNA-mRNA network was established. A prognostic risk model was developed to categorize patients according to the risk scores. Furthermore, the association between risk score and immune checkpoint genes including CD80, TNFSF9, CD276, CD274, LGALS9, and CD44 were significantly elevated in the high-risk group, while ICOSLG and ADORA2A were upregulated in the low-risk group.ConclusionsOur results may provide new clues for the prognosis and treatment of PDAC.     

Construction of miRNA-mRNA network for the identification of key biological markers and their associated pathways in IgA nephropathy by employing the integrated bioinformatics analysis

BackgroundAbout half-century ago, Immunoglobulin A nephropathy (IgAN) was discovered as a complicated disease with frequent clinical symptoms. Until now, exact mechanism underlying the pathogenesis of IgAN is poorly known. Therefore, current study was aimed to understand the molecular mechanism of IgAN by identifying the key miRNAs and their targeted hub genes. The key miRNAs might contribute to the diagnosis and therapy of IgAN, and could turn out to be a new star in the field of IgAN.MethodsThe microarray datasets were downloaded from Gene Expresssion Omnibus (GEO) database and analyzed using R package (LIMMA) in order to obtain differential expressed genes (DEGs). Then, the hub genes were identified using cytoHubba plugin of cytoscpae tool and other bioinformatics approaches including protein-protein interaction (PPI) network analysis, module analysis, and miRNA-hub gene network construction was also performed.ResultsA total of 348 DEGs were identified, of which 107 were upregulated genes and 241 were downregulated genes. Subsequently, the 12 overlapped genes were predicted from cytoHubba, and considered as hub genes. Moreover, a network among miRNA-hub genes was created to explore the correlation between the hub genes and their targeted miRNAs. Network construction ultimately lead to the identification of nine gene named FN1, EGR1, FOS, JUN, SERPINE1, MMP2, ATF3, MYC, and IL1B and one novel key miRNA namely, has-miR-144-3p as biomarker for diagnosis and therapy of IgAN.ConclusionThis study updates the information and yield a new perspective in context of understanding the pathogenesis and development of IgAN. In future, key miRNAs might be capable of improving the personalized detection and therapies for IgAN. In vivo and in vitro investigation of miRNAs and pathway interaction is essential to delineate the specific roles of the novel miRNAs, which may help to further reveal the mechanisms underlying IgAN.     

β-Carotene enhances the expression of inflammation-related genes and histone H3 K9 acetylation,K4 dimethylation,and K36 trimethylation around these genes in juvenile macrophage-like THP-1 cells

β-Carotene is converted into vitamin A in the body and can remove reactive oxygen species. However, it is still unclear whether β-carotene alters the expression levels of inflammation-related genes in macrophages and how this is regulated. In the present study, we investigated whether the administration of β-carotene under hyperglycemic conditions altered the expression level of inflammation-related genes and whether any observed differences were associated with changes in histone modifications in juvenile macrophage-like THP-1 cells. THP-1 cells (from a human monocytic leukemia cell line) were cultured in low glucose (5 mM), high glucose (25 mM), or high glucose (25 mM) + β-carotene (5 μM) media for 1 day, and mRNA expression levels of genes related to oxidative stress and inflammation, and histone modifications were determined by mRNA microarray and qRT-PCR analyses, and chromatin immunoprecipitation assays, respectively. The expression of inflammation-related genes, such as IL31RA, CD38, and NCF1B, and inflammation-associated signaling pathway genes, such as ITGAL, PRAM1, and CSF3R, were upregulated by β-carotene under high-glucose conditions. Under these conditions, histone H3 lysine 4 (K4) demethylation, H3K36 trimethylation, and H3K9 acetylation around the CD38, NCF1B, and ITGAL genes were higher in β-carotene-treated cells than in untreated cells. Treatment of juvenile macrophage-like THP-1 cells with β-carotene under these high glucose conditions induced the expression of inflammation-related genes, K9 acetylation, and K4 di- and K36 trimethylation of histone H3 around these genes.     

Analogs of nitrofuran antibiotics are potent GroEL/ES inhibitor pro-drugs

In two previous studies, we identified compound 1 as a moderate GroEL/ES inhibitor with weak to moderate antibacterial activity against Gram-positive and Gram-negative bacteria including Bacillus subtilis, methicillin-resistant Staphylococcus aureus, Klebsiella pneumonia, Acinetobacter baumannii, and SM101 Escherichia coli (which has a compromised lipopolysaccharide biosynthetic pathway making bacteria more permeable to drugs). Extending from those studies, we developed two series of analogs with key substructures resembling those of known antibacterials, nitroxoline (hydroxyquinoline moiety) and nifuroxazide/nitrofurantoin (bis-cyclic-N-acylhydrazone scaffolds). Through biochemical and cell-based assays, we identified potent GroEL/ES inhibitors that selectively blocked E. faecium, S. aureus, and E. coli proliferation with low cytotoxicity to human colon and intestine cells in vitro. Initially, only the hydroxyquinoline-bearing analogs were found to be potent inhibitors in our GroEL/ES-mediated substrate refolding assays; however, subsequent testing in the presence of an E. coli nitroreductase (NfsB) in situ indicated that metabolites of the nitrofuran-bearing analogs were potent GroEL/ES inhibitor pro-drugs. Consequently, this study has identified a new target of nitrofuran-containing drugs, and is the first reported instance of such a unique class of GroEL/ES chaperonin inhibitors. The intriguing results presented herein provide impetus for expanded studies to validate inhibitor mechanisms and optimize this antibacterial class using the respective GroEL/ES chaperonin systems and nitroreductases from E. coli and the ESKAPE bacteria.     

Comprehensive genomics and expression analysis of eceriferum (CER) genes in sunflower (Helianthus annuus)

Sunflower occupies the fourth position among oilseed crops the around the world. Eceriferum (CER) is an important gene family that plays critical role in very-long-chain fatty acids elongation and biosynthesis of epicuticular waxes under both biotic and abiotic stress conditions. The aim of present study was to investigate the effect of sunflower CER genes during drought stress condition. Thus, comparative analysis was undertaken for sunflower CER genes with Arabidopsis genome to determine phylogenetic relationship, chromosomal mapping, gene structures, gene ontology and conserved motifs. Furthermore, we subjected the sunflower cultivars under drought stress and used qRT-PCR analysis to explore the expression pattern of CER genes during drought conditions. We identified thirty-seven unevenly distributed CER genes in the sunflower genome. The phylogenetic analysis revealed that CER genes were grouped into seven clades in Arabidopsis, Helianthus annuus, and Gossypium hirsutum. Expression analysis showed that genes CER10 and CER60 were upregulated in sunflower during drought conditions, indicating that these genes are activated during drought stress. The results obtained will serve to characterize the CER gene family in sunflower and exploit the role of these genes in wax biosynthesis under limited water conditions.Key messageCuticular waxes protect the plants from drought stress, so we observed the expression of wax bio synthesis genes in recently sequences genome of Helianthus annuus. We observed that expression of wax biosynthesis genes CER10 and CER60 was upregulated when the plants were subjected to drought stress.     

Mendelian randomization analysis identified genes potentially pleiotropically associated with periodontitis

ObjectiveTo prioritize genes that were pleiotropically or potentially causally associated with periodontitis.MethodsWe applied the summary data-based Mendelian randomization (SMR) method integrating genome-wide association study (GWAS) for periodontitis and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with periodontitis. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analysis were done for participants of European and East Asian ancestries, separately.ResultsWe identified multiple genes showing pleiotropic association with periodontitis in participants of European ancestry and participants of East Asian ancestry. PDCD2 (corresponding probe: ILMN_1758915) was the top hit showing pleotropic association with periodontitis in the participants of European ancestry using CAGE eQTL data, and BX093763 (corresponding probe: ILMN_1899903) and AC104135.3 (corresponding probe: ENSG00000204792.2) were the top hits in the participants of East Asian ancestry using CAGE eQTL data and GTEx eQTL data, respectively.ConclusionWe identified multiple genes that may be involved in the pathogenesis of periodontitis in participants of European ancestry and participants of East Asian ancestry. Our findings provided important leads to a better understanding of the mechanisms underlying periodontitis and revealed potential therapeutic targets for the effective treatment of periodontitis.     

Heparanase promotes myeloma stemness and in vivo tumorigenesis

Heparanase is known to enhance the progression of many cancer types and is associated with poor patient prognosis. We recently reported that after patients with multiple myeloma were treated with high dose chemotherapy, the tumor cells that emerged upon relapse expressed a much higher level of heparanase than was present prior to therapy. Because tumor cells having stemness properties are thought to seed tumor relapse, we investigated whether heparanase had a role in promoting myeloma stemness. When plated at low density and grown in serum-free conditions that support survival and expansion of stem-like cells, myeloma cells expressing a low level of heparanase formed tumor spheroids poorly. In contrast, cells expressing a high level of heparanase formed significantly more and larger spheroids than did the heparanase low cells. Importantly, heparanase-low expressing cells exhibited plasticity and were induced to exhibit stemness properties when exposed to recombinant heparanase or to exosomes that contained a high level of heparanase cargo. The spheroid-forming heparanase-high cells had elevated expression of GLI1, SOX2 and ALDH1A1, three genes known to be associated with myeloma stemness. Inhibitors that block the heparan sulfate degrading activity of heparanase significantly diminished spheroid formation and expression of stemness genes implying a direct role of the enzyme in regulating stemness. Blocking the NF-κB pathway inhibited spheroid formation and expression of stemness genes demonstrating a role for NF-κB in heparanase-mediated stemness. Myeloma cells made deficient in heparanase exhibited decreased stemness properties in vitro and when injected into mice they formed tumors poorly compared to the robust tumorigenic capacity of cells expressing higher levels of heparanase. These studies reveal for the first time a role for heparanase in promoting cancer stemness and provide new insight into its function in driving tumor progression and its association with poor prognosis in cancer patients.     

The Endo-lysosomal System in Parkinson’s Disease: Expanding the Horizon

Parkinson’s disease (PD) is the second most common neurodegenerative disorder after Alzheimer’s disease, and its prevalence is increasing with age. A wealth of genetic evidence indicates that the endo-lysosomal system is a major pathway driving PD pathogenesis with a growing number of genes encoding endo-lysosomal proteins identified as risk factors for PD, making it a promising target for therapeutic intervention. However, detailed knowledge and understanding of the molecular mechanisms linking these genes to the disease are available for only a handful of them (e.g. LRRK2, GBA1, VPS35). Taking on the challenge of studying poorly characterized genes and proteins can be daunting, due to the limited availability of tools and knowledge from previous literature. This review aims at providing a valuable source of molecular and cellular insights into the biology of lesser-studied PD-linked endo-lysosomal genes, to help and encourage researchers in filling the knowledge gap around these less popular genetic players. Specific endo-lysosomal pathways discussed range from endocytosis, sorting, and vesicular trafficking to the regulation of membrane lipids of these membrane-bound organelles and the specific enzymatic activities they contain. We also provide perspectives on future challenges that the community needs to tackle and propose approaches to move forward in our understanding of these poorly studied endo-lysosomal genes. This will help harness their potential in designing innovative and efficient treatments to ultimately re-establish neuronal homeostasis in PD but also other diseases involving endo-lysosomal dysfunction.     

Exonic variants in multiple myeloma patients associated with relapsed/ refractory and response to bortezomib regimens

Novel treatment in multiple myeloma represented by proteasome inhibitors, immunomodulatory drugs and monoclonal antibodies have produced a deep response. However, relapses are possible, and all classes of drugs are refractory to patients. Next-generation sequencing has improved our understanding of the multiple myeloma genome related to drug resistance and has discovered many genomic variants. Therefore, this study was conducted to investigate new variants associated with drug resistance in MM patients who relapsed and refractory to bortezomib regimen and daratumumab treatment using next-generation sequencing for whole-exome sequencing. Peripheral blood samples were collected in EDTA tubes from six patients; four were in relapsed and refractory to bortezomib regimens and daratumumab; two patients responded to bortezomib regimens. Whole-exome sequencing was performed by the MGI-DNBSEQ-G400 instrument. We identified 21 variants in multiple myeloma patients. Seventeen variants were found in relapsed and refractory multiple myeloma in 11 genes (GNAQ, PMS1, CREB1, NSUNS2, PIK3CG, ROS1, PMS2, FIT4, KDM5A, STK11 and ZFHX3). And four variants were identified in two patients with response to bortezomib regimens in 4 genes (RAF1, CREB1, ZFHX3 and INSR). We have observed several genetic variants in many genes that may have been associated with the poor prognosis and poor response to treatment in these patients. These values should be further confirmed in large sample studies using the RNA-seq technique to identify genome expression.     

Construction and validation of a novel apoptosis-associated prognostic signature related to osteosarcoma metastasis and immune infiltration

BackgroundApoptosis played vital roles in the formation and progression of osteosarcoma. However, no studies elucidated the prognostic relationships between apoptosis-associated genes (AAGs) and osteosarcoma.MethodsThe differentially expressed genes associated with osteosarcoma metastasis and apoptosis were identified from GEO and MSigDB databases. The apoptosis-associated prognostic signature was established through univariate and multivariate cox regression analyses. The Kaplan–Meier (KM) survival curve, ROC curve and nomogram were constructed to investigate the predictive value of this signature. CIBERSORT algorithm and ssGSEA were used to explore the relationships between immune infiltration and AAG signature. The above results were validated in another GEO dataset and the expression of AAGs was also validated in osteosarcoma patient samples by immunohistochemistry.ResultsHSPB1 and IER3 were involved in AAG signature. In training and validation datasets, apoptosis-associated risk scores were negatively related to patient survival rates and the AAG signature was regarded as the independent prognostic factor. ROC and calibration curves demonstrated the signature and nomogram were reliable. GSEA revealed the signature related to immune-associated pathways. ssGSEA indicated that one immune cell and three immune functions were significantly dysregulated. The immunohistochemistry analyses of patients’ samples revealed that AAGs were significantly differently expressed between metastasis and non-metastasis osteosarcomas.ConclusionsThe present study identified and validated a novel apoptosis-associated prognostic signature related to osteosarcoma metastasis. It could serve as the potential biomarker and therapeutic targets for osteosarcoma in the future.     

Identification of chalcone synthase genes and their expression patterns reveal pollen abortion in cotton

Chalcone synthase (CHS) is a key enzyme and producing flavonoid derivatives as well play a vital roles in sustaining plant growth and development. However, the systematic and comprehensive analysis of CHS genes in island cotton (G. barbadense) has not been reported yet especially response to cytoplasmic male sterility (CMS). To fill this knowledge gap, a genome-wide investigation of CHS genes were studied in island cotton. A total of 20 GbCHS genes were identified and grouped into five GbCHSs. The gene structure analysis revealed that most of GbCHS genes consisted of two exons and one intron, and 20 motifs were identified. Twenty five pairs duplicated events (12 GbCHS genes) were identified including 23 segmental duplication pairs and two tandem duplication events, representing that GbCHS gene family amplification mainly owned to segmental duplication events and evolving slowly. Gene expression analysis exhibited that the GbCHS family genes presented a diversity expression patterns in various organs of cotton. Coupled with functional predictions and gene expression, the abnormal expression of GbCHS06, 10, 16 and 19 might be associated with pollen abortion of CMS line in island cotton. Conclusively, GbCHS genes exhibited diversity and conservation in many aspects, which will help to better understand functional studies and a reference for CHS research in island cotton and other plants.     

In-Depth Characterization of the Clostridioides difficile Phosphoproteome to Identify Ser/Thr Kinase Substrates

Clostridioides difficile is the leading cause of postantibiotic diarrhea in adults. During infection, the bacterium must rapidly adapt to the host environment by using survival strategies. Protein phosphorylation is a reversible post-translational modification employed ubiquitously for signal transduction and cellular regulation. Hanks-type serine/threonine kinases (STKs) and serine/threonine phosphatases have emerged as important players in bacterial cell signaling and pathogenicity. C. difficile encodes two STKs (PrkC and CD2148) and one phosphatase. We optimized a titanium dioxide phosphopeptide enrichment approach to determine the phosphoproteome of C. difficile. We identified and quantified 2500 proteins representing 63% of the theoretical proteome. To identify STK and serine/threonine phosphatase targets, we then performed comparative large-scale phosphoproteomics of the WT strain and isogenic ΔprkC, CD2148, Δstp, and prkC CD2148 mutants. We detected 635 proteins containing phosphorylated peptides. We showed that PrkC is phosphorylated on multiple sites in vivo and autophosphorylates in vitro. We were unable to detect a phosphorylation for CD2148 in vivo, whereas this kinase was phosphorylated in vitro only in the presence of PrkC. Forty-one phosphoproteins were identified as phosphorylated under the control of CD2148, whereas 114 proteins were phosphorylated under the control of PrkC including 27 phosphoproteins more phosphorylated in the ?stp mutant. We also observed enrichment for phosphothreonine among the phosphopeptides more phosphorylated in the Δstp mutant. Both kinases targeted pathways required for metabolism, translation, and stress response, whereas cell division and peptidoglycan metabolism were more specifically controlled by PrkC-dependent phosphorylation in agreement with the phenotypes of the ΔprkC mutant. Using a combination of approaches, we confirmed that FtsK was phosphorylated in vivo under the control of PrkC and that Spo0A was a substrate of PrkC in vitro. This study provides a detailed mapping of kinase–substrate relationships in C. difficile, paving the way for the identification of new biomarkers and therapeutic targets.