Chromosomal microarray analysis in fetuses with high-risk prenatal indications: A retrospective study in China.

OBJECTIVE: The present study aimed to determine the diagnostic value of prenatal chromosomal microarray analysis (CMA) for fetuses with several indications of being at high risk for various conditions. MATERIALS AND METHODS: This retrospective analysis included 1256 pregnancies that were prenatally evaluated due to high-risk indications using invasive CMA. The indications for invasive prenatal diagnosis mainly included ultrasound anomalies, high-risk for maternal serum screening (MSS), high-risk for non-invasive prenatal tests (NIPT), family history of genetic disorders or birth defects, and advanced maternal age (AMA). The rate of clinically significant genomic imbalances between the different groups was compared. RESULTS: The overall prenatal diagnostic yield was 98 (7.8%) of 1256 pregnancies. Clinically significant genomic aberrations were identified in 2 (1.5%) of 132 patients with non-structural ultrasound anomalies, 36 (12.7%) of 283 with structural ultrasound anomalies, 2 (4.5%) of 44 at high-risk for MSS, 38 (26.6%) of 143 at high-risk for NIPT, 11 (3.8%) of 288 with a family history, and 7 (2.1%) of 328 with AMA. Submicroscopic findings were identified in 29 fetuses, 19 of whom showed structural ultrasound anomalies. CONCLUSION: The diagnostic yields of CMA for pregnancies with different indications greatly varied. CMA could serve as a first-tier test for structural anomalies, especially multiple anomalies, craniofacial dysplasia, urinary defects, and cardiac dysplasia. Our results have important implications for genetic counseling.

Integrative single-cell transcriptome analysis reveals a subpopulation of fibroblasts associated with favorable prognosis of liver cancer patients.

Single-cell transcriptome analysis has provided detailed insights into the ecosystem of liver cancer. However, the changes of the cellular and molecular components of liver tumors in comparison with normal livers have not been described at single-cell level. Here, we performed an integrative single-cell analysis of both normal livers and liver cancers. Principal component analysis was firstly performed to delineate the cell lineages in liver tissues. Differential gene expression within major cell types were then analyzed between tumor and normal samples, thus resolved the cell type-specific molecular alterations in liver cancer development. Moreover, a comparison between liver cancer derived versus normal liver derived cell components revealed that two subpopulations of fibroblasts were exclusively expanded in liver cancer tissues. By further defining subpopulation-specific gene signatures, characterizing their spatial distribution in tumor tissues and investigating their clinical significance, we found that the SPARCL1 positive fibroblasts, representing a group of tumor vessel associated fibroblasts, were related to reduced vascular invasion and prolonged survival of liver cancer patients. Through establishing an in-vitro endothelial-to-mesenchymal transition model, we verified the conversion of the fetal liver sinusoidal endothelial cells into the fibroblast-like cells, demonstrating a possible endothelial cell origination of the SPARCL1 positive fibroblasts. Our study provides new insights into the cell atlas alteration, especially the expanded fibroblasts in liver cancers.

Dissecting transcriptional heterogeneity in primary gastric adenocarcinoma by single cell RNA sequencing.

OBJECTIVE: Tumour heterogeneity represents a major obstacle to accurate diagnosis and treatment in gastric adenocarcinoma (GA). Here, we report a systematic transcriptional atlas to delineate molecular and cellular heterogeneity in GA using single-cell RNA sequencing (scRNA-seq). DESIGN: We performed unbiased transcriptome-wide scRNA-seq analysis on 27 677 cells from 9 tumour and 3 non-tumour samples. Analysis results were validated using large-scale histological assays and bulk transcriptomic datasets. RESULTS: Our integrative analysis of tumour cells identified five cell subgroups with distinct expression profiles. A panel of differentiation-related genes reveals a high diversity of differentiation degrees within and between tumours. Low differentiation degrees can predict poor prognosis in GA. Among them, three subgroups exhibited different differentiation grade which corresponded well to histopathological features of Lauren's subtypes. Interestingly, the other two subgroups displayed unique transcriptome features. One subgroup expressing chief-cell markers (eg, LIPF and PGC) and RNF43 with Wnt/ß-catenin signalling pathway activated is consistent with the previously described entity fundic gland-type GA (chief cell-predominant, GA-FG-CCP). We further confirmed the presence of GA-FG-CCP in two public bulk datasets using transcriptomic profiles and histological images. The other subgroup specifically expressed immune-related signature genes (eg, LY6K and major histocompatibility complex class II) with the infection of Epstein-Barr virus. In addition, we also analysed non-malignant epithelium and provided molecular evidences for potential transition from gastric chief cells into MUC6+TFF2+ spasmolytic polypeptide expressing metaplasia. CONCLUSION: Altogether, our study offers valuable resource for deciphering gastric tumour heterogeneity, which will provide assistance for precision diagnosis and prognosis.

Susceptibility Factors of Stomach for SARS-CoV-2 and Treatment Implication of Mucosal Protective Agent in COVID-19.

Objectives: This work aims to study the gastrointestinal (GI) symptoms in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected patients and the susceptibility factors of the stomach for SARS-CoV-2. Materials and Methods: We investigated the SARS-CoV-2 susceptibility by analyzing the expression distribution of viral entry-associated genes, ACE2 and TMPRSS2, in single-cell RNA sequencing data derived from 12 gastric mucosa samples. We also analyzed the epidemiological, demographic, clinical, and laboratory data of 420 cases with SARS-CoV-2-caused coronavirus disease 2019 (COVID-19). Results: ACE2 and TMPRSS2 are specifically expressed in enterocytes which are mainly from gastric mucosa samples with Helicobacter pylori (H. pylori) infection history and intestinal metaplasia (IM). A total of 420 patients were surveyed, of which 62 were with and 358 were without GI symptoms. There is a significant difference in average hospital stay (p < 0.001) and cost (p < 0.001) between the two groups. Among 23 hospitalized patients including seven with upper GI symptoms and 16 with lower GI symptoms, six (85.7%) and five (31.3%) had H. pylori infection history, respectively (p = 0.03). Of 18 hospitalized patients with initial upper GI symptoms, none of the eight patients with mucosal protective agent therapy (e.g., sucralfate suspension gel, hydrotalcite tablets) had diarrhea subsequently, whereas six out of 10 patients without mucosal protective agent therapy had diarrhea subsequently (p = 0.01). Conclusion: IM and H. pylori infection history may be susceptibility factors of SARS-CoV-2, and the mucosal protective agent may be useful for the blockade of SARS-CoV-2 transmission from the stomach to the intestine.

Immune receptor repertoires in pediatric and adult acute myeloid leukemia.

BACKGROUND: Acute myeloid leukemia (AML), caused by the abnormal proliferation of immature myeloid cells in the blood or bone marrow, is one of the most common hematologic malignancies. Currently, the interactions between malignant myeloid cells and the immune microenvironment, especially T cells and B cells, remain poorly characterized. METHODS: In this study, we systematically analyzed the T cell receptor and B cell receptor (TCR and BCR) repertoires from the RNA-seq data of 145 pediatric and 151 adult AML samples as well as 73 non-tumor peripheral blood samples. RESULTS: We inferred over 225,000 complementarity-determining region 3 (CDR3) sequences in TCR α, ß, γ, and δ chains and 1,210,000 CDR3 sequences in B cell immunoglobulin (Ig) heavy and light chains. We found higher clonal expansion of both T cells and B cells in the AML microenvironment and observed many differences between pediatric and adult AML. Most notably, adult AML samples have significantly higher level of B cell activation and more secondary Ig class switch events than pediatric AML or non-tumor samples. Furthermore, adult AML with highly expanded IgA2 B cells, which might represent an immunosuppressive microenvironment, are associated with regulatory T cells and worse overall survival. CONCLUSIONS: Our comprehensive characterization of the AML immune receptor repertoires improved our understanding of T cell and B cell immunity in AML, which may provide insights into immunotherapies in hematological malignancies.

[Diagnosis of two cases from one family with Joubert syndrome caused by novel mutations of TCTN1 gene by whole exome sequencing].

OBJECTIVE: To explore the pathogenesis of two fetuses from one family affected with Joubert syndrome (JS). METHODS: Whole exome sequencing was employed to screen potential mutations in both fetuses. Suspected mutations were verified by Sanger sequencing. Impact of intronic mutations on DNA transcription was validated by cDNA analysis. RESULTS: Two novel TCTN1 mutations, c.342-8A>G and c.1494+1G>A, were identified in exons 2 and 12, respectively.cDNA analysis confirmed the pathogenic nature of both mutations with interference of normal splicing resulting in production of truncated proteins. CONCLUSION: The genetic etiology of the family affected with JS has been identified.Above findings have enriched the mutation spectrum of TCTN1gene and facilitated understanding of the genotype-phenotype correlation of JS.

Loss of PICH promotes chromosome instability and cell death in triple-negative breast cancer.

Triple-negative breast cancer (TNBC), defined by the lack of expression of estrogen, progesterone, and ERBB2 receptors, has the worst prognosis of all breast cancers. It is difficult to treat owing to a lack of effective molecular targets. Here, we report that the growth of TNBC cells is exceptionally dependent on PICH, a DNA-dependent ATPase. Clinical samples analysis showed that PICH is highly expressed in TNBC compared to other breast cancer subtypes. Importantly, its high expression correlates with higher risk of distal metastasis and worse clinical outcomes. Further analysis revealed that PICH depletion selectively impairs the proliferation of TNBC cells, but not that of luminal breast cancer cells, in vitro and in vivo. In addition, knockdown of PICH in TNBC cells induces the formation of chromatin bridges and lagging chromosomes in anaphase, frequently resulting in micronucleation or binucleation, finally leading to mitotic catastrophe and apoptosis. Collectively, our findings show the dependency of TNBC cells on PICH for faithful chromosome segregation and the clinical potential of PICH inhibition to improve treatment of patients with high-risk TNBC.

Predicting the regrowth of clinically non-functioning pituitary adenoma with a statistical model.

BACKGROUND: Compared with clinically functioning pituitary adenoma (FPA), clinically non-functioning pituitary adenoma (NFPA) lacks of detectable hypersecreting serum hormones and related symptoms which make it difficult to predict the prognosis and monitoring for postoperative tumour regrowth. We aim to investigate whether the expression of selected tumour-related proteins and clinical features could be used as tumour markers to effectively predict the regrowth of NFPA. METHOD: Tumour samples were collected from 295 patients with NFPA from Beijing Tiantan Hospital. The expression levels of 41 tumour-associated proteins were assessed using tissue microarray analyses. Clinical characteristics were analysed via univariate and multivariate logistic regression analyses. Logistic regression algorithm was applied to build a prediction model based on the expression levels of selected proteins and clinical signatures, which was then assessed in the testing set. RESULTS: Three proteins and two clinical signatures were confirmed to be significantly related to the regrowth of NFPA, including cyclin-dependent kinase inhibitor 2A (CDKN2A/p16), WNT inhibitory factor 1 (WIF1), tumour growth factor beta (TGF-ß), age and tumour volume. A prediction model was generated on the training set, which achieved a fivefold predictive accuracy of 81.2%. The prediction ability was validated on the testing set with an accuracy of 83.9%. The area under the receiver operating characteristic curves (AUC) for the signatures were 0.895 and 0.881 in the training and testing sets, respectively. CONCLUSION: The prediction model could effectively predict the regrowth of NFPA, which may facilitate the prognostic evaluation and guide early interventions.

Publisher Correction: Landscape of B cell immunity and related immune evasion in human cancers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Up-regulation of FAM64A promotes epithelial-to-mesenchymal transition and enhances stemness features in breast cancer cells.

FAM64A was found to be markedly up-regulated in tumor samples and associated with worse overall survival in multiple cancer types, including breast cancer. However, the functional significance of FAM64A in breast cancer remains largely unknown. In this study, we systematically investigated the expression of FAM64A in multiple public breast cancer datasets. We found that FAM64A is significantly positively correlated with tumor stemness index in breast cancer samples, corresponding with an advanced clinical grade, metastasis and unfavorable prognosis. In vitro experiments further showed an up-regulation of stemness genes after over-expressing FAM64A in breast cancer cells. FAM64A overexpression also promoted breast cancer cell proliferation, migration, accompanied by the activation of epithelial-to-mesenchymal transition (EMT). Besides, we identified a strong association of FAM64A expression with TP53 mutations in TCGA and three additional breast cancer datasets. In summary, our study revealed a novel function of FAM64A in promoting breast cancer stemness and EMT, suggesting that targeting of FAM64A may have therapeutic values in advanced breast cancer.

Transcriptome Alterations in Liver Metastases of Colorectal Cancer After Acquired Resistance to Cetuximab.

BACKGROUND/AIM: Cetuximab in combination with chemotherapy is recommended as first-line therapy for metastatic colorectal cancer (mCRC) with wild-type RAS. However, drug resistance to cetuximab exists widely in mCRC and reduces the prognosis of patients. Although some genomic alterations have been demonstrated to drive acquired resistance to cetuximab, the overall compendium of inherent molecular mechanisms is still incomplete. MATERIALS AND METHODS: Four liver metastasis biopsies were collected from two mCRC patients who were treated with cetuximab in combination with 5-fluororacil plus leucovorin and oxaliplatin (FOLFOX) regimen. RESULTS: Transcriptomic analysis revealed global gene expression alterations between paired samples prior to treatment and after acquired resistance. Further bioinformatics analysis discovered differentially expressed protein-coding genes/lncRNAs/miRNAs, potential miRNA-mRNA regulatory networks and lncRNA-mRNA competing endogenous RNA network, which may be potential biomarkers or play roles during the process of acquired resistance to cetuximab. CONCLUSION: Our study contributes to deciphering the molecular mechanisms of acquired resistance to cetuximab.

Landscape of B cell immunity and related immune evasion in human cancers.

Tumor-infiltrating B cells are an important component in the microenvironment but have unclear anti-tumor effects. We enhanced our previous computational algorithm TRUST to extract the B cell immunoglobulin hypervariable regions from bulk tumor RNA-sequencing data. TRUST assembled more than 30 million complementarity-determining region 3 sequences of the B cell heavy chain (IgH) from The Cancer Genome Atlas. Widespread B cell clonal expansions and immunoglobulin subclass switch events were observed in diverse human cancers. Prevalent somatic copy number alterations in the MICA and MICB genes related to antibody-dependent cell-mediated cytotoxicity were identified in tumors with elevated B cell activity. The IgG3-1 subclass switch interacts with B cell-receptor affinity maturation and defects in the antibody-dependent cell-mediated cytotoxicity pathway. Comprehensive pancancer analyses of tumor-infiltrating B cell-receptor repertoires identified novel tumor immune evasion mechanisms through genetic alterations. The IgH sequences identified here are potentially useful resources for future development of immunotherapies.

Expression of Bcl-2 is a favorable prognostic biomarker in lung squamous cell carcinoma.

Lung squamous cell carcinoma (LUSC) is the second major type of lung cancer globally. The majority of patients with LUSC are clinically diagnosed at the advanced stages, thus it is urgent to identify suitable prognostic markers for LUSC. B-cell lymphoma 2 (Bcl-2) has been widely studied in non-small cell lung cancer (NSCLC). However, the prognostic role of Bcl-2 in NSCLC remains conflicting and controversial, particularly for LUSC. Although certain studies have been performed to identify the prognostic value of Bcl-2, to the best of our knowledge, no study has investigated the prognostic role of Bcl-2 in LUSC specifically. The present study aimed to comprehensively evaluate the prognostic value of Bcl-2 in LUSC. Microarray data for LUSC were downloaded from public databases, including the Gene Expression Omnibus and The Cancer Genome Atlas. Microarray data of 901 patients with LUSC from 16 data sets were retrieved. The meta-z algorithm was applied and the combined z score was identified as -2.43, suggesting Bcl-2 is a favorable prognostic biomarker. Furthermore, immunohistochemical staining of Bcl-2 expression was performed in a tissue microarray of 72 patients with LUSC and survival analysis demonstrated that patients with high expression Bcl-2 exhibited significantly more improved overall survival rates compared with those with low Bcl-2 expression. Multivariate Cox regression revealed that high expression of Bcl-2 is an independent favorable prognostic factor (hazard ratio, 0.295; confidence interval, 0.097-0.904; P<0.05). Therefore, the results of the present study demonstrated that Bcl-2 is a favorable prognostic biomarker in LUSC.

Signaling protein signature predicts clinical outcome of non-small-cell lung cancer.

BACKGROUND: Non-small-cell lung cancer (NSCLC) is characterized by abnormalities of numerous signaling proteins that play pivotal roles in cancer development and progression. Many of these proteins have been reported to be correlated with clinical outcomes of NSCLC. However, none of them could provide adequate accuracy of prognosis prediction in clinical application. METHODS: A total of 384 resected NSCLC specimens from two hospitals in Beijing (BJ) and Chongqing (CQ) were collected. Using immunohistochemistry (IHC) staining on stored formalin-fixed paraffin-embedded (FFPE) surgical samples, we examined the expression levels of 75 critical proteins on BJ samples. Random forest algorithm (RFA) and support vector machines (SVM) computation were applied to identify protein signatures on 2/3 randomly assigned BJ samples. The identified signatures were tested on the remaining BJ samples, and were further validated with CQ independent cohort. RESULTS: A 6-protein signature for adenocarcinoma (ADC) and a 5-protein signature for squamous cell carcinoma (SCC) were identified from training sets and tested in testing sets. In independent validation with CQ cohort, patients can also be divided into high- and low-risk groups with significantly different median overall survivals by Kaplan-Meier analysis, both in ADC (31 months vs. 87 months, HR 2.81; P <  0.001) and SCC patients (27 months vs. not reached, HR 9.97; P <  0.001). Cox regression analysis showed that both signatures are independent prognostic indicators and outperformed TNM staging (ADC: adjusted HR 3.07 vs. 2.43, SCC: adjusted HR 7.84 vs. 2.24). Particularly, we found that only the ADC patients in high-risk group significantly benefited from adjuvant chemotherapy (P = 0.018). CONCLUSIONS: Both ADC and SCC protein signatures could effectively stratify the prognosis of NSCLC patients, and may support patient selection for adjuvant chemotherapy.

Analysis of PD1, PDL1, PDL2 expression and T cells infiltration in 1014 gastric cancer patients.

Although immune checkpoint blockade have demonstrated promising results, their effects on gastric cancer (GC) are under investigation. Understanding the clinical significance of PD1 and its ligands' expression, together with T cell infiltration might provide clues for biomarkers screening in GC immunotherapy. Immunohistochemistry were performed on a tissue microarray including 1,014 GC specimens using PD1, PDL1 and PDL2 antibodies. T cell markers CD3 and CD8 were also stained and quantified by automated image analysis. Correlation with clinical features and outcome were analyzed after controlling for potential confounders including EBV infection, HER2, C-met and PCNA expression. 37.8% of the cases showed membranous PD-L1 expression in tumor cells and 74.9% in infiltrating immune cells. PDL1 expression rate was rather higher in patients without metastasis, in EBV positive group and those with C-met and PCNA expression. GC patients with high level PDL1 expression exhibited better survival. GC Patients with higher T cell infiltration also showed elevated PDL1, PDL2 and PD1 expression and predict favorable outcome, indicating an adaptive immune resistance mechanism may exist. The group of patients infiltrated with lower density CD3+ T cells also without PDL1 expression in tumor cells predict the worst outcome in the subgroup of different PTNM stage, which may suggest an inactive immune status. These results highlights the need to assess both PDL1 expression in all tumor context and the characterization of the GC immune microenvironment.

Corrigendum: BCIP: A gene-centered platform for identifying potential regulatory genes in breast cancer.

This corrects the article DOI: 10.1038/srep45235.

Combining multi-dimensional data to identify key genes and pathways in gastric cancer.

Gastric cancer is an aggressive cancer that is often diagnosed late. Early detection and treatment require a better understanding of the molecular pathology of the disease. The present study combined data on gene expression and regulatory levels (microRNA, methylation, copy number) with the aim of identifying key genes and pathways for gastric cancer. Data used in this study was retrieved from The Cancer Genomic Atlas. Differential analyses between gastric cancer and normal tissues were carried out using Limma. Copy number alterations were identified for tumor samples. Bimodal filtering of differentially expressed genes (DEGs) based on regulatory changes was performed to identify candidate genes. Protein-protein interaction networks for candidate genes were generated by Cytoscape software. Gene ontology and pathway analyses were performed, and disease-associated network was constructed using the Agilent literature search plugin on Cytoscape. In total, we identified 3602 DEGs, 251 differentially expressed microRNAs, 604 differential methylation-sites, and 52 copy number altered regions. Three groups of candidate genes controlled by different regulatory mechanisms were screened out. Interaction networks for candidate genes were constructed consisting of 415, 228, and 233 genes, respectively, all of which were enriched in cell cycle, P53 signaling, DNA replication, viral carcinogenesis, HTLV-1 infection, and progesterone mediated oocyte maturation pathways. Nine hub genes (SRC, KAT2B, NR3C1, CDK6, MCM2, PRKDC, BLM, CCNE1, PARK2) were identified that were presumed to be key regulators of the networks; seven of these were shown to be implicated in gastric cancer through disease-associated network construction. The genes and pathways identified in our study may play pivotal roles in gastric carcinogenesis and have clinical significance.

A long non-coding RNA lncRNA-PE promotes invasion and epithelial-mesenchymal transition in hepatocellular carcinoma through the miR-200a/b-ZEB1 pathway.

Long non-coding RNAs have been revealed to play important roles in the progression of hepatocellular carcinoma. However, the detailed mechanisms underlying their activities are not fully understood. Using microarray technology, a number of long non-coding RNAs were previously identified to be aberrantly expressed in hepatocellular carcinoma. In this study, one of these long non-coding RNAs, designated lncRNA-PE (lncRNA promotes epithelial-mesenchymal transition), was further explored to study its expression profile and function. A cohort of human hepatocellular carcinoma tissue samples combined with benign controls and established human hepatocellular carcinoma cell lines were examined for the expression of lncRNA-PE. The biological functions of lncRNA-PE were examined by wound-healing and Transwell assays, which revealed that lncRNA-PE promotes cell invasion and migration. By detecting the level of epithelial-mesenchymal transition markers, lncRNA-PE was revealed to promote epithelial-mesenchymal transition in hepatocellular carcinoma cells. Further study suggested that lncRNA-PE downregulated miR-200a/b by repressing the primary transcript expression, enhanced ZEB1 expression, and promoted epithelial-mesenchymal transition of hepatocellular carcinoma cells. All these data imply that lncRNA-PE might play an important role in hepatocellular carcinoma development via the miR-200a/b-ZEB1 pathway.

Transcriptional response profiles of paired tumor-normal samples offer novel perspectives in pan-cancer analysis.

Both tumor and adjacent normal tissues are valuable in cancer research. Transcriptional response profiles represent the changes of gene expression levels between paired tumor and adjacent normal tissues. In this study, we performed a pan-cancer analysis based on the transcriptional response profiles from 633 samples across 13 cancer types. We obtained two interesting results. Using consensus clustering method, we characterized ten clusters with distinct transcriptional response patterns and enriched pathways. Notably, head and neck squamous cell carcinoma was divided in two subtypes, enriched in cell cycle-related pathways and cell adhesion-related pathways respectively. The other interesting result is that we identified 92 potential pan-cancer genes that were consistently upregulated across multiple cancer types. Knockdown of FAM64A or TROAP inhibited the growth of cancer cells, suggesting that these genes may promote tumor development and are worthy of further validations. Our results suggest that transcriptional response profiles of paired tumor-normal tissues can provide novel perspectives in pan-cancer analysis.

Quantitative proteomics by SWATH-MS reveals sophisticated metabolic reprogramming in hepatocellular carcinoma tissues.

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide, and understanding its molecular pathogenesis is pivotal to managing this disease. Sequential window acquisition of all theoretical mass spectra (SWATH-MS) is an optimal proteomic strategy to seek crucial proteins involved in HCC development and progression. In this study, a quantitative proteomic study of tumour and adjacent non-tumour liver tissues was performed using a SWATH-MS strategy. In total, 4,216 proteins were reliably quantified, and 338 were differentially expressed, with 191 proteins up-regulated and 147 down-regulated in HCC tissues compared with adjacent non-tumourous tissues. Functional analysis revealed distinct pathway enrichment of up- and down-regulated proteins. The most significantly down-regulated proteins were involved in metabolic pathways. Notably, our study revealed sophisticated metabolic reprogramming in HCC, including alteration of the pentose phosphate pathway; serine, glycine and sarcosine biosynthesis/metabolism; glycolysis; gluconeogenesis; fatty acid biosynthesis; and fatty acid ß-oxidation. Twenty-seven metabolic enzymes, including PCK2, PDH and G6PD, were significantly changed in this study. To our knowledge, this study presents the most complete view of tissue-specific metabolic reprogramming in HCC, identifying hundreds of differentially expressed proteins, which together form a rich resource for novel drug targets or diagnostic biomarker discovery.