Monitoring pharmacologically induced immunosuppression by immune repertoire sequencing to detect acute allograft rejection in heart transplant patients: a proof-of-concept diagnostic accuracy study

Christopher Vollmers, Iwijn De Vlaminck, Hannah A Valantine, Lolita Penland, Helen Luikart, Calvin Strehl, Garrett Cohen, Kiran K Khush, Stephen R Quake, Christopher Vollmers, Iwijn De Vlaminck, Hannah A Valantine, Lolita Penland, Helen Luikart, Calvin Strehl, Garrett Cohen, Kiran K Khush, Stephen R Quake

Abstract

Background: It remains difficult to predict and to measure the efficacy of pharmacological immunosuppression. We hypothesized that measuring the B-cell repertoire would enable assessment of the overall level of immunosuppression after heart transplantation.

Methods and findings: In this proof-of-concept study, we implemented a molecular-barcode-based immune repertoire sequencing assay that sensitively and accurately measures the isotype and clonal composition of the circulating B cell repertoire. We used this assay to measure the temporal response of the B cell repertoire to immunosuppression after heart transplantation. We selected a subset of 12 participants from a larger prospective cohort study (ClinicalTrials.gov NCT01985412) that is ongoing at Stanford Medical Center and for which enrollment started in March 2010. This subset of 12 participants was selected to represent post-heart-transplant events, with and without acute rejection (six participants with moderate-to-severe rejection and six without). We analyzed 130 samples from these patients, with an average follow-up period of 15 mo. Immune repertoire sequencing enables the measurement of a patient's net state of immunosuppression (correlation with tacrolimus level, r = -0.867, 95% CI -0.968 to -0.523, p = 0.0014), as well as the diagnosis of acute allograft rejection, which is preceded by increased immune activity with a sensitivity of 71.4% (95% CI 30.3% to 94.9%) and a specificity of 82.0% (95% CI 72.1% to 89.1%) (cell-free donor-derived DNA as noninvasive gold standard). To illustrate the potential of immune repertoire sequencing to monitor atypical post-transplant trajectories, we analyzed two more patients, one with chronic infections and one with amyloidosis. A larger, prospective study will be needed to validate the power of immune repertoire sequencing to predict rejection events, as this proof-of-concept study is limited to a small number of patients who were selected based on several criteria including the availability of a large number of samples and the absence or presence of rejection events.

Conclusions: If confirmed in larger, prospective studies, the method described here has potential applications in the tailored management of post-transplant immunosuppression and, more broadly, as a method for assessing the overall activity of the immune system.

Conflict of interest statement

CV, IDV, and SRQ have a patent S14-048 (STAN-1110PRV) 61/951,908: Provisional Application pending. HAV, HL, CS, GC, and KKK have nothing to declare.

Figures

Fig 1. Monitoring overall immunosuppression.
Fig 1. Monitoring overall immunosuppression.
Organ transplant recipients receive immunosuppressive therapy, which usually includes calcineurin inhibitors. These drugs inhibit T cell activation, leading to fewer B cells expressing class-switched and mutated IGH transcripts. We measured overall immunosuppression using immune repertoire sequencing (IGH-Seq). A schematic of the IGH-Seq protocol and data analysis is shown on the right. C_Region, primer pool specific to constant region; CDR3, complementary determining region 3; PBMC, peripheral blood mononuclear cell; UID, unique identifier; V_FR3, primer pool specific to V segments.
Fig 2. Using immune repertoire sequencing to…
Fig 2. Using immune repertoire sequencing to measure immunosuppression.
(A) Relative IGH sequence composition at increasing expression levels (molecules/sequence) is shown across all individuals (n = 12) in a histogram for days 1 and 180 post-transplant. Bar colors indicate the contribution of different immunoglobulin isotypes. Inserts show rescaled histograms. Highly expressed (HE) sequences (sequences represented by more than one molecule) are labeled in the insert. (B) ABS level (defined as the ratio of all highly expressed IgA, IgD, IgE, IgG, and mutated IgM sequences to the total number of molecules) as a function of tacrolimus trough level (n = 130) (average ± standard error of the mean). (C) Median ABS level (colors indicate isotype contribution), tacrolimus trough level, and total anellovirus load for all study patients (n = 12) during the first year after transplant.
Fig 3. Reduced level of overall immunosuppression…
Fig 3. Reduced level of overall immunosuppression correlates with acute rejection events.
ABS levels are shown for two individuals with (C and D) and two individuals without (A and B) moderate-to-severe acute rejection events (top). Colors indicate isotype contribution. cfdDNA levels are shown for the same individuals (bottom). Biopsy grades are shown as vertical blue lines (bottom). Thin lines indicate mild ACR events (grade 1R), and thick lines indicate moderate-to-severe ACR events (grade ≥ 2R).
Fig 4. ABS levels correlate well with…
Fig 4. ABS levels correlate well with cfdDNA and endomyocardial biopsy results.
(A) ABS levels plotted for all study participants. Acute rejection events as diagnosed by endomyocardial biopsy (ACR, violet; acute antibody-mediated rejection [AMR], green) or cfdDNA level (pink) are shown in color. (B) Samples are classified by rejection status diagnosed by biopsy as no evidence of rejection, mild ACR (biopsy grade 1R), or moderate-to-severe ACR (biopsy grade ≥ 2R), and by rejection status diagnosed by the cfdDNA assay. Samples are colored as in (A) if at any point the organ recipient experienced a serious rejection event. As biopsy grades are independent of cfdDNA measurements, there is overlap between samples in the elevated cfdDNA and biopsy-proven rejection bins. The groups were compared using a one-sided Mann-Whitney U test, which has high efficiency for both normally and non-normally distributed datasets. (C) Receiver operating characteristic (ROC) curves for ABS levels compared to cfdDNA (pink) and endomyocardial biopsy (violet) results. FPR, false rejection rate; TPR, true positive rate. (D) ABS levels measured before and after rejection diagnosis. The black line shows the single-exponent fit of ABS levels measured prior to rejection, y = a × (exp[b × t]), with best fit values (least squares) of a = 0.0317 and b = 0.00468. The gray line shows the single-exponent fit of ABS levels after diagnosis and during treatment of rejection: y = a × (exp[b× t]), with best fit values (least squares) of a = 0.0344 and b = −0.0324. Green lines show the median ABS levels in the covered time periods.
Fig 5. Reduced net immunosuppression correlates with…
Fig 5. Reduced net immunosuppression correlates with acute rejection.
ABS level is shown for three individuals: (A) a heart transplant recipient without moderate-to-severe rejection events, (B) a heart transplant recipient with infectious complications, and (C) a heart transplant recipient with AL amyloidosis. For the patient with AL amyloidosis (C), the approximate schedule of drug treatment is indicated at the top of the panel. Colors indicate isotype contribution. cfdDNA level is shown for the same individuals (bottom). Biopsy grades are shown as vertical blue lines (bottom): thin lines indicate mild rejection events (grade 1R), and thick lines indicate moderate-to-severe rejection events (grade ≥ 2R). SCT, stem cell transplant.

References

    1. Plosker GL, Foster RH. Tacrolimus: a further update of its pharmacology and therapeutic use in the management of organ transplantation. Drugs. 2000;59:323–389.
    1. Liu J, Farmer JD, Lane WS, Friedman J, Weissman I, Schreiber SL. Calcineurin is a common target of cyclophilin-cyclosporin A and FKBP-FK506 complexes. Cell. 1991;66:807–815.
    1. Heidt S, Roelen DL, Eijsink C, Eikmans M, Van Kooten C, Claas FHJ, et al. Calcineurin inhibitors affect B cell antibody responses indirectly by interfering with T cell help. Clin Exp Immunol. 2010;159:199–207. 10.1111/j.1365-2249.2009.04051.x
    1. Perez-Andres M, Paiva B, Nieto WG, Caraux A, Schmitz A, Almeida J, et al. Human peripheral blood B-cell compartments: a crossroad in B-cell traffic. Cytometry B Clin Cytom. 2010;78 (Suppl 1):S47–S60. 10.1002/cyto.b.20547
    1. Newell KA, Asare A, Kirk AD, Gisler TD, Bourcier K, Suthanthiran M, et al. Identification of a B cell signature associated with renal transplant tolerance in humans. J Clin Invest. 2010;120:1836–1847. 10.1172/JCI39933
    1. Glanville J, Kuo TC, von Budingen H-C, Guey L, Berka J, Sundar PD, et al. Naive antibody gene-segment frequencies are heritable and unaltered by chronic lymphocyte ablation. Proc Natl Acad Sci U S A. 2011;108:20066–20071. 10.1073/pnas.1107498108
    1. Vollmers C, Sit RV, Weinstein JA, Dekker CL, Quake SR. Genetic measurement of memory B-cell recall using antibody repertoire sequencing. Proc Natl Acad Sci U S A. 2013;110:13463–13468. 10.1073/pnas.1312146110
    1. Boyd SD, Marshall EL, Merker JD, Maniar JM, Zhang LN, Sahaf B, et al. Measurement and clinical monitoring of human lymphocyte clonality by massively parallel VDJ pyrosequencing. Sci Transl Med. 2009;1:12ra23
    1. Georgiou G, Ippolito GC, Beausang J, Busse CE, Wardemann H, Quake SR. The promise and challenge of high-throughput sequencing of the antibody repertoire. Nat Biotechnol. 2014;32:158–168. 10.1038/nbt.2782
    1. Jiang N, He J, Weinstein JA, Penland L, Sasaki S, He X-S, et al. Lineage structure of the human antibody repertoire in response to influenza vaccination. Sci Transl Med. 2013;5:171ra19 10.1126/scitranslmed.3004794
    1. Logan AC, Gao H, Wang C, Sahaf B, Jones CD, Marshall EL, et al. High-throughput VDJ sequencing for quantification of minimal residual disease in chronic lymphocytic leukemia and immune reconstitution assessment. Proc Natl Acad Sci U S A. 2011:108:21194–21199. 10.1073/pnas.1118357109
    1. De Vlaminck I, Khush KK, Strehl C, Kohli B, Luikart H, Neff NF, et al. Temporal response of the human virome to immunosuppression and antiviral therapy. Cell. 2013;155:1178–1187. 10.1016/j.cell.2013.10.034
    1. De Vlaminck I, Valantine HA, Snyder TM, Strehl C, Cohen G, Luikart H, et al. Circulating cell-free DNA enables noninvasive diagnosis of heart transplant rejection. Sci Transl Med. 2014;6:241ra77 10.1126/scitranslmed.3007803
    1. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–410.
    1. Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9:90–95.
    1. SciPy Developers. SciPy [computer program]. . Accessed 10 September 2015.
    1. Oliphant TE. Python for scientific computing. Comput Sci Eng. 2007;9:10–20.
    1. Van der Walt S, Colbert SC, Varoquaux G. The NumPy array: a structure for efficient numerical computation. Comput Sci Eng. 2011;13:22–30.
    1. Xia LC, Cram JA, Chen T, Fuhrman JA, Sun F. Accurate genome relative abundance estimation based on shotgun metagenomic reads. PLoS ONE. 2011;6:e27992 10.1371/journal.pone.0027992
    1. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12:2825–2830.
    1. Wilson EB. Probable inference, the law of succession, and statistical inference. J Am Stat Assoc. 1927;22:209.
    1. Julious SA. Two-sided confidence intervals for the single proportion: comparison of seven methods by Robert G. Newcombe, Statistics in Medicine 1998; 17:857–872. Stat Med. 2005;24:3383–3384.
    1. Nair N, Newell EW, Vollmers C, Quake SR, Morton JM, Davis MM, et al. High-dimensional immune profiling of total and rotavirus VP6-specific intestinal and circulating B cells by mass cytometry. Mucosal Immunol. 2015. April 22
    1. Tayebi H, Lapierre V, Saas P, Lienard A, Sutton L, Milpied N, et al. Enhanced activation of B cells in a granulocyte colony-stimulating factor-mobilized peripheral blood stem cell graft. Br J Haematol. 2001;114:698–700.
    1. Jensen GS, Poppema S, Mant MJ, Pilarski LM. Transition in CD45 isoform expression during differentiation of normal and abnormal B cells. Int Immunol. 1989;1:229–236.
    1. Hoffman FM. Outcomes and complications after heart transplantation: a review. J Cardiovasc Nurs. 2005;20:S31–S42.
    1. Giese T, Zeier M, Meuer S. Analysis of NFAT-regulated gene expression in vivo: a novel perspective for optimal individualized doses of calcineurin inhibitors. Nephrol Dial Transplant. 2004;19 (Suppl 4):iv55–iv60.
    1. Wieland E, Olbricht CJ, Süsal C, Gurragchaa P, Böhler T, Israeli M, et al. Biomarkers as a tool for management of immunosuppression in transplant patients. Ther Drug Monit. 2010;32:560–572. 10.1097/FTD.0b013e3181efb3d2
    1. Rossano JW, Denfield SW, Kim JJ, Price JF, Jefferies JL, Decker JA, et al. Assessment of the Cylex ImmuKnow cell function assay in pediatric heart transplant patients. J Hear Lung Transplant. 2009;28:26–31.
    1. Pham MX, Teuteberg JJ, Kfoury AG, Starling RC, Deng MC, Cappola TP, et al. Gene-expression profiling for rejection surveillance after cardiac transplantation. N Engl J Med. 2010;362:1890–1900. 10.1056/NEJMoa0912965
    1. Deng MC, Eisen HJ, Mehra MR, Billingham M, Marboe CC, Berry G, et al. Noninvasive discrimination of rejection in cardiac allograft recipients using gene expression profiling. Am J Transplant. 2006;6:150–160.

Source: PubMed

3
S'abonner