Abstract
Survival after lung transplantation lags that of other solid organ transplants. Long-term survival is hampered primarily due to chronic lung allograft dysfunction (CLAD) development. It remains elusive how (chronic) rejection is organized within the lung graft over time post-transplant. Using a model of orthotopic left lung transplantation in major mismatched mouse strains with daily immunosuppression, we aimed to study the spatiotemporal dynamics of (chronic) rejection, using micro-computed tomography imaging, flow cytometric analyses and spatial proteomics. Endothelial cells demonstrated early activation and destruction (day 7 post-transplant). The accompanying early inflammation at the vascular compartment, progressed towards aberrant tissue repair resulting in irreversible bronchovascular fibrosis and chronic graft dysfunction. We provide new insights in the spatiotemporal dynamics of (chronic) rejection with a vascular-oriented onset that may have future implications for diagnosis and treatment in clinical lung transplantation.
Graphical Abstract
Introduction
Lung transplantation is the only curative treatment for end-stage lung diseases. Despite advances in donor selection, surgical techniques, and immunosuppressive strategies, survival after lung transplantation lags that of other solid organ transplants with a median survival of only 6 years post-transplantation [–]. Beyond the first post-transplant year, survival is predominantly limited by the development of chronic lung allograft dysfunction (CLAD) [, , ]. CLAD is clinically defined by an irreversible decline in lung function and represents an umbrella term encompassing heterogeneous pathological processes. It is not exclusively defined by rejection, but rather as a broad term for irreversible loss of lung function. Three clinical phenotypes of CLAD are currently recognized, including, ‘bronchiolitis obliterans syndrome’ (BOS), ‘restrictive allograft syndrome’ (RAS) and a mixed phenotype []. These phenotypes reflect the complexity of CLAD from a clinical point of view. Importantly, this CLAD definition does not provide insights into the underlying drivers of graft injury [, ]. CLAD has historically been conceptualized as a predominantly airway-centered disease, with clinical and research efforts largely focused on bronchiolar pathology [–]. However, this airway-centric view contrasts with established paradigms of vascular chronic rejection that is well characterized in other solid organ transplants, where chronic rejection has been perceived as a vessel-centered process in which persistent immunological and endothelial activation induces thrombosis, intimal proliferation, and vascular remodeling act as hallmarks of chronic allograft vasculopathy [–].
The apparent disconnect between lung transplantation and other solid organs is therefore striking, particularly given the shared immunological principles governing allograft rejection. Considering transplant rejection can be described as an immune-mediated recognition of the donor graft HLA antigens as ‘non-self’, where both cellular and humoral mechanisms of the innate and adaptive immune system get involved [].
So, despite vascular chronic rejection being well characterized in other transplanted organs and some studies within the field of lung transplantation reporting perivascular inflammation histologically, its central role and the temporal and spatial relationships between vascular injury and graft dysfunction have never been investigated in the context of lung transplantation, and direct evidence in human lung allografts remains limited [, , , ]. This conceptual gap has contributed to an incomplete understanding of chronic lung allograft injury and may partially explain the lack of effective targeted therapies [, , ].
Importantly, the traditional dichotomous classification of rejection into “acute” and “chronic” forms complicates interpretation. In clinical practice, rejection cannot be regarded as a binary or purely temporal process. Lung transplant recipients may experience chronically ongoing antibody-mediated rejection, recurrent or persistent episodes of acute rejection, or overlapping immunopathological processes resulting in chronic graft injury []. Distinguishing “chronic” rejection as a pathophysiological mechanism from chronically ongoing forms of “acute” rejection is therefore essential [, –]. Controlled experimental models may act as a valuable tool to obtain essential and fundamental information on the spatiotemporal dynamics of vascular involvement in (chronic) rejection [–].
Building on these observations, we hypothesize that (chronic) rejection in lung transplantation may originate in the vascular compartment with a significant role for endothelial cells. We aim to map the temporal and spatial progression of vascular pathology in a murine lung transplant model to understand how vascular injury contributes to chronic graft dysfunction. We will use our murine model of rejection after orthotopic left lung transplantation and use advanced methodologies including multicolor flow cytometry, in vivo and ex vivo µCT imaging and spatial proteomics.
Materials and methods
Orthotopic mouse left lung transplantation and study design
Mice (male) were purchased from Janvier Labs (Le Genest-Saint-Isle, France) and transplanted between eight to 10 weeks of age. Recipients were C57BL/6N (H-2Kb) and donors were either C57BL/6N (H-2Kb) for isograft transplantation or BalbC (H-2Kd) for allograft transplantation. Serial sacrifice was performed at three different time points, on post-operative day 1, 7 and 35. All mice were housed in a conventional facility with individually ventilated cages (IVC) and received ad libitum standard chow and water. Orthotopic left lung transplantation was performed as previously described [22, , ]. After surgery, mice were allowed to recover on a heating pad overnight and received buprenorphine for the following 72 h. Both isografts and allografts received daily subcutaneous immunosuppression consisting of 10 mg/kg/day cyclosporine (Sandimmun®, Novartis, Vilvoorde, Belgium) and 1.6 mg/kg/day methylprednisolone (SoluMedrol®, Pfizer, Brussels, Belgium) started immediately after transplantation. The specifics of the study design are illustrated in Figure 1A. The experimental procedure was approved by the Ethical Committee for Animal Research at KU Leuven (P194/2019). To determine leukocyte and endothelial cell numbers and activation, flow cytometry was performed on the entire left lung graft at the three timepoints (n = 5–6/group/timepoint). Longitudinal, non-terminal in vivo µCT imaging was performed on day 1, 7 and 35 in the grafts (n = 6/group) of the flow cytometry experiment sacrificed day 35. Terminal exsanguination with retro-orbital blood collection was performed, prior to sacrifice for measuring the cyclosporine A trough level by an immunoassay (Dimension® RXL, Siemens Medical solutions, Diamond diagnostics, USA) (Supplementary Figure S1). The experiment with isografts and allografts was repeated for histological analyses (n = 4–5/group/timepoint). At the time of sacrifice, lungs were perfused with saline and inflated, followed by perfusion with 4% PFA at 4 °C for 24 h. The grafts were subsequently embedded in paraffine and processed into 4 µm thick sections which were stained with Hematoxylin-Eosin and Masson’s trichrome. Images were taken with an Olympus BX61 microscope. The percentage of collagen in a stained lung section with Masson’s trichrome was quantified using QuPath 0.3.2 with a personalized Pixel Classifier created with the ‘Train Pixel Classifier’ tool.
FIGURE 1
Flow cytometric analysis to quantify leukocyte and endothelial cell numbers and activation
At time of autopsy, the transplanted lung was removed and cells were isolated as previously described []. The lung was extracted and collected in RPMI buffer [RPMI GlutaMAX/FCS (5%)/ 1% penicillin/streptomycin/ 0.1% beta-mercaptoethanol], minced and incubated for 30 min at 37 °C in digestion medium consisting of 2 mg/mL collagenase D and 0.1 mg/mL DNase I in RPMI buffer. Subsequently, lung tissue was homogenized using a 20-gauge needle and new digestion medium was added, followed by a second incubation period of 15 min at 37 °C. Cells were washed and red blood cell lysed using 0.83% NH4Cl/10 mM Tris at 37 °C and passed through a 70 µm nylon cell strainer. After a last washing step, live leukocytes were counted in a Bürker chamber with trypan-blue staining.
For the flow cytometry 1.5-3 million lung cells were used. Cells were incubated with a viability dye, Zombie Aqua (Biolegend, San Diego, CA, USA) or Zombie UV (Biolegend), together with Mouse Fc block (MACS Miltenyi Biotec.) for 15 min at room temperature. After washing twice, cells were stained with a panel of monoclonal antibodies (Supplementary Table S1) dissolved in Brilliant stain buffer (BD Biosciences, Erembodegem, Belgium) for 20 min at 4 °C. After surface staining, 100,000 or 200,000 live single cells were analyzed per sample with a BD LSR Fortessa Flow cytometer (BD Biosciences). The flow cytometric panels used allow to quantify viable cells, leukocytes, endothelial cells numbers and activation status (Supplementary Table S1). For the activation status we used MHC1/2 (central for T cell activation), adhesion molecules ICAM-1 (present on leukocytes and endothelial cells) and VCAM-1 (present on endothelial cells). Data were analyzed with FlowJo v10 software (FlowJo LLC, Ashland, OR, USA) and cells were gated according to predefined gating strategies (Supplementary Figure S2). The absolute number of leukocytes and endothelial cells were calculated by the percentage of CD45+ and CD45− CD31+ cells, respectively, among the live cells multiplied by the total number of live cells determined with the Bürker chamber.
In vivo and ex vivo µCT imaging to identify and locate rejection
For in vivo µCT imaging free-breathing mice were scanned with a whole-body small animal µCT scanner (SkyScan 1278, Bruker µCT, Kontich, Belgium). Animals were anesthetized with isoflurane (2% in pure oxygen) and placed in supine position. Following scan parameters were used: 50 kVp X-ray source, 350 µA current, 1 mm aluminum X-ray filter and 150 ms exposure time per projection acquiring projections with 0.9° increments over a total angle of 220°. Acquisition resulted in a respiratory-weighted reconstructed 3D dataset with an isotropic voxel size of 50 µm in a total scan time of 3 min, associated with a radiation dose exposure of 69–89 mGy []. The images were reconstructed with the following parameters: smoothing of one, beam-hardening correction of 10%, post-alignment and ring artefact reduction were optimally set for each individual scan. Images were processed and calibrated to Hounsfield units (HU) as described before []. Image reconstruction, analysis and quantification were performed using NRecon, DataViewer and CTAn software provided by the manufacturer. Quantification of total lung volume and mean lung density was performed for a manually delineated region of interest (ROI), resulting in a volume of interest (VOI) on the transversal µCT images at end-expiration.
For ex vivo µCT imaging separate transplant allograft, isograft and non-transplant control mice were used (Figure 1A). The lungs were fixated as previously described and chemically dehydrated in a graded series of ethanol concentrations followed by hexamethyldisilazane (Sigma-Aldrich, Overijse, Belgium) submersion overnight. Dried lungs were scanned using a Skyscan 1272 µCT scanner (Bruker, Kontich, Belgium) with a resolution of 3.5 µm. 3D images were reconstructed using the NRecon software (version 1.7.0.4, Bruker microCT, Kontich, Belgium) and 3D segmentations of airway lumina, arteries and veins was performed semi-automatically by ITK-SNAP software [].
Bulk spatial GeoMx proteomics to validate vessel orientated rejection in lung sections
GeoMX spatial proteomic profiling was performed to determine the spatial localization of different immune-related proteins of interest. By selecting specific regions of interest based on the anatomical location, we performed a regional analysis subdividing airways, arteries, veins, and parenchyma (Figures 5A,B). Formalin-fixed and paraffin-embedded (FFPE) lung tissue sections (4 µm) were obtained from one isograft at each time point (day 1, 7 and 35) and three allografts at each timepoint. For each tissue section, three regions of interest (ROIs) were selected and paired for the parenchyma, bronchioles, arteries, and veins (total of 12 ROIs), based on the anatomical orientation in the H&E image. For the spatial proteomics, the GeoMx assay and platform was performed at the F Polverino lab in Baylor (Houston, USA). Tissue sections were incubated with 47 oligo-labelled primary antibodies of interest (Supplementary Table S2) and analyzed using the Nanostring GeoMx Digital Spatial Profiler as per manufacturer’s instructions. The detection antibodies comprised three fixed panels of the GeoMx assay (GeoMx immune cell profiling panel, GeoMx immune activation status panel, GeoMx immune cell typing panel).
Statistical analysis
Data analysis was performed using GraphPad statistical software (Prism, version 10, San Diego, CA, USA). To compare the different groups one-way ANOVA was used. To compare the different time points in isografts or allografts, a mixed effects model with Šidák multiple comparisons test was conducted. To compare isografts and allografts, a mixed effect model with Tukey’s multiple comparisons post hoc was used. A p-value of <0.05 was considered significant. Differential protein expression analyses (spatial proteomics) were performed to compare allografts and isografts in different compartments of the lung. Third quartile normalization was performed on the obtained expression matrix to allow comparison between samples. Of the 47 included proteins, 35 were assigned to a specific cluster (structural, hematopoietic, innate immunity, adaptive immunity, stromal and proliferation). To be able to pool multiple proteins per cluster, expression values were scaled to Z-scores for each protein separately. Statistical analysis on these differences was performed using mixed modeling in R (version 4.2.2), using the “lmer” function. The model used a full factorial design between graft (isograft and allograft), day (1, 7, 35) and tissue type (airway, artery, vein, parenchyma), with notation “graft*day*tissue”. Different mouse individuals were included as random effect. The least squares method (“emmeans” function) was used to perform pairwise comparisons between allograft and isograft.
Results
Lung morphology and density after murine orthotopic single lung transplantation
Isografts maintained normal macroscopic appearance throughout follow-up, whereas allografts progressively developed yellow discoloration with patchy dark areas by day 35 (Figure 1B). Longitudinal in vivo µCT revealed normal morphology in both groups on day 1; by day 7, allografts showed complete consolidation corresponding to atelectasis, while isografts remained aerated. On day 35, patchy dense regions persisted in allografts, contrasting with preserved isograft structure. Quantitative analysis confirmed these changes: mean lung density was significantly higher in allografts compared to isografts on day 35 (−87 ± 154 HU vs. −376 ± 27 HU, p = 0.016), while isografts showed a modest decline in density over time (p = 0.043) (Figure 1C). Lung volumes increased from day 7 to 35 in both groups (isograft p = 0.04; allograft p = 0.008) with no overall differences between both groups (Figure 1D). Histology confirmed these findings by preserved parenchymal architecture without cellular infiltrates or connective tissue deposition and fibrosis formation in isografts versus bronchovascular inflammation and fibrosis in allografts on day 35 (Figures 1E–G; Supplementary Figure S3).
Flow cytometry for leukocyte and endothelial cell numbers and activation status
The number of leukocytes (CD45+) significantly increased in the allograft left lung on day 7 (p = 0.037) but decreased again on day 35 (p = 0.020) (Figure 1H). These leukocytes significantly presented more ICAM-1 on their surface on day 7 (p = 0.0026) and 35 (p = 0.0045) compared to isografts, while no change in VCAM-1 was observed (Figures 1I,J).
Endothelial cell (CD31+) numbers remained stable in isografts. In contrast, allografts showed a significant reduction on day 7 compared to day 1 (p = 0.024), which only partially recovered by day 35 (p = 0.199) (Figure 2A). Endothelial cells showed to be activated, as VCAM-1 was significantly increased on day 7 compared to isografts (p = 0.047) and remained increased on day 35 versus day 1 (p = 0.049) (Figure 2B). No change in ICAM-1 expression by endothelial cells was observed (Figure 2C). Endothelial MHCI and MHCII expression were increased in allografts compared to isografts on day 7 (p = 0.045; p = 0.0002, respectively) and MHCII remained significantly increased on day 35 (p = 0.006) (Figures 2D,E).
FIGURE 2
Histological visualization of vessel-centered rejection
Isografts showed preserved vascular, airway, and pleural architecture, with minimal focal edema on day 1 that resolved by day 35. Allografts displayed a vascular driven pathology. On day 1, the multifocal intra-alveolar edema was accompanied by subtle perivascular changes (Figure 3). By day 7, dense cuff-like perivascular infiltrates dominated, composed of lymphocytes and histiocytes, frequently accompanied by pleural extension. On day 35, these infiltrates became organized and shifted in composition toward plasma cell–rich cuffs with fewer lymphocytes, paralleled by progressive perivascular fibrosis (Figure 3; Supplementary Figure S3). The parenchyma showed edema on day 1 that was largely resolved in isografts but more extensive in allografts (Figure 3). The airways remained preserved in isografts, whereas allografts demonstrated peribronchial fibrosis on day 35 (Figure 3; Supplementary Figure S3).
FIGURE 3
Bulk spatial GeoMx proteomics confirm vessel orientated rejection
2D-UMAP projection of all included samples, based on normalized protein expression, showed isografts at all timepoints clustering together with allograft day 1, while a second cluster consisted of allograft day 7 and 35 together (Supplementary Figure S4). Indeed, no significant differences were found comparing isografts and allografts on day 1. To continue, based on this clustering, data of allograft day 7 and 35 were pooled and this was compared to isograft (Figure 4C). We observed a higher expression of most of the proteins at the arterial regions in allografts compared to isografts. Particularly CD45 (pan leukocyte marker), CD4 and CD44 (T cell markers) were more present around the arteries of allografts compared to isografts (p < 0.001; p = 0.001; p = 0.002, respectively) and the other structures. Expression of CD19 (B cell marker) was also higher around the arteries of allografts compared to isografts (p = 0.006). Both innate and adaptive immune markers were more expressed in the allografts’ arteries compared to isografts and other structures. The expression of CD11c (dendritic cell marker) was increased in allografts in the parenchyma (p = 0.001) and around the arteries (p = 0.006) and airways (p = 0.029), compared to isografts. Protein expression of immune regulatory markers CD39, CD163, PD-L1 and FoxP3 decreased in the allograft’s parenchyma (p < 0.001; p < 0.001; p < 0.001; p < 0.001) and veins (p < 0.001; p < 0.001; p < 0.001; p = 0.002) compared to isografts. Expression of other adaptive immunity markers CD27, ICOS and CD40L (T cell activation markers) were significantly less detected in allografts’ vein compared to isografts (p < 0.001; p < 0.001; p < 0.001). EpCAM and PanCK (epithelial markers) were decreased in allografts’ airways (p < 0.001; p < 0.001) in comparison to isografts. Smooth muscle actin (SMA) was significantly lower in allografts’ airways (p < 0.001) and arteries (p < 0.001). However, fibronectin seemed to be higher in allografts, although not significant. Ki-67 (proliferation marker) was expressed more highly around the arteries and parenchyma of the allografts compared to the other structures, and to isografts (p = 0.023; p = 0.023).
FIGURE 4
3D visualization by ex vivo µCT imaging to locate rejection
Ex vivo µCT was performed on donor controls, isografts, and allografts on day 1, 7, and 35, enabled segmentation of the airway tree, pulmonary arteries, and veins (Figure 5). Donor control lungs (Balb/c and C57BL/6) and isografts demonstrated preserved morphology of both airway and vascular compartments at all time points. Allografts appeared normal on day 1, with intact airways and vasculature comparable to controls. In contrast, by day 35, allografts exhibited marked vascular pathology: arteries and veins were narrowed, collapsed, and embedded within dense fibrotic regions. Notably, airway structures appeared intact without evidence of distortion or malformation and airway obstruction was only found secondary and in severe stages of rejection.
FIGURE 5
Discussion
This study demonstrated a vasculature and endothelial cell centered onset of rejection in a mouse model of chronic lung rejection. Our findings indicate that endothelial cells are key in the immune response of acute rejection (day 7) and the fibrotic remodeling of chronic rejection (day 35) around the arteries, analogous to chronic vascular injury observed in other solid organ transplants.
We previously characterized the alloimmune response in this mouse transplant model with both innate and adaptive immune cells involved [22, ]. In this study we confirmed the temporal relation between immune activation on day 7 (with little fibrosis), and organized fibrosis on day 35 in the allografts. As such, day 7 may be regarded as acute rejection and day 35 may be regarded as chronic rejection. In the current study, we observed an early (day 7) loss of endothelial cells that coincided with this inflammatory peak of the acute phase of rejection. Moreover, early lymphocyte activation and leukocyte adhesion to the endothelium was enhanced by endothelial cell activation with increased expression of VCAM-1 and increased expression of ICAM-1 in leukocytes. The parallel increase in MHC molecules on endothelial cells (especially MHCII) further supports the leading role for the endothelium in orchestrating the immune response, by activating the adaptive immune system with its antigen-presenting capacity, directly activating T-helper cells. Jambusaria et al. demonstrated that lung endothelial cells express higher levels of immunoregulatory genes than those from the brain or heart, underscoring their unique immune-modulating capacity []. Carmeliet et al. further introduced the concept of “immunomodulatory endothelial cells” (IMECs), a specialized subset capable of engaging directly in immune signaling []. Histology and especially the spatial proteomics highlighted compartment-specific patterns of injury, confirming innate and adaptive organization specifically localized to arteries. Our findings suggest that endothelial cells function not only as targets of the immune response but also as active contributors to ongoing alloimmune activation within the graft. Since the alloantigenic stimulus persists, the immune response cannot resolve. Continuous endothelial activation and injury likely promote maladaptive repair processes, further resulting in perivascular fibrosis. This remodeling extends into adjacent parenchymal, pleural and airway compartments, providing a mechanistic link between early vascular injury and later graft dysfunction. Collectively, these observations position the vasculature as a central orchestrator of both the inflammatory and fibrotic phases of chronic lung allograft injury.
How this concept of a vessel-centered rejection process translates to the human clinical setting with CLAD (both RAS and BOS) is highly relevant. Acute rejection is histologically confirmed by the presence of lymphocytic perivascular infiltrates in the clinic, which is in parallel with the acute rejection observed on day 7 in the mice. In addition, some human studies highlighted the vascular contribution to CLAD, with reduced VEGF levels in bronchoalveolar lavage fluid of patients with RAS [], and vascular remodeling, including pulmonary arteriopathy and venopathy, has been described in both BOS and RAS [, ].
To extend our molecular and histological observations to the macrostructural level, a high-resolution ex vivo microCT was used enabling segmentation of the airway tree and pulmonary vascular compartments. Consistent with our molecular findings, vascular structures in allografts showed progressive narrowing, pruning, and remodeling by day 35, while airway architecture remained largely preserved, providing 3D confirmation that chronic graft injury in this mouse model primarily affects the vasculature rather than the airway compartment. This vascular involvement suggests that vascular injury precedes or drives subsequent parenchymal or airway remodeling which ultimately leads to loss of lung function. From a translational perspective, this also implies that, potentially, in vivo imaging approaches focusing on vascular alterations could serve as sensitive tools to detect CLAD progression in clinical lung transplantation at earlier stages, prior to the decrease in lung function. Advanced imaging modalities capable of assessing pulmonary perfusion or vascular integrity may therefore hold diagnostic value in identifying patients at risk for CLAD. Moreover, therapeutic strategies aimed at preserving endothelial integrity or modulating vascular inflammation might complement current immunosuppressive regimens and offer a new avenue to improve long-term graft survival.
The limitations of this study are the nature of the model, being a highly controlled setting to dissect alloimmune mechanisms, which cannot fully replicate the complexity of human lung transplantation. Immunosuppressive regimens also differ between species; we used cyclosporine A to moderate the immune response, whereas tacrolimus-based therapy is standard in clinical practice, which may affect the kinetics and nature of rejection. Only male mice were included, and potential sex-specific immune differences therefore remain unexplored. Furthermore, the number of isograft controls was limited, which may reduce statistical power for some comparisons. Mechanistic pathways underlying endothelial cell loss and its direct effects on fibroblast activation and fibrosis were not addressed here but represent an important direction for future work. Finally, this is a murine model of rejection, not CLAD. Modeling CLAD requires graft-specific functional measurements, and their absence is a limitation of this study. Despite these limitations, this model and study offered a unique opportunity to study rejection processes under reproducible and isolated conditions, minimizing confounding factors that are difficult to control in human transplantation.
In conclusion, the current study reveals that chronic vascular injury is a central feature of rejection also in the lung, preceding and potentially driving fibrotic remodeling. By identifying the endothelium as an active participant rather than a passive target, this study redefines lung rejection as a vessel-centered process. Targeting vascular inflammation and preserving endothelial integrity may thus hold the key to improving long-term lung allograft survival. This study invites to investigate how vascular rejection fits into the concept of CLAD.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by Ethical Committee for Animal Experimentation (ECD) of KU Leuven (Katholieke Universiteit Leuven) (P194/2019). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
Conceptualization: JK, CH, PEVdS, RV, and BMV; Transplant procedure: JK, CH, XJ, TH, HL, JVS, MM, MZ-O, CV, AOY, AB, LJC, and BMV; µCT imaging and visualization: JK, CH, XJ, AV, MM, MZ-O, YM, BT, and BMV; Flow cytometry: JK, CH, XJ, EP, FP, PEVdS, and BMV; Spatial GeoMx: JK, CH, PK, JRQ, FP, BMV, BO and LV; Histology: JK, CH, XJ, AV, GA, TH, MM, VG, AV, AOY, BMV, AS, and LV; writing – original draft preparation: JK, CH, AZ, RV, and BMV; writing – review and editing: JK, CH, XJ, TH, HL, EP, FP, HB, PK, AV, JV, MM, MZ-O, AZ, GA, VG, AV, CV, YM, BT, AOY, SB, LDS, LG, AB, JRQ, FP, PEVdS, LJC, RV, BMV, LV, BO, and AS. All authors contributed to the article and approved the submitted version.
Funding
The author(s) declared that financial support was received for this work and/or its publication. JK (1198920N), CH (1152225N), XJ (11PGT24N), PK (1120425N) and VG (11L9822N) are sponsored by a fundamental research grant from the FWO. EP is a recipient of the L’Oréal-Unesco Women for Sciences PhD fellowship. RV is supported by FWO as a senior clinical researcher (1803526N) and by a research grant (G060322N). PEVdS is funded by KU Leuven (C14/23/143) and an FWO research grant (G066723N). LJC is supported by an FWO research grant (G090922N). BMV is funded by KU Leuven (C16/19/005) and by an FWO research grant (G098025N).
Conflict of interest
The authors(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontierspartnerships.org/articles/10.3389/ti.2026.16293/full#supplementary-material
References
1.
VerledenGMGlanvilleARLeaseEDFisherAJCalabreseFCorrisPAet alChronic lung allograft dysfunction: definition, diagnostic criteria and approaches to treatment. A consensus report from the pulmonary council of the ISHLT. J Heart Lung Transpl (2019) 38(5):493–503. 10.1016/j.healun.2019.03.009
2.
GrahamCNWatsonCBarlevAStevensonMDharnidharkaVR. Mean lifetime survival estimates following solid organ transplantation in the US and UK. J Med Econ (2022) 25(1):230–7. 10.1080/13696998.2022.2033050
3.
ZajacovaADupontLJDe LeynPCeulemansLJVosR, Leuven Lung Transplant Group. Characteristics and outcomes of 1500 lung transplantations in the leuven lung transplant program: turning past lessons into tomorrow’s foundations. Transpl Int (2025) 38:15495. 10.3389/ti.2025.15495
4.
ChambersDCCherikhWSHarhayMOHayesDHsichEKhushKKet alThe international thoracic organ transplant registry of the international society for heart and lung transplantation: Thirty-sixth adult lung and heart–lung transplantation Report—2019; focus theme: donor and recipient size match. J Heart Lung Transpl (2019) 38(10):1042–55. 10.1016/j.healun.2019.08.001
5.
GauthierJMRuiz-PérezDLiWHachemRRPuriVGelmanAEet alDiagnosis, pathophysiology and experimental models of chronic lung allograft rejection. Transplantation (2018) 102(9):1459–66. 10.1097/TP.0000000000002250
6.
BosSVosRVan RaemdonckDEVerledenGM. Survival in adult lung transplantation: where are we in 2020?Curr Opin Organ Transpl (2020) 25(3):268–73. 10.1097/MOT.0000000000000753
7.
VanstapelAVerledenSVerbekenEBraubachPGoosTDe SadeleerLet alBeyond bronchiolitis obliterans: in-depth histopathologic characterization of bronchiolitis obliterans syndrome after lung transplantation. J Clin Med (2021) 11(1):111. 10.3390/jcm11010111
8.
MeyerKCRaghuGVerledenGMCorrisPAAuroraPWilsonKCet alAn international ISHLT/ATS/ERS clinical practice guideline: diagnosis and management of bronchiolitis Obliterans syndrome. Eur Respir J (2014) 44(6):1479–503. 10.1183/09031936.00107514
9.
VerledenSEVosRVanaudenaerdeBMVerledenGM. Chronic lung allograft dysfunction phenotypes and treatment. J Thorac Dis (2017) 9(8):2650–9. 10.21037/jtd.2017.07.81
10.
MitchellRNLibbyP. Vascular remodeling in transplant vasculopathy. Circ Res (2007) 100(7):967–78. 10.1161/01.RES.0000261982.76892.09
11.
LibbyPPoberJChronic rejection, 14 (2021). p. 387–97.
12.
NaesensMRoufosseCHaasMLefaucheurCMannonRBAdamBAet alThe banff 2022 kidney Meeting Report: reappraisal of microvascular inflammation and the role of biopsy-based transcript diagnostics. Am J Transpl (2024) 24(3):338–49. 10.1016/j.ajt.2023.10.016
13.
KlocMGhobrialRM. Chronic allograft rejection: a significant hurdle to transplant success. Burns Trauma (2014) 2(1):2321–3868. 10.4103/2321-3868.121646
14.
MehraMRCrespo-LeiroMGDipchandAEnsmingerSMHiemannNEKobashigawaJAet alInternational society for heart and Lung transplantation working formulation of a standardized nomenclature for cardiac allograft vasculopathy—2010. J Heart Lung Transpl (2010) 29(7):717–27. 10.1016/j.healun.2010.05.017
15.
BosSMilrossLFilbyAJVosRFisherAJ. Immune processes in the pathogenesis of chronic lung allograft dysfunction: identifying the missing pieces of the puzzle. Eur Respir Rev (2022) 31(165):220060. 10.1183/16000617.0060-2022
16.
BosSPradèrePBeeckmansHZajacovaAVanaudenaerdeBMFisherAJet alLymphocyte depleting and modulating therapies for chronic lung allograft dysfunction. Pharmacol Rev (2023) 75(6):1200–17. 10.1124/pharmrev.123.000834
17.
AunerSHillebrandCBoehmPMBoeckerJKorenDSchwarzSet alImpact of transient and persistent donor-specific antibodies in lung transplantation. Transpl Int (2024) 37:12774. 10.3389/ti.2024.12774
18.
LoupyALefaucheurC. Antibody-mediated rejection of solid-organ allografts. N Engl J Med (2018) 379(12):1150–60. 10.1056/NEJMra1802677
19.
BeryAIHachemRR. Antibody-mediated rejection after lung transplantation. Ann Transl Med (2020) 8(6):411. 10.21037/atm.2019.11.86
20.
BeeckmansHBosSVosRGlanvilleAR. Acute rejection and chronic lung allograft dysfunction. Clin Chest Med (2023) 44(1):137–57. 10.1016/j.ccm.2022.10.011
21.
LamaVNBelperioJAChristieJDEl-ChemalySFishbeinMCGelmanAEet alModels of lung transplant research: a consensus statement from the national heart, lung, and blood institute workshop. JCI Insight (2017) 2(9), e93121. 10.1172/jci.insight.93121
22.
HeiglTKaesJAelbrechtCSerréJYamadaYGeudensVet alThe nature of chronic rejection after lung transplantation: a murine orthotopic lung transplant study. Front Immunol (2024) 15:1369536. 10.3389/fimmu.2024.1369536
23.
ChiangCFChaoTTSuYFHsuCCChienCYChiuKCet alMetformin-treated cancer cells modulate macrophage polarization through AMPK-NF-κB signaling. Oncotarget (2017) 8(13):20706–18. 10.18632/oncotarget.14982
24.
JungraithmayrWDraenertAMarquardtKWederW. Ultrastructural changes in acute lung allograft rejection: novel insights from an animal study. The J Heart Lung Transplantation: The Official Publication Int Soc Heart Transplant (2012) 31(1):94–100. 10.1016/j.healun.2011.10.003
25.
KaesJPollenusEHooftCLiuHAelbrechtCCambierSet alThe immunopathology of pulmonary rejection after murine lung transplantation. Cells (2024) 13(3):241. 10.3390/cells13030241
26.
JinXKaesJVan SlambrouckJInciIArniSGeudensVet alA comprehensive review on the surgical aspect of lung transplant models in mice and rats. Cells (2022) 11(3):480. 10.3390/cells11030480
27.
PollenusEPhamTTVandermostenLPossemiersHKnoopsSOpdenakkerGet alCCR2 is dispensable for disease resolution but required for the restoration of leukocyte homeostasis upon experimental malaria-associated acute respiratory distress syndrome. Front Immunol (2021) 11:628643. 10.3389/fimmu.2020.628643
28.
BerghenNDekosterKMarienEDabinJHillenAWoutersJet alRadiosafe micro-computed tomography for longitudinal evaluation of murine disease models. Sci Rep (2019) 9(1):1–10. 10.1038/s41598-019-53876-x
29.
Vande VeldeGPoelmansJDe LangheEHillenAVanoirbeekJHimmelreichUet alLongitudinal micro-CT provides biomarkers of lung disease that can be used to assess the effect of therapy in preclinical mouse models, and reveal compensatory changes in lung volume. Dis Model Mech (2016) 9(1):91–8. 10.1242/dmm.020321
30.
YushkevichPAPivenJHazlettHCSmithRGHoSGeeJCet alUser-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. NeuroImage (2006) 31(3):1116–28. 10.1016/j.neuroimage.2006.01.015
31.
JambusariaAHongZZhangLSrivastavaSJanaATothPTet alEndothelial heterogeneity across distinct vascular beds during homeostasis and inflammation. eLife (2020) 9:e51413. 10.7554/eLife.51413
32.
AmersfoortJEelenGCarmelietP. Immunomodulation by endothelial cells — partnering up with the immune system?Nat Rev Immunol (2022) 22(9):576–88. 10.1038/s41577-022-00694-4
33.
VanstapelAGoldschmedingRBroekhuizenRNguyenTSacreasAKaesJet alConnective tissue growth factor is overexpressed in explant lung tissue and broncho-alveolar lavage in transplant-related pulmonary fibrosis. Front Immunol (2021) 12:661761. 10.3389/fimmu.2021.661761
Summary
Keywords
blood vessels, chronic lung allograft dysfunction (CLAD), chronic rejection, flow cytometry, lung transplantation, murine orthotopic left lung transplantation, spatial proteomics
Citation
Kaes J, Hooft C, Jin X, Heigl T, Liu H, Pollenus E, Prenen F, Beeckmans H, Kerckhof P, Vanstapel A, Van Slambrouck J, Mercurio M, Zapata-Ortega M, Aerts G, Geudens V, Vermaut A, Vanvuchelen L, Somers A, Verbeylen L, Vanluyten C, Barbarossa A, Mohamady Y, Özsoy B, Zajacova A, Tielemans B, Yildirim AO, Bos S, De Sadeleer L, Godinas L, Rojas Quintero J, Polverino F, Van den Steen PE, Ceulemans LJ, Vos R and Vanaudenaerde BM (2026) Blood vessels as primary site of rejection in murine lung transplantation. Transpl. Int. 39:16293. doi: 10.3389/ti.2026.16293
Received
23 January 2026
Revised
18 April 2026
Accepted
01 May 2026
Published
29 May 2026
Volume
39 - 2026
Updates
Copyright
© 2026 Kaes, Hooft, Jin, Heigl, Liu, Pollenus, Prenen, Beeckmans, Kerckhof, Vanstapel, Van Slambrouck, Mercurio, Zapata-Ortega, Aerts, Geudens, Vermaut, Vanvuchelen, Somers, Verbeylen, Vanluyten, Barbarossa, Mohamady, Özsoy, Zajacova, Tielemans, Yildirim, Bos, De Sadeleer, Godinas, Rojas Quintero, Polverino, Van den Steen, Ceulemans, Vos and Vanaudenaerde.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bart M. Vanaudenaerde, bart.vanaudenaerde@kuleuven.be
‡These authors have contributed equally to this work
ORCID: Janne Kaes, orcid.org/0000-0002-3326-3769; Charlotte Hooft, orcid.org/0000-0002-8471-3686; Xin Jin, orcid.org/0000-0002-3295-7179; Tobias Heigl orcid.org/0000-0002-9591-8345; Emilie Pollenus, orcid.org/0000-0002-2706-6259; Fran Prenen, orcid.org/0000-0002-5431-1715; Hanne Beeckmans, orcid.org/0000-0002-6176-0404; Pieterjan Kerckhof, orcid.org/0000-0002-3806-4478; Arno Vanstapel, orcid.org/0000-0003-2412-8001; Jan Van Slambrouck, orcid.org/0000-0002-7069-1535; Martina Mercurio, orcid.org/0009-0004-9282-3430; Marta Zapata-Ortega, orcid.org/0000-0003-0379-358X; Gitte Aerts, orcid.org/0000-0001-9348-8081; Vincent Geudens, orcid.org/0000-0002-4386-1746; Astrid Vermaut, orcid.org/0000-0002-4145-491X; Lise Vanvuchelen, orcid.org/0009-0001-8346-9693; Anse Somers, orcid.org/0009-0008-7441-3656; Lara Verbeylen, orcid.org/0009-0007-8060-0207; Cedric Vanluyten, orcid.org/0000-0002-0128-3351; Annalisa Barbarossa, orcid.org/0000-0002-2358-7284; Yousry Mohamady, orcid.org/0000-0003-0103-2517; Balin Özsoy, orcid.org/0000-0002-3976-2286; Andrea Zajacova, orcid.org/0000-0002-9691-2500; Birger Tielemans, orcid.org/0000-0002-1775-1797; Ali Onder Yildirim, orcid.org/0000-0003-1969-480X; Saskia Bos, orcid.org/0000-0002-5336-5914; Laurens De Sadeleer, orcid.org/0000-0003-1237-7369; Laurent Godinas, orcid.org/0000-0003-2214-5879; Joselyn Rojas Quintero, orcid.org/0000-0003-4994-075X; Francesca Polverino, orcid.org/0000-0001-9686-5698; Philippe E. Van den Steen, orcid.org/0000-0002-7334-9145; Laurens J. Ceulemans, orcid.org/0000-0002-4261-7100; Robin Vos, orcid.org/0000-0002-3468-9251; Bart M. Vanaudenaerde, orcid.org/0000-0001-6435-6901
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.