Skip to main content

Discovery and validation of novel protein markers in mucosa of portal hypertensive gastropathy



Portal hypertension induced esophageal and gastric variceal bleeding is the main cause of death among patients of decompensated liver cirrhosis. Therefore, a standardized, biomarker-based test, to make an early-stage non-invasive risk assessment of portal hypertension, is highly desirable. However, no fit-for-purpose biomarkers have yet been identified.


We conducted a pilot study consisting of 5 portal hypertensive gastropathy (PHG) patients and 5 normal controls, sampling the gastric mucosa of normal controls and PHG patients before and after endoscopic cyanoacrylate injection, using label-free quantitative (LFQ) mass spectrometry, to identify potential biomarker candidates in gastric mucosa from PHG patients and normal controls. Then we further used parallel reaction monitoring (PRM) to verify the abundance of the targeted protein.


LFQ analyses identified 423 significantly differentially expressed proteins. 17 proteins that significantly elevated in the gastric mucosa of PHG patients were further validated using PRM.


This is the first application of an LFQ-PRM workflow to identify and validate PHG–specific biomarkers in patient gastric mucosa samples. Our findings lay the foundation for comprehending the molecular mechanisms of PHG pathogenesis, and provide potential applications for useful biomarkers in early diagnosis and treatment.

Trial registration and ethics approval: Trial registration was completed (ChiCTR2000029840) on February 25, 2020. Ethics Approvals were completed on July 17, 2017 (NYSZYYEC20180003) and February 15, 2020 (NYSZYYEC20200005).

Peer Review reports


Liver cirrhosis is the end stage of chronic liver disease characterized by fibrosis and structurally abnormal nodules. Decompensation stage cirrhosis manifests as hepatic dysfunction and portal hypertension (PH) [1]. The PH-induced esophagogastric variceal bleeding (EGVB) was the main cause of death in patients with cirrhosis [1]. Therefore, it is of great significance to understand the pathogenesis and clinical manifestations of EGVB and develop effective prevention and treatment strategies accordingly to reduce the occurrence of EGVB and its related mortality.

According to Portal Hypertensive Bleeding in Cirrhosis: Risk Stratification, Diagnosis, and Management: 2016 Practice Guidance by the American Association for the Study of Liver Diseases [2], platelet count, ultrasound, transient elastography (TE), computed tomography (CT), or magnetic resonance imaging (MRI) are suggested as non-invasive methods to evaluate the severity of cirrhosis and PH. It also proposed that patients with compensated cirrhosis with liver stiffness < 20 kPa (determined by TE) and a platelet count > 150,000/mm3 were very unlikely to have high-risk varices (< 5%), and endoscopy could be safely avoided in them [3]. Patients with cirrhosis who do not meet the above criteria should undergo endoscopic examination to assess the presence of varicose veins. However, the imaging examination methods mentioned above have higher requirements on operating physicians and are subjective and restrictive to a certain extent. Currently, there is still a lack of non-invasive indicators for PH severity and objective quantitative assessment. Accordingly, there remains a clinical need for new biomarkers that achieve quantified, inexpensive, and minimally invasive testing to assess the severity of PH.

For patients with gastroesophageal varices (GOV) or EGVB, endoscopic cyanoacrylate injection eradicates varices with assistance of preoperative portal vein CT angiography and endoscopic ultrasound. However, our follow-up after 1 week of operation found that portal hypertensive gastropathy (PHG) became more severe, presenting as exacerbation of redness, point bleeding or a fine reticular pattern separating areas of raised edematous mucosa, the "mosaic pattern", which may attributes to the development of gastric capillary ectasia and collateral veins after obstruction of the esophageal varices [4]. Gastroscopy can provide an accurate PHG diagnosis and classification [4, 5]. Therefore, we sampled the gastric mucosa of the patients with compensated cirrhosis before and after endoscopic cyanoacrylate injection, and attempted to discover differentially expressed proteins (DEPs) which may appropriate for non-invasive diagnosis of PH severity.

Several studies have sought to determine disease biomarkers in gastric mucosa, with a view to the development of clinically relevant risk stratification and diagnostic tests [6, 7]. However, these studies didn’t provide sufficient resolution to differentiate the stage of PHG. The need to develop new and more effective non-invasive tests is therefore still obvious. Here, we report a pilot study applying label-free quantitative (LFQ) mass spectrometry to discover potential biomarker candidates, followed by targeted parallel reaction monitoring (PRM) to verify the PHG-related protein changes between decompensated liver cirrhosis patients and normal controls. This is the first use of LFQ-PRM for PHG biomarker discovery, and establishes the methodology required for subsequent studies to evaluate specific disease biomarkers of decompensated liver cirrhosis.


Clinical sample information

In this study, we evaluated the application of label-free mass spectrometry and PRM to identify and validate promising biomarkers for assessing the severity of PH in decompensated liver cirrhosis patients. The overall strategy and simplified technology roadmap are shown in Fig. 1a. In total, 29 gastric mucosa samples (n = 9 Nor group, n = 10 Ph group, n = 10 Inter group) were collected from individual patients for this study. The baseline characteristics of each group are presented in Table 1. Compared with healthy volunteers, the gastric mucosa of liver cirrhosis patients presents reddening, edematous, or even point bleeding and the appearance of the gastric mucosa after endoscopic cyanoacrylate injection were more severe (Fig. 1b–d). To provide consistency at point-of-collection, gastric mucosa samples were collected immediately during gastroscopy and processed identically.

Fig. 1
figure 1

a Workflow describing LFQ and PRM analysis. Proteins from samples were digested by trypsin and analyzed by LS–MS/MS. The differentially expressed proteins were analyzed by database researches and validated by PRM. Nor. group = healthy individuals, Ph. group = PHG patients and Inter. group = PHG patients after endoscopic cyanoacrylate injection for 1 week. b The gastroscopic photographs show normal appearance of esophagus and fundus, c severe esophageal varices and mild PHG and d esophageal varices alleviation and PHG exacerbation after endoscopic cyanoacrylate injection

Table 1 Baseline characteristics of PHG patients compared with controls

Values reported as mean ± SEM. Classification of PHG: Grade 0: no PHG. Grade 1: mild reddening, mucosa was congestive but there was no mosaic pattern. Grade 2: severe redness and a fine reticular pattern separating areas of raised edematous mucosa, the "mosaic pattern," or a fine red speckling were present. Grade 3: Point bleeding was recognized in the state of grade 2.

Identification of DEPs by LS–MS/MS Analysis

Using the LFQ workflow we identified a total of 6443 proteins (57,277 peptides) in the discovery cohort samples, with quantified a total of 5595 proteins. We utilized the quantitative ratio value of relative expression of protein and the corresponding t-test p value to screen the DEPs of normal controls and PHG patients before or after operation. A total of 423 proteins were found to be significantly differentially expressed under the selection criteria, of which 226 between Inter. group and Nor. group (113 up-regulated and 113 down-regulated), 19 between Inter. group and Ph. group (4 up-regulated and 15 down-regulated), 178 were found between Ph. group and Nor. group (83 up-regulated and 95 down-regulated) (Fig. 2a and Additional file 1). Out of 423 DEPs, 101 DEPs were found intersected among three group by comparison, which is shown by Venn diagram (Fig. 2b). The volcano plots based on the 423 DEPs are shown in Fig. 2c. All up-regulated and down-regulated DEPs are presented by heatmap, and their cluster analysis are performed (Fig. 3).

Fig. 2
figure 2

DEPs distribution in different comparative analysis groups between Inter versus Nor, Ph versus Nor and Inter versus Ph. a The DEPs distribution of each group is presented by bar chart. b A Venn diagram showing the distribution of DEPs. c DEPs distribution exhibits by volcano plots, which had fold change > 1.5 or < 1/1.5 and p < 0.05

Fig. 3
figure 3

Heatmap and cluster analysis of 423 DEPs in Inter versus Nor, Ph versus Nor and Inter versus Ph. Each line in the heatmap represents the mean fold-change of up-regulation (red) or down-regulation (green) of a protein's expression levels

Bioinformatics analysis of differentially expressed proteins

Gene Ontology (GO) level 2 analysis was performed to explain relevant biological processes (BP), molecular functions (MF), and cellular components (CC) of these differentially express proteins (Fig. 4a).

Fig. 4
figure 4

Gene Ontology (GO) level 2 analysis and subcellular location of DEPs. a Functional classification of DEPs for relevant biological processes, molecular functions, and cellular components. b Subcellular location of DEPs

In the Inter/Nor comparison, the top three biological processes associated with the 226 DEPs were: cellular process (187 proteins), biological regulation (154 proteins) and metabolic process (117 proteins). The most relevant cellular component was the cell (197 proteins), followed by organelle (162 proteins), and membrane (107 proteins). For the top three molecular functions, GO terms relating to binding (141 proteins) were the most highly represented, followed by catalytic activity (68 proteins), and molecular function regulator (26 proteins).

In the Inter /Ph comparison, the top three biological processes associated with the 19 DEPs were: cellular process (15 proteins), cellular component organization or biogenesis (10 proteins), and biological regulation (10 proteins). The most relevant cellular component was the cell (16 proteins), followed by organelle (14 proteins), and membrane (9 proteins). For the top three molecular functions, GO terms relating to binding (10 proteins) were the most highly represented, followed by catalytic activity (5 proteins), and molecular function regulator (2 proteins).

In the Ph /Nor comparison, the top three biological processes associated with the 178 DEPs were: cellular process (129 proteins), biological regulation (120 proteins), and response to stimulus (94 proteins). The most relevant cellular component was the cell (146 proteins), followed by organelle (117 proteins), and membrane (79 proteins). For the top three molecular functions, GO terms relating to binding (99 proteins) were the most highly represented, followed by catalytic activity (50 proteins), and molecular function regulator (24 proteins).

To further determine the location of DEPs, subcellular location analysis was performed (Fig. 4b). The top three subcellular locations in Inter/Nor groups were extracellular (75 proteins, 33.19%), nucleus (47 proteins, 20.8%) and cytoplasm (40 proteins, 17.7%). The top three subcellular locations in Inter/Ph groups were mitochondria (5 proteins, 26.32%), nucleus (5 proteins, 26.32%) and cytoplasm (4 proteins, 21.05%). The top three subcellular locations in Ph/Nor groups were extracellular (69 proteins, 38.76%), nucleus (32 proteins, 17.98%) and cytoplasm (24 proteins, 13.48%).

Pathway enrichment analysis of differentially expressed proteins

Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used as annotations for protein metabolic pathways. As it’s shown in Fig. 5, the gastric mucosa proteins up-regulated in PHG were involved in several pathways. In Ph/Nor comparison, hypertrophic cardiomyopathy (map05410), vascular smooth muscle contraction (map04270) and cardiac muscle contraction (map04260) were the top 3 enriched pathways. In Inter /Nor comparison, focal adhesion (map04510), protein digestion and absorption (map04974), and extracellular matrix (ECM)-receptor interaction (map04512) were the top 3 enriched pathways. It’s worth noting that focal adhesion (map04510) and vascular smooth muscle contraction (map04270) were both enriched in Ph/Nor groups (p = 0.03, p = 0.0001989, respectively) and Inter/Nor groups (p = 1.50 × 10–6, p = 0.01, respectively). However, no pathway was enriched in Inter/Ph groups.

Fig. 5
figure 5

KEGG pathway analysis and heatmap [50,51,52] of DEPs in Inter versus Nor, Ph versus Nor and Inter versus Ph

Protein interaction analysis (PPI)

We used the STRING database to build a regulatory network from DEPs participated in focal adhesion and vascular smooth contraction pathway in Inter/Nor group and Ph/Nor group to visualize the predicted interactions between DEPs in PHG (Fig. 6). In Inter/Nor group, the DEPs involved in focal adhesion include COL6A1, COL6A2, COL6A3, CAV2, TNC, ITGA1, MYLK, VWF, ACTN1, FLNA and FLNC. The DEPs involved in vascular smooth contraction include MYLK, ACTA2, ACTG2 and CALD1. In Ph/Nor group, the DEPs involved in focal adhesion include PDGFRB, TNC, CAV2 and ITGA1. The DEPs involved in vascular smooth contraction include AGT, ACTA2, ACTG2, CALD1 and CACNA1S.

Fig. 6
figure 6

STRING protein–protein interaction (PPI) network of validated interactions between DEPs enriched in focal adhesion and vascular smooth muscle contraction. They were created using STRING v11.0 software together with a ‘Homo sapiens’ database

PRM analysis for protein quantification

18 up-regulated DEPs validated by LS–MS/MS analysis related to liver diseases and gastrointestinal system were chosen for further study, their expression in gastric mucosa was measured by PRM. Due to limitation of the proteins’ properties and abundance, 17 of them were successfully quantified and 1 was not. The protein quantification results are shown in Table 2. The heatmap and boxplot are shown in Fig. 7. In accordance with LFQ, these proteins were significantly up-regulated in Ph group and Inter group compared with Nor group, the same trend was also observed in Inter group compared with Ph group.

Table 2 17 differentially expressed proteins validated by parallel reaction monitoring quantification (PRM)
Fig. 7
figure 7

Heatmap and box plot of 17 DEPs in PHG patients before and after endoscopic cyanoacrylate injection compared with controls by PRM analysis. Each line in the heatmap represents the mean fold-change of upregulation (red) or downregulation (green) of a protein's expression levels


Currently, the severity of GOV caused by PH is critical for the prognosis of patients with decompensated liver cirrhosis. Most of the need lies in the detection of liver cirrhosis patients at risk of occurrence and the progression of PH, and early diagnosis and treatment may greatly contribute to enhancing the survival of such patients. Decompensated liver cirrhosis presents not only as GOV but also PHG. The severity of swelling and reddening of gastric mucosa [8], interstitial vascular ectasia, spindle cell proliferation and fibrohyalinosis [9] reflect the degree of PHG and are associated with severity of liver disease. Monitoring changes in the gastric mucosa proteome may have an advantage in the diagnosis of decompensated liver cirrhosis. However, targeted proteomics, which identifies promising biomarkers of disease activity and organ involvement, is rarely involved in non-invasive method for detection and diagnosis of liver cirrhosis induced PH. Our group devoted to the research of proteomics in liver disease [10, 11]. To the best of our knowledge, no proteomic analysis of gastric mucosa has been reported for the discovery of biomarkers and to describe the involvement of related biological pathways in liver cirrhosis induced PH. In this study, we performed a comparative analysis of the gastric mucosa proteome to obtain insights into the quantitative assessment and prognosis in liver cirrhosis patients via LC–MS/MS and PRM. As a result, we found that 423 proteins were differentially expressed in gastric from PHG patients before and after endoscopic cyanoacrylate injection compared to that from normal controls. Due to the importance of pathogenesis and early diagnosis of PHG, we studied the characterization and function of some differentially expressed proteins in depth.

According to KEGG function classification analysis, several pathways have been proved to be enriched in Ph/Nor and Inter/Nor group, but only focal adhesion pathway and vascular smooth muscle contraction pathway are enriched in both Ph/Nor and Inter/Nor group while others are not, indicating these two pathways may play critical roles in the process of PHG deterioration. In the process of liver fibrogenesis, focal adhesion connects hepatic stellate cells (HSCs) and ECM, and provides a direct sensor to the integrity of interaction with the extracellular environment, causing cell adhesion and migration [12, 13]. The patients with liver cirrhosis exhibit a progression of vascular smooth muscle contraction dysfunction, which presents as an increased intrahepatic vasoconstriction and reduced systemic vascular resistance [14]. The progress of PHG is related to hemodynamic alterations in the gastric mucosa, presenting as hemangiectasis [15], abnormal formation of arteriovenous shunts and blood flow augmentation [16].

Among these DEPs, PDLIM4 is the only protein which alters among the three groups, indicating cytoskeleton organization and stress fiber formation [17] may relate to development of PHG. Moreover, 17 DEPs about liver diseases and gastrointestinal system were further selected and validated by PRM based on our large-scale LS–MS/MS study. LS–MS/MS quantification is an agreeing approach to provide thorough and relative fold changes of thousands of peptides across numerous samples and to discover differential expressed peptides, based on so-called data-dependent acquisition. However, as a trade-off to the in-depth analysis, the obtained fold changes can be underestimated [18]. There DEPs are up-regulated in PHG patients and even more after endoscopic cyanoacrylate injection, indicating their expression level is consistent with severity of PHG. Among these DEPs, MYLK, TNC, ITGA1, FLNA, COL6A1, COL6A2 and COL6A3 are enriched in focal adhesion pathway. MYLK participate in smooth muscle contraction pathway at the same time. Activation of MYLK triggers gastrointestinal smooth muscle contraction with an increase of intracellular Ca2+ [Ca2+]I. Compared with antrum muscles, fundus muscles have a higher level of MYLK expression [19], which may contribute to the mechanism of gastric motility. In addition, MYLK also mediates cell–cell and cell-matric adhesion [20]. TNC can promote arterial smooth muscle cells proliferation and ECM elements deposition after stimulation of inflammatory cytokines and growth factors in cerebral arteries [21] and pulmonary arteries [22] by enhancing focal adhesion kinase phosphorylation [23]. These mechanisms may contribute to TNC mediated vascular remodeling [24, 25] and participate in the pathogenesis of subarachnoid hemorrhage and pulmonary arterial hypertension. FLNA couples cell cytoskeleton to ECM and regulates integrin alpha-1 (ITGA1), thus affects the function of collagen and ECM remodeling [26,27,28,29]. With the development of liver fibrosis, FLNA can be activated in HSCs and may serve as a biomarker [30]. Integrin modulates angiogenesis, regulates inflammation, and translates biliary injury into stimulation of matrix-producing mesenchymal cells in the process of fibrosis development [31]. Thereinto, ITGA1 is a receptor for laminin and collagen, it is involved in cell adhesion and correlates with inflammatory activity and fibrosis [32]. Collagen IV contributes to ECM structure and expresses extracellularly. It is structured by a heterotrimer of the alpha1(VI), alpha2(VI), and alpha3(VI) chains which are encoded by COL6A1, COL6A2, COL6A3, respectively [33].

Other validated proteins, though not enriched in the two pathways mentioned above, also play a role in the process of liver and gastrointestinal diseases. CPT1A and ACSL3, which are responsible for fatty acid oxidation, have been proved to stimulate gastric cancer and hepatocellular carcinoma cell growth and migration [34,35,36]. MCAM (CD146) involves in actin cytoskeleton rearrangement, blood vessel endothelial cells remodeling and intercellular junction [37]. CSRP2 expresses specifically in HSCs in liver and mediates the transdifferentiation of HSCs into myofibroblasts with the activation of TGF-β [38]. Focal adhesion protein TGFB1I1 (also known as Hic5) is induced by TGF-βand regulates myofibroblast differentiation in liver [39]. TNS1 is also localized to focal adhesions and acts as a bridge linking ECM and the actin cytoskeleton [40]. EFEMP1 and FBLN5 are ECM proteins and highly express in portal fibroblasts, they have been proved to play a role in progressive liver fibrosis [41, 42]. SORBS1 involves in cytoskeleton organization and insulin signaling pathway in human hepatoma cell line [43]. PDLIM7 is composed of PDZ and LIM7 domain, the LIM7 domain is associated with progression of liver fibrosis in hepatitis C virus infected patients [44]. Although the association with PHG is not quite clear, this study identified these DEPs as potential biomarkers on gastric mucosa for liver cirrhosis induced PH. However, more studies are required to make sure their availability as predictive biomarkers and diagnostic targets.


In this study, we identified quantitative differences in expression of several proteins on gastric mucosa for PHG patients. This proteomic study demonstrates that a LFQ analysis can be useful in detection of novel predictive biomarkers, and 17 promising proteins are further validated by PRM analysis. Our findings will provide a prospect for non-invasive early diagnosis and prediction of severity of liver cirrhosis induced portal hypertension, but more mechanisms of the interaction of these DEPs and their effect on the progression of PHG require further research.


Clinical samples

From March 2020 to June 2020, decompensated liver cirrhosis patients meeting the criteria below were admitted to the Department of Digestive Medicine, Shenzhen Hospital of Southern Medical University. Inclusion criteria included: (1) an age range of 15–85 years, (2) computed tomography (CT), laboratory examination and endoscopy confirmed the diagnosis of PHG according to Portal Hypertensive Bleeding in Cirrhosis: Risk Stratification, Diagnosis, and Management: 2016 Practice Guidance by the American Association for the Study of Liver Diseases, (3) Moderate to severe esophageal varices have no bleeding, but there is an obvious risk of bleeding, endoscopic primary prevention treatment is needed, and (4) Cases with a history of esophageal varices rupture and bleeding and secondary endoscopic prophylaxis. Exclusion criteria included acute EGVB, presinus and retrosinus portal hypertension portal hypertension, esophageal or gastric submucosal tumor, hypertrophic gastritis, and patients who have endoscopic examination contraindications. All individuals were excluded helicobacter pylori infection by Carbon-13 urea breath test or histologic examination of gastric mucosa. 5 healthy individuals who underwent physical examinations and participated in another research of our group [5] with no history of alcoholic or viral hepatitis and helicobacter pylori infection were assigned to the normal control group (Nor group), 5 liver cirrhosis patients with gastroscopy confirmed PHG were assigned to PHG group (Ph group). Ph group patients underwent gastroscopy reexamination 1 week after endoscopic cyanoacrylate injection (Inter group). Gastric mucosa samples were collected during gastroscopy by a biopsy forceps from the middle part of greater curvature by gastroscopic biopsy in Department of Endoscopy, Shenzhen Hospital of Southern Medical University. All samples were assigned to LFQ analysis, 2 individuals of each group were randomly selected to participate in PRM verification. Besides, 4 healthy individuals and 5 liver cirrhosis patients meeting the criteria above were assigned to PRM analysis. The study design was approved by the Shenzhen Hospital of Southern Medical University Ethic Committee complying with Declaration of Helsinki, Operational Guideline for the Ethic Review of Biomedical Research involving Human Subject [(2016) No.11] and International Ethical Guidelines for Biomedical Research Involving Human Subjects, CIOMS. All participants have signed the written informed consent.

Protein extraction and LS–MS/MS analysis

Proteins extraction and LS–MS/MS analysis were performed according to the method of Zheng and Xie et al. [45, 46]. Gastric mucosa samples were grinded by liquid nitrogen into cell powder and then transferred to a 5-mL centrifuge tube. After that, four volumes of lysis buffer (8 M urea, 1% Protease Inhibitor Cocktail) were added to the cell powder, followed by sonication three times on ice using a high intensity ultrasonic processor (Scientz). (Note: For PTM experiments, inhibitors were also added to the lysis buffer, e.g. 3 μM TSA and 50 mM NAM for acetylation.) The remaining debris was removed by centrifugation at 12,000 g at 4 °C for 10 min. Finally, the supernatant was collected and the protein concentration was determined with BCA kit according to the manufacturer’s instructions.

For digestion, the protein solution was reduced with 5 mM dithiothreitol for 30 min at 56 °C and alkylated with 11 mM iodoacetamide for 15 min at room temperature in darkness. The protein sample was then diluted by adding 100 mM TEAB to urea concentration less than 2 M. Finally, trypsin was added at 1:50 trypsin-to-protein mass ratio for the first digestion overnight and 1:100 trypsin-to-protein mass ratio for a second 4 h-digestion.

The tryptic peptides were dissolved in solvent A (0.1% formic acid, 2% acetonitrile in water), directly loaded onto a home-made reversed-phase analytical column ( ReproSil-Pur Basic C18 column(1.9 μm, 100 μm i.d., 25 cm)). Peptides were separated with a gradient from 4 to 22% solvent B (0.1% formic acid in acetonitrile) over 70 min, 22–30% in 14 min and climbing to 80% in 3 min then holding at 80% for the last 3 min, all at a constant flow rate of 450 nL/min on a nanoElute UHPLC system (Bruker Daltonics).

The peptides were subjected to Capillary source followed by the timsTOF Pro (Bruker Daltonics) mass spectrometry. The electrospray voltage applied was 1.75 kV. Precursors and fragments were analyzed at the TOF detector, with a MS/MS scan range from 100 to 1700 m/z. The timsTOF Pro was operated in parallel accumulation serial fragmentation (PASEF) mode. Precursors with charge states 0 to 5 were selected for fragmentation, and 10 PASEF-MS/MS scans were acquired per cycle. The dynamic exclusion was set to 30 s.

Database search

The resulting MS/MS data were processed using MaxQuant search engine (v. Tandem mass spectra were searched against the human uniprot database [Homo_sapiens_9606_SP_20200509 (20,366 entries)] concatenated with reverse decoy database. Trypsin/P was specified as cleavage enzyme allowing up to 2 missing cleavages. The mass tolerance for precursor ions was set as 20 ppm in the first search and 20 ppm in Main search, and the mass tolerance for fragment ions was set as 20 ppm. Carbamidomethyl on Cys was specified as fixed modification, and acetylation on protein N-terminal and oxidation on Met were specified as variable modifications. FDR was adjusted to < 1% and minimum score for modified peptides was set > 40.

Bioinformatics analyses

Differentially expressed protein analysis

In this study, the quantitative values of each sample in three replicates were obtained by LFQ intensity. For each protein, LFQ intensity was normalized by mean of all samples. And then normalized LFQ intensity were taken as log2 transform (so that the data conforms to the normal distribution), and p value was calculated by the two-sample two-tailed T-test method. When p value < 0.05 and protein ratio > 1.5 was regarded as up-regulation. When p value < 0.05 and protein ratio < 1/1.5 was regarded as down-regulation.

GO classification [47]

Gene Ontology (GO) annotation proteome was derived from the UniProt-GOA database (www. Firstly, DEPs were mapped to GO IDs by protein accession. If some DEPs were not annotated by UniProt-GOA database, the InterProScan soft would be used to annotated protein’s GO functional based on the protein sequence alignment method. Then DEPs were classified by Gene Ontology annotation based on three categories: biological process, cellular component and molecular function. A bar plot graph was used to present GO terms by visualization R package “ggplot2” in RStudio.

KEGG pathway enrichment [48]

KEGG database was used pathways. Firstly, using KEGG online service tools KAAS to annotated protein’s KEGG database description. Then mapping the annotation result on the KEGG pathway database using KEGG online service tools KEGG mapper. DEPs enriched pathways were identified by a two-tailed Fisher’s exact test. The pathway with p value < 0.05 was considered significant. A bubble plot graph was used to present enriched pathway by visualization R package “ggplot2”. All calculation and visualization steps were performed in RStudio.

Protein–protein interaction network

All DEPs accessions were searched against the STRING database version 11.0 for protein–protein interactions. Only interactions between the proteins belonging to the searched data set were selected, thereby excluding external candidates. STRING defines a metric called “confidence score” to define interaction confidence; we fetched all interactions that had a confidence score > 0.7 (high confidence). Interaction network form STRING was visualized in CytoScape software.

Parallel reaction monitoring (PRM)

We use PRM analysis to verify protein expression levels according to the method of Zhang et al. [49]. The resulting MS data were processed using Skyline (v. Peptide settings: enzyme was set as Trypsin [KR/P], Max missed cleavage set as 0. The peptide length was set as 7–25, Variable modification was set as Carbamidomethyl on Cys and oxidation on Met, and max variable modifications was set as 3. Transition settings: precursor charges were set as 2, 3, ion charges were set as 1, ion types were set as b, y, p. The product ions were set as from ion 3 to last ion, the ion match tolerance was set as 0.02 Da.


All data were processed using Microsoft Excel 2019 (Microsoft Corporation, Redmond, WA, USA). Data were expressed as the mean ± SEM, and comparisons of differences between two groups were analyzed using a one-way analysis of variance (ANOVA). Count data are expressed as a ratio, and the χ2 test was used for comparisons of differences between groups. p < 0.05 was designated to be of statistical significance. Statistical tests were conducted using SPSS v26.0 (SPSS, Chicago, IL, USA).

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.



Biological processes


Cellular components


Computed tomography


Differentially expressed proteins


Esophagogastric variceal bleeding


Extracellular matrix


Gene ontology


Gastroesophageal varices


Hepatic stellate cells


Kyoto Encyclopedia of Genes and Genomes


Label-free quantitative


Molecular functions


Magnetic resonance imaging


Portal hypertension


Portal hypertensive gastropathy


Protein–protein interaction


Parallel reaction monitoring


Transient elastography


  1. Fukui H, Saito H, Ueno Y, Uto H, Obara K, Sakaida I, et al. Evidence-based clinical practice guidelines for liver cirrhosis 2015. J Gastroenterol. 2016;51(7):629–50.

    Article  CAS  PubMed  Google Scholar 

  2. Garcia-Tsao G, Abraldes J, Berzigotti A, Bosch J. Portal hypertensive bleeding in cirrhosis: Risk stratification, diagnosis, and management: 2016 practice guidance by the American Association for the study of liver diseases. Hepatology. 2017;65(1):310–35.

    Article  PubMed  Google Scholar 

  3. de Franchis R. Revising consensus in portal hypertension: report of the Baveno V consensus workshop on methodology of diagnosis and therapy in portal hypertension. J Hepatol. 2010;53(4):762–8.

    Article  PubMed  Google Scholar 

  4. Tanoue K, Hashizume M, Wada H, Ohta M, Kitano S, Sugimachi K. Effects of endoscopic injection sclerotherapy on portal hypertensive gastropathy: a prospective study. Gastrointest Endosc. 1992;38(5):582–5.

    Article  CAS  PubMed  Google Scholar 

  5. Zhu Y, Wang F, Zhou Y, Xia G, Dong L, He W, et al. Blue laser magnifying endoscopy in the diagnosis of chronic gastritis. Exp Ther Med. 2019;18(3):1993–2000.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Bella M, Casas M, Vergara M, Brullet E, Junquera F, Martínez-Bauer E, et al. Utility of histology for the diagnosis of portal hypertensive gastroenteropathy. Concordance between the endoscopic image and gastrointestinal biopsies. Role of the CD34 marker. Gastroenterol Hepatol. 2019;42(3):150–6.

    Article  PubMed  Google Scholar 

  7. Westerhoff M, Tretiakova M, Hovan L, Miller J, Noffsinger A, Hart J. CD61, CD31, and CD34 improve diagnostic accuracy in gastric antral vascular ectasia and portal hypertensive gastropathy: an immunohistochemical and digital morphometric study. Am J Surg Pathol. 2010;34(4):494–501.

    Article  PubMed  Google Scholar 

  8. Fontana R, Sanyal A, Mehta S, Doherty M, Neuschwander-Tetri B, Everson G, et al. Portal hypertensive gastropathy in chronic hepatitis C patients with bridging fibrosis and compensated cirrhosis: results from the HALT-C trial. Am J Gastroenterol. 2006;101(5):983–92.

    Article  PubMed  Google Scholar 

  9. Payen J, Calès P, Voigt J, Barbe S, Pilette C, Dubuisson L, et al. Severe portal hypertensive gastropathy and antral vascular ectasia are distinct entities in patients with cirrhosis. Gastroenterology. 1995;108(1):138–44.

    Article  CAS  PubMed  Google Scholar 

  10. Xu J, Zhu Y, Wang F, Zhou Y, Xia G, Xu W. ICMT contributes to hepatocellular carcinoma growth, survival, migration and chemoresistance via multiple oncogenic pathways. Biochem Biophys Res Commun. 2019;518(3):584–9.

    Article  CAS  PubMed  Google Scholar 

  11. Zhu Y, Xu J, Hu W, Wang F, Zhou Y, Xu W, et al. TFAM depletion overcomes hepatocellular carcinoma resistance to doxorubicin and sorafenib through AMPK activation and mitochondrial dysfunction. Gene. 2020;753:144807.

    Article  CAS  PubMed  Google Scholar 

  12. Carloni V, Romanelli R, Pinzani M, Laffi G, Gentilini P. Focal adhesion kinase and phospholipase C gamma involvement in adhesion and migration of human hepatic stellate cells. Gastroenterology. 1997;112(2):522–31.

    Article  CAS  PubMed  Google Scholar 

  13. Parsons C, Takashima M, Rippe R. Molecular mechanisms of hepatic fibrogenesis. J Gastroenterol Hepatol. 2007.

    Article  PubMed  Google Scholar 

  14. Kajita M, Murata T, Horiguchi K, Iizuka M, Hori M, Ozaki H. iNOS expression in vascular resident macrophages contributes to circulatory dysfunction of splanchnic vascular smooth muscle contractions in portal hypertensive rats. Am J Physiol Heart Circ Physiol. 2011;300(3):H1021–31.

    Article  CAS  PubMed  Google Scholar 

  15. Iwao T, Toyonaga A, Sumino M, Takagi K, Oho K, Nishizono M, et al. Portal hypertensive gastropathy in patients with cirrhosis. Gastroenterology. 1992;102(6):2060–5.

    Article  CAS  PubMed  Google Scholar 

  16. Albillos A, Colombato L, Enriquez R, Ng O, Sikuler E, Groszmann R. Sequence of morphological and hemodynamic changes of gastric microvessels in portal hypertension. Gastroenterology. 1992;102(6):2066–70.

    Article  CAS  PubMed  Google Scholar 

  17. Vallenius T, Scharm B, Vesikansa A, Luukko K, Schäfer R, Mäkelä T. The PDZ-LIM protein RIL modulates actin stress fiber turnover and enhances the association of alpha-actinin with F-actin. Exp Cell Res. 2004;293(1):117–28.

    Article  CAS  PubMed  Google Scholar 

  18. Karp N, Huber W, Sadowski P, Charles P, Hester S, Lilley K. Addressing accuracy and precision issues in iTRAQ quantitation. Mol Cell Proteomics MCP. 2010;9(9):1885–97.

    Article  CAS  PubMed  Google Scholar 

  19. Li W, Sasse K, Bayguinov Y, Ward S, Perrino BJJ. Contractile protein expression and phosphorylation and contractility of gastric smooth muscles from obese patients and patients with obesity and diabetes. J Diabetes Res. 2018;2018:8743874.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lin J, He Y, Chen L, Chen X, Zang S, Lin W. MYLK promotes hepatocellular carcinoma progression through regulating cytoskeleton to enhance epithelial-mesenchymal transition. Clin Exp Med. 2018;18(4):523–33.

    Article  CAS  PubMed  Google Scholar 

  21. Borel C, McKee A, Parra A, Haglund M, Solan A, Prabhakar V, et al. Possible role for vascular cell proliferation in cerebral vasospasm after subarachnoid hemorrhage. Stroke. 2003;34(2):427–33.

    Article  CAS  PubMed  Google Scholar 

  22. Morrell N, Adnot S, Archer S, Dupuis J, Jones P, MacLean M, et al. Cellular and molecular basis of pulmonary arterial hypertension. J Am Coll Cardiol. 2009;54:S20-31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zagzag D, Shiff B, Jallo G, Greco M, Blanco C, Cohen H, et al. Tenascin-C promotes microvascular cell migration and phosphorylation of focal adhesion kinase. Cancer Res. 2002;62(9):2660–8.

    CAS  PubMed  Google Scholar 

  24. Suzuki H, Kanamaru K, Suzuki Y, Aimi Y, Matsubara N, Araki T, et al. Tenascin-C is induced in cerebral vasospasm after subarachnoid hemorrhage in rats and humans: a pilot study. Neurol Res. 2010;32(2):179–84.

    Article  CAS  PubMed  Google Scholar 

  25. Wallner K, Sharifi B, Shah P, Noguchi S, DeLeon H, Wilcox J. Adventitial remodeling after angioplasty is associated with expression of tenascin mRNA by adventitial myofibroblasts. J Am Coll Cardiol. 2001;37(2):655–61.

    Article  CAS  PubMed  Google Scholar 

  26. Xu Y, Bismar T, Su J, Xu B, Kristiansen G, Varga Z, et al. Filamin A regulates focal adhesion disassembly and suppresses breast cancer cell migration and invasion. J Exp Med. 2010;207(11):2421–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Krebs K, Ruusmann A, Simonlatser G, Velling TJE. Expression of FLNa in human melanoma cells regulates the function of integrin α1β1 and phosphorylation and localisation of PKB/AKT/ERK1/2 kinases. Eur J Cell Biol. 2015;94(12):564–75.

    Article  CAS  PubMed  Google Scholar 

  28. Karimi A, Milewicz D. Structure of the elastin-contractile units in the thoracic aorta and how genes that cause thoracic aortic aneurysms and dissections disrupt this structure. Can J Cardiol. 2016;32(1):26–34.

    Article  PubMed  Google Scholar 

  29. Mohammadi H, Pinto V, Wang Y, Hinz B, Janmey P, McCulloch C. Filamin A mediates wound closure by promoting elastic deformation and maintenance of tension in the collagen matrix. J Invest Dermatol. 2015;135(11):2852–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zhang H, Chen F, Fan X, Lin C, Hao Y, Wei H, et al. Quantitative proteomic analysis on activated hepatic stellate cells reversion reveal STAT1 as a key regulator between Liver Fibrosis and recovery. Sci Rep. 2017;7:44910.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Patsenker E, Stickel F. Role of integrins in fibrosing liver diseases. Am J Physiol Gastrointestinal Liv Physiol. 2011;301(3):G425–34.

    Article  CAS  Google Scholar 

  32. Patsenker E, Popov Y, Stickel F, Schneider V, Ledermann M, Sägesser H, et al. Pharmacological inhibition of integrin alphavbeta3 aggravates experimental liver fibrosis and suppresses hepatic angiogenesis. Hepatology. 2009;50(5):1501–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Cescon M, Gattazzo F, Chen P, Bonaldo PJJ. Collagen VI at a glance. J Cell Sci. 2015;128(19):3525–31.

    Article  CAS  PubMed  Google Scholar 

  34. Wang L, Li C, Song Y, Yan Z. Inhibition of carnitine palmitoyl transferase 1A-induced fatty acid oxidation suppresses cell progression in gastric cancer. Arch Biochem Biophys. 2020;696:108664.

    Article  CAS  PubMed  Google Scholar 

  35. Xu A, Wang B, Fu J, Qin W, Yu T, Yang Z, et al. Diet-induced hepatic steatosis activates Ras to promote hepatocarcinogenesis via CPT1α. Cancer Lett. 2019;442:40–52.

    Article  CAS  PubMed  Google Scholar 

  36. Ndiaye H, Liu J, Hall A, Minogue S, Morgan M, Waugh M. Immunohistochemical staining reveals differential expression of ACSL3 and ACSL4 in hepatocellular carcinoma and hepatic gastrointestinal metastases. 2020. Biosci Rep.

  37. Ouhtit A, Gaur R, Abd Elmageed Z, Fernando A, Thouta R, Trappey A, et al. Towards understanding the mode of action of the multifaceted cell adhesion receptor CD146. Biochim Biophys Acta. 2009;1795(2):130–6.

    Article  CAS  PubMed  Google Scholar 

  38. Herrmann J, Borkham-Kamphorst E, Haas U, Van de Leur E, Fraga M, Esteller M, et al. The expression of CSRP2 encoding the LIM domain protein CRP2 is mediated by TGF-beta in smooth muscle and hepatic stellate cells. Biochem Biophys Res Commun. 2006;345(4):1526–35.

    Article  CAS  PubMed  Google Scholar 

  39. Varney S, Betts C, Zheng R, Wu L, Hinz B, Zhou J, et al. Hic-5 is required for myofibroblast differentiation by regulating mechanically dependent MRTF-A nuclear accumulation. J Cell Sci. 2016;129(4):774–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Chen H, Duncan I, Bozorgchami H, Lo S. Tensin1 and a previously undocumented family member, tensin2, positively regulate cell migration. Proc Natl Acad Sci USA. 2002;99(2):733–8.

    Article  CAS  PubMed  Google Scholar 

  41. Chen Y, Gilbert M, Grochowski C, McEldrew D, Llewellyn J, Waisbourd-Zinman O, et al. A genome-wide association study identifies a susceptibility locus for biliary atresia on 2p16.1 within the gene EFEMP1. PLoS Genet. 2018;14(8):e1007532.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Bracht T, Schweinsberg V, Trippler M, Kohl M, Ahrens M, Padden J, et al. Analysis of disease-associated protein expression using quantitative proteomics—fibulin-5 is expressed in association with hepatic fibrosis. J Proteome Res. 2015;14(5):2278–86.

    Article  CAS  PubMed  Google Scholar 

  43. Lin W, Huang C, Liu M, Chang H, Chen Y, Tai T, et al. Cloning, mapping, and characterization of the human sorbin and SH3 domain containing 1 (SORBS1) gene: a protein associated with c-Abl during insulin signaling in the hepatoma cell line Hep3B. Genomics. 2001;74(1):12–20.

    Article  CAS  PubMed  Google Scholar 

  44. Ibrahim M, Salama H, Abd El Rahman M, Dawood R, Bader El Din N, Salem H, et al. Three gene signature for predicting the development of hepatocellular carcinoma in chronically infected hepatitis C virus patients. J Interferon Cytokine Res. 2016;36(12):698–705.

    Article  CAS  PubMed  Google Scholar 

  45. Zheng X, Yang Q, Zhao L, Apaliya M, Zhang X, Zhang H. Crosstalk between proteins expression and lysine acetylation in response to patulin stress in Rhodotorula mucilaginosa. Sci Rep. 2017;7(1):13490.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Xie H, Wang Y, Ding Y, Qiu C, Sun L, Gai Z, et al. Global ubiquitome profiling revealed the roles of ubiquitinated proteins in metabolic pathways of tea leaves in responding to drought stress. Sci Rep. 2019;9(1):4286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Wu Y, Li J, Qin X, Sun S, Xiao Z, Dong X, et al. Proteome and microbiota analysis reveals alterations of liver-gut axis under different stocking density of Peking ducks. PLoS ONE. 2018;13(10):e0198985.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Ying J, Zhao J, Hou Y, Wang Y, Qiu J, Li Z, et al. Mapping the N-linked glycosites of rice (Oryza sativa L.) germinating embryos. PLoS ONE. 2017;12(3):e0173853.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Zhang P, Zhu Y, Luo X, Zhou S. Comparative proteomic analysis provides insights into the complex responses to Pseudoperonospora cubensis infection of cucumber (Cucumis sativus L.). Sci Rep. 2019;9(1):9433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucl Acids Res. 2020;49(D1):D545–51.

    Article  Google Scholar 

  51. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucl Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  Google Scholar 

  52. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank PTM Biolabs (Hangzhou, China) for proteomic quantification.


This work was supported by research grants provided by National Natural Science Foundation of China (Nos. 81800503 and 81974070), the Basic Research Subject of Science and Technology plan in Bao’ an District (No. 2019JD444).

Author information

Authors and Affiliations



YZ designed the project and collected the samples and prepared Fig. 2. WX drafted the work and prepared Figs. 1, 3, 4, 5, 6, and 7. WH collected the data and prepared Table 1. FW and YZ performed the analyses, evaluated the results and prepared Table 2. JGX and WG made substantial contributions to the conception. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jianguo Xu or Wei Gong.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for this study was obtained from the Shenzhen Hospital of Southern Medical University Ethic Committee (NYSZYYEC20180003 and NYSZYYEC20200005) and clinical study registration was completed (ChiCTR2000029840). All subjects gave informed consent, and all treatment and management measures were monitored by the attending physicians.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

423 differentially expressed proteins (DEPs) between normal controls and PHG patients before or after operation.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, Y., Xu, W., Hu, W. et al. Discovery and validation of novel protein markers in mucosa of portal hypertensive gastropathy. BMC Gastroenterol 21, 214 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: