Next Article in Journal
CDC42-IQGAP Interactions Scrutinized: New Insights into the Binding Properties of the GAP-Related Domain
Next Article in Special Issue
The Human Milk Microbiota Produces Potential Therapeutic Biomolecules and Shapes the Intestinal Microbiota of Infants
Previous Article in Journal
Molecular Biomarkers in Glioblastoma: A Systematic Review and Meta-Analysis
Previous Article in Special Issue
The Gut Microbiome in Depression and Potential Benefit of Prebiotics, Probiotics and Synbiotics: A Systematic Review of Clinical Trials and Observational Studies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 16
10.3390/ijms23168841
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovering Biomarkers for Non-Alcoholic Steatohepatitis Patients with and without Hepatocellular Carcinoma Using Fecal Metaproteomics

by
Svenja Sydor
1,
Christian Dandyk
2,
Johannes Schwerdt
3,
Paul Manka
1,
Dirk Benndorf
2,4,*,
Theresa Lehmann
2,
Kay Schallert
2,4,5,
Maximilian Wolf
6,
Udo Reichl
2,4,
Ali Canbay
1,
Lars P. Bechmann
1 and
Robert Heyer
2,4,5,6,*
1
Department of Internal Medicine, University Hospital Knappschaftskrankenhaus, Ruhr-University Bochum, In der Schornau 23–25, 44892 Bochum, Germany
2
Bioprocess Engineering, Otto von Guericke University Magdeburg, Universitätsplatz 2, 39106 Magdeburg, Germany
3
Data and Knowledge Engineering Group, Otto von Guericke University Magdeburg, Universitätsplatz 2, 39106 Magdeburg, Germany
4
Max Planck Institute for Dynamics of Complex Technical Systems, Bioprocess Engineering, Otto von Guericke University Magdeburg, Sandtorstraße 1, 39106 Magdeburg, Germany
5
Multidimensional Omics Analysis Group, Leibniz-Institut für Analytische Wissenschaften—ISAS—e.V., Bunsen-Kirchhoff-Straße 11, 44139 Dortmund, Germany
6
Faculty of Technology, Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(16), 8841; https://doi.org/10.3390/ijms23168841
Submission received: 1 July 2022 / Revised: 1 August 2022 / Accepted: 3 August 2022 / Published: 9 August 2022
(This article belongs to the Special Issue Human Gut Microbiome and Diet in Health and Disease)
Browse Figures
Versions Notes

Abstract

:
High-calorie diets lead to hepatic steatosis and to the development of non-alcoholic fatty liver disease (NAFLD), which can evolve over many years into the inflammatory form of non-alcoholic steatohepatitis (NASH), posing a risk for the development of hepatocellular carcinoma (HCC). Due to diet and liver alteration, the axis between liver and gut is disturbed, resulting in gut microbiome alterations. Consequently, detecting these gut microbiome alterations represents a promising strategy for early NASH and HCC detection. We analyzed medical parameters and the fecal metaproteome of 19 healthy controls, 32 NASH patients, and 29 HCC patients, targeting the discovery of diagnostic biomarkers. Here, NASH and HCC resulted in increased inflammation status and shifts within the composition of the gut microbiome. An increased abundance of kielin/chordin, E3 ubiquitin ligase, and nucleophosmin 1 represented valuable fecal biomarkers, indicating disease-related changes in the liver. Although a single biomarker failed to separate NASH and HCC, machine learning-based classification algorithms provided an 86% accuracy in distinguishing between controls, NASH, and HCC. Fecal metaproteomics enables early detection of NASH and HCC by providing single biomarkers and machine learning-based metaprotein panels.

1. Introduction

During the last few decades, the prevalence of obesity and metabolic syndrome has increased tremendously [1]. Due to this, non-alcoholic fatty liver disease (NAFLD) has emerged as one of the leading causes of chronic liver diseases worldwide [2]. The NAFLD spectrum ranges from visible fatty degeneration of the organ, but can progress to the inflammatory form of non-alcoholic steatohepatitis (NASH) [3,4]. Due to local and systemic inflammatory processes, this chronic inflammation carries the risk of further disease progression and cancer development [5]. Different mechanisms related to fat and glucose metabolism can promote fibrotic remodeling of the organ, the development of cirrhosis, and even hepatocellular carcinoma (HCC), whereas HCC is one of the most common cancers and causes of cancer-associated deaths worldwide [6,7]. To date, not all mechanisms affecting the progression of NAFLD have been elucidated. For example, fibrosis is a main driver of HCC development, whereas NASH can develop into HCC even without prior cirrhosis [8].
Timely diagnosis and progress prediction of NAFLD, NASH, and HCC is a major challenge, as it remains clinically inconspicuous for long periods and lacks appropriate biomarkers applicable for preventive medical checkups. Several factors and mechanisms that affect NAFLD’s progress are currently discussed, such as the presence of diabetes and metabolic syndrome, and the composition of intestinal microbiota [9]. We and others described the interaction between the gut and liver diseases within the enterohepatic circulation. Shifts in the composition of the intestinal microbiome affect bile acid composition and the formation of bioactive metabolites and substrates that can affect fatty acid and glucose metabolism, and thus, also promote NAFLD and its progress [10,11,12,13,14]. Based on sequencing methods, gut microbiota compositions have previously been described in different measures and NAFLD patient cohorts. However, most findings were descriptive, and the underlying mechanisms have not been fully elucidated. In contrast to monitoring the taxonomic composition by sequencing methods, metaproteomics detects the actual gene expression. Furthermore, fecal metaproteomics also reveals host proteins (e.g., from the immune system), indicating health status [15] or problems with food digestion [16].
Consequently, metaproteomics is a promising method for preventive, non-invasive medical checkups, but we still lack meaningful biomarkers. To identify the required biomarkers and to understand the correlation between the pathogenesis and the gut microbiome in NASH and NASH-derived HCC, we analyzed the human proteins and the taxonomic and functional gut microbiome composition through fecal metaproteomics.

2. Results

2.1. Characteristics of the Study Cohort

Eighty subjects including healthy controls (n = 19, 23.5 average age, and 23.0 BMI), NASH patients (n = 32, 53.3 average age, and 30.9 BMI), and HCC patients (n = 29, 67.8 average age, and 30.3 BMI), were included. The increasing age and BMI from healthy controls to NASH and HCC patients reflected the progression of disease over time due to elevated BMI. Demographic data of transient elastography (TE) and different serum parameters of individual patient groups are depicted in Supplementary Table S1.
Patients (NASH and HCC) showed significantly increased liver stiffness and hepatic fat accumulation compared to controls, as assessed by TE, including measurement of the controlled attenuation parameter (CAP). Tumor markers such as alpha-fetoprotein-Centaur (AFP), lectin-3-reactive alpha-fetoprotein (AFP-L3), and des-gamma-carboxyprothrombin (DCP) also showed a slight increase in NASH; however, this increase was below the established cutoff values for tumor diagnosis in NASH, but highly increased in HCC. Furthermore, NASH and HCC patients possessed elevated serum parameters of liver injury such as alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (AP), gamma-glutamyltransferase (γGT), glutamate dehydrogenase (GLDH), and lactate dehydrogenase (LDH). Serum levels of bilirubin were not changed when comparing the groups.
Total bile acids and their target, fibroblast growth factor 19 (FGF19), were increased in the serum of NASH and were even higher in HCC. Regarding the triglyceride (TAG) levels, a significantly higher concentration could be observed for NASH and HCC. Although not significantly, in NASH and HCC patients, low-density lipoprotein (LDL) cholesterol was increased, whereas high-density lipoprotein (HDL) cholesterol levels showed a significant decrease. As important metabolic mediators, we measured serum levels of adiponectin, glucagon-like peptide-1 (GLP1), and FGF21. GLP1 and FGF21 levels were significantly increased in NASH and HCC compared to controls, whereas adiponectin levels were elevated but without statistical significance. Cell death marker M65, apoptosis marker M30, and serum levels of pro-inflammatory cytokines interleukin 6 (IL6) and tumor necrosis factor alpha (TNFα) were significantly increased in NASH and HCC. This increase was accompanied by higher levels of the inflammation-associated protein lipoprotein-binding protein 1 (LBP1) (Supplementary Table S1).

2.2. Characterization of Fecal Metaproteomics

The metaproteomic analysis yielded an average of 19,221 ± 6752 identified spectra for each patient, which were assigned, on average, to 4218 ± 940 metaproteins and 193 ± 20 taxonomic families (Supplementary Tables S2 and S3). The taxonomic assignment of all identified spectra (Supplementary Figure S1) resulted in an average of 11.41 ± 2.51% bacterial spectra, 0.71 ± 0.20% archaeal spectra, 5.41 ± 1.79% eukaryotic spectra, and 0.32 ± 0.10% viral spectra. Furthermore, 30% of the identified spectra belonged to unknown protein entries from the metagenome, and for 36.5%, no specific superkingdom could be assigned due to overlapping protein identifications.
The set of eukaryotic spectra was divided into a fraction belonging to the host (Hominidae: 1.78 ± 1.01%) and a fraction belonging to photosynthetically active eukaryotes related to diet (e.g., Poaceae: 0.39 ± 0.27%). However, most eukaryotic spectra lack a sufficiently precise taxonomic assignment to assign them to one of these two groups without detailed functional interpretation. The most important microbial families were Bacillaceae, 2.01 ± 0.55%; Enterobacteriaceae, 1.70 ± 0.42%; Clostridiaceae, 0.84 ± 0.25%; Mycobacteriaceae, 0.40 ± 0.14%; and Pasteurellaceae, 0.37 ± 0.12%.
The main metabolic functions were summarized based on Biemann et al. 2021 [16], comprising human and microbial hydrolysis, microbial metabolisms and transporters, and host proteins derived from the intestinal barrier and the immune system (Supplementary Table S2).

2.3. Identification of Disease-Specific Metaprotein Patterns

To identify possible disease-specific patterns, we compared the human (family Hominidae) and bacterial metaproteins with the measured medical parameters and taxonomic composition at the family level using multilayer principal component analysis (PCA) (Figure 1) and an analysis of similarities (ANOSIM) (Supplementary Table S4).
The medical parameters allowed a separation of the healthy and diseased people, matching weak (0.12 < R < 0.20, p-value < 0.01) but significant differences found by ANOSIM for all three groups. Based on microbial metaproteins and families, only a separation between controls and the diseased was possible by the ANOSIM, and a weak trend in the PCAs was observable (Figure 1). In contrast, the ANOSIM of the human proteins enabled no significant differences between the three groups (p-values > 0.26). However, the human metaprotein profiles of the healthy individuals were closer together than those of the diseased patients.
Interlayer connections between the PCAs showed that the control samples varied less than the diseased samples across all layers. Analysis of the PCA loadings showed that NASH and HCC were correlated, among others, with increased blood fat levels, albumin levels, age, and a lower abundance of the families Clostridiaceae and Enterobacteriaceae (Figure 1). Furthermore, NASH and HCC patients’ feces contained more antibodies and metaproteins associated with the gut barrier and immune system (e.g., polymeric immunoglobulin receptor). Regarding microbial metaproteins, correlations with low-abundant metaproteins such as MP6 and a probable serine/threonine protein kinase, SPs1, were observable.

2.4. Significantly Altered Metaproteins, Taxonomies, and Functions

To identify potential diagnostic NASH and HCC biomarkers, we considered all metabolic functions from our summary and all metaproteins and families whose identified spectra abundance was above 0.01% (Supplementary Tables S3 and S5). In total, 15 of 40 functions and 9 of 34 families were significantly altered (p-value < 0.01). For the metaproteins, we increased the significance threshold to p < 10−5, revealing 25 of 961 changed metaproteins.
Among others, we observed in NASH and HCC patients more proteins assigned to Hominidae (spectral abundance (#SpecAb)_C: 1.08%, #SpecAb_NASH: 1.92%, #SpecAb_HCC: 2.07%) (Figure 2) and Thermotogaceae (#SpecAb_C: 0.11%, #SpecAb_NASH: 0.15%, #SpecAb_HCC: 0.18%; Figure 3), and less assigned to Enterobacteriaceae (#SpecAb_C: 1.92%, #SpecAb_NASH: 1.55%, #SpecAb_HCC: 1.73%), Clostridiaceae (#SpecAb_C: 1.00%, #SpecAb_NASH: 0.76%, #SpecAb_HCC: 0.83%), and Lactobacillaceae (#SpecAb_C: 0.16%, #SpecAb_NASH: 0.12%, #SpecAb_HCC: 0.14%) (Figure 3).
Furthermore, in NASH and HCC patients more proteins for the intestinal barrier (#SpecAb_C: 0.88%, #SpecAb_NASH: 1.66%, #SpecAb_HCC: 1.65%) and neutrophil granulocytes were detected (#SpecAb_C: 0.77%, #SpecAb_NASH: 1.36%, #SpecAb_HCC: 1.32%). Within the microbiome, we observed decreased microbial metabolism (e.g., butyrate fermentation; #SpecAb_C: 3.42%, #SpecAb_NASH: 1.58%, #SpecAb_HCC: 1.92%) and transporters (e.g., sugar transport; #SpecAb_C: 4.96%, #SpecAb_NASH: 2.67%, #SpecAb_HCC: 2.73%) in NASH and HCC patients. Exceptions were transporters for vitamin B12 (#SpecAb_C: 0.03%, #SpecAb_NASH: 0.05%, #SpecAb_HCC: 0.06%) and lactate fermentation (#SpecAb_C: 0.04%, #SpecAb_NASH: 0.07%, #SpecAb_HCC: 0.07%), being more abundant in NASH and HCC patients. Potential marker metaproteins for NASH and HCC were a decreased abundance of the sn-glycerol-3-phosphate import ATP-binding protein (#SpecAb_C: 0.55%%, #SpecAb_NASH: 0.20%, #SpecAb_HCC: 0.25%; unknown superkingdom) and ketol-acid reductoisomerase (NADP(+)) (#SpecAb_C: 0.51%, #SpecAb_NASH: 0.20%, #SpecAb_HCC: 0.27%), and increased abundances for the kielin/chordin-like protein (#SpecAb_C: 0.84%, #SpecAb_NASH: 2.96%, #SpecAb_HCC: 3.28%; class: Mammalia) and protein S100-A9 (#SpecAb_C: 0.09%, #SpecAb_NASH: 0.34%, #SpecAb_HCC: 0.39%; unknown superkingdom). Furthermore, we observed in NASH and HCC patients an increased abundance of the E3 ubiquitin ligase (#SpecAb_C: 0.04%, #SpecAb_NASH: 0.14%, #SpecAb_HCC: 0.20%) (Figure 3). Although this metaprotein was assigned to the fungal species Arthroderma otae, we would speculate that it belongs actually to the host, since it was much more abundant than other low-abundant metaproteins assigned to fungi.
Unfortunately, identified markers only enabled to separate between controls and diseased people, but not between NASH and HCC. The only exceptions were an increased abundance of Pasteurellaceae (#SpecAb_C: 0.33%, #SpecAb_NASH: 0.34%, #SpecAb_HCC: 0.42%%) and Pseudomonadaceae (#SpecAb_C: 0.26%%, #SpecAb_NASH: 0.26%%, #SpecAb_HCC: 0.32%%) in the feces of HCC patients and a decreased ratio between Firmicutes and Bacteriodetes in NASH samples (Figure 2).

2.5. Potential Biomarkers to Distinguish NASH and HCC from Controls

In the next step, we evaluated the performance of the significantly changed top ten metaproteins for separating between healthy and diseased patients (Supplementary Table S5). Therefore, we performed an ROC curve analysis and compared the area under the curve (Table 1). For the analyzed metaproteins, the area varied between 0.913 and 0.815, indicating a good classification. Exemplary, for the kielin/chordin-like protein, about 80% of the diseased people could be diagnosed with only 10% false positives. However, the number of false positives was still too high for routine diagnosis or preventive medical checkups.

2.6. Machine Learning-Based Biomarker Panels to Separate NASH from HCC and Controls

To improve the diagnosis of NASH and HCC, we developed machine learning-based classification algorithms (Figure 4A). Therefore, we ranked the normalized features according to their p-values derived from a t-test. Subsequently, we applied the wrapper technique on the top-ranked features to further reduce the set to the most relevant molecules for the classification task. Under several models, the diagonal linear discriminant analysis and logistic regression performed best and enabled a separation of controls from NASH or HCC at 99.98% and 100% using seven and five features, respectively. In contrast, the correct distinction between NASH and HCC was only possible in 86.4% of samples using ten features. Consequently, the correct classification of all three groups was only possible in 86.0% of all cases using eleven features (Table 2). Thereby, the NASH samples were either wrongly classified as healthy or HCC, and HCC patients were wrongly classified as NASH patients (Figure 4C). Misclassification of HCC as controls or other was not observed. Evaluating all selected features manually (Supplementary Table S6) revealed nucleophosmin as a promising biomarker. We observed a considerable overlap of identified molecules between both the machine learning and ROC curve analysis approaches, as described in Section 2.5. Ergo, both approaches complement each other in evidence. It was enriched in NASH and HCC compared to controls by a factor of 103, or rather, 129 (#SpecAb_C: 0.00%, #SpecAb_NASH: 0.026%, #SpecAb_HCC: 0.033%).

3. Discussion

NASH and HCC are severe liver diseases that progressively reduce liver function as the disease progresses, as observed through worsened liver function parameters, liver fibrosis, and liver damage [17]. The progressive carcinogenesis in HCC patients coincided with the increased tumor markers such as AFP-Centaur, AFP-L3, and DCP [18]. A primary cause of NASH and HCC is an unhealthy diet, which is reflected in a higher BMI, worse blood lipids, and sugar parameters. The unhealthy lifestyle leads over time, as observable with the higher age of the NASH and HCC patients, to fibrosis and tumor formation [19]. NASH- and HCC-induced liver alterations indeed result in changed bile acid production and liver protein expression, leading to changes in the gut microbiome via bile acid secretion from the gallbladder. Therefore, gut microbiome alterations may indicate NASH or HCC, or could even promote it via bile acid conjugation or ethanol production [20].
Although there was no clear separation of the metaproteome of NASH and HCC patients from the controls, differences in the metaprotein and family fingerprint were already observable in the PCA plots. Thereby, differences between healthy and diseased patients were bigger than between NASH and HCC, reflecting that HCC often develops from NASH [21,22]. The bigger variations in the PCA plots between the NASH and HCC patients than within the controls show how diseases possess different forms and severities.
Subsequent identification of potential marker functions, families, and metaproteins revealed multiple potential biomarkers enabling the separation between healthy and diseased patients. However, as the ROC plots and comparison of the area under the curve showed, the biomarker accuracy was insufficient for separating NASH from HCC. A better performance provided machine learning-based classification using five to seven metaproteins. Therefore, developing comprehensive clinical panels for the fecal metaproteome similar to a blood picture may represent a promising approach for NASH and HCC diagnosis. The particular advantage of fecal metaproteomics is that samples can be taken non-invasively at home and sent to clinical laboratories, making it a perfect approach for preventive clinical checkups.
The overlap between the presented top PCA loadings, significantly altered metaproteins, taxonomies, functions, and machine learning-based features appeared small. However, the latter groups were very similar since they were based on significantly changed metaproteins, but we focused within this manuscript on different aspects. We applied a smaller p-value cutoff for the manual selection and concentrated on high-abundant marker proteins with a potential clinical significance. For machine learning, the algorithms selected the smallest number of features, enabling the best separation of the groups. Furthermore, evaluation of the PCA loadings revealed that most loadings also differed significantly between controls, NASH, and HCC, but not all. This observation reflects that the PCA visualizes the variance of the samples, and the top loadings indicate the main differences. However, as observable in the PCA plots, NASH, HCC, and control samples were only weakly separated, also suggesting other differences within the microbiome.
Although we identified several promising marker metaproteins in our cohort, they required a clinical evaluation since they could also be linked with other diseases or with lifestyle. Elevated levels of kielin/chordin, an E3 ubiquitin ligase, and nucleophosmin 1 could be directly involved in the pathogenesis. NASH results in the accumulation of adipose cells in the liver, secreting tumor growth factor (TGFβ). TGFβ indeed is induced via endoplasmic reticulum (ER) stress and the unfolded protein response, apoptosis, and thus, fibrosis [23,24,25]. The E3 ubiquitin ligase is involved in protein degradation and nucleophosmin 1 in nucleic transport and ribosome biosynthesis. Upregulation of both would fit to enhance the unfolded protein response. Furthermore, both are described as enhanced in liver cancers [25,26]. Kielin/chordin indeed represses TGFß signaling, representing a protection mechanism against NASH. Soofi et al. showed that kielin/chordin knockout mice were more susceptible to developing hepatic steatosis and liver fibrosis [27]. Conversely, overexpression of kielin/chordin protected the mice’s liver from the effects of an excessively high-fat diet.
In contrast to these three disease-specific biomarkers, increased abundance of the family Hominidae combined with more metaproteins derived from the immune system (e.g., calprotectin from neutrophilic granulocytes) and the gut barrier reflected a worsened health status of NASH and HCC patients in contrast to the healthy control. However, this is not specific to NASH or HCC. For example, Lehmann et al. [15] observed higher abundances of calprotectin (protein S100-A9) in the feces of patients with inflammatory bowel disease, and Biemann et al. (2021) [16] showed that obese patients possess a systemic inflammation, higher abundances of protein S100-A9, and increased abundance of the family Thermotogaceae.
Similarly, obesity and diet also explain gut microbiome alterations, including decreased microbial transporters or less butyrate fermentation. Increased food uptake elevates the production and secretion of bile acids into the gut. Bile acids indeed possess antimicrobial properties, leading to an altered gut microbiome [28]. Although obesity is linked with an increased abundance of nutrients, the observed increase in vitamin B12 transporters suggests a higher competition for vitamins and a potential lack. In line with this, Voland et al. (2021) [29] reviewed the lack of vitamins in obese people and its impact on the gut microbiome, which could also explain the altered ratio of Firmicutes to Bacteriodetes in NASH patients [30]. On the contrary, there are also some research articles suggesting that the gut microbiome contributes to NASH and HCC by the production of toxic components such as alcohol, toxic bile acids, or inflammatory microbial metabolites and the activation of cancerogenesis-associated signaling pathways [10,31,32,33]. We observed some hints, such as an increase in NASH and HCC patients of probable serine/threonine protein kinase SPs1, which is associated with signaling [34], and of pyruvate decarboxylase isozyme 3 for ethanol production [35]. However, the evidence was insufficient, or the alteration was not significant.
The main shortcoming of the study was that no patient’s specific metagenomes were available. Therefore, and to keep this study comparable with previous studies [15,16], we selected the same metagenome database. However, specific metagenomes would increase the number of identified metaproteins, the taxonomic and functional protein annotation, and the protein grouping. Since our study focused on identifying fecal marker proteins, we included no liver or epithelial biopsies. Therefore, spatial protein assignment was impossible, but would be valuable for mechanistic studies. For example, the identified polymeric immunoglobulin receptor is upregulated in the liver of patients with liver fibrosis and liver cancer [36,37]. However, the observed increase is more likely caused by the degradation of the intestinal epithelia expressing the receptor in high amounts.
Analogous to all other NASH and HCC studies, our study suffers from the dependency of NASH and HCC on the cofounding factors of lifestyle, age, and co-morbidity with other diseases. For example, NASH is caused by a calorie-rich diet and is a metabolic syndrome. Thus, it often occurs together with diabetes and hypertension. In future studies, documentation of dietary habits of the patients would be useful since diet may influence the composition of the gut microbiota and may have an impact on the development of diseases, especially cancer [38].

4. Conclusions

In conclusion, we proposed, by fecal metaproteomics, several potential biomarkers enabling the separation of NASH and HCC patients from healthy people, presenting a valuable tool for preventive medical checkups. An even better diagnosis than single marker proteins, our findings provide machine learning-based biomarker panels.

5. Materials and Methods

5.1. Patient Recruitment and Sample Collection

This study was conducted based on a previous study by Sydor et al. (2020), analyzing potential links between the liver and the gut in NASH-related hepatocarcinogenesis. Therefore, they compared alterations of gut microbiota and mediators of bile acid signaling in the absence or presence of cirrhosis through analysis of feces and serum from patients with NASH and NASH-HCC and healthy volunteers [14].
The Ethics Committee (Institutional Review Board) of the University Hospital Essen (reference number: 14-6044-BO) approved the study, and all subjects provided informed written consent. The study protocol conformed to the ethical guidelines of the Declaration of Helsinki.
For the analysis, serum and fecal samples of subjects with NASH (n = 32), HCC (n = 29), and healthy controls (n = 19) without evidence of NAFLD were analyzed. The inclusion criteria for the study were the presence of NAFLD, specifically with appropriately confirmed NASH, and HCC based on NASH. Known chronic viral, toxic, hereditary, and immunologic liver diseases were considered exclusion criteria (e.g., HBV, HCV, autoimmune hepatitis, primary and secondary biliary cholangitis (PBC/PSC), Wilson’s disease, etc.).
Confounding factors of our study were age and BMI, which could not be mitigated since we wanted to focus on the development of NASH and HCC over time due to increased calorie uptake.
Diagnosis of NASH and HCC was performed as described before [14]. Patients with significant alcohol intake, as defined as consuming more than two standard drinks daily or more than six daily drinks on weekends for at least five years [39], were not considered for the study. The presence of HBV and HCV was excluded by seronegativity for HBV or HCV following standard laboratory tests. Healthy volunteers with a BMI below 30 and without NAFLD were selected as healthy controls.
All serum samples were collected in a fasted state in the morning and stored at −80 °C until measurement. The central laboratory of the University Hospital Essen evaluated by routine diagnostics the general clinical parameters, enzymes (ALT, AST, AP, γGT), total bile acids, and tumor markers (AFP, AFP-L3, and DCP). Fecal samples were collected from every patient in sterile tubes and stored at −80 °C.

5.2. Transient Elastography and Controlled Attenuation Parameter

Liver stiffness and the controlled attenuation parameter (CAP) to assess hepatic fat accumulation were measured using the Fibroscan® (Echosens, Paris, France), with samples taken from the subjects in a fasted state.

5.3. ELISA

Commercially available kits were used to measure serum levels of the overall cell death marker M65 and apoptosis marker M30 (TecoMedical, Sissach, Switzerland). Quantification of serum concentrations of FGF19, FGF21, GLP1, IL6, and TNFα was performed using the specific Quantikine ELISA Kit from R&D Systems (Minneapolis, MN, USA). Serum amounts of LBP1 were quantified using the LBP ELISA Kit (Hycult Biotech Uden, Uden, The Netherlands). All procedures were performed following the manufacturer’s instructions.

5.4. Fecal Sample Preparation for Metaproteomics

Proteins from approx. 100–200 mg stool samples were extracted by cell lysis and phenol extraction as described in Lehmann et al. [15]. After FASP digestion [40], LC-MS/MS analysis was performed using an UltiMate 3000 RSLCnano splitless liquid chromatography system coupled online to an Orbitrap Elite ™ hybrid ion trap, the Orbitrap-MS (both from Thermo Fisher Scientific, Bremen, Germany) using a 120 min gradient. All chemicals used were at least analysis quality and the solvents used were LC-MS/MS quality.

5.5. Data Handling

The MetaProteomeAnalyzer (version 3.1) [40] was used for protein identification, which included the search engines X! Tandem, OMSSA, and Mascot and the following parameters: enzyme trypsin, one missed cleavage, monoisotopic mass, carbamidomethylation (cysteine) as a fixed modification, oxidation (methionine) as a variable modification, ±10 ppm precursor and ±0.5 Da MS/MS fragment tolerance, 1 13 C, +2/+3 charged peptide ions, and a false detection rate of 1%. The used protein database was the UniProtKB/Swiss-Prot database (as of 16 January 2019) combined with a human gut microbiome database [41]. A BLAST search (NCBI-Blast version 2.2.31) against UniProtKB/Swiss-Prot was carried out for proteins that could not be annotated taxonomically or functionally. All BLAST hits with the best E-value that were at least below 10−4 were combined and used to annotate the protein identifications. Redundant homologous protein identifications were combined into a protein group (also referred to as metaprotein) if they had at least one peptide identification in common. Finally, all results were uploaded to PRIDE (Accession: PXD034175).

5.6. Statistical Analysis

Statistical analysis, including the multilayer PCA, ANOSIM, Kruskal–Wallis test, and violin plots [42], was carried out using R Statistics (version 1.2.5001) and Rstudio. For Krona plots, the provided Excel template by Ondov et al. [43], and for the ROC plots, the web service from Eng et al. (2014) [44] were used. Power analysis for metaproteomics using standard deviations of three spectra from a previous study [15] showed that for 20 samples per group and for proteins with an abundance of at least five spectra, a doubling of the spectra could be observed with a power of 0.993 and with a significance value below 0.01.

5.7. Development of a Biomarker Panel

A comprehensive software package using R and Java was used to develop machine learning-based classification algorithms. The software ranked the normalized features according to their t-test-based predictive power. We considered only features when 2/3 of the samples in at least one group had measurements (values above zero). Subsequently, different feature sets were identified by the wrapper method. The following machine learning models were used: linear discriminant analysis (LDA) [45], diagonal LDA [46], logistic regression [47], support-vector machine [48], random forest [49], extremely randomized trees [50], and k-nearest neighbors [51]. The evaluation was performed by averaging a 5-fold cross-validation on 100 repeats for the tree-like models and 10,000 repeats for the rest. Results were summarized in a confusion matrix and a clustergram using Ward linkage and Canberra distances.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23168841/s1.

Author Contributions

Conceptualization, R.H., S.S. and L.P.B.; methodology, R.H. and T.L.; software, K.S. and J.S.; investigation, C.D., R.H. and S.S.; resources, P.M.; writing—original draft preparation, S.S. and R.H.; writing—review and editing, D.B., L.P.B., U.R. and A.C.; visualization: J.S. and M.W.; supervision, R.H. and L.P.B.; project administration, R.H.; funding acquisition, R.H. and D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the German Federal Ministry of Education and Research (de.NBI network project MetaProtServ grant No. 031L0103) and the German Research Society (DFG, grant No. HE 8077/1-1).

Institutional Review Board Statement

The Ethics Committee (Institutional Review Board) of the University Hospital Essen (reference number: 14-6044-BO) approved the study. The study protocol conformed to the ethical guidelines of the Declaration of Helsinki.

Informed Consent Statement

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

Data Availability Statement

The datasets analyzed during the current study are accessible in the PRIDE Archive (Proteomics Identifications Database) under the accession number PXD034175. For reasons of data protection, patient data are presented in summarized form, but can be requested from the authors by personal inquiry.

Acknowledgments

We acknowledge Corina Siewert for laboratory support. This manuscript is dedicated to the memory of Professor Lars Bechmann, a passionate, excellent researcher and dear friend.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AFP/AFP-L3alpha-fetoprotein/lectin-3-reactive alpha-fetoprotein
ALTalanine aminotransferase
ANOSIManalysis of similarities
APalkaline phosphatase
ASTaspartate aminotransferase
BMIbody mass index
CAPcontrolled attenuation parameter
DCPdes-gamma-carboxyprothrombin
ERendoplasmic reticulum
FGF19/21fibroblast growth factor 19/21
γGTgamma-glutamyltransferase
GLDHglutamate dehydrogenase
GLP1glucagon-like peptide
HCChepatocellular carcinoma
HDLhigh-density lipoprotein
IL6interleukin 6
LBP1lipoprotein-binding protein 1
LDHlactate dehydrogenase
LDLlow-density lipoprotein
NAFLDnon-alcoholic fatty liver disease
NASHnon-alcoholic steatohepatitis
PCAprincipal component analysis
ROCreceiver operating characteristic
#SpecAbspectral abundance
TAGtriglyceride
TEtransient elastography
TGFβtumor growth factor beta
TNFαtumor necrosis factor alpha

References

  1. Younossi, Z.; Anstee, Q.M.; Marietti, M.; Hardy, T.; Henry, L.; Eslam, M.; George, J.; Bugianesi, E. Global Burden of NAFLD and NASH: Trends, Predictions, Risk Factors and Prevention. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 11–20. [Google Scholar] [CrossRef] [PubMed]
  2. Boyle, M.; Masson, S.; Anstee, Q.M. The Bidirectional Impacts of Alcohol Consumption and the Metabolic Syndrome: Cofactors for Progressive Fatty Liver Disease. J. Hepatol. 2018, 68, 251–267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yeh, M.M.; Brunt, E.M. Pathology of Nonalcoholic Fatty Liver Disease. Am. J. Clin. Pathol. 2007, 128, 837–847. [Google Scholar] [CrossRef] [PubMed]
  4. Powell, E.E.; Wong, V.W.-S.; Rinella, M. Non-Alcoholic Fatty Liver Disease. Lancet 2021, 397, 2212–2224. [Google Scholar] [CrossRef]
  5. Bugianesi, E. Non-Alcoholic Steatohepatitis and Cancer. Clin. Liver Dis. 2007, 11, 191–207. [Google Scholar] [CrossRef] [PubMed]
  6. Bechmann, L.P.; Hannivoort, R.A.; Gerken, G.; Hotamisligil, G.S.; Trauner, M.; Canbay, A. The Interaction of Hepatic Lipid and Glucose Metabolism in Liver Diseases. J. Hepatol. 2012, 56, 952–964. [Google Scholar] [CrossRef] [Green Version]
  7. El-Serag, H.B. Hepatocellular Carcinoma. N. Engl. J. Med. 2011, 365, 1118–1127. [Google Scholar] [CrossRef]
  8. Ertle, J.; Dechêne, A.; Sowa, J.-P.; Penndorf, V.; Herzer, K.; Kaiser, G.; Schlaak, J.F.; Gerken, G.; Syn, W.-K.; Canbay, A. Non-Alcoholic Fatty Liver Disease Progresses to Hepatocellular Carcinoma in the Absence of Apparent Cirrhosis. Int. J. Cancer 2011, 128, 2436–2443. [Google Scholar] [CrossRef]
  9. Friedman, S.L.; Neuschwander-Tetri, B.A.; Rinella, M.; Sanyal, A.J. Mechanisms of NAFLD Development and Therapeutic Strategies. Nat. Med. 2018, 24, 908–922. [Google Scholar] [CrossRef]
  10. Yoshimoto, S.; Loo, T.M.; Atarashi, K.; Kanda, H.; Sato, S.; Oyadomari, S.; Iwakura, Y.; Oshima, K.; Morita, H.; Hattori, M.; et al. Obesity-Induced Gut Microbial Metabolite Promotes Liver Cancer through Senescence Secretome. Nature 2013, 499, 97–101. [Google Scholar] [CrossRef]
  11. Mouzaki, M.; Wang, A.Y.; Bandsma, R.; Comelli, E.M.; Arendt, B.M.; Zhang, L.; Fung, S.; Fischer, S.E.; McGilvray, I.G.; Allard, J.P. Bile Acids and Dysbiosis in Non-Alcoholic Fatty Liver Disease. PLoS ONE 2016, 11, e0151829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Konturek, P.C.; Harsch, I.A.; Konturek, K.; Schink, M.; Konturek, T.; Neurath, M.F.; Zopf, Y. Gut-Liver Axis: How Do Gut Bacteria Influence the Liver? Med. Sci. 2018, 6, 79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Jiang, J.-W.; Chen, X.-H.; Ren, Z.; Zheng, S.-S. Gut Microbial Dysbiosis Associates Hepatocellular Carcinoma via the Gut-Liver Axis. Hepatobiliary Pancreat. Dis. Int. 2019, 18, 19–27. [Google Scholar] [CrossRef] [PubMed]
  14. Sydor, S.; Best, J.; Messerschmidt, I.; Manka, P.; Vilchez-Vargas, R.; Brodesser, S.; Lucas, C.; Wegehaupt, A.; Wenning, C.; Aßmuth, S.; et al. Altered Microbiota Diversity and Bile Acid Signaling in Cirrhotic and Noncirrhotic NASH-HCC. Clin. Transl. Gastroenterol. 2020, 11, e00131. [Google Scholar] [CrossRef]
  15. Lehmann, T.; Schallert, K.; Vilchez-Vargas, R.; Benndorf, D.; Püttker, S.; Sydor, S.; Schulz, C.; Bechmann, L.; Canbay, A.; Heidrich, B.; et al. Metaproteomics of Fecal Samples of Crohn’s Disease and Ulcerative Colitis. J. Proteom. 2019, 201, 93–103. [Google Scholar] [CrossRef] [PubMed]
  16. Biemann, R.; Buß, E.; Benndorf, D.; Lehmann, T.; Schallert, K.; Püttker, S.; Reichl, U.; Isermann, B.; Schneider, J.G.; Saake, G.; et al. Fecal Metaproteomics Reveals Reduced Gut Inflammation and Changed Microbial Metabolism Following Lifestyle-Induced Weight Loss. Biomolecules 2021, 11, 726. [Google Scholar] [CrossRef]
  17. Kupčová, V.; Fedelešová, M.; Bulas, J.; Kozmonová, P.; Turecký, L. Overview of the Pathogenesis, Genetic, and Non-Invasive Clinical, Biochemical, and Scoring Methods in the Assessment of NAFLD. Int. J. Environ. Res. Public Health 2019, 16, 3570. [Google Scholar] [CrossRef] [Green Version]
  18. Zhou, L.; Liu, J.; Luo, F. Serum Tumor Markers for Detection of Hepatocellular Carcinoma. World J. Gastroenterol. 2006, 12, 1175–1181. [Google Scholar] [CrossRef] [Green Version]
  19. Negro, F. Natural History of NASH and HCC. Liver Int. 2020, 40 (Suppl. S1), 72–76. [Google Scholar] [CrossRef] [Green Version]
  20. Bashiardes, S.; Shapiro, H.; Rozin, S.; Shibolet, O.; Elinav, E. Non-Alcoholic Fatty Liver and the Gut Microbiota. Mol. Metab. 2016, 5, 782–794. [Google Scholar] [CrossRef]
  21. Michelotti, G.A.; Machado, M.V.; Diehl, A.M. NAFLD, NASH and Liver Cancer. Nat. Rev. Gastroenterol. Hepatol. 2013, 10, 656–665. [Google Scholar] [CrossRef] [PubMed]
  22. Anstee, Q.M.; Reeves, H.L.; Kotsiliti, E.; Govaere, O.; Heikenwalder, M. From NASH to HCC: Current Concepts and Future Challenges. Nat. Rev. Gastroenterol. Hepatol. 2019, 16, 411–428. [Google Scholar] [CrossRef] [PubMed]
  23. Biernacka, A.; Dobaczewski, M.; Frangogiannis, N.G. TGF-β Signaling in Fibrosis. Growth Factors 2011, 29, 196–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kandeil, M.A.; Hashem, R.M.; Mahmoud, M.O.; Hetta, M.H.; Tohamy, M.A. Zingiber Officinale Extract and Omega-3 Fatty Acids Ameliorate Endoplasmic Reticulum Stress in a Nonalcoholic Fatty Liver Rat Model. J. Food Biochem. 2019, 43, e13076. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, Z.; Li, C.; Kang, N.; Malhi, H.; Shah, V.H.; Maiers, J.L. Transforming Growth Factor β (TGFβ) Cross-Talk with the Unfolded Protein Response Is Critical for Hepatic Stellate Cell Activation. J. Biol. Chem. 2019, 294, 3137–3151. [Google Scholar] [CrossRef] [Green Version]
  26. Sun, Y. E3 Ubiquitin Ligases as Cancer Targets and Biomarkers. Neoplasia 2006, 8, 645–654. [Google Scholar] [CrossRef] [Green Version]
  27. Soofi, A.; Wolf, K.I.; Emont, M.P.; Qi, N.; Martinez-Santibanez, G.; Grimley, E.; Ostwani, W.; Dressler, G.R. The Kielin/Chordin-like Protein (KCP) Attenuates High-Fat Diet-Induced Obesity and Metabolic Syndrome in Mice. J. Biol. Chem. 2017, 292, 9051–9062. [Google Scholar] [CrossRef] [Green Version]
  28. Tian, Y.; Gui, W.; Koo, I.; Smith, P.B.; Allman, E.L.; Nichols, R.G.; Rimal, B.; Cai, J.; Liu, Q.; Patterson, A.D. The Microbiome Modulating Activity of Bile Acids. Gut Microbes 2020, 11, 979–996. [Google Scholar] [CrossRef]
  29. Voland, L.; Le Roy, T.; Debédat, J.; Clément, K. Gut Microbiota and Vitamin Status in Persons with Obesity: A Key Interplay. Obes. Rev. 2022, 23, e13377. [Google Scholar] [CrossRef]
  30. Xu, Y.; Xiang, S.; Ye, K.; Zheng, Y.; Feng, X.; Zhu, X.; Chen, J.; Chen, Y. Cobalamin (Vitamin B12) Induced a Shift in Microbial Composition and Metabolic Activity in an in Vitro Colon Simulation. Front. Microbiol. 2018, 9, 2780. [Google Scholar] [CrossRef] [Green Version]
  31. Xie, G.; Wang, X.; Huang, F.; Zhao, A.; Chen, W.; Yan, J.; Zhang, Y.; Lei, S.; Ge, K.; Zheng, X.; et al. Dysregulated Hepatic Bile Acids Collaboratively Promote Liver Carcinogenesis. Int. J. Cancer 2016, 139, 1764–1775. [Google Scholar] [CrossRef] [PubMed]
  32. Miura, S.; Mitsuhashi, N.; Shimizu, H.; Kimura, F.; Yoshidome, H.; Otsuka, M.; Kato, A.; Shida, T.; Okamura, D.; Miyazaki, M. Fibroblast Growth Factor 19 Expression Correlates with Tumor Progression and Poorer Prognosis of Hepatocellular Carcinoma. BMC Cancer 2012, 12, 56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, Y.; Zhang, W.; Doughtie, A.; Cui, G.; Li, X.; Pandit, H.; Yang, Y.; Li, S.; Martin, R. Up-Regulation of Fibroblast Growth Factor 19 and Its Receptor Associates with Progression from Fatty Liver to Hepatocellular Carcinoma. Oncotarget 2016, 7, 52329–52339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Canova, M.J.; Molle, V. Bacterial Serine/Threonine Protein Kinases in Host-Pathogen Interactions. J. Biol. Chem. 2014, 289, 9473–9479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhu, L.; Baker, S.S.; Gill, C.; Liu, W.; Alkhouri, R.; Baker, R.D.; Gill, S.R. Characterization of Gut Microbiomes in Nonalcoholic Steatohepatitis (NASH) Patients: A Connection between Endogenous Alcohol and NASH. Hepatology 2013, 57, 601–609. [Google Scholar] [CrossRef]
  36. Yue, X.; Ai, J.; Xu, Y.; Chen, Y.; Huang, M.; Yang, X.; Hu, B.; Zhang, H.; He, C.; Yang, X.; et al. Polymeric Immunoglobulin Receptor Promotes Tumor Growth in Hepatocellular Carcinoma. Hepatology 2017, 65, 1948–1962. [Google Scholar] [CrossRef]
  37. Zhang, Y.; Lu, W.; Chen, X.; Cao, Y.; Yang, Z. A Bioinformatic Analysis of Correlations between Polymeric Immunoglobulin Receptor (PIGR) and Liver Fibrosis Progression. BioMed Res. Int. 2021, 2021, 5541780. [Google Scholar] [CrossRef]
  38. Samraj, A.N.; Pearce, O.M.T.; Läubli, H.; Crittenden, A.N.; Bergfeld, A.K.; Banda, K.; Gregg, C.J.; Bingman, A.E.; Secrest, P.; Diaz, S.L.; et al. A Red Meat-Derived Glycan Promotes Inflammation and Cancer Progression. Proc. Natl. Acad. Sci. USA 2015, 112, 542–547. [Google Scholar] [CrossRef] [Green Version]
  39. Dawson, D.A.; Li, T.-K.; Grant, B.F. A Prospective Study of Risk Drinking: At Risk for What? Drug Alcohol Depend. 2008, 95, 62–72. [Google Scholar] [CrossRef] [Green Version]
  40. Heyer, R.; Schallert, K.; Büdel, A.; Zoun, R.; Dorl, S.; Behne, A.; Kohrs, F.; Püttker, S.; Siewert, C.; Muth, T.; et al. A Robust and Universal Metaproteomics Workflow for Research Studies and Routine Diagnostics Within 24 h Using Phenol Extraction, FASP Digest, and the MetaProteomeAnalyzer. Front. Microbiol. 2019, 10, 1883. [Google Scholar] [CrossRef] [Green Version]
  41. Human Microbiome Jumpstart Reference Strains Consortium; Nelson, K.E.; Weinstock, G.M.; Highlander, S.K.; Worley, K.C.; Creasy, H.H.; Wortman, J.R.; Rusch, D.B.; Mitreva, M.; Sodergren, E.; et al. A Catalog of Reference Genomes from the Human Microbiome. Science 2010, 328, 994–999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kanai, M.; Maeda, Y.; Okada, Y. Grimon: Graphical Interface to Visualize Multi-Omics Networks. Bioinformatics 2018, 34, 3934–3936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Ondov, B.D.; Bergman, N.H.; Phillippy, A.M. Interactive Metagenomic Visualization in a Web Browser. BMC Bioinform. 2011, 12, 385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. ROC Analysis: Online ROC Curve Calculator. Available online: http://www.rad.jhmi.edu/jeng/javarad/roc/JROCFITi.html (accessed on 23 May 2022).
  45. Fisher, R.A. The Use of Multiple Measurements in Taxonomic Problems. Ann. Eugen. 1936, 7, 179–188. [Google Scholar] [CrossRef]
  46. Dudoit, S.; Fridlyand, J.; Speed, T.P. Comparison of Discrimination Methods for the Classification of Tumors Using Gene Expression Data. J. Am. Stat. Assoc. 2002, 97, 77–87. [Google Scholar] [CrossRef] [Green Version]
  47. Berkson, J. Application of the Logistic Function to Bio-Assay. J. Am. Stat. Assoc. 1944, 39, 357–365. [Google Scholar] [CrossRef]
  48. Vapnik, V. The Nature of Statistical Learning Theory; Springer: New York, NY, USA, 2000. [Google Scholar]
  49. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  50. Geurts, P.; Ernst, D.; Wehenkel, L. Extremely Randomized Trees. Mach. Learn. 2006, 63, 3–42. [Google Scholar] [CrossRef] [Green Version]
  51. Fix, E.; Hodges, J.L. Discriminatory Analysis—Nonparametric Discrimination: Consistency Properties. Int. Stat. Rev. 1989, 57, 238–247. [Google Scholar] [CrossRef]
Ijms 23 08841 g001 550
Figure 1. Multilayer PCA of clinical parameters, human and metaproteins, and family-level taxonomy for all samples. The individual PCAs were visualized on the left, and same samples were connected across all layers. The associated biplots with the top 10 loadings were shown on the right side. For human and microbial metaproteins, the metaproteins with at least 0.01% of the total spectral count were selected. In contrast, for families and clinical parameters, all 419 and 30 were chosen, respectively. For better readability, the top ten human and microbial loadings (metaproteins) were summarized in a table below the plot.
Figure 1. Multilayer PCA of clinical parameters, human and metaproteins, and family-level taxonomy for all samples. The individual PCAs were visualized on the left, and same samples were connected across all layers. The associated biplots with the top 10 loadings were shown on the right side. For human and microbial metaproteins, the metaproteins with at least 0.01% of the total spectral count were selected. In contrast, for families and clinical parameters, all 419 and 30 were chosen, respectively. For better readability, the top ten human and microbial loadings (metaproteins) were summarized in a table below the plot.
Ijms 23 08841 g001
Ijms 23 08841 g002 550
Figure 2. Abundance of the family Hominidae and the ratio of Bacteriodetes to Firmicutes. The abundance is based on the normalized abundance of identified spectra.
Figure 2. Abundance of the family Hominidae and the ratio of Bacteriodetes to Firmicutes. The abundance is based on the normalized abundance of identified spectra.
Ijms 23 08841 g002
Ijms 23 08841 g003 550
Figure 3. Summary of selected changed features. Significance was evaluated by the Kruskal–Wallis test using a p-value cutoff smaller than 0.01 for families and metabolic functions and smaller than 10−5 for metaproteins.
Figure 3. Summary of selected changed features. Significance was evaluated by the Kruskal–Wallis test using a p-value cutoff smaller than 0.01 for families and metabolic functions and smaller than 10−5 for metaproteins.
Ijms 23 08841 g003
Ijms 23 08841 g004 550
Figure 4. Machine learning-based sample classification between NASH, HCC, and controls. (A) shows the workflow of the software for feature selection, feature wrapping, and development of the classification algorithms. (B) shows display of the confusion matrix of the best-observed classifier (linear discriminant analysis.). Cross-validation ensured that the patient numbers in the confusion matrix were real and not natural numbers. The evaluation was performed by averaging a 5-fold cross-validation on 10,000 repeats. (C) shows the clustering of the data and their intrinsic similarity using Ward linkage and Canberra distances.
Figure 4. Machine learning-based sample classification between NASH, HCC, and controls. (A) shows the workflow of the software for feature selection, feature wrapping, and development of the classification algorithms. (B) shows display of the confusion matrix of the best-observed classifier (linear discriminant analysis.). Cross-validation ensured that the patient numbers in the confusion matrix were real and not natural numbers. The evaluation was performed by averaging a 5-fold cross-validation on 10,000 repeats. (C) shows the clustering of the data and their intrinsic similarity using Ward linkage and Canberra distances.
Ijms 23 08841 g004
Table
Table 1. Potential biomarker metaproteins between controls and diseased patients. We analyzed the ROC plot analysis for the most abundant ten metaproteins and summarized the area under the curve to evaluate metaprotein biomarkers.
Table 1. Potential biomarker metaproteins between controls and diseased patients. We analyzed the ROC plot analysis for the most abundant ten metaproteins and summarized the area under the curve to evaluate metaprotein biomarkers.
Metaproteins#SpecAbArea under Curve
Kielin/chordin-like protein2.568%0.893
Sn-glycerol-3-phosphate import ATP-binding protein0.303%0.868
Ketol-acid reductoisomerase (NADP(+))0.297%0.862
Protein S100-A90.296%0.815
Probable E3 ubiquitin ligase complex SCF0.135%0.839
30S ribosomal protein S30.120%0.879
Formate-tetrahydrofolate ligase 20.073%0.913
30S ribosomal protein S20.066%0.842
Acyl-CoA dehydrogenase, short-chain specific0.066%0.883
Glyceraldehyde-3-phosphate dehydrogenate0.063%0.905
Table
Table 2. Classification accuracy for NASH, HCC, and controls. Results were obtained by the given number of features and the algorithms by averaging a 5-fold cross-validation on 10,000 repeats.
Table 2. Classification accuracy for NASH, HCC, and controls. Results were obtained by the given number of features and the algorithms by averaging a 5-fold cross-validation on 10,000 repeats.
ComparisonAccuracyNumber of Features
NASH vs. Control0.99987 features
HCC vs. Control:15 features
HCC vs. NASH0.864010 features
HCC vs. NASH vs. Control0.8611 features
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sydor, S.; Dandyk, C.; Schwerdt, J.; Manka, P.; Benndorf, D.; Lehmann, T.; Schallert, K.; Wolf, M.; Reichl, U.; Canbay, A.; et al. Discovering Biomarkers for Non-Alcoholic Steatohepatitis Patients with and without Hepatocellular Carcinoma Using Fecal Metaproteomics. Int. J. Mol. Sci. 2022, 23, 8841. https://doi.org/10.3390/ijms23168841

AMA Style

Sydor S, Dandyk C, Schwerdt J, Manka P, Benndorf D, Lehmann T, Schallert K, Wolf M, Reichl U, Canbay A, et al. Discovering Biomarkers for Non-Alcoholic Steatohepatitis Patients with and without Hepatocellular Carcinoma Using Fecal Metaproteomics. International Journal of Molecular Sciences. 2022; 23(16):8841. https://doi.org/10.3390/ijms23168841

Chicago/Turabian Style

Sydor, Svenja, Christian Dandyk, Johannes Schwerdt, Paul Manka, Dirk Benndorf, Theresa Lehmann, Kay Schallert, Maximilian Wolf, Udo Reichl, Ali Canbay, and et al. 2022. "Discovering Biomarkers for Non-Alcoholic Steatohepatitis Patients with and without Hepatocellular Carcinoma Using Fecal Metaproteomics" International Journal of Molecular Sciences 23, no. 16: 8841. https://doi.org/10.3390/ijms23168841

APA Style

Sydor, S., Dandyk, C., Schwerdt, J., Manka, P., Benndorf, D., Lehmann, T., Schallert, K., Wolf, M., Reichl, U., Canbay, A., Bechmann, L. P., & Heyer, R. (2022). Discovering Biomarkers for Non-Alcoholic Steatohepatitis Patients with and without Hepatocellular Carcinoma Using Fecal Metaproteomics. International Journal of Molecular Sciences, 23(16), 8841. https://doi.org/10.3390/ijms23168841

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop