Next Article in Journal
Yes-Associated Protein (Yap) Is Up-Regulated in Heart Failure and Promotes Cardiac Fibroblast Proliferation
Previous Article in Journal
Diagnostics and Therapeutics in Targeting HER2 Breast Cancer: A Novel Approach
Previous Article in Special Issue
Muscle Proteomic Profile before and after Enzyme Replacement Therapy in Late-Onset Pompe Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plasma Proteomic Analysis in Morquio A Disease

by
José V. Álvarez
1,2,3,
Susana B. Bravo
4,
María Pilar Chantada-Vázquez
4,
Sofía Barbosa-Gouveia
2,
Cristóbal Colón
1,2,
Olalla López-Suarez
1,2,
Shunji Tomatsu
3,
Francisco J. Otero-Espinar
5,6,* and
María L. Couce
1,2,*
1
Department of Forensic Sciences, Pathology, Gynecology and Obstetrics, Pediatrics, Neonatology Service, Health Research Institute of Santiago de Compostela (IDIS), Hospital Clínico Universitario de Santiago de Compostela, CIBERER, MetabERN, 15706 Santiago de Compostela, Spain
2
Health Research Institute of Santiago de Compostela (IDIS), CIBERER, MetabERN, 15706 Santiago de Compostela, Spain
3
Skeletal Dysplasia Lab Nemours Biomedical Research Nemours/Alfred I, duPont Hospital for Children, 1600 Rockland Road, Wilmington, DE 19803, USA
4
Proteomic Platform, Health Research Institute of Santiago de Compostela (IDIS), Hospital ClínicoUniversitario de Santiago de Compostela, 15706 Santiago de Compostela, Spain
5
Paraquasil Platform, Health Research Institute of Santiago de Compostela (IDIS), Hospital Clínico Universitario de Santiago de Compostela, 15706 Santiago de Compostela, Spain
6
Department of Pharmacology, Pharmacy and Pharmaceutical Technology, School of Pharmacy, Campus Vida, University of Santiago de Compostela, 15872 Santiago de Compostela, Spain
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(11), 6165; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22116165
Submission received: 8 May 2021 / Revised: 2 June 2021 / Accepted: 4 June 2021 / Published: 7 June 2021
(This article belongs to the Special Issue Omics Technologies in Rare Diseases)

Abstract

:
Mucopolysaccharidosis type IVA (MPS IVA) is a lysosomal disease caused by mutations in the gene encoding the enzymeN-acetylgalactosamine-6-sulfate sulfatase (GALNS), and is characterized by systemic skeletal dysplasia due to excessive storage of keratan sulfate (KS) and chondroitin-6-sulfate in chondrocytes. Although improvements in the activity of daily living and endurance tests have been achieved with enzyme replacement therapy (ERT) with recombinant human GALNS, recovery of bone lesions and bone growth in MPS IVA has not been demonstrated to date. Moreover, no correlation has been described between therapeutic efficacy and urine levels of KS, which accumulates in MPS IVA patients. The objective of this study was to assess the validity of potential biomarkers proposed by other authors and to identify new biomarkers. To identify candidate biomarkers of this disease, we analyzed plasma samples from healthy controls (n=6) and from untreated (n=8) and ERT-treated (n=5, sampled before and after treatment) MPS IVA patients using both qualitative and quantitative proteomics analyses. The qualitative proteomics approach analyzed the proteomic profile of the different study groups. In the quantitative analysis, we identified/quantified 215 proteins after comparing healthy control untreated, ERT-treated MPSIVA patients. We selected a group of proteins that were dysregulated in MPS IVA patients. We identified four potential protein biomarkers, all of which may influence bone and cartilage metabolism: fetuin-A, vitronectin, alpha-1antitrypsin, and clusterin. Further studies of cartilage and bone samples from MPS IVA patients will be required to verify the validity of these proteins as potential biomarkers of MPS IVA.

1. Introduction

Mucopolysaccharidosis type IVA (MPS IVA, OMIM #253000) or Morquio A syndrome is an autosomal recessive disease caused by mutations in the GALNS gene. This molecular defect results in a deficiency in the enzyme N-acetylgalactosamine-6-sulfate sulfatase (GALNS, E.C.3.1.6.4) [1,2,3], which degrades the specific glycosaminoglycans(GAGs) keratan sulfate (KS) and chondroitin-6-sulfate (C6S) [4,5,6,7]. These GAGs accumulate in multiple tissues, mainly bone, cartilage, heart valves, and cornea, leading to progressive systemic skeletal dysplasia [8,9]. Significant non-skeletal manifestations, including respiratory disease, spinal cord compression, cardiac disease, impaired vision, hearing loss, and dental problems, have also been described in MPS IVA patients [10,11,12,13]. MPS IVA is more frequently associated with severe and extensive skeletal manifestations than the other MPS types. Specifically, hypermobility of the joints is a characteristic of MPS IVA that distinguishes this disease from other forms of MPS. Furthermore, cognitive involvement is generally absent in MPS IVA [14]. Among untreated patients, respiratory failure is the primary cause of death, which typically occurs during the second or third decades of life [15,16].
Currently, two therapies for MPS IVA are used in clinical practice: enzyme replacement therapy (ERT) with recombinant human GALNS; and hematopoietic stem cell transplantation (HSCT) [5,11]. ERT and HSCT are based on the cross-correction principle, where by lysosomal enzymes are taken up by the cells of deficient recipients and targeted to the lysosomes via the mannose-6-phosphate receptor. However, as described for other lysosomal storage diseases, ERT for MPS IVA suffers from several limitations, including the need for weekly infusions of 4–6 h; rapid clearance due to a short half-life (35 min in human, 2 min in mouse) [17,18,19]; a high price [19,20]; limited penetration of the avascular cartilage; and immunological responses against the infused enzyme [16,17,21]. Moreover, clinical trials have shown little improvements in terms of bone growth and skeletal dysplasia. There is no definitive therapy that markedly improves bone and cartilage lesions in MPS IVA [14,18,19,20]. Potential biomarkers proposed in recent studies include KS, C6S, and blood levels of collagen II [22,23]. Although clinical trials have studied urinary KS as a potential biomarker, it has not been demonstrated that decreases in this parameter correlate with clinical improvement. Urinary KS originates in the kidneys but does not reflect the involvement of bone or other relevant tissues in MPS IVA, and is considered a pharmacokinetic marker but not a surrogate biomarker. The identification of valid biomarkers of bone lesions to evaluate therapeutic efficacy thus remains an unmet need. While bone cartilage samples are the best samples in which to search for candidate biomarkers of MPS IVA, the objective of this study was to assess the validity of candidate plasma biomarkers proposed in previous studies and to search for new biomarkers. To this end, we carried out both quantitative and qualitative proteomic analyses of plasma samples from.
Within rare disease research, proteomics is a fast-growing field. Approaches based on liquid chromatography-tandem mass spectrometry (LC-MS/MS), which are used to search for biomarkers in other human diseases [24,25,26,27,28], are key tools used to study normal and affected tissues and body fluids (e.g., plasma) in rare diseases [29].
Disease states involve alterations in the expression levels, phosphorylation states, and post-transcriptional modifications of proteins. Therefore, proteomics techniques have been applied to the study of normal and pathological conditions to search for biological information that could further our understanding of pathophysiological mechanisms and identify novel disease biomarkers [26,28,30]. Specific disease biomarkers identified using proteomics techniques could help improve diagnostic and prognostic precision and serve as key tools for the evaluation of treatment outcomes [24,25,26,27,28,29,30,31]. We recently characterized new proteomic biomarker candidates in leukocytes of MPS IVA patients [32]. However, in the present study we sought to identify disease biomarkers in body fluids that can be more easily sampled and feasibly incorporated into a clinical setting to increase diagnostic accuracy in patients with this disease. We focused on potential biomarkers in plasma, which we identified using two proteomics techniques: qualitative LC-MS/MS and a quantitative assay that compared protein expression between controls and untreated or ERT-treated MPS IVA patients.

2. Results

Plasma samples were acquired from 13 MPS IVA patients and 6 controls and classified as follows: healthy controls; untreated MPS IVA patients; and MPS IVA patients who had been receiving weekly ERT infusions for a mean of 7 years. From these ERT-treated patients, samples were acquired 24 h before (ERT-a) and 24 h after (ERT-b) one of their weekly infusions (Figure 1). The demographic features of the participating MPS IVA patients have been previously described [32].
Qualitative and quantitative proteomic analyses of plasma proteins were performed on all patients and control samples.

2.1. Qualitative Analysis of Proteins in Plasma Samples

In this study, we identified all proteins in plasma samples using data-dependent acquisition (DDA)-LC-MS/MS. Subsequent analysis was limited to proteins identified with 1% error (i.e., those with a false discovery rate (FDR) < 1). To obtain the most representative proteins per group, we selected proteins identified in at least n-1 samples per group (where n is the number of samples per group) (Supplementary Table S1). The Venn diagram in Figure 2 depicts the distribution of proteins in each of the 4 groups, and indicates those that were common to multiple groups.
Once we determined the proteins characteristic of each group, we performed a functional analysis using the FunRich program. This analysis allows us to characterize the behavior of the proteins identified in each group. Comparison of ERT-a and ERT-b samples revealed a favorable evolution (i.e., values obtained in ERT-a or ERT-b were close to thoseobtained for control samples) for certain proteins. Next, we selected groups of proteins in which alterations are observed in MPS IVA, and analyzed their behavior in each study group. From these, we identified proteins involved in hemostasis, bone regeneration, reorganization of theextracellular matrix and collagen, the immune response, complement system, and oxidative stress.
Among proteins implicated in hemostasis (blood coagulation, fibrin clot formation, and blood coagulation intrinsic pathway), we observed that values obtained for ERT-a samples (i.e., those collected 6 days after the preceding enzyme infusion) were closer to those of control samples than to ERT-b samples (Figure 3a). This result suggests that blood levels of these proteins require several days to stabilize after ERT infusion. A similar effect was observed for proteins implicated in different aspects of hemostatic processes, including platelet activation, platelet aggregation, and platelet degranulation (Figure 3b).
Other proteins identified in the proteomic analysis were implicated in bone regeneration and collagen function. Analysis of the behavior of these proteins showed that ERT had a negative regulatory effect on some (including proteins involved in mineralization, ossification, and osteoblast differentiation) and a positive regulatory effect on others (bone resorption and regulation of bone mineralization) (Figure 4a). In the case of proteins that were negatively regulated, ERT required more than 24 h to affect protein expression. Conversely, in the case of positively regulated proteins (involved in bone resorption and osteoblast differentiation) plasma levels increased after treatment in the ERT-b group (Figure 4a).
Proteomic analysis of proteins involved in collagen function revealed values close to those of the control group for proteins related to a positive regulation of the collagen biosynthesis (Figure 4b) in the ERT-a group. Moreover, in the ERT-a study group, we observed decreases in the expression of proteins implicated in collagen biosynthesis and collagen fibril organization. The behavior of proteins involved in metabolic and catabolic processes was unaltered after ERT administration (data not shown).
A search for proteins that interact with collagen in the extracellular matrix identified 5 proteins; P02671, P02751 (fibronectin 1), P02675, P51884 (lumican), and P02679.
We detected a relatively high number of proteins involved in the immune response. The results of the functional analysis of these proteins are shown in Figure 5 and Figure 6. Figure 5a shows the results of functional enrichment analysis of proteins implicated in inflammation, including proteins related to interleukins (ILs) IL-1, IL-1β, IL-6, IL-8, and to interferon γ (TNF-γ). Proteins were grouped according to their participation in different processes, including the cellular response to interferon γ, the cellular response to IL-1, the cellular response to IL-6, positive regulation of IL-1 secretion, positive regulation of IL-1β secretion, negative regulation of IL-1β secretion, negative regulation of IL-6 production, and positive regulation of IL-8 biosynthesis. In all cases, we observed a reduction in the expression of the proteins of interest24 h after ERT administration (ERT-b), except for proteins involved in positive regulation of IL-1 secretion (no change observed after ERT). In UG we detected 1.79% of proteins related to IL-1, IL-1β, and TNF-γ. This corresponding value in the ERT-a group (1.25%) was closer to that observed in controls (1.23%). In the case of proteins related to IL-6 and IL-8 we detected values of 0.88% in UG,0.63% in ERT-a, and 0.61 in the control group.
The results of the functional enrichment analysis of proteins involved in biological processes related tothe immune response, immunoglobulin production, positive regulation of B cell activation, negative regulation inflammatory response, and leukocyte migrationare shown in Figure 5b. In all cases, ERT increased the expression of proteins involved in each of these biological processes, except for negative regulation of the inflammatory response.
Functional enrichment analysis was also performed to identify protein partner expression between groups of proteins related to inflammatory processes (biological processes relating to (i) the complement system, including the classical complement activation pathway, (ii) general complement activation, and (iii) the alternative complement activation pathway). The results (Figure 6) reveal a positive effect of ERT, although levels of protein expression levels did not reach those observed in the control group (Figure 6a). Figure 6b shows the results of the functional enrichment analysis of proteins involved in the oxidative stress response. Treatment (ERT-a) resulted in protein levels similar to those observed in the control group for proteins involved in both complement activation and oxidative stress.

2.2. SWATH Quantitative Analysisof Proteins in Plasma Samples

Quantitative analysis of proteins using SWATH revealed significant alterations in the expression of multiple proteins in the ERT-a and ERT-b groups relative to the control group (Table 1).

Quantitative Analysis Using SWATH

Proteins that were upregulated or down regulated in each group may be relevant in MPS IVA disease. Differences in protein expression in the UG may reveal candidate biomarkers of disease. Proteins upregulated or downregulated in the ERT-a and ERT-b groups can be considered as markers of disease evolution as they were altered before treatment. We identified certain proteins that were downregulated by more than 100 fold, including AMBP, ASM3B, KV117, and ZBTB4 (highlighted in dark red in Figure 7 and Table 2).
Of the proteins that were the upregulated and downregulated in the UG, ERT-a, and ERT-b groups relative to controls, 31% were related to immune system function. It is important to point out that most of the proteins for which changes in expression were observed in the study groups were down regulated. In addition to upregulated protein in the UG, ERT-a, and ERT-b groups relative to CG, we alsodetected2 immunoglobulins (KV320 and LV321) that were not present in the CG.
Some of these downregulated proteins (highlighted in red in Table 2 and in the heatmap in Figure 8) may be of particular importance in MPS IVA.
In order to shown the results in a better way, we represent a heatmap using the proteins up and down regulated shown in Table 2. Moreover, the proteins were classified into different functions and were highline in different colors. Among these proteins, those shown in Figure 9 and Figure 10 stand out.
The SWATH analysis identified proteins that were downregulated in the UG and the ERT-a and b groups relative to controls, or were absent entirely. These included FA10 and SAA4, which were absent in all patient groups, and RET4 and C1QB, which were absent in UG and in UG and ERT-a, respectively (see Figure 10 and Figure 11).
In the UG group, we detected changes in expression in 4 proteins with respect to the 3 other groups (Figure 12 and Table 3).
Because these proteins were likely dysregulated as a consequence of the disease, they may include potential disease biomarkers that are either up-or down regulated in MPS IVA.
Theproteinsfor which changes in expression were detected in the UG included severalimmunoglobulins. This is unsurprising, given the important inflammatory component of MPS. Moreover, proteins that were dysregulated in this group included some previously proposed as biomarkers of MPS (e.g., A1AT[alpha 1 antitrypsin]) [33].We also observed dysregulation of proteins involved in inflammation (S10A9), binding to collagen (C1QC), and prothrombin (THRB), a thyroid hormone receptor that mediates the biological activities of thyroid hormone [34]. Proteins that were downregulated in the UG included C1QB (in line with previously reports, Figure 9) and the serine protease inhibitor A2AP (Figure 10 and Table 4).
The most important dysregulated protein identified in this analysis was A1AT. Analysis of A1AT levels in individuals in each of the different groups (Figure 13) revealed significant upregulation in the UG compared with each of the other groups. A1AT levels in the ERT-a and ERT-b groups were comparable to those observed in controls. This protein may constitute not only a biomarker of MPSIVA, but also an indicator of disease progression, potentially providing a better means of monitoring treatment responses in MPS IVA.
While treatment resulted in changes in the expression of proteins in both the ERT-a and ERT-b groups (Table 4), none of these proteins were common to both groups.
In this analysis, we identified proteins that were expressed in only one of the treated groups (ERT-a or ERT-b) and were absent from all other groups. SAA1 and APOF weredetected exclusively in the ERT-a group and TPIS, QSOX1, and KV108exclusively in the ERT-b group. Changes in the expression of these proteins can be considered a consequence of ERT. More detailed analyses will be required to determine their role in the treatment response. GELS was upregulated in the ERT-b group compared with controls (Figure 11). The fact that this protein was upregulated in the untreated group, and even more so in the ERT-b group, suggests that its expression is increased as a consequence of the disease process.

3. Discussion

Quality surrogate biomarkers to monitor disease progression and treatment response are currently lacking in MPS in general, and in Morquio A disease in particular. Currently, the diagnosis of MPS IVA is established by GAG assay, enzymatic activity analysis, and genetic studies. To date, the most useful biomarker for MPS IVA is GAG levels in urine and blood. High levels of KS and C6S are indicative of MPS IVA disease. However, decreases in urinary KS levels after ERT do not imply that treatment is working as urinary GAG originates from the kidneys and does not correlate with bone lesion severity in MPS IVA [4]. Blood KS could constitute a more appropriate readout of therapeutic effect, as it is directly derived from bone. Another problem is that GAG levels in MPS IVA patients decrease naturally with age and therefore may be normalized or nearly normalized without treatment due to destruction of the cartilage in which KS is predominantly synthesized [35]. While a recent study reported that the disaccharide chondroitin sulfate can be used to monitor treatment with ERT, to date this biomarker has only been studied in a small number of samples from MPS IVA patients. [36].
Proteomic studies can provide greater insight into the pathological mechanisms underlying diseases and help identify candidate biomarkers. Sample collection is complicated in MPS IVA, in which bone is the main tissue affected. Proteomic analyses of fibroblasts and leukocytes from MPS IVA patients published by our group [37] have demonstrated correction of the levels of certain dysregulated proteins following ERT. In the present study, we focused our analysis on the effects of MPS IVA and ERT on plasma proteins implicated in bone physiology, inflammatory processes, and oxidative stress. To this end, we performed both qualitative and quantitative analyses. Qualitative analyses revealed changes in the expression of proteins associated with coagulation, bone regeneration, collagen reorganization, inflammation, and oxidative stress. ERT reversed the dysregulation of plasma proteins involved in hemostasis and oxidative stress, but had no such effect on proteins implicated in inflammation or the immune response. The inflammatory component of MPS IVA results from GAG accumulation within tissues that activate Toll-like receptors and inflammation pathways [38]. Anti-inflammatory drugs can inhibit the activity of altered cytokines and interleukins in MPS IVA [4]. In our study population, expression levels of certain dysregulated ILs (e.g., IL1 and IL1β) were not corrected after ERT. ILs have been described by other authors as possible candidate biomarkers for MPS IVA [23]. The combination of anti-inflammatory drugs (infliximab) and ERT has been tested in rats with MPS IVA, in which analysis of cultured articular chondrocytes showed that this therapy restored collagen IIA1 expression and reduced the expression of apoptotic markers [39]. However, this treatment only suppressed hyperplasia of synovial cells. Micro-computed tomography (CT) analysis revealed no significant positive effects on bone microarchitecture and no histological improvement in bone growth plates was observed. In our study, ERT resulted in only partial correction of the expression of dysregulated proteins implicated in the immune response. It is not clear if the possible alterations in the immune response that can affect patients with MPS IVA are due to the neutralizing antibodies that they produce against the ERT enzyme and there are no studies to clarify this.
Our quantitative analysis examined potential alterations in the expression of 215 proteins included in the SWATH library. Using this method, we detected proteins that showed changes in expression all other groups relative to healthy controls. In line with the findings of the qualitative analysis, several of the dysregulated proteins identified using this approach are involved in bone metabolism. Two such proteins, VTNC and FINC, participate in remodeling of the extracellular matrix (ECM) and in osseous integration mechanisms [40]. VTNC is a glycoprotein that circulates in blood and is also found in the ECM. This protein can bind GAGs, collagen, plasminogen, and the urokinase receptor. VTNC regulates the proteolytic degradation of this matrix. Moreover, VTNC binds to the complement system, heparin, and thrombin-antithrombin III complexes, facilitating the regulation of the immune response and clot formation [41].
Fetuin A (FETUA) [42] is a family member structurally related to the cystatin-like protein domain. Cystatin can inhibit papain, caspases, and cathepsins. FETUA inhibits insulin receptor autophosphorylation and plays a key role in insulin-dependent metabolism. Fetuin-A is also involved in osteogenesis and bone resorption and in the control of calcium salt precipitation in blood [43,44]. This protein accumulates in the bone matrix and inhibits transforming growth factor, necessary for bone mineralization [39]. Fetuin-A is also an anti-inflammatory mediator that participates in macrophage deactivation. Decreases in plasma levels of this protein may occur in response to inflammatory processes, such as those that occur in MPS IVA [45]. A proteomic analysis of pig bone samples from Procopio et al. presents many of the proteins located in different regions of the bone, such as in the proximal area of the tibia where the FA10, VTNC, SAMP, CRP, and FETUA proteins were found [46]. The remaining protein expression found in the study was located in the metatarsus, the midshaft of the tibia and the proximal femur. The authors concluded that three proteins, A1T1, chromogranin A, and FETUA, could be considered biological age markers in the tibia [46]. To identify candidate protein biomarkers of MPS IVA, Martell et al. performed quantitative multiplexed assays using plasma samples from 78 MPS IVA patients and 58 healthy controls [33]. The authors identified three candidate proteins: alpha-1-antitrypsin, lipoprotein(a), and serum amyloid. In our analysis, we found that A1T1 expression was increased in untreated MPS IVA patients and decreased in ERT-treated patients.
Another protein identified in our study was serum amyloid A 4 (SAA4). The SAA protein family has been described as important clinical markers of inflammation and acts as a precursor of elements that are involved in secondary reactive amyloidosis [47]. Moreover, SAA family proteins are known to contribute to cellular cholesterol homeostasis, modulate intracellular calcium levels, and promote signaling cascades [48,49,50]. SAA proteins also participate in lipid transport and metabolism and are implicated in atherosclerosis. SAA proteins share several common functions with cytokines, including cell–cell communication and feedback in inflammatory, immunological, neoplastic, and protective pathways [51]. Retinol-binding protein 4 (RET4) is the main transporter of retinol in the bloodstream, whose metabolites may contribute to osteogenesis, chondrogenesis, and adipogenesis. RET4 is mainly produced in the liver (but is also expressed in several extrahepatic tissues, including the limbs), where it mobilizes hepatic retinol stores [52].
Another protein of interest detected in the present study is CLUS. This protein helps to prevent cell stress-induced apoptosis and participates in the aggregation of blood plasma proteins [53,54,55,56,57], along with other proteins such as APP and APOC2. CLUS can also inhibit the formation of amyloid fibrils.

Study Limitations

The SWATH method is a specific, reproducible, and sensitive approach that allows relative or absolute quantification of proteins. However, some limitations of this approach should be noted. First, the number of proteins identified/quantified is primarily limited by the composition of the spectral library. In this study, the library was generated from a pool of proteins for each group of study, and consisted of 215 proteins, each with a FDR < 1%. In contrast to the selected reaction monitoring (SRM) technique, in which only three transitions are quantified, the SWATH technique quantifies seven transitions are for each peptide. Therefore, 10 peptides per protein and seven fragments per peptide must be identified in order to extract the peak areas necessary for quantification. Despite increasing the precision of the technique, this aspect also constitutes a limitation, as proteins in the library for which less than 10 peptides are identified will not be quantified.
Another limitation of our study is that normalization of the expression of proteins in plasma does not necessarily correlate with normalization in chondrocytes in the vascular region. Enzyme is easily taken up by hepatic proteins, restoring dysregulated metabolic pathways in these cells while the infused enzyme circulates in the blood. It is critical to understand whether the blood/plasma proteins are dysregulated in bone and cartilage in MPS IVA patients in order to identify potential diagnostic biomarkers of disease severity. Moreover, when a bone-penetrating drug becomes available, surrogate disease, progression disease, or treatment outcome biomarker will be essential.
The selected candidate proteins must be studied using techniques that are used daily in clinical practice to demonstrate their validity. However, the low number of patients with this disease makes this analysis difficult.

4. Materials and Methods

4.1. Study Workflow

The graphical abstract depicts the workflow used for the qualitative (DDA-LC-MS/MS) and quantitative (SWATH-MS) analysis of proteins in plasma samples from healthy controls and MPSIVA patients (treated or untreated) (Figure 14).

4.2. Samples

After obtaining written informed consent, blood samples from MPS IVA patients were acquired for proteomic analysis of plasma from three Spanish centers of reference for MPS IVA: eight from untreated MPS IVA patients; and five from MPS IVA patients who had been undergoing ERT for a mean duration of 7 years, and from whom samples were taken 24 h before (ERT-a group) and 24 h after (ERT-b group) one of their weekly ERT infusions. All patients presented the classical phenotype as is shown in the previous study [32]. Another 6 blood samples were acquired from 6 healthy donors (Figure 15). This study was approved (22 May 2017) by the local Ethics Committee of Santiago-Lugo (Reference number 2017/298).

4.3. Protein Extraction

Plasma samples were depleted in order to eliminate the most abundant proteins prior to the proteomic analysis. Leukocytes were sonicated to rupture the cell membrane and centrifuged for 10 min at 10,000 rpm and 4 °C. Protein extracts from the supernatant were recovered and subsequently frozen at −20 °C.

4.4. Enzyme Activity Test

An enzyme activity test was used to analyze the GALNS activity in plasma samples [58]. Enzymatic activity in samples from ERT-a patients was lower than reference values (1.8–20 nmol/h/mg protein) established by the laboratory to diagnose MPSIVA.

4.5. Proteomic Analysis

Proteomic analysis was performed as described previously by our group [32]. Plasma samples were first depleted as previously described [59]. A restrictive analysis was performed to select only proteins with a p-value < 0.05 and a fold change >1.5 as possible candidate biomarkers. A heat map was generated using http://www.heatmapper.ca/expression/, accessed on 1 May 2021.

5. Conclusions

Morquio A is a rare disease. Diagnosis can be complicated in the absence of clinical signs and symptoms, and bone impact therapy remains undeveloped. The discovery of effective biomarkers for rare diseases is essential to improve diagnosis and treatment. Based on our findings, we propose plasma levels of FETUA, VINT, A1AT, and CLUS as potential candidate biomarkers for MPS IVA. These proteins may influence bone and cartilage metabolism but more studies are needed. Other proteins identified in this study appear to provide an indication of the effects produced as a result of bone alteration. In particular, further investigation of the role of A1AT in bone is warranted, as several studies have reported dysregulation of this protein in MPS IVA.

Supplementary Materials

Author Contributions

Contributions were provided by the authors in the statements: conceptualization, J.V.Á., S.B.B., F.J.O.-E. and M.L.C.; methodology, S.B.B., J.V.Á., M.P.C.-V., S.B.B. and C.C.; validation, J.V.Á. and S.B.B.; sample collection. O.L.-S., and M.L.C.; formal analysis, J.V.Á., S.B.B., S.B.-G. and C.C.; writing—original draft preparation, J.V.Á., S.B.B.; writing—review and editing, M.L.C., F.J.O.-E. and S.T.; visualization, S.B.B.; supervision, M.L.C., F.J.O.-E. and S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the local Ethics Committee of Santiago-Lugo (data approval 22 May 2017 Reference number 2017/298).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available in this article and this Supplementary Material.

Acknowledgments

Thanks to all participating patients.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Matalon, R.; Arbogast, B.; Justice, P.; Brandt, I.K.; Dorfman, A. Morquio’s syndrome: Deficiency of a chondroitin sulfate N-acetylhexosamine sulfate sulfatase. Biochem. Biophys Res. Commun. 1974, 61, 759–765. [Google Scholar] [CrossRef]
  2. Singh, J.; Di Ferrante, N.; Niebes, P.; Tavella, D. N-acetylgalactosamine-6-sulfate sulfatase in man. Absence of the enzyme in Morquio disease. J. Clin. Investig. 1976, 57, 1036–1040. [Google Scholar] [CrossRef] [PubMed]
  3. Di Ferrante, N.; Ginsberg, L.C.; Donnelly, P.V.; Di Ferrante, D.T.; Caskey, C.T. Deficiencies of glucosamine-6-sulfate or galactosamine-6-sulfate sulfatases are responsible for different mucopolysaccharidoses. Science 1978, 199, 79–81. [Google Scholar] [CrossRef] [PubMed]
  4. Tomatsu, S.; Yasuda, E.; Patel, P.; Ruhnke, K.; Shimada, T.; Mackenzie, W.G.; Mason, R.; Thacker, M.M.; Theroux, M.; Montaño, A.M.; et al. Morquio A syndrome: Diagnosis and current and future therapies. Pediatr. Endocrinol. Rev. PER 2014, 12 (Suppl. 1), 141–151. [Google Scholar] [PubMed]
  5. Hendriksz, C.J.; Berger, K.I.; Giugliani, R.; Harmatz, P.; Kampmann, C.; Mackenzie, W.G.; Raiman, J.; Villarreal, M.S.; Savarirayan, R. International Guidelines for the Management and Treatment of Morquio A Syndrome. Am. J. Med. Genet. Part A 2015, 167, 11–25. [Google Scholar] [CrossRef] [PubMed]
  6. Khan, S.; Alméciga-Díaz, C.J.; Sawamoto, K.; Mackenzie, W.G.; Theroux, M.C.; Pizarro, C.; Mason, R.W.; Orii, T.; Tomatsu, S. Mucopolysaccharidosis IVA and glycosaminoglycans. Mol. Genet. Metab. 2017, 120, 78–95. [Google Scholar] [CrossRef] [Green Version]
  7. Peracha, H.; Sawamoto, K.; Averill, L.; Kecskemethy, H.; Theroux, M.; Thacker, M.; Nagao, K.; Pizarro, C.; Mackenzie, W.; Kobayashi, H.; et al. Diagnosis and prognosis of Mucopolysaccharidosis IVA. Mol. Genet. Metab. 2018, 125, 18–37. [Google Scholar] [CrossRef] [PubMed]
  8. Melbouci, M.; Mason, R.W.; Suzuki, Y.; Fukao, T.; Orii, T.; Tomatsu, S. Growth impairment in mucopolysaccharidoses. Mol. Genet. Metab. 2018, 124, 1–10. [Google Scholar] [CrossRef]
  9. Tomatsu, S.; Montaño, A.M.; Oikawa, H.; Smith, M.; Barrera, L.; Chinen, Y.; Thacker, M.M.; Mackenzie, W.G.; Suzuki, Y.; Orii, T. Mucopolysaccharidosis type IVA (Morquio A disease): Clinical review and current treatment. Curr. Pharm. Biotechnol. 2011, 12, 931–945. [Google Scholar] [CrossRef] [PubMed]
  10. Lavery, C.; Hendriksz, C. Mortality in patients with Morquio syndrome A. JIMD Rep. 2015, 15, 59. [Google Scholar]
  11. Hendriksz, C.; Santra, S.; Jones, S.A.; Geberhiwot, T.; Jesaitis, L.; Long, B.; Qi, Y.; Hawley, S.M.; Decker, C. Safety, immunogenicity, and clinical outcomes in patients with Morquio A syndrome participating in 2 sequential open-label studies of elosulfasealfa enzyme replacement therapy (MOR-002/MOR-100), representing 5 years of treatment. Mol. Genet. Metab. 2018, 123, 479–487. [Google Scholar] [CrossRef] [PubMed]
  12. Hendriksz, C.J. Elosulfasealfa (BMN 110) for the treatment of mucopolysaccharidosis IVA (Morquio A Syndrome). Expert Rev. Clin. Pharmacol. 2016, 9, 1521–1532. [Google Scholar] [CrossRef] [Green Version]
  13. Akyol, M.U.; Alden, T.D.; Amartino, H.; Ashworth, J.; Belani, K.; Berger, K.I.; Borgo, A.; Braunlin, E.; Eto, Y.; Gold, J.I.; et al. Recommendations for the management of MPS IVA: Systematic evidence- and consensus-based guidance. Orphanet J. Rare Dis. 2019, 14, 137. [Google Scholar] [CrossRef] [Green Version]
  14. Tomatsu, S.; Sawamoto, K.; Shimada, T.; Bober, M.B.; Kubaski, F.; Yasuda, E.; Mason, R.W.; Khan, S.; Alméciga-Díaz, C.J.; Barrera, L.A.; et al. Enzyme replacement therapy for treating mucopolysaccharidosis type IVA (Morquio A syndrome): Effect and limitations. Expert Opin. Orphan Drugs 2015, 3, 1279–1290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tomatsu, S.; Montaño, A.M.; Dung, V.C.; Ohashi, A.; Oikawa, H.; Oguma, T.; Orii, T.; Barrera, L.; Sly, W.S. Enhancement of drug delivery: Enzyme-replacement therapy for murine Morquio A syndrome. Mol. Ther. 2010, 18, 1094–1102. [Google Scholar] [CrossRef] [PubMed]
  16. Sawamoto, K.; Stapleton, M.; Alméciga-Díaz, C.J.; Espejo-Mojica, A.J.; Losada, J.C.; Suarez, D.A.; Tomatsu, S. Therapeutic Options for Mucopolysaccharidoses: Current and Emerging Treatments. Drugs 2019, 79, 1103–1134. [Google Scholar] [CrossRef]
  17. Long, B.; Tompkins, T.; Decker, C.; Jesaitis, L.; Khan, S.; Slasor, P.; Harmatz, P.; O’Neill, C.A.; Schweighardt, B. Long-term Immunogenicity of Elosulfase Alfa in the Treatment of Morquio A Syndrome: Results From MOR-005, a Phase III Extension Study. Clin. Ther. 2017, 39, 118–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Doherty, C.; Stapleton, M.; Piechnik, M.; Mason, R.W.; Mackenzie, W.G.; Yamaguchi, S.; Kobayashi, H.; Suzuki, Y.; Tomatsu, S. Effect of enzyme replacement therapy on the growth of patients with Morquio A. J. Hum. Genet. 2019, 64, 625–635. [Google Scholar] [CrossRef] [PubMed]
  19. Tomatsu, S.; Averill, L.W.; Sawamoto, K.; Mackenzie, W.G.; Bober, M.B.; Pizarro, C.; Goff, C.J.; Xie, L.; Orii, T.; Theroux, M. Obstructive airway in Morquio A syndrome, the past, the present and the future. Mol. Genet. Metab. 2016, 117, 150–156. [Google Scholar] [CrossRef] [Green Version]
  20. Tomatsu, S.; Sawamoto, K.; Alméciga-Díaz, C.J.; Shimada, T.; Bober, M.B.; Chinen, Y.; Yabe, H.; Montaño, A.M.; Giugliani, R.; Kubaski, F.; et al. Impact of enzyme replacement therapy and hematopoietic stem cell transplantation in patients with Morquio A syndrome. Drug Des. Dev. Ther. 2015, 9, 1937–1953. [Google Scholar] [CrossRef] [Green Version]
  21. Qi, Y.; Musson, D.G.; Schweighardt, B.; Tompkins, T.; Jesaitis, L.; Shaywitz, A.J.; Yang, K.; O’Neill, C.A. Pharmacokinetic and Pharmacodynamic Evaluation of Elosulfase Alfa, an Enzyme Replacement Therapy in Patients with Morquio A Syndrome. Clin. Pharmacokinet. 2014, 53, 1137–1147. [Google Scholar] [CrossRef] [Green Version]
  22. De Franceschi, L.; Roseti, L.; Desando, G.; Facchini, A.; Grigolo, B. A molecular and histological characterization of cartilage from patients with Morquio syndrome. Osteoarthr. Cartil. 2007, 15, 1311–1317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Fujitsuka, H.; Sawamoto, K.; Peracha, H.; Mason, R.W.; Mackenzie, W.; Kobayashi, H.; Yamaguchi, S.; Suzuki, Y.; Orii, K.; Orii, T.; et al. Biomarkers in patients with mucopolysaccharidosis type II and IV. Mol. Genet. Metab. Rep. 2019, 19, 100455. [Google Scholar] [CrossRef] [PubMed]
  24. Tan, H.T.; Chung, M.C.M. Label-Free Quantitative Phosphoproteomics Reveals Regulation of Vasodilator-Stimulated Phosphoprotein upon Stathmin-1 Silencing in a Pair of Isogenic Colorectal Cancer Cell Lines. Proteomics 2018, 18, e1700242. [Google Scholar] [PubMed]
  25. Braccia, C.; Espinal, M.P.; Pini, M.; De Pietri Tonelli, D.; Armirotti, A. A new SWATH ion library for mouse adult hippocampal neural stem cells. Data Brief 2018, 18, 1–8. [Google Scholar] [CrossRef] [PubMed]
  26. Jayabalan, N.; Lai, A.; Nair, S.; Guanzon, D.; Scholz-Romero, K.; Palma, C.; McIntyre, H.D.; Lappas, M.; Salomon, C. Quantitative proteomics by SWATH-MS suggest an association between circulating exosomes and maternal metabolic changes in gestational diabetes mellitus. Proteomics 2019, 19, e1800164. [Google Scholar] [CrossRef] [Green Version]
  27. Jamwal, R.; Barlock, B.J.; Adusumalli, S.; Ogasawara, K.; Simons, B.L.; Akhlaghi, F. Multiplex and Label-Free Relative Quantification Approach for Studying Protein Abundance of Drug Metabolizing Enzymes in Human Liver Microsomes using SWATH-MS. J. Proteome Res. 2017, 16, 4134–4143. [Google Scholar] [CrossRef] [PubMed]
  28. Zhu, Y.; Zhu, J.; Lu, C.; Zhang, Q.; Xie, W.; Sun, P.; Dong, X.; Yue, L.; Sun, Y.; Yi, X.; et al. Identification of Protein Abundance Changes in Hepatocellular Carcinoma Tissues Using PCT–SWATH. Proteom. Clin. Appl. 2018, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Li, H.; Mao, Y.; Xiong, Y.; Zhao, H.H.; Shen, F.; Gao, X.; Yang, P.; Liu, X.; Fu, D. A Comprehensive Proteome Analysis of Peripheral Blood Mononuclear Cells (PBMCs) to Identify Candidate Biomarkers of Pancreatic Cancer. Cancer Genom. Proteom. 2019, 16, 81–89. [Google Scholar] [CrossRef] [Green Version]
  30. Chang, R.Y.; Etheridge, N.; Nouwens, A.S.; Dodd, P.R. SWATH analysis of the synaptic proteome in Alzheimer’s disease. Neurochem. Int. 2015, 87, 1–12. [Google Scholar] [CrossRef] [PubMed]
  31. Gregorich, Z.R.; Ge, Y. Top-down proteomics in health and disease: Challenges and opportunities. Proteomics 2014, 14, 1195–1210. [Google Scholar] [CrossRef] [PubMed]
  32. Álvarez, J.V.; Bravo, S.B.; Chantada-Vazquez, M.P.; Colón, C.; De Castro, M.J.; Morales, M.; Vitoria, I.; Tomatsu, S.; Otero-Espinar, F.J.; Couce, M.L. Characterization of New Proteomic Biomarker Candidates in Mucopolysaccharidosis Type IVA. Int. J. Mol. Sci. 2020, 22, 226. [Google Scholar] [CrossRef]
  33. Martell, L.; Lau, K.; Mei, M.; Burnett, V.; Decker, C.; Foehr, E.D. Biomarker analysis of Morquio syndrome: Identification of disease state and drug responsive markers. Orphanet J. Rare Dis. 2011, 6, 84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zhang, A.; Sieglaff, D.H.; York, J.P.; Suh, J.H.; Ayers, S.D.; Winnier, G.E.; Kharitonenkov, A.; Pin, C.; Zhang, P.; Webb, P.; et al. Thyroid hormone receptor regulates most genes independently of fibroblast growth factor 21 in liver. J. Endocrinol. 2015, 224, 289–301. [Google Scholar] [CrossRef] [Green Version]
  35. Shimada, T.; Tomatsu, S.; Mason, R.W.; Yasuda, E.; Mackenzie, W.G.; Hossain, J.; Shibata, Y.; Montaño, A.M.; Kubaski, F.; Giugliani, R.; et al. Di-sulfated Keratan Sulfate as a Novel Biomarker for Mucopolysaccharidosis II, IVA, and IVB. JIMD Rep. 2015, 21, 1–13. [Google Scholar]
  36. Chin, S.J.; Saville, J.T.; McDermott, B.K.; Zankl, A.; Fletcher, J.M.; Fuller, M. Chondroitin sulfate disaccharide is a specific and sensitive biomarker for mucopolysaccharidosis type IVA. JIMD Rep. 2020, 55, 68–74. [Google Scholar] [CrossRef]
  37. Álvarez, J.V.; Bravo, S.B.; García-Vence, M.; De Castro, M.J.; Luzardo, A.; Colón, C.; Tomatsu, S.; Otero-Espinar, F.J.; Couce, M.L. Proteomic Analysis in Morquio A Cells Treated with Immobilized Enzymatic Replacement Therapy on Nanostructured Lipid Systems. Int. J. Mol. Sci. 2019, 20, 4610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Algahim, M.F.; Almassi, G.H. Current and emerging management options for patients with Morquio A syndrome. Ther. Clin. Risk Manag. 2013, 9, 45–53. [Google Scholar]
  39. Eliyahu, E.; Wolfson, T.; Ge, Y.; Jepsen, K.J.; Schuchman, E.H.; Simonaro, C.M. Anti-TNF-Alpha Therapy Enhances the Effects of Enzyme Replacement Therapy in Rats with Mucopolysaccharidosis Type VI. PLoS ONE 2011, 6, e22447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Barallobre-Barreiro, J.; Oklu, R.; Lynch, M.; Fava, M.; Baig, F.; Yin, X.; Barwari, T.; Potier, D.N.; Albadawi, H.; Jahangiri, M.; et al. Extracellular matrix remodelling in response to venous hypertension: Proteomics of human varicose veins. Cardiovasc. Res. 2016, 110, 419–430. [Google Scholar] [CrossRef]
  41. De Boer, H.C.; Preissner, K.T.; Bouma, B.N.; de Groot, P.G. Binding of vitronectin-thrombin-antithrombin III complex to human endothelial cells is mediated by the heparin binding site of vitronectin. J. Biol. Chem. 1992, 267, 2264–2268. [Google Scholar] [CrossRef]
  42. Gangneux, C.; Daveau, M.; Hiron, M.; Derambure, C.; Papaconstantinou, J.; Salier, J.P. The inflammation-induced down-regulation of plasma Fetuin-A (alpha2HS-Glycoprotein) in liver results from the loss of interaction between long C/EBP isoforms at two neighbouring binding sites. Nucleic Acids Res. 2003, 31, 5957–5970. [Google Scholar] [CrossRef] [Green Version]
  43. Brylka, L.; Jahnen-Dechent, W. The role of fetuin-A in physiological and pathological mineralization. Calcif. Tissue Int. 2013, 93, 355–364. [Google Scholar] [CrossRef] [PubMed]
  44. Jahnen-Dechent, W.; Heiss, A.; Schäfer, C.; Ketteler, M. Fetuin-A regulation of calcified matrix metabolism. Circ. Res. 2011, 108, 1494–1509. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, H.; Sama, A.E. Anti-inflammatory role of Fetuin-A in Injury and Infection. Curr. Mol. Med. 2012, 12, 625. [Google Scholar] [CrossRef] [Green Version]
  46. Procopio, N.; Chamberlain, A.T.; Buckley, M. Intra- and Interskeletal Proteome Variations in Fresh and Buried Bones. J. Proteome Res. 2017, 16, 2016–2029. [Google Scholar] [CrossRef]
  47. Liao, F.; Lusis, A.J.; Berliner, J.A.; Fogelman, A.M.; Kindy, M.; de Beer, M.C.; de Beer, F.C. Serum amyloid A protein family. Differential induction by oxidized lipids in mouse strains. Arterioscler. Thromb. J. Vasc. Biol. 1994, 14, 1475–1479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Badolato, R.; Johnston, J.A.; Wang, J.M.; McVicar, D.; Xu, L.L.; Oppenheim, J.J.; Kelvin, D.J. Serum amyloid A induces calcium mobilization and chemotaxis of human monocytes by activating a pertussis toxin-sensitive signaling pathway. J. Immunol. 1995, 155, 4004–4010. [Google Scholar] [PubMed]
  49. Artl, A.; Marsche, G.; Lestavel, S.; Sattler, W.; Malle, E. Role of serum amyloid A during metabolism of acute-phase HDL by macrophages. Arterioscler. Thromb. Vasc. Biol. 2000, 20, 763–772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Baranova, I.N.; Vishnyakova, T.G.; Bocharov, A.V.; Kurlander, R.; Chen, Z.; Kimelman, M.L.; Remaley, A.T.; Csako, G.; Thomas, F.; Eggerman, T.L. Serum amyloid A binding to CLA-1 (CD36 and LIMPII analogous-1) mediates serum amyloid A protein-induced activation of ERK1/2 and p38 mitogen-activated protein kinases. J. Biol. Chem. 2005, 280, 8031–8040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Sack, G.H. Serum amyloid—A review. Mol. Med. 2018, 24, 46. [Google Scholar] [CrossRef]
  52. Hatfield, J.T.; Anderson, P.J.; Powell, B.C. Retinol-binding protein 4 is expressed in chondrocytes of developing mouse long bones: Implications for a local role in formation of the secondary ossification center. Histochem. Cell Biol. 2013, 139, 727–734. [Google Scholar] [CrossRef] [PubMed]
  53. Abdallah, B.M.; Alzahrani, A.M.; Kassem, M. Secreted Clusterin protein inhibits osteoblast differentiation of bone marrow mesenchymal stem cells by suppressing ERK1/2 signaling pathway. Bone 2018, 110, 221–229. [Google Scholar] [CrossRef] [PubMed]
  54. Poon, S.; Easterbrook-Smith, S.B.; Rybchyn, M.S.; Carver, J.A.; Wilson, M.R. Clusterin is an ATP-independent chaperone with very broad substrate specificity that stabilizes stressed proteins in a folding-competent state. Biochemistry 2000, 39, 15953–15960. [Google Scholar] [CrossRef] [PubMed]
  55. Poon, S.; Rybchyn, M.S.; Easterbrook-Smith, S.B.; Carver, J.A.; Pankhurst, G.J.; Wilson, M.R. Mildly Acidic pH Activates the Extracellular Molecular Chaperone Clusterin. J. Biol. Chem. 2002, 277, 39532–39540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Stewart, E.M.; Aquilina, J.A.; Easterbrook-Smith, S.B.; Murphy-Durland, D.; Jacobsen, C.; Moestrup, S.; Wilson, M.R. Effects of Glycosylation on the Structure and Function of the Extracellular Chaperone Clusterin. Biochemistry 2007, 46, 1412–1422. [Google Scholar] [CrossRef]
  57. Wyatt, A.R.; Wilson, M.R. Identification of Human Plasma Proteins as Major Clients for the Extracellular Chaperone Clusterin. J. Biol. Chem. 2010, 285, 3532–3539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Camelier, M.V.; Burin, M.G.; De Mari, J.; Vieira, T.A.; Marasca, G.; Giugliani, R. Practical and reliable enzyme test for the detection of mucopolysaccharidosis IVA (Morquio Syndrome type A) in dried blood samples. Clin. Chim. Acta 2011, 412, 1805–1808. [Google Scholar] [CrossRef] [Green Version]
  59. Del PilarChantada-Vázquez, M.; López, A.C.; Vence, M.G.; Vázquez-Estévez, S.; Acea-Nebril, B.; Calatayud, D.G.; Jardiel, T.; Bravo, S.B.; Núñez, C. Proteomic investigation on bio-corona of Au, Ag and Fe nanoparticles for the discovery of triple negative breast cancer serum protein biomarkers. J. Proteom. 2020, 212, 103581. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study groups: controls and MPS IVA patients.CG, control group; MPS IVA, mucopolysaccharidosis type IVA; ERT-a, MPS IVA patients sampled before ERT; ERT-b, MPS IVA patients sampled 24 h after ERT.
Figure 1. Study groups: controls and MPS IVA patients.CG, control group; MPS IVA, mucopolysaccharidosis type IVA; ERT-a, MPS IVA patients sampled before ERT; ERT-b, MPS IVA patients sampled 24 h after ERT.
Ijms 22 06165 g001
Figure 2. Venn diagram showing the distribution of proteins identified in the qualitative proteomic analysis across each of the 4 study groups.
Figure 2. Venn diagram showing the distribution of proteins identified in the qualitative proteomic analysis across each of the 4 study groups.
Ijms 22 06165 g002
Figure 3. Functional enrichment analysis using FunRich. Proteins identified in each the different groupswere submitted to molecular functional analysis according to their involvement inhemostasis (a) and platelet function (b). CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group. Proteins listed for qualitative proteomic analysis appear in Supplementary Table S2.
Figure 3. Functional enrichment analysis using FunRich. Proteins identified in each the different groupswere submitted to molecular functional analysis according to their involvement inhemostasis (a) and platelet function (b). CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group. Proteins listed for qualitative proteomic analysis appear in Supplementary Table S2.
Ijms 22 06165 g003
Figure 4. Functional enrichment analysis using FunRich. Proteins identified in each the different groups were submitted to molecular functional analysis according to their involvement in: (a) processes that occur at the bone level; (b) collagen function (activity/processing/biosynthesis); (c) collagen-containing extracellular matrix. CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group. Proteins included in this analysis, identified in the qualitative proteomic analysis, are listed in the Supplementary Table S2.
Figure 4. Functional enrichment analysis using FunRich. Proteins identified in each the different groups were submitted to molecular functional analysis according to their involvement in: (a) processes that occur at the bone level; (b) collagen function (activity/processing/biosynthesis); (c) collagen-containing extracellular matrix. CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group. Proteins included in this analysis, identified in the qualitative proteomic analysis, are listed in the Supplementary Table S2.
Ijms 22 06165 g004
Figure 5. Functional enrichment analysis using FunRich. Proteins identified in each the different groups were submitted to molecular functional analysis according to their involvement in the immune response: proteinsrelated tointerleukins (ILs) IL-1, IL-1β, IL-6, IL-8, and interferon γ (a); proteins involved in immunoglobulin production, positive regulation of B-cell activation, negative regulation of the inflammatory response, and leukocyte migration (b).
Figure 5. Functional enrichment analysis using FunRich. Proteins identified in each the different groups were submitted to molecular functional analysis according to their involvement in the immune response: proteinsrelated tointerleukins (ILs) IL-1, IL-1β, IL-6, IL-8, and interferon γ (a); proteins involved in immunoglobulin production, positive regulation of B-cell activation, negative regulation of the inflammatory response, and leukocyte migration (b).
Ijms 22 06165 g005
Figure 6. Functional enrichment analysis using FunRich. Proteins identified in each of the different groups were submitted to molecular functional analysis according to their involvement in immune processes: proteins involved in the classical or alternative complement activation pathways (a); proteins involved in the response to oxidative stress (b).
Figure 6. Functional enrichment analysis using FunRich. Proteins identified in each of the different groups were submitted to molecular functional analysis according to their involvement in immune processes: proteins involved in the classical or alternative complement activation pathways (a); proteins involved in the response to oxidative stress (b).
Ijms 22 06165 g006
Figure 7. Volcano plot showing the results of the SWATH analysis of plasma protein expression. Data are represented as the results obtained for the control group relative to all other groups (UG, ERT-a, and ERT-b). Graph shows the negative base 10 log of the p-value plotted against the base 2 log of the fold change for each protein included in the SWATH analysis. Changes were considered significant at p < 0.05 and fold change >1.5 compared with control group.
Figure 7. Volcano plot showing the results of the SWATH analysis of plasma protein expression. Data are represented as the results obtained for the control group relative to all other groups (UG, ERT-a, and ERT-b). Graph shows the negative base 10 log of the p-value plotted against the base 2 log of the fold change for each protein included in the SWATH analysis. Changes were considered significant at p < 0.05 and fold change >1.5 compared with control group.
Ijms 22 06165 g007
Figure 8. Heatmap depicting the expression profiles of different proteins (also shown in Table 2, Table 3 and Table 4) as determined by SWATH quantitative analysis. The relative abundance of each protein is reflected by the colors in the heatmap, based on the z-score of the protein’s normalized peak area. The dendrogram represents the biological replicates of the UG, ERT-a, ERT-b, and CG groups. Proteins, listed in the sidebar on the right, are color coded according to their function, as follows: apolipoprotein-related proteins, dark green; blood coagulation, purple; calcium binding, dark orange; carrier proteins, light green; cytoskeletal proteins, light pink; immune system, grey; metabolic processes, yellow; metabolic interconversion, light blue; other functions, white; protein binding, pink; transport, light grey.
Figure 8. Heatmap depicting the expression profiles of different proteins (also shown in Table 2, Table 3 and Table 4) as determined by SWATH quantitative analysis. The relative abundance of each protein is reflected by the colors in the heatmap, based on the z-score of the protein’s normalized peak area. The dendrogram represents the biological replicates of the UG, ERT-a, ERT-b, and CG groups. Proteins, listed in the sidebar on the right, are color coded according to their function, as follows: apolipoprotein-related proteins, dark green; blood coagulation, purple; calcium binding, dark orange; carrier proteins, light green; cytoskeletal proteins, light pink; immune system, grey; metabolic processes, yellow; metabolic interconversion, light blue; other functions, white; protein binding, pink; transport, light grey.
Ijms 22 06165 g008
Figure 9. Box plot depicting expression levels of the6 most relevant plasma proteins (ALBU, HBA, FETUA, APOC1, APOC2, VTNC, and GPX3) for which downregulation was observed in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median protein expression value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 9. Box plot depicting expression levels of the6 most relevant plasma proteins (ALBU, HBA, FETUA, APOC1, APOC2, VTNC, and GPX3) for which downregulation was observed in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median protein expression value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Ijms 22 06165 g009
Figure 10. Box plot depicting expression levels of the 5 most relevant plasma proteins (RET4, SAA4, FA10, FINC, and C1QB) for which very low expression was observed in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 10. Box plot depicting expression levels of the 5 most relevant plasma proteins (RET4, SAA4, FA10, FINC, and C1QB) for which very low expression was observed in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Ijms 22 06165 g010
Figure 11. Box plot depictingthe4 plasma proteins (A2GL, CLUS, GELS, ITIH2) with the highest expression levels in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 11. Box plot depictingthe4 plasma proteins (A2GL, CLUS, GELS, ITIH2) with the highest expression levels in the UG and ERT-a and ERT-b groups relative to the control group (CG). Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05; ** p < 0.01; *** p < 0.001.
Ijms 22 06165 g011
Figure 12. Volcano plot showing the results of the SWATH analysis of plasma protein expression. Data are represented as the results obtained for the untreated group (UG) relative to all other groups (CG, ERT-a, and ERT-b). Graph shows the negative base 10 log of the p-value plotted against the base 2 log of the fold change for each protein included in the SWATH analysis. Changes were considered significant at p < 0.05 and fold change >1.5.
Figure 12. Volcano plot showing the results of the SWATH analysis of plasma protein expression. Data are represented as the results obtained for the untreated group (UG) relative to all other groups (CG, ERT-a, and ERT-b). Graph shows the negative base 10 log of the p-value plotted against the base 2 log of the fold change for each protein included in the SWATH analysis. Changes were considered significant at p < 0.05 and fold change >1.5.
Ijms 22 06165 g012
Figure 13. Box plot depicting ATA1 levels in each of the study groups. Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05.
Figure 13. Box plot depicting ATA1 levels in each of the study groups. Each data point represents the median value obtained in an individual sample. The line inside the box represents the median of all values obtained. The upper and lower limits of the box represent the first and third quartiles. Whiskers represent the minimum and maximum values within 1.5 times the interquartile range. Any data points not included between the whiskers are considered outliers. * p < 0.05.
Ijms 22 06165 g013
Figure 14. Graphical abstract showing study workflow. Plasma samples were acquired from 3 distinct groups: untreated MPS IVA patients; MPS IVA patients that received ERT (from whom samples were taken 1 day before and 1 day after ERT), and healthy blood donors. Blood was separated to obtain plasma, the samples were lysed, and proteomic analyses were performed. On completion of these two distinct proteomic analyses, a bioinformatics analysis was performed to obtain more information about the identified proteins.
Figure 14. Graphical abstract showing study workflow. Plasma samples were acquired from 3 distinct groups: untreated MPS IVA patients; MPS IVA patients that received ERT (from whom samples were taken 1 day before and 1 day after ERT), and healthy blood donors. Blood was separated to obtain plasma, the samples were lysed, and proteomic analyses were performed. On completion of these two distinct proteomic analyses, a bioinformatics analysis was performed to obtain more information about the identified proteins.
Ijms 22 06165 g014
Figure 15. Enzyme activity of GALNS. CG, control group; ERT-a, MPS IVA patients sampled 1 day before ERT; ERT-b, MPS IVA patients sampled 24 h after ERT; UG, untreated group.
Figure 15. Enzyme activity of GALNS. CG, control group; ERT-a, MPS IVA patients sampled 1 day before ERT; ERT-b, MPS IVA patients sampled 24 h after ERT; UG, untreated group.
Ijms 22 06165 g015
Table 1. Changes in plasma protein expression detected by SWATH in any of the patient groups (UG, ERT-a, or ERT-b) relative to the control group.
Table 1. Changes in plasma protein expression detected by SWATH in any of the patient groups (UG, ERT-a, or ERT-b) relative to the control group.
Group ComparisonProteins with Altered Expression*Upregulated ProteinsDownregulated Proteins
UG vs CG14122
ERT-a vs CG85(2 only in ERT-a)1
ERT-b vs CG1410(2 only in ERT-b)1
CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group. * p < 0.05 and fold change >1.5 compared with the control group.
Table 2. Changes in plasma protein expression detected in all patient groups (UG, ERT-a, and ERT-b) relative to the control group. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Table 2. Changes in plasma protein expression detected in all patient groups (UG, ERT-a, and ERT-b) relative to the control group. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Uniprot CodeProtein
Code
Protein Namep-ValueFold Change
Downregulated in UG, ERT-a, and ERT-b
P01860IGHG3Immunoglobulin heavy constant gamma 33.30 × 10−5Only CG
P02741CRPC-reactive protein8.93 × 10−3Only CG
P02760AMBPProtein AMBP7.75 × 10−50.0008
Q92485ASM3BAcid sphingomyelinase-like phosphodiesterase 3b4.80 × 10−60.0021
P01599KV117Immunoglobulin kappa variable 1–176.69 × 10−40.0027
Q9P1Z0ZBTB4Zinc finger and BTB domain-containing protein 42.81 × 10−70.0086
P01780HV307Immunoglobulin heavy variable 3–78.23 × 10−60.0333
P22352GPX3Glutathione peroxidase 31.31 × 10−40.0349
P35542SAA4Serum amyloid A-4 protein8.01 × 10−50.0368
P01766HV313Immunoglobulin heavy variable 3–134.50 × 10−20.0402
P04004VTNCVitronectin1.00 × 10−60.0417
P00742FA10Coagulation factor X4.71 × 10−90.0535
P02042HBDHemoglobin subunit delta2.79 × 10−50.0565
A0A075B6J9LV218Immunoglobulin lambda variable 2–188.30 × 10−40.0588
P02753RET4Retinol-binding protein 46.31 × 10−30.1103
P06312KV401Immunoglobulin kappa variable 4–12.12 × 10−60.1292
P02654APOC1Apolipoprotein C-I7.01 × 10−60.1515
P02765FETUAAlpha-2-HS-glycoprotein6.60 × 10−50.1527
P69905HBAHemoglobin subunit alpha1.86 × 10−50.1585
P02655APOC2Apolipoprotein C-II8.27 × 10−50.1642
P00915CAH1Carbonicanhydrase 11.37 × 10−20.1984
P02751FINCFibronectin9.79 × 10−30.2304
A0A0B4J1Y8LV949Immunoglobulin lambda variable 9–496.70 × 10−30.2353
P02746C1QBComplement C1q subcomponent subunit B1.68 × 10−30.2381
P00738HPTHaptoglobin3.00 × 10−30.2445
P01857IGHG1Immunoglobulin heavy constant gamma 12.37 × 10−20.2519
P01834IGKCImmunoglobulin kappa constant4.67 × 10−40.2591
P68871HBBHemoglobin subunit beta1.01 × 10−30.2681
P01859IGHG2Immunoglobulin heavy constant gamma 23.60 × 10−40.2959
P29622KAINKallistatin1.99 × 10−30.3344
P02743SAMPSerum amyloid P-component1.91 × 10−20.3509
P02768ALBUSerum albumin1.92 × 10−50.3509
P02787TRFESerotransferrin4.11 × 10−20.3571
P02749APOHBeta-2-glycoprotein 11.26 × 10−20.3584
P02679FIBGFibrinogen gamma chain5.71 × 10−30.3597
P02763A1AG1Alpha-1-acid glycoprotein 12.43 × 10−20.4115
P05543THBGThyroxine-binding globulin4.82 × 10−40.4348
Upregulated in UG, ERT-a, and ERT-b
P01619KV320Immunoglobulin kappa variable 3–200.0005Not found CG
P80748LV321Immunoglobulin lambda variable 3-210.0179Not found CG
B9A064IGLL5Immunoglobulin lambda-like polypeptide 50.000244.0529
P01615KVD28Immunoglobulin kappa variable 2D–280.047427.4725
P0DOY2IGLC2Immunoglobulin lambda constant 20.021726.1780
P35858ALSInsulin-like growth factor-binding protein complex acid labile subunit0.029910.8460
P10909CLUSClusterin0.001010.1523
P02750A2GLLeucine-rich alpha-2-glycoprotein0.00049.9800
P01023A2MGAlpha-2-macroglobulin0.00004.7170
P00734THRBProthrombin0.01324.3178
P06396GELSGelsolin0.00573.3223
P19823ITIH2Inter-alpha-trypsin inhibitor heavy chain H20.00233.1536
Q96IY4CBPB2Carboxypeptidase B20.01212.8977
Q08380LG3BPGalectin-3-binding protein0.03612.5880
P08185CBGCorticosteroid-binding globulin0.03252.2763
A0A0A0MS15HV349Immunoglobulin heavy variable 3–490.04822.1418
A0A0C4DH38HV551Immunoglobulin heavy variable 5–510.01652.1345
Only CG: These proteins were not detected in theUG, ERT-a, or ERT-b groups, and were only found in CG. Not found CG: these proteins were not found in CG, and were only detected in the UG, ERT-a, or ERT-b groups.CG, control group; ERT-a, patients sampled before ERT; ERT-b, patients sampled 24 h after ERT; UG, untreated group.
Table 3. Changes in protein expression in the untreated group relative to all other groups. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Table 3. Changes in protein expression in the untreated group relative to all other groups. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Uniprot CodeProtein CodeProtein Namep-ValueFold Change
Upregulated in Untreated Group
P06702S10A9Protein S100-A90.003617.4468
P01615KVD28Immunoglobulin kappa variable 2D-280.00085.6619
P80748LV321Immunoglobulin lambda variable 3-210.00074.6098
P01594KV133Immunoglobulin kappa variable 1-330.01823.5490
P02652APOA2Apolipoprotein A-II0.01972.8151
P01591IGJImmunoglobulin J chain0.03952.4479
P63261ACTGActin, cytoplasmic 20.03042.2374
P00734THRBProthrombin0.02652.0770
B9A064IGLL5Immunoglobulin lambda-like polypeptide 50.04201.9341
P01714LV319Immunoglobulin lambda variable 3-190.03131.7496
P01009A1ATAlpha-1-antitrypsin0.00141.6875
P02747C1QCComplement C1q subcomponent subunit C0.02231.6522
Downregulated in Untreated Group
P02746C1QBComplement C1q subcomponent subunit B0.03290.1549
P08697A2APAlpha-2-antiplasmin0.02750.4938
Table 4. Changes in protein expression in the ERT-a or ERT-b groups relative to all other groups. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Table 4. Changes in protein expression in the ERT-a or ERT-b groups relative to all other groups. Protein expression is considered altered in cases in which p < 0.05 and the fold change >1.5 (upregulated proteins) or <0.6 (downregulated proteins).
Uniprot CodeProtein CodeProtein Namep-ValueFold Change
Upregulated in ERT-a
P0DJI8SAA1Serum amyloid A-1 protein0.0487Only ERT-a
Q13790APOFApolipoprotein F0.0487Only ERT-a
A0A0C4DH68KV224Immunoglobulin kappa variable 2-240.01096.1248
P05156CFAIComplement factor I0.00655.2409
A0A075B6I0LV861Immunoglobulin lambda variable 8-610.03863.1547
P43251BTDBiotinidase0.03012.8996
P04196HRGHistidine-rich glycoprotein0.02241.5688
Downregulated in ERT-a
P02747C1QCComplement C1q subcomponent subunit C0.04650.5069
Uniprot CodeProtein CodeProtein Namep-ValueFold Change
Upregulated in ERT-b
P60174TPISTriosephosphateisomerase0.0487Only ERT-b
O00391QSOX1Sulfhydryl oxidase 10.0487Only ERT-b
A0A0C4DH67KV108Immunoglobulin kappa variable 1-80.0028Only ERT-b
P0C0L5CO4BComplement C4-B0.002917.9592
O14791APOL1Apolipoprotein L10.04802.2700
P20742PZPPregnancyzone protein0.04182.1171
P06396GELSGelsolin0.00802.0649
P00751CFABComplement factor B0.01842.0058
P19823ITIH2Inter-alpha-trypsin inhibitor heavy chain H20.00931.9167
P01023A2MGAlpha-2-macroglobulin0.00931.8574
P08697A2APAlpha-2-antiplasmin0.02981.8165
Q14624ITIH4Inter-alpha-trypsin inhibitor heavy chain H40.03891.7428
P05546HEP2Heparin cofactor 20.04211.6111
Downregulated in ERT-b
P02766TTHYTransthyretin0.02240.2358
Only ERT-a: these proteins wereupregulated in the ERT-a group with an FC value of infinity, meaning that they wereabsentin all other groups (UG, CG, and ERT-b). Only ERT-b: These proteins wereupregulated in the ERT-b group with aFC value of infinity, meaning that they wereabsent in all other groups (UG, CG, and ERT-a).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Álvarez, J.V.; Bravo, S.B.; Chantada-Vázquez, M.P.; Barbosa-Gouveia, S.; Colón, C.; López-Suarez, O.; Tomatsu, S.; Otero-Espinar, F.J.; Couce, M.L. Plasma Proteomic Analysis in Morquio A Disease. Int. J. Mol. Sci. 2021, 22, 6165. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22116165

AMA Style

Álvarez JV, Bravo SB, Chantada-Vázquez MP, Barbosa-Gouveia S, Colón C, López-Suarez O, Tomatsu S, Otero-Espinar FJ, Couce ML. Plasma Proteomic Analysis in Morquio A Disease. International Journal of Molecular Sciences. 2021; 22(11):6165. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22116165

Chicago/Turabian Style

Álvarez, José V., Susana B. Bravo, María Pilar Chantada-Vázquez, Sofía Barbosa-Gouveia, Cristóbal Colón, Olalla López-Suarez, Shunji Tomatsu, Francisco J. Otero-Espinar, and María L. Couce. 2021. "Plasma Proteomic Analysis in Morquio A Disease" International Journal of Molecular Sciences 22, no. 11: 6165. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22116165

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