Relationship between donor fraction cell-free DNA and clinical rejection in heart transplantation.

BACKGROUND: Clinical rejection (CR) defined as decision to treat clinically suspected rejection with change in immunotherapy based on clinical presentation with or without diagnostic biopsy findings is an important part of care in heart transplantation. We sought to assess the utility of donor fraction cell-free DNA (DF cfDNA) in CR and the utility of serial DF cfDNA in CR patients in predicting outcomes of clinical interest. METHODS: Patients with heart transplantation were enrolled in two sequential, multi-center, prospective observational studies. Blood samples were collected for surveillance or clinical events. Clinicians were blinded to the results of DF cfDNA. RESULTS: A total of 835 samples from 269 subjects (57% pediatric) were included for this analysis, including 28 samples associated with CR were analyzed. Median DF cfDNA was 0.43 (IQR 0.15, 1.36)% for CR and 0.10 (IQR 0.07, 0.16)% for healthy controls (p < .0001). At cutoff value of 0.13%, the area under curve (AUC) was 0.82, sensitivity of 0.86, specificity of 0.67, and negative predictive value of 0.99. There was serial decline in DF cfDNA post-therapy, however, those with cardiovascular events (cardiac arrest, need for mechanical support or death) showed significantly higher levels of DF cfDNA on Day 0 (2.11 vs 0.31%) and Day 14 (0.51 vs 0.22%) compared to those who did not have such an event (p < .0001). CONCLUSION: DF cfDNA has excellent agreement with clinical rejection and, importantly, serial measurement of DF cfDNA predict clinically significant outcomes post treatment for rejection in these patients.

Early changes in cell-free DNA levels in newly transplanted heart transplant patients.

Heart transplantation is a well-established therapy for end-stage heart failure in children and young adults. The highest risk of graft loss occurs in the first 60 days post-transplant. Donor fraction of cell-free DNA is a highly sensitive marker of graft injury. Changes in cell-free DNA levels have not previously been studied in depth in patients early after heart transplant. A prospective study was conducted among heart transplant recipients at a single pediatric heart center. Blood samples were collected from children and young adult transplant patients at three time points within 10 days of transplantation. DF and total cell-free DNA levels were measured using a targeted method (myTAIHEART ). In 17 patients with serial post-transplant samples, DF peaks in the first 2 days after transplant (3.5%, [1.9-10]%) and then declines toward baseline (0.27%, [0.19-0.52]%) by 6-9 days. There were 4 deaths in the first year among the 10 patients with complete sample sets, and 3 out of 4 who died had a late rise or blunted decline in donor fraction. Patients who died trended toward an elevated total cell-free DNA at 1 week (41.5, [34-65] vs 13.6, [6.2-22] P = .07). Donor fraction peaks early after heart transplant and then declines toward baseline. Patients without sustained decline in donor fraction and/or elevated total cell-free DNA at 1 week may have worse outcomes.

Impact of MYH6 variants in hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a clinically and anatomically severe form of congenital heart disease (CHD). Although prior studies suggest that HLHS has a complex genetic inheritance, its etiology remains largely unknown. The goal of this study was to characterize a risk gene in HLHS and its effect on HLHS etiology and outcome. We performed next-generation sequencing on a multigenerational family with a high prevalence of CHD/HLHS, identifying a rare variant in the α-myosin heavy chain (MYH6) gene. A case-control study of 190 unrelated HLHS subjects was then performed and compared with the 1000 Genomes Project. Damaging MYH6 variants, including novel, missense, in-frame deletion, premature stop, de novo, and compound heterozygous variants, were significantly enriched in HLHS cases (P < 1 × 10-5). Clinical outcomes analysis showed reduced transplant-free survival in HLHS subjects with damaging MYH6 variants (P < 1 × 10-2). Transcriptome and protein expression analyses with cardiac tissue revealed differential expression of cardiac contractility genes, notably upregulation of the ß-myosin heavy chain (MYH7) gene in subjects with MYH6 variants (P < 1 × 10-3). We subsequently used patient-specific induced pluripotent stem cells (iPSCs) to model HLHS in vitro. Early stages of in vitro cardiomyogenesis in iPSCs derived from two unrelated HLHS families mimicked the increased expression of MYH7 observed in vivo (P < 1 × 10-2), while revealing defective cardiomyogenic differentiation. Rare, damaging variants in MYH6 are enriched in HLHS, affect molecular expression of contractility genes, and are predictive of poor outcome. These findings indicate that the etiology of MYH6-associated HLHS can be informed using iPSCs and suggest utility in future clinical applications.

Nuclear Cell-Free DNA Predicts Adverse Events After Pediatric Cardiothoracic Surgery.

BACKGROUND: Preoperative risk stratification in cardiac surgery includes patient and procedure factors that are used in clinical decision-making. Despite these tools, unidentified factors contribute to variation in outcomes. Identification of latent physiologic risk factors may strengthen predictive models. Nuclear cell-free DNA (ncfDNA) increases with tissue injury and drops to baseline levels rapidly. The goal of this investigation is to measure and to observe ncfDNA kinetics in children undergoing heart operations with cardiopulmonary bypass (CPB), linking biomarkers, organ dysfunction, and outcomes. METHODS: This is a prospective observational study of 116 children <18 years and >3 kg undergoing operations with CPB. Plasma ncfDNA samples were collected and processed in a stepwise manner at predefined perioperative time points. The primary outcome measure was occurrence of postoperative cardiac arrest or extracorporeal membrane oxygenation. RESULTS: Data were available in 116 patients (median age, 0.9 years [range, 0-17.4 years]; median weight, 7.8 kg [range, 3.2-98 kg]). The primary outcome was met in 6 of 116 (5.2%). Risk of primary outcome was 2% with ncfDNA <20 ng/mL and 33% with ncfDNA >20 ng/mL (odds ratio, 25; CI, 3.96-158; P = .001). Elevated ncfDNA was associated with fewer hospital-free days (P < .01). CONCLUSIONS: This study analyzes ncfDNA kinetics in children undergoing operations with CPB for congenital heart disease. Elevated preoperative ncfDNA is strongly associated with postoperative arrest and extracorporeal membrane oxygenation. Further studies are needed to validate this technology as a tool to predict morbidity in children after cardiac surgical procedures.

Human gene copy number spectra analysis in congenital heart malformations.

The clinical significance of copy number variants (CNVs) in congenital heart disease (CHD) continues to be a challenge. Although CNVs including genes can confer disease risk, relationships between gene dosage and phenotype are still being defined. Our goal was to perform a quantitative analysis of CNVs involving 100 well-defined CHD risk genes identified through previously published human association studies in subjects with anatomically defined cardiac malformations. A novel analytical approach permitting CNV gene frequency "spectra" to be computed over prespecified regions to determine phenotype-gene dosage relationships was employed. CNVs in subjects with CHD (n = 945), subphenotyped into 40 groups and verified in accordance with the European Paediatric Cardiac Code, were compared with two control groups, a disease-free cohort (n = 2,026) and a population with coronary artery disease (n = 880). Gains (≥200 kb) and losses (≥100 kb) were determined over 100 CHD risk genes and compared using a Barnard exact test. Six subphenotypes showed significant enrichment (P ≤ 0.05), including aortic stenosis (valvar), atrioventricular canal (partial), atrioventricular septal defect with tetralogy of Fallot, subaortic stenosis, tetralogy of Fallot, and truncus arteriosus. Furthermore, CNV gene frequency spectra were enriched (P ≤ 0.05) for losses at: FKBP6, ELN, GTF2IRD1, GATA4, CRKL, TBX1, ATRX, GPC3, BCOR, ZIC3, FLNA and MID1; and gains at: PRKAB2, FMO5, CHD1L, BCL9, ACP6, GJA5, HRAS, GATA6 and RUNX1. Of CHD subjects, 14% had causal chromosomal abnormalities, and 4.3% had likely causal (significantly enriched), large, rare CNVs. CNV frequency spectra combined with precision phenotyping may lead to increased molecular understanding of etiologic pathways.

Elevated nuclear and mitochondrial cell-free deoxyribonucleic acid measurements are associated with death after infant cardiac surgery.

OBJECTIVES: Mortality rates following pediatric cardiac surgery with cardiopulmonary bypass have declined over decades, but have plateaued in recent years. This is in part attributable to persistent issues with postoperative global inflammation and myocardial dysfunction, commonly manifested by systemic inflammatory response syndrome and low cardiac output syndrome, respectively. Quantified cell-free DNA (cfDNA), of nuclear or mitochondrial origin, has emerged as a biomarker for both inflammation and myocardial injury. Recent data suggest that nuclear cfDNA (ncfDNA) may quantify inflammation, whereas mitochondrial cfDNA (mcfDNA) may correlate with the degree of myocardial injury. We hypothesize that threshold levels of ncfDNA and mcfDNA can be established that are sensitive and specific for postoperative mortality mediated through independent pathways, and that association will be enhanced with combined analysis. METHODS: Prospective observational study of infants younger than age 1 year undergoing planned surgery with cardiopulmonary bypass. The study received institutional review board approval. Samples were drawn before skin incision, immediately after completion of cardiopulmonary bypass, and subsequently at predetermined intervals postoperatively. Association of early postoperative ncfDNA and mcfDNA levels with mortality were assessed by logistic regression with cut-points chosen by receiving operating characteristic curve exploration. RESULTS: Data were available in 59 patients. Median age and weight were 122 days (interquartile range, 63-154 days) and 4.9 kg (interquartile range, 3.9-6.2 kg). Median STAT category was 3 (interquartile range, 1-4). The primary outcome of death was met in 3 out of 59 (5%). Combined analysis of ncfDNA and mcfDNA levels at 12 hours after the initiation of cardiopulmonary bypass with death at a threshold of 50 ng/mL ncfDNA and 17 copies/µL mcfDNA yielded 100% sensitivity and negative predictive value. The specificity (91%) and positive predictive value (38%) increased through combined analysis compared with univariate analysis. Combined analysis exhibited high specificity (93%) and negative predictive value (78%) for prolonged (>30 postoperative days) hospitalization. CONCLUSIONS: Combined analysis of early postoperative ncfDNA and mcfDNA can stratify risk of mortality and prolonged hospitalization following infant cardiac surgery. Evaluation of both ncfDNA and mcfDNA to identify states of generalized inflammation and myocardial injury may allow for targeted interventions and improved outcomes.

Total Cell-Free DNA Predicts Death and Infection Following Pediatric and Adult Heart Transplantation.

BACKGROUND: Elevated total cell-free DNA (TCF) concentration has been associated with critical illness in adults and elevated donor fraction (DF), the ratio of donor specific cell-free DNA to total cell-free DNA present in the recipient's plasma, is associated with rejection following cardiac transplantation. This study investigates relationships between TCF and clinical outcomes after heart transplantation. METHODS: A prospective, blinded, observational study of 87 heart transplantation recipients was performed. Samples were collected at transplantation, prior to endomyocardial biopsy, during treatment for rejection, and at hospital readmissions. Longitudinal clinical data were collected and entered into a RedCAP (Vanderbilt University) database. TCF and DF levels were correlated with endomyocardial biopsy and angiography results, as well as clinical outcomes. Logistic regression for modeling and repeated measures analysis using generalized linear modeling was used. The standard receiver operating characteristic curve, hazard ratios, and odds ratios were calculated. RESULTS: There were 257 samples from 87 recipients analyzed. TCF greater than 50 ng/mL were associated with increased mortality (P = .01, area under the curve 0.93, sensitivity 0.44, specificity 0.97) and treatment for infection (P < .005, area under the curve 0.68, sensitivity 0.45, specificity 0.96). Increased DF was not correlated with treatment for infection. DF was associated with rejection and cardiac allograft vasculopathy (P < .001), but TCF was not. CONCLUSIONS: TCF elevation predicted death and treatment for infection. DF elevation predicted histopathologic acute rejection and cardiac allograft vasculopathy. Surveillance of TCF and DF levels may inform treatment after heart transplantation.

Cell-free DNA donor fraction analysis in pediatric and adult heart transplant patients by multiplexed allele-specific quantitative PCR: Validation of a rapid and highly sensitive clinical test for stratification of rejection probability.

Lifelong noninvasive rejection monitoring in heart transplant patients is a critical clinical need historically poorly met in adults and unavailable for children and infants. Cell-free DNA (cfDNA) donor-specific fraction (DF), a direct marker of selective donor organ injury, is a promising analytical target. Methodological differences in sample processing and DF determination profoundly affect quality and sensitivity of cfDNA analyses, requiring specialized optimization for low cfDNA levels typical of transplant patients. Using next-generation sequencing, we previously correlated elevated DF with acute cellular and antibody-mediated rejection (ACR and AMR) in pediatric and adult heart transplant patients. However, next-generation sequencing is limited by cost, TAT, and sensitivity, leading us to clinically validate a rapid, highly sensitive, quantitative genotyping test, myTAIHEART®, addressing these limitations. To assure pre-analytical quality and consider interrelated cfDNA measures, plasma preparation was optimized and total cfDNA (TCF) concentration, DNA fragmentation, and DF quantification were validated in parallel for integration into myTAIHEART reporting. Analytical validations employed individual and reconstructed mixtures of human blood-derived genomic DNA (gDNA), cfDNA, and gDNA sheared to apoptotic length. Precision, linearity, and limits of blank/detection/quantification were established for TCF concentration, DNA fragmentation ratio, and DF determinations. For DF, multiplexed high-fidelity amplification followed by quantitative genotyping of 94 SNP targets was applied to 1168 samples to evaluate donor options in staged simulations, demonstrating DF call equivalency with/without donor genotype. Clinical validation studies using 158 matched endomyocardial biopsy-plasma pairs from 76 pediatric and adult heart transplant recipients selected a DF cutoff (0.32%) producing 100% NPV for ≥2R ACR. This supports the assay's conservative intended use of stratifying low versus increased probability of ≥2R ACR. myTAIHEART is clinically validated for heart transplant recipients ≥2 months old and ≥8 days post-transplant, expanding opportunity for noninvasive transplant rejection assessment to infants and children and to all recipients >1 week post-transplant.

Novel KLHL26 variant associated with a familial case of Ebstein's anomaly and left ventricular noncompaction.

BACKGROUND: Ebstein's anomaly (EA) is a rare congenital heart disease of the tricuspid valve and right ventricle. Patients with EA often manifest with left ventricular noncompaction (LVNC), a cardiomyopathy. Despite implication of cardiac sarcomere genes in some cases, very little is understood regarding the genetic etiology of EA/LVNC. Our study describes a multigenerational family with at least 10 of 17 members affected by EA/LVNC. METHODS: We performed echocardiography on all family members and conducted exome sequencing of six individuals. After identifying candidate variants using two different bioinformatic strategies, we confirmed segregation with phenotype using Sanger sequencing. We investigated structural implications of candidate variants using protein prediction models. RESULTS: Exome sequencing analysis of four affected and two unaffected members identified a novel, rare, and damaging coding variant in the Kelch-like family member 26 (KLHL26) gene located on chromosome 19 at position 237 of the protein (GRCh37). This variant region was confirmed by Sanger sequencing in the remaining family members. KLHL26 (c.709C > T p.R237C) segregates only with EA/LVNC-affected individuals (FBAT p < .05). Investigating structural implications of the candidate variant using protein prediction models suggested that the KLHL26 variant disrupts electrostatic interactions when binding to part of the ubiquitin proteasome, specifically Cullin3 (CUL3), a component of E3 ubiquitin ligase. CONCLUSION: In this familial case of EA/LVNC, we have identified a candidate gene variant, KLHL26 (p.R237C), which may have an important role in ubiquitin-mediated protein degradation during cardiac development.

The inferred cardiogenic gene regulatory network in the mammalian heart.

Cardiac development is a complex, multiscale process encompassing cell fate adoption, differentiation and morphogenesis. To elucidate pathways underlying this process, a recently developed algorithm to reverse engineer gene regulatory networks was applied to time-course microarray data obtained from the developing mouse heart. Approximately 200 genes of interest were input into the algorithm to generate putative network topologies that are capable of explaining the experimental data via model simulation. To cull specious network interactions, thousands of putative networks are merged and filtered to generate scale-free, hierarchical networks that are statistically significant and biologically relevant. The networks are validated with known gene interactions and used to predict regulatory pathways important for the developing mammalian heart. Area under the precision-recall curve and receiver operator characteristic curve are 9% and 58%, respectively. Of the top 10 ranked predicted interactions, 4 have already been validated. The algorithm is further tested using a network enriched with known interactions and another depleted of them. The inferred networks contained more interactions for the enriched network versus the depleted network. In all test cases, maximum performance of the algorithm was achieved when the purely data-driven method of network inference was combined with a data-independent, functional-based association method. Lastly, the network generated from the list of approximately 200 genes of interest was expanded using gene-profile uniqueness metrics to include approximately 900 additional known mouse genes and to form the most likely cardiogenic gene regulatory network. The resultant network supports known regulatory interactions and contains several novel cardiogenic regulatory interactions. The method outlined herein provides an informative approach to network inference and leads to clear testable hypotheses related to gene regulation.