Skip to main content

Native T1 mapping detects both acute clinical rejection and graft dysfunction in pediatric heart transplant patients

Abstract

Background

Cardiovascular magnetic resonance (CMR) is emerging as an important tool for cardiac allograft assessment. Native T1 mapping may add value in identifying rejection and in assessing graft dysfunction and myocardial fibrosis burden. We hypothesized that CMR native T1 values and features of textural analysis of T1 maps would identify acute rejection, and in a secondary analysis, correlate with markers of graft dysfunction, and with fibrosis percentage from endomyocardial biopsy (EMB).

Methods

Fifty cases with simultaneous EMB, right heart catheterization, and 1.5 T CMR with breath-held T1 mapping via modified Look-Locker inversion recovery (MOLLI) in 8 short-axis slices and subsequent quantification of mean and peak native T1 values, were performed on 24 pediatric subjects. A single mid-ventricular slice was used for image texture analysis using nine gray-level co-occurrence matrix features. Digital quantification of Masson trichrome stained EMB samples established degree of fibrosis. Markers of graft dysfunction, including serum brain natriuretic peptide levels and hemodynamic measurements from echocardiography, catheterization, and CMR were collated. Subjects were divided into three groups based on degree of rejection: acute rejection requiring new therapy, mild rejection requiring increased ongoing therapy, and no rejection with no change in treatment. Statistical analysis included student’s t-test and linear regression.

Results

Peak and mean T1 values were significantly associated with acute rejection, with a monotonic trend observed with increased grade of rejection. Texture analysis demonstrated greater spatial heterogeneity in T1 values, as demonstrated by energy, entropy, and variance, in cases requiring treatment. Interestingly, 2 subjects who required increased therapy despite low grade EMB results had abnormal peak T1 values. Peak T1 values also correlated with increased BNP, right-sided filling pressures, and capillary wedge pressures. There was no difference in histopathological fibrosis percentage among the 3 groups; histopathological fibrosis did not correlate with T1 values or markers of graft dysfunction.

Conclusion

In pediatric heart transplant patients, native T1 values identify acute rejection requiring treatment and may identify graft dysfunction. CMR shows promise as an important tool for evaluation of cardiac grafts in children, with T1 imaging outperforming biopsy findings in the assessment of rejection.

Introduction

Heart transplant remains a life-saving intervention for pediatric patients with advanced heart failure. Despite significant advancements in immunosuppressive therapy, cardiac allograft rejection remains an important cause of mortality and morbidity after transplantation [1, 2]. Allograft rejection and graft failure are often clinically insidious, yet early detection significantly improves outcomes [3]. Identification of allograft rejection is complex; the current gold standard, direct pathologic evaluation of the myocardium via endomyocardial biopsy (EMB), frequently results in false negatives [4,5,6] due to random sampling of the accessible regions of the right ventricle [7]. Further, EMB requires an invasive procedure for tissue procurement, which while generally safe, involves a risk of cardiac perforation, arrhythmia, effusion, and tricuspid valve injury [5, 8, 9]. Other noninvasive methods of identifying rejection, including echocardiography [10, 11], serum markers [12], and clinical signs and symptoms also carry limitations in detecting rejection but have been used for diagnosis of late graft dysfunction. Characteristics of late graft dysfunction can be diagnosed by hemodynamic measurements through echocardiography and cardiac catheterization [13, 14], and myocardial fibrosis noted on histological evaluation [15].

Cardiovascular magnetic resonance (CMR) offers a diagnostic advantage in detection of myocardial pathology due to its ability to characterize the entire myocardium for evidence of fibrosis or edema using quantitative T1 and T2 weighted techniques [16]. Native T1 values increase with processes that result in myocardial edema and fibrosis [17]. In adult heart transplant recipients, T1 and T2 mapping techniques have been shown to be more sensitive for the detection of acute allograft rejection than EMB, with excellent negative predictive capacity [18,19,20,21].

T1 and T2 mapping techniques are complex pulse sequences which require systematic use and careful interpretation to yield accurate data, particularly in the pediatric population [22]. Typical spatial and temporal resolution values favor larger hearts with slower heart rates to minimize cardiac motion artifacts and provide enough time for T1 recovery [23], making the application of parametric mapping in pediatric populations challenging. However, the use of parametric mapping has been successfully demonstrated in other conditions in the pediatric population [22], including myocarditis [24], hypertrophic cardiomyopathy [25], Duchenne muscular dystrophy myocardial disease [26], iron deposition cardiomyopathy [27], and anthracycline cardiotoxicity [28].

Another frequent challenge of using CMR parametric mapping is the variability of native T1 values, in particular between different field strengths [29], vendors, and individual scanners [16]. This is compounded by lack of normative pediatric T1 data to understand and identify disease states. For this reason, there is a consensus to use institutional reference ranges [16, 17], making comparability of measurements across patients and institutions challenging [23]. We hypothesized that in addition to changes in native T1 value, T1 maps from patients with rejection would demonstrate changes in spatial variability of T1 voxel intensities. Texture analysis is a computational imaging analysis method that quantifies spatial heterogeneity using the relative intensity differences between neighboring voxels [30, 31], independent of the actual voxel intensities. Texture analysis thus overcomes the lack of normative values available and has been demonstrated to have utility in other conditions, including myocarditis [32, 33], dilated cardiomyopathy [34], and hypertrophic cardiomyopathy [35].

There is a clear need for a non-invasive and more accurate method of detecting acute allograft rejection and late graft dysfunction in pediatric heart transplant patients. The overall aim of this study was to determine the utility of CMR native T1 values in (1) identifying acute graft rejection and (2) assessing graft dysfunction in pediatric heart transplant patients. In secondary analysis, we aimed to compare textural features of parametric maps in cases of rejection versus those without rejection and to assess the relationship of native T1 mapping and histological myocardial fibrosis burden.

Methods

Clinical procedure

In this IRB-approved study, with consent/assent as appropriate, heart transplant patients referred for clinically-indicated EMB were prospectively enrolled to undergo noncontrast-CMR at 1.5 T (MAGNETOM Aera, Siemens Healthineers, Erlangen, Germany) followed by cardiac catheterization, with EMB in an adjoining biplane fluoroscopy suite. CMR included right heart catheterization, standard volumetry cines, phase contrast imaging, and breath-held native T1 mapping using Modified Look-Locker Inversion recovery (MOLLI) in eight short-axis slices. The following parameters were used for MOLLI T1 map acquisition: field of view 360 mm × 307 mm, percent phase field of view 50–80% based on patient’s body habitus, and matrix size 1.4 × 1.4 mm. For patients with a RR interval > 700, echo time (TE) was 1.12 ms, repetition time (TR) was 2.80 ms, and flip angle of 35 degrees; for patients with RR interval < 700, TE was 1.06 ms, TR was 1.93 ms, and flip angle of 70 degrees. Slice thickness and slice skip also varied based on patient’s size, ranging from 4–8 mm to 0–2 mm respectively. Measurements of volume, function, and cardiac output were performed using standard offline software (Medis Medical Imaging, AJ Leiden, Netherlands). Catheterization included collection of standard oximetry and hemodynamic data and EMB. If clinically indicated coronary angiography was performed according to surveillance standard.

Clinical data were recorded, including patient demographics, transplant history, rejection history, and serum brain natriuretic peptide (BNP) levels. Hemodynamic measurements from echocardiography (left ventricular (LV) ejection fraction, mitral E/e’), catheterization (right atrial mean pressure, right ventricle (RV) systolic pressure, RV end diastolic pressure (RVEDP), main pulmonary artery mean pressure, and pulmonary capillary wedge pressure), and CMR (LV and RV end diastolic volume) were also included. Serum BNP and hemodynamic measurements listed above were considered markers of graft dysfunction for the purpose of analysis.

Rejection analysis

Evidence of cellular rejection and antibody mediated rejection on EMB were graded based on the International Society of Heart and Lung Transplantation (ISHLT) guidelines [36] by a surgical pathologist following standard clinical practices. Clinically, rejection was defined based on treatment plan created by transplant team following the catheterization/EMB procedure. Of note, CMR T1 measurements and analyses did not affect treatment decision, as these are not reported with standard clinical data at our institution. Per our institutional protocol, those cases with Grade 2R or above acute cellular rejection on biopsy received new rejection treatment including intravenous steroids, thymoglobulin, immunoglobulins, etc.. Those cases with Grade 0R or 1R acute cellular rejection on biopsy but no hemodynamic compromise did not receive any treatment or modification to immunosuppressive therapy. Cases with Grade 0R or 1R rejection, with hemodynamic compromise, on echocardiogram, catheterization, or CMR, received treatment per transplant team using other data including biopsy, history of rejection, clinical symptoms, and donor specific antibody (DSA) results, without impact of CMR T1 mapping data. Cases were divided into one of three outcome groups (A, B, C), based on their degree of rejection therapy recommended by their physicians, who were blinded to T1 mapping data. Group A included cases with no rejection and thus no changes made to their treatment regimen. Group B included cases with some evidence of rejection which required augmentation of maintenance treatment regimen or initiation of oral steroids under pulse dosing. Group C included cases with significant evidence of acute rejection, who received new rejection treatment such as intravenous (IV) or oral steroids at pulse dosing, thymoglobulin, or intravenous immunoglobulin (IVIG).

Parametric map analysis

Native T1 parametric maps were deidentified and analyzed using additional offline software (OsiriX, Bernex, Switzerland). The middle 6 slices were used (2 basal, 2 mid, and 2 apical slices), to minimize through-plane motion artifacts associated with the most basal and apical slices. Based on the American Heart Association 17-segment model for myocardial segmentation [37], 16 total regions of interest (ROIs) were generated, as demonstrated in Fig. 1. Per lab standard [26], ROIs were traced in the septal and lateral walls only to avoid partial volume effect with fat and lungs in the anterior wall and partial volume effect with the diaphragm and stomach in the inferior wall commonly found in the pediatric population. ROIs were traced by a blinded reviewer using the “middle-third” technique to avoid artifacts and blood pool at the endocardial border, per lab standard [26, 38], yielding 16 segmental average regional voxel T1 values. Global mean and peak native T1 values were quantified using these 16 segmental T1 values.

Fig. 1
figure 1

Location representation of regions of interest (ROIs) based on the American Heart Association 17-segment model for myocardial segmentation. [37] Shading indicates percentage of native T1 values in that segment that were abnormal (> 1050). In basal slices, ROIs were generated on two septal segments each (segment 2 and 3) and a lateral wall region (combined segments 5 and 6). In mid slices, ROIs were generated on two septal segments each (segment 8 and 9) and a lateral wall region (combined segments 11 and 12). In apical slices, ROIs were generated on the septal (segment 14) and lateral (segment 16) segments. Of note, per lab standard, anterior and inferior segments were not analyzed to avoid partial volume effect

A second reviewer traced ROIs on 20% of cases (10 cases) to measure interobserver variability.

Given the myriad of factors contribution to local variation in T1 values, our institution created and maintains a local normal T1 database [16], which has both general norms (900–1050 ms) and normal values per body surface area (BSA) quartile. Per internal local studies on healthy control patients, normal T1 values were considered 900–1050 ms, with values above 1050 ms considered to be abnormal (Additional file 1: Table S2).

Texture analysis

LV myocardium was semiautomatically segmented using Otsu thresholding from T1 maps with offline software (Seg3D, Scientific Computing and Imaging Institute, Salt Lake City, Utah, USA). Blood pool and epicardial myocardium were carefully excluded. Image texture analysis (Fig. 2) was performed on a single mid-ventricular slice, using a previously-described texture toolbox [39] in MATLAB (The MathWorks Inc., Natick, Massachusetts, USA). Segmented maps were normalized to 256 intensity levels. A 256 × 256 Gy-level co-occurrence matrix (GLCM) was created, quantifying the frequency with which a pixel of a given intensity neighbors a pixel of another given pixel intensity. The GLCM is constructed such that each entry (i, j) in the GLCM quantifies the number of times that a pixel of intensity i neighbors a pixel of intensity j. An image with more homogeneous pixel intensity is thus more diagonally dominant than one with pixel heterogeneity. To quantify image heterogeneity, nine GLCM-based texture features were computed for each slice as previously described: [30, 40,41,42] energy, contrast, entropy, homogeneity, correlation, sum average, variance, dissimilarity, and autocorrelation (Additional file 1: Table S1).

Fig. 2
figure 2

Demonstrated is the clinical treatment algorithm at our institution. Patients undergo transthoracic echocardiogram, right heart catheterization, and endomyocardial biopsy at each surveillance encounter. Coronary angiography is performed if indicated. Of note, the 50 cases included this study also underwent CMR per research protocol and right heart catheterization was performed under CMR guidance

EMB fibrosis percentage analysis

EMB specimens underwent staining with Masson’s trichrome following a standard laboratory staining protocol. Whole slide scanning was performed using Tissue Scope (Huron Technologies, Ontario, Canada). Using ImageScope (Aperio, Vista, California, USA), endocardial collagen and artifacts were manually excluded. A pixel counting algorithm was developed by setting a colorimetric threshold specific to blue stained collagen. Using this algorithm, percentage tissue fibrosis was calculated per biopsy. Surgical pathologists served as expert readers, with involvement in the creation and tuning of the percentage fibrosis tool.

Statistical analysis

Interclass correlation coefficient was used to verify interobserver variability. Student’s t-test was performed on native T1 values and parametric map texture features of cases with and without rejection, as defined by clinical rejection (Group A vs Groups B + C above) and by biopsy grade (Grade 0R vs Grade 1R or greater); ANOVA analysis was also performed across the three treatment groups. A receiver operating characteristic curve was created for mean T1 values, peak T1 values, and significant texture analysis features. Pearson correlation was performed on T1 mean and peak values against markers of graft dysfunction, which included serum BNP and hemodynamic data from echocardiography, catheterization, and CMR. All statistical analyses were performed using GraphPad Prism version 9.0.0 for Windows (GraphPad Software, San Diego, California, USA).

Results

Patient cohort

Twenty-four pediatric patients in 50 cases (12.2 ± 4.6 years, graft age 5.3 ± 4.1 years, 36% female) underwent study procedures for surveillance (70%), suspected rejection (4%), and follow-up of prior rejection (26%). Thirty-seven cases (74%) were in Group A (no rejection, no therapy changes), 6 cases (12%) in Group B (mild rejection, minor therapy changes), and 7 cases (14%) in Group C (rejection present, major therapy changes). For analysis purposes, groups B + C represented those cases with clinical rejection. There was no difference in mean age, graft age, or graft ischemia time between Group A and Group B + C (Table 1). LV ejection fraction on echocardiogram was lower in Group B + C than Group A (60.8 ± 3.4%vs 64.3 ± 5.2%, p = 0.009), though both were normal (Table 2). A higher percentage of cases in Group B + C, compared to Group A, were positive for DSAs (53.8% vs 20%, p = 0.051).

Table 1 Patient cohort: demographics
Table 2 Indication for treatment

Similarly, a higher percentage of cases in Group B + C had biopsies with cellular rejection grade > 0R compared to Group A (84.6% vs 16.2%, p = 0.001).

Two patients in Group B + C who had Grade 0R biopsies and abnormal hemodynamics, prompted modification of immunosuppressive therapy (Table 2). Of the six cases in Group A with Grade 1R biopsies, all had normal hemodynamics on echocardiogram, catheterization, and CMR; further all six cases were receiving more frequent surveillance due to history of rejection and continued to show improvement during these cases. There were no cases of antibody mediated rejection among the 50 cases.

T1 parametric map analysis

Intraclass correlation coefficient of 20% of cases, with regions of interest traced by 2 reviewers, was 0.829 for global mean T1 and 0.830 for peak T1.

T1 correlation with clinical rejection

A monotonic, increasing trend was noted in both mean and peak T1 values, with increasing degree of rejection (Figs. 3, 4), with peak T1 values in the abnormal range (T1 values greater than 1050) for Group B and C. Area under the curve (AUC) of receiver operating characteristic curve of mean T1 values was 0.746 (p = 0.007) and of peak T1 values was 0.730 (p = 0.012) (Fig. 5). Notably, ROC analysis demonstrated 100% sensitivity at peak T1 values > 1050 ms, which is the cutoff of normal vs abnormal T1 values at our institution based on internal studies.

Fig. 3
figure 3

Representation of texture analysis of two example images that have the same pixel intensity mean, standard deviation, and distribution but differ in texture features. Image A has higher energy with lower entropy and variance than Image B

Fig. 4
figure 4

ANOVA analysis between three groups of global mean T1 and peak T1 values yielded p value of 0.0005 and 0.033 respectively. T-test analysis between group A (no change in treatment) versus combined groups B (maintenance treatment augmentation) and C (new treatment initiated) of global mean T1 and peak T2 values yielded pvalue of 0.007 and 0.016 respectively

Fig. 5
figure 5

Receive Operative Characteristics Curve demonstrating sensitivity and specificity of mean and peak T1 values between groups with and without clinical rejection requiring treatment. AUC for mean T1 values was 0.746 (p = 0.007) and for peak T1 values was 0.730 (p = 0.012)

Interestingly, T1 values did not differ between cases with Grade 0R and those with Grade 1R or higher cellular rejection [mean T1: 1018 ± 28 ms vs 1029 ± 45 ms respectively (p = 0.388); peak T1: 1067 ± 44 ms vs 1075 ± 44 ms respectively (p = 0.547)]. Of note, 2 of 13 cases (15%) in group B + C had biopsies with grade 0R rejection, were treated for clinical rejection, and both were found to have elevated peak T1 values.

T1 correlation with graft function

Mean T1 values demonstrated moderate correlation with BNP (r = 0.59) and mean pulmonary artery pressure (r = 0.33). Peak T1 values demonstrated moderate correlation with BNP (r = 0.52), right atrial pressure (r = 0.40), mean pulmonary artery pressure (r = 0.46), RVEDP (r = 0.36), and average pulmonary capillary wedge pressure (r = 0.33). T1 values did not correlate with other hemodynamic markers (Table 3). When cases with clinical rejection were removed, leaving only those cases with possible graft dysfunction, peak T1 correlated with BNP moderately (r = 0.54) but with no other markers.

Table 3 Correlation of native T1 values and markers of graft dysfunction

T1 correlation with fibrosis percentage

Fibrosis percentage, calculated from random EMB of the RV, did not demonstrate a difference between cases of clinical rejection (Group B + C), versus cases without rejection (Group A) (8.2% ± 3.1 vs 7.4% ± 5.2, p = 0.703). Fibrosis percentage did not correlate with hemodynamic markers listed above or with mean or peak T1 values.

Native T1 texture analysis

Image texture analysis demonstrated that mid-ventricular short axis slices from Group A differed from groups B and C in three of nine texture features computed. Energy, a measure of image homogeneity, was higher in Group A than Groups B + C (p = 0.033). Entropy and variance, which reflect randomness and heterogeneity, were higher in Groups B + C than in Group A (p = 0.008 and p = 0.001 respectively). ROC analysis of these texture features demonstrated AUC of 0.750 (p = 0.016) for energy, 0.779 (p = 0.007) for entropy, and 0.831 (p = 0.002) for variance (Fig. 6). Other texture features did not differ between the groups.

Fig. 6
figure 6

Receive Operative Characteristics Curve demonstrating sensitivity and specificity of 3 textural features noted to show significant differences between groups with and without clinical rejection requiring treatment. Noted is an AUC of 0.750 (p = 0.016) for energy, 0.779 (p = 0.007) for entropy, and 0.831 (p = 0.002) for variance

Discussion

This prospective study demonstrates that mean and peak native T1 values correlate with both degree of clinical rejection and graft dysfunction in pediatric heart transplant patients. Further, no correlation between native T1 values and biopsy findings was demonstrated, indicating that native T1 mapping holds promise to identify clinically important changes in the heart transplant population, possibly outperforming random EMB of the RV. CMR is emerging as a useful tool for rejection screening in pediatric transplant patients, particularly when biopsy negative clinical rejection is present. In addition, texture analysis of parametric T1 maps identified additional promising trends for detection of rejection using imaging heterogeneity. Similar utility of native T1 values has been noted in adult studies, which demonstrated T1 and T2 mapping as useful tools for detecting rejection and therefore reducing the need for invasive EMB. Specifically, Imran et al. demonstrated an excellent negative predictive value of T1 mapping for cardiac allograft rejection [20]. Similarly, Vermes et al., Dolan et al., and Usman et al., have demonstrated higher T2 levels and higher extracellular volume fraction (ECV) in patients with acute rejection [18, 19, 43]. Pediatric studies have noted higher native T1 values in all pediatric heart transplant patients compared to controls without transplants [44] and higher native T2 values in cases of allograft rejection [38]. However, in other studies in which rejection was defined solely by biopsy grade, no difference in T1 values was demonstrated [45]. Using a larger population and whole-heart imaging, we have demonstrated a difference in T1 values between cases of rejection and cases without rejection, when rejection is defined clinically by necessity of treatment rather than exclusively by biopsy grade. Further, due to significant inter-center variability in the frequency of EMB for rejection surveillance, prior studies have investigated the yield of EMB and its impact on outcomes in the pediatric population [46, 47]. Studies investigating frequency of EMB both in the first year post-transplantation [47] when rejection risk is highest as well as long-term routine surveillance [46], demonstrated similar outcomes between high intensity and low intensity protocols. Adult studies have supported discontinuing routine EMB for long-term surveillance [48, 49].

Though texture analysis has been demonstrated in other pediatric conditions such as myocarditis [32, 33] and cardiomyopathies [34, 35], it has not previously been applied in pediatric heart transplant population and we therefore sought to do so. Texture analysis allows for the evaluation of heterogeneity in relative voxel intensities rather than using absolute voxel intensities [30, 31], offering the advantage of overcoming lack of normal T1 values between institutions given differences in field strengths, vendors, and individual scanners [16].

Diastolic dysfunction impacts prognosis of various etiologies of chronic heart disease [50], including graft failure [51]. Several studies have demonstrated a relationship between diastolic dysfunction and myocardial fibrosis as measured by histological analysis [15, 51] and late post-gadolinium myocardial enhancement CMR [50]. In our study, random EMB fibrosis percentage did not correlate with T1 values or with markers of graft dysfunction. Further, fibrosis analysis of random EMB also did not differ between cases of clinical rejection. It is unclear if this reflects a limitation in the ability of T1 values to assess for myocardial fibrosis burden; however, it is more likely that fibrosis percentage measurements from randomly sampled EMB tissue samples do not accurately reflect whole heart myocardial fibrosis burden. Non-invasive markers, such as BNP and echocardiography mitral E/e’, have also been shown to correlate with RVEDP and pulmonary capillary wedge pressure [52], markers of graft dysfunction. Though native T1 values showed moderate correlation with hemodynamic markers and BNP when patients in acute clinical rejection were included, this correlation was not found when looking only at patients without rejection. Perhaps, native T1 values, while showing promise in identifying clinical rejection and in assessing graft dysfunction in patients with rejection, may be limited in their ability to assess for more chronic graft changes that occur, including those mediated by fibrosis.

Limitations

There are limitations to this work; this is a single center study including a relatively small cohort which includes longitudinal, repeat encounters. We recognize that though this is a prospective study, the definition of clinical rejection is retrospectively based on decision to treat. However, these treatment groups are in line with the clinical treatment algorithm at our institution and therefore we find are a reliable outcome measure to compare. Also of note is that only septal and lateral segments were used for T1 value measurements. Given the nature of T1 mapping, we preferred to analyze only highly reliable data at the cost of excluding certain segments. In the pediatric population, particularly with patients as small as 15 kg, the anterior and inferior segments are subject to partial volume effect from epicardial fat and lungs for anterior segments and diaphragm and stomach for inferior segments, leading to unreliable data. Therefore, these segments were not included in the analysis, as is consistent with lab standards [26]. It is also important to note that very few patients underwent this procedure due to clinical concern for rejection, but rather for screening. Further, screening protocols and treatment protocols differ between institutions as standardized guidelines do not exist; therefore, further investigation at a multi-institutional level is required.

Despite these limitations, we find it promising that CMR serves as a noninvasive screening tool during surveillance encounters and may be used to identify those patients that may be at higher risk of rejection and therefore require further evaluation. Transplant rejection surveillance remains a multi-faceted approach, including assessment of clinical presentation, echocardiography, catheterization hemodynamics, and serum markers. Our team sought to demonstrate the possibility of a synergistic value of a combined CMR and EMB protocol for evaluating patients. In our population, two patients had Grade 0R biopsies but required treatment due to abnormal hemodynamics; these same patients demonstrated peak T1 values > 1050 ms, considered abnormal in per institutional normal values, further emphasizing the utility of CMR. ROC analysis demonstrates 100% sensitivity at peak T1 values > 1050 ms, demonstrating the possibility of a model in which patients requiring transplant rejection surveillance undergo CMR including parametric mapping, right heart catheterization, and hemodynamic measurements as initial screening. For those patients with abnormal findings, EMB would be performed and also inform treatment decision. In our analysis, if this model were used, all patients who required treatment would have underwent EMB and the appropriate treatment. Of the 37 patients who did not require treatment, 15 patients would have been saved from invasive EMB. Further, at our institution where CMR guided right heart catheterization is available and routinely performed for hemodynamic data, these patients are exposed to radiation only during coronary angiography, typically only performed at annual surveillance encounters. CMR, with its ability to perform radiation free evaluations that allow for hemodynamic assessment, particularly with the use of CMR guided right heart catheterization, and for entire myocardium assessment for fibrosis and edema, may be a promising tool in the pediatric heart transplant population. Further investigation is needed in the other ways in which CMR may serve as a useful tool in transplant screening, including the potential for guiding EMB. Preliminary work completed by our group shows parametric mapping patterns of T1/T2 elevations, or hotspots [53], in myocardial diseases such as rejection, reinforcing that these pathological changes are not uniform in nature. This has been suspected in the past given the false negative rate of random EMB [4,5,6]. Evaluation of these hotspots using guided EMB may provide further insight.

Conclusion

CMR native T1 parametric mapping demonstrates utility in identifying rejection and assessing for graft dysfunction in pediatric heart transplant patients, possibly beyond random endomyocardial biopsy. Further work is needed to determine how CMR can best fit into the current clinical multi-faceted approach to transplant rejection.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

BNP:

Brain natriuretic peptide

BSA:

Body surface area

CMR:

Cardiovascular magnetic resonance

DSA:

Donor specific antibodies

ECV:

Extracellular volume fraction

EMB:

Endomyocardial biopsy

GLCM:

Gray-level co-occurrence matrix

IRB:

Institutional review board

ISHLT:

International society of heart and lung transplantation

IV:

Intravenous

IVIG:

Intravenous immunoglobulin

LV:

Left ventricle/left ventricular

MOLLI:

Modified look-locker inversion

RA:

Right atrial

ROI:

Region of interest

RV:

Right ventricle/right ventricular

RVEDP:

Right ventricle end diastolic pressure

References

  1. Chin C, Naftel DC, Singh TP, Blume ED, Luikart H, Bernstein D, Gamberg P, Kirklin JK, Morrow WR. Risk factors for recurrent rejection in pediatric heart transplantation: a multicenter experience. J Heart Lung Transplant. 2004;23(2):178–85.

    Article  PubMed  Google Scholar 

  2. Dipchand AI. Current state of pediatric cardiac transplantation. Ann Cardiothorac Surg. 2018;7(1):31–55.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ameduri RK, Zheng J, Schechtman KB, Hoffman TM, Gajarski RJ, Chinnock R, Naftel DC, Kirklin JK, Dipchand AI, Canter CE. Has late rejection decreased in pediatric heart transplantation in the current era? A multi-institutional study. J Heart Lung Transplant. 2012;31(9):980–6.

    Article  PubMed  Google Scholar 

  4. Nakhleh R, Jones J, Goswitz J, Anderson E, Titus J. Correlation of endomyocardial biopsy findings with autopsy findings in human cardiac allografts. J Heart Lung Transplant. 1992;May-Jun(11):479–85.

  5. From AM, Maleszewski JJ, Rihal CS. Current status of endomyocardial biopsy. Mayo Clin Proc. 2011;86(11):1095–102.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Asimacopoulos EP, Garbern JC, Gauvreau K, Blume ED, Daly KP, Singh TP. Association of clinical rejection versus rejection on protocol biopsy with cardiac allograft vasculopathy in pediatric heart transplant recipients. Transplantation. 2020;104(1):e31–7.

    Article  PubMed  Google Scholar 

  7. Laser JA, Fowles RE, Mason JW. Endomyocardial biopsy. Cardiovasc Clin. 1985;15(1):141–63.

    CAS  PubMed  Google Scholar 

  8. Baraldi-Junkins C, Levin H, Kasper E, Rayburn B, Herskowitz A, Baughman K. Complications of endomyocardial biopsy in heart transplant patients. J Heart Lung Transplant. 1993;12(1 Pt 1):63–7.

    CAS  PubMed  Google Scholar 

  9. Daly KP, Marshall AC, Vincent JA, Zuckerman WA, Hoffman TM, Canter CE, Blume ED, Bergersen L. Endomyocardial biopsy and selective coronary angiography are low-risk procedures in pediatric heart transplant recipients: Results of a multicenter experience. J Heart Lung Transplant. 2012;31(4):398–409.

    Article  PubMed  Google Scholar 

  10. Rosenthal DN, Chin C, Nishimura K, Perry SB, Robbins RC, Reitz B, Bernstein D, Feinstein JA. Identifying cardiac transplant rejection in children: diagnostic utility of echocardiography, right heart catheterization and endomyocardial biopsy data. J Heart Lung Transplant. 2004;23(3):323–9.

    Article  PubMed  Google Scholar 

  11. Aggarwal S, Blake J, Sehgal S. Right ventricular dysfunction as an echocardiographic measure of acute rejection following heart transplantation in children. Pediatr Cardiol. 2017;38(3):442–7.

    Article  PubMed  Google Scholar 

  12. De Vlaminck I, Valantine HA, Snyder TM, Strehl C, Cohen G, Luikart H, Neff NF, Okamoto J, Bernstein D, Weisshaar D, Quake SR, Khush KK. Circulating cell-free DNA enables noninvasive diagnosis of heart transplant rejection. Sci Transl Med. 2014;6(241).

  13. Perens G, Li F, Meier S, Kaur R, Alejos JC, Fishbein M. Clinico-pathologic findings in end-stage pediatric heart transplant grafts. Pediatr Transplant. 2009;13(7):887–91.

    Article  PubMed  Google Scholar 

  14. Koestenberger M, Schweintzger S, Pocivalnik M, Kurath-Koller S, Grangl G. Relevance of right ventricular outflow tract velocity time integral (RVOT VTI) and tricuspid regurgitation velocity/RVOT VTI ratio determination in children following heart transplantation. Pediatr Cardiol. 2018;39(3):642–3.

    Article  PubMed  Google Scholar 

  15. Feingold B, Picarsic J, Lesniak A, Popp BA, Wood-Trageser MA, Demetris AJ. Late graft dysfunction after pediatric heart transplantation is associated with fibrosis and microvasculopathy by automated, digital whole-slide analysis. J Heart Lung Transplant. 2017;36(12):1336–43.

    Article  PubMed  Google Scholar 

  16. Messroghli DR, Moon JC, Ferreira VM, Grosse-Wortmann L, He T, Kellman P, Mascherbauer J, Nezafat R, Salerno M, Schelbert EB, Taylor AJ, Thompson R, Ugander M, van Heeswijk RB, Friedrich MG. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: a consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imaging (EACVI). J Cardiovasc Magn Reson. 2017;19(1):75.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Haaf P, Garg P, Messroghli DR, Broadbent DA, Greenwood JP, Plein S. Cardiac T1 mapping and extracellular volume (ECV) in clinical practice: a comprehensive review. J Cardiovasc Magn Reson. 2017;18(1):89.

    Article  Google Scholar 

  18. Vermes E, Pantaléon C, Auvet A, Cazeneuve N, Machet MC, Delhommais A, Bourguignon T, Aupart M, Brunereau L. Cardiovascular magnetic resonance in heart transplant patients: diagnostic value of quantitative tissue markers: T2 mapping and extracellular volume fraction, for acute rejection diagnosis. J Cardiovasc Magn Reson. 2018;20(1):59.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Usman AA, Taimen K, Wasielewski M, McDonald J, Shah S, Giri S, Cotts W, McGee E, Gordon R, Collins JD, Markl M, Carr JC. Cardiac magnetic resonance T2 mapping in the monitoring and follow-up of acute cardiac transplant rejection: a pilot study. Circ Cardiovasc Imaging. 2012;5(6):782–90.

    Article  PubMed  Google Scholar 

  20. Imran M, Wang L, McCrohon J, Yu C, Holloway C, Otton J, Huang J, Stehning C, Moffat KJ, Ross J, Puntmann VO, Vassiliou VS, Prasad S, Kotlyar E, Keogh A, et al. Native T1 mapping in the diagnosis of cardiac allograft rejection. JACC Cardiovasc Imaging. 2019;12(8):1618–28.

    Article  PubMed  Google Scholar 

  21. Butler CR, Savu A, Bakal JA, Toma M, Thompson R, Chow K, Wang H, Kim DH, Mengel M, Haykowsky M, Pearson GJ, Kaul P, Paterson I. Correlation of cardiovascular magnetic resonance imaging findings and endomyocardial biopsy results in patients undergoing screening for heart transplant rejection. J Heart Lung Transplant. 2015;34(5):643–50.

    Article  PubMed  Google Scholar 

  22. Rao S, Tseng SY, Pednekar A, Siddiqui S, Kocaoglu M, Fares M, Lang SM, Kutty S, Christopher AB, Olivieri LJ, Taylor MD, Alsaied T. Myocardial parametric mapping by cardiac magnetic resonance imaging in pediatric cardiology and congenital heart disease. Circ Cardiovasc Imaging. 0(0):CIRCIMAGING.120.012242.

  23. Riesenkampff E, Messroghli DR, Redington AN, Grosse-Wortmann L. Myocardial T1 mapping in pediatric and congenital heart disease. Circ Cardiovasc Imaging. 2015;8(2).

  24. Cornicelli MD, Rigsby CK, Rychlik K, Pahl E, Robinson JD. Diagnostic performance of cardiovascular magnetic resonance native T1 and T2 mapping in pediatric patients with acute myocarditis. J Cardiovasc Magn Reson. 2019;21(1):40.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Parekh K, Markl M, Deng J, de Freitas RA, Rigsby CK. T1 mapping in children and young adults with hypertrophic cardiomyopathy. Int J Cardiovasc Imaging. 2017;33(1):109–17.

    Article  PubMed  Google Scholar 

  26. Olivieri LJ, Kellman P, McCarter RJ, Cross RR, Hansen MS, Spurney CF. Native T1 values identify myocardial changes and stratify disease severity in patients with Duchenne muscular dystrophy. J Cardiovasc Magn Reson. 2017;18(1):72.

    Article  Google Scholar 

  27. Westwood M, Anderson LJ, Firmin DN, Gatehouse PD, Charrier CC, Wonke B, Pennell DJ. A single breath-hold multiecho T2* cardiovascular magnetic resonance technique for diagnosis of myocardial iron overload. J Magn Reson Imaging. 2003;18(1):33–9.

    Article  PubMed  Google Scholar 

  28. Tham EB, Haykowsky MJ, Chow K, Spavor M, Kaneko S, Khoo NS, Pagano JJ, Mackie AS, Thompson RB. Diffuse myocardial fibrosis by T1-mapping in children with subclinical anthracycline cardiotoxicity: relationship to exercise capacity, cumulative dose and remodeling. J Cardiovasc Magn Reson. 2013;15(1):48.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Messroghli DR, Plein S, Higgins DM, Walters K, Jones TR, Ridgway JP, Sivananthan MU. Human myocardium: single-Breath-hold MR T1 mapping with high spatial resolution—reproducibility study. Radiology. 2006;238(3):1004–12.

    Article  PubMed  Google Scholar 

  30. Haralick RM, Shanmugam K, Dinstein I. Textural features for image classification. IEEE Trans Syst Man Cybern. 1973;SMC3(6):610–21.

    Article  Google Scholar 

  31. Amadasun M, King R. Textural features corresponding to textural properties. IEEE Trans Syst Man Cybern. 1989;19(5):1264–74.

    Article  Google Scholar 

  32. Doshi A, Olivieri L, Trayanova N. Abstract 12520: texture analysis of cardiac magnetic resonance native T1 maps characterizes the proarrhythmic substrate in pediatric myocarditis. Circulation. 2019;140(Suppl_1):A12520–A12520.

    Google Scholar 

  33. Baessler B, Luecke C, Lurz J, Klingel K, Das A, von Roeder M, de Waha-Thiele S, Besler C, Rommel K-P, Maintz D, Gutberlet M, Thiele H, Lurz P. Cardiac MRI and texture analysis of myocardial T1 and T2 maps in myocarditis with acute versus chronic symptoms of heart failure. Radiology. 2019;292(3):608–17.

    Article  PubMed  Google Scholar 

  34. El-Rewaidy H, Neisius U, Nakamori S, Ngo L, Rodriguez J, Manning WJ, Nezafat R. Characterization of interstitial diffuse fibrosis patterns using texture analysis of myocardial native T1 mapping. Lionetti V, ed. PLOS ONE. 2020;15(6):e0233694.

  35. Neisius U, El-Rewaidy H, Nakamori S, Rodriguez J, Manning WJ, Nezafat R. Radiomic analysis of myocardial native T1 imaging discriminates between hypertensive heart disease and hypertrophic cardiomyopathy. JACC Cardiovasc Imaging. 2019;12(10):1946–54.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Stewart S, Winters GL, Fishbein MC, Tazelaar HD, Kobashigawa J, Abrams J, Andersen CB, Angelini A, Berry GJ, Burke MM, Demetris AJ, Hammond E, Itescu S, Marboe CC, McManus B, et al. Revision of the 1990 working formulation for the standardization of nomenclature in the diagnosis of heart rejection. J Heart Lung Transplant. 2005;24(11):1710–20.

    Article  PubMed  Google Scholar 

  37. American Heart Association Writing Group on Myocardial Segmentation and Registration for Cardiac Imaging, Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a Statement for Healthcare Professionals From the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002;105(4):539–42.

    Article  Google Scholar 

  38. Sethi N, Doshi A, Doshi T, Cross R, Cronin I, Amin E, Kanter J, Scheel J, Khan S, Campbell-Washburn A, Olivieri L. Quantitative cardiac magnetic resonance T2 imaging offers ability to non-invasively predict acute allograft rejection in children. Cardiol Young. 2020;30(6):852–9.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Vallières M, Freeman CR, Skamene SR, El Naqa I. A radiomics model from joint FDG-PET and MRI texture features for the prediction of lung metastases in soft-tissue sarcomas of the extremities. Phys Med Biol. 2015;60(14):5471–96.

    Article  PubMed  Google Scholar 

  40. Assefa D, Keller H, Ménard C, Laperriere N, Ferrari RJ, Yeung I. Robust texture features for response monitoring of glioblastoma multiforme on T1-weighted and T2-FLAIR MR images: a preliminary investigation in terms of identification and segmentation: GBM texture features on T1- and T2-weighted MR images. Med Phys. 2010;37(4):1722–36.

    Article  PubMed  Google Scholar 

  41. Thibault G, Fertil B, Navarro C, Pereira S, Lévy N, Sequeira J, Mari J-L. Texture indexes and gray level size zone matrix application to cell nuclei classification. In: 10th International Conference on Pattern Recognition and Information Processing. 2009.

  42. Aerts HJWL, Velazquez ER, Leijenaar RTH, Parmar C, Grossmann P, Carvalho S, Bussink J, Monshouwer R, Haibe-Kains B, Rietveld D, Hoebers F, Rietbergen MM, Leemans CR, Dekker A, Quackenbush J, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5(1):4006.

    Article  CAS  PubMed  Google Scholar 

  43. Dolan RS, Rahsepar AA, Blaisdell J, Suwa K, Ghafourian K, Wilcox JE, Khan SS, Vorovich EE, Rich JD, Anderson AS, Yancy CW, Collins JD, Carr JC, Markl M. Multiparametric cardiac magnetic resonance imaging can detect acute cardiac allograft rejection after heart transplantation. JACC Cardiovasc Imaging. 2019;12(8):1632–41.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Husain N, Watanabe K, Berhane H, Gupta A, Markl M, Rigsby CK, Robinson JD. Multi-parametric cardiovascular magnetic resonance with regadenoson stress perfusion is safe following pediatric heart transplantation and identifies history of rejection and cardiac allograft vasculopathy. J Cardiovasc Magn Reson. 2021;23(1):135.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Greenway SC, Dallaire F, Kantor PF, Dipchand AI, Chaturvedi RR, Warade M, Riesenkampff E, Yoo S-J, Grosse-Wortmann L. Magnetic resonance imaging of the transplanted pediatric heart as a potential predictor of rejection. World J Transplant. 2016;6(4):751.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Peng DM, Ding VY, Hollander SA, Khalapyan T, Dykes JC, Rosenthal DN, Almond CS, Sakarovitch C, Desai M, McElhinney DB. Long-term surveillance biopsy: Is it necessary after pediatric heart transplant? Pediatr Transplant. 2019;23(1): e13330.

    Article  PubMed  Google Scholar 

  47. Evers PD, Jorgensen N, Hong B, Albers E, Kemna M, Friedland-Little J, Boucek RJ, Law Y. Outcomes of low-intensity biopsy surveillance for rejection in paediatric cardiac transplantation. Cardiol Young. 2019;29(7):910–6.

    Article  PubMed  Google Scholar 

  48. Chi N-H, Chou N-K, Tsao C-I, Huang S-C, Wu I-H, Yu H-Y, Chen Y-S, Wang S-S. Endomyocardial biopsy in heart transplantation: schedule or event? Transpl Proc. 2012;44(4):894–6.

    Article  Google Scholar 

  49. Shah KB, Flattery MP, Smallfield MC, Merinar G, Tang DG, Sheldon EH, Thacker LR, Kasirajan V, Cooke RH, Hess ML. Surveillance endomyocardial biopsy in the modern era produces low diagnostic yield for cardiac allograft rejection. Transplantation. 2015;99(8):e75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Moreo A, Ambrosio G, De Chiara B, Pu M, Tran T, Mauri F, Raman SV. Influence of myocardial fibrosis on left ventricular diastolic function: noninvasive assessment by cardiac magnetic resonance and echo. Circ Cardiovasc Imaging. 2009;2(6):437–43.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Kobashigawa JA, Itagaki BK, Razi RR, Patel JK, Chai W, Kawano MA, Goldstein Z, Kittleson MM, Fishbein MC. Correlation between myocardial fibrosis and restrictive cardiac physiology in patients undergoing retransplantation. Clin Transplant. 2013;27(6):E679–84.

    Article  PubMed  Google Scholar 

  52. DiPietro A, Sexson-Tejtel SK, Pignatelli R, Cui Q, Cabrera AG, Price J, Mitchell D, Dreyer WJ, Jeewa A, Denfield S. Noninvasive predictors of diastolic dysfunction in children after cardiac transplant. J Heart Lung Transplant. 2013;32(4):S190.

    Article  Google Scholar 

  53. Patel D, Contento J, Cleveland V, Hamman K, Irizarry S, Downing T, Kanter J, Berger III J, Cross R, Christopher A, Chow K. Abstract: Accuracy of free-breathing SASHA in identifying T1 and T2 elevations in pediatric orthotopic heart transplant patients. SCMR Scientific Sessions 2022. 2022.

Download references

Acknowledgements

Children’s National Interventional Cardiac Magnetic Resonance team, Children’s National Stimulating Access to Research in Residency (CNStARR).

Funding

This publication was supported by Award Number R01-HL144494 and R38HL143585 from the NIH National Heart Lung and Blood Institute. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the National Heart Lung and Blood Institute or the National Institutes of Health.

Author information

Authors and Affiliations

Authors

Contributions

DP performed data collection, analysis, and write-up of this manuscript. NG performed digital analysis of histological slides for fibrosis percentage and assisted with data collection. JM performed data collection during cases. AD performed texture analysis. KH performed patient recruitment, consent, and data input during cases. CR performed histological examination of endomyocardial biopsies, providing transplant rejection grading, and was involved in the creation of the fibrosis percentage measurement tool. AR was involved in the creation of the fibrosis percentage measurement tool and served as consultant regarding assessing fibrosis quantitatively using digital histological analysis. RC performed cardiac magnetic resonance and parametric mapping collection during cases and was involved in design of the study. JK performed endomyocardial biopsy attainment and catheterization hemodynamic data collection and was involved in design of the study. JB served as expert consultant regarding pediatric heart transplant patients and often was involved in deciding treatment plans for particular patients; was also involved in design of the study. LO is the principal investigator of this study, involved in design and execution of the study including performing cardiac magnetic resonance and parametric mapping collection during cases. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Devika P. Richmann.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Children’s National Institutional Review Board. Informed consent/assent was obtained prior to each case.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Definition and calculations of gray-level co-occurrence matrix texture features [30, 4042].

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Richmann, D.P., Gurijala, N., Mandell, J.G. et al. Native T1 mapping detects both acute clinical rejection and graft dysfunction in pediatric heart transplant patients. J Cardiovasc Magn Reson 24, 51 (2022). https://doi.org/10.1186/s12968-022-00875-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12968-022-00875-z

Keywords