Dynamic Role of trans Regulation of Gene Expression in Relation to Complex Traits.

Identifying causal genetic variants and understanding their mechanisms of effect on traits remains a challenge in genome-wide association studies (GWASs). In particular, how genetic variants (i.e., trans-eQTLs) affect expression of remote genes (i.e., trans-eGenes) remains unknown. We hypothesized that some trans-eQTLs regulate expression of distant genes by altering the expression of nearby genes (cis-eGenes). Using published GWAS datasets with 39,165 single-nucleotide polymorphisms (SNPs) associated with 1,960 traits, we explored whole blood gene expression associations of trait-associated SNPs in 5,257 individuals from the Framingham Heart Study. We identified 2,350 trans-eQTLs (at p < 10-7); more than 80% of them were found to have cis-associated eGenes. Mediation testing suggested that for 35% of trans-eQTL-trans-eGene pairs in different chromosomes and 90% pairs in the same chromosome, the disease-associated SNP may alter expression of the trans-eGene via cis-eGene expression. In addition, we identified 13 trans-eQTL hotspots, affecting from ten to hundreds of genes, suggesting the existence of master genetic regulators. Using causal inference testing, we searched causal variants across eight cardiometabolic traits (BMI, systolic and diastolic blood pressure, LDL cholesterol, HDL cholesterol, total cholesterol, triglycerides, and fasting blood glucose) and identified several cis-eGenes (ALDH2 for systolic and diastolic blood pressure, MCM6 and DARS for total cholesterol, and TRIB1 for triglycerides) that were causal mediators for the corresponding traits, as well as examples of trans-mediators (TAGAP for LDL cholesterol). The finding of extensive evidence of genome-wide mediation effects suggests a critical role of cryptic gene regulation underlying many disease traits.

Integrating genetic, transcriptional, and biological information provides insights into obesity.

OBJECTIVE: Indices of body fat distribution are heritable, but few genetic signals have been reported from genome-wide association studies (GWAS) of computed tomography (CT) imaging measurements of body fat distribution. We aimed to identify genes associated with adiposity traits and the key drivers that are central to adipose regulatory networks. SUBJECTS: We analyzed gene transcript expression data in blood from participants in the Framingham Heart Study, a large community-based cohort (n up to 4303), as well as implemented an integrative analysis of these data and existing biological information. RESULTS: Our association analyses identified unique and common gene expression signatures across several adiposity traits, including body mass index, waist-hip ratio, waist circumference, and CT-measured indices, including volume and quality of visceral and subcutaneous adipose tissues. We identified six enriched KEGG pathways and two co-expression modules for further exploration of adipose regulatory networks. The integrative analysis revealed four gene sets (Apoptosis, p53 signaling pathway, Proteasome, Ubiquitin-mediated proteolysis) and two co-expression modules with significant genetic variants and 94 key drivers/genes whose local networks were enriched with adiposity-associated genes, suggesting that these enriched pathways or modules have genetic effects on adiposity. Most identified key driver genes are involved in essential biological processes such as controlling cell cycle, DNA repair, and degradation of regulatory proteins are cancer related. CONCLUSIONS: Our integrative analysis of genetic, transcriptional, and biological information provides a list of compelling candidates for further follow-up functional studies to uncover the biological mechanisms underlying obesity. These candidates highlight the value of examining CT-derived and central adiposity traits.

A whole-blood transcriptome meta-analysis identifies gene expression signatures of cigarette smoking.

Cigarette smoking is a leading modifiable cause of death worldwide. We hypothesized that cigarette smoking induces extensive transcriptomic changes that lead to target-organ damage and smoking-related diseases. We performed a meta-analysis of transcriptome-wide gene expression using whole blood-derived RNA from 10,233 participants of European ancestry in six cohorts (including 1421 current and 3955 former smokers) to identify associations between smoking and altered gene expression levels. At a false discovery rate (FDR) <0.1, we identified 1270 differentially expressed genes in current vs. never smokers, and 39 genes in former vs. never smokers. Expression levels of 12 genes remained elevated up to 30 years after smoking cessation, suggesting that the molecular consequence of smoking may persist for decades. Gene ontology analysis revealed enrichment of smoking-related genes for activation of platelets and lymphocytes, immune response, and apoptosis. Many of the top smoking-related differentially expressed genes, including LRRN3 and GPR15, have DNA methylation loci in promoter regions that were recently reported to be hypomethylated among smokers. By linking differential gene expression with smoking-related disease phenotypes, we demonstrated that stroke and pulmonary function show enrichment for smoking-related gene expression signatures. Mediation analysis revealed the expression of several genes (e.g. ALAS2) to be putative mediators of the associations between smoking and inflammatory biomarkers (IL6 and C-reactive protein levels). Our transcriptomic study provides potential insights into the effects of cigarette smoking on gene expression in whole blood and their relations to smoking-related diseases. The results of such analyses may highlight attractive targets for treating or preventing smoking-related health effects.

A meta-analysis of gene expression signatures of blood pressure and hypertension.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

Transcriptome Assessment of Erythema Migrans Skin Lesions in Patients With Early Lyme Disease Reveals Predominant Interferon Signaling.

Background: The most common clinical manifestation of early Lyme disease is the erythema migrans (EM) skin lesion that develops at the tick bite site typically between 7 and 14 days after infection with Borreliella burgdorferi. The host-pathogen interactions that occur in the skin may have a critical role in determining outcome of infection. Methods: Gene arrays were used to characterize the global transcriptional alterations in skin biopsy samples of EM lesions from untreated adult patients with Lyme disease in comparison to controls. Results: The transcriptional pattern in EM biopsies consisted of 254 differentially regulated genes (180 induced and 74 repressed) characterized by the induction of chemokines, cytokines, Toll-like receptors, antimicrobial peptides, monocytoid cell activation markers, and numerous genes annotated as interferon (IFN)-inducible. The IFN-inducible genes included 3 transcripts involved in tryptophan catabolism (IDO1, KMO, KYNU) that play a pivotal role in immune evasion by certain other microbial pathogens by driving the differentiation of regulatory T cells. Conclusions: This is the first study to globally assess the human skin transcriptional response during early Lyme disease. Borreliella burgdorferi elicits a predominant IFN signature in the EM lesion, suggesting a potential mechanism for spirochetal dissemination via IDO1-mediated localized immunosuppression.

Messenger RNA and MicroRNA transcriptomic signatures of cardiometabolic risk factors.

BACKGROUND: Cardiometabolic (CM) risk factors are heritable and cluster in individuals. We hypothesized that CM risk factors are associated with multiple shared and unique mRNA and microRNA (miRNA) signatures. We examined associations of mRNA and miRNA levels with 6 CM traits: body mass index, HDL-cholesterol and triglycerides, fasting glucose, and systolic and diastolic blood pressures through cross-sectional analysis of 2812 Framingham Heart Study who had whole blood collection for RNA isolation for mRNA and miRNA expression studies and who consented to genetic research. We excluded participants taking medication for hypertension, dyslipidemia, or diabetes. We measured mRNA (n = 17,318; using the Affymetrix GeneChip Human Exon 1.0 ST Array) and miRNA (n = 315; using qRT-PCR) expression in whole blood. We used linear regression for mRNA analyses and a combination of linear and logistic regression for miRNA analyses. We conducted miRNA-mRNA coexpression and gene ontology enrichment analyses to explore relations between pleiotropic miRNAs, mRNA expression, and CM trait clustering. RESULTS: We identified hundreds of significant associations between mRNAs, miRNAs, and individual CM traits. Four mRNAs (FAM13A, CSF2RB, HIST1H2AC, WNK1) were associated with all 6 CM traits (FDR < 0.001) and four miRNAs (miR-197-3p, miR-328, miR-505-5p, miR-145-5p) were associated with four CM traits (FDR < 0.05). Twelve mRNAs, including WNK1, that were coexpressed with the four most pleiotropic miRNAs, were also miRNA targets. mRNAs coexpressed with pleiotropic miRNAs were enriched for RNA metabolism (miR-505-5p), ubiquitin-dependent protein catabolism (miR-197-3p, miR-328) and chromatin assembly (miR-328). CONCLUSIONS: We identified mRNA and miRNA signatures of individual CM traits and their clustering. Implicated transcripts may play causal roles in CM risk or be downstream consequences of CM risk factors on the transcriptome. Studies are needed to establish whether or not pleiotropic circulating transcripts illuminate causal pathways for CM risk.

Integrative DNA, RNA, and protein evidence connects TREML4 to coronary artery calcification.

Coronary artery calcification (CAC) is a heritable and definitive morphologic marker of atherosclerosis that strongly predicts risk for future cardiovascular events. To search for genes involved in CAC, we used an integrative transcriptomic, genomic, and protein expression strategy by using next-generation DNA sequencing in the discovery phase with follow-up studies using traditional molecular biology and histopathology techniques. RNA sequencing of peripheral blood from a discovery set of CAC cases and controls was used to identify dysregulated genes, which were validated by ClinSeq and Framingham Heart Study data. Only a single gene, TREML4, was upregulated in CAC cases in both studies. Further examination showed that rs2803496 was a TREML4 cis-eQTL and that the minor allele at this locus conferred up to a 6.5-fold increased relative risk of CAC. We characterized human TREML4 and demonstrated by immunohistochemical techniques that it is localized in macrophages surrounding the necrotic core of coronary plaques complicated by calcification (but not in arteries with less advanced disease). Finally, we determined by von Kossa staining that TREML4 colocalizes with areas of microcalcification within coronary plaques. Overall, we present integrative RNA, DNA, and protein evidence implicating TREML4 in coronary artery calcification. Our findings connect multimodal genomics data with a commonly used clinical marker of cardiovascular disease.

Subgingival microbial communities in Leukocyte Adhesion Deficiency and their relationship with local immunopathology.

Leukocyte Adhesion Deficiency I (LAD-I) is a primary immunodeficiency caused by single gene mutations in the CD18 subunit of ß2 integrins which result in defective transmigration of neutrophils into the tissues. Affected patients suffer from recurrent life threatening infections and severe oral disease (periodontitis). Microbial communities in the local environment (subgingival plaque) are thought to be the triggers for inflammatory periodontitis, yet little is known regarding the microbial communities associated with LAD-I periodontitis. Here we present the first comprehensive characterization of the subgingival communities in LAD-I, using a 16S rRNA gene-based microarray, and investigate the relationship of this tooth adherent microbiome to the local immunopathology of periodontitis. We show that the LAD subgingival microbiome is distinct from that of health and Localized Aggressive Periodontitits. Select periodontitis-associated species in the LAD microbiome included Parvimonas micra, Porphyromonas endodontalis, Eubacterium brachy and Treponema species. Pseudomonas aeruginosa, a bacterium not typically found in subgingival plaque is detected in LAD-I. We suggest that microbial products from LAD-associated communities may have a role in stimulating the local inflammatory response. We demonstrate that bacterial LPS translocates into the lesions of LAD-periodontitis potentially triggering immunopathology. We also show in in vitro assays with human macrophages and in vivo in animal models that microbial products from LAD-associated subgingival plaque trigger IL-23-related immune responses, which have been shown to dominate in patient lesions. In conclusion, our current study characterizes the subgingival microbial communities in LAD-periodontitis and supports their role as triggers of disease pathogenesis.

Integromic analysis of genetic variation and gene expression identifies networks for cardiovascular disease phenotypes.

BACKGROUND: Cardiovascular disease (CVD) reflects a highly coordinated complex of traits. Although genome-wide association studies have reported numerous single nucleotide polymorphisms (SNPs) to be associated with CVD, the role of most of these variants in disease processes remains unknown. METHODS AND RESULTS: We built a CVD network using 1512 SNPs associated with 21 CVD traits in genome-wide association studies (at P≤5×10(-8)) and cross-linked different traits by virtue of their shared SNP associations. We then explored whole blood gene expression in relation to these SNPs in 5257 participants in the Framingham Heart Study. At a false discovery rate <0.05, we identified 370 cis-expression quantitative trait loci (eQTLs; SNPs associated with altered expression of nearby genes) and 44 trans-eQTLs (SNPs associated with altered expression of remote genes). The eQTL network revealed 13 CVD-related modules. Searching for association of eQTL genes with CVD risk factors (lipids, blood pressure, fasting blood glucose, and body mass index) in the same individuals, we found examples in which the expression of eQTL genes was significantly associated with these CVD phenotypes. In addition, mediation tests suggested that a subset of SNPs previously associated with CVD phenotypes in genome-wide association studies may exert their function by altering expression of eQTL genes (eg, LDLR and PCSK7), which in turn may promote interindividual variation in phenotypes. CONCLUSIONS: Using a network approach to analyze CVD traits, we identified complex networks of SNP-phenotype and SNP-transcript connections. Integrating the CVD network with phenotypic data, we identified biological pathways that may provide insights into potential drug targets for treatment or prevention of CVD.

Sex- and age-interacting eQTLs in human complex diseases.

Many complex human diseases exhibit sex or age differences in gene expression. However, the presence and the extent of genotype-specific variations in gene regulation are largely unknown. Here, we report results of a comprehensive analysis of expression regulation of genetic variation related to 11,672 complex disease-associated SNPs as a function of sex and age in whole-blood-derived RNA from 5254 individuals. At false discovery rate <0.05, we identified 14 sex- and 10 age-interacting expression quantitative trait loci (eQTLs). We show that these eQTLs are also associated with many sex- or age-associated traits. These findings provide important context regarding the regulation of phenotypes by genotype-environment interaction.

Integrative network analysis reveals molecular mechanisms of blood pressure regulation.

Genome-wide association studies (GWAS) have identified numerous loci associated with blood pressure (BP). The molecular mechanisms underlying BP regulation, however, remain unclear. We investigated BP-associated molecular mechanisms by integrating BP GWAS with whole blood mRNA expression profiles in 3,679 individuals, using network approaches. BP transcriptomic signatures at the single-gene and the coexpression network module levels were identified. Four coexpression modules were identified as potentially causal based on genetic inference because expression-related SNPs for their corresponding genes demonstrated enrichment for BP GWAS signals. Genes from the four modules were further projected onto predefined molecular interaction networks, revealing key drivers. Gene subnetworks entailing molecular interactions between key drivers and BP-related genes were uncovered. As proof-of-concept, we validated SH2B3, one of the top key drivers, using Sh2b3(-/-) mice. We found that a significant number of genes predicted to be regulated by SH2B3 in gene networks are perturbed in Sh2b3(-/-) mice, which demonstrate an exaggerated pressor response to angiotensin II infusion. Our findings may help to identify novel targets for the prevention or treatment of hypertension.

Dissecting the roles of microRNAs in coronary heart disease via integrative genomic analyses.

OBJECTIVE: The roles of microRNAs (miRNAs) in coronary heart disease (CHD) have not been well characterized. This study sought to systematically characterize the complex genomic architecture of CHD by integrating whole blood miRNA and mRNA expression with genetic variation in 186 CHD cases and 186 controls. APPROACH AND RESULTS: At false discovery rate <0.2, 15 miRNAs were differentially expressed between CHD cases and controls. To explore regulatory mechanisms, we integrated miRNA and mRNA expression with genome-wide genotype data to investigate miRNA and mRNA associations and relationships of genetic variation with miRNAs. We identified a large number of correlated miRNA-mRNA pairs and genetic loci that seem to regulate miRNA levels. Subsequently, we explored the relationships of these complex molecular associations with CHD status. We identified a large difference in miRNA-mRNA associations between CHD cases and controls, as demonstrated by a significantly higher proportion of inversely correlated miRNA-mRNA pairs in cases versus controls (80% versus 30%; P<1×10(-16)), suggesting a genome-wide shift in the regulatory structure of the transcriptome in CHD. The differentially coexpressed miRNA-mRNA pairs showed enrichment for CHD risk genetic variants affecting both miRNA and mRNA expression levels, implicating a putatively causal role in CHD. Furthermore, 3 miRNAs (miR-1275, miR-365a-3p, and miR-150-5p) were associated with an mRNA coexpression module that was causally linked to CHD and reflected the dysregulation of B-cell centered immune function. CONCLUSIONS: Our results provide novel evidence that miRNAs are important regulators of biological processes involved in CHD via genetic control and via their tight coexpression with mRNAs.

A systematic heritability analysis of the human whole blood transcriptome.

Genome-wide expression quantitative trait locus (eQTL) mapping may reveal common genetic variants regulating gene expression. In addition to mapping eQTLs, we systematically evaluated the heritability of the whole blood transcriptome in 5,626 participants from the Framingham Heart Study. Of all gene expression measurements, about 40 % exhibit evidence of being heritable [hgeneExp(2) > 0, (p < 0.05)], the average heritability was estimated to be 0.13, and 10 % display hgeneExp(2) > 0.2. To identify the role of eQTLs in promoting phenotype differences and disease susceptibility, we investigated the proportion of cis/trans eQTLs in different heritability categories and discovered that genes with higher heritability are more likely to have cis eQTLs that explain large proportions of variance in the expression of the corresponding genes. Single cis eQTLs explain 0.33-0.53 of variance in transcripts on average, whereas single trans eQTLs only explain 0.02-0.07. The top cis eQTLs tend to explain more variance in the corresponding gene when its hgeneExp(2) is greater. Taking body mass index (BMI) as a case study, we cross-linked cis/trans eQTLs with both GWAS SNPs and differentially expressed genes for BMI. We discovered that BMI GWAS SNPs in 16p11.2 (e.g., rs7359397) are associated with several BMI differentially expressed genes in a cis manner (e.g. SULT1A1, SPNS1, and TUFM). These BMI signature genes explain a much larger proportion of variance in BMI than do the GWAS SNPs. Our results shed light on the impact of eQTLs on the heritability of the human whole blood transcriptome and its relations to phenotype differences.

Circadian changes in long noncoding RNAs in the pineal gland.

Long noncoding RNAs (lncRNAs) play a broad range of biological roles, including regulation of expression of genes and chromosomes. Here, we present evidence that lncRNAs are involved in vertebrate circadian biology. Differential night/day expression of 112 lncRNAs (0.3 to >50 kb) occurs in the rat pineal gland, which is the source of melatonin, the hormone of the night. Approximately one-half of these changes reflect nocturnal increases. Studies of eight lncRNAs with 2- to >100-fold daily rhythms indicate that, in most cases, the change results from neural stimulation from the central circadian oscillator in the suprachiasmatic nucleus (doubling time = 0.5-1.3 h). Light exposure at night rapidly reverses (halving time = 9-32 min) levels of some of these lncRNAs. Organ culture studies indicate that expression of these lncRNAs is regulated by norepinephrine acting through cAMP. These findings point to a dynamic role of lncRNAs in the circadian system.

Whole blood gene expression and interleukin-6 levels.

BACKGROUND: Circulating interleukin-6 levels increase with advancing age and are a risk factor for various diseases and mortality. The characterization of gene expression profiles associated with interleukin-6 levels might suggest important molecular events underlying its regulation. METHODS AND RESULTS: We studied the association of transcriptional profiles with interleukin-6 levels in 2422 participants from the Framingham Heart Study Offspring Cohort using Affymetrix Human Exon 1.0 ST Array. We identified 4139 genes that were significantly associated with interleukin-6 levels (FDR<0.05) after adjusting for age, sex and blood cell components. We then replicated 807 genes in the InCHIANTI study with 694 participants. Many of the top genes are involved in inflammation-related pathways or erythrocyte function, including JAK/Stat signaling pathway and interleukin-10 signaling pathway. CONCLUSION: We identified and replicated 807 genes that were associated with circulating interleukin-6 levels. Future characterization of interleukin-6 regulation networks may facilitate the identification of additional potential targets for treating inflammation-related diseases.

Identification of candidate genes involved in coronary artery calcification by transcriptome sequencing of cell lines.

BACKGROUND: Massively-parallel cDNA sequencing (RNA-Seq) is a new technique that holds great promise for cardiovascular genomics. Here, we used RNA-Seq to study the transcriptomes of matched coronary artery disease cases and controls in the ClinSeq® study, using cell lines as tissue surrogates. RESULTS: Lymphoblastoid cell lines (LCLs) from 16 cases and controls representing phenotypic extremes for coronary calcification were cultured and analyzed using RNA-Seq. All cell lines were then independently re-cultured and along with another set of 16 independent cases and controls, were profiled with Affymetrix microarrays to perform a technical validation of the RNA-Seq results. Statistically significant changes (p < 0.05) were detected in 186 transcripts, many of which are expressed at extremely low levels (5-10 copies/cell), which we confirmed through a separate spike-in control RNA-Seq experiment. Next, by fitting a linear model to exon-level RNA-Seq read counts, we detected signals of alternative splicing in 18 transcripts. Finally, we used the RNA-Seq data to identify differential expression (p < 0.0001) in eight previously unannotated regions that may represent novel transcripts. Overall, differentially expressed genes showed strong enrichment (p = 0.0002) for prior association with cardiovascular disease. At the network level, we found evidence for perturbation in pathways involving both cardiovascular system development and function as well as lipid metabolism. CONCLUSIONS: We present a pilot study for transcriptome involvement in coronary artery calcification and demonstrate how RNA-Seq analyses using LCLs as a tissue surrogate may yield fruitful results in a clinical sequencing project. In addition to canonical gene expression, we present candidate variants from alternative splicing and novel transcript detection, which have been unexplored in the context of this disease.

Gene expression signatures of coronary heart disease.

OBJECTIVE: To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study. APPROACH AND RESULTS: A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell-related pathways with CHD. CONCLUSIONS: Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

A systems biology framework identifies molecular underpinnings of coronary heart disease.

OBJECTIVE: Genetic approaches have identified numerous loci associated with coronary heart disease (CHD). The molecular mechanisms underlying CHD gene-disease associations, however, remain unclear. We hypothesized that genetic variants with both strong and subtle effects drive gene subnetworks that in turn affect CHD. APPROACH AND RESULTS: We surveyed CHD-associated molecular interactions by constructing coexpression networks using whole blood gene expression profiles from 188 CHD cases and 188 age- and sex-matched controls. Twenty-four coexpression modules were identified, including 1 case-specific and 1 control-specific differential module (DM). The DMs were enriched for genes involved in B-cell activation, immune response, and ion transport. By integrating the DMs with gene expression-associated single-nucleotide polymorphisms and with results of genome-wide association studies of CHD and its risk factors, the control-specific DM was implicated as CHD causal based on its significant enrichment for both CHD and lipid expression-associated single-nucleotide polymorphisms. This causal DM was further integrated with tissue-specific Bayesian networks and protein-protein interaction networks to identify regulatory key driver genes. Multitissue key drivers (SPIB and TNFRSF13C) and tissue-specific key drivers (eg, EBF1) were identified. CONCLUSIONS: Our network-driven integrative analysis not only identified CHD-related genes, but also defined network structure that sheds light on the molecular interactions of genes associated with CHD risk.