Background: Pulmonary hypertension is a condition in which blood pressure in the lung's arterial circulation is abnormally elevated, defined by a mean pulmonary artery pressure of 20 mmHg or greater measured by right heart catheterization. In patients who undergo lung transplantation, the removed lung (the explant) shows characteristic changes in the blood vessel walls that reflect years of elevated pressure -- including thickening of vessel walls, abnormal cell proliferation, and in severe cases, complex vascular malformations. These changes are visible on standard hematoxylin and eosin (H&E) stained microscopy slides, but evaluating them systematically across an entire slide is time-consuming and subjective. No artificial intelligence system for detecting pulmonary hypertension from explant pathology slides has been previously described.
What we did: We applied a machine learning method called attention-based multiple instance learning to whole-slide digital images of lung transplant explants from 18 patients (10 with pulmonary hypertension, 8 without). The AI system was trained to classify each case as pulmonary hypertension or normal using only the slide images and the clinical pressure measurement as ground truth -- no hand-drawn annotations or region labels were required. We used a leave-one-out cross-validation strategy, meaning for each patient the model was trained on the other 17 patients and then tested on that patient's slides, giving an unbiased estimate of how the model would perform on a new case.
Results: The model correctly classified 16 of 18 cases at a probability cutoff of 0.5, and 17 of 18 at the statistically optimal cutoff. The area under the receiver operating characteristic curve (a standard measure of discriminative ability ranging from 0.5 for chance to 1.0 for perfect) was 0.975. Sensitivity (the proportion of true pulmonary hypertension cases identified) was 90%, and specificity (the proportion of true non-PH cases correctly excluded) was 100%. There were no false positives -- every case the model called pulmonary hypertension actually was.
Conclusion: A foundation model artificial intelligence system can detect pulmonary hypertension from lung transplant explant slides with strong accuracy, using only the whole slide and the clinical diagnosis -- no specialized region-by-region annotation. These results warrant validation in a larger, multi-center cohort.
Pulmonary hypertension (mean pulmonary artery pressure, or mPAP, of 20 mmHg or greater at rest, by right heart catheterization) substantially worsens the prognosis of patients with chronic lung disease, including interstitial lung disease (ILD), chronic obstructive pulmonary disease (COPD), and connective tissue disorders. At the time of lung transplantation, the explanted lung provides a uniquely comprehensive tissue sample capturing the full burden of disease across multiple lobes and decades of pathological change.
The blood vessels in these explants show a spectrum of changes that correlate with pulmonary artery pressure: medial (muscular wall) hypertrophy, intimal (inner layer) proliferation and fibrosis, and in severe cases, plexiform lesions -- complex, disorganized vascular growths that represent end-stage vascular remodeling. Prior work by Seki and colleagues (2019), from the Wallace and Saggar research groups at USC and Cedars-Sinai, demonstrated that the presence of capillary proliferation specifically in preserved alveolated lung -- not in fibrotic or destroyed regions -- was the strongest histologic predictor of pulmonary hypertension, with a c-statistic of 0.869. This is an important distinction: the predictive signal is not where one might intuitively look (the obviously diseased regions), but in the relatively preserved parenchyma. Any AI model trained to detect pulmonary hypertension should therefore not restrict itself to grossly abnormal tissue zones.
Separately, Dotan and colleagues (2020) showed that severe histologic vascular changes can be present even when the measured pulmonary artery pressure is not dramatically elevated, suggesting that tissue-based assessment may identify disease severity that hemodynamic measurements alone miss.
Foundation models for digital pathology are large neural networks pretrained on millions of pathology image patches across many tissue types and diseases. They produce rich numerical representations of image patches that capture a wide range of visual features without being trained specifically for any one task. We used phikon-v2, a publicly available pathology foundation model, to extract these representations from every tissue patch across each whole slide, then trained a lightweight classifier on top of those representations to predict pulmonary hypertension at the case level.
Eighteen consecutive lung transplant explant cases from USC Keck Medical Center were included in this analysis, provided by Dr. William D. Wallace. Ten patients had pulmonary hypertension (mean pulmonary artery pressure 26 to 51 mmHg, measured by right heart catheterization at the time of transplant listing). Eight patients had no pulmonary hypertension and served as controls. Slides were digitized at 20x magnification using an Aperio scanner, producing 128 whole-slide images in total. The number of slides per patient ranged from 3 to 10 (median 7). Some slides contain expert annotations of anatomical compartments; these annotations were not used in this analysis and are reserved for a planned follow-up analysis.
All 128 slides were screened for pen/grossing ink, tissue fold, focus artifact, and dark spot using an in-house gated-attention multiple instance learning quality control classifier (Track B; pen AUROC 0.813 on an independent BRAF cohort). No pen or grossing ink artifact was identified in the WDW cohort. All slides were retained for analysis.
Each whole-slide image was divided into non-overlapping 256x256-pixel square patches at 20x magnification. Patches containing predominantly background (glass or white space) were excluded using tissue masking. The remaining tissue patches were each processed through phikon-v2, a pathology foundation model based on a Vision Transformer Large (ViT-L) architecture that produces a 1,024-dimensional numerical feature vector for each patch. These vectors were stored for each slide along with the spatial coordinates of each patch within the slide. Feature extraction was performed on an NVIDIA RTX 5090 graphics processing unit and required one overnight run to process all 128 slides.
We used attention-based multiple instance learning (ABMIL), a framework designed for whole-slide image analysis in which the final prediction is made at the case (patient) level rather than at the patch level. The model learns to assign an importance weight (attention score) to each image patch, then combines all patches from a patient's slides into a single case-level prediction using those weights as a weighted average. Patches the model finds informative receive high attention scores; uninformative patches are down-weighted. For each patient, all patches from all available slides were combined into one bag of feature vectors fed to the model together.
The classifier itself consists of: (1) a patch encoder that projects the 1,024-dimensional features to a 512-dimensional space using a linear layer, layer normalization, and a nonlinear activation; (2) a gated attention mechanism using parallel tanh and sigmoid pathways that produces a scalar attention weight for each patch; (3) a weighted aggregation step that produces one 512-dimensional vector per patient; and (4) a two-layer classification head that produces the final probability of pulmonary hypertension. The model was trained using the AdamW optimizer, a cosine learning rate schedule, and binary cross-entropy loss. Early stopping was applied when training loss did not improve for 20 consecutive training epochs, with a maximum of 60 epochs. During training, patches were randomly subsampled to a maximum of 8,192 per patient to manage memory; at test time, all patches were used.
With only 18 patients, standard multi-fold cross-validation would leave too few patients in each test set to provide meaningful evaluation. We therefore used leave-one-out cross-validation: for each patient, a fresh model was trained from scratch on all 17 other patients, and that model was then applied to the held-out patient whose slides it had never seen. This produced one predicted probability of pulmonary hypertension for each of the 18 patients. These 18 probabilities were then used to compute the area under the receiver operating characteristic curve (AUROC), sensitivity, specificity, and classification accuracy. The Youden index was used to identify the probability threshold that maximized the sum of sensitivity and specificity. Total runtime: 77 minutes on an RTX 5090 graphics processing unit.
The model produced highly confident predictions for most cases. Fifteen of 18 patients received predicted probabilities below 0.005 or above 0.999, indicating strong separation between the two groups. The two cases the model missed have clear explanations (see Discussion).
| Case | True Diagnosis | Pulmonary Artery Pressure | No. of Slides | AI Probability of PH | AI Prediction | Outcome |
|---|---|---|---|---|---|---|
| 1 | PH | 35 mmHg | 5 | 0.999 | PH | Correct |
| 3 | PH | 51 mmHg | 6 | 0.999 | PH | Correct |
| 60 | Normal | Normal (no PH) | 5 | 0.001 | Normal | Correct |
| 71 | PH | 26 mmHg Borderline mPAP; mildest PH in cohort | 4 | 0.001 | Normal | Incorrect |
| 103 | PH | 32 mmHg | 10 | 1.000 | PH | Correct |
| 118 | Normal | Normal (no PH) | 10 | 0.001 | Normal | Correct |
| 126 | PH | 43 mmHg | 10 | 1.000 | PH | Correct |
| 128 | PH | 49 mmHg | 10 | 1.000 | PH | Correct |
| 161 | PH | 29 mmHg | 6 | 1.000 | PH | Correct |
| 166 | PH | 46 mmHg | 10 | 1.000 | PH | Correct |
| 185 | PH | 49 mmHg | 10 | 1.000 | PH | Correct |
| 213 | Normal | Normal (no PH) | 6 | 0.004 | Normal | Correct |
| 265 | Normal | Normal (no PH) | 10 | 0.001 | Normal | Correct |
| 319 | Normal | Normal (no PH) | 8 | 0.000 | Normal | Correct |
| 331 | Normal | Normal (no PH) | 5 | 0.000 | Normal | Correct |
| 344 | Normal | Normal (no PH) | 5 | 0.001 | Normal | Correct |
| 353 | PH | 40 mmHg Only 3 slides available | 3 | 0.051 | Normal | Incorrect |
| 463 | Normal | Normal (no PH) | 5 | 0.000 | Normal | Correct |
This study provides the first demonstration that a foundation model artificial intelligence system can detect pulmonary hypertension from lung transplant explant whole-slide images with clinically meaningful accuracy, using only slide images and the clinical hemodynamic diagnosis -- no region-by-region annotation by a pathologist was required.
The striking feature of the results is not just the overall accuracy but the extreme confidence of most individual predictions. Probabilities were clustered near 0 (control-like) or near 1 (PH-like), with very little in between. This bimodal distribution suggests the model has learned a global texture or architectural signature that differs qualitatively between pulmonary hypertension and non-PH explants, rather than detecting rare focal lesions. This is consistent with the histological understanding of PH, in which vascular remodeling is a diffuse, multifocal process involving many vessels throughout the lung.
The perfect specificity (zero false positives) is particularly noteworthy from a clinical standpoint. A screening tool that produces no false alarms -- where every case it calls pulmonary hypertension actually has it -- is highly practical. This likely reflects the distinct and consistent global appearance of PH tissue versus non-PH tissue even in the setting of other underlying lung disease.
This patient has the mildest pulmonary hypertension in the cohort. A mean pulmonary artery pressure of 26 mmHg is above the diagnostic threshold of 20 mmHg, but only just. At this level of disease, the histologic vascular changes are likely early and subtle -- modest muscular hypertrophy of small vessels, minimal intimal thickening, no plexiform lesions. The model produced a probability of 0.0008, essentially identical to its scores for normal control cases. This may represent a genuine limitation at the mild end of the disease spectrum rather than a failure of the system. Additionally, this patient had only 4 slides available -- the fewest of any pulmonary hypertension patient.
This patient has moderate pulmonary hypertension, well above threshold. The model's probability was 0.051 -- much higher than for the other missed case and just below the decision boundary, suggesting the model was detecting something. However, this patient had only 3 slides -- the fewest of any case in the entire cohort. Three slides represent a fraction of the tissue available from a typical lung explant. The most likely explanation is insufficient tissue sampling: the vascular changes were not adequately represented in the three available digitized sections. Retrieving additional archived slides from pathology blocks for this patient would be the most informative next step.
All eight control patients (those without pulmonary hypertension) were correctly classified, with probabilities ranging from 0.000 to 0.004. Importantly, some of these patients had significant underlying lung disease (interstitial lung disease, fibrosis) that produces secondary vascular changes. The model nonetheless distinguished the absence of primary pulmonary hypertension remodeling from these secondary changes, suggesting it has learned a specific pattern associated with hemodynamically confirmed pulmonary hypertension rather than generic lung injury.
Seki and colleagues demonstrated that the pathologist-assessed presence of capillary proliferation in preserved (non-fibrotic) alveolated lung had a c-statistic of 0.869 for predicting pulmonary hypertension -- a strong result from careful manual assessment. The Computational Pathology Kit (CPK), requiring no manual assessment at all, achieved an AUROC of 0.975 in a leave-one-out analysis. The attention maps shown in the Attention Maps tabs of this report provide visual evidence about which tissue regions drove the AI predictions, and may confirm or extend the finding that preserved parenchyma (rather than fibrotic or destroyed zones) is where the discriminative signal lives.
The cohort size of 18 cases is small. Leave-one-out cross-validation provides an unbiased estimate but the resulting confidence intervals are wide. All cases are from a single institution, a single scanner, and a single pathology processing protocol. The pulmonary hypertension etiology varies across cases (Group 3 pulmonary hypertension associated with lung disease being the most common in this transplant population) and this variation has not been analyzed separately. The model has not been tested on cases from any other institution. Prospective external validation is essential before any clinical conclusions can be drawn.
No publicly available labeled pulmonary hypertension lung transplant explant whole-slide image dataset exists (as of May 2026). The most realistic path to external validation is via expansion of the USC/Keck cohort or collaboration with other high-volume lung transplant programs. A multi-institutional dataset of 50–90 cases would provide the sample size needed for publication-grade validation.
Attention-based multiple instance learning applied to phikon-v2 foundation model features achieves an AUROC of 0.975, sensitivity of 90%, specificity of 100%, and zero false positives for detecting pulmonary hypertension from lung transplant explant whole-slide images in leave-one-out cross-validation across 18 cases. The two misclassifications have straightforward clinico-pathologic explanations: the mildest disease in the cohort, and the case with the fewest available slides. The system requires no pathologist annotation and can process a full case in minutes once features are extracted. External validation in a multi-center cohort is the clear next step.
Each map shows the distribution of the model's attention weights across the tissue patches in the highest-attention slide for that patient, from the leave-one-out analysis (meaning the model had never seen this patient's slides during training). Brighter regions received higher attention from the model when making its prediction. These are not pathology annotations -- they show where the model focused, which can then be compared to the tissue features a pathologist sees at those locations. See the next tab for side-by-side H&E comparisons on representative cases.
Red bars = pulmonary hypertension cases, teal bars = controls. Dashed line = 0.5 decision boundary. The two cases labeled MISS are the ones the model incorrectly classified.
Six representative patients are shown: two confidently correct pulmonary hypertension cases, two confidently correct controls, and the two cases the model missed. For each patient, the slide with the highest mean attention score was selected. Two views are provided: a full-slide overview at approximately 2.5x magnification, and a zoomed crop (~10x) centered on the region of the slide where the model focused most. The left half of each image is the unmodified hematoxylin and eosin stain; the right half is the same image with the model's attention scores overlaid in color (inferno colormap: black = low attention, yellow/white = high attention).
This presentation is intended to allow pathologic interpretation of what the model is responding to. The critical question is whether the model's attention concentrates in regions with known pulmonary hypertension histologic features -- particularly in preserved alveolated tissue rather than fibrotic or destroyed zones -- consistent with the Seki et al. (2019) findings.
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay
Full slide view (~2.5x) — Left: H&E | Right: same slide with attention overlay (brighter = higher model focus)
Zoomed region (~10x) showing area of highest model attention — Left: H&E | Right: attention overlay