Phosphoproteome of Pristionchus pacificus provides insights into architecture of signaling networks in nematode models

Multi-Omics Analysis Revealed a Significant Alteration of Critical Metabolic Pathways Due to Sorafenib-Resistance in Hep3B Cell Lines

Hepatocellular carcinoma (HCC) is the second prominent cause of cancer-associated death worldwide. Usually, HCC is diagnosed in advanced stages, wherein sorafenib, a multiple target tyrosine kinase inhibitor, is used as the first line of treatment. Unfortunately, resistance to sorafenib is usually encountered within six months of treatment. Therefore, there is a critical need to identify the underlying reasons for drug resistance. In the present study, we investigated the proteomic and metabolomics alterations accompanying sorafenib resistance in hepatocellular carcinoma Hep3B cells by employing ultra-high-performance liquid chromatography quadrupole time of flight mass spectrometry (UHPLC-QTOF-MS). The Bruker Human Metabolome Database (HMDB) library was used to identify the differentially abundant metabolites through MetaboScape 4.0 software (Bruker). For protein annotation and identification, the Uniprot proteome for Homo sapiens (Human) database was utilized through MaxQuant. The results revealed that 27 metabolites and 18 proteins were significantly dysregulated due to sorafenib resistance in Hep3B cells compared to the parental phenotype. D-alanine, L-proline, o-tyrosine, succinic acid and phosphatidylcholine (PC, 16:0/16:0) were among the significantly altered metabolites. Ubiquitin carboxyl-terminal hydrolase isozyme L1, mitochondrial superoxide dismutase, UDP-glucose-6-dehydrogenase, sorbitol dehydrogenase and calpain small subunit 1 were among the significantly altered proteins. The findings revealed that resistant Hep3B cells demonstrated significant alterations in amino acid and nucleotide metabolic pathways, energy production pathways and other pathways related to cancer aggressiveness, such as migration, proliferation and drug-resistance. Joint pathway enrichment analysis unveiled unique pathways, including the antifolate resistance pathway and other important pathways that maintain cancer cells' survival, growth, and proliferation. Collectively, the results identified potential biomarkers for sorafenib-resistant HCC and gave insights into their role in chemotherapeutic drug resistance, cancer initiation, progression and aggressiveness, which may contribute to better prognosis and chemotherapeutic outcomes.

Comparative Phosphoproteomic Analysis of Sporulated Oocysts and Tachyzoites of Toxoplasma gondii Reveals Stage-Specific Patterns

Toxoplasma gondii is an obligate intracellular protozoan of severe threat to humans and livestock, whose life history harbors both gamic and apogamic stages. Chinese 1 (ToxoDB#9) was a preponderant genotype epidemic in food-derived animals and humans in China, with a different pathogenesis from the strains from the other nations of the world. Posttranslational modifications (PTMs) of proteins were critical mediators of the biology, developmental transforms, and pathogenesis of protozoan parasites. The phosphoprotein profiling and the difference between the developmental phases of T. gondii, contributing to development and infectivity, remain unknown. A quantitative phosphoproteomic approach using IBT integrated with TiO₂ affinity chromatography was applied to identify and analyze the difference in the phosphoproteomes between the sporulated oocysts and the tachyzoites of the virulent ToxoDB#9 (PYS) strain of T. gondii. A total of 4058 differential phosphopeptides, consisting of 2597 upregulated and 1461 downregulated phosphopeptides, were characterized between sporulated the oocysts and tachyzoites. Twenty-one motifs extracted from the upregulated phosphopeptides contained 19 serine motifs and 2 threonine motifs (GxxTP and TP), whereas 16 motifs identified from downregulated phosphopeptides included 13 serine motifs and 3 threonine motifs (KxxT, RxxT, and TP). Beyond the traditional kinases, some infrequent classes of kinases, including Ab1, EGFR, INSR, Jak, Src and Syk, were found to be corresponding to motifs from the upregulated and downregulated phosphopeptides. Remarkable functional properties of the differentially expressed phosphoproteins were discovered by GO analysis, KEGG pathway analysis, and STRING analysis. S8GFS8 (DNMT1-RFD domain-containing protein) and S8F5G5 (Histone kinase SNF1) were the two most connected peptides in the kinase-associated network. Out of these, phosphorylated modifications in histone kinase SNF1 have functioned in mitosis and interphase of T. gondii, as well as in the regulation of gene expression relevant to differentiation. Our study discovered a remarkable difference in the abundance of phosphopeptides between the sporulated oocysts and tachyzoites of the virulent ToxoDB#9 (PYS) strain of T. gondii, which may provide a new resource for understanding stage-specific differences in PTMs and may enhance the illustration of the regulatory mechanisms contributing to the development and infectivity of T. gondii.

Quantitative phosphoproteomics reveals the effect of baeSR and acrB genes on protein phosphorylation in Salmonella enterica serovar Typhimurium

The BaeSR two-component system and the AcrB efflux pump are closely associated with Salmonella resistance to antibiotics. However, the relationship between the two-component system, efflux pumps and protein phosphorylation of Salmonella is poorly understood. In this study, S. typhimurium ciprofloxacin-resistant strain CR, baeSR gene deletion strain CR^ΔbaeSR, acrB gene deletion strain CR^ΔacrB, and double gene deletion strain CR^{ΔbaeSRΔacrB} were used to explore phosphorylated proteins with significant difference, based on non-marker, quantitative phosphorylation modified proteomics technique. Consequently, 363 phosphosites of 213 phosphoproteins were identified in the four strains. More than 70% of the phosphosites were serine phosphorylation. In the CR^ΔbaeSR/CR, CR^ΔacrB/CR and CR^{ΔbaeSRΔacrB}/CR comparison groups, 36, 37 and 49 phosphosites were significantly altered, respectively. Bioinformatic analysis revealed that the main enrichment pathways of these differentially phosphorylated proteins were metabolic pathways, biosynthesis of antibiotics, phosphotransferase system (PTS), ABC transporters, and lipopolysaccharide biosynthesis. Furthermore, 21 differentially phosphorylated proteins were identified to be associated with antibiotic resistance. These results suggest that the BaeSR two-component system and the AcrB efflux pump affect the phosphorylation of proteins in S. typhimurium and may influence the drug resistance and virulence of S. typhimurium by affecting protein phosphorylation, providing a new idea to explore the mechanism of drug resistance in Salmonella.

Elucidating the molecular and developmental biology of parasitic nematodes: Moving to a multiomics paradigm

In the past two decades, significant progress has been made in the sequencing, assembly, annotation and analyses of genomes and transcriptomes of parasitic worms of socioeconomic importance. This progress has somewhat improved our knowledge and understanding of these pathogens at the molecular level. However, compared with the free-living nematode Caenorhabditis elegans, the areas of functional genomics, transcriptomics, proteomics and metabolomics of parasitic nematodes are still in their infancy, and there are major gaps in our knowledge and understanding of the molecular biology of parasitic nematodes. The information on signalling molecules, molecular pathways and microRNAs (miRNAs) that are known to be involved in developmental processes in C. elegans and the availability of some molecular resources (draft genomes, transcriptomes and some proteomes) for selected parasitic nematodes provide a basis to start exploring the developmental biology of parasitic nematodes. Indeed, some studies have identified molecules and pathways that might associate with developmental processes in related, parasitic nematodes, such as Haemonchus contortus (barber's pole worm). However, detailed information is often scant and 'omics resources are limited, preventing a proper integration of 'omic data sets and comprehensive analyses. Moreover, little is known about the functional roles of pheromones, hormones, signalling pathways and post-transcriptional/post-translational regulations in the development of key parasitic nematodes throughout their entire life cycles. Although C. elegans is an excellent model to assist molecular studies of parasitic nematodes, its use is limited when it comes to explorations of processes that are specific to parasitism within host animals. A deep understanding of parasitic nematodes, such as H. contortus, requires substantially enhanced resources and the use of integrative 'omics approaches for analyses. The improved genome and well-established in vitro larval culture system for H. contortus provide unprecedented opportunities for comprehensive studies of the transcriptomes (mRNA and miRNA), proteomes (somatic, excretory/secretory and phosphorylated proteins) and lipidomes (e.g., polar and neutral lipids) of this nematode. Such resources should enable in-depth explorations of its developmental biology at a level, not previously possible. The main aims of this review are (i) to provide a background on the development of nematodes, with a particular emphasis on the molecular aspects involved in the dauer formation and exit in C. elegans; (ii) to critically appraise the current state of knowledge of the developmental biology of parasitic nematodes and identify key knowledge gaps; (iii) to cover salient aspects of H. contortus, with a focus on the recent advances in genomics, transcriptomics, proteomics and lipidomics as well as in vitro culturing systems; (iv) to review recent advances in our knowledge and understanding of the molecular and developmental biology of H. contortus using an integrative multiomics approach, and discuss the implications of this approach for detailed explorations of signalling molecules, molecular processes and pathways likely associated with nematode development, adaptation and parasitism, and for the identification of novel intervention targets against these pathogens. Clearly, the multiomics approach established recently is readily applicable to exploring a wide range of interesting and socioeconomically significant parasitic worms (including also trematodes and cestodes) at the molecular level, and to elucidate host-parasite interactions and disease processes.

The developmental phosphoproteome of Haemonchus contortus

Protein phosphorylation plays essential roles in many cellular processes. Despite recent progress in the genomics, transcriptomics and proteomics of socioeconomically important parasitic nematodes, there is scant phosphoproteomic data to underpin molecular biological discovery. Here, using the phosphopeptide enrichment-based LC-MS/MS and data-independent acquisition (DIA) quantitation, we characterised the first developmental phosphoproteome of the parasitic nematode Haemonchus contortus - one of the most pathogenic parasites of ruminant livestock. Totally, 1804 phosphorylated proteins with 4406 phosphorylation sites ('phosphosites') from different developmental stages/sexes were identified. Bioinformatic analyses of quantified 'phosphosites' exhibited distinctive stage- and sex-specific patterns during development, and identified a subset of phosphoproteins proposed to play crucial roles in processes such as spindle positioning, signal transduction and kinase activity. A sequence-based comparison of the phosphoproteome of H. contortus with those of two free-living nematode species (Caenorhabditis elegans and Pristionchus pacificus) suggested a limited number of common protein phosphorylation events among these species. Our findings infer active roles for protein phosphorylation in the adaptation of a parasitic nematode to a constantly changing external environment. The phosphoproteomic data set for H. contortus provides a basis to better understand phosphorylation and associated biological processes (e.g., regulation of signal transduction), and might enable the discovery of novel anthelmintic targets. SIGNIFICANCE: Here, we report the first phosphoproteome for a socioeconomically parasitic nematode (Haemonchus contortus). This phosphoproteome exhibits distinctive patterns during development, suggesting active roles of post-translational modification in the parasite's adaptation to changing environments within and outside of the host animal. This work sheds a light on the developmental phosphorylation in a parasitic nematode, and could enable the discovery of novel interventions against major pathogens.

iTRAQ-Based Global Phosphoproteomics Reveals Novel Molecular Differences Between Toxoplasma gondii Strains of Different Genotypes

To gain insights into differences in the virulence among T. gondii strains at the post-translational level, we conducted a quantitative analysis of the phosphoproteome profile of T. gondii strains belonging to three different genotypes. Phosphopeptides from three strains, type I (RH strain), type II (PRU strain) and ToxoDB#9 (PYS strain), were enriched by titanium dioxide (TiO2) affinity chromatography and quantified using iTRAQ technology. A total of 1,441 phosphopeptides, 1,250 phosphorylation sites and 759 phosphoproteins were detected. In addition, 392, 298, and 436 differentially expressed phosphoproteins (DEPs) were identified in RH strain when comparing RH/PRU strains, in PRU strain when comparing PRU/PYS strains, and in PYS strain when comparing PYS/RH strains, respectively. Functional characterization of the DEPs using GO, KEGG, and STRING analyses revealed marked differences between the three strains. In silico kinase substrate motif analysis of the DEPs revealed three (RxxS, SxxE, and SxxxE), three (RxxS, SxxE, and SP), and five (SxxE, SP, SxE, LxRxxS, and RxxS) motifs in RH strain when comparing RH/PRU strains, in PRU strain when comparing PRU/PYS, and in PYS strain when comparing PYS/RH strains, respectively. This suggests that multiple overrepresented protein kinases including PKA, PKG, CKII, IKK, and MAPK could be involved in such a difference between T. gondii strains. Kinase associated network analysis showed that ROP5, ROP16, and cell-cycle-associated protein kinase CDK were the most connected kinase peptides. Our data reveal significant changes in the abundance of phosphoproteins between T. gondii genotypes, which explain some of the mechanisms that contribute to the virulence heterogeneity of this parasite.

Dauer signalling pathway model for Haemonchus contortus

BACKGROUND

Signalling pathways have been extensively investigated in the free-living nematode Caenorhabditis elegans, but very little is known about these pathways in parasitic nematodes. Here, we constructed a model for the dauer-associated signalling pathways in an economically highly significant parasitic worm, Haemonchus contortus.

METHODS

Guided by data and information available for C. elegans, we used extensive genomic and transcriptomic datasets to infer gene homologues in the dauer-associated pathways, explore developmental transcriptomic, proteomic and phosphoproteomic profiles in H. contortus and study selected molecular structures.

RESULTS

The canonical cyclic guanosine monophosphate (cGMP), transforming growth factor-β (TGF-β), insulin-like growth factor 1 (IGF-1) and steroid hormone signalling pathways of H. contortus were inferred to represent a total of 61 gene homologues. Compared with C. elegans, H. contortus has a reduced set of genes encoding insulin-like peptides, implying evolutionary and biological divergences between the parasitic and free-living nematodes. Similar transcription profiles were found for all gene homologues between the infective stage of H. contortus and dauer stage of C. elegans. High transcriptional levels for genes encoding G protein-coupled receptors (GPCRs), TGF-β, insulin-like ligands (e.g. ins-1, ins-17 and ins-18) and transcriptional factors (e.g. daf-16) in the infective L3 stage of H. contortus were suggestive of critical functional roles in this stage. Conspicuous protein expression patterns and extensive phosphorylation of some components of these pathways suggested marked post-translational modifications also in the L3 stage. The high structural similarity in the DAF-12 ligand binding domain among nematodes indicated functional conservation in steroid (i.e. dafachronic acid) signalling linked to worm development.

CONCLUSIONS

Taken together, this pathway model provides a basis to explore hypotheses regarding biological processes and regulatory mechanisms (via particular microRNAs, phosphorylation events and/or lipids) associated with the development of H. contortus and related nematodes as well as parasite-host cross talk, which could aid the discovery of new therapeutic targets.

Advances in kinome research of parasitic worms - implications for fundamental research and applied biotechnological outcomes

Protein kinases are enzymes that play essential roles in the regulation of many cellular processes. Despite expansions in the fields of genomics, transcriptomics and bioinformatics, there is limited information on the kinase complements (kinomes) of most eukaryotic organisms, including parasitic worms that cause serious diseases of humans and animals. The biological uniqueness of these worms and the draft status of their genomes pose challenges for the identification and classification of protein kinases using established tools. In this article, we provide an account of kinase biology, the roles of kinases in diseases and their importance as drug targets, and drug discovery efforts in key socioeconomically important parasitic worms. In this context, we summarise methods and resources commonly used for the curation, identification, classification and functional annotation of protein kinase sequences from draft genomes; review recent advances made in the characterisation of the worm kinomes; and discuss the implications of these advances for investigating kinase signalling and developing small-molecule inhibitors as new anti-parasitic drugs.

Comparative Genomics of Gene Loss and Gain in Caenorhabditis and Other Nematodes

Nematodes, such as Caenorhabditis elegans, form one of the most species-rich animal phyla. By now more than 30 nematode genomes have been published allowing for comparative genomic analyses at various different time-scales. The majority of a nematode's gene repertoire is represented by either duplicated or so-called orphan genes of unknown origin. This indicates the importance of mechanisms that generate new genes during the course of evolution. While it is certain that nematodes have acquired genes by horizontal gene transfer from various donors, this process only explains a small portion of the nematode gene content. As evolutionary genomic analyses strongly support that most orphan genes are indeed protein-coding, future studies will have to decide, whether they are result from extreme divergence or evolved de novo from previously noncoding sequences. In this contribution, I summarize several studies investigating gene loss and gain in nematodes and discuss the strengths and weaknesses of individual approaches and datasets. These approaches can be used to ask nematode-specific questions such as associated with the evolution of parasitism or with switches in mating systems, but also can complement studies in other animal phyla like vertebrates and insects to broaden our general view on genome evolution.

A temporal proteome dynamics study reveals the molecular basis of induced phenotypic resistance in Mycobacterium smegmatis at sub-lethal rifampicin concentrations

In the last 40 years only one new antitubercular drug has been approved, whilst resistance to current drugs, including rifampicin, is spreading. Here, we used the model organism Mycobacterium smegmatis to study mechanisms of phenotypic mycobacterial resistance, employing quantitative mass spectrometry-based proteomics to investigate the temporal effects of sub-lethal concentrations of rifampicin on the mycobacterial proteome at time-points corresponding to early response, onset of bacteriostasis and early recovery. Across 18 samples, a total of 3,218 proteins were identified from 31,846 distinct peptides averaging 16,250 identified peptides per sample. We found evidence that two component signal transduction systems (e.g. MprA/MprB) play a major role during initial mycobacterial adaptive responses to sub-lethal rifampicin and that, after dampening an initial SOS response, the bacteria supress the DevR (DosR) regulon and also upregulate their transcriptional and translational machineries. Furthermore, we found a co-ordinated dysregulation in haeme and mycobactin synthesis. Finally, gradual upregulation of the M. smegmatis-specific rifampin ADP-ribosyl transferase was observed which, together with upregulation of transcriptional and translational machinery, likely explains recovery of normal growth. Overall, our data indicates that in mycobacteria, sub-lethal rifampicin triggers a concerted phenotypic response that contrasts significantly with that observed at higher antimicrobial doses.

First insights into the nature and evolution of antisense transcription in nematodes

BACKGROUND

The development of multicellular organisms is coordinated by various gene regulatory mechanisms that ensure correct spatio-temporal patterns of gene expression. Recently, the role of antisense transcription in gene regulation has moved into focus of research. To characterize genome-wide patterns of antisense transcription and to study their evolutionary conservation, we sequenced a strand-specific RNA-seq library of the nematode Pristionchus pacificus.

RESULTS

We identified 1112 antisense configurations of which the largest group represents 465 antisense transcripts (ASTs) that are fully embedded in introns of their host genes. We find that most ASTs show homology to protein-coding genes and are overrepresented in proteomic data. Together with the finding, that expression levels of ASTs and host genes are uncorrelated, this indicates that most ASTs in P. pacificus do not represent non-coding RNAs and do not exhibit regulatory functions on their host genes. We studied the evolution of antisense gene pairs across 20 nematode genomes, showing that the majority of pairs is lineage-specific and even the highly conserved vps-4, ddx-27, and sel-2 loci show abundant structural changes including duplications, deletions, intron gains and loss of antisense transcription. In contrast, host genes in general, are remarkably conserved and encode exceptionally long introns leading to unusually large blocks of conserved synteny.

CONCLUSIONS

Our study has shown that in P. pacificus antisense transcription as such does not define non-coding RNAs but is rather a feature of highly conserved genes with long introns. We hypothesize that the presence of regulatory elements imposes evolutionary constraint on the intron length, but simultaneously, their large size makes them a likely target for translocation of genomic elements including protein-coding genes that eventually end up as ASTs.

Are orphan genes protein-coding, prediction artifacts, or non-coding RNAs?

Background

Current genome sequencing projects reveal substantial numbers of taxonomically restricted, so called orphan genes that lack homology with genes from other evolutionary lineages. However, it is not clear to what extent orphan genes are real, genomic artifacts, or represent non-coding RNAs.

Results

Here, we use a simple set of assumptions to test the nature of orphan genes. First, a sequence that is transcribed is considered a real biological entity. Second, every sequence that is supported by proteome data or shows a depletion of non-synonymous substitutions is a protein-coding gene. Using genomic, transcriptomic and proteomic data for the nematode Pristionchus pacificus, we show that between 4129–7997 (42–81 %) of predicted orphan genes are expressed and 3818–7545 (39–76 %) of orphan genes are under negative selection. In three cases that exhibited strong evolutionary constraint but lacked expression evidence in 14 RNA-seq samples, we could experimentally validate the predicted gene structures. Comparing different data sets to infer selection on orphan gene clusters, we find that the presence of a closely related genome provides the most powerful resource to robustly identify evidence of negative selection. However, even in the absence of other genomic data, the availability of paralogous sequences was enough to show negative selection in 8–10 % of orphan genes.

Conclusions

Our study shows that the great majority of previously identified orphan genes in P. pacificus are indeed protein-coding genes. Even though this work represents a case study on a single species, our approach can be transferred to genomic data of other non-model organisms in order to ascertain the protein-coding nature of orphan genes.