Patients
The study included a total of 15 fresh serum specimens, which were divided into three groups: SONFH, TONFH and CK group, with 5 independent individuals in each group (Table 1). All the female patients and healthy volunteers were derived from a clinical laboratory in the First Affiliated Hospital of Guangxi University of Chinese Medicine (Nanning, China) from April 2014 to June 2016. The systemic lupus erythematosus combined with SONFH and TONFH patients were selected according to the definite diagnosis by at least three orthopedists (Supplementary Table 1). Serum was distinguished from blood moderately drawn from bodies and then stored at -80 ℃ until detection.
Table 1
Sample information of experiments
Group
|
Number of females
|
Age (years, Median[range])
|
DDA
|
WB
|
DDA
|
WB
|
CK
|
5
|
8
|
43.0 (32 ~ 46)(28 ~ 46)
|
42.0
|
SONFH
|
5
|
10
|
44.6 (39–53)(23 ~ 65)
|
43.5
|
TONFH
|
5
|
5
|
45.4 (30 ~ 61)(30 ~ 61)
|
45.4
|
Protein extraction and digestion
For protein extraction, 100 µl of the serum samples (each for 5 adult females) were employed for high abundant protein removal based on the ProteoPrep Blue Albumin and IgG Depletion Kit13. The protein eluent was precipitated with cold acetone for 3 h at -20 ℃. After centrifugation at 4 ℃ at 12000 g for 10 min, the protein deposit was redissolved in lysis buffer (8 M urea, 100 mM triethylammonium bicarbonate (TEAB)). The protein concentration was determined by a Modified Bradford Protein Assay Kit according to the manufacturer’s instructions. For digestion, 100 µg protein of each sample was first reduced with 10 mM DTT at 37 ℃ for 60 min and then alkylated with 25 mM iodoacetamide (IAM) at room temperature for 45 min in darkness. The urea concentration of the protein sample was diluted to less than 2 M by adding 100 mM TEAB. The protein pool of each sample was digested with Sequencing Grade Modified Trypsin with the ratio of protein: trypsin = 50:1 mass ratio at 37 ℃ overnight and 100:1 for continuous digestion for 4 h.
Protein isobaric labelling and sample cleanup
After trypsin digestion, peptides were desalted by Strata X SPE column and vacuum-dried. Peptides were reconstituted in 100 µL 100 mM TEAB and processed according to the protocol of manufacturer’s 6-plex tandem mass tags (TMT) kit40. Next, one unit of TMT reagent was added to the peptide solution after thawed and dissolved in 41 µL isopropanol. The peptide mixtures were incubated for 1 h at room temperature, then pooled and dried by vacuum centrifugation. The dried and labeled peptides were reconstituted with HPLC solution A (2% ACN, pH 10) and then fractionated into fractions by high-pH reverse-phase HPLC using Waters Bridge Peptide BEH C18 (130 Å, 3.5 µm, 4.6×250 mm). Briefly, peptides were first separated into 60 fractions at a speed of 0.5 ml/min over 88 min with a gradient of 2–98% acetonitrile in pH 10. Until then, the peptides were combined into 20 fractions and dried by vacuum centrifugation. After which, the peptide fractions were desalted by Ziptip C18 according to manufacturer’s instructions and finally dried under vacuum before being kept at -20 ℃ until MS analyses were performed.
High-resolution LC-MS/MS analysis
The mass spectrometry survey was performed by Nano LC 1000 LC-MS/MS using a Proxeon EASY-nLC 1000 coupled to an LTQ-Orbitrap Elite. Trypsin digestion fractions were reconstituted in 0.1% formic acid (FA) and directly loaded onto a reversed-phase pre-column (Acclaim PepMap®100 C18, 3µm, 100Å, 75µm× 2cm) delivering at 5 µL/min in 100% solvent A (0.1 M acetic acid in water). After that, peptides eluted from the trap column were loaded onto a reversed-phase analytical column (Acclaim PepMap® RSLC C18, 2 µm, 100Å, 50 µm × 15 cm) which gradient was comprised of an increase from 10–35% solvent B (0.1% FA in 98% ACN) over 60 min, 35–50% in 10 min and climbing to 100% in 5 min at a constant flow rate of 250 nl/min on an EASY-nLC 1000 system. The eluent was sprayed via NSI source at the 1.8 kV electrospray voltage and then analyzed by MS/MS in LTQ-Orbitrap Elite. The mass spectrometer was operated in data-dependent mode, and automatically switches between MS and MS/MS. Full-scan MS spectra (from m/z 350 to 1800) were acquired in the Orbitrap with a resolution of 60,000. Ion fragments were detected in the Orbitrap at a resolution of 15,000. And the 20 most intense precursors were selected for subsequent decision tree-based ion trap higher energy collision-induced dissociation (HCD) fragmentation at the collision energy of 38% above a threshold ion count of 300 in the MS survey scan with 30.0s dynamic exclusion. The full width at half maximum (FHMW) at 400 m/z using an automatic gain control (AGC) setting of 1e6 ions. The fixed first mass was set as 100 m/z.
Data processing
The resulting MS/MS raw data were converted to mgf format profile with the software Proteome Discoverer (version 1.3, Thermo Scientific). Generated peak lists were searched against the Homo Sapiens database (Taxon identifier: 9606, including 154527 protein sequences) which consists of reviewed Swissprot database combined with unreviewed TrEMBL database. Mascot 2.3.0 is used to supplement frequently observed contaminants. Trypsin/P was chosen as the enzyme and two missed cleavages were allowed. Carbamidomethylation on cysteine was set as a fixed modification, and oxidation on methionine, acetylation in N-term were set as variable modifications. The peptide used in the search had a mass tolerance of 20 ppm, a production tolerance of 0.02 DA and a false discovery rate (FDR) of 1%. For quantitation, proteins were required to contain at least two quantitated unique peptides. Proteins with fold changes > 1.2 or < 0.83 and significance t-test p-value < 0.05 between two compared groups were considered to be differentially expressed.
Then proteins were classified by Gene Ontology annotation based on three categories: biological process, cellular component and molecular function. Gene Ontology (GO) annotation proteome was derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/). KEGG signaling pathway is part of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database41–43, which is a reference database for pathway mapping. KEGG Pathway analyses of identified proteins were extracted using the Search pathway tool in the KEGG Mapper platform (http://www.genome.jp/kegg/mapper.html). InterProScan, a sequence analysis application was also used for protein domain annotation based on the protein sequence alignment algorithm and the InterPro domain database44.
A two-tailed Fisher’s exact test was employed to test the GO, KEGG pathway and Domain enrichment of the differential expression protein for all identified proteins. Correction for multiple hypothesis testing was carried out under standard FDR control methods, and pathways with a corrected p-value < 0.05 were considered the most significant.
Different proteome clustering methods based on expression and function enrichment were used to explore potential relationships between different protein groups in special protein functions. Firstly, all the protein groups obtained after functional enrichment analysis along with their p values were collected. Secondly, these categories enriched were sorted in at least one protein group with a p-value < 0.05. This filtered p-value matrix was transformed by the function x = − log10 (p value). Thirdly, z-transformed applies to x values for each functional category and z scores were clustered by one-way hierarchical clustering (Euclidean distance, average linkage clustering). Finally, the expression pattern of the protein was completed using the timeclust function in the R package ‘TCseq’, and a heatmap was generated by the software Heml 1.0.3 45. We also present the STRING network of some known/novel protein-protein interactions as an evidence view by using the STRING Version 9.0 (Search Tool for the Retrieval of Interacting Genes/Proteins), a database of physical and functional interactions (http://string-db.org/) 46.
Western blotting
Equal amounts of proteins (30 µg) were separated by 12% SDS-PAGE and electrophoretically transferred to a polyvinylidene fluoride (PVDF) membrane (Millipore, MA, USA). Membranes were blocked for 4 h at 4 ℃ with 5% non-fat milk in 1 × TBST buffer and rinsed 3 times (10 min each) with 1 × TBST buffer. Next, the membranes were incubated with primary antibodies: Anti-C reactive protein (1:1000, Abcam, ab32412), Anti-alpha-1-acid glycoprotein 1 (1:1000, Abcam, ab134160), Anti-SAA1 (1:1000, Abcam, ab201660), Anti-dopamine beta-hydroxylase (1:2000, Abcam, ab96615) at 4 ℃ overnight and then washed with 1 × TBST buffer for 3 times. After that, the secondary antibodies was be in progress: Peroxidase AffiniPure Goat Anti-Mouse IgG (H + L) (1:20000, Jackson, 115-035-003) and Peroxidase AffiniPure Goat Anti-Rabbit IgG (H + L) (1:20000, Jackson, 111-035-003) for 4 h at 4 ℃. After washing 3 times, the membranes were developed using Western blotting ECL Substrate (Bio-Rad) detection methods (Solarbio, PE0010) and scanned with the ChemiDoc MP imaging system using Image Lab software (Bio-Rad).