Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Altered contractility, Ca2+ transients, and cell morphology seen in a patient-specific iPSC-CM model of Ebstein's anomaly with left ventricular noncompaction.

Patients with two congenital heart diseases (CHDs), Ebstein's anomaly (EA) and left ventricular noncompaction (LVNC), suffer higher morbidity than either CHD alone. The genetic etiology and pathogenesis of combined EA/LVNC remain largely unknown. We investigated a familial EA/LVNC case associated with a variant (p.R237C) in the gene encoding Kelch-like protein 26 (KLHL26) by differentiating induced pluripotent stem cells (iPSCs) generated from affected and unaffected family members into cardiomyocytes (iPSC-CMs) and assessing iPSC-CM morphology, function, gene expression, and protein abundance. Compared with unaffected iPSC-CMs, CMs containing the KLHL26 (p.R237C) variant exhibited aberrant morphology including distended endo(sarco)plasmic reticulum (ER/SR) and dysmorphic mitochondria and aberrant function that included decreased contractions per minute, altered calcium transients, and increased proliferation. Pathway enrichment analyses based on RNASeq data indicated that the "structural constituent of muscle" pathway was suppressed, whereas the "ER lumen" pathway was activated. Taken together, these findings suggest that iPSC-CMs containing this KLHL26 (p.R237C) variant develop dysregulated ER/SR, calcium signaling, contractility, and proliferation.NEW & NOTEWORTHY We demonstrate here that iPSCs derived from patients with Ebstein's anomaly and left ventricular noncompaction, when differentiated into cardiomyocytes, display significant structural and functional changes that offer insight into disease pathogenesis, including altered ER/SR and mitochondrial morphology, contractility, and calcium signaling.

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.

GSEPD: a Bioconductor package for RNA-seq gene set enrichment and projection display.

BACKGROUND: RNA-seq, wherein RNA transcripts expressed in a sample are sequenced and quantified, has become a widely used technique to study disease and development. With RNA-seq, transcription abundance can be measured, differential expression genes between groups and functional enrichment of those genes can be computed. However, biological insights from RNA-seq are often limited by computational analysis and the enormous volume of resulting data, preventing facile and meaningful review and interpretation of gene expression profiles. Particularly, in cases where the samples under study exhibit uncontrolled variation, deeper analysis of functional enrichment would be necessary to visualize samples' gene expression activity under each biological function. RESULTS: We developed a Bioconductor package rgsepd that streamlines RNA-seq data analysis by wrapping commonly used tools DESeq2 and GOSeq in a user-friendly interface and performs a gene-subset linear projection to cluster heterogeneous samples by Gene Ontology (GO) terms. Rgsepd computes significantly enriched GO terms for each experimental condition and generates multidimensional projection plots highlighting how each predefined gene set's multidimensional expression may delineate samples. CONCLUSIONS: The rgsepd serves to automate differential expression, functional annotation, and exploratory data analyses to highlight subtle expression differences among samples based on each significant biological function.

Human genotyping and an experimental model reveal NPR-C as a possible contributor to morbidity in coarctation of the aorta.

Coarctation of the aorta (CoA) is a common congenital cardiovascular (CV) defect characterized by a stenosis of the descending thoracic aorta. Treatment exists, but many patients develop hypertension (HTN). Identifying the cause of HTN is challenging because of patient variability (e.g., age, follow-up duration, severity) and concurrent CV abnormalities. Our objective was to conduct RNA sequencing of aortic tissue from humans with CoA to identify a candidate gene for mechanistic studies of arterial dysfunction in a rabbit model of CoA devoid of the variability seen with humans. We present the first known evidence of natriuretic peptide receptor C (NPR-C; aka NPR3) downregulation in human aortic sections subjected to high blood pressure (BP) from CoA versus normal BP regions (validated to PCR). These changes in NPR-C, a gene associated with BP and proliferation, were replicated in the rabbit model of CoA. Artery segments from this model were used with human aortic endothelial cells to reveal the functional relevance of altered NPR-C activity. Results showed decreased intracellular calcium ([Ca2+]i) activity to C-type natriuretic peptide (CNP). Normal relaxation induced by CNP and atrial natriuretic peptide was impaired for aortic segments exposed to elevated BP from CoA. Inhibition of NPR-C (M372049) also impaired aortic relaxation and [Ca2+]i activity. Genotyping of NPR-C variants predicted to be damaging revealed that rs146301345 was enriched in our CoA patients, but sample size limited association with HTN. These results may ultimately be used to tailor treatment for CoA based on mechanical stimuli, genotyping, and/or changes in arterial function.

Ciliopathy variant burden and developmental delay in children with hypoplastic left heart syndrome.

PURPOSE: To test the hypothesis that patients with hypoplastic left heart syndrome (HLHS) and developmental delay will have a higher average summative C-score in ciliopathy genes than patients with HLHS without developmental delay. METHODS: Ciliopathy gene variant burden was determined utilizing a summative C-score for 14 ciliopathy genes in children with HLHS (n = 24). Mean summative C-scores were compared between children with and without developmental delay. Genome-wide randomizing gene sets were evaluated as a scoring control. RESULTS: Children with developmental delay had a mean summative C-score of 4.05 in ciliopathy genes as compared to a mean summative C-score of 2.02 for children without developmental delay. This difference in means was higher than 99.1% (empirical P value <0.01) of 2 million random lists of 14 genes. CONCLUSION: Genetically complex disorders such as ciliopathies can be assessed to determine phenotypic risk with summative C-score in appropriately chosen gene sets. If these results are replicated in subsequent cohorts, a diagnostic gene panel could identify risk for developmental delay and other ciliopathy-related comorbidities in infants with congenital heart disease.Genet Med advance online publication 27 October 2016Genetics in Medicine (2016); doi:10.1038/gim.2016.167.

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.

Activin-A and Bmp4 levels modulate cell type specification during CHIR-induced cardiomyogenesis.

The use of human pluripotent cell progeny for cardiac disease modeling, drug testing and therapeutics requires the ability to efficiently induce pluripotent cells into the cardiomyogenic lineage. Although direct activation of the Activin-A and/or Bmp pathways with growth factors yields context-dependent success, recent studies have shown that induction of Wnt signaling using low molecular weight molecules such as CHIR, which in turn induces the Activin-A and Bmp pathways, is widely effective. To further enhance the reproducibility of CHIR-induced cardiomyogenesis, and to ultimately promote myocyte maturation, we are using exogenous growth factors to optimize cardiomyogenic signaling downstream of CHIR induction. As indicated by RNA-seq, induction with CHIR during Day 1 (Days 0-1) was followed by immediate expression of Nodal ligands and receptors, followed later by Bmp ligands and receptors. Co-induction with CHIR and high levels of the Nodal mimetic Activin-A (50-100 ng/ml) during Day 0-1 efficiently induced definitive endoderm, whereas CHIR supplemented with Activin-A at low levels (10 ng/ml) consistently improved cardiomyogenic efficiency, even when CHIR alone was ineffective. Moreover, co-induction using CHIR and low levels of Activin-A apparently increased the rate of cardiomyogenesis, as indicated by the initial appearance of rhythmically beating cells by Day 6 instead of Day 8. By contrast, co-induction with CHIR plus low levels (3-10 ng/ml) of Bmp4 during Day 0-1 consistently and strongly inhibited cardiomyogenesis. These findings, which demonstrate that cardiomyogenic efficacy is improved by optimizing levels of CHIR-induced growth factors when applied in accord with their sequence of endogenous expression, are consistent with the idea that Nodal (Activin-A) levels toggle the entry of cells into the endodermal or mesodermal lineages, while Bmp levels regulate subsequent allocation into mesodermal cell types.

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.