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.

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.