Proteomic analysis of urine reveals biomarkers for the diagnosis and phenotyping of abdominal-type Henoch-Schonlein purpura.

BACKGROUND: Abdominal-type Henoch-Schonlein purpura (HSP) is a common refractory disease in children. Currently, no specific diagnostic biomarker is available for HSP.
METHODS: Children with abdominal type HSP were first diagnosed with three syndromes using Chinese traditional medicine. The urinary proteomes among the three syndromes of patients with abdominal type HSP and healthy controls were compared using two label-free proteomics quantifications, including data-dependent acquisition and data-independent acquisition.
RESULTS: For the comparison between patients with abdominal type HSP and healthy children, a total of 75 differential urinary proteins were identified by determining the overlap of the two experiments. The ingenuity pathway analysis (IPA) analysis showed that these differential proteins were correlated with the pathogenesis of abdominal type HSP. Of these, 37 proteins were distributed in 13 solid tissues as tissue-enriched proteins. Monitoring changes in these proteins might help us detect uncommon clinical manifestations of HSP. Patients with abdominal type HSP can be further distinguished into three syndromes based on the urine proteome. Finally, a panel of six urinary proteins (P25774, P09417, Q7Z5L0, P60900, P14550 and P09668) was constructed for both the diagnosis and phenotyping of abdominal type HSP.
CONCLUSIONS: Urinary protein biomarkers for the diagnosis and phenotyping of abdominal type HSP were identified, which will contribute to the personalized treatment of patients with abdominal type HSP. 2021 Translational Pediatrics. All rights reserved.

Introduction

Henoch-Schonlein purpura (HSP) is the most common systemic vasculitis in children, with an incidence of 10–20 cases per 100,000 per year (1). The dominant clinical features include cutaneous purpura, arthritis, abdominal pain, gastrointestinal bleeding, and nephritis (2). Three-quarters of children with HSP develop abdominal symptoms, including diarrhea, abdominal pain, vomiting, gastrointestinal bleeding and hematochezia, and even intussusception and intestinal perforation (3,4). These patients can be classified as abdominal-type HSP. Abdominal type HSP is an acute and severe disease in pediatric patients, and severe cases can be secondary to intestinal obstruction and intussusception.

Currently, there are no specific diagnostic laboratory abnormalities for HSP. The serum IgA concentration is increased in 50% of patients, and serum C3 or C4 concentrations is low in a few patients. Other laboratory studies are useful only to exclude other conditions that may resemble HSP (5-8). The diagnosis of HSP is mainly based on the clinical manifestation. However, the abdominal symptoms often precede purpura in children with abdominal type HSP, and the clinical manifestation lacks specificity, which means that abdominal type HSP cannot be diagnosed early and is easily misdiagnosed. Diagnostic biomarkers for abdominal type HSP urgently need to be developed.

Information about the effective treatment for HSP is limited. Anti-infective, anti-allergic and symptomatic treatment are now dominant. Corticosteroid therapy has recently emerged as an effective treatment for abdominal pain (9). However, there is no evidence that corticosteroid therapy is effective in treating purpura, shortening the duration of the disease, or preventing recurrence (5). In addition, long-term use will lead to more serious adverse reactions. To reduce the incidence of adverse reactions and improve efficacy, corticosteroid therapy is often combined with Chinese herbs to treat abdominal type HSP in China.

Several studies have shown that corticosteroid therapy combined with Chinese herbs is more effective and safer than corticosteroid therapy alone (10,11). However, it is interesting that corticosteroid therapy is not combined with the same Chinese herbs for all patients with abdominal type HSP. In China, patients with abdominal type HSP will be divided into different syndromes according their clinical symptoms, such as the location and appearance of purpura, fever, cough, changes in stool characteristics, changes in tongue coating and pulse manifestations. Based on these different manifestations, three subtypes of HSP can be distinguished by expert doctors and then patients will be received significantly different drug prescriptions.

This inspired us to further divide abdominal type HSP into different types based on molecular phenotype. To achieve this aim, we adopted systems biology research methods. Urine is one of the main excreted human body fluids. Due to the lack of homeostatic mechanisms, it exhibits substantial changes in the body. Therefore, we hypothesized that the differences among different syndromes of patients with abdominal type HSP would be reflected in the urine. Proteins are one of the main components of urine. The urinary proteome is mainly composed of plasma proteins that pass through the glomerular barrier and proteins shed by cells within the urogenital system (12). It can reflect the changes in both systemic and local conditions of the body. With the development of MS techniques, more than 6,000 proteins have been identified in the normal human urinary proteome (13). Therefore, biomarker discovery through urine proteomics has become a very popular field and is also of great significance for early diagnosis, for clarifying the classification problem, and for better understanding the pathogenesis of abdominal type HSP.

In this study, urine samples were collected from children with abdominal type HSP with different syndromes under the guidance of traditional Chinese medicine. The urine proteomes were compared with the healthy controls by quantitative proteomics techniques to investigate whether patients with abdominal type HSP or the syndromes can be distinguished based on the urine proteome profile. Some differential urinary proteins with the potential for diagnosis and phenotyping of abdominal type HSP were identified, which might help guide drug therapy for abdominal type HSP in children.

We present the following article in accordance with the MDAR reporting checklist (available at http://dx.doi.org/10.21037/tp-20-317).

Methods

Participants and sample collection

Patients with abdominal type HSP and healthy controls were recruited from Beijing Children’s Hospital. These children were all aged from 1 to 10 years. Abdominal type HSP was diagnosed using the international criteria for HSP with abdominal symptoms. Children with abdominal type HSP were further diagnosed with three syndromes [wind-heat (WH) syndrome, damp-poison (DP) syndrome, and spleen-deficiency (SD) syndrome] based on traditional Chinese medicine by experienced Chinese medicine doctors (14). Patients with hematuria, severe proteinuria, and other diseases were excluded. Urine samples were collected from patients with abdominal type HSP and healthy controls. This study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by the local ethics committee of Beijing Children’s Hospital of Capital Medical University (No. 2016-91) and informed consent was taken from all the patients.

Classification standard of abdominal type HSP syndromes based on traditional Chinese medicine

WH syndrome

Patients have fever, slight chills, cough, and loss of appetite. Purpura occurs in the lower body, especially the lower limbs and buttocks. They are often symmetrical, bright red, present as papules or erythema, and vary in size and shape. They can be fused into a patch or itchy. Facial swelling, abdominal pain and bloody stool can also be observed. A red tongue with a thin greasy coating could be observed, and a superficial and rapid pulse may be felt (14).

DP syndrome

Skin of the lower extremities is full of purple freckle spots. Some symptoms can be observed, such as abdominal pain, bad breath, poor appetite, gastric fullness, bleeding gums, loose stool that is dark red or brown-purple or with roundworms. A red tongue with yellow coating can be observed, and a slippery pulse may be felt (14).

SD syndrome

The purpura is often pale red or recurrent. Patients have some symptoms, such as weight loss, sallow complexion, lassitude, general weakness, poor appetite, spontaneous sweating, low urine output, loose stool, abdominal pain, or even edema of the whole body or lower extremities. A pale tongue with a thin white coating could be observed, and a thready weak pulse or deep weak pulse may be felt (14).

Urine sample preparation

Midstream first morning urine was collected from these subjects with no serial measurements. All the urine samples of patients were collected when they first visit the hospital before treatment. After collection, the urine was immediately centrifuged at 5,000 g for 30 min at 4 °C to remove the cell debris. Protein from individual urine samples (30 mL) was precipitated with ethanol (90 mL) overnight at −20 °C followed by centrifugation at 10,000 g for 30 min to remove the supernatant. The precipitate was then resuspended in lysis buffer (8 M urea, 2 M thiourea, 50 mM Tris and 25 mM DTT). The protein amount of each sample was measured using the Bradford method. Then, 100 µg of protein from each individual sample was prepared using the filter-aided sample preparation method (FASP) for tryptic digestion (15). The peptide concentration was measured using the BCA method. Finally, 1 µg of each peptide sample was used for the subsequent liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis.

LC-MS/MS

Digested peptides were dissolved in 0.1% formic acid and loaded onto a trap column (75 µm × 2 cm, 3 µm, C18, 100 Å). The elute was separated on a reverse-phase C18 analytical column (50 µm × 250 mm, 2 µm, C18, 100 Å) using the EASY-nLC 1200 HPLC system (Thermo Fisher Scientific, Waltham, MA, USA). Then, the peptides were analyzed with an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA).

In the data-dependent acquisition (DDA) experiment, the elution was set to 5–30% buffer B (0.1% formic acid in 99.9% acetonitrile; flow rate of 0.3 μL/min) for 90 min. The MS parameters were set as follows: the full MS scan was acquired from 150 to 2,000 m/z with a resolution of 120,000, the MS/MS scan was acquired from 110 to 2,000 m/z with a resolution of 30,000 in Orbitrap, MS data were acquired in high-sensitivity mode, and 30% higher-energy collision dissociation (HCD) collision energy and charge-state screening (+2 to +7) were used.

In the data-independent acquisition (DIA) experiment, the LC settings were the same as in the DDA experiment to maintain retention time stability. The calibration kit (iRT kit from Biognosys, Switzerland) reagent was spiked at a concentration of 1:20 v/v in all samples. The iRT standard peptides were spiked into the urinary peptides for spectral library generation, and the following single DIA runs were used to calibrate the retention time of extracted peptide peaks. First, to generate the spectral library, all 39 samples were pooled as a mixture and fractionated by a high-pH reversed-phase peptide fractionation kit with 10 fractions. These ten fractions from the spin column were analyzed with a mass spectrometer in DDA mode. Then, the mass spectrometer was changed to the DIA mode. A variable isolation window DIA method with 45 windows was developed. The full scan was set at a resolution of 60,000 over a m/z range of 350 to 1,400, followed by DIA scans with a resolution of 30,000. The HCD collision energy was set to 30%. The maximal injection time was 50 ms for the full scan and DIA scans. The precursor ion number was equalized in each isolation window according to the precursor m/z distribution of the pooled sample. During the DIA analysis, a pooled sample was inserted as the quality control (QC).

Data processing

For the DDA experiment, the raw MS files were searched against the Swiss-Prot Human database using Mascot software (version 2.6.1, Matrix Science, London, UK). The parent ion tolerance was set at 10 ppm, and the fragment ion mass tolerance was set to 0.05 Da. A maximum of two missed cleavage sites in the trypsin digestion was allowed. Carbamidomethylation of cysteines was set as a fixed modification, and the oxidation of methionine was considered a variable modification. Label-free quantitation of the proteomic data was performed using Scaffold software (version 4.8.4, Proteome Software Inc., Portland, OR, USA). Peptide identifications were accepted with ≥2 unique peptides at a 1.0% false discovery rate (FDR) using the Scaffold Local FDR algorithm. Protein identifications were also accepted at an FDR of less than 1.0%. Spectral counts in different samples were compared after normalization with the total spectra according to previously described procedures (16).

For the DIA experiment, first, for library generation, the raw DDA MS files from the ten fractions were processed using Proteome Discoverer (Version 2.1, Thermo, USA) software with SEQUEST HT against the SwissProt database (Homo sapiens, year 201803) appended with the iRT fusion protein sequence. The search allowed two missed cleavage sites in trypsin digestion. Carbamidomethyl (C) was specified as the fixed modification. Oxidation (M) was specified as the variable modification. The parent ion mass tolerances were set to 10 ppm, and fragment ion mass tolerance was set to 0.02 Da. The Q value (FDR) cutoff at the precursor and protein levels was 1%. Finally, a total of 3,350 protein groups, 19,684 peptide groups and 61,401 peptide-spectrum matches (PSMs) were identified. Then, the results were imported to Spectronaut Pulsar X (Biognosys AG, Switzerland) software to generate the spectra library (17).

Then, all the raw DIA files were imported into Spectronaut Pulsar X as default settings. The optimal XIC extraction window was determined according to the iRT calibration strategy. The mass tolerance strategy was set to dynamic based on the extensive mass calibration. The cross-run normalization was set to the local normalization based on the local regression (18) Protein inference was performed using the implemented IDPicker algorithm (19). The Q value was set to 0.01 for data filtering. The sum peak areas of the respective fragment ions in MS2 were used for the quantification of peptide intensities. For later analysis, we required the proteins to be removed if the coefficient of variation (CV) in mixed samples were more than 30% or the missing values were present in more than 12 samples. Then, the missing values of protein abundance were filled with the K-nearest neighbors (KNN) method.

Statistical analysis

Average normalized spectral counts were used for the following DDA statistical analysis, and the normalized peak areas were used for DIA statistical analysis. Data were analyzed using SPSS statistical software (version 22.0). For comparisons between two groups, unpaired t-tests were used. For comparisons between more than two groups, analysis of variance (ANOVA) followed by a least significant difference (LSD) post-hoc test was used. Differences were considered significant when P<0.05.

To identify potential urine proteins to differentiate patients with abdominal type HSP from healthy controls, differential proteins in both DDA and DIA experiments were identified according to the following criteria: fold change ≥1.5 or ≤0.67. The above two criteria must be met at the same time. To identify potential urine proteins to differentiate the three disease syndromes, differential urinary proteins in patients with abdominal type HSP were further selected based on the following criteria: (I) ANOVA P<0.05 for the three disease syndromes; and (II) P value <0.05, change ratio >1.5 in any of these three disease syndromes when compared to each other.

Then, principal component analysis (PCA) and orthogonal projection to latent structures discriminant analysis (O-PLS-DA) were conducted in SIMCA 14.1 (Umetrics, Umea, Sweden) based on the selected biomarkers above. The biomarkers with the best predictive value to discriminate between the patients with abdominal type HSP and controls or between the abdominal type HSP with three syndromes were determined using area under the receiver operating characteristic (ROC) curve (AUC) analysis followed by calculations of sensitivity and specificity.

Biomarker annotation and functional analysis

The information for the urinary disease biomarkers was acquired from the Urinary Protein Biomarker Database (UPBD) (http://upbd.bmicc.cn/biomarker/web/indexdb). The information for the tissue specificity of the proteins was acquired from The Tissue Atlas (https://www.proteinatlas.org/humanproteome/tissue), which contains information regarding the expression profiles of human genes at both the mRNA and protein levels. The classification and criteria are as follows. Tissue enriched: at least four-fold higher mRNA level in a particular tissue compared to any other tissues. Group enriched: at least four-fold higher average mRNA level in a group of 2–5 tissues compared to any other tissue. Tissue enhanced: at least four-fold higher mRNA level in a particular tissue compared to the average level in all other tissues. Ingenuity pathway analysis (IPA) was used to determine the functional characterization and potential functional interaction for the differential proteins.

Results

Clinical characteristics of patients with abdominal type HSP and study design

A total of 20 subjects were included in the DDA experiment, including 5 for each symptom group and 5 for the healthy control group. A total of 39 subjects were included in the DIA experiment, including 11 children with WH syndrome, 10 children with DP syndrome, 9 children with SD syndrome and 9 healthy controls. The clinical characteristics of the patients are shown in and table online: https://cdn.amegroups.cn/static/public/tp-20-317-1.pdf. According to , based on traditional biochemical indicators, the three syndromes of patients with abdominal type HSP could not be distinguished.

Table 1

Clinical characteristics of subjects in this study

Variable	DDA method (n=20)						DIA method (n=39)
	HC	WH	DP	SD	P value‡		HC	WH	DP	SD	P value‡
Number	5	5	5	5	–		9	11	10	9	–
Sex (M/F)	2/3	3/2	4/1	3/2	–		6/3	4/7	8/2	6/3	–
Age	7.4±2.6	8.8±1.2	7.8±1.1	6.4±0.5	–		7.4±3.2	8.6±0.7	7.9±0.7	7.6±0.8	–
WBC (109/L)	4–10†	7.8±0.9	10.3±2.2	8.7±1.9	0.782		4–10†	8.39±0.68	8.80±0.90	9.10±1.20	0.927
N%	35–65†	68.2±3.3	72.9±6.4	70.2±3.1	0.577		35–65†	71.6±2.5	70.8±3.4	69.4±3.0	0.937
PLT (109/L)	100–400†	392.2±51.5	467.2±55.2	403.0±33.5	0.490		100–400†	406.9±26.7	464.9±34.5	417.4±28.5	0.527
FIB (g/L)	2–4†	5.5±1.9	2.8±0.2	2.6±0.4	0.079		2–4†	4.3±0.9	3.3±0.3	3.0±0.3	0.519
D-D (mg/L)	0–0.243†	0.92±0.48	1.66±0.97	0.59±0.22	0.976		0–0.243†	1.32±0.36	1.11±0.36	0.48±0.12	0.424
IgA (g/L)	0.6–2.2†	2.76±0.25	2.46±0.34	1.55±0.19*	0.019		0.6–2.2†	2.51±0.16	2.77±0.37	2.27±0.43	0.457

Values are given as the mean ± SD. †, the normal reference values of healthy children; ‡, comparison of the three syndromes was performed using the Kruskal-Wallis test, followed by post hoc multiple comparisons using the corrected Dunn’s test; *, P value <0.05 between WH and SD. HC, healthy control; WH, wind-heat syndrome; DP, damp-poison syndrome; SD, spleen-deficiency syndrome; WBC, white blood cell; N%, neutrophilic granulocyte percentage; PLT, blood platelet; FIB, fibrinogen; D-D, D-Dimer; IgA, immunoglobin A.

In this study, the urinary proteomes among three syndromes of patients with abdominal type HSP and healthy controls were compared. The flow chart of the study design is shown in . First, an experiment was initiated in which label-free proteomics quantifications with DDA were used, and 20 urine samples were included (five samples in each group). These subjects were matched by age and sex. Then, a larger sample size experiment was performed that used label-free proteomics quantifications with DIA and included a total of 39 urine samples, that is, 11 patients with WH syndrome in the WH group, 10 patients with DP syndrome in the DP group, 9 patients with SD syndrome in the SD group and 9 healthy children in the control group. The consistent comparison results from both experiments were selected for subsequent biomarker discovery analysis.

Figure 1

Flow chart of the study design. DDA, data-dependent acquisition; DIA, data-independent acquisition; PCA, principal component analysis; O-PLS-DA, orthogonal projection to latent structures discriminant analysis; ROC, receiver operating characteristic.

Identification of the potential diagnostic biomarkers for abdominal type HSP

To identify the potential diagnostic biomarkers for abdominal type HSP, the urinary proteome of abdominal type HSP children was compared with that of healthy children. In the DDA proteomics quantification experiment, a total of 1,292 urinary protein groups were identified in 20 urine samples. Finally, a total of 174 differential urine proteins were identified between abdominal type HSP and healthy children with a ratio >1.5 and P value <0.05 (table online: https://cdn.amegroups.cn/static/public/tp-20-317-2.pdf). Among these differential proteins, 84 urine proteins were upregulated, and 90 urine proteins were downregulated.

In the DIA proteomics quantification experiment, a total of 2,366 urinary proteins were identified in 39 urine samples. Specifically, 1,627 proteins were accepted when the CV values of protein abundance in QC samples and missing values in all samples were both less than 30%. Finally, 340 urinary proteins changed significantly in patients with abdominal type HSP compared with healthy controls (table online: https://cdn.amegroups.cn/static/public/tp-20-317-3.pdf). Among these differential proteins, 154 urine proteins were upregulated, and 186 urine proteins were downregulated.

The potential diagnostic biomarkers for abdominal type HSP were identified by comparing the differential proteins identified by the DDA experiment and DIA experiment. Finally, a total of 75 differentially expressed proteins, including 46 commonly upregulated proteins and 29 commonly downregulated proteins, were identified (table online: https://cdn.amegroups.cn/static/public/tp-20-317-4.pdf). These proteins have the potential to differentiate patients with abdominal type HSP from healthy children.

Annotation and functional analysis of the potential diagnostic biomarkers for abdominal type HSP

Whether the potential diagnostic biomarkers of abdominal type HSP found in this study have ever been reported as biomarkers for other diseases was investigated. The UPBD comprehensively collects information on urinary biomarkers or biomarker candidates from published literature. The 75 proteins were searched from the UPBD one by one. Ten proteins have been reported as candidate biomarkers for twenty-one diseases ().

Table 2

Ten potential biomarkers for abdominal type HSP were reported as candidate biomarkers for other diseases

Protein ID	Disease	Pmid	The potential reason why the abdominal type HSP shared biomarkers
P00450	Other diseases of intestines (K55–K64)	21915437	Intestinal involvement is the clinical characteristic of the patients included in this study. The most common complications of abdominal type HSP include intussusception, perforation, obstruction or gastrointestinal hemorrhage
P01019	Henoch-Schonlein purpura nephritis	21854508	Henoch-Schonlein purpura, also referred to as IgA vasculitis (IgAV), is characterized by immunoglobulin A1 (IgA1)-dominant immune deposits (20). HSP frequently affects kidneys, which is called HSP nephritis. It seems to share pathological mechanisms with IgA nephropathy (21). The renal involvement of HSP might share some biomarkers with other kidney diseases such as acute kidney failure, dent disease, kidney transplant rejection, and renal cell carcinoma. In addition, several biomarkers of other disease were shared with biomarkers of IgA glomerulonephritis, which is associated with pathological progression of HSP
P01019	IgA glomerulonephritis	21366514; 25523477
P00450	IgA glomerulonephritis	21595033
P04217	IgA glomerulonephritis	25957429
P01019	Chronic kidney disease	24065527; 24664631
P04217	Nephrotic syndrome	21591266
P07686	Proteinuria	21265931
P07686	Acute kidney failure	12584277; 17267747
P01019	Acute kidney failure	27538426
P07686	Nephropathy induced by other drugs, medicaments and biological substances	15967208
P01019	Dent disease	15140760
P07195	Dent disease	15140760
P07737	Dent disease	15140760
P01011	Kidney transplant rejection	17331118
P39059	Complications of kidney transplant	22253069
P00450	Renal cell carcinoma	23511837
P01019	Hydronephrosis	23772991	Ureteral obstruction and ureteritis are the main manifestations of IgAV involving the ureter, which will cause hydronephrosis. The involvement most commonly occurred at the uretero-pelvic junction (22,23)
P01019	Congenital occlusion of ureteropelvic junction	20639044
P04406	Congenital occlusion of ureteropelvic junction	20639044
P31151	Congenital occlusion of ureteropelvic junction	20639044
P00450	Diabetic nephropathy	17327332;26608305	HSP is a small vessel vasculitis. Microvasculitis is also one of major complications of diabetes (24). Vascular damage results in weakening of walls of the blood vessels and leakage of protein into the tissues, and an inflammatory response might happen in the pathological process of both diseases (24,25). In addition, many diabetes biomarkers were found to be shared with biomarkers of kidney diseases, which can also be involved in the progression of HSP
P01011	Diabetic nephropathy	17327332
P04217	Diabetic nephropathy	17327332
P00450	Type 2 diabetes mellitus	15111541
P00450	Type 2 diabetes mellitus with diabetic nephropathy	22536212
P07686	Type 2 diabetes mellitus	15111541
P01019	Type 2 diabetes mellitus	26617876
P04217	Type 1 diabetes mellitus	22678621
P01019	Cardiorenal syndrome	25722365	Cardiac and renal involvement can be mutually induced. In addition to kidney involvement, cardiac involvement may be common in patients with HSP. IgA and C3 deposited in the walls of intramyocardial vessel might account for its pathogenesis (26)In addition, P01019 is shared with biomarkers of Henoch-Schonlein purpura nephritis, which is associated with pathological progression of HSP
P01011	Acute appendicitis	19556024	Patients with HSP can also develop acute appendicitis (27). In addition, P01011 is shared with biomarkers of kidney involvement disease
P07737	Malignant neoplasm of bladder	22065568	HSP may affect the bladder, which leads to bladder wall hematomas and urinary retention (28). In addition, many bladder biomarkers were found to be shared with biomarkers of kidney diseases, which can also be involved in the progression of HSP
P04217	Malignant neoplasm of bladder	17518487
P39059	Malignant neoplasm of bladder	23389364
P07195	Malignant neoplasm of bladder	21496341

For speculating target organ damage, the expression levels of protein biomarkers in different tissues are very useful. The Tissue Atlas contains information regarding the expression data of human genes from 44 normal human tissue types at both the mRNA and protein levels. The 75 potential biomarkers were searched from the Tissue Atlas, and their expression-enriched or enhanced tissue types were acquired. A total of 24 tissues were identified (). The biomarker expression of enriched or enhanced tissue may be an organ that abdominal type HSP can affect. We speculated that the more biomarkers there were, the greater the possibility.

Figure 2

Annotation and functional characterization of differential proteins. (A) The distribution of the tissue-enriched proteins. (B,C,D) Disease and biofunctions analysis, including immune and inflammatory response class (B), lipid and carbohydrate metabolism class (C), and organismal injury and abnormalities class (D). (E,F) Differential proteins involved in two networks. Red: upregulated proteins; Green: downregulated proteins.

Then, functional characterization of the differentially expressed proteins was performed using IPA software. The IPA analysis revealed that the differential proteins were linked to three classes of disease and biofunctions, including the immune and inflammatory response class (), lipid and carbohydrate metabolism class (), and organismal injury and abnormalities class (). In addition, the differentially expressed proteins formed two tightly connected networks, which were associated with cardiovascular disease: organismal injury () and abnormal humoral immune response ().

Investigation of whether abdominal type HSP can be distinguished into different syndromes based on the urine proteome

In this study, patients with abdominal type HSP were divided into three different syndromes guided by traditional Chinese medicine: WH syndrome, DP syndrome and SD syndrome). The differential proteins among the three disease syndromes were investigated based on ANOVA.

Differential urinary proteins were selected based on the following criteria: (I) ANOVA P<0.05 for the three disease syndromes; and (II) P value <0.05, change ratio >1.5 in any of these three disease syndromes when compared to each other (post-hoc analysis). Finally, 91 and 73 differential proteins were acquired based on DDA and DIA quantification experimental data, respectively (tables online: https://cdn.amegroups.cn/static/public/tp-20-317-5.pdf). PCA and O-PLS-DA were performed based on differential urinary proteins. The results are shown in , indicating that the abdominal type HSP can be distinguished as three different syndromes based on urinary proteins.

Figure 3

PCA and O-PLS-DA analysis of differential proteins among three disease syndromes of abdominal type HSP. (A) PCA analysis with 91 differential proteins identified by DDA quantification; (B) PCA analysis with 73 differential proteins identified by DIA quantification; (C) O-PLS-DA analysis with 91 differential proteins identified by DDA quantification; (D) O-PLS-DA analysis with 73 differential proteins identified by DIA quantification. Green represents wind-hot syndrome, blue represents damp-poison syndrome, and red represents spleen-deficiency syndrome. WH, wind-heat; DP, damp-poison; SD, spleen-deficiency; PCA, principal component analysis; O-PLS-DA, orthogonal projection to latent structures discriminant analysis; HSP, Henoch-Schonlein purpura; DDA, data-dependent acquisition; DIA, data-independent acquisition.

Construction of the urinary protein biomarker panel for the diagnosis and phenotyping of abdominal type HSP

Only proteins with consistent changes between the DDA experiment and the DIA experiment were selected to determine the phenotyping and diagnosis urinary protein biomarkers. When the abdominal type HSP patients were compared with controls, we found 75 proteins with consistent changes between the two experiments. When the three syndromes were compared, a total of 21 proteins were found to have consistent changes by comparing the results of the DDA experiment and the DIA experiment, including two proteins for the comparison of WH and DP syndrome, three proteins for the comparison of WH and SD syndrome, and sixteen proteins for the comparison of DP and SD syndrome (table online: https://cdn.amegroups.cn/static/public/tp-20-317-6.pdf).

To evaluate the diagnostic performance of these candidate biomarkers, which were differentially expressed among the three disease syndromes, ROC curves were plotted. As shown in , a panel of three urinary proteins (P25774 + P09417 + Q7Z5L0) showed excellent sensitivity and specificity for the diagnosis of abdominal type HSP, with an AUC value of 0.966.

Figure 4

Diagnostic performance of biomarker panels for distinguishing abdominal type HSP and healthy controls and distinguishing three disease syndromes. The x-axis represents the diagnostic sensitivity of the biomarker panel, and the y-axis represents its diagnostic specificity. (A) panel for distinguishing abdominal type HSP from healthy controls; (B) panel for distinguishing wind-heat syndrome and damp-poison syndrome; (C) panel for distinguishing wind-heat syndrome and spleen-deficiency syndrome; (D) panel for distinguishing damp-poison syndrome and spleen-deficiency syndrome. P25774: cathepsin S; P09417: dihydropteridine reductase; Q7Z5L0: vitelline membrane outer layer protein 1 homolog; P14550: aldo-keto reductase family 1 member A1; P60900: proteasome subunit alpha type-6; P09668: pro-cathepsin H. CI, confidence interval; WH, wind-heat; DP, damp-poison; SD, spleen-deficiency; HSP, Henoch-Schonlein purpura.

Several urine proteins could be used for distinguishing the three disease syndromes. As shown in , a protein panel (P14550 + P60900) could differentiate WH and DP syndromes with an AUC value of 0.96. Another protein panel of P25774 and P09668 could differentiate WH and SD syndromes with an AUC value of 0.83 (). One urine protein (P09417) could differentiate DP syndrome and SD syndrome with an AUC value of 0.956 ().

These results indicated that a minimum of six urinary proteins (P25774 + P09417 + Q7Z5L0 + P60900 + P14550 + P09668) is required to form a urinary protein panel, which has the potential for the diagnosis of abdominal type HSP and the three syndromes of abdominal type HSP.

Discussion

In the present study, we collected urine samples from healthy children and abdominal type HSP patients with three syndromes according to the guidance of traditional Chinese medicine. The urine proteome was analyzed by two quantitative proteomics experiments. The consistent findings of the two experiments were presented in this study. First, the proteome of the abdominal type HSP patients was compared with that of the controls. By determining the overlap of the different proteins in the two experiments, 75 disease-related proteins were identified. The annotation and function of these 75 proteins were analyzed.

Abdominal type HSP is Henoch-Schonlein purpura, which involves the gastrointestinal tract. HSP is also referred to as IgA vasculitis (IgAV). It is characterized by immunoglobulin A1 (IgA1)-dominant immune deposits affecting small vessels and often involves the skin, gastrointestinal tract, joints, and kidneys (20). Due to the involvement of small blood vessels throughout the body, this disease can cause a variety of symptoms in different organs.

Ten of 75 proteins have been reported as candidate biomarkers for twenty-one diseases. Most proteins have been reported as candidate biomarkers for more than one disease. This suggests that similar pathological damage may exist among these diseases, and the differential proteins can reflect the changes caused by such pathological damage. In this case, we could not determine which disease was the primary cause of the pathological damage. However, according to the multisystem manifestations of HSP, we believe that most of the 21 diseases may have similar pathological damage. From the point of view of pathological damage, the reasoning for each protein being annotated to different disease biomarkers is shown in .

Based on the analysis in , we believe it is not surprising that abdominal type HSP shares biomarkers with most of the 21 diseases. In turn, this increases the reliability of these proteins as abdominal type HSP biomarkers. At the same time, many kidney diseases share the same biomarkers as abdominal type HSP, which suggests that in addition to presenting abdominal symptoms, abdominal type HSP patients may also have some degree of renal damage. It should be noted that the renal clinical indicators of these patients were normal. Therefore, these injuries should be minor, and urine biomarkers are very sensitive and can reflect renal changes in very early stages.

Currently, uncommon manifestations of HSP pose diagnostic difficulties. Some children may lack clinical symptoms in the early stage. Alternatively, the signs of uncommon manifestations may be nonspecific and masked by common complications. Therefore, the indicators of organ damage in the early stage are very useful for the diagnosis of uncommon clinical manifestations of HSP. Since urine biomarkers can reflect changes in organs in the early stage, the organ origin of the 75 differentially expressed proteins was analyzed. Because only a few proteins show strict tissue-specific expression in a single tissue or organ, we can only ascertain the tissues or organs that express the protein at a significantly elevated level. Finally, thirty-seven of the 75 proteins were found to be distributed in 13 solid tissues as tissue-enriched proteins. All 13 solid tissues can be involved in the progression of HSP. The analysis results are shown in . In the future, monitoring changes in these tissue-enriched proteins over a longer time period and in larger populations may help us identify biomarkers that can predict specific organ damage by HSP.

Table 3

Tissue-enriched protein biomarkers and tissue involvement evidence in the pathogenesis of abdominal type HSP

Differential proteins found in this study	Number of proteins	Enriched or enhanced tissue	Evidence of tissue involvement in the pathogenesis of abdominal type HSP
P01019; P61457; P00450; P01011; P09467; P04217; P52758; P17174; P09417; P34896; P09210; Q9Y2S2; Q5R3I4; P30046; P05062	15	Liver	Hepatobiliary involvement may happen in children with HSP (29). In addition, cases of HSP with primary biliary cirrhosis (PBC) have also been reported, and autoantigen-specific IgA antibody seems to play an important role in their pathogenesis (30)
P52758; P34896; O75309; P09210; O94760; Q9Y2S2; P05062; P07195; Q07075	9	Kidney	HSP frequently affects kidneys, which is called HSP nephritis (31). It seems to share pathological mechanisms with IgA nephropathy (21)
Q8TAB3; P09417; P21802; A6NL88; O94910; Q9HC56; Q9BXN2	7	Brain	Central nervous system (CNS) involvement may present as a complication of HSP, which include central nervous system (CNS) injury, cerebral vasculitis, cerebral hemorrhage, posterior reversible encephalopathy syndrome (PRES), seizures, and peripheral nervous system injury (32)
P05062; P29323; Q5R3I4; Q07075	4	Intestine	The small intestine is the most frequently involved site in the GI tract by HSP due to its predilection toward ischemic injury (25)
P04406; P21695; P17174	3	Skeletal muscle	HSP with muscle involvement, presenting as myositis (33)
P09467; Q01151; Q9BXN2	3	Lung	Pulmonary involvement may be a complication of HSP and diffuse alveolar hemorrhage (DAH) is the most frequent clinical presentation (34)
P17174; P07195	2	Heart muscle	Cardiac involvement may be common in patients with HSP. IgA and C3 deposited in the walls of intramyocardial vessel might account for its pathogenesis (26)
P09210	1	Pancreas	Acute pancreatitis may present as a complication of HSP (35)
P48745; Q8NHP8; P80370	3	Adrenal gland	The adrenal gland may be involved in the progression of HSP (36)
Q9H9P2; Q8IUL8	2	Testis	The testis may be involved in the progression of HSP (37)
Q13145; Q8IUL8	2	Ovary	The ovary may be involved in the progression of HSP (38)
P31151	1	Esophagus	The esophagus may be involved in the progression of HSP (39)
Q6FHJ7	1	Retina	Bilateral central retinal artery occlusion may be a complication of HSP (40)

HSP, Henoch-Schonlein purpura; GI, gastrointestinal.

The IPA analysis showed that these differential proteins were correlated with the pathogenesis of abdominal type HSP (). Macrophages, neutrophils, the immune system, and the inflammatory response were first shown to be linked with differentially expressed proteins. HSP is an immune-complex mediated vasculitis in childhood. Immune complex deposits result in necrosis of the wall of small- and medium-sized arteries with infiltration of tissue by neutrophils and deposition of nuclear fragments. This process is called leukocytoclastic vasculitis (LCV) (25). Neutrophils can contribute to vasculitis by increasing the release of neutrophil extracellular traps (41). Macrophages have also been found to play a decisive role in the pathogenesis of vasculitis (42). Lipid and carbohydrate metabolism were also shown to be linked with differentially expressed proteins, which might be because cellular metabolism is important in determining the immune responsiveness of lymphocytes (43). Finally, the IPA results also showed that these differentially expressed proteins can reflect multiorgan abnormalities, which was consistent with the results obtained from the tissue origin analysis of the 75 proteins.

Second, whether abdominal type HSP can be distinguished into different syndromes based on the urine proteome was investigated. PCA showed that although they could not be completely separated, there was a tendency to distinguish among the three syndromes either in DDA or DIA quantification experiments. O-PLS-DA showed that the abdominal type HSP can be classified into three syndromes based on the selected differential proteins. SD syndrome and DP syndrome were more easily separated. In the theory of traditional Chinese medicine, SD syndrome is a deficiency syndrome, while DP syndrome is an excessive syndrome. They represent opposite pathogeneses, so there will be significant differences in the proteome between the two groups.

Third, the urinary protein biomarker panel for the diagnosis and phenotyping of abdominal type HSP was constructed. As most biomarkers of HSP are not specific and originate from inflammation, immune response or vasculitis, we therefore constructed a biomarker panel to improve the diagnosis performance of HSP and facilitate its future clinical application. After using ROC analysis to evaluate the diagnostic effects of candidate biomarkers, we found that at least six differential proteins (P25774 + P09417 + Q7Z5L0 + P60900 + P14550 + P09668) are needed. After searching the UPBD (44), none of these six proteins in this biomarker panel had been reported as urinary disease biomarkers in previous studies.

This study is an attempt to investigate the markers of the molecular phenotype of the disease based on the urine proteome with the help of clues obtained from clinical treatment experience. The protein panel might be useful not only for the diagnosis of abdominal type HSP but also for distinguishing the three syndromes of abdominal type HSP. Furthermore, because the three disease syndromes respond differently to treatment, this biomarker panel could advantageously be used for patient stratification. Although the protein panel was constructed by the common proteins found in two urine proteome quantitative comparison experiments, it can still only be regarded as a potential biomarker, and large-scale clinical verification is needed in future experiments.

In China, to obtain better clinical outcomes, many diseases are divided into several syndromes based on traditional Chinese medicine and receive different treatments. Although use of traditional Chinese medicine theory has been controversial, all these results can be regarded as indicators when studying disease phenotypes based on modern multiomics strategies. Specifically, urine can accumulate systematic changes in the body; thus, urine proteomics has great clinical application prospects in the classification of syndromes of diseases, especially for verifying these clues. In the future, urine proteome biomarkers will be increasingly used in the classification of disease syndromes and will guide precise treatment.

The article’s supplementary files as