Skip to main content
  • Research article
  • Open access
  • Published:

Messenger RNA and MicroRNA transcriptomic signatures of cardiometabolic risk factors

Abstract

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.

Background

Metabolic risk factors cluster in individuals and their presence is associated with increased risk for type II diabetes mellitus (T2DM) and cardiovascular disease (CVD) [1, 2]. Genome-wide association studies (GWAS) have identified hundreds of loci associated with cardiometabolic (CM) risk factors including body mass index (BMI), lipid levels, glucose levels, T2DM, and blood pressure [15]. In more recent years, circulating mRNA and miRNA transcriptional patterns have been linked to CVD and CM phenotypes [69].

Despite having identified molecular associations with CM risk factors at the genetic, mRNA, and miRNA levels, multidimensional interrelations of these molecular elements and how they interact to influence susceptibility to CM risk factors and CVD risk remain unknown. Recent studies have shown that circulating mRNA and miRNA transcript levels are heritable quantitative traits that are partly under genetic control [10]. The molecular underpinnings of complex CM diseases may be explained in part by genetic variation, mRNA and miRNA expression, and by miRNA–mRNA interaction [11]. Therefore, integrative analyses that incorporate multidimensional genomic data are necessary to investigate and characterize complex changes in the regulatory machinery and their effects on biological functions and complex CM phenotypes [6, 7].

With the goal of generating new insights into potential gene regulatory factors responsible for the clustering of CM risk factors, we examined the relations of circulating mRNA and miRNA expression levels to six CM traits: BMI, plasma lipid levels (HDL cholesterol [HDL-C], triglycerides [TG]), fasting glucose levels, and systolic (SBP) and diastolic (DBP) blood pressure. A high-throughput, unbiased approach was used to detect novel relationships among mRNAs and miRNAs across multiple CM traits since clustering of these traits is frequently seen clinically.

Methods

Study sample

The Framingham Heart Study (FHS) is a prospective, community-based observational study of CVD and its risk factors. We included participants from the FHS Offspring and Third Generation cohorts [12, 13]. We focused on 725 Offspring cohort participants (examination 8, 2005–2008) and 2087 Third Generation cohort participants (examination 2, 2008–2011) 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. Venous blood samples were obtained after an overnight fast and samples were stored using methods that maintain RNA stability (http://www.preanalytix.com/products/blood/RNA/paxgene-blood-rna-tube) [14]. All participants gave informed consent. The Boston University Medical Center Institutional Review Board approved FHS examination protocols and University of Massachusetts Medical School Review Board approved the miRNA and RNA profiling protocols.

Risk factor definitions

FHS participants had a physician-administered history and physical examination including anthropometric measurements and a laboratory evaluation focused on CVD and its risk factors. Blood pressure was measured twice by a physician with the participant seated; the average of both measurements was used to calculate SBP and DBP. BMI was calculated by dividing the weight in kilograms by the square of height in meters (kg/m2). Fasting plasma glucose was measured using a hexokinase reagent kit (A-gent glucose test, Abbott Laboratories, Inc., South Pasadena, CA); the intra-assay CV was <3% [15]. Venous blood samples were collected in 0.1% EDTA tubes and plasma was separated by centrifugation. Plasma lipids levels were measured before freezing. Triglycerides were measured using an automated enzymatic assay. HDL-C was measured after dextran sulfate magnesium precipitation [16].

mRNA expression profiling

Whole blood was collected in PAXgene (QIAGEN, Valencia, CA) tubes from each study participant after an overnight fast and stored at −80°C. RNA was extracted from whole blood using the PAXgene Blood RNA System Kit according to published methods [14]. RNA quality was validated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA); A NanoDrop ND-1000 spectrophotometer was used to quantify RNA concentration (NanoDrop Technologies, Wilmington, DE). NuGEN’s WT-Ovation Pico RNA Amplification System was used to amplify 50 ng of total RNA, which was then labeled according to established protocols with FL-Ovation cDNA Biotin Module V2 (NuGEN, San Carlos, CA) [11].

miRNA expression profiling

The same RNA sample was used for miRNA isolation. The high throughput Gene Expression and Biomarker Core Laboratory at the University of Massachusetts Medical School profiled 346 miRNAs isolated from whole blood (RNA isolation was performed by Asuragen, Inc, Austin, TX) in 2445 FHS Offspring and 3245 Third Generation cohort participants using TaqMan chemistry-based assays (Additional file 1: Methods). The initial miRNA list encompassed all TaqMan miRNA assays (774) available at the start of the study. If a miRNA was not expressed in any of 550 randomly selected samples, this miRNA was not examined in the entire cohort. Using this method, we restricted miRNA quantification in the overall cohort to the aforementioned 346 miRNAs. Methods used for cDNA conversion, preamplification, and quantification are reported in the Additional file 1: Methods. Among 70 replicate samples, >95% of the data points had coefficients of variation <10% (mean ~4%). As described previously, miRNA expression was quantified using cycle threshold (Ct), where higher Ct values reflect lower miRNA expression [14]. We analyzed 315 miRNAs that were expressed in at least 100 people.

Statistical analyses

We report descriptive statistics using counts and percentages for binary variables and means ± standard deviations (SD) for continuous variables. We modeled mRNAs and miRNAs as response variables versus each CM risk factor, adjusting for age and sex. mRNA expression was quantified by log-2 transformed expression intensities. These models also adjusted for RNA processing variables and differential cell counts (imputed). RNA processing variables included isolation batch, RNA quality, concentration, and 260/280 ratio (ratio of absorbance at 260 and 280nm using a spectrophotometer) [14]. Differential cell counts (white blood count, and percent lymphocytes, monocytes, eosinophils, basophils) were imputed from mRNA expression values via partial least square (PLS) prediction, with cross-validated prediction accuracy estimates ranging from 0.25–0.89. We conducted modeling in the full sample and in two random, equal subsets (discovery and validation, Additional file 1), which preserved intact pedigrees. For each metabolic risk factor, we applied false discovery rate (FDR) calculations in the whole sample, the discovery set, and the validation set.

mRNA analyses

The robust multichip analysis (RMA) algorithm [17] was applied using Affymetrix Power Tools (APT) for generation of signal values to yield an initially normalized dataset using log-2 transformed expression intensities. For a detailed description on data quality control and normalization, please refer to Joehanes et al. [18] Further adjustment of this dataset by Affymetrix quality control parameters (all_probeset_mean, all_probeset_stdev, neg_control_mean, neg_control_stdev, pos_control_mean, pos_control_stdev, all_probeset_rle_mean, all_probeset_mad_residual_mean, mm_mean), the first principal component of the dataset (PC1), batch (Batch_Lump), and a factor accounting for the non-random layout of probesets on the array (ProbesetGroupDiff) yielded a final normalized dataset. A pedigree-based mixed-effects model implemented in the R package "pedigreemm (version 1.0-4)" was applied to this final normalized dataset for analysis of differential gene expression. Age, sex, measured or imputed blood cell counts (RBC, WBC, PLT, LY_PER, MO_PER, EO_PER, BA_PER, Retic_Per) were included in the model as covariates.

miRNA analyses

We used pedigree-based linear mixed-effects models (R package “lmekin”) to analyze continuous miRNA values (i.e., when Ct < 27) and logistic regression models to analyze binary values (i.e., Ct < 27 versus Ct > = 27). This was necessary because expression of miRNAs was not universal, but varied from 1.5–99.9% among miRNAs. The observed Ct values generally did not have a truncated normal distribution, which precluded Tobit modeling [19]. Furthermore, imputation produced extreme bi-modal distributions and was not acceptable for response data in linear modeling. Therefore, we employed an adaptive approach. For a given miRNA, if at least 90% of participants expressed it, we used the linear-model p value; if <10% of participants expressed it, we used the logistic-model p value; if between 10% and 90% of participants expressed the miRNA, we combined results from the two models. Specifically, we added their X 2 statistics and we calculated the p-value from the distribution of a X 2 variate with two degrees of freedom.

miRNA-mRNA coexpression analysis

The coexpression analysis was performed on FHS samples for which miRNA and mRNA data were both available (N = 5626). Linear mixed models (R package “lmekin”) were used to conduct pairwise coexpression analyses for all profiled mRNAs (dependent variable, N = 17,318) and 280 miRNAs (independent variable) expressed in >200 samples, with fixed effects including age, sex, technical covariates, imputed cell types, surrogate variables (SV), and a random effect to account for family structure. As described above, the mRNA expression was quantified by log-2 transformed expression intensities. miRNA expression used Ct values with higher values reflected lower expression levels of miRNAs. Adjustment was made for technical covariates (11 for mRNA expression and 4 for miRNA expression). Surrogate variables (SVs) were computed from the mRNA expression data using the R package “SVA,” and 51 SVs associated with at least 1 miRNA at Bonferroni corrected P < 1.7 × 10−4 (0.05/280) were included in the statistical model. We chose SVs that correlated with at least 1 miRNA to adjust for additional hidden effects in the mRNA expression measurements that might have affected miRNA-mRNA correlations. The Benjamini-Hochberg method [20] was used to compute the false discovery rate (FDR). The significant miRNA-mRNA coexpression pairs were selected using FDR < 0.05.

miRNA target databases

For the top four multiple-trait associated miRNAs, we used OmniSearch to search for computationally predicted and experimentally validated targets. OmniSearch is a semantics-based integration system to search miRNA targets. OmniSearch provides access to multiple prediction databases, including miRDB, TargetScan, and miRanda, and experimentally validated including miRTarBase. The predicted miRNA-mRNA target pairs from these databases were compared with the miRNA-mRNA coexpression pairs identified in the current study. The predicted miRNA-mRNA target pairs from these databases were compared with the miRNA-mRNA coexpression pairs identified in the current study.

Pathway and gene ontology enrichment analysis

Co-expressed mRNAs for each miRNA were combined as a set and classified using Gene Ontology (GO) databases to identify potentially relevant biological processes. Fisher’s exact test was used to calculate enrichment p values of the overlapped genes in comparison with the number of co-expressed mRNAs and the number of genes in each GO biological process terms. Because of many GO terms duplicated, we only used 825 unique GO biological process terms as suggested by MsigDB [21]. We used a Bonferroni adjusted statistical significance threshold of p < 0.05/825 = 6.0 × 10−5.

All statistical analyses were performed using SAS software version 9.2 (SAS Institute Inc., Cary, NC, USA) or R software version 3.1.1 (R Foundation for Statistical Computing, Vienna, Austria).

Results

Demographic and CM risk factor characteristics of the 2812 participants are shown in Table 1. The study sample consisted of middle-aged (49 ± 12 years) participants, a slight majority of whom were women (59%).

Table 1 Framingham Heart Study Offspring and Generation 3 Study participant characteristicsa

Circulating mRNAs in association with metabolic traits

Each CM trait was associated with multiple mRNA gene transcripts at FDR <0.05. TG was associated with the greatest number of circulating mRNAs (N = 5049), followed by BMI (N = 4826), HDL-C (N = 1768), DBP (N = 1499), SBP (N = 1019), and glucose (N = 1014).

Numerous transcripts were associated with multiple CM traits (Table 2 and Fig. 1). As shown in Fig. 1, genes associated with BMI shared associations with other CM traits, most notably DBP, TG, and HDL-C. Genes associated with both BMI and TG and those associated with BMI and DBP were positively correlated, whereas genes associated with both BMI and HDL-C were, in general, inversely correlated. SBP and DBP shared a large number of associated genes and the directionality of shared transcripts was concordant.

Table 2 Eighteen mRNAsa with greatest pleiotropy across metabolic traits using a cut-off of FDR <0.001 to define significance of association
Fig. 1
figure 1

Similarity networks for whole blood miRNA (panel A) and mRNA (panel B) across cardiometabolic traits

Four circulating transcripts were associated with all six CM risk factors at an FDR <0.001 (Table 2), including several genes with known associations with cardiovascular and/or pulmonary disease: FAM13A (family with sequence similarity 13, member a) [22], CSF2RB (colony stimulating factor 2 receptor, beta) [23], HIST1H2AC (histone cluster 1, H2ac), and WNK1 (WNK lysine deficient protein kinase 1) [4, 24]. An additional 14 transcripts were associated with five of the six CM traits (Table 2).

We conducted an analysis of all four of the most pleiotropic genes to examine relations to ‘energy metabolism’ or ‘cardiomyocyte function’ by searching these key words using the GeneRif database (https://www.ncbi.nlm.nih.gov/gene/about-generif), and found none of the four aforementioned mRNAs to have reported functions related to cardiomyocyte function or energy metabolism.

The top ten gene transcripts associated with each CM trait (BMI, HDL-C, TG, glucose, and SBP and DBP) are shown in Additional file 1: Tables S3–S8. The directionality and strength of association of the top four mRNAs are graphically depicted in Additional file 2: Figure S1B.

Circulating miRNAs in association with cardiometabolic traits

Five CM traits were associated with multiple miRNAs (FDR < 0.05). TG was associated with the greatest number of miRNAs (N = 150), followed by DBP (N = 112), BMI (N = 99), SBP (N = 4). HDL-C and glucose were not associated with any circulating miRNAs at this FDR threshold.

Fifty miRNAs were associated with three or more CM traits and four were associated with four CM traits (miR-197-3p, miR-328, miR-505-5p, miR-145-5p) (Table 3 and Fig. 1). As shown in Fig. 1, miRNAs associated with BMI shared associations with other CM traits, most notably DBP and TG. miRNAs associated with both BMI and TG, and those associated with BMI and DBP, were positively correlated (reflected by line color, Fig. 1).

Table 3 Four miRNAsa with greatest pleiotropy across metabolic traits using a cut-off of FDR <0.05 to define significance of association

The top ten miRNAs associated with each CM trait (BMI, HDL-C, TG, glucose, and SBP and DBP) are shown in Additional file 1: Tables S9–S14. The directionality and strength of association of the top four miRNAs are graphically depicted in Additional file 2: Figure S1A.

Replication results

Results of the separate analyses of discovery and validation sets revealed a high degree of concordance for mRNA results across all traits (Additional file 1: Table S1). Only two genes, ARRDC3 and CAPN2, which were associated with glucose in the discovery set, failed to validate (at FDR <0.05) in our validation set; their validation set FDR values were 0.06. In contrast, owing to reduced power to detect miRNA-trait associations, far less concordance was noted between the discovery and replication sets for miRNA-trait associations (Additional file 1: Table S2).

Coexpression analysis

We identified coexpressed mRNAs for the four miRNAs that were associated with four CM traits (Table 3). These highly pleiotropic miRNAs were associated with a large number of mRNAs (1109 mRNAs in total; 396 coexpressed mRNAs for miR-505-5p, 241 for miR-197-3p, 177 for miR-145-5p, and 649 for miR-328). Notably, similar patterns of associations across traits were seen for mRNAs and miRNAs (Fig. 2).

Fig. 2
figure 2

Coexpression network analysis

Among the 1109 coexpressed mRNAs, 807 mRNAs were associated with at least one CM trait at FDR <0.1, and 249 were associated with more than three traits at FDR <0.1. A less restrictive FDR threshold was used in these analyses in light of the relatively low number of coexpressed mRNAs and miRNAs with CM trait associations included in this model. The coexpressed mRNAs for the four highly pleiotropic miRNAs were highly enriched for associations with CM traits (enrichment P < 1 × 10−32 by hypergeometric test). Figure 3 shows the miRNA-mRNA coexpression network for the four highly pleiotropic miRNAs and their coexpressed mRNAs (n = 249) that were each associated with at least four CM traits. Among the coexpressed mRNAs, 17 mRNAs were also miRNA targets reported by at least one miRNA target database, including WNK1 for miR-197-3p.

Fig. 3
figure 3

Pleotropic RNAs and miRNAs with heterogeneous effect directions across cardiometabolic traits. a Within trait associations showing heterogenous effect directions for pleiotropic mRNAs and miRNAs. b Heat map showing strength and directionality of associations between pleiotropic mRNA and miRNAs

Gene ontology enrichment analysis (Additional file 1: Table S15) revealed that the co-expressed mRNAs for miR-505-5p were enriched for genes involved in RNA metabolism (P = 3.5 × 10−5). The coexpressed mRNAs for miR-197-3p and miR-328 were enriched for cellular macromolecule catabolism (P = 7.67 × 10−7) and ubiquitin-dependent protein catabolism (P = 1.86 × 10−10). The coexpressed mRNAs for miR-328 were enriched for DNA packaging and chromatin assembly (P = 1.73 × 10−5). miRNAs coexpressed with miR-145 did not show significant enrichment for GO terms.

Discussion

In a large, community-based cohort, we identified distinct as well as shared circulating transcriptomic signatures for CM risk factors. Four mRNAs were associated with all 6 CM risk factors (FDR <0.0001) and four miRNAs were associated with four CM risk factors (FDR <0.05). miRNAs associated with the greatest number of CM traits were coexpressed with many of the mRNAs that associated with multiple CM traits, including validated miR targets. Results of GO analyses revealed enrichment for processes relevant to regulation of gene expression and protein levels. These findings are consistent with causal roles of the implicated genes in CM risk. Alternatively, they may be due to downstream consequences of metabolic syndrome on the transcriptomic landscape. Functional studies are warranted to explore the mechanistic role of altered mRNA and miRNA expression in the pathogenesis of CVD and its CM risk factors.

To our knowledge, no prior study has performed a large-scale analysis of circulating miRNA and mRNA expression across CM traits. Several small studies, however, have examined mRNAs or miRNAs in relation to CM risk factors or disease. In the Young Finns Study [25], which included 71 participants, several circulating miRs were found to be associated with components of metabolic syndrome, including glucose and lipids. The authors also showed that down-regulated targets of two miRs, miR-1207-5p and miR-129-2-3p, were enriched in PI3K and MAPK pathways and that eight of 12 enriched pathways were downregulated in individuals with metabolic syndrome. Reflecting the translational relevance of our findings and validity of our approach, as discussed below, several of the most pleiotropic miRNAs and mRNAs identified in our analyses have been previously related to individual CM traits [26].

Circulating mRNAs associated with multiple metabolic traits

Genes with the greatest pleiotropy across multiple CM traits (i.e. associated with all six CM traits at FDR <0.001) included FAM13A (family with sequence similarity 13, member A), CSF2RB (colony stimulating factor 2 receptor, beta), HIST1H2AC (histone cluster 1, H2ac), and WNK1 (WNK lysine deficient protein kinase 1). FAM13A encodes a Rho GTPase activating protein involved in signal transduction. Variants in FAM13A, as well as PARK2 and RGS6, have been associated with chronic lung disease in prior genome-wide association studies (GWAS) [27, 28]. In another recent GWAS involving over 180,000 participants, a SNP intronic to FAM13A was reported to be associated with HDL-C (p = 4 × 10−12) [22].

CSF2RB encodes a protein that is a common subunit to the three type I cytokine receptors (granulocyte-macrophage colony stimulating factor receptor, as well as the interleukin-3 and interleukin-5 receptors). As it has long been known that inflammatory cells, including macrophages, play important roles in the pathophysiology of atherogenesis, it is perhaps not surprising that recent animal work has demonstrated that CSF2RB expression affects monocyte and macrophage number and function in atherosclerotic lesions [23].

HIST1H2AC (histone cluster 1, H2ac) encodes Histone H2A type IC, one of the four core histones responsible for nucleosome structure in eukaryotic cells. Histones and other DNA-modifying/chromatin remodeling proteins play important roles as mediators of age-related DNA change and have been associated with cardiovascular risk factors, e.g., T2DM, and diseases, including atherosclerosis, myocardial infarction, and heart failure [29]. A GWAS involving 17,000 participants identified SNPs intronic to HIS1H2AC (including rs806971) that were associated with type 1 diabetes mellitus (p = 1.2 × 10−10) [30].

The WNK1 (WNK lysine deficient protein kinase 1) protein is a serine/threonine protein kinase that plays a role in angiogenesis associated with VEGF signaling. Overexpression of WNK1 has also been linked to hypertension and hyperkalemia through alterations in sodium and potassium handling [4]. GWAS have linked two SNPs intronic to, or near, WNK1 with TG levels as well as stroke risk [24].

Another gene exhibiting significant pleiotropy (5 traits; Table 2), ATP-binding cassette G1 (ABCG1), has been associated with total cholesterol levels in GWAS [24]. Lymphocyte-specific protein 1 (LSP1), which was associated with all metabolic traits except glucose (Table 2) harbors variants that are associated with both SBP and DBP [4].

Circulating miRNAs associated with multiple cardiometabolic traits

miRNAs miR-505-5p, miR-197-3p, miR-145-5p, and miR-328 exhibited significant associations with BMI, SBP, DBP, and TG. miR-505-5p targets SRSF1 (Serine/arginine-rich splicing factor 1). SRSF1, in turn, regulates endoglin, vascular endothelial growth factor A, and tissue factor, and controls a molecular senescence program in endothelial cells, leading to age-dependent vascular pathologies [31]. As shown in Table 3, the top gene target of miR-505 is major histocompatibility complex, class I (MR1). Recently, a novel susceptibility locus in MR1 has been associated with coronary disease, likely as a result of dysregulated endothelial function and atherogenesis [31]. As with miR-505-5p, miR-145-5p is highly expressed in smooth muscle cells and controls smooth muscle cell differentiation and function, particularly in the context of metabolic syndrome [32]. Circulating miR-145 deletion results in impaired vascular contractility and differential expression of miR-145 in peripheral blood mononuclear cells has been shown to relate to hypertension [33]. Both miRNAs have strong signals that they are implicated in vascular function and susceptibility to cardiovascular disease.

In another recent study, circulating levels of miR-197 were associated with dyslipidemia in participants with metabolic syndrome and miR-197 levels correlation tightly to body mass index (p = 0.029) [34]. In contrast to the other pleiotropic miRNAs, miR-328 is highly expressed in platelets and has been associated with atrial fibrillation and cardiac hypertrophy, likely through SERCA2A dependent signaling pathways [35]. Notably, variants in the primary gene target of miR-328, transcription factor 7-like 2 gene (TCF7L2) [36], have been associated with T2DM in several genome wide association studies [37]. Prior associations and mRNA targets suggest that miRNAs 197 and 328 may play important roles in the regulation of gene networks influencing body mass and susceptibility to T2DM.

miRNA-miRNA coexpression analyses and cardiometabolic traits

To better understand the molecular mechanisms underlying relations among circulating miRNAs, mRNAs, and CM traits, we conducted miRNA-miRNA coexpression analyses. Figure 3 displays the network of genes coexpressed with the most highly pleiotropic miRNAs. Seventeen mRNAs, including WNK1, that were coexpressed with the four most highly pleiotropic miRNAs, were each associated with multiple CM traits, and were found to be miR targets. We did not otherwise see extensive overlap between miRNA-mRNA targets and miRNA-mRNA coexpression. This may relate to the fact that miRNAs influence many non-traditional mRNA targets, and the effect of a miRNA on a single protein-coding gene target may be too small to be detected. Our findings do suggest, however, that key miRNAs (“hub” miRNAs, such as miR-197-3p) are shared across CM traits and co-express with genes previously associated with CM traits, e.g. WNK1 [38]. We also identified notable associations between miR-145-5p with ARHGEF12 (or LARG), a gene encoding Rho guanine nucleotide exchange factor 12, a vascular smooth muscle signaling protein that is required to develop salt-induced hypertension [39, 40]. Another notable finding was the association between miR-197-3p with SLC4A1, which encodes Band 3 anion transport protein. SLC4A1 has been associated in GWAS with hypertension and its expression in the kidney is altered in animal models with altered sodium absorption [41, 42]. These findings suggest that miRNA-mRNA coexpression pairs may influence vascular phenotypes.

Gene ontology enrichment analysis (Additional file 1: Table S15) revealed that coexpressed mRNAs for the most pleiotropic miRNAs were enriched for RNA metabolism (miR-505-5p), ubiquitin-dependent protein catabolism (miR-197-3p and miR-328), and chromatin assembly (miR-328). The enriched GO terms (e.g., RNA metabolism, protein catabolism, and chromatin assembly) are all relevant to controlling gene expression and protein levels. For example, several genes involved in chromatin assembly (HIST1H4E, HIST1H4B) were highly co-expressed and are known to be related to several cardiovascular diseases [33].

Although our analyses revealed specific miRNAs and mRNAs associated and coexpressed with multiple CM traits, the global effect of these associations is likely to be complex. While “master regulation” may occur in specific settings, these data suggest that a cluster of gene expression changes is contributing to the many relevant pathways found in complex CM systems. It is well known that an individual miRNA can target multiple genes and each protein-coding gene can be regulated by several miRNAs but this complexity is compounded by the fact that most existing studies are performed with single miRNAs, limiting the interpretation of intricate observations. However, the unbiased approach of this study is a strength as it presents the complex findings as potential starting points for future mechanistic investigation.

Study Limitations Transcriptomic signatures may vary by cell type and patterns of mRNA and miRNA expression are known to differ between cell types. Since CM risk factors may influence white cell lineage differentiation, miRNA and mRNA levels in adults with CM traits may reflect differential leukocyte development. Nevertheless, in contrast to many prior analyses focusing on the circulating transcriptome, all observed associations were adjusted for white blood cell counts. Utilization of whole blood derived RNA for the analyses in our study does not provide detailed information on the specific cellular RNA source. We have previously observed that plasma derived extracellular miRNA and blood miRNA levels are often divergent, suggesting that distinct biological sources of RNA may reflect different biological processes and disease associations [43]. This conclusion is strengthened by our previous observations of both concordance and divergence amongst different blood sources of miRNA [44].

We excluded participants receiving several medications, including statins. Although this may have introduced bias (e.g., less severe CM phenotypes), this bias would likely have biased our results toward the null and does not threaten the validity of our findings. Finally, the FHS participants are largely middle-age adults of European ancestry. Generalizability to other younger individuals or those from other racial groups is uncertain.

Conclusions

We found multiple circulating mRNAs and miRNAs that were associated with individual CM phenotypes and with their clustering. Our work supports the hypothesis that circulating transcriptomic patterns can be identified for CM traits and can be used to identify pathways involved in development and progression of CVD and its risk factors.

Abbreviations

CM:

Cardiometabolic

CV:

coefficient of variation

CVD:

cardiovascular disease

DBP:

diastolic blood pressure

GWAS:

genome-wide association studies

HDL-C:

HDL cholesterol

miRNA:

microRNA

mRNA:

messenger RNA

SBP:

systolic

T2DM:

type II diabetes mellitus

TG:

triglycerides

References

  1. O’Donnell CJ, Nabel EG. Genomics of cardiovascular disease. N Engl J Med. 2011;365(22):2098–109.

    Article  PubMed  Google Scholar 

  2. Ejebe KG, Döring A, Dominiczak AF, Demissie S, Deloukas P, de Geus EJ, de Faire U, Crawford G, Collins FS, Chen YD, Caulfield MJ, Campbell H, Burtt NP, Bonnycastle LL, Boomsma DI, Boekholdt SM, Bergman RN, Barroso I, Bandinelli S, Ballantyne CM, Assimes TL, Quertermous T, Altshuler D, Seielstad M, Wong TY, Tai ES, Feranil AB, Kuzawa CW, Adair LS, Taylor Jr HA, Borecki IB, Gabriel SB, Wilson JG, Holm H, Thorsteinsdottir U, Gudnason V, Krauss RM, Mohlke KL, Ordovas JM, Munroe PB, Kooner JS, Tall AR, Hegele RA, Kastelein JJ, Schadt EE, Rotter JI, Boerwinkle E, Strachan DP, Mooser V, Stefansson K, Reilly MP, Samani NJ, Schunkert H, Cupples LA, Sandhu MS, Ridker PM, Rader DJ, van Duijn CM, Peltonen L, Abecasis GR, Boehnke M, Kathiresan S. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466(7307):707–13.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Morris AP, Voight BF, Teslovich TM, Ferreira T, Segrè AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A, Prokopenko I, Kang HM, Dina C, Esko T, Fraser RM, Kanoni S, Kumar A, Lagou V, Langenberg C, Luan J, Lindgren CM, Müller-Nurasyid M, Pechlivanis S, Rayner NW, Scott LJ, Wiltshire S, Yengo L, Kinnunen L, Rossin EJ, Raychaudhuri S, Johnson AD, Dimas AS, Loos RJ, Vedantam S, Chen H, Florez JC, Fox C, Liu CT, Rybin D, Couper DJ, Kao WH, Li M, Cornelis MC, Kraft P, Sun Q, van Dam RM, Stringham HM, Chines PS, Fischer K, Fontanillas P, Holmen OL, Hunt SE, Jackson AU, Kong A, Lawrence R, Meyer J, Perry JR, Platou CG, Potter S, Rehnberg E, Robertson N, Sivapalaratnam S, Stančáková A, Stirrups K, Thorleifsson G, Tikkanen E, Wood AR, Almgren P, Atalay M, Benediktsson R, Bonnycastle LL, Burtt N, Carey J, Charpentier G, Crenshaw AT, Doney AS, Dorkhan M, Edkins S, Emilsson V, Eury E, Forsen T, Gertow K, Gigante B, Grant GB, Groves CJ, Guiducci C, Herder C, Hreidarsson AB, Hui J, James A, Jonsson A, Rathmann W, Klopp N, Kravic J, Krjutškov K, Langford C, Leander K, Lindholm E, Lobbens S, Männistö S, Mirza G, Mühleisen TW, Musk B, Parkin M, Rallidis L, Saramies J, Sennblad B, Shah S, Sigurðsson G, Silveira A, Steinbach G, Thorand B, Trakalo J, Veglia F, Wennauer R, Winckler W, Zabaneh D, Campbell H, van Duijn C, Uitterlinden AG, Hofman A, Sijbrands E, Abecasis GR, Owen KR, Zeggini E, Trip MD, Forouhi NG, Syvänen AC, Eriksson JG, Peltonen L, Nöthen MM, Balkau B, Palmer CN, Lyssenko V, Tuomi T, Isomaa B, Hunter DJ, Qi L, Wellcome Trust Case Control Consortium; Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) Investigators; Genetic Investigation of ANthropometric Traits (GIANT) Consortium; Asian Genetic Epidemiology Network–Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium, Shuldiner AR, Roden M, Barroso I, Wilsgaard T, Beilby J, Hovingh K, Price JF, Wilson JF, Rauramaa R, Lakka TA, Lind L, Dedoussis G, Njølstad I, Pedersen NL, Khaw KT, Wareham NJ, Keinanen-Kiukaanniemi SM, Saaristo TE, Korpi-Hyövälti E, Saltevo J, Laakso M, Kuusisto J, Metspalu A, Collins FS, Mohlke KL, Bergman RN, Tuomilehto J, Boehm BO, Gieger C, Hveem K, Cauchi S, Froguel P, Baldassarre D, Tremoli E, Humphries SE, Saleheen D, Danesh J, Ingelsson E, Ripatti S, Salomaa V, Erbel R, Jöckel KH, Moebus S, Peters A, Illig T, de Faire U, Hamsten A, Morris AD, Donnelly PJ, Frayling TM, Hattersley AT, Boerwinkle E, Melander O, Kathiresan S, Nilsson PM, Deloukas P, Thorsteinsdottir U, Groop LC, Stefansson K, Hu F, Pankow JS, Dupuis J, Meigs JB, Altshuler D, Boehnke M, McCarthy MI, DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 2012;44(9):981–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. International Consortium for Blood Pressure Genome-Wide Association Studies, Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, Chasman DI, Smith AV, Tobin MD, Verwoert GC, Hwang SJ, Pihur V, Vollenweider P, O'Reilly PF, Amin N, Bragg-Gresham JL, Teumer A, Glazer NL, Launer L, Zhao JH, Aulchenko Y, Heath S, Sõber S, Parsa A, Luan J, Arora P, Dehghan A, Zhang F, Lucas G, Hicks AA, Jackson AU, Peden JF, Tanaka T, Wild SH, Rudan I, Igl W, Milaneschi Y, Parker AN, Fava C, Chambers JC, Fox ER, Kumari M, Go MJ, van der Harst P, Kao WH, Sjögren M, Vinay DG, Alexander M, Tabara Y, Shaw-Hawkins S, Whincup PH, Liu Y, Shi G, Kuusisto J, Tayo B, Seielstad M, Sim X, Nguyen KD, Lehtimäki T, Matullo G, Wu Y, Gaunt TR, Onland-Moret NC, Cooper MN, Platou CG, Org E, Hardy R, Dahgam S, Palmen J, Vitart V, Braund PS, Kuznetsova T, Uiterwaal CS, Adeyemo A, Palmas W, Campbell H, Ludwig B, Tomaszewski M, Tzoulaki I, Palmer ND, CARDIoGRAM consortium; CKDGen Consortium; KidneyGen Consortium; EchoGen consortium; CHARGE-HF consortium, Aspelund T, Garcia M, Chang YP, O'Connell JR, Steinle NI, Grobbee DE, Arking DE, Kardia SL, Morrison AC, Hernandez D, Najjar S, McArdle WL, Hadley D, Brown MJ, Connell JM, Hingorani AD, Day IN, Lawlor DA, Beilby JP, Lawrence RW, Clarke R, Hopewell JC, Ongen H, Dreisbach AW, Li Y, Young JH, Bis JC, Kähönen M, Viikari J, Adair LS, Lee NR, Chen MH, Olden M, Pattaro C, Bolton JA, Köttgen A, Bergmann S, Mooser V, Chaturvedi N, Frayling TM, Islam M, Jafar TH, Erdmann J, Kulkarni SR, Bornstein SR, Grässler J, Groop L, Voight BF, Kettunen J, Howard P, Taylor A, Guarrera S, Ricceri F, Emilsson V, Plump A, Barroso I, Khaw KT, Weder AB, Hunt SC, Sun YV, Bergman RN, Collins FS, Bonnycastle LL, Scott LJ, Stringham HM, Peltonen L, Perola M, Vartiainen E, Brand SM, Staessen JA, Wang TJ, Burton PR, Soler Artigas M, Dong Y, Snieder H, Wang X, Zhu H, Lohman KK, Rudock ME, Heckbert SR, Smith NL, Wiggins KL, Doumatey A, Shriner D, Veldre G, Viigimaa M, Kinra S, Prabhakaran D, Tripathy V, Langefeld CD, Rosengren A, Thelle DS, Corsi AM, Singleton A, Forrester T, Hilton G, McKenzie CA, Salako T, Iwai N, Kita Y, Ogihara T, Ohkubo T, Okamura T, Ueshima H, Umemura S, Eyheramendy S, Meitinger T, Wichmann HE, Cho YS, Kim HL, Lee JY, Scott J, Sehmi JS, Zhang W, Hedblad B, Nilsson P, Smith GD, Wong A, Narisu N, Stančáková A, Raffel LJ, Yao J, Kathiresan S, O'Donnell CJ, Schwartz SM, Ikram MA, Longstreth Jr WT, Mosley TH, Seshadri S, Shrine NR, Wain LV, Morken MA, Swift AJ, Laitinen J, Prokopenko I, Zitting P, Cooper JA, Humphries SE, Danesh J, Rasheed A, Goel A, Hamsten A, Watkins H, Bakker SJ, van Gilst WH, Janipalli CS, Mani KR, Yajnik CS, Hofman A, Mattace-Raso FU, Oostra BA, Demirkan A, Isaacs A, Rivadeneira F, Lakatta EG, Orru M, Scuteri A, Ala-Korpela M, Kangas AJ, Lyytikäinen LP, Soininen P, Tukiainen T, Würtz P, Ong RT, Dörr M, Kroemer HK, Völker U, Völzke H, Galan P, Hercberg S, Lathrop M, Zelenika D, Deloukas P, Mangino M, Spector TD, Zhai G, Meschia JF, Nalls MA, Sharma P, Terzic J, Kumar MV, Denniff M, Zukowska-Szczechowska E, Wagenknecht LE, Fowkes FG, Charchar FJ, Schwarz PE, Hayward C, Guo X, Rotimi C, Bots ML, Brand E, Samani NJ, Polasek O, Talmud PJ, Nyberg F, Kuh D, Laan M, Hveem K, Palmer LJ, van der Schouw YT, Casas JP, Mohlke KL, Vineis P, Raitakari O, Ganesh SK, Wong TY, Tai ES, Cooper RS, Laakso M, Rao DC, Harris TB, Morris RW, Dominiczak AF, Kivimaki M, Marmot MG, Miki T, Saleheen D, Chandak GR, Coresh J, Navis G, Salomaa V, Han BG, Zhu X, Kooner JS, Melander O, Ridker PM, Bandinelli S, Gyllensten UB, Wright AF, Wilson JF, Ferrucci L, Farrall M, Tuomilehto J, Pramstaller PP, Elosua R, Soranzo N, Sijbrands EJ, Altshuler D, Loos RJ, Shuldiner AR, Gieger C, Meneton P, Uitterlinden AG, Wareham NJ, Gudnason V, Rotter JI, Rettig R, Uda M, Strachan DP, Witteman JC, Hartikainen AL, Beckmann JS, Boerwinkle E, Vasan RS, Boehnke M, Larson MG, Järvelin MR, Psaty BM, Abecasis GR, Chakravarti A, Elliott P, van Duijn CM, Newton-Cheh C, Levy D, Caulfield MJ, Johnson T. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011;478(7367):103–9.

    Article  Google Scholar 

  5. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, Lango Allen H, Lindgren CM, Luan J, Mägi R, Randall JC, Vedantam S, Winkler TW, Qi L, Workalemahu T, Heid IM, Steinthorsdottir V, Stringham HM, Weedon MN, Wheeler E, Wood AR, Ferreira T, Weyant RJ, Segrè AV, Estrada K, Liang L, Nemesh J, Park JH, Gustafsson S, Kilpeläinen TO, Yang J, Bouatia-Naji N, Esko T, Feitosa MF, Kutalik Z, Mangino M, Raychaudhuri S, Scherag A, Smith AV, Welch R, Zhao JH, Aben KK, Absher DM, Amin N, Dixon AL, Fisher E, Glazer NL, Goddard ME, Heard-Costa NL, Hoesel V, Hottenga JJ, Johansson A, Johnson T, Ketkar S, Lamina C, Li S, Moffatt MF, Myers RH, Narisu N, Perry JR, Peters MJ, Preuss M, Ripatti S, Rivadeneira F, Sandholt C, Scott LJ, Timpson NJ, Tyrer JP, van Wingerden S, Watanabe RM, White CC, Wiklund F, Barlassina C, Chasman DI, Cooper MN, Jansson JO, Lawrence RW, Pellikka N, Prokopenko I, Shi J, Thiering E, Alavere H, Alibrandi MT, Almgren P, Arnold AM, Aspelund T, Atwood LD, Balkau B, Balmforth AJ, Bennett AJ, Ben-Shlomo Y, Bergman RN, Bergmann S, Biebermann H, Blakemore AI, Boes T, Bonnycastle LL, Bornstein SR, Brown MJ, Buchanan TA, Busonero F, Campbell H, Cappuccio FP, Cavalcanti-Proença C, Chen YD, Chen CM, Chines PS, Clarke R, Coin L, Connell J, Day IN, den Heijer M, Duan J, Ebrahim S, Elliott P, Elosua R, Eiriksdottir G, Erdos MR, Eriksson JG, Facheris MF, Felix SB, Fischer-Posovszky P, Folsom AR, Friedrich N, Freimer NB, Fu M, Gaget S, Gejman PV, Geus EJ, Gieger C, Gjesing AP, Goel A, Goyette P, Grallert H, Grässler J, Greenawalt DM, Groves CJ, Gudnason V, Guiducci C, Hartikainen AL, Hassanali N, Hall AS, Havulinna AS, Hayward C, Heath AC, Hengstenberg C, Hicks AA, Hinney A, Hofman A, Homuth G, Hui J, Igl W, Iribarren C, Isomaa B, Jacobs KB, Jarick I, Jewell E, John U, Jørgensen T, Jousilahti P, Jula A, Kaakinen M, Kajantie E, Kaplan LM, Kathiresan S, Kettunen J, Kinnunen L, Knowles JW, Kolcic I, König IR, Koskinen S, Kovacs P, Kuusisto J, Kraft P, Kvaløy K, Laitinen J, Lantieri O, Lanzani C, Launer LJ, Lecoeur C, Lehtimäki T, Lettre G, Liu J, Lokki ML, Lorentzon M, Luben RN, Ludwig B, MAGIC, Manunta P, Marek D, Marre M, Martin NG, McArdle WL, McCarthy A, McKnight B, Meitinger T, Melander O, Meyre D, Midthjell K, Montgomery GW, Morken MA, Morris AP, Mulic R, Ngwa JS, Nelis M, Neville MJ, Nyholt DR, O'Donnell CJ, O'Rahilly S, Ong KK, Oostra B, Paré G, Parker AN, Perola M, Pichler I, Pietiläinen KH, Platou CG, Polasek O, Pouta A, Rafelt S, Raitakari O, Rayner NW, Ridderstråle M, Rief W, Ruokonen A, Robertson NR, Rzehak P, Salomaa V, Sanders AR, Sandhu MS, Sanna S, Saramies J, Savolainen MJ, Scherag S, Schipf S, Schreiber S, Schunkert H, Silander K, Sinisalo J, Siscovick DS, Smit JH, Soranzo N, Sovio U, Stephens J, Surakka I, Swift AJ, Tammesoo ML, Tardif JC, Teder-Laving M, Teslovich TM, Thompson JR, Thomson B, Tönjes A, Tuomi T, van Meurs JB, van Ommen GJ, Vatin V, Viikari J, Visvikis-Siest S, Vitart V, Vogel CI, Voight BF, Waite LL, Wallaschofski H, Walters GB, Widen E, Wiegand S, Wild SH, Willemsen G, Witte DR, Witteman JC, Xu J, Zhang Q, Zgaga L, Ziegler A, Zitting P, Beilby JP, Farooqi IS, Hebebrand J, Huikuri HV, James AL, Kähönen M, Levinson DF, Macciardi F, Nieminen MS, Ohlsson C, Palmer LJ, Ridker PM, Stumvoll M, Beckmann JS, Boeing H, Boerwinkle E, Boomsma DI, Caulfield MJ, Chanock SJ, Collins FS, Cupples LA, Smith GD, Erdmann J, Froguel P, Grönberg H, Gyllensten U, Hall P, Hansen T, Harris TB, Hattersley AT, Hayes RB, Heinrich J, Hu FB, Hveem K, Illig T, Jarvelin MR, Kaprio J, Karpe F, Khaw KT, Kiemeney LA, Krude H, Laakso M, Lawlor DA, Metspalu A, Munroe PB, Ouwehand WH, Pedersen O, Penninx BW, Peters A, Pramstaller PP, Quertermous T, Reinehr T, Rissanen A, Rudan I, Samani NJ, Schwarz PE, Shuldiner AR, Spector TD, Tuomilehto J, Uda M, Uitterlinden A, Valle TT, Wabitsch M, Waeber G, Wareham NJ, Watkins H, Procardis C, Wilson JF, Wright AF, Zillikens MC, Chatterjee N, McCarroll SA, Purcell S, Schadt EE, Visscher PM, Assimes TL, Borecki IB, Deloukas P, Fox CS, Groop LC, Haritunians T, Hunter DJ, Kaplan RC, Mohlke KL, O'Connell JR, Peltonen L, Schlessinger D, Strachan DP, van Duijn CM, Wichmann HE, Frayling TM, Thorsteinsdottir U, Abecasis GR, Barroso I, Boehnke M, Stefansson K, North KE, McCarthy MI, Hirschhorn JN, Ingelsson E, Loos RJ. Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. 2010;42(11):937–48.

    Article  CAS  PubMed Central  Google Scholar 

  6. Huan T, Zhang B, Wang Z, Joehanes R, Zhu J, Johnson AD, Ying S, Munson PJ, Raghavachari N, Wang R, Liu P, Courchesne P, Hwang SJ, Assimes TL, McPherson R, Samani NJ, Schunkert H, Coronary ARteryDIsease Genome wide Replication and Meta-analysis (CARDIoGRAM) Consortium, International Consortium for Blood Pressure GWAS (ICBP), Meng Q, Suver C, O'Donnell CJ, Derry J, Yang X, Levy D. A systems biology framework identifies molecular underpinnings of coronary heart disease. Arterioscler Thromb Vasc Biol. 2013;33(6):1427–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Huan T, Joehanes R, Schurmann C, Schramm K, Pilling LC, Peters MJ, Mägi R, DeMeo D, O'Connor GT, Ferrucci L, Teumer A, Homuth G, Biffar R, Völker U, Herder C, Waldenberger M, Peters A, Zeilinger S, Metspalu A, Hofman A, Uitterlinden AG, Hernandez DG, Singleton AB, Bandinelli S, Munson PJ, Lin H, Benjamin EJ, Esko T, Grabe HJ, Prokisch H, van Meurs JB, Melzer D, Levy D. Whole-Blood Transcriptome Meta-Analysis Identifies Gene Expression Signatures of Cigarette Smoking. Hum Mol Genet. 2016; doi:10.1093/hmg/ddw288

  8. Thum T, Galuppo P, Wolf C, Fiedler J, Kneitz S, van Laake LW, Doevendans PA, Mummery CL, Borlak J, Haverich A, Gross C, Engelhardt S, Ertl G, Bauersachs J. MicroRNAs in the human heart: a clue to fetal gene reprogramming in heart failure. Circulation. 2007;116(3):258–67.

    Article  CAS  PubMed  Google Scholar 

  9. Small EM, Frost RJA, Olson EN. MicroRNAs add a new dimension to cardiovascular disease. Circulation. 2010;121(8):1022–32.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Emilsson V, Thorleifsson G, Zhang B, Leonardson AS, Zink F, Zhu J, Carlson S, Helgason A, Walters GB, Gunnarsdottir S, Mouy M, Steinthorsdottir V, Eiriksdottir GH, Bjornsdottir G, Reynisdottir I, Gudbjartsson D, Helgadottir A, Jonasdottir A, Jonasdottir A, Styrkarsdottir U, Gretarsdottir S, Magnusson KP, Stefansson H, Fossdal R, Kristjansson K, Gislason HG, Stefansson T, Leifsson BG, Thorsteinsdottir U, Lamb JR, Gulcher JR, Reitman ML, Kong A, Schadt EE, Stefansson K. Genetics of gene expression and its effect on disease. Nature. 2008;452(7186):423–8.

    Article  CAS  PubMed  Google Scholar 

  11. Huan T, Esko T, Peters MJ, Pilling LC, Schramm K, Schurmann C, Chen BH, Liu C, Joehanes R, Johnson AD, Yao C, Ying SX, Courchesne P, Milani L, Raghavachari N, Wang R, Liu P, Reinmaa E, Dehghan A, Hofman A, Uitterlinden AG, Hernandez DG, Bandinelli S, Singleton A, Melzer D, Metspalu A, Carstensen M, Grallert H, Herder C, Meitinger T, Peters A, Roden M, Waldenberger M, Dörr M, Felix SB, Zeller T, International Consortium for Blood Pressure GWAS (ICBP), Vasan R, O'Donnell CJ, Munson PJ, Yang X, Prokisch H, Völker U, van Meurs JB, Ferrucci L, Levy D. A meta-analysis of gene expression signatures of blood pressure and hypertension. PLoS Genet. 2015;11(3):e1005035.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Kannel WB, Feinleib M, McNamara PM, Garrison RJ, Castelli WP. An investigation of coronary heart disease in families. The Framingham offspring study. Am J Epidemiol. 1979;110(3):281–90.

    CAS  PubMed  Google Scholar 

  13. Splansky GL, Corey D, Yang Q, Atwood LD, Cupples LA, Benjamin EJ, D'Agostino Sr RB, Fox CS, Larson MG, Murabito JM, O'Donnell CJ, Vasan RS, Wolf PA, Levy D. The Third Generation Cohort of the National Heart, Lung, and Blood Institute’s Framingham Heart Study: design, recruitment, and initial examination. Am J Epidemiol. 2007;165(11):1328–35.

    Article  PubMed  Google Scholar 

  14. McManus DD, Lin H, Tanriverdi K, Quercio M, Yin X, Larson MG, Ellinor PT, Levy D, Freedman JE, Benjamin EJ. Relations between circulating microRNAs and atrial fibrillation: data from the Framingham Offspring Study. Heart Rhythm. 2014;11(4):663–9.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Liu E, McKeown NM, Newby PK, Meigs JB, Vasan RS, Quatromoni PA, D'Agostino RB, Jacques PF. Cross-sectional association of dietary patterns with insulin-resistant phenotypes among adults without diabetes in the Framingham Offspring Study. Br J Nutr. 2009;102(4):576–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ingelsson E, Massaro JM, Sutherland P, Jacques PF, Levy D, D'Agostino RB, Vasan RS, Robins SJ. Contemporary trends in dyslipidemia in the Framingham Heart Study. Arch Intern Med. 2009;169(3):279–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Katz S, Irizarry RA, Lin X, Tripputi M, Porter MW. A summarization approach for Affymetrix GeneChip data using a reference training set from a large, biologically diverse database. BMC Bioinforma. 2006;7:464.

    Article  Google Scholar 

  18. Joehanes R, Ying S, Huan T, Johnson AD, Raghavachari N, Wang R, Liu P, Woodhouse KA, Sen SK, Tanriverdi K, Courchesne P, Freedman JE, O'Donnell CJ, Levy D, Munson PJ. Gene expression signatures of coronary heart disease. Arterioscler Thromb Vasc Biol. 2013;33(6):1418–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Twisk J, Rijmen F. Longitudinal tobit regression: a new approach to analyze outcome variables with floor or ceiling effects. J Clin Epidemiol. 2009;62(9):953–8.

    Article  PubMed  Google Scholar 

  20. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995;57(1):289–300.

    Google Scholar 

  21. Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Global Lipids Genetics Consortium, Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, Beckmann JS, Bragg-Gresham JL, Chang HY, Demirkan A, Den Hertog HM, Do R, Donnelly LA, Ehret GB, Esko T, Feitosa MF, Ferreira T, Fischer K, Fontanillas P, Fraser RM, Freitag DF, Gurdasani D, Heikkilä K, Hyppönen E, Isaacs A, Jackson AU, Johansson A, Johnson T, Kaakinen M, Kettunen J, Kleber ME, Li X, Luan J, Lyytikäinen LP, Magnusson PK, Mangino M, Mihailov E, Montasser ME, Müller-Nurasyid M, Nolte IM, O'Connell JR, Palmer CD, Perola M, Petersen AK, Sanna S, Saxena R, Service SK, Shah S, Shungin D, Sidore C, Song C, Strawbridge RJ, Surakka I, Tanaka T, Teslovich TM, Thorleifsson G, Van den Herik EG, Voight BF, Volcik KA, Waite LL, Wong A, Wu Y, Zhang W, Absher D, Asiki G, Barroso I, Been LF, Bolton JL, Bonnycastle LL, Brambilla P, Burnett MS, Cesana G, Dimitriou M, Doney AS, Döring A, Elliott P, Epstein SE, Eyjolfsson GI, Gigante B, Goodarzi MO, Grallert H, Gravito ML, Groves CJ, Hallmans G, Hartikainen AL, Hayward C, Hernandez D, Hicks AA, Holm H, Hung YJ, Illig T, Jones MR, Kaleebu P, Kastelein JJ, Khaw KT, Kim E, Klopp N, Komulainen P, Kumari M, Langenberg C, Lehtimäki T, Lin SY, Lindström J, Loos RJ, Mach F, McArdle WL, Meisinger C, Mitchell BD, Müller G, Nagaraja R, Narisu N, Nieminen TV, Nsubuga RN, Olafsson I, Ong KK, Palotie A, Papamarkou T, Pomilla C, Pouta A, Rader DJ, Reilly MP, Ridker PM, Rivadeneira F, Rudan I, Ruokonen A, Samani N, Scharnagl H, Seeley J, Silander K, Stancáková A, Stirrups K, Swift AJ, Tiret L, Uitterlinden AG, van Pelt LJ, Vedantam S, Wainwright N, Wijmenga C, Wild SH, Willemsen G, Wilsgaard T, Wilson JF, Young EH, Zhao JH, Adair LS, Arveiler D, Assimes TL, Bandinelli S, Bennett F, Bochud M, Boehm BO, Boomsma DI, Borecki IB, Bornstein SR, Bovet P, Burnier M, Campbell H, Chakravarti A, Chambers JC, Chen YD, Collins FS, Cooper RS, Danesh J, Dedoussis G, de Faire U, Feranil AB, Ferrières J, Ferrucci L, Freimer NB, Gieger C, Groop LC, Gudnason V, Gyllensten U, Hamsten A, Harris TB, Hingorani A, Hirschhorn JN, Hofman A, Hovingh GK, Hsiung CA, Humphries SE, Hunt SC, Hveem K, Iribarren C, Järvelin MR, Jula A, Kähönen M, Kaprio J, Kesäniemi A, Kivimaki M, Kooner JS, Koudstaal PJ, Krauss RM, Kuh D, Kuusisto J, Kyvik KO, Laakso M, Lakka TA, Lind L, Lindgren CM, Martin NG, März W, McCarthy MI, McKenzie CA, Meneton P, Metspalu A, Moilanen L, Morris AD, Munroe PB, Njølstad I, Pedersen NL, Power C, Pramstaller PP, Price JF, Psaty BM, Quertermous T, Rauramaa R, Saleheen D, Salomaa V, Sanghera DK, Saramies J, Schwarz PE, Sheu WH, Shuldiner AR, Siegbahn A, Spector TD, Stefansson K, Strachan DP, Tayo BO, Tremoli E, Tuomilehto J, Uusitupa M, van Duijn CM, Vollenweider P, Wallentin L, Wareham NJ, Whitfield JB, Wolffenbuttel BH, Ordovas JM, Boerwinkle E, Palmer CN, Thorsteinsdottir U, Chasman DI, Rotter JI, Franks PW, Ripatti S, Cupples LA, Sandhu MS, Rich SS, Boehnke M, Deloukas P, Kathiresan S, Mohlke KL, Ingelsson E, Abecasis GR. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45(11):1274–83.

    Article  Google Scholar 

  23. Wang M, Subramanian M, Abramowicz S, Murphy AJ, Gonen A, Witztum J, Welch C, Tabas I, Westerterp M, Tall AR. Interleukin-3/granulocyte macrophage colony-stimulating factor receptor promotes stem cell expansion, monocytosis, and atheroma macrophage burden in mice with hematopoietic ApoE deficiency. Arterioscler Thromb Vasc Biol. 2014;34(5):976–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ikram MA, Seshadri S, Bis JC, Murphy AJ, Gonen A, Witztum J, Welch C, Tabas I, Westerterp M, Tall AR. Genomewide association studies of stroke. N Engl J Med. 2009;360(17):1718–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Raitoharju E, Seppälä I, Oksala N, Lyytikäinen LP, Raitakari O, Viikari J, Ala-Korpela M, Soininen P, Kangas AJ, Waldenberger M, Klopp N, Illig T, Leiviskä J, Loo BM, Hutri-Kähönen N, Kähönen M, Laaksonen R, Lehtimäki T. Blood microRNA profile associates with the levels of serum lipids and metabolites associated with glucose metabolism and insulin resistance and pinpoints pathways underlying metabolic syndrome: the cardiovascular risk in Young Finns Study. Mol Cell Endocrinol. 2014;391(1–2):41–9.

    Article  CAS  PubMed  Google Scholar 

  26. Santovito D, De Nardis V, Marcantonio P, et al. Plasma exosome microRNA profiling unravels a new potential modulator of adiponectin pathway in diabetes: effect of glycemic control. J Clin Endocrinol Metab. 2014;99(9):E1681–5.

    Article  CAS  PubMed  Google Scholar 

  27. Cho MH, Boutaoui N, Klanderman BJ, Mandolini C, Paganelli C, Vitale E, Buttitta F, Bucci M, Mezzetti A, Consoli A, Cipollone F. Variants in FAM13A are associated with chronic obstructive pulmonary disease. Nat Genet. 2010;42(3):200–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lee S, She J, Deng B, Kim J, de Andrade M, Na J, Sun Z, Wampfler JA, Cunningham JM, Wu Y, Limper AH, Aubry MC, Wendt C, Biterman P, Yang P, Lou Z. Multiple-level validation identifies PARK2 in the development of lung cancer and chronic obstructive pulmonary disease. Oncotarget. 2016;7(28):44211–23.

    PubMed  PubMed Central  Google Scholar 

  29. Schiano C, Vietri MT, Grimaldi V, Picascia A, De Pascale MR, Napoli C. Epigenetic-related therapeutic challenges in cardiovascular disease. Trends Pharmacol Sci. 2015;36(4):226–35.

    Article  CAS  PubMed  Google Scholar 

  30. Stankov K, Benc D, Draskovic D. Genetic and epigenetic factors in etiology of diabetes mellitus type 1. Pediatrics. 2013;132(6):1112–22.

    Article  PubMed  Google Scholar 

  31. Davies RW, Wells GA, Stewart AFR, Erdmann J, Shah SH, Ferguson JF, Hall AS, Anand SS, Burnett MS, Epstein SE, Dandona S, Chen L, Nahrstaedt J, Loley C, König IR, Kraus WE, Granger CB, Engert JC, Hengstenberg C, Wichmann HE, Schreiber S, Tang WH, Ellis SG, Rader DJ, Hazen SL, Reilly MP, Samani NJ, Schunkert H, Roberts R, McPherson R. A genome-wide association study for coronary artery disease identifies a novel susceptibility locus in the major histocompatibility complex. Circ Cardiovasc Genet. 2012;5(2):217–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Riches K, Alshanwani AR, Warburton P, O'Regan DJ, Ball SG, Wood IC, Turner NA, Porter KE. Elevated expression levels of miR-143/5 in saphenous vein smooth muscle cells from patients with Type 2 diabetes drive persistent changes in phenotype and function. J Mol Cell Cardiol. 2014;74:240–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Kontaraki JE, Marketou ME, Zacharis EA, Parthenakis FI, Vardas PE. Differential expression of vascular smooth muscle-modulating microRNAs in human peripheral blood mononuclear cells: novel targets in essential hypertension. J Hum Hypertens. 2014;28(8):510–6.

    Article  CAS  PubMed  Google Scholar 

  34. Karolina DS, Tavintharan S, Armugam A, Sepramaniam S, Pek SL, Wong MT, Lim SC, Sum CF, Jeyaseelan K. Circulating miRNA profiles in patients with metabolic syndrome. J Clin Endocrinol Metab. 2012;97(12):E2271–6.

    Article  CAS  PubMed  Google Scholar 

  35. Li C, Li X, Gao X, Zhang R, Zhang Y, Liang H, Xu C, Du W, Zhang Y, Liu X, Ma N, Xu Z, Wang L, Chen X, Lu Y, Ju J, Yang B, Shan H. MicroRNA-328 as a regulator of cardiac hypertrophy. Int J Cardiol. 2014;173(2):268–76. doi:10.1016/j.ijcard.2014.02.035.

    Article  PubMed  Google Scholar 

  36. Chang Y-C, Chang T-J, Jiang Y-D, Kuo SS, Lee KC, Chiu KC, Chuang LM. Association study of the genetic polymorphisms of the transcription factor 7-like 2 (TCF7L2) gene and type 2 diabetes in the Chinese population. Diabetes. 2007;56(10):2631–7.

    Article  CAS  PubMed  Google Scholar 

  37. Marzi C, Huth C, Kolz M, Grallert H, Meisinger C, Wichmann HE, Rathmann W, Herder C, Illig T. Variants of the transcription factor 7-like 2 gene (TCF7L2) are strongly associated with type 2 diabetes but not with the metabolic syndrome in the MONICA/KORA surveys. Horm Metab Res. 2007;39(1):46–52.

    Article  CAS  PubMed  Google Scholar 

  38. Wilson FH, Disse-Nicodème S, Choate KA, Ishikawa K, Nelson-Williams C, Desitter I, Gunel M, Milford DV, Lipkin GW, Achard JM, Feely MP, Dussol B, Berland Y, Unwin RJ, Mayan H, Simon DB, Farfel Z, Jeunemaitre X, Lifton RP. Human hypertension caused by mutations in WNK kinases. Science. 2001;293(5532):1107–12.

    Article  CAS  PubMed  Google Scholar 

  39. Chiu WC, Juang JM, Chang SN, Wu CK, Tsai CT, Tseng CD, Tseng YZ, Su MJ, Chiang FT. Differential baseline expression and angiotensin II-stimulation of leukemia-associated RhoGEF in vascular smooth muscle cells of spontaneously hypertensive rats. Int J Nanomedicine. 2012;7:5929–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Wirth A, Benyó Z, Lukasova M, Leutgeb B, Wettschureck N, Gorbey S, Orsy P, Horváth B, Maser-Gluth C, Greiner E, Lemmer B, Schütz G, Gutkind JS. Offermanns pp. G12–G13-LARG-mediated signaling in vascular smooth muscle is required for salt-induced hypertension. Nat Med. 2008;14(1):64–8.

    Article  CAS  PubMed  Google Scholar 

  41. Schultheis PJ, Clarke LL, Meneton P, Miller ML, Soleimani M, Gawenis LR, Riddle TM, Duffy JJ, Doetschman T, Wang T, Giebisch G, Aronson PS, Lorenz JN, Shull GE. Renal and intestinal absorptive defects in mice lacking the NHE3 Na+/H+ exchanger. Nat Genet. 1998;19(3):282–5.

    Article  CAS  PubMed  Google Scholar 

  42. Kokubo Y, Tomoike H, Tanaka C, Banno M, Okuda T, Inamoto N, Kamide K, Kawano Y, Miyata T. Association of sixty-one non-synonymous polymorphisms in forty-one hypertension candidate genes with blood pressure variation and hypertension. Hypertens Res Off J Jpn Soc Hypertens. 2006;29(8):611–9.

    Article  CAS  Google Scholar 

  43. Shah R, Tanriverdi K, Levy D, Larson M, Gerstein M, Mick E, Rozowsky J, Kitchen R, Murthy V, Mikalev E, Freedman JE. Discordant Expression of Circulating microRNA from Cellular and Extracellular Sources. PLoS One. 2016;11(4), e0153691.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Ward JA, Esa N, Pidikiti R, Freedman JE, Keaney JF, Tanriverdi K, Vitseva O, Ambros V, Lee R, McManus DD. Circulating Cell and Plasma microRNA Profiles Differ between Non-ST-Segment and ST-Segment-Elevation Myocardial Infarction. Fam Med Med Sci Res. 2013;2(2):108.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

DDM’s time was supported by KL2RR031981, 1R01HL126911-01A1, 1R15HL121761-01A1, and 1UH2TR000921-02 from the National Heart, Lung and Blood Institute of the National Institutes of Health, Bethesda, MA.

Funding

The Framingham Heart Study is funded by National Institutes of Health contract N01-HC-25195. The laboratory work for this investigation was funded by the Division of Intramural Research, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD and by NIH contract N01-HC-25195. The analytical component of this project was funded by the Division of Intramural Research, National Heart, Lung, and Blood Institute, and the Center for Information Technology, National Institutes of Health, Bethesda, MD. DDM’s time was supported by KL2RR031981, 1R01HL126911-01A1, 1R15HL121761-01A1, and 1UH2TR000921-02 from the National Heart, Lung and Blood Institute of the National Institutes of Health, Bethesda, MA.

Availability of data and materials

dbGAP. All data used herein are available online in dbGaP (http://www.ncbi.nlm.nih.gov/gap; accession number phs000007).

Authors’ contributions

KT, PM, and JE carried out the RNA and miRNA profiling. DM drafted the manuscript. JR, TH, ML, SL, RJ, and PM participated in the design of the study and performed the statistical analysis. DL and DM conceived of the study, and participated in its design and coordination and VM, RS helped to draft the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interest.

Ethics approval and consent to participate

All participants provided informed consent for genetic research and ethics/human subjects research approval provided by the Boston University Medical School Institutional Review Board.

Deposition of data

RNA sequences are deposited in dbGAP.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to David D. McManus.

Additional files

Additional file 1:

Supplemental data. (DOCX 179 kb)

Additional file 2:

Supplemental figures. (ZIP 249 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

McManus, D.D., Rong, J., Huan, T. et al. Messenger RNA and MicroRNA transcriptomic signatures of cardiometabolic risk factors. BMC Genomics 18, 139 (2017). https://doi.org/10.1186/s12864-017-3533-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-3533-9

Keywords