- Research
- Open access
- Published:
Integrated transcriptomics and histopathology approach identifies a subset of rejected donor livers with potential suitability for transplantation
BMC Genomics volume 25, Article number: 437 (2024)
Abstract
Background
Liver transplantation is an effective treatment for liver failure. There is a large unmet demand, even as not all donated livers are transplanted. The clinical selection criteria for donor livers based on histopathological evaluation and liver function tests are variable. We integrated transcriptomics and histopathology to characterize donor liver biopsies obtained at the time of organ recovery. We performed RNA sequencing as well as manual and artificial intelligence-based histopathology (10 accepted and 21 rejected for transplantation).
Results
We identified two transcriptomically distinct rejected subsets (termed rejected-1 and rejected-2), where rejected-2 exhibited a near-complete transcriptomic overlap with the accepted livers, suggesting acceptability from a molecular standpoint. Liver metabolic functional genes were similarly upregulated, and extracellular matrix genes were similarly downregulated in the accepted and rejected-2 groups compared to rejected-1. The transcriptomic pattern of the rejected-2 subset was enriched for a gene expression signature of graft success post-transplantation. Serum AST, ALT, and total bilirubin levels showed similar overlapping patterns. Additional histopathological filtering identified cases with borderline scores and extensive molecular overlap with accepted donor livers.
Conclusions
Our integrated approach identified a subset of rejected donor livers that are likely suitable for transplantation, demonstrating the potential to expand the pool of transplantable livers.
Introduction
Liver transplantation is a crucial intervention for acute [1,2,3] and chronic liver failure [4], addressing a range of liver disease conditions [5,6,7,8,9,10]. While there has been a significant improvement in post-transplant survival rates, there is an increasing demand for donor livers available for transplants. Recent metrics from the Organ Procurement and Transplantation Network (OPTN) in the United States indicate an 11.18% increase in the number of deceased donors, reaching 8,788 donors compared to the previous year. However, there remains an unmet demand for liver transplantation because not all donated livers are transplanted (7.2% are not transplanted per OPTN [11] ). Additionally, 12,053 recipients were added to the waitlist for liver transplantation in the previous year (as of October 27, 2023 [11] ).
Current clinical selection criteria for donor livers include liver function tests to evaluate the levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST) [12, 13], total bilirubin, blood type compatibility between the donor and recipient [14, 15], absence of viral infections such as Hepatitis and HIV [16, 17], assessment of alcohol and drug-induced liver damage and histopathological evaluation of liver biopsy for steatosis, fibrosis, and necrosis [18, 19]. These comprehensive criteria aid in determining the suitability of donor livers for transplantation [20]. However, variability in the application of the clinical selection criteria between transplant centers results in the rejection of a fraction of donor livers for transplantation [21]. Furthermore, rejected donor livers are highly variable, with a range of differences in several parameters, including tissue steatosis, fibrosis, donor age, and underlying conditions. The wide range of variability among rejected donor livers could potentially lead to the rejection of a fraction of rejected donor livers that might be suitable for transplantation. For instance, livers from donation after circulatory death (DCD) show high reperfusion injury in contrast to livers from donation after brain death (DBD), which reduces graft success post-transplantation [22]. There is variability in the perfusion time between DBD and DCD livers [23], and the difference in the perfusion time influences the downstream functions impacting graft success [24, 25].
To address the variability between transplant centers and within the pool of rejected donor livers, a comprehensive characterization of donor livers beyond histopathology is warranted. Transcriptomics approach provides a global view of the underlying molecular processes and enables deeper insights into the pathophysiological mechanisms of liver function and regeneration. Previous studies have focused on characterizing the transcriptomic profiles of the liver to predict the post-transplantation outcomes, including the effects of ischemia-reperfusion injury [26] and initial graft function [27]. However, to date, no study has focused on the transcriptomic analysis of the donor livers for assessing the similarities and differences between the accepted and rejected livers, which can yield new information for augmenting the current criteria to evaluate the suitability for transplantation.
Here, we applied an integrated molecular and histopathological approach to characterize donor livers that were either accepted or rejected for transplantation. We evaluated the highly variable rejected group for differences and overlap with the accepted donor livers. Our analysis aimed to identify a subset, if any, of rejected donor livers that shares molecular and histopathological features with accepted donor livers.
Materials and methods
Ethical approval and sample collection
All research was conducted in accordance with both the Declarations of Helsinki and Istanbul. Informed Consent to Research was obtained by the Gift of Life Foundation to collect liver biopsy samples from the deceased donors. As per the regulations outlined in 45 CFR 46.101 (b) [4], the current study falls under Exemption 4, and does not require additional Institutional Review Board approval. A Certification for Protected Health Information of Decedents for Research was obtained from the Privacy Officer of the Office of University Counsel and Corporate Compliance Department of Thomas Jefferson University Hospital (TJUH). The informed consent was obtained by the Gift of Life Foundation from the deceased donor`s next of kin to collect liver biopsy samples for research purposes. Clinical information regarding patient demographics, including age, sex, race, body mass index, AST, ALT, and total bilirubin levels was collected (Table 1). Tissue samples were categorized into two groups: (1) the accepted group (n = 10), corresponding to donor livers accepted for transplantation by the TJUH Liver Transplant Center, and (2) the rejected group (n = 22), consisting of donor livers rejected for transplantation by multiple transplant centers (Table S1). Liver wedge biopsy tissue samples were obtained at the time of organ recovery. A portion of the tissue was flash-frozen and another was fixed in paraformaldehyde.
RNA extraction and sequencing
Total RNA was extracted from flash-frozen liver tissue samples using the miRNeasy Mini Kit (217,004; Qiagen, Hilden, Germany), and the extracted RNA was stored at -80 °C. The purity and quality of the RNA were assessed using a Bioanalyzer system (Agilent Technologies, California, United States). RNA sequencing was performed at Thomas Jefferson University Genomics Core Facility. Libraries were prepared using the Illumina Stranded Total RNA Prep with Ribo-Zero (20,040,525; Illumina, California, United States) and sequenced on a NovaSeq 6000 platform (Illumina). Paired-end sequencing was performed with read lengths of 2 × 100 bp, and each sample had a sequencing depth of at least 50 million reads. Fastq generation was performed using the Basespace cloud-based analysis platform (Illumina).
Data processing and quality control
We used the nf-core RNA-seq pipeline (version 3.7; doi: https://doi.org/10.5281/zenodo.1400710) within the Nextflow framework (version 22.04.0) [28]. This pipeline facilitated fastq preprocessing, quality assessment, transcriptome alignment using STAR aligner [29], and transcript quantification using salmon [30]. Alignments were performed against the GRCh38 genome using Gencode v40/ENSEMBL gene annotations [31]. One sample was an outlier with total gene counts less than two standard deviations below the mean across all 32 samples and was excluded from the analysis. The RNAseq data files have been deposited in the NCBI Gene Expression Omnibus under the accession number GSE243887. The normalized gene expression values are available in Table S2.
Differential gene expression analysis
The raw gene counts were normalized to counts per million (CPM), and genes with the sum of expression values across all samples less than 10 CPM were filtered out. The DESeq2 package (version 1.38.2) was used for differential gene expression analysis [32]. The raw filtered count data was normalized using the variance stabilizing transformation (vst) function from the DESeq2 package. Genes with two-fold up or down-regulation and padj < 0.01 were selected for downstream functional pathway enrichment analysis. Dimension reduction was performed using the umap R package [33]. K-means clustering with k = 2 was performed on the vst normalized data for accepted and rejected groups.
Estimating sensitivity of the statistical test
We performed a power analysis using the RNASeqPower package in R to assess the statistical confidence of our study based on the study sample size [34]. We estimated sensitivity across various comparison groups: accepted vs. rejected, accepted vs. rejected-1, rejected-2 vs. rejected-1, and accepted vs. rejected-2 groups.
Pathway enrichment and network analysis
Functional annotation and gene ontology analysis were conducted using the DAVID tool version 2021 [35]. Protein-protein interaction (PPI) networks were constructed using the Cytoscape String App [36], and highly interconnected regions within the network were identified using the Molecular Complex Detection (MCODE) app [37].
Comparison with transcriptomics of graft failure post-transplantation
We used the results from an earlier study that developed a set of classifier genes (78 genes) that are differentially regulated in livers displaying initial poor graft function (IPGF) post-perfusion and transplantation [27]. There were 69 downregulated and 7 upregulated genes in the non-IPGF vs. IPGF group (Table 5 of Defamie et al., 2008 [27]). We performed Gene Set Enrichment Analysis using fgsea package in R to evaluate whether the non-IPGF downregulated genes were enriched in the present differential gene expression comparisons (rejected-2 vs. rejected-1 and accepted vs. rejected-1). Differentially expressed genes in our study were ranked in decreasing order of log2 fold change. Enrichment scores were calculated based on the distribution of non-IPGF classifier genes in the rank-ordered differentially expressed gene list. Since the non-IPGF upregulated genes had only 4 genes, it was not statistically meaningful to evaluate enrichment.
Histopathological staining
We preserved 26 of 31 donor liver biopsies for independent histopathological evaluation (n = 7 accepted, n = 7 rejected-1, n = 12 rejected-2). Liver biopsy tissue samples (∼2 mm diameter) were locked in Tissue Embedding Cassettes and fixed in 4% paraformaldehyde (PFA) (Electron Microscopy Sciences, Hatfield, PA, United States) for 24 h on a shaker. The cassettes with the biopsy samples were washed three times in distilled water and then placed in 70% ethanol (Decon Labs, Pennsylvania, United States) at 4 °C until paraffin embedding. Histopathology processing was performed by the Pathology Core at Thomas Jefferson University. The sections were deparaffinized in xylene (534,056; Sigma Aldrich, Massachusetts, United States) and rehydrated in a series of decreasing ethanol concentrations, followed by distilled water. Serial sections were stained with Hematoxylin & Eosin (H&E) and Masson’s Trichrome.
H&E staining: Liver tissue sections were first incubated in Haematoxylin (SL100; StatLab, Texas, United States) for 2 min at room temperature and then washed under running cold water for 1 min. The sections were next incubated in acid water (SL103; StatLab) for 1 min, washed under running cold water for 1 min, and then incubated in Bluing Reagent (SL102; StatLab) for 1 min, followed by two washes in deionized water. Slides were then incubated in 95% Alcohol for 1 min, followed by incubation in Eosin (SL104; StatLab) for 2 min.
Trichrome staining: The sections were incubated in Mordant in Bouin’s Fixative (s129; Poly Scientific R&D Corp, New York, United States) for 1 h at room temperature. The sections were first gently rinsed in running tap water for 10–15 min and then rinsed once in distilled water. Next, the sections were incubated in Weigert’s Iron Hematoxylin Working Solution (s216b; Poly Scientific R&D Corp) for 10 min, gently washed in warm water for 10–15 min and then rinsed in distilled water (10 dips) for three changes of distilled water. The sections were incubated in Biebrich Scarlet – Acid Fuchsin (s125; Poly Scientific R&D Corp) for 5 min and rinsed twice in distilled water. The sections were then differentiated in Phosphotungstic – Phosphomolybdic Acid (s225; Poly Scientific R&D Corp.) for 12 min and then placed in Aniline Blue Solution (s116; Poly Scientific R&D Corp) for 20 min. The sections were rinsed twice with distilled water for 15 s each. Next, the sections were incubated in 1% Aqueous Acetic Acid (s100; Poly Scientific R&D Corp) for 3 min. Lastly, the sections were dehydrated in increasing ethanol concentrations, followed by xylene, and then coverslipped.
Histopathological scoring for liver steatosis, fibrosis, and necrosis
Whole slide images (WSIs) were acquired at ×20 magnification using an Aperio ImageScope digital slide scanner (Leica Biosystems, Wetzlar, Germany). Histopathological evaluation of steatosis, fibrosis, necrosis, and hepatocyte ballooning was performed by a gastrointestinal pathologist. In parallel, the percentage of steatosis and tissue collagen were assessed using pre-trained AI-based image analysis applications in Visiopharm software (Visiopharm Corporation, Colorado, United States). The areas of collagenous and steatotic tissues were divided by the total tissue area to calculate the respective percentages shown in Table 2.
Results
Identification of molecularly defined human donor liver subsets
We applied an integrated approach, including evaluation of clinical data, histopathology, and transcriptomics, to identify a subset of rejected donor livers that may be suitable for transplantation (Fig. 1A). We obtained liver biopsies from deceased donors at the time of organ recovery, either accepted or rejected for liver transplantation. Information on donor demographics and standard serological tests is provided in Table 1. We performed RNA sequencing of donor liver biopsies to characterize the global gene expression profiles and enriched functional pathways. In parallel, a histopathological evaluation was performed to assess steatosis, fibrosis and necrosis.
To characterize the molecular differences between the accepted and rejected donor livers, as well as the variability within the rejected group, we performed dimension reduction using Uniform Manifold Approximation and Projection (UMAP) on the normalized transcriptomics data on the expression of 18,033 genes. The accepted and rejected groups were not distinctly separated within the broader molecular heterogeneity of the donor livers (Fig. 1B). We examined the expression of select genes critical to liver function and homeostasis, including albumin (ALB), glutamate pyruvate transaminase gene (GPT) encoding ALT, glutamic-oxaloacetic transaminase gene (GOT) encoding AST, cytochrome P450 Family 2 Subfamily E Member 1 (CYP2E1), and solute Carrier Family 2 Member 2 (SLC2A2). We observed an overlapping gene expression between the accepted and rejected groups (Fig. 1C). We found a high variability in gene expression in the rejected group, with a fraction of the rejected donor livers and accepted livers showing similar expression levels of key genes (Fig. 1C).
Molecular profiling indicates extensive gene expression overlap between accepted and rejected donor livers
Differential gene expression analysis between the accepted and rejected liver donor groups yielded only 55 of the 18,033 genes as statistically significant (padj<0.01, estimated sensitivity of 87.95%, Fig. 2A). The corresponding list of differential gene expression fold changes is included in Supplementary Table S3. Visualization of the expression profiles of these 55 genes suggests that a subset of rejected donor livers exhibit similar gene expression profiles to that of the accepted group. Dimension reduction analysis and classification based on these 55 genes identified two molecularly distinguishable subsets within rejected donor livers (rejected-1 and rejected-2 groups hereon) (Fig. 2B). Of the 21 rejected donor livers, 9 were classified as rejected-1 and 12 as rejected-2. The rejected-2 subset clustered closely with the accepted group, while the rejected-1 subset was further away from the accepted group. Analysis of the selected genes involved in drug metabolism (CYP2C19), fatty acid metabolism (FABP5, ALDH1L2), cell adhesion (SDCBP2), cytoskeletal organization (ITGA2), inflammation, and hepatic stellate cell activation (S100A6) illustrated the extensive molecular overlap between the accepted and rejected-2 groups (Fig. 2C).
A subset of rejected donor livers share broad transcriptomic features with the accepted group
We determined the significantly differentially expressed genes (DEGs) between three groups: accepted vs. rejected-1, accepted vs. rejected-2, and rejected-2 vs. rejected-1 (padj < 0.01, two-fold up or down). Our analysis revealed a total of 1942 DEGs between the accepted and rejected-1 groups (765 upregulated genes and 1177 downregulated genes, estimated sensitivity of 79.62%), and 1821 DEGs between the rejected-2 and rejected-1 groups (751 upregulated genes and 1070 downregulated genes, estimated sensitivity of 83.42%, Fig. 3A). There were 2 DEGs between the accepted and rejected-2 donor livers (2 upregulated genes and no downregulated genes, estimated sensitivity 93.46%), indicating extensive and transcriptome-wide molecular overlap between the accepted and the rejected-2 subset. The corresponding lists of differential gene expression fold changes for the three pairwise comparisons are included in Supplementary Table S4. Dimension reduction analysis based on the DEGs in either the accepted and rejected-2 groups compared to the rejected-1 group revealed the rejected-1 group was separate from the other two groups (accepted and rejected-1). There was an extensive overlap between the accepted and the rejected 2 groups (Fig. 3B). A comparison of the differential gene expression profiles in the accepted and rejected-2 groups suggests that the majority of the genes show similar expression patterns (Fig. 3C). The two rejected subsets had similar proportions of DBD and DCD cases (Table 1 - rejected-2: 58% DCD; rejected-1: 55% DCD). Overall, our findings suggest that the rejected-2 subset within the highly variable rejected group shared many molecular features with the accepted donor livers, indicating potential similarities in their underlying molecular states and, by implication, in their functional capacities in multiple pathways critical to liver pathophysiology.
Shared molecular features between accepted and a subset of rejected donor livers contribute to metabolic functions and extracellular matrix organization
Pathway enrichment analysis revealed that the upregulated genes in accepted and rejected-2 were mainly enriched for lipid metabolism, fatty acid metabolism, cholesterol biosynthesis, cholesterol homeostasis, electron transport chain, bile acid biosynthesis, and transmembrane transport (False Discovery Rate < 5%) (Fig. 4A). The downregulated genes in accepted and rejected-2 were mainly enriched for extracellular liver matrix organization, integrin-mediated signaling, collagen biosynthesis, collagen fibril organization, cell adhesion, cell migration, and cell differentiation (Fig. 4B).
Our analysis thus far is based on comparing the rejected-2 group with the accepted donor livers for similarity of transcriptomics, suggesting similar functional capacities. As an additional comparative analysis for the suitability of transplantation of the rejected-2 group, we considered a published data set on transcriptomically derived classifier genes for predicting successful grant function post-transplantation [27]. The classifier consists of 69 down-regulated genes and 9 up-regulated genes in non-IPGF vs. IPGF. We tested the hypothesis that genes with an expression signature predictive of post-transplantation success are enriched in the rejected-2 subset relative to the rejected-1 group. We found that non-IPGF vs. IPGF classifier genes showed a robust signature in our transcriptomics data, and out of the 69 downregulated genes in the non-IPGF group that showed successful graft function post-transplantation, 47 genes were also downregulated in the rejected-2 and accepted groups compared to the rejected-1 group in our study (Supplementary Fig. 1A, B).
We constructed protein-protein interaction networks for the genes overlapping between the accepted and rejected-2 groups. The most interconnected module of the upregulated genes contained 14 proteins (nodes) and 71 interactions (edges) corresponding to fatty acid metabolism, lipid metabolism, amino acid metabolism, and glucose transport (Fig. 4C). The most interconnected module of the downregulated genes contained 30 proteins and 310 interactions corresponding to ECM organization, cell adhesion, cell-matrix interactions, tissue integrity, and remodeling (Fig. 4D). These findings highlight the extensive overlap in transcriptomic features spanning system-wide liver functions between the accepted and rejected-2 groups of donor livers. In particular, the molecular characterization and pathway enrichment analysis indicated that the high expression of liver metabolic functional genes and low expression of fibrotic genes are likely to be similar in the accepted and rejected-2 donor livers, with implications on the suitability for transplantation of the latter group.
Integrated Molecular profiling, Histopathological analysis, and Clinical factors predict a cohort of rejected donor livers as potentially suitable for transplantation
We identified a subset within the rejected group, rejected-2, that is molecularly similar to accepted donor livers. We filtered the rejected-2 subset based on histopathological analysis to assess micro- and macrosteatosis, fibrosis, necrosis, and hepatocyte ballooning (Table 2; Fig. 5). As expected, the accepted group showed no steatosis or minimal steatosis, no fibrosis, or perisinusoidal fibrosis (Fig. 5A, B). In contrast, donor livers within the rejected-2 subset displayed a range of tissue steatosis, varying from no steatosis to 40–70% steatosis. Fibrosis in the rejected-2 donor livers ranged from no fibrosis, perisinusoidal fibrosis, portal, and periportal fibrosis to bridging fibrosis (Fig. 5A, B). Conversely, the rejected-1 subset showed steatosis exceeding 33% in most biopsy samples, with some exhibiting more than 50% steatosis. Furthermore, most liver biopsies in the rejected-1 subset demonstrated fibrosis, ranging from perisinusoidal fibrosis and portal & periportal fibrosis to bridging fibrosis and even cirrhosis. Additionally, some donor livers in the rejected-1 group show high necrosis and hepatocyte ballooning, whereas there was no necrosis in the accepted donor livers (Table 2). There was no necrosis in most of the rejected-2 group except for two donor livers with mild to moderate necrosis. We compared AST, ALT, and total bilirubin levels across the three liver groups. Consistent with the transcriptomics findings, there were no statistical differences between the accepted and rejected-2 groups, whereas there were significant differences in the AST, ALT, and total bilirubin levels compared to the rejected-1 group (Fig. 5C). Histopathological filtering of the molecularly identified rejected-2 subset revealed 5 of the 12 donor livers as potentially transplantable (highlighted in Fig. 5B). These findings highlight the potential utility of incorporating transcriptomic analysis alongside current histopathological selection criteria to aid in decision-making regarding the transplantability of donor livers.
Discussion
We implemented an integrated transcriptomics and histopathological approach to characterize deceased donor livers that were either accepted or rejected for transplantation (Fig. 6). We performed bulk RNA sequencing to characterize the global gene expression profiles. Additionally, we conducted a histopathological evaluation through manual assessment by a gastrointestinal pathologist and AI-based image analysis. Gene expression analysis revealed high variability within rejected donor livers, with a fraction of rejected donor livers showing extensive transcriptomic overlap with the accepted group. Liver functional tests for serum AST, ALT, and total bilirubin showed similar overlapping patterns. These results suggest functional similarities across critical liver pathophysiological processes related to metabolic functions and ECM organization between some rejected donor livers and the accepted group, pointing to the potential suitability of transplantation of these livers. Importantly, the transcriptomic pattern of this subset of potentially suitable livers was enriched for a gene expression signature of graft success post-transplantation. Additional filtering of the rejected group by histopathological evaluation identified a select subset of donor livers that are likely suitable for transplantation.
Guidelines for assessing the suitability of donor livers for transplantation are lacking, leading to considerable variability among transplant centers [21]. Our approach captures the variability of selection criteria across transplant centers. This variability affects the rejected group as livers considered borderline by one center may be acceptable for transplantation by another. The wide range of gene expression at least partly reflects the variability in the selection criteria that led to rejection for transplantation. The accepted group was less variable in gene expression, likely reflecting the consistency of the selection criteria within a transplant center (TJUH). Together, the transcriptomic consistency of the accepted group and the relatively higher molecular heterogeneity of the rejected group enabled us to identify a subset of livers likely suitable for transplantation based on the molecular state of the donor liver at the time of organ recovery. Our results are further supported by the enrichment of a published transcriptomics signature linked to post-transplantation graft failure or success [27]. Follow-up studies can build on our findings using perfusion and functional testing of livers identified as potentially transplantable through an integrated transcriptomics and histopathological approach.
Biomarker assay development is largely focused on predicting post-liver transplant success [26, 27]. There are no studies on a diagnostic biomarker-based assay to support decision-making in selecting donor livers during organ recovery. Our current transcriptomics data provides a resource to identify such predictive biomarkers. The underlying transcriptomic profiles of the livers typically accepted for transplantation reflect preserved functions crucial for liver homeostasis and show no evidence of tissue damage. Additionally, linking histopathological features to underlying molecular characteristics corresponding to functions critical to liver pathophysiology may help predict the suitability of donor livers solely based on histopathological analysis. Such an approach has the potential to overcome the existing challenges of obtaining molecular profiling data within the critical time frame between organ recovery and transplantation and also serves as a cost-effective alternative to expensive multigene molecular assays. For example, in the case of renal transplantation, a molecular assay based on gene expression signature was developed to complement histopathological evaluation [38]. A combined use of molecular and histopathological evaluation has shown promise in overcoming the considerable variability across centers in evaluating renal transplant suitability of donor kidneys [38,39,40]. Recent efforts through the VITTAL study demonstrated that a large fraction (71%) of otherwise rejected livers can be transplanted after normothermic perfusion with 100% graft and patient survival for 90 days and significant patient survival over 5 years post-transplantation [41]. These results are consistent with our transcriptomic findings that a significant fraction of currently rejected livers retain a molecular state similar to that of accepted livers, demonstrating the potential for recovering liver functional capacity and leading to successful graft function post-transplantation.
There is an opportunity to build on our transcriptomics results to define molecular states at the time of organ recovery to predict post-transplant outcomes. Including molecular profiling and defining molecular states of the marginal livers can help guide the parameters for the downstream methods of organ transplant. While the marginal livers could expand the pool of transplantable livers, there is a higher chance of primary graft failure (PGF) compared with standard criteria donor (SCD) grafts [42]. Extrapolating from our results on the molecular overlap between accepted and a subset of rejected livers, the transplantation success of some of the marginal livers could potentially be attributed to the similarity of the underlying molecular functionalities of the marginal livers with that of livers typically regarded as acceptable for transplantation. Importantly, this similarity of underlying molecular state may not be readily apparent in the current clinical practice of histopathological assessment of the marginal livers. Biomarker evaluation guided by our transcriptomics results and potential integration with histopathology may improve decision-making in evaluating marginal livers for transplantation. Our study includes both DBD and DCD livers in all sample groups. The rejected-1 and rejected-2 groups included both DCD and DBD cases, demonstrating that our results in subgrouping the rejected donors are not biased by the underlying cause of death. We applied an integrated transcriptomics and histopathological approach and identified the subset within the rejected group. Our data presents an opportunity for leveraging transcriptomic information to discover molecular biomarkers that could play a crucial role in the selection of DCD livers for transplantation.
Conclusions
The present study utilized an integrated transcriptomic and histopathological approach to characterize transplant donor livers from biopsies of deceased donor livers either accepted or rejected for liver transplantation at the time of organ collection. Our integrated approach led us to identify potentially transplantable livers within the rejected group, which leads to exciting translational opportunities for our findings, pointing to an opportunity to expand the pool of donor livers that could be used in transplantation. Our findings hold the potential to address the variability observed between transplant centers in the selection of donor livers for transplantation, where the current evaluation primarily relies on histopathology, and there is variability in the application of the selection of donor organs for transplantation. In future endeavors, we aim to construct a comprehensive landscape of morphological features in donor liver samples and integrate these findings with our transcriptomics data to strengthen the ability to predict the suitability of the donor liver for transplantation. Moreover, we intend to leverage our integrated approach to predict post-transplant success. Our results demonstrate a proof of principle for the utility of molecular characterization and histopathological assessment in evaluating donor livers for transplantation.
Data availability
RNA sequencing data are available on NCBI Gene Expression Omnibus: GEO accession ID GSE243887.
Abbreviations
- ALB:
-
Albumin
- ALT:
-
Alanine aminotransferase
- APP:
-
Application Protocol Packages
- AST:
-
Aminotransferase enzyme
- CPM:
-
Counts per million
- CYP2E1:
-
Cytochrome P450 Family 2 Subfamily E Member 1
- DBD:
-
Donation after brain death
- DCD:
-
Donation after circulatory death
- ECD:
-
Extended criteria donor
- GOT:
-
Glutamic-oxaloacetic transaminase gene
- GPT:
-
Glutamate pyruvate transaminase gene
- H&E:
-
Hematoxylin & Eosin
- MASLD:
-
Metabolic dysfunction-associated steatotic liver disease
- MCODE:
-
Molecular Complex Detection
- MELD:
-
Model for End-Stage Liver Disease
- OPTN:
-
Organ Procurement and Transplantation Network
- PFA:
-
Paraformaldehyde
- PGF:
-
primary graft failure
- PPI:
-
Protein-protein interaction
- SCD:
-
Standard criteria donor
- SLC2A2:
-
Carrier Family 2 Member 2
- TJUH:
-
Thomas Jefferson University Hospital
- UMAP:
-
Uniform manifold approximation and projection
- WSI:
-
Whole slide images
References
Russo MW, Galanko JA, Shrestha R, Fried MW, Watkins P. Liver transplantation for acute liver failure from drug induced liver injury in the United States. Liver Transpl. 2004;10(8):1018–23.
O’Grady J. Timing and benefit of liver transplantation in acute liver failure. J Hepatol. 2014;60(3):663–70.
Donnelly MC, Hayes PC, Simpson KJ. The changing face of liver transplantation for acute liver failure: Assessment of current status and implications for future practice. Liver Transpl. 2016;22(4):527–35.
Prince MI, Hudson M. Liver transplantation for chronic liver disease: advances and controversies in an era of organ shortages. Postgrad Med J. 2002;78(917):135–41.
Van Thiel DH. Liver transplantation for alcoholics with Terminal Liver Disease. Alcohol Health Res World. 1996;20(4):261–5.
Mathurin P, Lucey MR. Liver transplantation in patients with alcohol-related liver disease: current status and future directions. Lancet Gastroenterol Hepatol. 2020;5(5):507–14.
Tsoulfas G, Goulis I, Giakoustidis D, Akriviadis E, Agorastou P, Imvrios G, et al. Hepatitis C and liver transplantation. Hippokratia. 2009;13(4):211–5.
Crespo G, Mariño Z, Navasa M, Forns X. Viral hepatitis in liver transplantation. Gastroenterology. 2012;142(6):1373–e13831.
Santopaolo F, Lenci I, Milana M, Manzia TM, Baiocchi L. Liver transplantation for hepatocellular carcinoma: where do we stand? World J Gastroenterol. 2019;25(21):2591–602.
Haydon GH, Neuberger J. Liver transplantation of patients in end-stage cirrhosis. Baillieres Best Pract Res Clin Gastroenterol. 2000;14(6):1049–73.
OPTN Metrics. [cited 2023 Nov 2]. https://insights.unos.org/OPTN-metrics/.
Radunz S, Paul A, Nowak K, Treckmann JW, Saner FH, Mathé Z. Liver transplantation using donor organs with markedly elevated liver enzymes: how far can we go? Liver Int. 2011;31(7):1021–7.
Kaltenbach MG, Harhay MO, Abt PL, Goldberg DS. Trends in deceased donor liver enzymes prior to transplant: the impact on graft selection and outcomes. Am J Transpl. 2020;20(1):213–9.
Lee JH, Mintz PD. Graft versus host anti-rho(D) following minor Rh-incompatible orthotopic liver transplantation. Am J Hematol. 1993;44(3):168–71.
Kluger MD, Guarrera JV, Olsen SK, Brown RS, Emond JC, Cherqui D. Safety of blood group A2-to-O liver transplantation: an analysis of the United Network of Organ Sharing database. Transplantation. 2012;94(5):526–31.
Garzoni C, Ison MG. Uniform definitions for donor-derived infectious disease transmissions in solid organ transplantation. Transplantation. 2011;92(12):1297–300.
Nam H, Nilles KM, Levitsky J, Ison MG. Donor-derived viral infections in liver transplantation. Transplantation. 2018;102(11):1824–36.
Markin RS, Wisecarver JL, Radio SJ, Stratta RJ, Langnas AN, Hirst K, et al. Frozen section evaluation of donor livers before transplantation. Transplantation. 1993;56(6):1403–9.
Abdelsameea E, Alsebaey A, Ehsan N, Elrefaey A, Aboelela K, Ibrahim E-S, et al. Liver biopsy is still needed in liver transplantation recipients: a single center experience. Egypt Liver J. 2020;10(1):18.
Reddy MS, Varghese J, Venkataraman J, Rela M. Matching donor to recipient in liver transplantation: relevance in clinical practice. World J Hepatol. 2013;5(11):603–11.
Schaffhausen CR, Bruin MJ, Chu S, Fu H, McKinney WT, Schladt D, et al. Tool to aid patients in selecting a liver transplant center. Liver Transpl. 2020;26(3):337–48.
Wertheim JA, Petrowsky H, Saab S, Kupiec-Weglinski JW, Busuttil RW. Major challenges limiting liver transplantation in the United States. Am J Transpl. 2011;11(9):1773–84.
Haque O, Yuan Q, Uygun K, Markmann JF. Evolving utilization of donation after circulatory death livers in liver transplantation: the day of DCD has come. Clin Transpl. 2021;35(3):e14211.
de Vera ME, Lopez-Solis R, Dvorchik I, Campos S, Morris W, Demetris AJ, et al. Liver transplantation using donation after cardiac death donors: long-term follow-up from a single center. Am J Transpl. 2009;9(4):773–81.
Croome KP, Taner CB. The changing landscapes in DCD liver transplantation. Curr Transpl Rep. 2020;7(3):194–204.
Hoyer DP, Swoboda S, Treckmann JW, Benkö T, Paul A, Brocke-Ahmadinejad N. Transcriptomic profiles of human livers undergoing rewarming machine perfusion before transplantation-first insights. Funct Integr Genomics. 2021;21(3–4):367–76.
Defamie V, Cursio R, Le Brigand K, Moreilhon C, Saint-Paul MC, Laurens M, et al. Gene expression profiling of human liver transplants identifies an early transcriptional signature associated with initial poor graft function. Am J Transpl. 2008;8(6):1221–36.
Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J, Wilm A, et al. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol. 2020;38(3):276–8.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.
Frankish A, Diekhans M, Ferreira A-M, Johnson R, Jungreis I, Loveland J, et al. GENCODE reference annotation for the human and mouse genomes. Nucleic Acids Res. 2019;47(D1):D766–73.
Love MI, Huber W, Anders S. Moderated estimation of Fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Yang Y, Sun H, Zhang Y, Zhang T, Gong J, Wei Y, et al. Dimensionality reduction by UMAP reinforces sample heterogeneity analysis in bulk transcriptomic data. Cell Rep. 2021;36(4):109442.
Hart SN, Therneau TM, Zhang Y, Poland GA, Kocher J-P. Calculating sample size estimates for RNA sequencing data. J Comput Biol. 2013;20(12):970–8.
Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022;50(W1):W216–21.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape stringapp: network analysis and visualization of proteomics data. J Proteome Res. 2019;18(2):623–32.
Bader GD, Hogue CWV. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Halloran PF, Famulski KS, Reeve J. Molecular assessment of disease states in kidney transplant biopsy samples. Nat Rev Nephrol. 2016;12(9):534–48.
Schutter R, Sanders J-SF, Ramspek CL, Crop MJ, Bemelman FJ, Christiaans MHL, et al. Considerable variability among transplant nephrologists in judging deceased donor kidney offers. Kidney Int Rep. 2023;8(10):2008–16.
Madill-Thomsen K, Perkowska-Ptasińska A, Böhmig GA, Eskandary F, Einecke G, Gupta G, et al. Discrepancy analysis comparing molecular and histology diagnoses in kidney transplant biopsies. Am J Transpl. 2020;20(5):1341–50.
Mergental H, Laing RW, Kirkham AJ, Clarke G, Boteon YL, Barton D, et al. Discarded livers tested by normothermic machine perfusion in the VITTAL trial: secondary end points and 5-year outcomes. Liver Transpl. 2024;30(1):30–45.
Huang C, Chen Z, Wang T, He X, Chen M, Ju W. A marginal liver graft with hyperbilirubinemia transplanted successfully by ischemia-free liver transplantation. Ann Transl Med. 2021;9(5):425.
Funding
Research was supported by the National Institute on Alcohol Abuse and Alcoholism R01 AA018873 and T32 AA007463 to RV and JH, F31 AA030214 to AM, and the Gift of Life Donor Program, the Transplant Foundation.
Author information
Authors and Affiliations
Contributions
J.B.H. and R.V. developed the concept and designed the research; A.P.S., J.B.H. and R.V. designed the study protocol. A.S. performed the experiments; A.S. and A.M. analyzed the data; A.S., A.M., A.P.S., J.B.H., and R.V. interpreted the results of the experiments; A.S., A.M., and R.V. prepared the figures; J.W. and M.A. performed histopathological scoring; B.K.B collected the samples and collated clinical data; A.S. drafted the initial manuscript; J.B.H., A.P.S., and R.V. edited and revised the manuscript; all authors approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study used deidentified liver tissue samples and data obtained from deceased individuals. As per the regulations outlined in 45 CFR 46.101 (b) [4], the current study falls under Exemption 4, and does not require additional Institutional review board approval. The need for approval was waived by the Thomas Jefferson University Institutional review board. A Certification for Protected Health Information of Decedents for Research was obtained from the Privacy Officer of the Office of University Counsel and Corporate Compliance Department of Thomas Jefferson University Hospital (TJUH) for the use of deidentified data in this study. The informed consent was obtained by the Gift of Life Foundation from the deceased donor`s next of kin to collect liver biopsy samples for research purposes.
Consent for publication
Not Applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Srivastava, A., Manchel, A., Waters, J. et al. Integrated transcriptomics and histopathology approach identifies a subset of rejected donor livers with potential suitability for transplantation. BMC Genomics 25, 437 (2024). https://doi.org/10.1186/s12864-024-10362-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10362-7