Systematic analysis of the phosphoproteome and kinase-substrate networks in the mouse testis

A Review of Machine Learning and Algorithmic Methods for Protein Phosphorylation Site Prediction

Post-translational modifications (PTMs) have key roles in extending the functional diversity of proteins and, as a result, regulating diverse cellular processes in prokaryotic and eukaryotic organisms. Phosphorylation modification is a vital PTM that occurs in most proteins and plays a significant role in many biological processes. Disorders in the phosphorylation process lead to multiple diseases, including neurological disorders and cancers. The purpose of this review is to organize this body of knowledge associated with phosphorylation site (p-site) prediction to facilitate future research in this field. At first, we comprehensively review all related databases and introduce all steps regarding dataset creation, data preprocessing, and method evaluation in p-site prediction. Next, we investigate p-site prediction methods, which are divided into two computational groups: algorithmic and machine learning (ML). Additionally, it is shown that there are basically two main approaches for p-site prediction by ML: conventional and end-to-end deep learning methods, both of which are given an overview. Moreover, this review introduces the most important feature extraction techniques, which have mostly been used in p-site prediction. Finally, we create three test sets from new proteins related to the released version of the database of protein post-translational modifications (dbPTM) in 2022 based on general and human species. Evaluating online p-site prediction tools on newly added proteins introduced in the dbPTM 2022 release, distinct from those in the dbPTM 2019 release, reveals their limitations. In other words, the actual performance of these online p-site prediction tools on unseen proteins is notably lower than the results reported in their respective research papers.

STK33 phosphorylates fibrous sheath protein AKAP3/4 to regulate sperm flagella assembly in spermiogenesis

Spermatogenesis defects are important for male infertility; however, the etiology and pathogenesis are still unknown. Herein, we identified two loss of function mutations of STK33 in 7 individuals with non-obstructive azoospermia. Further functional studies of these frameshift and nonsense mutations revealed that Stk33^-/KI male mice were sterile, and Stk33^-/KI sperm were abnormal with defects in the mitochondrial sheath (MS), fibrous sheath (FS), outer dense fiber (ODF) and axoneme. Stk33^KI/KI male mice were subfertile and had oligoasthenozoospermia. Differential phosphoproteomic analysis and in vitro kinase assay identified novel phosphorylation substrates of STK33, fibrous sheath components AKAP3 and AKAP4, whose expression levels decreased in testis after deletion of Stk33. STK33 regulated the phosphorylation of AKAP3/4, affected the assembly of fibrous sheath in the sperm, and played an essential role in spermiogenesis and male infertility.

Deep Learning in Phosphoproteomics: Methods and Application in Cancer Drug Discovery

Protein phosphorylation is a key post-translational modification (PTM) that is a central regulatory mechanism of many cellular signaling pathways. Several protein kinases and phosphatases precisely control this biochemical process. Defects in the functions of these proteins have been implicated in many diseases, including cancer. Mass spectrometry (MS)-based analysis of biological samples provides in-depth coverage of phosphoproteome. A large amount of MS data available in public repositories has unveiled big data in the field of phosphoproteomics. To address the challenges associated with handling large data and expanding confidence in phosphorylation site prediction, the development of many computational algorithms and machine learning-based approaches have gained momentum in recent years. Together, the emergence of experimental methods with high resolution and sensitivity and data mining algorithms has provided robust analytical platforms for quantitative proteomics. In this review, we compile a comprehensive collection of bioinformatic resources used for the prediction of phosphorylation sites, and their potential therapeutic applications in the context of cancer.

Integrated analysis of SR-like protein kinases Sky1 and Sky2 links signaling networks with transcriptional regulation in Candida albicans

Fungal infections are a major global health burden where Candida albicans is among the most common fungal pathogen in humans and is a common cause of invasive candidiasis. Fungal phenotypes, such as those related to morphology, proliferation and virulence are mainly driven by gene expression, which is primarily regulated by kinase signaling cascades. Serine-arginine (SR) protein kinases are highly conserved among eukaryotes and are involved in major transcriptional processes in human and S. cerevisiae. Candida albicans harbors two SR protein kinases, while Sky2 is important for metabolic adaptation, Sky1 has similar functions as in S. cerevisiae. To investigate the role of these SR kinases for the regulation of transcriptional responses in C. albicans, we performed RNA sequencing of sky1Δ and sky2Δ and integrated a comprehensive phosphoproteome dataset of these mutants. Using a Systems Biology approach, we study transcriptional regulation in the context of kinase signaling networks. Transcriptomic enrichment analysis indicates that pathways involved in the regulation of gene expression are downregulated and mitochondrial processes are upregulated in sky1Δ. In sky2Δ, primarily metabolic processes are affected, especially for arginine, and we observed that arginine-induced hyphae formation is impaired in sky2Δ. In addition, our analysis identifies several transcription factors as potential drivers of the transcriptional response. Among these, a core set is shared between both kinase knockouts, but it appears to regulate different subsets of target genes. To elucidate these diverse regulatory patterns, we created network modules by integrating the data of site-specific protein phosphorylation and gene expression with kinase-substrate predictions and protein-protein interactions. These integrated signaling modules reveal shared parts but also highlight specific patterns characteristic for each kinase. Interestingly, the modules contain many proteins involved in fungal morphogenesis and stress response. Accordingly, experimental phenotyping shows a higher resistance to Hygromycin B for sky1Δ. Thus, our study demonstrates that a combination of computational approaches with integration of experimental data can offer a new systems biological perspective on the complex network of signaling and transcription. With that, the investigation of the interface between signaling and transcriptional regulation in C. albicans provides a deeper insight into how cellular mechanisms can shape the phenotype.

CDK6 activity in a recurring convergent kinase network motif

In humans, more than 500 kinases phosphorylate ~15% of all proteins in an emerging phosphorylation network. Convergent local interaction motifs, in which ≥two kinases phosphorylate the same substrate, underlie feedback loops and signal amplification events but have not been systematically analyzed. Here, we first report a network-wide computational analysis of convergent kinase-substrate relationships (cKSRs). In experimentally validated phosphorylation sites, we find that cKSRs are common and involve >80% of all human kinases and >24% of all substrates. We show that cKSRs occur over a wide range of stoichiometries, in many instances harnessing co-expressed kinases from family subgroups. We then experimentally demonstrate for the prototypical convergent CDK4/6 kinase pair how multiple inputs phosphorylate the tumor suppressor retinoblastoma protein (RB) and thereby hamper in situ analysis of the individual kinases. We hypothesize that overexpression of one kinase combined with a CDK4/6 inhibitor can dissect convergence. In breast cancer cells expressing high levels of CDK4, we confirm this hypothesis and develop a high-throughput compatible assay that quantifies genetically modified CDK6 variants and inhibitors. Collectively, our work reveals the occurrence, topology, and experimental dissection of convergent interactions toward a deeper understanding of kinase networks and functions.

Computational systems approach towards phosphoproteomics and their downstream regulation

Protein phosphorylation plays an essential role in modulating cell signalling and its downstream transcriptional and translational regulations. Until recently, protein phosphorylation has been studied mostly using low-throughput biochemical assays. The advancement of mass spectrometry (MS)-based phosphoproteomics transformed the field by enabling measurement of proteome-wide phosphorylation events, where tens of thousands of phosphosites are routinely identified and quantified in an experiment. This has brought a significant challenge in analysing large-scale phosphoproteomic data, making computational methods and systems approaches integral parts of phosphoproteomics. Previous works have primarily focused on reviewing the experimental techniques in MS-based phosphoproteomics, yet a systematic survey of the computational landscape in this field is still missing. Here, we review computational methods and tools, and systems approaches that have been developed for phosphoproteomics data analysis. We categorise them into four aspects including data processing, functional analysis, phosphoproteome annotation and their integration with other omics, and in each aspect, we discuss the key methods and example studies. Lastly, we highlight some of the potential research directions on which future work would make a significant contribution to this fast-growing field. We hope this review provides a useful snapshot of the field of computational systems phosphoproteomics and stimulates new research that drives future development.

Using transcriptomics, proteomics and phosphoproteomics as new approach methodology (NAM) to define biological responses for chemical safety assessment

Omic-based technologies are of particular interest and importance for hazard identification and health risk characterization of chemicals. Their application in the new approach methodologies (NAMs) anchored on cellular toxicity pathways is based on the premise that any apical health endpoint change must be underpinned by some alterations at the omic levels. In the present study we examined the cellular responses to two chemicals, caffeine and coumarin, by generating and integrating multi-omic data from multi-dose and multi-time point transcriptomic, proteomic and phosphoproteomic experiments. We showed that the methodology presented here was able to capture the complete chain of events from the first chemical-induced changes at the phosphoproteome level, to changes in gene expression, and lastly to changes in protein abundance, each with vastly different points of departure (PODs). In HepG2 cells we found that the metabolism of lipids and general cellular stress response to be the dominant biological processes in response to caffeine and coumarin exposure, respectively. The phosphoproteomic changes were detected early in time, at very low doses and provided a fast, adaptive cellular response to chemical exposure with 7-37-fold lower points of departure comparing to the transcriptomics. Changes in protein abundance were found much less frequently than transcriptomic changes. While challenges remain, our study provides strong and novel evidence supporting the notion that these three omic technologies can be used in an integrated manner to facilitate a more complete understanding of pathway perturbations and POD determinations for risk assessment of chemical exposures.

Global phosphoproteomic analysis identified key kinases regulating male meiosis in mouse

Meiosis, a highly conserved process in organisms from fungi to mammals, is subjected to protein phosphorylation regulation. Due to the low abundance of phosphorylation, there is a lack of systemic characterization of phosphorylation regulation of meiosis in mammals. Using the phosphoproteomic approach, we profiled large-scale phosphoproteome of purified primary spermatocytes undergoing meiosis I, and identified 14,660 phosphorylation sites in 4419 phosphoproteins. Kinase-substrate phosphorylation network analysis followed by in vitro meiosis study showed that CDK9 was essential for meiosis progression to metaphase I and had enriched substrate phosphorylation sites in proteins involved in meiotic cell cycle. In addition, histones and epigenetic factors were found to be widely phosphorylated. Among those, HASPIN was found to be essential for male fertility. Haspin knockout led to misalignment of chromosomes, apoptosis of metaphase spermatocytes and a decreased number of sperm by deregulation of H3T3ph, chromosomal passenger complex (CPC) and spindle assembly checkpoint (SAC). The complicated protein phosphorylation and its important regulatory functions in meiosis indicated that in-depth studies of phosphorylation-mediated signaling could help us elucidate the mechanisms of meiosis.

Homeodomain-interacting protein kinase HIPK4 regulates phosphorylation of manchette protein RIMBP3 during spermiogenesis

Nonobstructive azoospermia (NOA) is the most serious form of spermatogenesis abnormalities in male infertility. Genetic factors are important to consider as elements leading to NOA. Although many pathogenic genes have been reported, the causative genes of NOA for many patients are still unknown. In this study, we found ten point mutations in the gene encoding homeodomain-interacting protein kinase 4 (HIPK4) in patients with NOA, and using in vitro studies, we determined a premature termination point mutation (p. Lys490∗, c.1468A>T) that can cause decreased expression of HIPK4. Our phosphoproteomic analysis of Hipk4^−/− testes revealed phosphorylation of multiple proteins regulated by HIPK4 during spermiogenesis. We also confirmed that a substrate of HIPK4 with four downregulated phosphorylation sites matching the xSPx motif is the known manchette-related protein RIMS-binding protein 3, which is required for sperm head morphogenesis. Therefore, we conclude HIPK4 regulates the phosphorylation of manchette protein RIMS-binding protein 3 and plays essential roles in sperm head shaping and male fertility.

Sperm Phosphoproteome: Unraveling Male Infertility

Infertility affects approximately 15% of couples worldwide of childbearing age, and in many cases the etiology of male infertility is unknown. The current standard evaluation of semen is insufficient to establish an accurate diagnosis. Proteomics techniques, such as phosphoproteomics, applied in this field are a powerful tool to understand the mechanisms that regulate sperm functions such as motility, which is essential for successful fertilization. Among the post-translational modifications of sperm proteins, this review summarizes, from a proteomic perspective, the updated knowledge of protein phosphorylation, in human spermatozoa, as a relevant molecular mechanism involved in the regulation of sperm physiology. Specifically, the role of sperm protein phosphorylation in motility and, consequently, in sperm quality is highlighted. Additionally, through the analysis of published comparative phosphoproteomic studies, some candidate human sperm phosphoproteins associated with low sperm motility are proposed. Despite the remarkable advances in phosphoproteomics technologies, the relatively low number of studies performed in human spermatozoa suggests that phosphoproteomics has not been applied to its full potential in studying male infertility yet. Therefore, further studies will improve the application of this procedure and overcome the limitations, increasing the understanding of regulatory mechanisms underlying protein phosphorylation in sperm motility and, consequently, in male fertility.

Recent progress of proteomic analysis on spermatogenesis

Testis, the only organ responsible for generating sperm, is by far the organ with the largest variety of proteins and tissue-specific proteins in humans. In testis, spermatogenesis is a multi-step complex process well-accepted that protein and mRNA are decoupled in certain stages of spermatogenesis. With the fast development of mass spectrometry-based proteomics, it is possible to systemically study protein abundances and modifications in testis and sperm to help us understand the molecular mechanisms of spermatogenesis. This review provides an overview of the recent progress of proteomics analysis on spermatogenesis, including protein expression and multiple PTMs, such as phosphorylation, glycosylation, ubiquitylation, and acetylation.

Sperm acquire epididymis-derived proteins through epididymosomes

STUDY QUESTION

Are epididymosomes implicated in protein transfer from the epididymis to spermatozoa?

SUMMARY ANSWER

We characterized the contribution of epididymal secretions to the sperm proteome and demonstrated that sperm acquire epididymal proteins through epididymosomes.

WHAT IS KNOWN ALREADY

Testicular sperm are immature cells unable to fertilize an oocyte. After leaving the testis, sperm transit along the epididymis to acquire motility and fertilizing abilities. It is well known that marked changes in the sperm proteome profile occur during epididymal maturation. Since the sperm is a transcriptional and translational inert cell, previous studies have shown that sperm incorporate proteins, RNA and lipids from extracellular vesicles (EVs), released by epithelial cells lining the male reproductive tract.

STUDY DESIGN, SIZE, DURATION

We examined the contribution of the epididymis to the post-testicular maturation of spermatozoa, via the production of EVs named epididymosomes, released by epididymal epithelial cells. An integrative analysis using both human and mouse data was performed to identify sperm proteins with a potential epididymis-derived origin. Testes and epididymides from adult humans (n = 9) and adult mice (n = 3) were used to experimentally validate the tissue localization of four selected proteins using high-resolution confocal microscopy. Mouse epididymal sperm were co-incubated with carboxyfluorescein succinimidyl ester (CFSE)-labeled epididymosomes (n = 4 mice), and visualized using high-resolution confocal microscopy.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Adult (12-week-old) C57BL/CBAF1 wild-type male mice and adult humans were used for validation purposes. Testes and epididymides from both mice and humans were obtained and processed for immunofluorescence. Mouse epididymal sperm and mouse epididymosomes were obtained from the epididymal cauda segment. Fluorescent epididymosomes were obtained after labeling the epididymal vesicles with CFSE dye followed by epididymosome isolation using a density cushion. Immunofluorescence was performed following co-incubation of sperm with epididymosomes in vitro. High-resolution confocal microscopy and 3D image reconstruction were used to visualize protein localization and sperm-epididymosomes interactions.

MAIN RESULTS AND THE ROLE OF CHANCE

Through in silico analysis, we first identified 25 sperm proteins with a putative epididymal origin that were conserved in both human and mouse spermatozoa. From those, the epididymal origin of four sperm proteins (SLC27A2, EDDM3B, KRT19 and WFDC8) was validated by high-resolution confocal microscopy. SLC27A2, EDDM3B, KRT19 and WFDC8 were all detected in epithelial cells lining the human and mouse epididymis, and absent from human and mouse seminiferous tubules. We found region-specific expression patterns of these proteins throughout the mouse epididymides. In addition, while EDDM3B, KRT19 and WFDC8 were detected in both epididymal principal and clear cells (CCs), SLC27A2 was exclusively expressed in CCs. Finally, we showed that CFSE-fluorescently labeled epididymosomes interact with sperm in vitro and about 12-36% of the epididymosomes contain the targeted sperm proteins with an epididymal origin.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

The human and mouse sample size was limited and our results were descriptive. The analyses of epididymal sperm and epididymosomes were solely performed in the mouse model due to the difficulties in obtaining epididymal luminal fluid human samples. Alternatively, human ejaculated sperm and seminal EVs could not be used because ejaculated sperm have already contacted with the fluids secreted by the male accessory sex glands, and seminal EVs contain other EVs in addition to epididymosomes, such as the abundant prostate-derived EVs.

WIDER IMPLICATIONS OF THE FINDINGS

Our findings indicate that epididymosomes are capable of providing spermatozoa with a new set of epididymis-derived proteins that could modulate the sperm proteome and, subsequently, participate in the post-testicular maturation of sperm cells. Additionally, our data provide further evidence of the novel role of epididymal CCs in epididymosome production. Identifying mechanisms by which sperm mature to acquire their fertilization potential would, ultimately, lead to a better understanding of male reproductive health and may help to identify potential therapeutic strategies to improve male infertility.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Spanish Ministry of Economy and Competitiveness (Ministerio de Economía y Competividad; fondos FEDER 'una manera de hacer Europa' PI13/00699 and PI16/00346 to R.O.; and Sara Borrell Postdoctoral Fellowship, Acción Estratégica en Salud, CD17/00109 to J.C.), by National Institutes of Health (grants HD040793 and HD069623 to S.B., grant HD104672-01 to M.A.B.), by the Spanish Ministry of Education, Culture and Sports (Ministerio de Educación, Cultura y Deporte para la Formación de Profesorado Universitario, FPU15/02306 to F.B.), by a Lalor Foundation Fellowship (to F.B. and M.A.B.), by the Government of Catalonia (Generalitat de Catalunya, pla estratègic de recerca i innovació en salut, PERIS 2016-2020, SLT002/16/00337 to M.J.), by Fundació Universitària Agustí Pedro i Pons (to F.B.), and by the American Society for Biochemistry and Molecular Biology (PROLAB Award from ASBMB/IUBMB/PABMB to F.B.). Confocal microscopy and transmission electron microscopy was performed in the Microscopy Core facility of the Massachusetts General Hospital (MGH) Center for Systems Biology/Program in Membrane Biology which receives support from Boston Area Diabetes and Endocrinology Research Center (BADERC) award DK57521 and Center for the Study of Inflammatory Bowel Disease grant DK43351. The Zeiss LSM800 microscope was acquired using an NIH Shared Instrumentation Grant S10-OD-021577-01. The authors have no conflicts of interest to declare.

Phosphoproteomics of ATR signaling in mouse testes

The phosphatidylinositol 3' kinase (PI3K)-related kinase ATR is crucial for mammalian meiosis. ATR promotes meiotic progression by coordinating key events in DNA repair, meiotic sex chromosome inactivation (MSCI), and checkpoint-dependent quality control during meiotic prophase I. Despite its central roles in meiosis, the ATR-dependent meiotic signaling network remains largely unknown. Here, we used phosphoproteomics to define ATR signaling events in testes from mice following chemical and genetic ablation of ATR signaling. Quantitative analysis of phosphoproteomes obtained after germ cell-specific genetic ablation of the ATR activating 9-1-1 complex or treatment with ATR inhibitor identified over 14,000 phosphorylation sites from testes samples, of which 401 phosphorylation sites were found to be dependent on both the 9-1-1 complex and ATR. Our analyses identified ATR-dependent phosphorylation events in crucial DNA damage signaling and DNA repair proteins including TOPBP1, SMC3, MDC1, RAD50, and SLX4. Importantly, we identified ATR and RAD1-dependent phosphorylation events in proteins involved in mRNA regulatory processes, including SETX and RANBP3, whose localization to the sex body was lost upon ATR inhibition. In addition to identifying the expected ATR-targeted S/T-Q motif, we identified enrichment of an S/T-P-X-K motif in the set of ATR-dependent events, suggesting that ATR promotes signaling via proline-directed kinase(s) during meiosis. Indeed, we found that ATR signaling is important for the proper localization of CDK2 in spermatocytes. Overall, our analysis establishes a map of ATR signaling in mouse testes and highlights potential meiotic-specific actions of ATR during prophase I progression.

Cyclin-dependent kinase 7 is essential for spermatogenesis by regulating retinoic acid signaling pathways and the STAT3 molecular pathway

Spermatogenesis is a complex process that requires precise regulation. Phosphorylation plays a role in spermatogenesis by regulating protein structure and activity. This study focused on cyclin-dependent kinase 7 (CDK7), and explored its function and molecular mechanisms in spermatogenesis in vitro in a cell line and in vivo in a mouse model. Inhibition of CDK7 activity affected spermatogonia proliferation and differentiation, and we found that CDK7 regulates retinoic acid (RA)-mediated c-KIT expression to play a role in spermatogonia. Then, we demonstrated that inhibition of CDK7 affected meiosis initiation, DNA repair, and synaptonemal complex formation in meiosis progression, and CDK7 played this role by regulating RA-mediated STRA8 and REC8 signaling pathways. Moreover, inhibition of CDK7 impacted spermatid differentiation and resulted in decreased counts, decreased motility, and increased head deformity of sperm. We demonstrated that CDK7 affects germ cell apoptosis and sperm motility by activating STAT3 and that STAT3 further regulates Cortactin expression to influence the nuclear elongation, chromatin condensation, and acrosome formation of sperm. Additionally, EP300 was identified as another potential target phosphorylated by CDK7 that participates in chromatin condensation. Our results demonstrated the important role of CDK7 in all key aspects of spermatogenesis, potentially providing an effective target for clinical diagnosis and pathogenesis.

Quantitative phosphoproteomics reveals GSK3A substrate network is involved in the cryodamage of sperm motility

During sperm cryopreservation, the most significant phenotype of cryodamage is the decrease in sperm motility. Several proteomics studies have already been performed to search for key regulators at the protein level. However, sperm functions are known to be highly regulated by phosphorylation signaling. Here, we constructed a quantitative phosphoproteome to investigate the expression change of phosphorylated sites during sperm cryopreservation. A total of 3107 phosphorylated sites are identified and 848 of them are found to be significantly differentially expressed (DE). Bioinformatics analysis showed that the corresponding genes of these regulated sites are highly associated with sperm motility, providing a connection between the molecular basis and the phenotype of cryodamage. We then performed kinase enrichment analysis and successfully identified glycogen synthase kinase-3α (GSK3A) as the key kinase that may play an important role in the regulation of sperm motility. We further constructed a GSK3A centric network that could help us better understand the molecular mechanism of cryodamage in sperm motility. Finally, we also verified that GSK3A was abnormally activated during this process. The presented phosphoproteome and functional associations provide abundant research resources for us to learn the regulation of sperm functions, as well as to optimize the cryoprotectant for sperm cryopreservation.

ZDHHC19 Is Dispensable for Spermatogenesis, but Is Essential for Sperm Functions in Mice

Spermatogenesis is a complicated process involving mitotically proliferating spermatogonial cells, meiotically dividing spermatocytes, and spermatid going through maturation into spermatozoa. The post-translational modifications of proteins play important roles in this biological process. S-palmitoylation is one type of protein modifications catalyzed by zinc finger Asp-His-His-Cys (ZDHHC)-family palmitoyl S-acyltransferases. There are 23 mammalian ZDHHCs that have been identified in mouse. Among them, Zdhhc19 is highly expressed in adult testis. However, the in vivo function of Zdhhc19 in mouse spermatogenesis and fertility remains unknown. In this study, we knocked out the Zdhhc19 gene by generating a 2609 bp deletion from exon 3 to exon 6 in mice. No differences were found in testis morphology and testis/body weight ratios upon Zdhhc19 deletion. Spermatogenesis was not disrupted in Zdhhc19 knockout mice, in which properly developed TRA98+ germ cells, SYCP3+ spermatocytes, and TNP1+ spermatids/spermatozoa were detected in seminiferous tubules. Nevertheless, Zdhhc19 knockout mice were male infertile. Zdhhc19 deficient spermatozoa exhibited multiple defects including abnormal morphology of sperm tails and heads, decreased motility, and disturbed acrosome reaction. All of these led to the inability of Zdhhc19 mutant sperm to fertilize oocytes in IVF assays. Taken together, our results support the fact that Zdhhc19 is a testis enriched gene dispensable for spermatogenesis, but is essential for sperm functions in mice.

Comprehensive bioinformatics analysis reveals kinase activity profiling associated with heart failure

Heart failure is a complex clinical syndrome originating from cardiac injury, which leads to considerable morbidity and mortality. Among the dynamic molecular adaptations occurring in heart failure development, aggravation of the disease is often attributed to global or local abnormality of the kinase. Therefore, the overall monitoring of kinase activity is indispensable. In this study, a bioinformatics analysis method was developed to conduct deep mining of transcriptome and phosphoproteome in failing heart tissue. A total of 982 differentially expressed genes and 9781 phosphorylation sites on 3252 proteins were identified. Via upstream regulator relations and kinase-substrate relations, a dendrogram of kinases can be constructed to monitor its abnormality. The results show that, on the dendrogram, the distribution of kinases demonstrated complex kinase activity changes and certain rules that occur during heart failure. Finally, we also identified the hub kinases in heart failure and verified the expression of these kinases by reverse-transcription polymerase chain reaction and Western blot analysis. In conclusion, for the first time, we have systematically analyzed the differences in kinases during heart failure and provided an unprecedented breadth of multi-omics data. These results can bring about a sufficient data foundation and novel research perspectives.

Integrated analysis of phosphoproteome and ubiquitylome in epididymal sperm of buffalo (Bubalus bubalis)

In mammals, sperm need to mature in the epididymis to gain fertilization competency. However, the molecular mechanism underlying buffalo sperm maturation remains elusive. Exploring sperm physiology at the posttranslational modification (PTM) level could help to develop our understanding of these mechanisms. Protein phosphorylation and ubiquitination are major PTMs in the regulation of many biological processes. In the present study, to our knowledge, we report the first phosphoproteome and ubiquitylome of sperm collected from the caput, corpus, and cauda segments of the epididymis using liquid chromatography-mass spectrometry combined with affinity purification. In total, 647 phosphorylation sites in 294 proteins and 1063 ubiquitination sites in 446 proteins were characterized. Some of these proteins were associated with cellular developmental processes and energy metabolic pathways. Interestingly, 84 proteins were both phosphorylated and ubiquitinated, simultaneously. Some of these proteins were involved in, for example, spermatogenesis, reproduction, and spermatid development. Taken together, these data provide a theoretical basis for further functional analysis of phosphorylation and ubiquitination in epididymal sperm of buffalo and other mammals, and serve as an important resource for exploring the physiological mechanism underlying sperm maturation.

Novel targets identified by integrated proteomic and phosphoproteomic analysis in spermatogenesis of swamp buffalo (Bubalus bubalis)

To understand mechanisms of spermatogenesis, the proteome and the phosphoproteome in prepubertal and pubertal swamp buffalo (Bubalus bubalis) testes were analyzed using tandem mass tag (TMT) coupled with liquid chromatography-tandem mass spectrometry (LC–MS/MS). In prepubertal testes, 80 proteins were overexpressed, 148 proteins were underexpressed, and 139 and 142 protein sites had higher and lower phosphorylation, respectively, compared to the levels in pubertal testes. Several of these proteins were associated with reproductive processes such as sexual reproduction, spermatogenesis, fertilization, and spermatid development. In particular, outer dense fiber protein 1 (ODF1), protein maelstrom homolog (MAEL), actin-like protein 7B (ACTL7B), tyrosine-(Y)-phosphorylation regulated (CABYR), and tripartite motif containing 36 (TRIM36) were upregulated with age at both the proteome and phosphoproteome levels. Combining proteome and phosphoproteome analysis can be effectively applied to study the protein/phosphorylation patterns of buffalo testes. These data provide new regulatory candidates and evidence for a complex network in spermatogenesis in buffalo testes, and serve as an important resource for exploring the physiological mechanism of spermatogenesis in mammals.

In Silico Tools and Phosphoproteomic Software Exclusives

Proteomics and phosphoproteomics have been emerging as new dimensions of omics. Phosphorylation has a profound impact on the biological functions and applications of proteins. It influences everything from intrinsic activity and extrinsic executions to cellular localization. This post-translational modification has been subjected to detailed study and has been an object of analytical curiosity with the advent of faster instrumentation. The major strength of phosphoproteomic research lies in the fact that it gives an overall picture of the workforce of the cell. Phosphoproteomics gives deeper insights into understanding the mechanism behind development and progression of a disease. This review for the first time consolidates the list of existing bioinformatics tools developed for phosphoproteomics. The gap between development of bioinformatics tools and their implementation in clinical research is highlighted. The challenge facing progress is ideally believed to be the interdisciplinary arena this field of research is associated with. For meaningful solutions and deliverables, these tools need to be implemented in clinical studies for obtaining answers to pharmacodynamic questions, saving time, costs and energy. This review hopes to invoke some thought in this direction.

The Protein Phosphorylation Landscape of Mouse Spermatids during Spermiogenesis

The characteristic tadpole shape of sperm is formed from round spermatids via spermiogenesis, a process which results in dramatic morphological changes in the final stage of spermatogenesis in the testis. Protein phosphorylation, as one of the most important post-translational modifications, can regulate spermiogenesis; however, the phosphorylation events taking place during this process have not been systematically analyzed. In order to better understand the role of phosphorylation in spermiogenesis, large-scale phosphoproteome profiling is performed using IMAC and TiO₂ enrichment. In total, 13 835 phosphorylation sites, in 4196 phosphoproteins, are identified in purified mouse spermatids undergoing spermiogenesis in two biological replicates. Overall, 735 testis-specific proteins are identified to be phosphorylated, and are expressed at high levels during spermiogenesis. Gene ontology analysis shows enrichment of the identified phosphoproteins in terms of histone modification, cilium organization, centrosome and the adherens junction. Further characterization of the kinase-substrate phosphorylation network demonstrates enrichment of phosphorylation substrates related to the regulation of spermiogenesis. This global protein phosphorylation landscape of spermiogenesis shows wide phosphoregulation across a diverse range of processes during spermiogenesis and can help to further characterize the process of sperm generation. All MS data are available via ProteomeXchange with the identifier PXD011890.

Human Testis Phosphoproteome Reveals Kinases as Potential Targets in Spermatogenesis and Testicular Cancer

Spermatogenesis is a complex cell differentiation process that includes marked genetic, cellular, functional and structural changes. It requires tight regulation, because disturbances in any of the spermatogenic processes would lead to fertility deficiencies as well as disorders in offspring. To increase our knowledge of signal transduction during sperm development, we carried out a large-scale identification of the phosphorylation events that occur in the human male gonad. Metal oxide affinity chromatography using TiO2 combined with LC-MS/MS was conducted to profile the phosphoproteome of adult human testes with full spermatogenesis. A total of 8187 phosphopeptides derived from 2661 proteins were identified, resulting in the most complete report of human testicular phosphoproteins to date. Phosphorylation events were enriched in proteins functionally related to spermatogenesis, as well as to highly active processes in the male gonad, such as transcriptional and translational regulation, cytoskeleton organization, DNA packaging, cell cycle and apoptosis. Moreover, 174 phosphorylated kinases were identified. The most active human protein kinases in the testis were predicted both by the number of phosphopeptide spectra identified and the phosphorylation status of the kinase activation loop. The potential function of cyclin-dependent kinase 12 (CDK12) and p21-activated kinase 4 (PAK4) has been explored by in silico, protein-protein interaction analysis, immunodetection in testicular tissue, and a functional assay in a human embryonal carcinoma cell line. The colocalization of CDK12 with Golgi markers suggests a potential crucial role of this protein kinase during sperm formation. PAK4 has been found expressed in human spermatogonia, and a role in embryonal carcinoma cell response to apoptosis has been observed. Together, our protein discovery analysis confirms that phosphoregulation by protein kinases is highly active in sperm differentiation and opens a window to detailed characterization and validation of potential targets for the development of drugs modulating male fertility and tumor behavior.

Brain-derived neurotrophic factor: A steroidogenic regulator of Leydig cells

The brain-derived neurotrophic factor (BDNF) was first recognized for its roles in the peripheral and central nervous systems, and its complex functions on mammalian organs have been extended constantly. However, to date, little is known about its effects on the male reproductive system, including the steroidogenesis of mammals. The purpose of this study was to elucidate the effects of BDNF on testosterone generation of Leydig cells and the underlying mechanisms. We found that BDNF-induced proliferation of TM3 Leydig cells via upregulation of proliferating cell nuclear antigen ( Pcna) and promoted testosterone generation as a result of upregulation of steroidogenic acute regulatory protein ( Star), 3b-hydroxysteroid dehydrogenase ( Hsd3b1), and cytochrome P450 side-chain cleavage enzyme ( Cyp11a1) both in primary Leydig cells and TM3 Leydig cells, which were all attenuated in Bdnf knockdown TM3 Leydig cells. Furthermore, the possible mechanism of testosterone synthesis was explored in TM3 Leydig cells. The results showed that BDNF enhanced extracellular signal-regulated protein kinases 1 and 2 (ERK1/2) phosphorylation, and the effect was disrupted by Bdnf deletion. Moreover, PD98059, a potent selective inhibitor of ERK1/2 activation, compromised BDNF-induced testosterone generation and upregulation of Star, Hsd3b1, and Cyp11a1. The Bdnf knockdown assay, on the other hand, indicated the autocrine effect of BDNF on steroidogenesis in TM3 Leydig cells. On the basis of these results, we concluded that BDNF, acting as an autocrine factor, induced testosterone generation as a result of the upregulation of Star, Hsd3b1, and Cyp11a1 via stimulation of the ERK1/2 pathway.

Clinically Relevant Post-Translational Modification Analyses-Maturing Workflows and Bioinformatics Tools

Post-translational modifications (PTMs) can occur soon after translation or at any stage in the lifecycle of a given protein, and they may help regulate protein folding, stability, cellular localisation, activity, or the interactions proteins have with other proteins or biomolecular species. PTMs are crucial to our functional understanding of biology, and new quantitative mass spectrometry (MS) and bioinformatics workflows are maturing both in labelled multiplexed and label-free techniques, offering increasing coverage and new opportunities to study human health and disease. Techniques such as Data Independent Acquisition (DIA) are emerging as promising approaches due to their re-mining capability. Many bioinformatics tools have been developed to support the analysis of PTMs by mass spectrometry, from prediction and identifying PTM site assignment, open searches enabling better mining of unassigned mass spectra-many of which likely harbour PTMs-through to understanding PTM associations and interactions. The remaining challenge lies in extracting functional information from clinically relevant PTM studies. This review focuses on canvassing the options and progress of PTM analysis for large quantitative studies, from choosing the platform, through to data analysis, with an emphasis on clinically relevant samples such as plasma and other body fluids, and well-established tools and options for data interpretation.

Loss of TIGAR Induces Oxidative Stress and Meiotic Defects in Oocytes from Obese Mice

Maternal obesity has been reported to impair oocyte quality in mice, however, the underlying mechanisms remain unclear. In the present study, by conducting a comparative proteomic analysis, we identified a reduced expression of TIGAR (TP53-induced glycolysis and apoptosis regulator) protein in ovulated oocytes from high-fat diet (HFD)-fed mice. Specific depletion of TIGAR in mouse oocytes results in the marked elevation of reactive oxygen species (ROS) levels and the failure of meiotic apparatus assembly. Importantly, forced expression of TIGAR in HFD oocytes not only attenuates ROS production, but also partly prevents spindle disorganization and chromosome misalignment during meiosis. Meantime, we noted that TIGAR knockdown in oocytes induces a strong activation of autophagy, whereas overexpression of TIGAR significantly reduces the LC3 accumulation in HFD oocytes. By anti-oxidant treatment, we further demonstrated that such an autophagic response is dependent on the TIGAR-controlled ROS production. In summary, our data indicate a role for TIGAR in modulating redox homeostasis during oocyte maturation, and uncover that loss of TIGAR is a critical pathway mediating the effects of maternal obesity on oocyte quality.

Phosphoproteomics-Based Profiling of Kinase Activities in Cancer Cells

Cellular signaling, predominantly mediated by phosphorylation through protein kinases, is found to be deregulated in most cancers. Accordingly, protein kinases have been subject to intense investigations in cancer research, to understand their role in oncogenesis and to discover new therapeutic targets. Despite great advances, an understanding of kinase dysfunction in cancer is far from complete.A powerful tool to investigate phosphorylation is mass-spectrometry (MS)-based phosphoproteomics, which enables the identification of thousands of phosphorylated peptides in a single experiment. Since every phosphorylation event results from the activity of a protein kinase, high-coverage phosphoproteomics data should indirectly contain comprehensive information about the activity of protein kinases.In this chapter, we discuss the use of computational methods to predict kinase activity scores from MS-based phosphoproteomics data. We start with a short explanation of the fundamental features of the phosphoproteomics data acquisition process from the perspective of the computational analysis. Next, we briefly review the existing databases with experimentally verified kinase-substrate relationships and present a set of bioinformatic tools to discover novel kinase targets. We then introduce different methods to infer kinase activities from phosphoproteomics data and these kinase-substrate relationships. We illustrate their application with a detailed protocol of one of the methods, KSEA (Kinase Substrate Enrichment Analysis). This method is implemented in Python within the framework of the open-source Kinase Activity Toolbox (kinact), which is freely available at http://github.com/saezlab/kinact/ .

Changes in testicular function proteins and sperm acrosome status in rats treated with valproic acid

Valproic acid (VPA), an anti-epileptic drug, reduces testosterone levels and sperm quality. However, the degree to which testosterone levels and sperm quality are decreased under VPA treatment needs to be clarified. The aim of the present study was to investigate the testicular proteins involved in testosterone synthesis and spermatogenesis, histopathology and sperm acrosome status in VPA-treated rats. Adult rats were divided into control and experimental groups (n = 8 in each). Rats in the experimental group were treated with 500 mg kg–1, i.p., VPA for 10 consecutive days. Expression of Ki-67, tyrosine phosphorylated proteins and testicular steroidogenic proteins was examined. As expected, VPA-treated rats exhibited adverse changes in almost all reproductive parameters, particularly an increase in precocious acrosome reactions, compared with the control group. In addition, fibrosis of the tunica albuginea and tubule basement membrane was observed in testes from VPA-treated rats. Moreover, the expression of testicular Ki-67, cholesterol side-chain cleavage enzyme (P450scc) and phosphorylated proteins (41, 51 and 83 kDa) was decreased significantly in VPA-treated rats compared with control. In contrast, the expression of steroidogenic acute regulatory proteins in the VPA-treated group was significantly higher than in the control group. In conclusion, VPA treatment changes the expression of testicular proteins responsible for spermatogenesis and testosterone production, resulting in male infertility.

dbPAF: an integrative database of protein phosphorylation in animals and fungi

Protein phosphorylation is one of the most important post-translational modifications (PTMs) and regulates a broad spectrum of biological processes. Recent progresses in phosphoproteomic identifications have generated a flood of phosphorylation sites, while the integration of these sites is an urgent need. In this work, we developed a curated database of dbPAF, containing known phosphorylation sites in H. sapiens, M. musculus, R. norvegicus, D. melanogaster, C. elegans, S. pombe and S. cerevisiae. From the scientific literature and public databases, we totally collected and integrated 54,148 phosphoproteins with 483,001 phosphorylation sites. Multiple options were provided for accessing the data, while original references and other annotations were also present for each phosphoprotein. Based on the new data set, we computationally detected significantly over-represented sequence motifs around phosphorylation sites, predicted potential kinases that are responsible for the modification of collected phospho-sites, and evolutionarily analyzed phosphorylation conservation states across different species. Besides to be largely consistent with previous reports, our results also proposed new features of phospho-regulation. Taken together, our database can be useful for further analyses of protein phosphorylation in human and other model organisms. The dbPAF database was implemented in PHP + MySQL and freely available at http://dbpaf.biocuckoo.org.

IKAP: A heuristic framework for inference of kinase activities from Phosphoproteomics data

MOTIVATION

Phosphoproteomics measurements are widely applied in cellular biology to detect changes in signalling dynamics. However, due to the inherent complexity of phosphorylation patterns and the lack of knowledge on how phosphorylations are related to functions, it is often not possible to directly deduce protein activities from those measurements. Here, we present a heuristic machine learning algorithm that infers the activities of kinases from Phosphoproteomics data using kinase-target information from the PhosphoSitePlus database. By comparing the estimated kinase activity profiles to the measured phosphosite profiles, it is furthermore possible to derive the kinases that are most likely to phosphorylate the respective phosphosite.

RESULTS

We apply our approach to published datasets of the human cell cycle generated from HeLaS3 cells, and insulin signalling dynamics in mouse hepatocytes. In the first case, we estimate the activities of 118 at six cell cycle stages and derive 94 new kinase-phosphosite links that can be validated through either database or motif information. In the second case, the activities of 143 kinases at eight time points are estimated and 49 new kinase-target links are derived.

AVAILABILITY AND IMPLEMENTATION

The algorithm is implemented in Matlab and be downloaded from github. It makes use of the Optimization and Statistics toolboxes. https://github.com/marcel-mischnik/IKAP.git.

CONTACT

marcel.mischnik@gmail.com

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.