Age-associated microRNA expression in human peripheral blood is associated with all-cause mortality and age-related traits

Recent studies provide evidence of correlations of DNA methylation and expression of protein‐coding genes with human aging. The relations of microRNA expression with age and age‐related clinical outcomes have not been characterized thoroughly. We explored associations of age with whole‐blood microRNA expression in 5221 adults and identified 127 microRNAs that were differentially expressed by age at P < 3.3 × 10⁻⁴ (Bonferroni‐corrected). Most microRNAs were underexpressed in older individuals. Integrative analysis of microRNA and mRNA expression revealed changes in age‐associated mRNA expression possibly driven by age‐associated microRNAs in pathways that involve RNA processing, translation, and immune function. We fitted a linear model to predict ‘microRNA age’ that incorporated expression levels of 80 microRNAs. MicroRNA age correlated modestly with predicted age from DNA methylation (r = 0.3) and mRNA expression (r = 0.2), suggesting that microRNA age may complement mRNA and epigenetic age prediction models. We used the difference between microRNA age and chronological age as a biomarker of accelerated aging (Δage) and found that Δage was associated with all‐cause mortality (hazards ratio 1.1 per year difference, P = 4.2 × 10⁻⁵ adjusted for sex and chronological age). Additionally, Δage was associated with coronary heart disease, hypertension, blood pressure, and glucose levels. In conclusion, we constructed a microRNA age prediction model based on whole‐blood microRNA expression profiling. Age‐associated microRNAs and their targets have potential utility to detect accelerated aging and to predict risks for age‐related diseases.

MicroRNA Signature of Cigarette Smoking and Evidence for a Putative Causal Role of MicroRNAs in Smoking-Related Inflammation and Target Organ Damage

BACKGROUND

Cigarette smoking increases risk for multiple diseases. MicroRNAs regulate gene expression and may play a role in smoking-induced target organ damage. We sought to describe a microRNA signature of cigarette smoking and relate it to smoking-associated clinical phenotypes, gene expression, and lung inflammatory signaling.

METHODS AND RESULTS

Expression profiling of 283 microRNAs was conducted on whole blood-derived RNA from 5023 Framingham Heart Study participants (54.0% women; mean age, 55±13 years) using TaqMan assays and high-throughput reverse transcription quantitative polymerase chain reaction. Associations of microRNA expression with smoking status and associations of smoking-related microRNAs with inflammatory biomarkers and pulmonary function were tested with linear mixed effects models. We identified a 6-microRNA signature of smoking. Five of the 6 smoking-related microRNAs were associated with serum levels of C-reactive protein or interleukin-6; miR-1180 was associated with pulmonary function measures at a marginally significant level. Bioinformatic evaluation of smoking-associated genes coexpressed with the microRNA signature of cigarette smoking revealed enrichment for immune-related pathways. Smoking-associated microRNAs altered expression of selected inflammatory mediators in cell culture gain-of-function assays.

CONCLUSIONS

We characterized a novel microRNA signature of cigarette smoking. The top microRNAs were associated with systemic inflammatory markers and reduced pulmonary function, correlated with expression of genes involved in immune function, and were sufficient to modulate inflammatory signaling. Our results highlight smoking-associated microRNAs and are consistent with the hypothesis that smoking-associated microRNAs serve as mediators of smoking-induced inflammation and target organ damage. These findings call for further mechanistic studies to explore the diagnostic and therapeutic use of smoking-related microRNAs.

Small RNA-seq during acute maximal exercise reveal RNAs involved in vascular inflammation and cardiometabolic health: brief report

Exercise improves cardiometabolic and vascular function, although the mechanisms remain unclear. Our objective was to demonstrate the diversity of circulating extracellular RNA (ex-RNA) release during acute exercise in humans and its relevance to exercise-mediated benefits on vascular inflammation. We performed plasma small RNA sequencing in 26 individuals undergoing symptom-limited maximal treadmill exercise, with replication of our top candidate miRNA in a separate cohort of 59 individuals undergoing bicycle ergometry. We found changes in miRNAs and other ex-RNAs with exercise (e.g., Y RNAs and tRNAs) implicated in cardiovascular disease. In two independent cohorts of acute maximal exercise, we identified miR-181b-5p as a key ex-RNA increased in plasma after exercise, with validation in a separate cohort. In a mouse model of acute exercise, we found significant increases in miR-181b-5p expression in skeletal muscle after acute exercise in young (but not older) mice. Previous work revealed a strong role for miR-181b-5p in vascular inflammation in obesity, insulin resistance, sepsis, and cardiovascular disease. We conclude that circulating ex-RNAs were altered in plasma after acute exercise target pathways involved in inflammation, including miR-181b-5p. Further investigation into the role of known (e.g., miRNA) and novel (e.g., Y RNAs) RNAs is warranted to uncover new mechanisms of vascular inflammation on exercise-mediated benefits on health.NEW & NOTEWORTHY How exercise provides benefits to cardiometabolic health remains unclear. We performed RNA sequencing in plasma during exercise to identify the landscape of small noncoding circulating transcriptional changes. Our results suggest a link between inflammation and exercise, providing rich data on circulating noncoding RNAs for future studies by the scientific community.

Peer Review Practices for Evaluating Biomedical Research Grants: A Scientific Statement From the American Heart Association

The biomedical research enterprise depends on the fair and objective peer review of research grants, leading to the distribution of resources through efficient and robust competitive methods. In the United States, federal funding agencies and foundations collectively distribute billions of dollars annually to support biomedical research. For the American Heart Association, a Peer Review Subcommittee is charged with establishing the highest standards for peer review. This scientific statement reviews the current literature on peer review practices, describes the current American Heart Association peer review process and those of other agencies, analyzes the strengths and weaknesses of American Heart Association peer review practices, and recommends best practices for the future.

Cross-sectional relations of whole-blood miRNA expression levels and hand grip strength in a community sample

MicroRNAs (miRNAs) regulate gene expression with emerging data suggesting miRNAs play a role in skeletal muscle biology. We sought to examine the association of miRNAs with grip strength in a community-based sample. Framingham Heart Study Offspring and Generation 3 participants (n = 5668 54% women, mean age 55 years, range 24, 90 years) underwent grip strength measurement and miRNA profiling using whole blood from fasting morning samples. Linear mixed-effects regression modeling of grip strength (kg) versus continuous miRNA 'Cq' values and versus binary miRNA expression was performed. We conducted an integrative miRNA-mRNA coexpression analysis and examined the enrichment of biologic pathways for the top miRNAs associated with grip strength. Grip strength was lower in women than in men and declined with age with a mean 44.7 (10.0) kg in men and 26.5 (6.3) kg in women. Among 299 miRNAs interrogated for association with grip strength, 93 (31%) had FDR q value < 0.05, 54 (18%) had an FDR q value < 0.01, and 15 (5%) had FDR q value < 0.001. For almost all miRNA-grip strength associations, increasing miRNA concentration is associated with increasing grip strength. miR-20a-5p (FDR q 1.8 × 10-6 ) had the most significant association and several among the top 15 miRNAs had links to skeletal muscle including miR-126-3p, miR-30a-5p, and miR-30d-5p. The top associated biologic pathways included metabolism, chemokine signaling, and ubiquitin-mediated proteolysis. Our comprehensive assessment in a community-based sample of miRNAs in blood associated with grip strength provides a framework to further our understanding of the biology of muscle strength.

Ideal Cardiovascular Health, Cardiovascular Remodeling, and Heart Failure in Blacks: The Jackson Heart Study

BACKGROUND

The lifetime risk of heart failure (HF) is higher in the black population than in other racial groups in the United States.

METHODS AND RESULTS

We measured the Life's Simple 7 ideal cardiovascular health metrics in 4195 blacks in the JHS (Jackson Heart Study; 2000-2004). We evaluated the association of Simple 7 metrics with incident HF and left ventricular structure and function by cardiac magnetic resonance (n=1188). Mean age at baseline was 54.4 years (65% women). Relative to 0 to 2 Simple 7 factors, blacks with 3 factors had 47% lower incident HF risk (hazard ratio [HR], 0.53; 95% confidence interval [CI], 0.39-0.73; P<0.0001); and those with ≥4 factors had 61% lower HF risk (HR, 0.39; 95% CI, 0.24-0.64; P=0.0002). Higher blood pressure (HR, 2.32; 95% CI, 1.28-4.20; P=0.005), physical inactivity (HR, 1.65; 95% CI, 1.07-2.55; P=0.02), smoking (HR, 2.04; 95% CI, 1.43-2.91; P<0.0001), and impaired glucose control (HR, 1.76; 95% CI, 1.34-2.29; P<0.0001) were associated with incident HF. The age-/sex-adjusted population attributable risk for these Simple 7 metrics combined was 37.1%. Achievement of ideal blood pressure, ideal body mass index, ideal glucose control, and nonsmoking was associated with less likelihood of adverse cardiac remodeling by cardiac magnetic resonance.

CONCLUSIONS

Cardiovascular risk factors in midlife (specifically elevated blood pressure, physical inactivity, smoking, and poor glucose control) are associated with incident HF in blacks and represent targets for intensified HF prevention.

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.

Extracellular RNAs Are Associated With Insulin Resistance and Metabolic Phenotypes

OBJECTIVE

Insulin resistance (IR) is a hallmark of obesity and metabolic disease. Circulating extracellular RNAs (ex-RNAs), stable RNA molecules in plasma, may play a role in IR, though most studies on ex-RNAs in IR are small. We sought to characterize the relationship between ex-RNAs and metabolic phenotypes in a large community-based human cohort.

RESEARCH DESIGN AND METHODS

We measured circulating plasma ex-RNAs in 2,317 participants without diabetes in the Framingham Heart Study (FHS) Offspring Cohort at cycle 8 and defined associations between ex-RNAs and IR (measured by circulating insulin level). We measured association between candidate ex-RNAs and markers of adiposity. Sensitivity analyses included individuals with diabetes. In a separate cohort of 90 overweight/obese youth, we measured selected ex-RNAs and metabolites. Biology of candidate microRNAs was investigated in silico.

RESULTS

The mean age of FHS participants was 65.8 years (56% female), with average BMI 27.7 kg/m²; participants in the youth cohort had a mean age of 15.5 years (60% female), with mean BMI 33.8 kg/m². In age-, sex-, and BMI-adjusted models across 391 ex-RNAs in FHS, 18 ex-RNAs were associated with IR (of which 16 were microRNAs). miR-122 was associated with IR and regional adiposity in adults and IR in children (independent of metabolites). Pathway analysis revealed metabolic regulatory roles for miR-122, including regulation of IR pathways (AMPK, target of rapamycin signaling, and mitogen-activated protein kinase).

CONCLUSIONS

These results provide translational evidence in support of an important role of ex-RNAs as novel circulating factors implicated in IR.

Thrombosis and platelets: an update

Haemostasis and thrombosis are complex, multifactorial processes. There is an evolving understanding of the mechanisms influencing vascular occlusion and the role of inflammation and immunity. Despite major advances in elucidating the mechanistic pathways mediating platelet function and thrombosis, challenges in the treatment of vascular occlusive diseases persist. Pharmacological advances have greatly affected thrombotic outcomes, but this has led to the unwanted side effect of bleeding. Detailed assessment of the impact of non-thrombotic diseases on haemostasis and thrombosis is necessary to better evaluate thrombotic risk and establish optimal treatment. This review will focus on recent advances in understanding the contribution of evolving risk factors to thrombosis.

Stroke and Circulating Extracellular RNAs

Supplemental Digital Content is available in the text.

Background and Purpose—

There is increasing interest in extracellular RNAs (ex-RNAs), with numerous reports of associations between selected microRNAs (miRNAs) and a variety of cardiovascular disease phenotypes. Previous studies of ex-RNAs in relation to risk for cardiovascular disease have investigated small numbers of patients and assayed only candidate miRNAs. No human studies have investigated links between novel ex-RNAs and stroke.

Methods—

We conducted unbiased next-generation sequencing using plasma from 40 participants of the FHS (Framingham Heart Study; Offspring Cohort Exam 8) followed by high-throughput polymerase chain reaction of 471 ex-RNAs. The reverse transcription quantitative polymerase chain reaction included 331 of the most abundant miRNAs, 43 small nucleolar RNAs, and 97 piwi-interacting RNAs in 2763 additional FHS participants and explored the relations of ex-RNAs and prevalent (n=63) and incident (n=51) stroke and coronary heart disease (prevalent=286, incident=69).

Results—

After adjustment for multiple cardiovascular disease risk factors, 7 ex-RNAs were associated with stroke prevalence or incidence; there were no ex-RNA associated with prevalent or incident coronary heart disease. Statistically significant ex-RNA associations with stroke were specific, with no overlap between prevalent and incident events.

Conclusions—

This is the largest study of ex-RNAs in relation to stroke using an unbiased approach in an observational cohort and the first large study to examine human small noncoding RNAs beyond miRNAs. These results demonstrate that when studied in a large observational cohort, extracellular miRNAs are associated with stroke risk.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3533-9) contains supplementary material, which is available to authorized users.

Integrated genome-wide analysis of expression quantitative trait loci aids interpretation of genomic association studies

Background

Identification of single nucleotide polymorphisms (SNPs) associated with gene expression levels, known as expression quantitative trait loci (eQTLs), may improve understanding of the functional role of phenotype-associated SNPs in genome-wide association studies (GWAS). The small sample sizes of some previous eQTL studies have limited their statistical power. We conducted an eQTL investigation of microarray-based gene and exon expression levels in whole blood in a cohort of 5257 individuals, exceeding the single cohort size of previous studies by more than a factor of 2.

Results

We detected over 19,000 independent lead cis-eQTLs and over 6000 independent lead trans-eQTLs, targeting over 10,000 gene targets (eGenes), with a false discovery rate (FDR) < 5%. Of previously published significant GWAS SNPs, 48% are identified to be significant eQTLs in our study. Some trans-eQTLs point toward novel mechanistic explanations for the association of the SNP with the GWAS-related phenotype. We also identify 59 distinct blocks or clusters of trans-eQTLs, each targeting the expression of sets of six to 229 distinct trans-eGenes. Ten of these sets of target genes are significantly enriched for microRNA targets (FDR < 5%). Many of these clusters are associated in GWAS with multiple phenotypes.

Conclusions

These findings provide insights into the molecular regulatory patterns involved in human physiology and pathophysiology. We illustrate the value of our eQTL database in the context of a recent GWAS meta-analysis of coronary artery disease and provide a list of targeted eGenes for 21 of 58 GWAS loci.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1142-6) contains supplementary material, which is available to authorized users.

Cardiometabolic disease in South Asians: A global health concern in an expanding population

Cardiovascular disease (CVD) is one of the main causes of mortality and morbidity worldwide. As an emerging population, South Asians (SAs) bear a disproportionately high burden of CVD relative to underlying classical risk factors, partly attributable to a greater prevalence of insulin resistance and diabetes and distinct genetic and epigenetic influences. While the phenotypic distinctions between SAs and other ethnicities in CVD risk are becoming increasingly clear, the biology of these conditions remains an area of active investigation, with emerging studies involving metabolism, genetic variation and epigenetic modifiers (e.g., extracellular RNA). In this review, we describe the current literature on prevalence, prognosis and CVD risk in SAs, and provide a landscape of translational research in this field toward ameliorating CVD risk in SAs.

Challenges and Opportunities in Linking Long Noncoding RNAs to Cardiovascular, Lung, and Blood Diseases

The new millennium heralds an unanticipated surge of genomic information, most notably an expansive class of long noncoding RNAs (lncRNAs). These transcripts, which now outnumber all protein-coding genes, often exhibit the same characteristics as mRNAs (RNA polymerase II-dependent, 5' methyl-capped, multiexonic, polyadenylated); yet, they do not encode for stable, well-conserved proteins. Elucidating the function of all relevant lncRNAs in heart, vasculature, lung, and blood is essential for generating a complete interactome in these tissues. This is particularly evident because an increasing number of investigators perform RNA-sequencing experiments where, typically, annotated lncRNAs exhibit impressive changes in gene expression. How does one go about evaluating an lncRNA when the sequence of the transcript lends no insight into how it may function within a cell type? Here, we provide a brief overview for the rational study of lncRNAs.

Discordant Expression of Circulating microRNA from Cellular and Extracellular Sources

MicroRNA (miRNA) expression has rapidly grown into one of the largest fields for disease characterization and development of clinical biomarkers. Consensus is lacking in regards to the optimal sample source or if different circulating sources are concordant. Here, using miRNA measurements from contemporaneously obtained whole blood- and plasma-derived RNA from 2391 individuals, we demonstrate that plasma and blood miRNA levels are divergent and may reflect different biological processes and disease associations.

Comparison of RNA isolation and associated methods for extracellular RNA detection by high-throughput quantitative polymerase chain reaction

MicroRNAs (miRNAs) are small noncoding RNA molecules that function in RNA silencing and posttranscriptional regulation of gene expression. miRNAs in biofluids are being used for clinical diagnosis as well as disease prediction. Efficient and reproducible isolation methods are crucial for extracellular RNA detection. To determine the best methodologies for miRNA detection from plasma, the performance of four RNA extraction kits, including an in-house kit, were determined with miScript miRNA assay technology; all were measured using a high-throughput quantitative polymerase chain reaction (qPCR) platform (BioMark System) with 90 human miRNA assays. In addition, the performances of complementary DNA (cDNA) and preamplification kits for TaqMan miRNA assays and miScript miRNA assays were compared using the same 90 miRNAs on the BioMark System. There were significant quantification cycle (Cq) value differences for the detection of miRNA targets between isolation kits. cDNA, preamplification, and qPCR performances were also varied. In summary, this study demonstrates differences among RNA isolation methods as measured by reverse transcription (RT)-qPCR. Importantly, differences were also noted in cDNA and preamplification performance using TaqMan and miScript. The in-house kit performed better than the other three kits. These findings demonstrate significant variability between isolation and detection methods for low-abundant miRNA detection from biofluids.

Characterization of the platelet transcriptome by RNA sequencing in patients with acute myocardial infarction

Transcripts in platelets are largely produced in precursor megakaryocytes but remain physiologically active as platelets translate RNAs and regulate protein/RNA levels. Recent studies using transcriptome sequencing (RNA-seq) characterized the platelet transcriptome in limited number of non-diseased individuals. Here, we expand upon these RNA-seq studies by completing RNA-seq in platelets from 32 patients with acute myocardial infarction (MI). Our goals were to characterize the platelet transcriptome using a population of patients with acute MI and relate gene expression to platelet aggregation measures and ST-segment elevation MI (STEMI) (n = 16) vs. non-STEMI (NSTEMI) (n = 16) subtypes. Similar to other studies, we detected 9565 expressed transcripts, including several known platelet-enriched markers (e.g. PPBP, OST4). Our RNA-seq data strongly correlated with independently ascertained platelet expression data and showed enrichment for platelet-related pathways (e.g. wound response, hemostasis, and platelet activation), as well as actin-related and post-transcriptional processes. Several transcripts displayed suggestively higher (FBXL4, ECHDC3, KCNE1, TAOK2, AURKB, ERG, and FKBP5) and lower (MIAT, PVRL3, and PZP) expression in STEMI platelets compared to NSTEMI. We also identified transcripts correlated with platelet aggregation to TRAP (ATP6V1G2, SLC2A3), collagen (CEACAM1, ITGA2), and ADP (PDGFB, PDGFC, ST3GAL6). Our study adds to current platelet gene expression resources by providing transcriptome-wide analyses in platelets isolated from patients with acute MI. In concert with prior studies, we identify various genes for further study in regards to platelet function and acute MI. Future platelet RNA-seq studies examining more diverse sets of healthy and diseased samples will add to our understanding of platelet thrombotic and non-thrombotic functions.

Extracellular RNAs: development as biomarkers of human disease

Ten ongoing studies designed to test the possibility that extracellular RNAs may serve as biomarkers in human disease are described. These studies, funded by the NIH Common Fund Extracellular RNA Communication Program, examine diverse extracellular body fluids, including plasma, serum, urine and cerebrospinal fluid. The disorders studied include hepatic and gastric cancer, cardiovascular disease, chronic kidney disease, neurodegenerative disease, brain tumours, intracranial haemorrhage, multiple sclerosis and placental disorders. Progress to date and the plans for future studies are outlined.

The role of circulating platelet transcripts

Our understanding of platelets, anucleate cells with a traditional role in hemostasis and inflammation, has developed greatly over the last decade. Platelets' role in the systemic response of the body to vascular injury, inflammation, and infection has expanded as has our understanding of their importance to the body's regulation of these processes. One recently explored mechanism by which platelets regulate the body's inflammatory and immune response is through its endogenous RNA. Platelets' messenger RNA (mRNAs) and microRNA (miRNAs) profiles have been shown to reflect disease and disease risk factors and have been correlated with select human clinical phenotypes. Developing an understanding of platelet transcripts in the circulation elucidates how platelets function in both their traditional thrombotic role and non-traditional functions and may have widespread implications in several fields including thrombosis, infection, cancer, and systemic inflammation.

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.