Combined High-Throughput Proteomics and Random Forest Machine-Learning Approach Differentiates and Classifies Metabolic, Immune, Signaling and ECM Intra-Tumor Heterogeneity of Colorectal Cancer

Colorectal cancer (CRC) is a frequent, worldwide tumor described for its huge complexity, including inter-/intra-heterogeneity and tumor microenvironment (TME) variability. Intra-tumor heterogeneity and its connections with metabolic reprogramming and epithelial-mesenchymal transition (EMT) were investigated with explorative shotgun proteomics complemented by a Random Forest (RF) machine-learning approach. Deep and superficial tumor regions and distant-site non-tumor samples from the same patients (n = 16) were analyzed. Among the 2009 proteins analyzed, 91 proteins, including 23 novel potential CRC hallmarks, showed significant quantitative changes. In addition, a 98.4% accurate classification of the three analyzed tissues was obtained by RF using a set of 21 proteins. Subunit E1 of 2-oxoglutarate dehydrogenase (OGDH-E1) was the best classifying factor for the superficial tumor region, while sorting nexin-18 and coatomer-beta protein (beta-COP), implicated in protein trafficking, classified the deep region. Down- and up-regulations of metabolic checkpoints involved different proteins in superficial and deep tumors. Analogously to immune checkpoints affecting the TME, cytoskeleton and extracellular matrix (ECM) dynamics were crucial for EMT. Galectin-3, basigin, S100A9, and fibronectin involved in TME-CRC-ECM crosstalk were found to be differently variated in both tumor regions. Different metabolic strategies appeared to be adopted by the two CRC regions to uncouple the Krebs cycle and cytosolic glucose metabolism, promote lipogenesis, promote amino acid synthesis, down-regulate bioenergetics in mitochondria, and up-regulate oxidative stress. Finally, correlations with the Dukes stage and budding supported the finding of novel potential CRC hallmarks and therapeutic targets.

A top-down proteomic approach reveals a salivary protein profile able to classify Parkinson's disease with respect to Alzheimer's disease patients and to healthy controls

Parkinson's disease (PD) is a complex neurodegenerative disease with motor and non-motor symptoms. Diagnosis is complicated by lack of reliable biomarkers. To individuate peptides and/or proteins with diagnostic potential for early diagnosis, severity and discrimination from similar pathologies, the salivary proteome in 36 PD patients was investigated in comparison with 36 healthy controls (HC) and 35 Alzheimer's disease (AD) patients. A top-down platform based on HPLC-ESI-IT-MS allowed characterizing and quantifying intact peptides, small proteins and their PTMs (overall 51). The three groups showed significantly different protein profiles, PD showed the highest levels of cystatin SA and antileukoproteinase and the lowest of cystatin SN and some statherin proteoforms. HC exhibited the lowest abundance of thymosin β4, short S100A9, cystatin A, and dimeric cystatin B. AD patients showed the highest abundance of α-defensins and short oxidized S100A9. Moreover, different proteoforms of the same protein, as S-cysteinylated and S-glutathionylated cystatin B, showed opposite trends in the two pathological groups. Statherin, cystatins SA and SN classified accurately PD from HC and AD subjects. α-defensins, histatin 1, oxidized S100A9, and P-B fragments were the best classifying factors between PD and AD patients. Interestingly statherin and thymosin β4 correlated with defective olfactory functions in PD patients. All these outcomes highlighted implications of specific proteoforms involved in the innate-immune response and inflammation regulation at oral and systemic level, suggesting a possible panel of molecular and clinical markers suitable to recognize subjects affected by PD.

Combined Salivary Proteome Profiling and Machine Learning Analysis Provides Insight into Molecular Signature for Autoimmune Liver Diseases Classification

Autoimmune hepatitis (AIH) and primary biliary cholangitis (PBC) are autoimmune liver diseases that target the liver and have a wide spectrum of presentation. A global overview of quantitative variations on the salivary proteome in presence of these two pathologies is investigated in this study. The acid-insoluble salivary fraction of AIH and PBC patients, and healthy controls (HCs), was analyzed using a gel-based bottom-up proteomic approach combined with a robust machine learning statistical analysis of the dataset. The abundance of Arginase, Junction plakoglobin, Desmoplakin, Hexokinase-3 and Desmocollin-1 decreased, while that of BPI fold-containing family A member 2 increased in AIHp compared to HCs; the abundance of Gelsolin, CD14, Tumor-associated calcium signal transducer 2, Clusterin, Heterogeneous nuclear ribonucleoproteins A2/B1, Cofilin-1 and BPI fold-containing family B member 2 increased in PBCp compared to HCs. The abundance of Hornerin decreased in both AIHp and PBCp with respect to HCs and provided an area under the ROC curve of 0.939. Machine learning analysis confirmed the feasibility of the salivary proteome to discriminate groups of subjects based on AIH or PBC occurrence as previously suggested by our group. The topology-based functional enrichment analysis performed on these potential salivary biomarkers highlights an enrichment of terms mostly related to the immune system, but also with a strong involvement in liver fibrosis process and with antimicrobial activity.

Salivary Proteomics Reveals Significant Changes in Relation to Alzheimer's Disease and Aging

BACKGROUND

Aging is a risk factor for several pathologies as Alzheimer's disease (AD). Great interest exists, therefore, in discovering diagnostic biomarkers and indicators discriminating biological aging and health status. To this aim, omic investigations of biological matrices, as saliva, whose sampling is easy and non-invasive, offer great potential.

OBJECTIVE

Investigate the salivary proteome through a statistical comparison of the proteomic data by several approaches to highlight quali-/quantitative variations associated specifically either to aging or to AD occurrence, and, thus, able to classify the subjects.

METHODS

Salivary proteomic data of healthy controls under-70 (adults) and over-70 (elderly) years old, and over-70 AD patients, obtained by liquid chromatography/mass spectrometry, were analyzed by multiple Mann-Whitney test, Kendall correlation, and Random-Forest (RF) analysis.

RESULTS

Almost all the investigated proteins/peptides significantly decreased in relation to aging in elderly subjects, with or without AD, in comparison with adults. AD subjects exhibited the highest levels of α-defensins, thymosin β4, cystatin B, S100A8 and A9. Correlation tests also highlighted age/disease associated differences. RF analysis individuated quali-/quantitative variations in 20 components, as oxidized S100A8 and S100A9, α-defensin 3, P-B peptide, able to classify with great accuracy the subjects into the three groups.

CONCLUSION

The findings demonstrated a strong change of the salivary protein profile in relation to the aging. Potential biomarkers candidates of AD were individuated in peptides/proteins involved in antimicrobial defense, innate immune system, inflammation, and in oxidative stress. RF analysis revealed the feasibility of the salivary proteome to discriminate groups of subjects based on age and health status.

Predictive validity for the clinical diagnosis of a new parent questionnaire, the CABI, compared with CBCL

BACKGROUND

The questionnaires completed by the parents give a first general information on the behavioral problems of the child-adolescent, as a useful orientation to the clinical evaluation. The Child and Adolescent Behavior Inventory (CABI) is a 75-item parent questionnaire, which explores a large number of problem areas. The study of its predictive validity for the clinical diagnosis, in comparison with the Diagnostic and Statistical Manual of Mental Disorders (DSM)-oriented scales of the Child Behavior Checklist (CBCL), can assess whether its use may be advantageous.

MATERIAL AND METHODS

Parents/caregivers of 462 children and adolescents responded to both CABI and CBCL as a preliminary routine investigation. The results were compared with those of diagnoses obtained after the completion of the usual clinical procedure.

RESULTS

Accuracy values (probability of correct classification) resulted high for both instruments and significantly better for CABI anxiety and attention-deficit hyperactivity disorder (ADHD) scales, and for CBCL oppositional defiant disorder (ODD) and conduct disorder (CD) scales; no significant difference was found for depression scales. All the areas under the curve (AUC) of the receiver operating characteristic analysis reached excellent values, suggesting a very good predictive ability of the five scales of the two instruments. The comparison of AUC showed the CABI's anxiety and ADHD scales to give significantly higher values than those of CBCL, indicating that these two scales have a better predictive ability.

CONCLUSION

The study indicates a very good comparative (vs CBCL) and predictive validity of the CABI, suggesting an advantage in the use of this shorter questionnaire, available for free use both for clinical practice and supposedly for screening and epidemiological evaluations.

Molecular Signature and Mechanisms of Hepatitis D Virus-Associated Hepatocellular Carcinoma

There is limited data on the molecular mechanisms whereby hepatitis D virus (HDV) promotes liver cancer. Therefore, serum and liver specimens obtained at the time of liver transplantation from well-characterized patients with HDV-HCC (n = 5) and with non-HCC HDV cirrhosis (n = 7) were studied using an integrated genomic approach. Transcriptomic profiling was performed using laser capture-microdissected (LCM) malignant and nonmalignant hepatocytes, tumorous and nontumorous liver tissue from patients with HDV-HCC, and liver tissue from patients with non-HCC HDV cirrhosis. HDV-HCC was also compared with hepatitis B virus (HBV) HBV-HCC alone, and hepatitis C virus (HCV) HCV-HCC. HDV malignant hepatocytes were characterized by an enrichment of upregulated transcripts associated with pathways involved in cell-cycle/DNA replication, damage, and repair (Sonic Hedgehog, GADD45, DNA-damage-induced 14-3-3σ, cyclins and cell-cycle regulation, cell cycle: G₂-M DNA-damage checkpoint regulation, and hereditary breast cancer). Moreover, a large network of genes identified functionally relate to DNA repair, cell cycle, mitotic apparatus, and cell division, including 4 cancer testis antigen genes, attesting to the critical role of genetic instability in this tumor. Besides being overexpressed, these genes were also strongly coregulated. Gene coregulation was high not only when compared with nonmalignant hepatocytes, but also to malignant hepatocytes from HBV-HCC alone or HCV-HCC. Activation and coregulation of genes critically associated with DNA replication, damage, and repair point to genetic instability as an important mechanism of HDV hepatocarcinogenesis. This specific HDV-HCC trait emerged also from the comparison of the molecular pathways identified for each hepatitis virus-associated HCC. Despite the dependence of HDV on HBV, these findings suggest that HDV and HBV promote carcinogenesis by distinct molecular mechanisms.Implications: This study identifies a molecular signature of HDV-associated hepatocellular carcinoma and suggests the potential for new biomarkers for early diagnostics. Mol Cancer Res; 16(9); 1406-19. ©2018 AACR.

Infection with hepatitis C virus depends on TACSTD2, a regulator of claudin-1 and occludin highly downregulated in hepatocellular carcinoma

Entry of hepatitis C virus (HCV) into hepatocytes is a complex process that involves numerous cellular factors, including the scavenger receptor class B type 1 (SR-B1), the tetraspanin CD81, and the tight junction (TJ) proteins claudin-1 (CLDN1) and occludin (OCLN). Despite expression of all known HCV-entry factors, in vitro models based on hepatoma cell lines do not fully reproduce the in vivo susceptibility of liver cells to primary HCV isolates, implying the existence of additional host factors which are critical for HCV entry and/or replication. Likewise, HCV replication is severely impaired within hepatocellular carcinoma (HCC) tissue in vivo, but the mechanisms responsible for this restriction are presently unknown. Here, we identify tumor-associated calcium signal transducer 2 (TACSTD2), one of the most downregulated genes in primary HCC tissue, as a host factor that interacts with CLDN1 and OCLN and regulates their cellular localization. TACSTD2 gene silencing disrupts the typical linear distribution of CLDN1 and OCLN along the cellular membrane in both hepatoma cells and primary human hepatocytes, recapitulating the pattern observed in vivo in primary HCC tissue. Mechanistic studies suggest that TACSTD2 is involved in the phosphorylation of CLDN1 and OCLN, which is required for their proper cellular localization. Silencing of TACSTD2 dramatically inhibits HCV infection with a pan-genotype effect that occurs at the level of viral entry. Our study identifies TACSTD2 as a novel regulator of two major HCV-entry factors, CLDN1 and OCLN, which is strongly downregulated in malignant hepatocytes. These results provide new insights into the complex process of HCV entry into hepatocytes and may assist in the development of more efficient cellular systems for HCV propagation in vitro.

Contour recognition of complex leaf shapes

The leaf shape is an important taxonomic character. Compared to the classic morphological leaf features such as veins, margin indentations, sinuses, etc., the shape is simpler to obtain by using the 'magic wand' or other contouring tools that are available in most of imaging applications. The only exception is when leaves develop large lobes that get in touch or overlap each other, as the presence of hidden or closed portions of the leaf border precludes the application of automatic methods and forces the leaf contour to be traced manually. This is a time consuming and relatively accurate operation that, nevertheless, can not be avoided, as overlapping lobes are characteristic features of the leaves of several plant species and varieties. The method described in the paper overcomes this problem as it allows the leaf contour to be achieved even in the presence of touching or overlapping lobes. The method involves three steps: (1) the acquisition of leaf images using a transilluminator, (2) a two-level image segmentation that allows all leaf components (blade, overlapping lobes and closed sinuses) to be represented in a single binary image, and (3) the contouring and concatenation of all binary outlines in a single, self-intersecting closed contour that reproduces accurately the leaf shape. The method can be extended to acquire the shape of leaves of herbarium specimens, that are often overlapped but can not be easily handled and repositioned because of their extreme fragility and relevant taxonomic value.

Vegetative propagation of Plukenetia polyadenia by cuttings: effects of leaf area and indole-3-butyric acid concentration

The seeds of Plukenetia polyadenia have high levels of unsaturated fatty acids and are used as medicine and food for native people in the Peruvian and Brazilian Amazon. The objective of this study was to develop a method for vegetative propagation of Plukenetia polyadenia by rooting of cuttings. The experiment was laid out in a randomized complete block design with 12 treatments and 3 replications of 8 cuttings, in a 3 × 4 factorial arrangement. The factors were: 3 levels of leaf area (25, 50 and 75%) and 3 indole-3-butyric acid - IBA concentrations (9.84, 19.68 and 29.52mM) and a control without IBA. Data were submitted to analysis of variance and means were compared by Tukey test at 5% probability. Our results show that the use of cuttings with 50% of leaf area and treatment with 29.52mM of IBA induced high percentages of rooting (93%) and the best root formation. Vegetative propagation of Plukenetia polyadenia by cuttings will be used as a tool to conserve and propagate germplasm in breeding programs.

Integrated ordination of miRNA and mRNA expression profiles

Background

Several studies have investigated miRNA and mRNA co-expression to identify regulatory networks at the transcriptional level. A typical finding of these studies is the presence of both negative and positive miRNA-mRNA correlations. Negative correlations are consistent with the expected, faster degradation of target mRNAs, whereas positive correlations denote the existence of feed-forward regulations mediated by transcription factors. Both mechanisms have been characterized at the molecular level, although comprehensive methods to represent miRNA-mRNA correlations are lacking. At present, genome-wide studies are able to assess the expression of more than 1000 mature miRNAs and more than 35,000 well-characterized human genes. Even if studies are generally restricted to a small subset of genes differentially expressed in specific diseases or experimental conditions, the number of potential correlations remains very high, and needs robust multivariate methods to be conveniently summarized by a small set of data.

Results

Nonparametric Kendall correlations were calculated between miRNAs and mRNAs differentially expressed in livers of patients with acute liver failure (ALF) using normal livers as controls. Spurious correlations due to the histopathological composition of samples were removed by partial correlations. Correlations were then transformed into distances and processed by multidimensional scaling (MDS) to map the miRNA and mRNA relationships. These showed: (a) a prominent displacement of miRNA and mRNA clusters in ALF livers, as compared to control livers, indicative of gene expression dysregulation; (b) a clustering of mRNAs consistent with their functional annotations [CYP450, transcription factors, complement, proliferation, HLA class II, monocytes/macrophages, T cells, T-NK cells and B cells], as well as a clustering of miRNAs with the same seed sequence; and (c) a tendency of miRNAs and mRNAs to populate distinct regions of the MDS plot. MDS also allowed to visualize the network of miRNA-mRNA target pairs.

Conclusions

Different features of miRNA and mRNA relationships can be represented as thematic maps within the framework of MDS obtained from pairwise correlations. The symmetric distribution of positive and negative correlations between miRNA and mRNA expression suggests that miRNAs are involved in a complex bidirectional molecular network, including, but not limited to, the inhibitory regulation of miRNA targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1971-9) contains supplementary material, which is available to authorized users.

Primary Graft Dysfunction After Zero-Mismatch Kidney Transplantation Secondary to Early Biopsy-Proven Acute Cell-Mediated Rejection: Case Report

We report a case of primary renal allograft dysfunction and early acute cell-mediated rejection after a 12/12 HLA antigen zero-mismatch (0MM) transplant. The recipient was a 40-year-old white man who was highly allosensitized, with a calculated panel reactive antibody score of 100%. In posteroperative day 1 the recipient remained anuric and underwent dialysis because of hyperkalemia. Graft biopsy showed early acute cellular rejection, Banff grade 2B. No evidence of antibody-mediated rejection was observed. To our knowledge, this case is the 1st to report early cell-mediated rejection after 12/12 HLA antigen 0MM kidney transplantation. This case suggests that highly sensitized candidates are at high immunologic risk even in the context of 0MM kidney transplantation.

A Landscape Analysis to Understand Orientation of Honey Bee (Hymenoptera: Apidae) Drones in Puerto Rico

Honey bees [Apis mellifera L. (Apidae, Hymenoptera)] show spatial learning behavior or orientation, in which animals make use of structured home ranges for their daily activities. Worker (female) orientation has been studied more extensively than drone (male) orientation. Given the extensive and large flight range of drones as part of their reproductive biology, the study of drone orientation may provide new insight on landscape features important for orientation. We report the return rate and orientation of drones released at three distances (1, 2, and 4 km) and at the four cardinal points from an apiary located in Gurabo, Puerto Rico. We used high-resolution aerial photographs to describe landscape characteristics at the releasing sites and at the apiary. Analyses of variance were used to test significance among returning times from different distances and directions. A principal components analysis was used to describe the landscape at the releasing sites and generalized linear models were used to identify landscape characteristics that influenced the returning times of drones. Our results showed for the first time that drones are able to return from as far as 4 km from the colony. Distance to drone congregation area, orientation, and tree lines were the most important landscape characteristics influencing drone return rate. We discuss the role of landscape in drone orientation.

Association of G1359A polymorphism of the cannabinoid receptor gene (CNR1) with macronutrient intakes in obese females

BACKGROUND

The endogenous cannabinoid system plays a role in metabolic aspects of body weight and feeding behaviour. A polymorphism (1359 G/A) (rs1049353) of the CB1 gene was reported as a common polymorphism in the Caucasian population. The present study aimed to investigate the association of the polymorphism (G1359A) of the CB1 receptor gene on macronutrient intake in females with obesity.

METHODS

A sample of 896 females was analysed. A bioimpedance measurement, a blood pressure measurement, a serial assessment of nutritional intake with 3 days of written food records, and a biochemical analysis were all performed. The genotype of the CNR1 receptor gene polymorphism (rs1049353) was studied.

RESULTS

Five hundred and sixteen patients (57.6%) had the genotype G1359G (non-A carriers) and 380 (42.4%) patients had G1359A (328 patients, 36.6%) or A1359A (52 patients, 5.8%) (A carriers). Triglycerides and high-density lipoprotein (HDL) cholesterol levels were higher in A non-A allele carriers than non-A allele carriers. The intakes of dietary cholesterol and saturated fat for the upper tertile (T3) compared to the baseline tertile were inversely associated with the CB1-R 1359 G/A polymorphism [odds ratio (OR) = 0.59; 95% confidence interval (CI) = 0.30-0.92 and OR = 0.66; 95% CI = 0.39-0.91, respectively]. These data were observed in the second tertile (T2) (OR = 0.61; 95% CI = 0.29-0.94 and OR = 0.58; 95% CI = 0.31-0.90, respectively).

CONCLUSIONS

The present study reports an association of the A allele with a better lipid profile (triglycerides and HDL cholesterol) than non-A allele carriers. In addition, this polymorphism is associated with a specific macronutrient intake, as well as with low cholesterol and fat saturated intakes.

Viral expression and molecular profiling in liver tissue versus microdissected hepatocytes in hepatitis B virus-associated hepatocellular carcinoma

Background

The molecular mechanisms whereby hepatitis B virus (HBV) induces hepatocellular carcinoma (HCC) remain elusive. We used genomic and molecular techniques to investigate host-virus interactions by studying multiple areas of the same liver from patients with HCC.

Methods

We compared the gene signature of whole liver tissue (WLT) versus laser capture-microdissected (LCM) hepatocytes along with the intrahepatic expression of HBV. Gene expression profiling was performed on up to 17 WLT specimens obtained at various distances from the tumor center from individual livers of 11 patients with HCC and on selected LCM samples. HBV markers in liver and serum were determined by real-time polymerase chain reaction (PCR) and confocal immunofluorescence.

Results

Analysis of 5 areas of the liver showed a sharp change in gene expression between the immediate perilesional area and tumor periphery that correlated with a significant decrease in the intrahepatic expression of HB surface antigen (HBsAg). The tumor was characterized by a large preponderance of down-regulated genes, mostly involved in the metabolism of lipids and fatty acids, glucose, amino acids and drugs, with down-regulation of pathways involved in the activation of PXR/RXR and PPARα/RXRα nuclear receptors, comprising PGC-1α and FOXO1, two key regulators critically involved not only in the metabolic functions of the liver but also in the life cycle of HBV, acting as essential transcription factors for viral gene expression. These findings were confirmed by gene expression of microdissected hepatocytes. Moreover, LCM of malignant hepatocytes also revealed up-regulation of unique genes associated with cancer and signaling pathways, including two novel HCC-associated cancer testis antigen genes, NUF2 and TTK.

Conclusions

Integrated gene expression profiling of whole liver tissue with that of microdissected hepatocytes demonstrated that HBV-associated HCC is characterized by a metabolism switch-off and by a significant reduction in HBsAg. LCM proved to be a critical tool to validate gene signatures associated with HCC and to identify genes that may play a role in hepatocarcinogenesis, opening new perspectives for the discovery of novel diagnostic markers and therapeutic targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-014-0230-1) contains supplementary material, which is available to authorized users.

Identification of microRNAs specifically expressed in hepatitis C virus-associated hepatocellular carcinoma

Although several studies have investigated the association of miRNAs with hepatocellular carcinoma (HCC), the data published so far are not concordant. A reason for these discrepancies may be the fact that most studies used the nontumorous tissue surrounding the HCC lesion as a control, which is almost invariably affected by cirrhosis or chronic hepatitis, as well as other pathological conditions such as hepatitis B virus (HBV) or hepatitis C virus (HCV) infection. Moreover, HCC is often analyzed as a single group regardless of the different viral etiologies. The miRNAs differentially expressed in HCV-related HCC were investigated by comparing the tumorous tissues to a wide range of liver specimens, including healthy livers obtained from liver donors and patients who underwent liver resection for angioma, in addition to tissues from various acute and chronic liver diseases, including HCV-related cirrhosis not associated with HCC, HCV-related cirrhosis associated with HCC and HBV-associated acute liver failure. The whole set of 2,226 human miRNAs were examined, including 1,121 pre-miRNAs and 1,105 mature miRNAs, available in a microarray platform. Stringent statistical methods were applied to reduce the risk of false discoveries to less than 1%. These data identified 18 miRNAs exclusively expressed in HCV-associated HCC, characterized by high specificity and selectivity versus all other liver diseases and healthy conditions and connected into a regulatory network pivoting on p53, phosphatase and tensin homolog and all-trans retinoic acid signaling.

Decompressive craniectomy and brain death prevalence and mortality: 8-year retrospective review

Decompressive craniectomy (DC) is a surgical practice that has been used since the late 19th century. The cerebral blood flow increase after the performance of a DC can delay and even prevent the development of cerebral circulatory arrest and brain death (BD). We aimed to determine the prevalence of BD, the use of DC, and the evolution to BD with versus without DC. This retrospective, observational, cross-sectional study was performed in a single high-intensity center in Argentina from January 2003 to December 2010. Inclusion criteria were all patients with Glasgow Coma Score of at most 7 on admission or during their stay in the intensive care units. Exclusion criteria were patients with incomplete data. In cases of death, we assessed whether they fulfilled BD criteria or if the cause of death was a cardiac arrest (CA). The 698 patients considered for analysis showed a 60% (n = 418) global mortality rate. The causes were: CA (n = 270); BD (n = 108) and others considered to be "undefined," namely not assessed completely for the diagnosis of BD (n = 40). According to diagnosis category, traumatic brain injury (TBI) was largest (nearly 50%). The DC group (n = 206) showed significant differences regarding sex and diagnosis category versus no DC group. Mortality was significantly lower in this group (48% versus 65%, P < .001). No significant differences were observed comparing causes of death (CA, BD, or undefined). The use of DC did not influence the frequency of BD development (24% versus 26%, P = .72). The average DC rate was 30% and of BD 16%. The prevalence of DC and better survival were recorded compared with subjects without DC. The prevalence of BD was lower than expected in accordance with national registries; however, among our group, DC did not seem to modify the evolution to BD.

Liver regeneration signature in hepatitis B virus (HBV)-associated acute liver failure identified by gene expression profiling

Introduction

The liver has inherent regenerative capacity via mitotic division of mature hepatocytes or, when the hepatic loss is massive or hepatocyte proliferation is impaired, through activation of hepatic stem/progenitor cells (HSPC). The dramatic clinical course of acute liver failure (ALF) has posed major limitations to investigating the molecular mechanisms of liver regeneration and the role of HSPC in this setting. We investigated the molecular mechanisms of liver regeneration in 4 patients who underwent liver transplantation for hepatitis B virus (HBV)-associated ALF.

Methods and Findings

Gene expression profiling of 17 liver specimens from the 4 ALF cases and individual specimens from 10 liver donors documented a distinct gene signature for ALF. However, unsupervised multidimensional scaling and hierarchical clustering identified two clusters of ALF that segregated according to histopathological severity massive hepatic necrosis (MHN; 2 patients) and submassive hepatic necrosis (SHN; 2 patients). We found that ALF is characterized by a strong HSPC gene signature, along with ductular reaction, both of which are more prominent in MHN. Interestingly, no evidence of further lineage differentiation was seen in MHN, whereas in SHN we detected cells with hepatocyte-like morphology. Strikingly, ALF was associated with a strong tumorigenesis gene signature. MHN had the greatest upregulation of stem cell genes (EpCAM, CK19, CK7), whereas the most up-regulated genes in SHN were related to cellular growth and proliferation. The extent of liver necrosis correlated with an overriding fibrogenesis gene signature, reflecting the wound-healing process.

Conclusion

Our data provide evidence for a distinct gene signature in HBV-associated ALF whose intensity is directly correlated with the histopathological severity. HSPC activation and fibrogenesis positively correlated with the extent of liver necrosis. Moreover, we detected a tumorigenesis gene signature in ALF, emphasizing the close relationship between liver regeneration and liver cancer.