Abstract
Lubricin or proteoglycan-4 (PRG-4) is a mucinous glycoprotein lubricating the cartilage and maintaining normal tissue function and cell-homeostasis. Altered O-glycoforms of lubricin has been found in osteoarthritis (OA) synovial fluid (SF). Here, we utilized a Fluorescent Immuno-Lectin Assay (FILA) to measure the levels of lubricin glycoforms in plasma/SF and their potential as biomarkers for disease. Five different lectins were used in the assay: Macrophage Galactose-type lectin (MGL), Sambucus Nigra Agglutinin (SNA), Maackia Amurensis Agglutinin (MAA), Peanut Agglutinin (PNA), and Galectin-3 (Gal-3). Our results showed that the levels of lubricin glycoforms in late-stage knee OA plasma (typically <10 μg/ml) are 1-3 order of magnitude lower than in OA SF (typically> 100 μg/ml). Furthermore, plasma lubricin glycans displayed higher level of sialylation, while SF derived lubricin displayed higher level of Tn-antigens. In OA patients we found decreased level in plasma of SNA binding lubricin (p=0.0023) compared to controls. This lectin binds 6 linked sialic acid. In addition, we found that lubricin glycoforms correlated both with Body Mass Index (BMI) and age, especially in regards of sialylation as measured by both MAA and SNA. Our data suggest that glycosylation of lubricin is different comparing SF with plasma. Moreover, the glycosylation of plasma lubricin is altered in OA patients compared controls.
INTRODUCTION
Osteoarthritis (OA) is a degenerative joint disease where the protective cartilage that cushions the ends of the bones wears down and the underlying bone changes due to imbalance in cartilage catabolism1,2. The OA condition affects joints such as hands, knees, hips, and/or spine1,3. The seriousness of the disease is illustrated by that according to the World Health Organization (WHO) 9.6% of all 60+ years old men and 18% of all 60+ years old women are estimated to have symptomatic OA, worldwide4. The development of OA is believed to start as an initial alteration in pericellular matrix and softening of the extracellular matrix, with an increased mechanical stress of the chondrocyte and increased secretion of proteolytic enzyme degrading type II collagen and proteoglycans5-7. It is this process that eventually leads to a macroscopically detectable degradation of cartilage8,9. In addition to this pronounced degradation of the cartilage, OA also involves a low-grade joint inflammation, which manifests as an increased local level of inflammatory cytokines, chemokines, and growth factors10. At the present time, the diagnosis of OA is limited by radiographic evaluation, clinical examination, and detection of signs/symptoms such as crepitus, morning stiffness of joints, joint pain, joint instability, and bony enlargement (osteophytes)3,11. There are no approved molecular diagnostics for OA. In an early non-chronic stage of OA, there is a potential that diagnostics will be able to identify individuals that will benefit from future disease modifying osteoarthritis drugs (DMOADs) or can be selected to less invasive therapy compared to total joint replacement surgery11,12. Hence, there is a need to use modern type of omics methodology to identify biomarker candidates that could aid in identifying subgroups of patients that are at risk for developing OA and would benefit from precision medicine targeting their specific phenotype of the disease. We and others have pursued the use of glycomics for stratification of OA patients, specifically looking into the O-glycome of the lubricin in synovial fluid (SF)13. The SF provides the lubrication to the joint giving an almost friction free motion. There are mainly three types of biomolecules that have been implicated being involved in this lubrication: phospholipids14, hyaluronan (HA), and glycoproteins (lubricin15 and aggrecan16).
Lubricin (Figure 1) is the key O-glycoprotein responsible for boundary lubrication of diarthrodial joints and for the integrity of articular cartilage15,17. Lubricin has earned the attention for its involvement in joint degrading diseases such as OA and rheumatoid arthritis18. In the articular cartilage, lubricin is expressed by synoviocytes located on synovial membrane19 and chondrocytes embedded in cartilage matrix20. In addition to the joints, lubricin is widely distributed in the body, being expressed also by cells in the liver, heart, lungs, brain, prostate, small intestine, joints and tendon21,22. The expression of lubricin by various cells, alternative splicing, and post-translational modifications give rise to different isoforms of lubricin23. In accordance with lubricin’s wide distribution, the protein has been implicated in other non-joint associate diseases, where its regulation in plasma and liver in a murine sepsis model is in line with the role of lubricin as an immune- and inflammatory-regulator24,25.
Lubricin is heavily glycosylated in its mucin-like domain with O-linked glycans26, that counts for > 50% of its mass27,28. Lubricin in its largest splice form is translated into 1404 amino acid residues, where the STP-rich domain consists of several imperfect repeats of the amino acid sequence KEPAPTTP13,29. The glycosylation of lubricin occurs on its STP-rich domain, that have shown to harbor more than 160 binding sites for O-linked glycosylation including sialylated/non sialylated core-1 O-glycans (Galβ1-3GalNAcα1-), and core-2 structures such as NeuAcα2-3Galβ1-3(NeuAcα2-3Galβ1-4GlcNAcβ1-6)GalNAcα1-28,29. In OA, there is an increase of truncated structures such as Tn-antigen (GalNAcα1-Ser/Thr) in SF lubricin10. Lubricin with truncated glycosylation acts as a signaling molecule that stimulates synoviocytes influencing the level of inflammatory cytokines, chemokines, and growth factors10.
In this report we set out to investigate the different glycoforms of lubricin found in SF10 and plasma, in respect to OA. In order to detect differences in glycosylation we developed and validated a Fluorescent Immuno-Lectin Assay (FILA), to ensure that we had the sensitivity required.
Results
Overview results
We selected five lectins including Macrophage Galactose-type lectin (MGL, GalNAcα-1, Tn-antigen), Sambucus Nigra Agglutinin (SNA, NeuAcα2-6, only found as NeuAcα2-6GalNAc in lubricin)23, Maackia Amurensis Agglutinin (MAA, NeuAcα2-3Galβ1-), Peanut Agglutinin (PNA,Galβ1-3GalNAcα1-, T-antigen, core 1) and Galectin-3 (Gal-3, Galβ1-4GlcNAc-, N-acetyllactosamine type 2), corresponding to glycoepitopes previously found on lubricin13, and attached them covalently to fluorescent nanoparticles. These lectin-nanoparticles were then used to measure the amount of glycoepitopes on lubricin captured by mAbs against lubricin (Figure 2A). We validated this FILA assay based on recommended tests such as spike-and-recovery, dilution-recovery, and inter-/intra-plate variation30,31. In addition, we benchmarked the FILA assay using SF from OA patients and compared it to a quantitative mass spectrometric Selected Reaction Monitoring (SRM) method for individual glycans from lubricin (Figure 2B-C). The validated method was then used for comparing lubricin glycosylation in plasma versus SF as well as comparing lubricin glycosylation in plasma from OA-patients and controls.
Validation of FILA for screening of plasma lubricin glycosylation
For the validation of the lectin immunoassay we performed a spike-and-recovery test using a low (1x), medium (3x), and high spike (9x) (Table 1). Overall, for the range tested and for all lectins we found a good recovery (92%). In general, we found that the recovery was lower the higher the spike concentration, with an average of 117% (1x) versus 72% (9x), indicating a suppression of signal for higher concentration.
Extreme high value of recovery for 1x was found for MAA (189%), and extreme low values of recovery were found for PNA and MGL with high suppression of the 9x spike of 58% and 47%, respectively. To test the linearity of dilution, we performed an 8 times 1:2 serial dilution of plasma samples starting from 1:150 for all lectins except MGL (1:3 instead) and we found a recovery of 88% - 122%. Intra-plate reproducibility was found to be slightly lower (14% CV) compared to Inter-plate (19%) (Table 1).
Benchmarking FILA with quantitative glyco SRM-MS
Lubricin glycans in SF are predominantly made up of simple O-linked glycans, consisting of Tn antigens, core 1 and core 2 (Galβ1-3[GlcNAcβ1-6]GalNAc with/without sialic acid and sulfate. The whole glycan set would comprise of a total of 15 structures, where three pairs are sequence isomers23. This relatively small number allows screening of the O-linked glycans released from lubricin to obtain quantitative information by mass spectrometry using LC separation in combination with SRM as we have shown previously23 (Figure 2B). We used LC-SRM to benchmark the FILA for lubricin glycoforms (Figure 2C). For this we used SF from OA patients, and isolated lubricin by ion exchange chromatography before the glycans were released. The variation of lubricin glycan structure found in OA patient’s synovial fluid was used to correlate the response between LC-SRM and the FILA. We usually found a good agreement between the two assays (Figure 2C), with a regression for all SRM/FILA above 0.6. The lowest regression was found for SNA, that also displayed the lowest variation in regards of level between the different samples, both found by SRM and FILA. The low correlation was probably a reflection of that the level of SNA-glycoform did not vary enough between our samples, superseded by the precision of the measurement.
Lubricin glycoforms in OA patients’ plasma are different compared their synovial lubricin glycoforms
We performed the FILA using the five lectins and selected matched plasma/SF from 19 OA patients (Figure 3A). What was obvious from the analysis was that for all glycoforms, the level in SF was 1-3 order of magnitudes higher compared to plasma. For most of the glycoforms the levels in plasma were found to be below 10 μg/ml, while in SF glycoforms were usually detected at above 100 μg/ml, except for the glycoform detected by Gal-3 that was slightly lower. The highest difference in the glycoform levels between SF and plasma were found using the MGL-lectin (Tn-antigen), where the levels in plasma were less than 1 μg/ml. To appreciate differences in lubricin glycosylation between SF and plasma, we generated the ratio of all glycoforms detected with MAA, PNA, SNA and Gal-3 against glycoforms detected by MGL (Figure 3B). We found that for all epitopes, the ratio was higher in plasma compared to SF, and thus highlighting that the glycosylation of plasma derived lubricin is markedly different to SF lubricin.
Plasma lubricin glycoforms in OA and control
To identify differences in lubricin glycoforms in plasma between controls and OA patients, we selected non-obese (BMI<30) OA patients due for Total Knee Replacement from two orthopedic clinics in Sweden (N=183) and controls (N=38). The selection criteria for controls were individuals with no self-reported history of recuring joint problems and not diagnosed with having previous OA. Females were overrepresented in OA patients (59%) reflecting the increased incidence of OA in females in the population. However, no significance was detected in the gender, age and BMI distribution between our OA patients and controls.
Using FILA, SNA, MAA, PNA and Gal-3, the concentration of glycoforms of lubricin was found to be around 10 μg/ml, while the MGL glycoform was found to be 1-4 order of magnitudes lower (Figure 4A). The statistic evaluation shows that plasma lubricin levels of lectin-binding glycoforms (MAA, Gal-3, PNA, and MGL), in OA and controls were not significant showing p-values ranging between 0.20 to 0.74 (Supplementary Table 3). Sialyl-Tn-antigen and related structure with 6 linked sialic acid detected by SNA lectin showed a decrease in OA (p=0.0023). We also performed the calculation of glycoform ratio dividing the individual lectin concentration using FILA with the MGL lectin. Also using this measure, we found that the ration between SNA/MGL was significantly different (p=0.0178), while the other lectin ratios were not (Figure 4B). This data suggest that OA patients display a decreased level of plasma lubricin glycoforms containing epitopes for SNA compared to controls.
Lubricin glycoforms in plasma correlate with increased BMI and age
We also addressed if there was a correlation between lubricin glycoforms and BMI. Selecting both obese and non-obese patients and controls from the biobank, significant positive correlations between BMI and lubricin glycoforms was detected using lectins Gal-3 (p=0.016), MAA (p=0.0033), and SNA (p=0.0018), but not for MGL and PNA (Figure 5). These data indicates that BMI/overweight is associated with increased sialylation (MAA and SNA) and complexity (Type 2 N-acetyllactosamine elongation, Gal-3) of lubricin.
We also correlated lubricin glycoforms in plasma with age. Significant negative correlations between age and lubricin glycoforms in our OA patients was detected using lectins MAA (p=0.043), and SNA (p=0.0047) and positive correlation for MGL (p=0.017). Gal-3 and PNA did not show statistical correlation with age (Figure 6).
Interestingly, the positive impact of BMI in regards of sialylation, measured both by MAA and SNA appeared to be counteracted by age, since there was a slight negative correlation between both SNA and MAA and age.
Discussion
We validated the lectin FILA according to selectivity, accuracy and precision and dilution linearity32, and cross validated the FILA with SRM of oligosaccharides from lubricin in synovial fluid using patient material. Overall and for individual lectins we could identify that we were within or close to recommendation (±20%) for spike and recovery, dilution recovery and inter-/intra-plate CV32. The data also showed good agreement between SRM and FILA. These validation data provided us with confidence that our FILA assay was providing a good picture of lubricin glycosylation in plasma.
Our results demonstrated that both the level and the glycosylation of lubricin from plasma and SF is different. Our FILA data suggested levels of lubricin of >100 μg/ml in SF and < 10 μg/ml in plasma. These concentrations are consistent with what have been reported in human SF and plasma previously33-35. The glycosylation of SF lubricin has been thoroughly characterized, demonstrating predominant level of Tn-antigen and un-sialylated/mono-sialylated core 1 glycans and lower levels of disialylated structures and core-2 glycans13,23, consistent with the FILA data. The FILA data from plasma suggest that lubricin has higher level of sialylation, both 3- and 6-linked relative to SF lubricin (Figure 3). The relative increased binding to PNA (core 1, T-antigen) and Gal-3 (core 2, Type-2 N-acetyllactosamine) is on the expense of a relative decrease of Tn-antigen in plasma compared to SF (MGL-binding) (Figure 3 B). Both core-1 and core-2 structures in plasma lubricin are then the target for the increased sialylation found in plasma compared to SF. The difference in glycosylation of lubricin between the two fluids suggest that the lubricin in plasma are mainly made up of lubricin secreted from different tissue origin(s) than SF lubricin. The cells that are synthesizing SF lubricin are expected to be found in the joint area such as synoviocytes and chondrocytes, while plasma lubricin would originate from high expressing tissue such as liver36 that will secrete the protein into the blood. Due to the high content of Tn-antigen on SF lubricin shown here, it would be predicted that Tn-containing lubricin escaping the joint capsule into the blood would have a short biological half-life, quickly removed by the Asialoglycoprotein receptor (ASGPR) expressed in liver37 . Hence, it is possible that lubricin SF glycoforms escaping the synovial areas may have different fate as they start to circulate in the plasma compared to glycoforms with less Tn antigens and higher sialylation synthesized by other cells.
The change of plasma glycosylation in OA with decreased amount of SNA epitopes is quite different compared to OA regarding SF lubricin where the glycans are shown to be more truncated and expected to have more Tn antigens compared to controls10. The origin of this lubricin SNA depleted glycoform found in the OA-cohort ‘s plasma is unknown and whether it is caused by alterations in catabolic and/or anabolic processes. Interestingly, the data from the comparison between SF and plasma lubricin (Figure 3B) indicates that SF lubricin has a relative lower content of SNA epitopes compared to plasma lubricin. The small change in SNA binding indicates that this alteration of lubricin in plasma could be included in a combined assay of biomarkers for OA diagnostic/prognostic, but it would not suffice by itself.
Lubricin level in plasma has been suggested to correlate with Joint Space Narrowing (JSN)38, but so far, a specific glycoform of lubricin has not been shown to be related to OA. Further differentiation between plasma and SF derived lubricin would be possible with additional molecular selection. It has been shown that the splice-form of lubricin produced by hepatocytes is shorter compared to the SF form, with difference in the deletion of a peptide segment in the non-glycosylated region22,39. Identifying splice form of lubricin found in plasma may be an additional criterion to trace the origin of lubricin in plasma. The biomedical background of the potential change of glycoforms in OA plasma, needs further investigation, whether it is originating from a local inflammation of the joint, or it is due to a more profound systemic alteration of the health status in our patients that were all from the late stage of OA.
The correlation between lubricin lectin binding in plasma and BMI, suggests that metabolic syndrome could a main driver of altered O-linked glycosylation of lubricin in plasma, especially the increase of sialic acid with BMI. A potential connection between lubricin and BMI is a systemic inflammation. In this case the inflammation is associated with overweight40,41. In other causes of cause of inflammation, such as an an acute sepsis model in mice, hepatic lubricin displays the highest upregulation of all proteins42. This may suggest that lubricin has an acute phase protein like behavior. The influence of O-linked glycosylation remodeling in hepatic acute phase glycoproteins has not been widely studied. However, remodeling of N-glycosylation of acute phase proteins has been reported in both chronic and acute inflammation43.
In regards of aging, the plasma lubricin glycosylation was indicated to display decreased sialylation (MAA and SNA) with a concomitant upregulation of the Tn antigen (MGL). Interestingly SNA (p=0,0023) and MAA (p=0.20) was the lectins that displayed the highest significance in OA versus control plasma. This indicates that aging and OA manifested similarly in regards of plasma lubricin glycoform alteration, in the opposite direction to BMI, where sialylation appeared to increase with increased BMI. The relation of O-linked glycosylation and age in plasma has not been widely studied. Studies on plasma N-linked glycosylation and age shows that the effect varies depending on protein (and potentially origin)44.
We conclude that using glyco-analytics, we have been able to trace that lubricin in plasma and SF are of different , since there are significant differences in the glycosylation. Hence, the glycosylation alteration of lubricin as we found in OA SF10,23would have limited impact on the glycosylation of lubricin in plasma, and we only found glycoforms detected by SNA to be altered in OA plasma. Our data also suggest that lubricin in plasma is altered by age and in metabolic syndrome manifested as increased BMI. Lubricin has been indicated to be involved in both atherosclerosis45 and in high-fat-diet induced glycose intolerance46. Changed glycosylation of lubricin in these diseases are likely to contribute to the disease progress.
Materials and Methods
Ethics declarations and clinical samples
Plasma samples and their controls were from a biobank maintained by Sahlgrenska University hospital (Gothenburg, Sweden). All patients and participants had given informed consent. Ethical permission for patient screening was obtained from the ethical board in Gothenburg, Sweden (ethical application 172-15). Analysis of OA biobank samples has been approved by REK (470919) and handling of personal data by NSD (431575). The study was performed using samples from 183 OA patients and 41 controls. The information about the patients is summarized in Table 2.
MS-FILA cross validation of lubricin glycosylation
Lubricin from OA patients (n = 20) were isolated from SF samples by anion exchange chromatography as previously described23. Lubricin containing fractions were concentrated using 100 kD spin filters (0.5 ml Amicon Ultra, Millipore) and salt exchange was performed with 3 x 0.5 ml 0.1M NH4HCO3, followed by speedvac to dryness. The standard mixture contained four standards: GalNAc (Tn) was from Sigma-Aldrich (St Louis MO, US), Galβ1-3GalNAc (T) from Dextra (Reading, UK), and NeuAcα2-3Galβ1-3GalNAc (ST) and Sulfo Lea (HSO3-3Galβ1-3(Fucα1-4)GlcNAc) were purchased from Carbosynth (Berkshire, UK). The O-linked oligosaccharides were released as alditols by reductive β-elimination in 100 μl sodium borohydride (1.0 M) in sodium hydroxide (0.10 M) at 50°C overnight followed by cleanup using 150 μl cation exchange media (AG50WX8, Biorad, Hercules, CA, US) on top of C18 SPE columns (Strata C18-E, 100 mg, Phenomenex, Torrance, CA, US). The oligosaccharides were dried in the speedvac, followed by repeated additions of 5 x 50 μl methanol and dried in speedvac, to evaporate borate methyl esters. Glycan standards were reduced to alditols at same conditions as above for 3 hours or overnight. The oligosaccharides were dissolved in 100 μl H2O, followed by injection (8 μl) on a Waters UPLC-MS/MS (Acquity Xevo TQ-S) triple quadrupole mass spectrometer. Oligosaccharides were separated on porous graphitized carbon columns (100 x 2.1 mm, 3 µm particles, Hypercarb, ThermoFisher Scientific, Waltham, MA) kept at 25°C. The gradient (36 min) consisted of 0-20 min 0-40% B (A: 10 mM ammonium bicarbonate, B: 80% acetonitrile in 10 mM ammonium bicarbonate), 20-23 min 40-100% B, 24-26 min wash with 1% HAc, then equilibration 26-36 min with 100% A. The flow rate was kept at 150 µl/min. The ESI capillary was kept on 2.5 kV. The source was at 150°C, and the cone at 40V. The samples were injected twice to cover a large range of glycans. The instrument was run in positive mode the first 4.35 min for analysis of GalNAcol (elution time 4 min). Positive mode transitions were the same for injections analyses, and covered the following six transitions at collision energy (CE)=30% and with a dwell time of 0.052 ms per transition during every cycle: m/z 224 to 182 ([M+H]+-C2H2O), m/z 224 to 206 ([M+H]+-H20) and also the corresponding C13 isotope transitions: m/z 225 to 183, m/z 225 to 207, m/z 246 to 182 ([M+Na]+-C2H2O) and m/z 246 to 206 ([M+Na]+-H2O). During the second period of the chromatographic run the instrument was run in negative mode (4.35-20 min), covering nine transitions per run and at a dwell time of 0.032 ms: m/z 384.2 to 204.1 (CE=15); m/z 464.1 to 241.1 (CE=40); m/z 464.1 to 302.1 (CE=40); m/z 587.2 to 384.2 (CE=20); m/z 587.2 to 407.2 (CE=20); m/z 610.2 to 464.2 (CE=40) ; m/z 667.2 to 241.1 (CE=45); m/z 667.2 to 505.2 (CE=45); m/z 675.2 to 290.1 (CE=30); m/z 749.3 to 569.2 (CE=30); m/z 755.2 to 465.2 (CE=50); 829.2 to 667.2 (CE=50); m/z 878.3 to 290.1 (CE=40); m/z 966.3 to 290.1 (CE=50);m/z 1040.4 to 290.1 (CE=45) and m/z 1331.5 to 290.1 (CE=60). Collision energies were optimized using syringe infusions of saliva mucin oligosaccharides, fetuin and porcine gastric mucin oligosaccharides. A dilution series (1:2) of the standard mixture was made (0.10-90 pmol on column), with analyses of 1-2 standards between every three samples. Standards were analyzed using the same method as for samples, although 2 μl were injected.
Biotinylation of antibodies
Sulfo N-hydroxysuccinimide Biotin (sulfo-NHS-biotin) (Thermos Fisher, Waltham, Massachusetts, US) was used for the biotinylating of antibodies. In-house mouse IgG antibodies 1E12 (0.42 mg, for plasma-SF comparison) and 14G10 (0.71 mg, for plasma OA-Control comparison) against mucin domain of lubricin was biotinylated in 2.5 ml of PBS with 20-fold molar excess of sulfo-NHS-biotin compared to antibody. The biotinylation reaction was allowed to proceed for two hours. Excess biotin was removed using a PD-10 column (Cytiva, Little Chalfont, UK) according to manufacturer’s recommendation.
Lectin coating of Eu3+ nanoparticles
Lectins (Vector labs, Newark, CA, US) were coated covalently attached to 97 nm Fluoro-Max Dyed Carboxylate-modified Microparticles (Thermo Fisher) (20 µl). The particles were applied to 1 Nanosep 300kDa Omega (Thermo Fisher) centrifugal filter. The Nanosep tube was centrifuged until dry at 8000 rpm for five minutes. Filter was washed twice with 200 µl of conjugation buffer (2-(N-morpholino)ethanesulfonic acid (MES) buffer, 50 mM, pH 6.0) and finally resuspended adding 150 µl of conjugation buffer followed by sonication with 40 pulses using an ultrasonic probe. Particles was activated adding sulfo-NHS (230 mM dissolved in conjugation buffer, 7.4 µl) and 1-ethyl-3-(3’dimethylaminopropyl)-carbodiimide (Sigma-Aldrich, 52 mM dissolved in conjugation buffer 10.4 µl) and agitated for 15 min at room-temperature using vortex. Lectins (0.5 mg/ml dissolved in PBS, 75 µl) were added to 50 µl of particle solution and the volume was adjusted to a total of 156 µl using 500 mM MES (4.5 µl), 3.5 µl 3.0 M NaCl, 17 µl water and 6.5 µl carbonate buffer 0. 50 M, pH 9.8 and incubated and repeatedly vortexed for another 30 min. The solution was stored over night at +4 °C in rotary mixer with the exception for SNA that was just stored in +4 °C without mixing. The reaction mixture was removed on the next day using one NanoSep 300 and centrifugated (8000 rpm) for five min followed by washing with storage buffer twice (100 μl) (25 mM Tris, 150 mM NaCl, 0.1% NaN3 pH 7.8) and centrifugation at 8000 rpm for five min. Storage buffer (100 µl) was added on the filter and the lectin coated particles was resuspend by ultra-sonication.
The fluorescent lectin-immune assay (FILA)
Streptavidin coated micro-plates (SA96-plates, Kaivogen, Turku, Finland) were prewashed 1x with washing buffer (Kaivogen, Turku, Finland). Biotinylated mAb 1E12 and 14G10 targeting lubricin was added to each well (25 ng/well, 25 μl) and plates were incubated for one hour in RT. After incubation the wells were washed 2x by washing buffer and plasma samples were added in triplicates. Each sample was diluted (1:150 and 1:400 or 1:5 and 1:10 for MGL) in red buffer. In addition, a serial dilution (1.3×10-1 µg/ml up to 2.5×10-4 µg/ml) of recombinant lubricin-(rhPRG4) stock conc.1300 µg/ml) were added to the plate in triplicates. Plates were then incubated for one hour at room temperature. After the incubation all plates were washed 2x and lectin Eu3+-nanoparticles (1×107 particle/well, 25 µl). The plates were incubated for one hour at RT and washed 6x with washing buffer. The fluorescence was read using VictorNivo multimode plate reader (PerkinElmer). The measurement operation was time-resolved fluorescence (TRF) with excitation 320 nm, and emission 615 nm wavelength. The concentration of lubricin glycoforms was calculated using MyAssays analysis software (MyAssays Ltd, Brighton, UK).
Validation of immune assay
For the validation of immune assay both spike and recovery test and test for the linearity of dilutions were performed. In addition, a system-control was used for calculation of intra- and inter-assay coefficient of variations (CVs) according to recommendation31.
For spike and recovery test, diluted plasma was used adding of known amounts of rhPRG4 (spike) to selected samples. Plasma diluted was spiked with low (0.55 µg/ml), medium (1.7 µg/ml), and high (5.0 µg/ml) concentrations using rhPRG4. For recovery 1:150, 1:300, 1:600, 1:1200 diluted plasma was used and added to the plate in triplicates. In case of MGL selected plasma samples were spiked with low (0.3 µg/ml), medium (1.0 µg/ml), and high (2.9 µg/ml) concentrations using rhPRG4 and for dilution recovery 1:3, 1:6, 1:12, 1:24 diluted plasma were used. Intra-plate CV% was determined for identical samples in triplicates from the same plate. Inter-plate CV% was determined from same system-control sample in different plates.
Statistics
In total 224 samples were used for analyzing plasma lubricin glycoforms and to determine the differences in the subgroups: OA (n=183) and controls (n=41). Statistical analyses were performed using statistical software: SPSS (IBM, Armonk, NY, US) and GraphPad Prism 9 (Dotmatics, CA, US). Continuous variables were summarized as mean and standard deviation or median and minimum-maximum. Assumptions such as outliers and normal distribution (Shapiro-wilk test) were considered, where outliers were kept. All data were normalized and non-parametric Mann-Whitney U-test (unpaired) or Wilcoxon matched pairs were used to determine if there were differences in glycoforms of lubricin between subgroups (for each glycoform separately). The significant level was considered as p<0.05 (two-tails).
Quantitative analysis of samples and standard curves were carried out regarding different types of lectins. For the calculation of MGL concentration in samples linear regression was used and for other lectins four parameter logistic (4PL).
Supplementary material
Supplementary Data 1: Raw data from FILA for all individuals
Supplementary Data 2: Raw data from FILA for individuals BMI<30
Supplementary Tables: Table 1: Comparison of FILA levels between OA and control plasma lubricin glycoforms in individuals with BMI<30 and for all individuals. Table 2: Correlation between lubricin glycoforms and BMI. Table 3: Correlation between lubricin glycoforms and age.
Data availability statement
Analytical data about individual patient FILA results are included in the Supplementary Material. Inquiries about additional patient meta data can be directed to the corresponding author.
Author Contributions
ARA and NGK planned, and ARA performed the validation of the FILA. ARA planned and executed the OA and control FILA experiments and performed the statistical evaluation of the results. JH, KG and KP adopted the FILA for lubricin glycosylation analysis and JH performed the plasma-SF comparison. KT and HR performed SRM experiments. OR, LIB and TE were responsible for providing control and patient sample material. TS and GJ contributed to the design and provision of rhPRG4 for the analysis and antibody generation. ARA wrote the initial draft, and NK completed the writing with input from all contributors. NGK designed the research plan, coordinated, and directed the research, and performed the overall statistical evaluation. NK and GJ generated funding for the research. All authors critically reviewed and approved the submitted version.
Conflict of interests/ Competing interests
NGK authored a patent involving the use of lubricin for diagnostics assigned to the company Lynxon AB, and NGK holds equity in Lynxon AB. GJ and TS authored patents related to rhPRG4. GJ and TS hold equity in Lubris BioPharma LLC. TS was also paid consultancy fees by Lubris BioPharma. The remaining authors declare no competing interests. The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Sofia Grindberg and Paula Therese Kelly Petterssons at Danderyd’s Hospital and Lotta Falkendahl at the University of Gothenburg are acknowledged for their assistance in collecting samples. Dr Shan Huang is acknowledged in guidance in validation of the FILA. Lubris BioPharma is acknowledged for the generous gift of rhPRG4. We also gratefully acknowledge financial contribution to NGK from the Swedish state under the agreement between the Swedish government and the county council, the ALF-agreement (ALFGBG-722391), the Swedish Research Council (621-2013-5895), Petrus and Augusta Hedlund’s foundation (M-2016-0353), and AA insurance research fund (dnr 150150). GJ acknowledge funding from the National Institute of Health R01AR067748. ARA gratefully acknowledge OsloMet for providing a stipend to perform this project. Sponsors had no role in the study design, collection, analysis, and interpretation of data, nor were they involved in the writing of the manuscript and the decision to submit it for publication.