Next Article in Journal
Western Diet Induced Remodelling of the Tongue Proteome
Next Article in Special Issue
H2B Type 1-K Accumulates in Senescent Fibroblasts with Persistent DNA Damage along with Methylated and Phosphorylated Forms of HMGA1
Previous Article in Journal
Macrophage Proteome Analysis at Different Stages of Mycobacterium avium Subspecies paratuberculosis Infection Reveals a Mechanism of Pathogen Dissemination
Previous Article in Special Issue
Small Mass but Strong Information: Diagnostic Ions Provide Crucial Clues to Correctly Identify Histone Lysine Modifications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Protamine Characterization by Top-Down Proteomics: Boosting Proteoform Identification with DBSCAN

by
Gianluca Arauz-Garofalo
1,
Meritxell Jodar
2,3,
Mar Vilanova
1,
Alberto de la Iglesia Rodriguez
2,
Judit Castillo
2,
Ada Soler-Ventura
2,
Rafael Oliva
2,3,*,
Marta Vilaseca
1,* and
Marina Gay
1,*
1
Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac, 10, 08028 Barcelona, Spain
2
Molecular Biology of Reproduction and Development Research Group, Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Fundació Clínic per a la Recerca Biomèdica (FCRB), Faculty of Medicine and Health Sciences, University of Barcelona, 08036 Barcelona, Spain
3
Biochemistry and Molecular Genetics Service, Hospital Clínic de Barcelona, 08036 Barcelona, Spain
*
Authors to whom correspondence should be addressed.
Submission received: 5 March 2021 / Revised: 25 April 2021 / Accepted: 27 April 2021 / Published: 30 April 2021
(This article belongs to the Special Issue Latest Discoveries and Proteomics-Based Methods in Histone Biology)

Abstract

:
Protamines replace histones as the main nuclear protein in the sperm cells of many species and play a crucial role in compacting the paternal genome. Human spermatozoa contain protamine 1 (P1) and the family of protamine 2 (P2) proteins. Alterations in protamine PTMs or the P1/P2 ratio may be associated with male infertility. Top-down proteomics enables large-scale analysis of intact proteoforms derived from alternative splicing, missense or nonsense genetic variants or PTMs. In contrast to current gold standard techniques, top-down proteomics permits a more in-depth analysis of protamine PTMs and proteoforms, thereby opening up new perspectives to unravel their impact on male fertility. We report on the analysis of two normozoospermic semen samples by top-down proteomics. We discuss the difficulties encountered with the data analysis and propose solutions as this step is one of the current bottlenecks in top-down proteomics with the bioinformatics tools currently available. Our strategy for the data analysis combines two software packages, ProSight PD (PS) and TopPIC suite (TP), with a clustering algorithm to decipher protamine proteoforms. We identified up to 32 protamine proteoforms at different levels of characterization. This in-depth analysis of the protamine proteoform landscape of normozoospermic individuals represents the first step towards the future study of sperm pathological conditions opening up the potential personalized diagnosis of male infertility.

1. Introduction

Protamines are the main nuclear proteins in the sperm cells of many species, and they play a crucial role in the correct packaging of the paternal DNA [1,2,3]. During the final stage of spermatogenesis (called spermiogenesis), the chromatin of immature germ cells undergoes marked remodeling, in which histones are sequentially replaced by specific histone variants, transition proteins and, finally, protamines [1,4,5,6,7]. Both the small size of protamines (20 to 55 amino acids) and their extremely rich sequence in positively charged arginine residues allow the formation of highly condensed toroidal complexes, together with the negatively charged paternal DNA [8,9]. In addition, the high number of cysteine residues in mammalian protamines promotes the generation of multiple intra- and inter-molecular disulfide bonds and zinc bridges between protamines, thereby stabilizing the compact toroidal nucleoprotamine complex [10,11]. Several functions have been proposed for the nucleoprotamine domain, including the protection of the paternal genome during transit from the male to the oocyte and contribution to obtaining the required hydrodynamic shape for mature sperm functionality [3,4,12,13,14,15,16,17]. Of note, the replacement of histones by protamines during spermatogenesis is not complete, resulting in 85–95% of sperm DNA packed by protamines and the remaining 5–15% attached to histones [18]. Several studies have reported that the distribution of the nucleohistone and nucleoprotamine complexes along sperm chromatin is not random. This distinctive chromatin structure has been proposed as a specific epigenetic mark of sperm that could regulate gene expression once the oocyte has been fertilized [19,20,21].
The crucial role of protamines in the compaction of the paternal genome has given rise to considerable interest in the analysis of protamine alterations in the andrology field. Mammalian spermatozoa contain two types of protamines: protamine 1 (P1) and the family of protamine 2 (P2) proteins. P1 is expressed in all species of protamine-containing vertebrates, while the P2 family is present only in some mammalian species, including the human and mouse. P2 is formed by the proteolysis of the P2 precursor (pre-P2), resulting in the HP2, HP3 and HP4 components, of which HP2 is the most abundant [4] (Figure 1).
Several groups have established similar protein levels of P1 and P2 in normal spermatozoa, with an approximate 1:1 proportion, and alterations of the P1/P2 ratio have been associated with male infertility and poor reproductive outcomes [22,23,24]. However, controversial results have been reported in the literature and the P1/P2 ratio shows a wide range in fertile men [25]. The P1/P2 ratio has been commonly assessed by polyacrylamide gel electrophoresis (PAGE) and protein staining [26,27,28,29]. However, this common approach does not provide information on whether a specific alteration of the P1/P2 ratio is due to the deregulation of P1, the P2 family, specific P2 components, or all of three. In addition, the emerging hypothesis that alterations in both histone variants and histone post-translational modifications (PTMs) might lead to male infertility [30] suggests that, in a similar way, alterations in protamine PTMs are also associated with male infertility. Recently, we suggested that the presence of human sperm protamine proteoforms, including truncations or proteoforms containing PTMs, explain the variability observed in the P1/P2 ratio assessed by conventional PAGE [30,31]. Therefore, we propose that more sophisticated methodologies such as mass spectrometry (MS) could help decipher the different protamine proteoforms present in mature sperm cells and their impact on male fertility [31]. However, the peculiar physicochemical properties of protamines (high arginine content with P1 containing 24 Arg residues out of 51 and HP2 32 Arg residues out of 54) and the high similarity between P2 components make them unsuitable for characterization by bottom-up proteomic strategies. The main reason is that conventional trypsin digestion results in extremely small protamine peptides that are lost during their processing. Additionally, the intra- and inter-protamine disulfide bonds and zinc bridges that form due to the high number of cysteine residues in the protamine sequences (6 Cys in P1 and 5 Cys in HP2) require methods for reduction and alkylation to avoid multiple cysteine oxidations. Of note, the low molecular weight of protamines (P1: 6687 Da and P2: 12,912 Da) is a big advantage for analyzing these molecules by top-down proteomics.
Top-down proteomics is a powerful MS-based technology that allows direct analysis of all proteoforms, whether derived from alternative splicing, missense or nonsense genetic variants or PTMs. This approach offers unique and precise information on proteoforms that is complementary to that achieved using bottom-up proteomics. The top-down analysis of MS-derived data is challenging both due to the deconvolution of the spectra and subsequent searches in databases for correct proteoform assignment. Unlike bottom-up proteomics, where robust solutions are found for the identification, characterization and quantification of peptides and proteins, data analysis in top-down proteomics is not yet fully established. In this regard, it is still one of the bottlenecks for top-down proteomic strategies. One of the main pitfalls is the lack of automated validation of the final IDs, which requires manual inspection of the assignments.
There are two main approaches for top-down data analysis and proteoform identification: (a) analysis based on a search in annotated databases, where all the possible proteoforms with defined PTMs are included in a certain position, like in ProSight PD (PS) [32], and (b) analysis based on spectral alignment algorithms [33], which matches the tandem mass spectra of the proteins with the unannotated database, without prior knowledge of the type of PTMs or their location, like in TopPIC suite (TP) [34,35].
The data analysis workflow that we describe herein combines two distinct software packages (referred to herein as nodes), namely PS and TP. To avoid redundant identifications on the same precursor, our pipeline is based on MS1 clustering and harmonization of the two datasets obtained. This new strategy allowed us to identify 13 proteoforms of P1 and 19 proteoforms of P2 in human semen samples. We envisioned that this strategy could be applied to other highly modified proteins such as histones.

2. Materials and Methods

2.1. Sample Collection

Human semen samples (n = 2) were obtained from patients undergoing routine semen analysis at the Assisted Reproduction Unit at the Clinic Institute of Gynecology, Obstetrics and Neonatology, at the Hospital Clínic of Barcelona, Spain. All patients gave signed informed consent in accordance with the Declaration of Helsinki. The ejaculates were collected by masturbation into sterile containers after 3–5 days of sexual abstinence. Seminal parameters were evaluated using the automatic semen analysis system CASA (Proiser, Paterna, Spain). After seminal analysis, all samples were classified as normozoospermic according to World Health Organization guidelines [36]. Subsequently, sperm was purified through 50% Puresperm® density gradient separation (NidaCon International AB, Gothenburg, Sweden), following the manufacturer’s instructions.

2.2. Purification of Human Sperm Protamine

The protamine-enriched fraction of spermatozoa was purified as described in [27]. except that histones and other basic proteins weakly associated with DNA were removed from the sperm by 0.5 M HCl, vortexing and incubation for 10 min at 37 °C and centrifugation at 2000× g for 20 min at 4 °C (three times) prior to suspending the sperm pellet with 0.5% Triton X-100, 20 mM (pH8) Trizma® base HCl, and 2 mM MgCl2. Samples were then processed as described in [27]. For the top-down MS approach, protamine-enriched fractions (equivalent to 15 million spermatozoa) were dried at room temperature overnight under a fume hood and kept at 20 °C. Before liquid chromatography MS analysis (LC-MS/MS), samples were reconstituted in 50 mM NH4HCO3, reduced with 10 mM dithiothreitol (DTT) for 45 min at 56 °C, alkylated for 30 min in the dark with 50 mM iodoacetamide (IAA) and further desalted using PolyLC C18 filter tips (PolyLC INC., Columbia, MD, USA).

2.3. Protamine Quantification

Purified extracts of protamines were quantified as described in [27]. Briefly, an aliquot of protamines for each sample was separated in an acid-urea PAGE together with increasing amounts of a protamine standard obtained as described elsewhere [37]. Approximate amounts of P1 and P2 were calculated from protamine optical density measurements and the quantification obtained with the protamine standard, using Quantity One 1-D Analysis Software (Bio-rad, Hercules, CA, USA).

2.4. Top-Down Proteomics

2.4.1. LC-MS

We analyzed two biological replicates, each with two technical replicates. Samples were reconstituted in an aqueous solution of 3% acetonitrile (ACN) and 1% formic acid (FA) and injected to the LC-MS/MS system.
LC-MS/MS analysis of intact protamines was conducted on a Dionex Ultimate 3000 coupled to an Orbitrap Fusion Lumos mass spectrometer (both from Thermo Scientific, San Jose, CA, USA). Samples were loaded onto C18 trap columns (100 µm × 2 cm Acclaim PepMap100, 5 µm, 100 Å (Thermo Scientific) at a flow rate of 15 µL/min. Protein content was separated using C18 analytical columns (Acclaim PepMap® RSLC 75 µm × 50 cm, nanoViper, C18, 2 µm, 100Å (Thermo Scientific) and a linear gradient from 3 to 15% B in 30 min at a flow rate of 250 nl/min (A = 0.1% FA in water, B = 0.1% FA in ACN). The Advion TriVersa NanoMate (Advion Ltd., Harlow, United Kingdom) was used as the nanospray interface with the direct on-line coupling mode with a spray voltage of 1.60 kV. The mass spectrometer was operated in data-dependent acquisition (DDA) mode. Survey MS scans were acquired in the orbitrap with the resolution (defined at 200 m/z) set at 120k. The top N (most intense) ions per scan were fragmented by electron transfer dissociation (ETD) and detected in the orbitrap at 120k. The ion count target value was 400,000 for the survey scan and 1,000,000 for the ETD-fragmented MS/MS scan. Target ions already selected for MS/MS were dynamically excluded for 30 s. RF lens was tuned to 30%. The minimal signal required to trigger the MS to MS/MS switch was set at 500,000.

2.4.2. Proteoform Search

We devised a twin search approach for our top-down proteomics workflow (Figure 2). This strategy encompasses the bioinformatic stools TP and PS.
The first node used for top-down searches was TP v1.4.2, including its tool TopFD to deconvolute the spectra. This node can identify proteoforms of an entire proteome with remarkable speed by searching against a FASTA database without prior knowledge of possible modifications or existing PTMs. Therefore, TP is particularly suitable when dealing with highly modified proteins such as histones or protamines.
We included the following PTMs in the parameters of the TP search: oxidation (in M), phosphorylation (in STY), acetylation (in KS), methylation, dimethylation (both in KR) and trimethylation (in K) for being the major PTMs found in histones and therefore also expected in protamines. We included the following N-terminal forms: none, N-terminal methionine excision (NME), NME_Acetylation and M_Acetylation, which referred to proteoforms with or without the initial methionine non-acetylated or acetylated. We also defined precursor and fragment mass tolerances of 15 ppm, and maximum mass shifts of 500 Da. We only considered Proteoform Spectrum Matches (PrSMs) with FDR < 1% for downstream data integration. We used a FASTA file containing P1 and P2 sequences as a database.
The second search was done with Proteome Discoverer v2.5.0.400 (Thermo Fisher Scientific) and its PS v4.0.0.228 node. In contrast with TP, PS requires a database with all theoretical proteoforms explicitly listed. The Protein Annotation option from the PS Database Manager allows such a database to be built from a typical *.xml file available in SwissProt.
It is crucial to mention here the existing bioinformatic bottleneck despite current computing power. As the number of PTMs for a given protein sequence increases, the number of entries in the database rapidly explodes. We can compute the number of theoretical proteoforms for a given database, N, using the following expression:
N = j = 1 P i = 1 n m i + 1 j
where P is the number of protein sequences in the database, n the number of amino acids for each sequence, and m the number of possible modifications for each amino acid on each sequence [38]. Given the PTMs typically found in protamines that we previously mentioned, this equation returns around NP1 ≈ 4.4 × 1015 and NP2 ≈ 3.9 × 1022 theoretical proteoforms for P1 and P2, respectively, which gives NNP2 ≈ 1022 in total. These values would imply an extremely long computational search time. To overcome this issue, we launched a preliminary PS search excluding those PTMs on consecutively repeated arginine residues, but keeping PTMs in selected arginines from the N-terminal, C-terminal and from the middle of the sequence. As a result, the number of theoretical proteoforms in this first search was reduced by 15 orders of magnitude. Since this preliminary search did not find any acetylations, methylations, dimethylations or trimethylations, we launched a second three-tier search (Annotated Proteoform, Subsequent Search and Annotated Proteoform) using a database with only phosphorylations (in STY), oxidations (in M) and acetylation in N-terminal in all the nodes. Using this strategy, we further reduced the number of theoretical proteoforms by 2 orders of magnitude, adding up a total of 17 orders of magnitude and thereby giving a much more reasonable search time. Precursor mass tolerances were 2.2 Da/10 ppm/500 Da for Annotated Proteoform, Subsequent Search and Annotated Proteoformrespectively, and fragment mass tolerances 10 ppm in all cases. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [39]. partner repository with the dataset identifier PXD024405.

2.4.3. Data Analysis

The PrSM output format used by PS and TP is highly heterogeneous. The column name used by each node to designate the same concept usually varies (for example: “Spectrum File” in PS and “Data file name” in TP). Even the unit used to express the same magnitude differs (for example: TP expresses retention time in minutes while PS uses seconds). We computationally addressed this issue with a custom Python code [40]. (available on the GitHub repository: https://github.com/MSPCF/2021_Proteomes, accessed on 27 April 2021). To this end, we leveraged some of Python’s most popular packages such as Pandas [41], NumPy [42], Biopython [43], Plotnine [44], and Seaborn [45].
Once this harmonized dataset with all PrSMs was ready, we filtered out those identifications with retention times beyond 42 min for being associated with the column wash. We did not apply any PrSMs pre-filter based on PrSM quality (such as the MS1 intensity or the E-value), but this might be convenient depending on the overall PrSM dataset quality. After pre-filtering, we performed the PrSM clustering by means of a density-based spatial clustering algorithm of applications with noise (DBSCAN) [46,47]. implemented within the Scikit-learn Python package [48]. This clustering algorithm is particularly suitable in our scenario because it does not require pre-definition of the number of clusters expected (a priori, we do not know the number of proteoforms within a sample). In addition, DBSCAN has just two main hyperparameters to work with: maximum distance between two samples to accept both as belonging to the same neighborhood (ε), and minimum number of samples in a neighborhood to be accepted as an independent cluster (nmin). In our scenario, “samples” and “independent clusters” would correspond to PrSMs and proteoforms, respectively. Regarding the “distance”, it is important to mention that we simply passed the standard score (commonly called z-score) of the deconvoluted PrSM masses to the clustering algorithm. Finally, we optimized the DBSCAN hyperparameters ε and nmin by interactive exploratory data analysis with a custom Python code leveraging the Bokeh Python package [49], getting ε = 1.2 × 10−3 (which in our filtered harmonized dataset corresponds to 1.28 Da) and nmin = 10.

2.4.4. Proteoform Annotation

Once clustering was done, proteoforms were annotated in a semi-automatic manner comprising two major steps. For a given cluster, we first checked to see if most of its PrSMs pointed to a particular protein, usually any full-length or truncated mature or immature protamine form or other sperm-related protein. In this case, we pre-assigned what we call a “cluster-protein” to this cluster and then computed its corresponding cluster-protein mass shift as follows:
Δ m C T = m C m T ,
where mT is the theoretical mass of the fully carbamidomethylated cluster-protein, and mC the mean mass of all PrSMs within this cluster. When a predominant cluster-protein was not found, we manually attempted to assign a suitable protein consistent with the particular mC value. Secondly, and depending on the precise ΔmCT we manually assigned a suitable PTM set. The proteoform annotation table is provided in the supplementary material (Table S1), as well as the associated harmonized PrSMs dataset, including the clustering outcome and the proteform annotation (Table S2).

3. Results and Discussion

We analyzed two biological replicates (A and B), with two technical replicates each (R01 and R02 for A, R03 and R04 for B). The results presented here are all the data combined with the final aim of characterizing the protamine proteoforms in the samples.
In the following sections we will use the five-level classification system for proteoform identifications proposed by Smith et al. [50]. We will refer to levels 1, 2A, 2B and 3. Level 1 describes proteoforms identified with full knowledge of their gene origin, amino acid sequence and, when present, the identification and localization of their PTMs. Level 2A refers to those proteoforms with a defined amino acid sequence from a known gene with all the PTMs identified but with incomplete localization. In level 2B, the PTMs are not fully characterized but their localization is complete. Level 3 refers to those proteoforms that are certain for two of the four ambiguity classes: gene, sequence, PTM identification and PTM localization. In this manuscript we refer to level 3 as proteoforms with a gene and sequence defined but with an unknown modification and localization.

3.1. Overview of the DBSCAN Approach

The goal of our DBSCAN approach was to integrate and harmonize the PS and TP datasets by PrSM deconvoluted experimental masses. We assumed that a proteoform at level 1 without PTMs, or level 2A or level 3 with an equivalent mass shift, always had the same mass and therefore could be annotated as the same proteoform. In this first step, we were not able to localize modification sites and reach level 1 for proteoforms with PTMs, but we localized PTMs in a second step considering only those PrSMs identified with a high confidence level (Section 3.4). The power of our workflow for top-down proteomic data analysis lies in combining multiple nodes to achieve high-quality PrSMs and a single harmonized dataset.
From the PS and TP searches, a total of 2838 PrSMs passed our filters (56.6% for PS and 43.4% for TP) and were further submitted for clustering. The DBSCAN algorithm successfully clustered 2578 of the PrSMs (55.9% for PS and 44.1% for TP) into 27 well differentiated mass clusters, leaving 260 unclustered PrSMs (63.8% for PS and 36.2% for TP).
We identified the following four main protamine groups: P1 proteoforms, comprising clusters from P|000 to P|008 (7029 Da to 7311 Da); mature P2 forms (HP2, 3 and 4) from P|009 to P|013 (7538 Da to 8123 Da); intermediate P2 forms (HPS1, 2 and HPI2) from P|014 to P|020 (9002 Da to 10,654 Da); and pre-P2 forms from P|021 to P|025 (13,195 Da to 13,516 Da) (Table 1). Cluster P|026 remained unidentified and could presumably be ascribed to a sperm-related protein other than protamines. In addition to several phosphorylated proteoforms (up to +80 Da × 4), we also identified some presenting an unrecognized +61 Da mass shift modification (up to +61 Da × 3), usually in combination with other P1 phosphorylated proteoforms. We found this +61 Da mass shift in P1, HP2, HP3, HP4, HPS1 and HPS2. We had already identified this +61 Da mass shift in previous studies [31]. We cannot explain this +61 Da mass shift yet but propose that these protamines exhibit some metal ion-binding sites for Zn2+ or Cu+ [51,52]. Further verification is required using other MS-based techniques such as inductively coupled plasma MS. Finally, we also observed the formation of the cyclic pyroglutamyl residue (pyro-Glu) from the N-terminal glutamine of HPS1 [53].
We used different projections to visualize the clustering outcome and better optimize its ε and nmin parameters with our Bokeh-based web application. In particular, we inspected the following three projections: retention time vs. deconvoluted experimental mass (Figure 3a, Figure S1), log2 intensity vs. deconvoluted experimental mass (Figure S2a), and retention time vs. log2 intensity (Figure S2b). This strategy was crucial for gaining greater insight into how the DBSCAN algorithm worked, as well as the differences between the PS and TP nodes.
The scatterplot of the retention time vs. the deconvoluted experimental masses was particularly useful. We used colors to differentiate clusters. In this plot, it is easy to see proteoforms with the same mass clustered or to localize PrSMs that did not cluster. The four main protamine groups mentioned above per mass region can be visualized in the same plot, grouped in mass regions: (a) P1: 7000 to 7500 Da; (b) mature P2: 7500 to 8200 Da; (c) Intermediate P2: from 9000 to 11,000 Da; and (c) pre-P2: from 13,000 to 13,600 Da. Each region had the canonical sequences and the subsequent proteoforms with different modifications.
Using the DBSCAN approach, we identified proteoforms at levels 1, 2A and 3 of the classification system [51]. We reached level 1 in proteoforms without PTMs (truncated and non-truncated forms) and in the case of HPS1 +PyroGlu, level 2A in proteoforms with known modifications (phosphorylations), and level 3 with proteoforms with the +61 mass shift.

3.2. Advantages of the DBSCAN Approach vs. PS and TP

When handling PS and TP outputs, we encountered some common difficulties and issues, such as the different output tables generated or incorrectly assigned proteoforms. With the DBSCAN approach, we overcame some of these issues and integrated, simplified and corrected the PS and TP outputs.
We first examined the results at the amino acid sequence level without considering known or unknown modifications. TP identified 65 sequences, of which 41 were truncated at the N-terminal, 5 at the C-terminal and 4 in both the N and C-terminal regions. PS found 14 sequences, of which 9 were truncated at the N-terminus (Table S2). These numbers, based only on the amino acid sequence level, already highlighted the high dispersion between the nodes. When we introduced modifications, the dispersion between results was even higher, thereby hindering comparison of the two datasets. However, after the cluster analysis using the PS and TP data, most of these sequences fell into the 8 sequences previously described in Uniprot (Figure 1).
Going deeper into the data analysis, we observed that PS and TP sometimes failed to assign a proteoform from a given PrSM. We identified two main common cases of incorrect PrSM assignment: (a) PrSMs with very different deconvoluted masses recognized as the same proteoform (which should not be possible for any classification level), and (b) PrSMs with a very similar deconvoluted mass recognized as different proteoforms. The first misassignment case (from now on “converging”) was usually found in PS while the second one (from now on “diverging”) was found in both PS and TP. Figure 3 illustrates these two misassignment scenarios by zooming in to highlight the mass region where HPS1 and HPS2 are located (Figure 3b). For example, PS considers those PrSMs around 9000 Da and 9060 Da to be the same proteoform (Figure 3b, green data points). Another example, TP does not report those PrSMs clustered around 9283 Da as the same proteoform (Figure 3b, purple and magenta data points). Note how our DBSCAN approach nicely groups those PrSMs with similar masses into well-defined clusters (proteoforms), in striking contrast to how PS and TP perform (Figure 3c).

3.2.1. Converging and Diverging Misassignment of PrSMs

We considered converging misassignment to be those cases in which PrSMs with very different experimental deconvoluted masses (i.e., belonging to different proteoforms according to our clustering approach) are identified by the node as the same proteoform. In general, for a given TP proteoform ID we found just one proteoform according to our clustering approach, meaning that TP did not mix proteoforms with different deconvoluted masses under the same proteoform ID. On the contrary, we observed that PS “converges” into the same proteoform ID PrSMs belonging to different proteoforms according to our clustering approach. In fact, PS classified under the same proteoform ID PrSMs with deconvoluted experimental mass differences up to 980 Da. This PS issue usually occurred when there were no identified modifications, as in the +61 mass shift case.
A clear example of this situation can be observed in Figure 3b. Note how PS converged as HPS2 (theoretical mass 9002.75 Da) PrSMs with quite different deconvoluted masses (9002 Da and 9064 Da). We consider that PrSMs around 9003 Da (cluster P|015) are correctly associated with HPS2, but those around 9064 Da (cluster P|016) should be linked to a new HPS2 +61 Da proteoform because they are clustered apart (Table 2a). We manually validated this hypothesis with the help of ProSight Lite [54], by precisely localizing such +61 Da mass shifts in the same spectrum identified as HPS2 by PS (Table 2b). Note how both the P-Score and the number of matching fragments improved once we assigned this + 61 Da mass shift to the C-46 residue in HPS2 (the first decreases from 1.1 × 10−91 to 7.6 × 10−155 and the second increases from 62 to 95) (Table 2a,b). Similarly, PS also “converged” as HPS1 (9300.91 Da), PrSMs with clearly different deconvoluted masses (this time with values around 9301 Da and 9362 Da). Again, since those PrSMs were clustered apart (P|019 and P|020, respectively), we assumed that PrSMs around 9362 Da correspond to a new HPS1 +61 Da proteoform (Figure 3b).
We considered diverging misassignment to be those cases in which PrSMs with very similar experimental deconvoluted masses (i.e., belonging to the same proteoform according to our clustering approach) were identified as different proteoforms. Both PS and TP tended to make this kind of mistake. Diverging misassignment translated into the identification of multiple proteoforms (according PS or TP) within one of our mass clusters. In particular, the average number of proteoforms per cluster according to PS and TP was 2.81 and 2.96, with extreme cases in which we found up to 9 and 5 proteoforms per cluster, respectively.
A clear example of this situation can be seen in Figure 3b (for PS and TP nodes) and Table 1. Note how the PS and TP nodes “diverged” those PrSMs with deconvoluted experimental mass around 9283 Da (cluster P|018) in up to nine and three distinct proteoforms, respectively. We consider that all these PrSMs should be assigned to an HPS1 proteoform with a pyroglutamic acid at the glutamine located in its N-terminus (−18.010565 Da). Similarly, as in the previous Section 3.2.1, we manually validated our hypothesis using ProSight Lite by including and localizing this pyroglutamic acid (Table 2d–f). Note how for two particular PrSMs, which PS originally assigned to HPS1 + 2 phosphorylations (Table 2c) and HPS2 + 1 phosphorylation + 1 N-term acetylation (Table 2e), both the P-Score and the number of matching fragments improved once we localized the pyroglutamic acid at HPS1 N-terminus (Table 2d,f) (the first decreases from 7.9 × 10-38 to 2.8 × 10−103 and from 9.7 × 10−45 to 1.8 × 10−118, respectively, and the second increases from 35 to 75 and from 40 to 85, respectively).

3.2.2. DBSCAN Approach Limitations

The DBSCAN approach was able to merge, harmonize and correct the PS and TP results but had some limitations. The three main limitations were (a) the requirement for previous knowledge of the sample to perform the final semi-manual annotation; (b) the assumption that the same mass meant the same proteoform when two amino acid sequences (with or without PTMs) may in fact have the same mass; and (c) some of the identified PrSMs remained unclustered and were discarded.
At the present stage of development, proteoform assignment to each cluster is semi-automatic and requires a minimum degree of manual validation and previous knowledge of the sample. We therefore believe that this approach can easily be applied to simple samples even when they have a large number of modifications, such as protamines or histones. However, we still have to check its viability in complex samples.
One of the strengths of the DBSCAN approach was also a limitation. We assumed that one mass corresponded to one proteoform without the modification localized. This assumption corrected and simplified the data, but could also lead to incorrect assignments. It is possible for two different proteoforms to have the same mass even if they have different amino acid sequences. To address this limitation and try to distinguish this kind of situation, we intend to add the retention time as a new parameter in the DBSCAN algorithm in a future version of the workflow. This modification would better differentiate and monitor proteforms of the same mass.
As stated in Section 3.1, the DBSCAN algorithm successfully clustered 2838 of the PrSMs, leaving 260 unclustered. Moreover, only one (P|026)—comprising an additional 28 PrSMs—out of 27 clusters was unknown. As a result, our approach left a total of 288 (10.1%) unclassified PrSMs (unclustered/unidentified). However, we observed that unclassified PrSMs showed overall poorer quality than classified ones (clustered/identified).
In particular, unclassified PrSMs had a lower mean E-value and intensity than classified PrSMs for both PS and TP (Figure S3). The mean of the -log10 of the E-value was 73.4 (PS) and 26.6 (TP) for the classified PrSMs and 32.3 (PS) and 11.1 (TP) for the unclassified ones. The mean of z-score normalized log2 intensities was 0.10 (PS) and 0.18 (TP) for the clustered PrSMs and −0.77 (PS) and −1.91 (TP) for the unclustered ones. The loss of 10.1% of all PrSMs is therefore not particularly concerning since it implies that low confidence matches are discarded.
The DBSCAN approach combines the strengths of the PS and TP nodes, in particular the high spectra identification rate of PS and the ability to identify proteoforms with an arbitrary mass shift of TP. The clustering approach also overcomes the weaknesses of the two nodes, like the two kinds of proteoform misassignment highlighted in Section 3.2.1.

3.3. Protamine Proteoform Landscape and Protamine Ratios

Our DBSCAN-based approach allowed us to identify a total of 26 protamine proteoforms. The protamine forms involved in this proteoform landscape were P1, pre-P2, HPI2, HPS1, HPS2, HP2, HP3 and HP4. Among the PTMs contributing to the landscape, we found phosphorylation, N-terminal cyclic pyroglutamyl, N-terminal glutamine loss and the +61 Da modification, which was discussed in Section 3.1.
To gain greater insight into this rich protamine proteoform landscape, we performed scatterplots for the cluster-protein mass shift associated with each protamine proteoform cluster (ΔmCT in Equation (2)) vs. the corresponding protamine forms so that each data point represents a different 1, 2A- or 3-level proteoform (Figure 4). We generated this scatterplot at three levels: global, by technical replicate, and by biological replicate (Figure 4, panels a, b and c, respectively).
Note how 22 out of 25 proteoforms appearing in biological replicate A are found within its two technical replicates R01 and R02 (exceptions are P1 +1ph +2(61), pre-P2 +4ph and HP3 +1(61)). Similarly, 24 out of 26 proteoforms appearing in biological replicate B are found within its two technical replicates R03 and R04 (exceptions are P1 +2ph +2(61) and HP3 +1(61)). The only proteoform not found in the two biological replicates was P1 +2ph +2(61). By revisiting Table 1, we checked that all these exceptions corresponded to clusters with limited PrSM counts. In particular: P1 +1ph +2(61) with 20 + 29 (TP + PS) counts; pre-P2 +4ph with 6 + 8; HP3 +1(61) with 13 + 17; and P1 +2ph +2(61) with 5 + 7. After inspecting the protamine proteoform landscapes at all levels, we saw that the performance of our characterization was acceptable for intertechnical replicate variability.
As an advantage over classical protamine ratio assessment through gel-staining approaches, our strategy enables computation of the abundance ratio between any protamine proteoforms. The ratio between P1 and P2 mature proteoforms (HP2, HP3, and HP4) is the ratio most widely used for the study of male infertility [3]. Alterations of the P2 family seem to be more frequent than those of P1, and the incomplete processing of the pre-P2 has also been linked to male infertility [55,56,57]. Therefore, the ratio of the different protamine proteoforms between immature proteoforms of P2, including pre-P2 and intermediate proteoforms (HPI2, HPS1, HPS2) and mature forms of P2, is highly relevant in the context of understanding the chromatin maturity state. We provided barplots for these two ratios (P1/P2-mature and P2-immature/P2-mature) at three levels: global, by technical replicate, and by biological replicate (Figure 4, panels d, e and f, respectively). To better monitor the contribution of each node in each ratio computation, we split them using a color scale. The P1/P2-mature and P2-immature/P2-mature ratios fall around 1 for our two biological replicates and their corresponding technical replicates. Note how PS gives slightly greater ratios than TP. However, this systematic discrepancy is smaller than associated standard deviations.

3.4. Modification Site Assignment

The DBSCAN cluster approach allowed us to identify modified proteoforms at level 2A or 3. We reached level 2A with phosphorylated proteoforms, but we did not localize the p-site. We reached level 2B classification for HPS2 +61 Da shown in Table 2c. However, this localization was done manually. Undertaking this task for all PrSMs would have been time-consuming, so we accepted that proteoforms modified with +61 Da were identified at level 3.
To reach level 1 of proteoform identification for phosphorylated proteoforms, we analyzed those modified PrSMs according to PS and TP. We took into account PrSMs with a 2A classification level that matched our cluster assignment and filtered those with a high localization confidence level (PTM localization probability > 75% in TP and C-Score > 40 in PS [58].). A total of 111 p-sites passed the filter (33% of the total), 63 for PS and 48 for TP. We manually checked two spectra for P1 +ph and P2 +2ph and their localization sites to be confident of the automatic assignment by the softwares (Figure 5, panel a). Both nodes found P1 and pre-P2 phosphorylated (P1 monophosphorylated and pre-P2 with one, two, three and four p-sites). PS and TP mainly localized the phosphorylation of P1 in site S-11 with 20 and 23 PrSMs, respectively. TP also found PrSMs with the phosphorylation site in S-13, Y-16, S-22 and Y-44. While S-13 was identified with 11 PrSMs, the other three sites were identified with only one PrSM each. However, PS and TP showed a distinct behavior when localizing p-sites in pre-P2. TP found pre-P2 monophosphorylated in S-59, S-13 and S-51 and PS found pre-P2 bi-phosphorylated in S-8 and S-10, tri-phosphorylated in Y-4, S-8 and S-10 and tetra-phosphorylated in Y-4, S-8, S-10 and T-98 (Figure 5, panel b). Given the above observations, we concluded that TP tends to localize mono-phosphorylations while PS performs better at localizing multiple sites.
Protamine phosphorylation sites have been studied in humans and other mammalian species through conventional approaches [15,59,60,61,62]. For instance, protamine 1 phosphorylation at S11 has been described in humans and is conserved among other species, as well as at S9 and S29 residues [24,63,64,65]. Regarding pre-P2, phosphorylation sites have been identified at S8, S10, S37, T47, S51, S59 and S73 [24,63,64,65]. Our results not only are consistent with these previous evidences, but also have allowed the identification of novel proteoforms in the human sperm [31]. These results support this new DBSCAN-based strategy combining PS and TP as a robust approach for identifying site-specific mono- and poly-phosphorylated human protamines.

4. Conclusions

Our clustering approach combined the strengths, and overcame the weaknesses, of two popular top-down proteomics bioinformatics tools: ProSight PD (PS) and TopPIC Suite (TP). We showed the usefulness of the DBSCAN algorithm to solve some of the misassignments found when analyzing top-down database searches provided by these two nodes. From a qualitative viewpoint, and despite the limited number of biological replicates available (n = 2), our proteoform characterization showed good reproducibility across individuals since only a level 3 proteoform was not identified in both biological replicates, namely, P1 +2ph +2(61). With this strategy, we identified 13 proteoforms of P1 and 19 proteoforms of P2, 11 of them with the p-site localized. Our results highlighted the capacity of our clustering approach to characterize purified small proteins with a high number of truncated forms or PTMs. The possibility of examining alterations through studying individual protamine proteoforms, including the assessment of specific modified proteoforms, is a first step towards a personalized diagnosis.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/proteomes9020021/s1, Figure S1: Scatterplots of the intensity vs. the deconvoluted experimental mass and the retention time for PS and TP by cluster, Figure S2: Scatterplots of the retention time vs. the deconvoluted experimental mass by proteoform ID. Figure S3: Boxplot of -log10 E-value and z-score of log2 intensity vs. Node by outcome. Table S1: Protamine proteform cluster annotation table. Table S2: Harmonized PrSMs dataset including clustering outcome and proteform annotation. The Python code we implemented to perform the data analysis and to generate some figures (Figures 3–5; Supplementary Figures S1–S3) and tables appearing on this paper (Table 1; Supplementary Table S1 and S2) is available on the following GitHub repository: https://github.com/MSPCF/2021_Proteomes.

Author Contributions

Conceptualization, G.A.-G., R.O., M.V. (Marta Vilaseca), and M.G.; methodology, A.S.-V., J.C., A.d.l.I.R., M.J., and M.V. (Mar Vilanova); software, G.A.-G. and M.G.; validation, G.A.-G. and M.G.; formal analysis, G.A.-G. and M.G.; investigation, A.S.-V., J.C., A.d.l.I.R., M.J., and M.V. (Mar Vilanova); resources, R.O. and M.V. (Marta Vilaseca); data curation, G.A.-G. and M.G.; writing—original draft preparation, G.A.-G. and M.G.; writing—review and editing, G.A.-G., M.J., M.V. (Mar Vilanova), A.d.l.I.R., J.C., A.S.-V., R.O., M.V. (Marta Vilaseca), and M.G.; visualization, G.A.-G. and M.G.; supervision, M.J., M.G.; project administration, R.O. and M.V. (Marta Vilaseca); funding acquisition, R.O. and M.V. (Marta Vilaseca). All authors have read and agreed to the published version of the manuscript.

Funding

IRB Barcelona Mass Spectrometry and Proteomics Core Facility is a principal unit of Proteored, PRB3-ISCIII, supported by grant PRB3 (IPT17/0019 - ISCIII-SGEFI/ERDF. GA was cofinanced from 01/01/2018-31/12/2021 with grant reference: PT17/0019/0022 (Plataforma de Proteómica, Genotipado y Líneas celulares, PRB3) in the frame of the “Plataforma de apoyo a la Investigación en Ciencias y Tecnología de la Salud de la convocatoria de 2017 de la Acción Estratégica en Salud 2013-2016 of Instituto de Salud Carlos III (ISCIII)”. The Proteomics unit is granted for the set up and development of the Biomedical Proteomics Platform (IRB Barcelona’s Unit) for biomedical and translational Research. Co-financed by the European Regional Development Fund (ERDF), in the framework of the 2014–2020 ERDF Operational Programme in Catalonia, with EUR 519,500. Reference: IU16-015983. Supported by the Spanish Ministerio de Ciencia e Innovación; fondos FEDER ‘una manera de hacer Europa’, Proyectos de Investigación en Salud (ISCIII), PI16/00346 and PI20/00936 to R.O, and by Sara Borrell Postdoctoral Fellowship, Acción Estratégica en Salud, CD17/00109 to J.C.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Clinical Research Ethics Committee of the Hospital Clínic de Barcelona (Barcelona, Spain, protocol approval code 2012/7942).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [39] partner repository with the dataset identifier PXD024405.

Acknowledgments

We would like to thank the whole developing team behind TopPIC Suite. We are particularly grateful to Xiaowen Liu and In Kwon Choi for their extremely efficient support. We thank Mick Greer, Joseph Greer and Neil Kelleher for their support with ProSight PD. Thanks to Tanya Yates from the IRB Editorial support. Finally, we are grateful to the whole IRB Mass Spec & Proteomics team for their professionality and effort.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oliva, R.; Dixon, G.H. Vertebrate protamine genes and the histone-to-protamine replacement reaction. Prog. Nucleic Acid Res. Mol. Biol. 1991, 40, 25–94. [Google Scholar] [CrossRef] [PubMed]
  2. Balhorn, R. The protamine family of sperm nuclear proteins. Genome Biol. 2007, 8, 227. [Google Scholar] [CrossRef]
  3. Castillo, J.; Simon, L.; de Mateo, S.; Lewis, S.; Oliva, R. Protamine/DNA ratios and DNA damage in native and density gradient centrifuged sperm from infertile patients. J. Androl. 2011, 32, 324–332. [Google Scholar] [CrossRef] [Green Version]
  4. Oliva, R. Protamines and male infertility. Hum. Reprod. Update 2006, 12, 417–435. [Google Scholar] [CrossRef] [Green Version]
  5. Mezquita, C. Chromatin composition, structure and function in spermatogenesis. Revis. Biol. Celular 1985, 5, 1–124. [Google Scholar]
  6. Kimmins, S.; Sassone-Corsi, P. Chromatin remodelling and epigenetic features of germ cells. Nature 2005, 434, 583–589. [Google Scholar] [CrossRef]
  7. Oliva, R.; Krawetz, S.A.; GHD Consortium. Gordon Henry Dixon. 25 March 1930–24 July 2016. Biogr. Mems Fell. R. Soc. 2020. [Google Scholar] [CrossRef]
  8. Hud, N.V.; Downing, K.H.; Balhorn, R. A constant radius of curvature model for the organization of DNA in toroidal condensates. Proc. Natl. Acad. Sci. USA 1995, 92, 3581–3585. [Google Scholar] [CrossRef] [Green Version]
  9. Ward, W.S. Function of sperm chromatin structural elements in fertilization and development. Mol. Hum. Reprod. 2010, 16, 30–36. [Google Scholar] [CrossRef] [Green Version]
  10. Balhorn, R.; Corzett, M.; Mazrimas, J.A. Formation of intraprotamine disulfides in vitro. Arch. Biochem. Biophys. 1992, 296, 384–393. [Google Scholar] [CrossRef]
  11. Björndahl, L.; Kvist, U. Human sperm chromatin stabilization: A proposed model including zinc bridges. Mol. Hum. Reprod. 2010, 16, 23–29. [Google Scholar] [CrossRef] [Green Version]
  12. Oliva, R.; Castillo, J. Proteomics and the genetics of sperm chromatin condensation. Asian J. Androl. 2011, 13, 24–30. [Google Scholar] [CrossRef] [Green Version]
  13. Zini, A.; Agarwal, A. Sperm Chromatin: Biological and Clinical Applications in Male Infertility and Assisted Reproduction; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; ISBN 9781441968579. [Google Scholar]
  14. Lüke, L.; Tourmente, M.; Roldan, E.R.S. Sexual Selection of Protamine 1 in Mammals. Mol. Biol. Evol. 2016, 33, 174–184. [Google Scholar] [CrossRef] [Green Version]
  15. Kasinsky, H.E.; Eirín-López, J.M.; Ausió, J. Protamines: Structural complexity, evolution and chromatin patterning. Protein Pept. Lett. 2011, 18, 755–771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Castillo, J.; Estanyol, J.M.; Ballescá, J.L.; Oliva, R. Human sperm chromatin epigenetic potential: Genomics, proteomics, and male infertility. Asian J. Androl. 2015, 17, 601–609. [Google Scholar] [CrossRef] [PubMed]
  17. Castillo, J.; Jodar, M.; Oliva, R. The contribution of human sperm proteins to the development and epigenome of the preimplantation embryo. Hum. Reprod. Update 2018, 24, 535–555. [Google Scholar] [CrossRef] [Green Version]
  18. Gatewood, J.M.; Cook, G.R.; Balhorn, R.; Bradbury, E.M.; Schmid, C.W. Sequence-specific packaging of DNA in human sperm chromatin. Science 1987, 236, 962–964. [Google Scholar] [CrossRef] [PubMed]
  19. Arpanahi, A.; Brinkworth, M.; Iles, D.; Krawetz, S.A.; Paradowska, A.; Platts, A.E.; Saida, M.; Steger, K.; Tedder, P.; Miller, D. Endonuclease-sensitive regions of human spermatozoal chromatin are highly enriched in promoter and CTCF binding sequences. Genome Res. 2009, 19, 1338–1349. [Google Scholar] [CrossRef] [Green Version]
  20. Hammoud, S.S.; Nix, D.A.; Zhang, H.; Purwar, J.; Carrell, D.T.; Cairns, B.R. Distinctive chromatin in human sperm packages genes for embryo development. Nature 2009, 460, 473–478. [Google Scholar] [CrossRef] [Green Version]
  21. Castillo, J.; Amaral, A.; Azpiazu, R.; Vavouri, T.; Estanyol, J.M.; Ballescà, J.L.; Oliva, R. Genomic and proteomic dissection and characterization of the human sperm chromatin. Mol. Hum. Reprod. 2014, 20, 1041–1053. [Google Scholar] [CrossRef] [Green Version]
  22. Jodar, M.; Oliva, R. Protamine alterations in human spermatozoa. Adv. Exp. Med. Biol. 2014, 791, 83–102. [Google Scholar] [CrossRef] [PubMed]
  23. Ni, K.; Spiess, A.-N.; Schuppe, H.-C.; Steger, K. The impact of sperm protamine deficiency and sperm DNA damage on human male fertility: A systematic review and meta-analysis. Andrology 2016, 4, 789–799. [Google Scholar] [CrossRef]
  24. Barrachina, F.; Soler-Ventura, A.; Oliva, R.; Jodar, M. Sperm Nucleoproteins (Histones and Protamines). In A Clinician’s Guide to Sperm DNA and Chromatin Damage; Springer: Cham, Switzerland, 2018; pp. 31–51. [Google Scholar]
  25. Nanassy, L.; Liu, L.; Griffin, J.; Carrell, D.T. The clinical utility of the protamine 1/protamine 2 ratio in sperm. Protein Pept. Lett. 2011, 18, 772–777. [Google Scholar] [CrossRef]
  26. De Yebra, L.; Oliva, R. Rapid analysis of mammalian sperm nuclear proteins. Anal. Biochem. 1993, 209, 201–203. [Google Scholar] [CrossRef]
  27. Soler-Ventura, A.; Castillo, J.; de la Iglesia, A.; Jodar, M.; Barrachina, F.; Ballesca, J.L.; Oliva, R. Mammalian Sperm Protamine Extraction and Analysis: A Step-By-Step Detailed Protocol and Brief Review of Protamine Alterations. Protein Pept. Lett. 2018, 25, 424–433. [Google Scholar] [CrossRef]
  28. Bogle, O.A.; Kumar, K.; Attardo-Parrinello, C.; Lewis, S.E.M.; Estanyol, J.M.; Ballescà, J.L.; Oliva, R. Identification of protein changes in human spermatozoa throughout the cryopreservation process. Andrology 2017, 5, 10–22. [Google Scholar] [CrossRef] [PubMed]
  29. De Mateo, S.; Ramos, L.; van der Vlag, J.; de Boer, P.; Oliva, R. Improvement in chromatin maturity of human spermatozoa selected through density gradient centrifugation. Int. J. Androl. 2011, 34, 256–267. [Google Scholar] [CrossRef]
  30. Wang, T.; Gao, H.; Li, W.; Liu, C. Essential Role of Histone Replacement and Modifications in Male Fertility. Front. Genet. 2019, 10, 962. [Google Scholar] [CrossRef]
  31. Soler-Ventura, A.; Gay, M.; Jodar, M.; Vilanova, M.; Castillo, J.; Arauz-Garofalo, G.; Villarreal, L.; Ballescà, J.L.; Vilaseca, M.; Oliva, R. Characterization of Human Sperm Protamine Proteoforms through a Combination of Top-Down and Bottom-Up Mass Spectrometry Approaches. J. Proteome Res. 2020, 19, 221–237. [Google Scholar] [CrossRef]
  32. LeDuc, R.D.; Taylor, G.K.; Kim, Y.-B.; Januszyk, T.E.; Bynum, L.H.; Sola, J.V.; Garavelli, J.S.; Kelleher, N.L. ProSight PTM: An integrated environment for protein identification and characterization by top-down mass spectrometry. Nucleic Acids Res. 2004, 32, W340–W345. [Google Scholar] [CrossRef] [PubMed]
  33. Frank, A.M.; Pesavento, J.J.; Mizzen, C.A.; Kelleher, N.L.; Pevzner, P.A. Interpreting top-down mass spectra using spectral alignment. Anal. Chem. 2008, 80, 2499–2505. [Google Scholar] [CrossRef] [PubMed]
  34. Kou, Q.; Xun, L.; Liu, X. TopPIC: A software tool for top-down mass spectrometry-based proteoform identification and characterization. Bioinformatics 2016, 32, 3495–3497. [Google Scholar] [CrossRef] [PubMed]
  35. Kou, Q.; Wu, S.; Tolic, N.; Paša-Tolic, L.; Liu, Y.; Liu, X. A mass graph-based approach for the identification of modified proteoforms using top-down tandem mass spectra. Bioinformatics 2017, 33, 1309–1316. [Google Scholar] [CrossRef] [Green Version]
  36. Cooper, T.G.; Noonan, E.; von Eckardstein, S.; Auger, J.; Baker, H.W.G.; Behre, H.M.; Haugen, T.B.; Kruger, T.; Wang, C.; Mbizvo, M.T.; et al. World Health Organization reference values for human semen characteristics. Hum. Reprod. Update 2010, 16, 231–245. [Google Scholar] [CrossRef] [PubMed]
  37. Mengual, L.; Ballescá, J.L.; Ascaso, C.; Oliva, R. Marked differences in protamine content and P1/P2 ratios in sperm cells from percoll fractions between patients and controls. J. Androl. 2003, 24, 438–447. [Google Scholar] [CrossRef]
  38. Aebersold, R.; Agar, J.N.; Jonathan Amster, I.; Baker, M.S.; Bertozzi, C.R.; Boja, E.S.; Costello, C.E.; Cravatt, B.F.; Fenselau, C.; Garcia, B.A.; et al. How many human proteoforms are there? Nat. Chem. Biol. 2018, 14, 206–214. [Google Scholar] [CrossRef] [Green Version]
  39. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef]
  40. 3.9.1 Documentation. Available online: https://docs.python.org/3/ (accessed on 30 December 2020).
  41. McKinney, W. Data Structures for Statistical Computing in Python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 3–8 July 2010; pp. 56–61. [Google Scholar]
  42. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef]
  43. Cock, P.J.A.; Antao, T.; Chang, J.T.; Chapman, B.A.; Cox, C.J.; Dalke, A.; Friedberg, I.; Hamelryck, T.; Kauff, F.; Wilczynski, B.; et al. Biopython: Freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 2009, 25, 1422–1423. [Google Scholar] [CrossRef]
  44. Kibirige, H.; Lamp, G.; Katins, J.; Gdowding; Austin, O.; Matthias, K.; Funnell, T.; Finkernagel, F.; Arnfred, J.; Blanchard, D.; et al. has2k1/Plotnine; v0.7.1; Zenodo: Genève, Switzerland, 2020. [Google Scholar] [CrossRef]
  45. Waskom, M.; Gelbart, M.; Botvinnik, O.; Ostblom, J.; Hobson, P.; Lukauskas, S.; Gemperline, D.C.; Augspurger, T.; Halchenko, Y.; Warmenhoven, J.; et al. Mwaskom/Seaborn; v0.11.1; Zenodo: Genève, Switzerland, 2020. [Google Scholar] [CrossRef]
  46. Schubert, E.; Sander, J.; Ester, M.; Kriegel, H.P.; Xu, X. DBSCAN Revisited, Revisited. ACM Trans. Database Syst. 2017, 42, 1–21. [Google Scholar] [CrossRef]
  47. Ester, M.; Kriegel, H.-P.; Sander, J.; Xu, X. A Density-Based Algorithm for Discovering Clusters in Large Spatial Databases with Noise; Association for the Advancement of Artificial Intelligence: Palo Alto, CA, USA, 1996. [Google Scholar]
  48. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  49. Van de Ven, B. Bokeh. Available online: https://bokeh.org/ (accessed on 29 January 2021).
  50. Smith, L.M.; Thomas, P.M.; Shortreed, M.R.; Schaffer, L.V.; Fellers, R.T.; LeDuc, R.D.; Tucholski, T.; Ge, Y.; Agar, J.N.; Anderson, L.C.; et al. A five-level classification system for proteoform identifications. Nat. Methods 2019, 16, 939–940. [Google Scholar] [CrossRef] [PubMed]
  51. Bianchi, F.; Rousseaux-Prevost, R.; Sautiere, P.; Rousseaux, J. P2 protamines from human sperm are zinc -finger proteins with one CYS2/HIS2 motif. Biochem. Biophys. Res. Commun. 1992, 182, 540–547. [Google Scholar] [CrossRef]
  52. Roskamp, K.W.; Azim, S.; Kassier, G.; Norton-Baker, B.; Sprague-Piercy, M.A.; Miller, R.J.D.; Martin, R.W. Human γS-Crystallin-Copper Binding Helps Buffer against Aggregation Caused by Oxidative Damage. Biochemistry 2020, 59, 2371–2385. [Google Scholar] [CrossRef]
  53. Liu, H.; Ponniah, G.; Zhang, H.-M.; Nowak, C.; Neill, A.; Gonzalez-Lopez, N.; Patel, R.; Cheng, G.; Kita, A.Z.; Andrien, B. In vitro and in vivo modifications of recombinant and human IgG antibodies. MAbs 2014, 6, 1145–1154. [Google Scholar] [CrossRef] [Green Version]
  54. Fellers, R.T.; Greer, J.B.; Early, B.P.; Yu, X.; LeDuc, R.D.; Kelleher, N.L.; Thomas, P.M. ProSight Lite: Graphical software to analyze top-down mass spectrometry data. Proteomics 2015, 15, 1235–1238. [Google Scholar] [CrossRef] [PubMed]
  55. De Yebra, L.; Ballescà, J.L.; Vanrell, J.A.; Bassas, L.; Oliva, R. Complete selective absence of protamine P2 in humans. J. Biol. Chem. 1993, 268, 10553–10557. [Google Scholar] [CrossRef]
  56. Torregrosa, N.; Domínguez-Fandos, D.; Camejo, M.I.; Shirley, C.R.; Meistrich, M.L.; Ballescà, J.L.; Oliva, R. Protamine 2 precursors, protamine 1/protamine 2 ratio, DNA integrity and other sperm parameters in infertile patients. Hum. Reprod. 2006, 21, 2084–2089. [Google Scholar] [CrossRef] [PubMed]
  57. De Mateo, S.; Gázquez, C.; Guimerà, M.; Balasch, J.; Meistrich, M.L.; Ballescà, J.L.; Oliva, R. Protamine 2 precursors (Pre-P2), protamine 1 to protamine 2 ratio (P1/P2), and assisted reproduction outcome. Fertil. Steril. 2009, 91, 715–722. [Google Scholar] [CrossRef]
  58. LeDuc, R.D.; Fellers, R.T.; Early, B.P.; Greer, J.B.; Thomas, P.M.; Kelleher, N.L. The C-score: A Bayesian framework to sharply improve proteoform scoring in high-throughput top down proteomics. J. Proteome Res. 2014, 13, 3231–3240. [Google Scholar] [CrossRef]
  59. Gusse, M.; Sautière, P.; Bélaiche, D.; Martinage, A.; Roux, C.; Dadoune, J.P.; Chevaillier, P. Purification and characterization of nuclear basic proteins of human sperm. Biochim. Biophys. Acta 1986, 884, 124–134. [Google Scholar] [CrossRef]
  60. Pruslin, F.H.; Imesch, E.; Winston, R.; Rodman, T.C. Phosphorylation state of protamines 1 and 2 in human spermatids and spermatozoa. Gamete Res. 1987, 18, 179–190. [Google Scholar] [CrossRef]
  61. Arkhis, A.; Martinage, A.; Sautiere, P.; Chevaillier, P. Molecular structure of human protamine P4 (HP4), a minor basic protein of human sperm nuclei. Eur. J. Biochem. 1991, 200, 387–392. [Google Scholar] [CrossRef] [PubMed]
  62. Papoutsopoulou, S.; Nikolakaki, E.; Chalepakis, G.; Kruft, V.; Chevaillier, P.; Giannakouros, T. SR protein-specific kinase 1 is highly expressed in testis and phosphorylates protamine 1. Nucleic Acids Res. 1999, 27, 2972–2980. [Google Scholar] [CrossRef] [Green Version]
  63. Chirat, F.; Arkhis, A.; Martinage, A.; Jaquinod, M.; Chevaillier, P.; Sautière, P. Phosphorylation of human sperm protamines HP1 and HP2: Identification of phosphorylation sites. Biochim. Biophys. Acta 1993, 1203, 109–114. [Google Scholar] [CrossRef]
  64. Pirhonen, A.; Linnala-Kankkunen, A.; Mäenpää, P.H. Identification of phosphoseryl residues in protamines from mature mammalian spermatozoa. Biol. Reprod. 1994, 50, 981–986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Hornbeck, P.V.; Zhang, B.; Murray, B.; Kornhauser, J.M.; Latham, V.; Skrzypek, E. PhosphoSitePlus, 2014: Mutations, PTMs and recalibrations. Nucleic Acids Res. 2015, 43, D512–D520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Human protamine amino acid sequences. Protamine 1 is translated as a mature form containing 51 residues (P1). In contrast, protamine 2 is translated as a precursor protein containing 102 residues (pre-P2), which is processed by proteolysis in multiple intermediate forms (HPI2 (22–102), HPS1 (34–102) and HPS2 (37–102)) and mature forms (HP2 (46–102), HP3 (49–102) and HP4 (45–102)).
Figure 1. Human protamine amino acid sequences. Protamine 1 is translated as a mature form containing 51 residues (P1). In contrast, protamine 2 is translated as a precursor protein containing 102 residues (pre-P2), which is processed by proteolysis in multiple intermediate forms (HPI2 (22–102), HPS1 (34–102) and HPS2 (37–102)) and mature forms (HP2 (46–102), HP3 (49–102) and HP4 (45–102)).
Proteomes 09 00021 g001
Figure 2. Protamine characterization workflow. We purified the protamine content of two biological replicates (A, B) and injected four technical replicates (R01, R02, R03, R04) into the LC-MS instrument. We then used two bioinformatic tools (ProSight PD (Thermo Scientific, San Jose, CA, USA) and TopPIC Suite (Indianapolis, IN, USA)) to identify proteoform spectrum matches. The whole data analysis pipeline downstream was implemented entirely using Python as a programming language.
Figure 2. Protamine characterization workflow. We purified the protamine content of two biological replicates (A, B) and injected four technical replicates (R01, R02, R03, R04) into the LC-MS instrument. We then used two bioinformatic tools (ProSight PD (Thermo Scientific, San Jose, CA, USA) and TopPIC Suite (Indianapolis, IN, USA)) to identify proteoform spectrum matches. The whole data analysis pipeline downstream was implemented entirely using Python as a programming language.
Proteomes 09 00021 g002
Figure 3. Scatterplots of the retention time (RT) vs. the deconvoluted experimental mass (m) for PS (Prosight PD) and TP (TopPIC Suite). (a) Deconvoluted masses from 6000 Da to 14,000 Da. The sequential color scale highlights the clusters obtained sorted by increasing experimental mass, while unclustered PrSMs appear in grey. (b,c) Deconvoluted masses from 8900 Da to 9400 Da corresponding to clusters that we associated with HPS1 and HPS2 proteoforms (the region highlighted with a rectangle in panel (a)). The color scale in panel (b) denotes each proteoform according to PS and TP, while the color scale in panel (c) denotes each proteoform according to our clustering approach. (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid).
Figure 3. Scatterplots of the retention time (RT) vs. the deconvoluted experimental mass (m) for PS (Prosight PD) and TP (TopPIC Suite). (a) Deconvoluted masses from 6000 Da to 14,000 Da. The sequential color scale highlights the clusters obtained sorted by increasing experimental mass, while unclustered PrSMs appear in grey. (b,c) Deconvoluted masses from 8900 Da to 9400 Da corresponding to clusters that we associated with HPS1 and HPS2 proteoforms (the region highlighted with a rectangle in panel (a)). The color scale in panel (b) denotes each proteoform according to PS and TP, while the color scale in panel (c) denotes each proteoform according to our clustering approach. (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid).
Proteomes 09 00021 g003aProteomes 09 00021 g003b
Figure 4. Protamine proteoform landscapes and ratios. (ac) Scatterplot of the mass shift vs. the protamine form at global, biological replicate and technical replicate levels, respectively. The sequential color scale in these three scatterplots highlights the particular set of PTMs. The data point size is proportional to the number of PrSMs hitting the corresponding protamine proteoform. (df) Barplots with P1/P2-mature and P2-immature/P2-mature protamine ratios according to each node at global, biological replicate and technical replicate levels, respectively. Error bars correspond to standard deviations of associated technical replicates (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid).
Figure 4. Protamine proteoform landscapes and ratios. (ac) Scatterplot of the mass shift vs. the protamine form at global, biological replicate and technical replicate levels, respectively. The sequential color scale in these three scatterplots highlights the particular set of PTMs. The data point size is proportional to the number of PrSMs hitting the corresponding protamine proteoform. (df) Barplots with P1/P2-mature and P2-immature/P2-mature protamine ratios according to each node at global, biological replicate and technical replicate levels, respectively. Error bars correspond to standard deviations of associated technical replicates (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid).
Proteomes 09 00021 g004
Figure 5. Protamine proteoform phosphosite localization analysis. (a) MS2 deconvoluted spectra and their fragment map of P1 +ph (up) and P2 +2ph (down) showing all detected c and z ions and p-site localizations (phosphorylated amino acids are marked in blue). (b) Barplots report the number of PrSMs hitting the corresponding localization site. The faceting grid breaks down the phosphosite count for P1 and pre-P2 protamine forms (left and right), as well as for each node (up and down). Color scale denotes the number of phosphorylations coexisting in the same sequence (ph = phosphorylation).
Figure 5. Protamine proteoform phosphosite localization analysis. (a) MS2 deconvoluted spectra and their fragment map of P1 +ph (up) and P2 +2ph (down) showing all detected c and z ions and p-site localizations (phosphorylated amino acids are marked in blue). (b) Barplots report the number of PrSMs hitting the corresponding localization site. The faceting grid breaks down the phosphosite count for P1 and pre-P2 protamine forms (left and right), as well as for each node (up and down). Color scale denotes the number of phosphorylations coexisting in the same sequence (ph = phosphorylation).
Proteomes 09 00021 g005
Table 1. Cluster annotation summary. We provide the cluster-to-proteoform correspondence in columns “Cluster” and “Proteoform”. Column “mT (Da)” gives the theoretical mass of the fully carbamidomethylated cluster-protein (for HP4 the mass is 8062.32 Da). Columns “mC (Da) (mean)” and “mC (Da) (std)” give the mean mass and the mass standard deviation of all PrSMs comprising each cluster. Column “ΔmCT (Da)” gives the cluster-protein mass shift returned by Equation (2). We also report proteoform and PrSM counts according to PS and TP nodes in columns “PF count” and “PrSM count”, respectively (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid, PF = proteoform, PrSMs = proteoform spectrum matches, ProSight = PS, TopPIC = TP).
Table 1. Cluster annotation summary. We provide the cluster-to-proteoform correspondence in columns “Cluster” and “Proteoform”. Column “mT (Da)” gives the theoretical mass of the fully carbamidomethylated cluster-protein (for HP4 the mass is 8062.32 Da). Columns “mC (Da) (mean)” and “mC (Da) (std)” give the mean mass and the mass standard deviation of all PrSMs comprising each cluster. Column “ΔmCT (Da)” gives the cluster-protein mass shift returned by Equation (2). We also report proteoform and PrSM counts according to PS and TP nodes in columns “PF count” and “PrSM count”, respectively (ph = phosphorylation, Q = glutamine, pyroGlu = pyroglutamic acid, PF = proteoform, PrSMs = proteoform spectrum matches, ProSight = PS, TopPIC = TP).
ClusterProteoformmT (Da)mC (Da)(mean)mC (Da)(std)ΔmCT (Da)PF CountPrSM Count
PSTPPSTP
P|000P17029.67029.451.39−0.1425247206
P|001P1 + 1(61)-7090.921.9261.325415697
P|002P1 + 1ph-7108.861.5879.27357174
P|003P1 + 2(61)-7151.591.56122637146
P|004P1 + 1ph + 1(61)-7170.391.67140.79235132
P|005P1 + 3(61)-7212.911.85183.32436441
P|006P1 + 1ph + 2(61)-7231.821.66202.23622920
P|007P1 + 1ph + 3(61)-7293.522.73263.92454229
P|008P1 + 2ph + 2(61)-7310.921.43281.323275
P|009HP37539.077538.541.26−0.53227253
P|010HP3 + 1(61)-7599.631.2060.56121713
P|011HP27933.287932.731.08−0.5514142119
P|012HP2 + 1(61)-7994.031.3560.75116451
P|013HP4 + 1(61)-8122.660.4960.342 16
P|014HPS29002.759002.091.10−0.66236954
P|015HPS2 + 1(61)-9063.41.9960.66233110
P|016HPS1 − 1Q-9172.081.09−128.83312319
P|017HPS1 +1Pyroglu-9283.321.27−17.59932621
P|018HPS19300.919301.171.070.25137970
P|019HPS1 + 1(61)-9363.591.5162.67322318
P|020HPI210,654.4610,653.971.01−0.491166
P|021pre-P213,196.8213,195.971.67−0.85145851
P|022pre-P2 + 1ph-13,275.941.8079.12342722
P|023pre-P2 + 2ph-13,355.921.54159.1152829
P|024pre-P2 + 3ph-13,435.641.59238.82331417
P|025pre-P2 + 4ph-13,515.51.19318.682286
P|026Unknown-13,873.220.78- 2 28
UnclusteredUnknown----578916694
Table 2. Converging and diverging missassigments examples. Row (a) show an example of converging missassigment and row (b) its corresponding manual reassignation after cluster inspection: HPS2 +1(61) proteoform is more likely than HPS2. Rows (c) and (d) along with (e) and (f) provide two examples of diverging missassigments and their associated corrections. For both examples, HPS1 +1pyroGlu is more likely than HPS2 +2ph or HPS1 +1AcNterm, respectively. Check ProSight Lite fragmentation pattern, mass shift and remaining scoring columns for details (ph = phosphorylation, pyroGlu = pyroglutamic acid). Columns “me” and “mt” give the experimental and theoretical masses, respectively. “Δm” gives the difference between me and mt. “P-Score” is the score as provided by ProSight Lite. “%Res. Cleav. is the percentage of inter-residues observed. “Match. Frag.” are the number of assigned matching fragments.
Table 2. Converging and diverging missassigments examples. Row (a) show an example of converging missassigment and row (b) its corresponding manual reassignation after cluster inspection: HPS2 +1(61) proteoform is more likely than HPS2. Rows (c) and (d) along with (e) and (f) provide two examples of diverging missassigments and their associated corrections. For both examples, HPS1 +1pyroGlu is more likely than HPS2 +2ph or HPS1 +1AcNterm, respectively. Check ProSight Lite fragmentation pattern, mass shift and remaining scoring columns for details (ph = phosphorylation, pyroGlu = pyroglutamic acid). Columns “me” and “mt” give the experimental and theoretical masses, respectively. “Δm” gives the difference between me and mt. “P-Score” is the score as provided by ProSight Lite. “%Res. Cleav. is the percentage of inter-residues observed. “Match. Frag.” are the number of assigned matching fragments.
ProSight Lite Fragmentation PatternFormMod.me (Da)mt (Da)Δm
(Da)
P-Score%Res. Cleav.Match Frag.
(a) Originally assigned by PS
Proteomes 09 00021 i001
BioRep A—TechRep 02—Scan 349
HPS2-9063.839002.7561.0841.10 × 10−917462
(b) Manual reassignation
Proteomes 09 00021 i002
BioRep A—TechRep R02—Scan 349
HPS2+1(61)9063.839063.760.0767.60·× 10−1558295
(c) Originally assigned by PS
Proteomes 09 00021 i003
BioRep A—TechRep R02—Scan 539
HPS2+2ph9282.899162.68120.2137.90·× 10−385235
(d) Manual reassignation
Proteomes 09 00021 i004
BioRep A—TechRep R02—Scan 539
HPS1+1pyroGlu9282.899283.89−0.9922.80·× 10−1037475
(e) Originally assigned by PS
Proteomes 09 00021 i005
BioRep B—TechRep R04—Scan 545
HPS1+1ph
+1AcNterm
9283.879422.89−139.0229.70·× 10−455940
(f) Manual reassignation
Proteomes 09 00021 i006
BioRep B—TechRep R04—Scan 545
HPS1+1pyroGlu9283.879283.89−0.0191.18·× 10−1187984
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arauz-Garofalo, G.; Jodar, M.; Vilanova, M.; de la Iglesia Rodriguez, A.; Castillo, J.; Soler-Ventura, A.; Oliva, R.; Vilaseca, M.; Gay, M. Protamine Characterization by Top-Down Proteomics: Boosting Proteoform Identification with DBSCAN. Proteomes 2021, 9, 21. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020021

AMA Style

Arauz-Garofalo G, Jodar M, Vilanova M, de la Iglesia Rodriguez A, Castillo J, Soler-Ventura A, Oliva R, Vilaseca M, Gay M. Protamine Characterization by Top-Down Proteomics: Boosting Proteoform Identification with DBSCAN. Proteomes. 2021; 9(2):21. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020021

Chicago/Turabian Style

Arauz-Garofalo, Gianluca, Meritxell Jodar, Mar Vilanova, Alberto de la Iglesia Rodriguez, Judit Castillo, Ada Soler-Ventura, Rafael Oliva, Marta Vilaseca, and Marina Gay. 2021. "Protamine Characterization by Top-Down Proteomics: Boosting Proteoform Identification with DBSCAN" Proteomes 9, no. 2: 21. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop