Whole blood gene expression and atrial fibrillation: the Framingham Heart Study

Peripheral blood RNA biomarkers for cardiovascular disease from bench to bedside: a position paper from the EU-CardioRNA COST action CA17129

Despite significant advances in the diagnosis and treatment of cardiovascular diseases, recent calls have emphasized the unmet need to improve precision-based approaches in cardiovascular disease. Although some studies provide preliminary evidence of the diagnostic and prognostic potential of circulating coding and non-coding RNAs, the complex RNA biology and lack of standardization have hampered the translation of these markers into clinical practice. In this position paper of the CardioRNA COST action CA17129, we provide recommendations to standardize the RNA development process in order to catalyse efforts to investigate novel RNAs for clinical use. We list the unmet clinical needs in cardiovascular disease, such as the identification of high-risk patients with ischaemic heart disease or heart failure who require more intensive therapies. The advantages and pitfalls of the different sample types, including RNAs from plasma, extracellular vesicles, and whole blood, are discussed in the sample matrix, together with their respective analytical methods. The effect of patient demographics and highly prevalent comorbidities, such as metabolic disorders, on the expression of the candidate RNA is presented and should be reported in biomarker studies. We discuss the statistical and regulatory aspects to translate a candidate RNA from a research use only assay to an in-vitro diagnostic test for clinical use. Optimal planning of this development track is required, with input from the researcher, statistician, industry, and regulatory partners.

Genetic and Metabolic Determinants of Atrial Fibrillation in a General Population Sample: The CHRIS Study

Atrial fibrillation (AF) is a supraventricular arrhythmia deriving from uncoordinated electrical activation with considerable associated morbidity and mortality. To expand the limited understanding of AF biological mechanisms, we performed two screenings, investigating the genetic and metabolic determinants of AF in the Cooperative Health Research in South Tyrol study. We found 110 AF cases out of 10,509 general population individuals. A genome-wide association scan (GWAS) identified two novel loci (p-value < 5 × 10⁻⁸) around SNPs rs745582874, next to gene PBX1, and rs768476991, within gene PCCA, with genotype calling confirmed by Sanger sequencing. Risk alleles at both SNPs were enriched in a family detected through familial aggregation analysis of the phenotype, and both rare alleles co-segregated with AF. The metabolic screening of 175 metabolites, in a subset of individuals, revealed a 41% lower concentration of lysophosphatidylcholine lysoPC a C20:3 in AF cases compared to controls (p-adj = 0.005). The genetic findings, combined with previous evidence, indicate that the two identified GWAS loci may be considered novel genetic rare determinants for AF. Considering additionally the association of lysoPC a C20:3 with AF by metabolic screening, our results demonstrate the valuable contribution of the combined genomic and metabolomic approach in studying AF in large-scale population studies.

New biomarkers from multiomics approaches: improving risk prediction of atrial fibrillation

Atrial fibrillation (AF) is a common cardiac arrhythmia leading to many adverse outcomes and increased mortality. Yet the molecular mechanisms underlying AF remain largely unknown. Recent advances in high-throughput technologies make large-scale molecular profiling possible. In the past decade, multiomics studies of AF have identified a number of potential biomarkers of AF. In this review, we focus on the studies of multiomics profiles with AF risk. We summarize recent advances in the discovery of novel biomarkers for AF through multiomics studies. We also discuss limitations and future directions in risk assessment and discovery of therapeutic targets for AF.

Transport of L-Arginine Related Cardiovascular Risk Markers

L-arginine and its derivatives, asymmetric and symmetric dimethylarginine (ADMA and SDMA) and L-homoarginine, have emerged as cardiovascular biomarkers linked to cardiovascular outcomes and various metabolic and functional pathways such as NO-mediated endothelial function. Cellular uptake and efflux of L-arginine and its derivatives are facilitated by transport proteins. In this respect the cationic amino acid transporters CAT1 and CAT2 (SLC7A1 and SLC7A2) and the system y⁺L amino acid transporters (SLC7A6 and SLC7A7) have been most extensively investigated, so far, but the number of transporters shown to mediate the transport of L-arginine and its derivatives is constantly increasing. In the present review we assess the growing body of evidence regarding the function, expression, and clinical relevance of these transporters and their possible relation to cardiovascular diseases.

Relations between plasma microRNAs, echocardiographic markers of atrial remodeling, and atrial fibrillation: Data from the Framingham Offspring study

BACKGROUND

Circulating microRNAs may reflect or influence pathological cardiac remodeling and contribute to atrial fibrillation (AF).

OBJECTIVE

The purpose of this study was to identify candidate plasma microRNAs that are associated with echocardiographic phenotypes of atrial remodeling, and incident and prevalent AF in a community-based cohort.

METHODS

We analyzed left atrial function index (LAFI) of 1788 Framingham Offspring 8 participants. We quantified expression of 339 plasma microRNAs. We examined associations between microRNA levels with LAFI and prevalent and incident AF. We constructed pathway analysis of microRNAs' predicted gene targets to identify molecular processes involved in adverse atrial remodeling in AF.

RESULTS

The mean age of the participants was 66 ± 9 years, and 54% were women. Five percent of participants had prevalent AF at the initial examination and 9% (n = 157) developed AF over a median 8.6 years of follow-up (IQR 8.1-9.2 years). Plasma microRNAs were associated with LAFI (N = 73, p<0.0001). Six of these plasma microRNAs were significantly associated with incident AF, including 4 also associated with prevalent AF (microRNAs 106b, 26a-5p, 484, 20a-5p). These microRNAs are predicted to regulate genes involved in cardiac hypertrophy, inflammation, and myocardial fibrosis.

CONCLUSIONS

Circulating microRNAs 106b, 26a-5p, 484, 20a-5p are associated with atrial remodeling and AF.

Insight into atrial fibrillation through analysis of the coding transcriptome in humans

Atrial fibrillation is the most common sustained cardiac arrhythmia in humans, and its prevalence continues to increase because of the aging of the world population. Much still needs to be learned about the molecular pathways involved in the development and the persistence of the disease. Analysis of the transcriptome of cardiac tissue has provided valuable insight into diverse aspects of atrial remodeling, in particular concerning electrical remodeling-related to ion channels-and structural remodeling identified by dysregulation of processes linked to inflammation, fibrosis, oxidative stress, and thrombogenesis. The huge amount of data produced by these studies now represents a valuable source for the identification of novel potential therapeutic targets. In addition, the shift from cardiac tissue to peripheral blood as a substrate for transcriptome analysis revealed this strategy as a promising tool for improved diagnosis and therefore better patient care.

Integrative Omics Approach to Identifying Genes Associated With Atrial Fibrillation

Rationale: GWAS (Genome-Wide Association Studies) have identified hundreds of genetic loci associated with atrial fibrillation (AF). However, these loci explain only a small proportion of AF heritability. Objective: To develop an approach to identify additional AF-related genes by integrating multiple omics data. Methods and Results: Three types of omics data were integrated: (1) summary statistics from the AFGen 2017 GWAS; (2) a whole blood EWAS (Epigenome-Wide Association Study) of AF; and (3) a whole blood TWAS (Transcriptome-Wide Association Study) of AF. The variant-level GWAS results were collapsed into gene-level associations using fast set-based association analysis. The CpG-level EWAS results were also collapsed into gene-level associations by an adapted SNP-set Kernel Association Test approach. Both GWAS and EWAS gene-based associations were then meta-analyzed with TWAS using a fixed-effects model weighted by the sample size of each data set. A tissue-specific network was subsequently constructed using the NetWAS (Network-Wide Association Study). The identified genes were then compared with the AFGen 2018 GWAS that contained more than triple the number of AF cases compared with AFGen 2017 GWAS. We observed that the multiomics approach identified many more relevant AF-related genes than using AFGen 2018 GWAS alone (1931 versus 206 genes). Many of these genes are involved in the development and regulation of heart- and muscle-related biological processes. Moreover, the gene set identified by multiomics approach explained much more AF variance than those identified by GWAS alone (10.4% versus 3.5%). Conclusions: We developed a strategy to integrate multiple omics data to identify AF-related genes. Our integrative approach may be useful to improve the power of traditional GWAS, which might be particularly useful for rare traits and diseases with limited sample size.

Whole Blood Gene Expression Associated With Clinical Biological Age

Background

Biologic age may better reflect an individual's rate of aging than chronologic age.

Methods

We conducted a transcriptome-wide association study with biologic age estimated with clinical biomarkers, which included: systolic blood pressure, forced expiratory volume at 1 second (FEV1), total cholesterol, fasting glucose, C-reactive protein, and serum creatinine. We assessed the association between the difference between biologic age and chronologic age (∆age) and gene expression in whole blood measured using the Affymetrix Human Exon 1.0st Array.

Results

Our discovery sample included 2,163 participants from the Framingham Offspring cohort (mean age 67 ± 9 years, 55% women). A total of 481 genes were significantly associated with ∆age (p < 2.8 × 10-6). Among them, 415 genes were validated (p < .05/481 = 1.0 × 10-4) in 2,946 participants from the Framingham Third Generation cohort (mean age 46 ± 9 years, 53% women). Many of the significant genes were involved in the ubiquitin-mediated proteolysis pathway. The replication in 414 Rotterdam Study participants (mean age 59 ± 8, 52% women) found 104 of 415 validated genes reached nominal significance (p < .05).

Conclusion

We identified and validated 415 genes associated with ∆age in a community-based cohort. Future functional characterization of the biologic age-related gene network may identify targets to test for interventions to delay aging in older adults.

Investigation of leukocyte RHO/ROCK gene expressions in patients with non-valvular atrial fibrillation

Atrial fibrillation (AF) is an arrhythmia caused by disorganized electrical activity in the atria, and it is an important cause of mortality and morbidity. There is a limited data about Rho/Rho-kinase (ROCK) pathway contribute to AF development. The aim of the present study was to elucidate leukocyte RHO/ROCK gene expressions in patients with non-valvular AF (NVAF). A total of 37 NVAF patients and 47 age and sex-matched controls were included in this study. mRNA was extracted from leukocytes, and real-time polymerase chain reaction was used for gene expression analysis. A marked increase in ROCK1 and ROCK2 gene expressions in patients with NVAF was observed (P<0.0001). The present study detected significant elevations in RHOBTB2, RND3 (RHOE), RHOC, RHOG, RHOH, RAC3, RHOB, RHOD, RHOV, RHOBTB1, RND2, RND1 and RHOJ gene expressions (P<0.01). However, there were marked decreases in CDC42, RAC2, and RHOQ gene expressions in patients with NVAF. No significant modifications were seen in the other Rho GTPase proteins RHOA, RAC1, RHOF, RHOU and RHOBTB3. To the best of our knowledge, the present study is the first to provide data that gene expression of leukocyte RHO/ROCK may contribute to the NVAF pathogenesis through activated leukocytes, which promotes the immune or inflammatory cascade.

Transcriptome-wide association study of inflammatory biologic age

Chronic low grade inflammation is a fundamental mechanism of aging. We estimated biologic age using nine biomarkers from diverse inflammatory pathways and we hypothesized that genes associated with inflammatory biological age would provide insights into human aging. In Framingham Offspring Study participants at examination 8 (2005 to 2008), we used the Klemera-Doubal method to estimate inflammatory biologic age and we computed the difference (∆Age) between biologic age and chronologic age. Gene expression in whole blood was measured using the Affymetrix Human Exon 1.0 ST Array. We used linear mixed effect models to test associations between inflammatory ∆Age and gene expression (dependent variable) adjusting for age, sex, imputed cell counts, and technical covariates. Our study sample included 2386 participants (mean age 67A±9 years, 55% women). There were 448 genes significantly were associated with inflammatory ∆Age (P<2.8x10^-6), 302 genes were positively associated and 146 genes were negatively associated. Pathway analysis among the identified genes highlighted the NOD-like receptor signaling and ubiquitin mediated proteolysis pathways. In summary, we identified 448 genes that were significantly associated with inflammatory biologic age. Future functional characterization may identify molecular interventions to delay aging and prolong healthspan in older adults.

Whole blood gene expression and white matter Hyperintensities

Background

White matter hyperintensities (WMH) are an important biomarker of cumulative vascular brain injury and have been associated with cognitive decline and an increased risk of dementia, stroke, depression, and gait impairments. The pathogenesis of white matter lesions however, remains uncertain. The characterization of gene expression profiles associated with WMH might help uncover molecular mechanisms underlying WMH.

Methods

We performed a transcriptome-wide association study of gene expression profiles with WMH in 3248 participants from the Framingham Heart Study using the Affymetrix Human Exon 1.0 ST Array.

Results

We identified 13 genes that were significantly associated with WMH (FDR < 0.05) after adjusting for age, sex and blood cell components. Many of these genes are involved in inflammation-related pathways.

Conclusion

Thirteen genes were significantly associated with WMH. Our study confirms the hypothesis that inflammation might be an important factor contributing to white matter lesions. Future work is needed to explore if these gene products might serve as potential therapeutic targets.

Electronic supplementary material

The online version of this article (10.1186/s13024-017-0209-5) contains supplementary material, which is available to authorized users.

Leukocyte TRP channel gene expressions in patients with non-valvular atrial fibrillation

Atrial fibrillation (AF) is the most common arrhythmia in clinical practice and is a major cause of morbidity and mortality. The upregulation of TRP channels is believed to mediate the progression of electrical remodelling and the arrhythmogenesis of the diseased heart. However, there is limited data about the contribution of the TRP channels to development of AF. The aim of this study was to investigate leukocyte TRP channels gene expressions in non-valvular atrial fibrillation (NVAF) patients. The study included 47 NVAF patients and 47 sex and age matched controls. mRNA was extracted from blood samples, and real-time polymerase chain reaction was performed for gene expressions by using a dynamic array system. Low levels of TRP channel expressions in the controls were markedly potentiated in NVAF group. We observed marked increases in MCOLN1 (TRPML1), MCOLN2 (TRPML2), MCOLN3 (TRPML3), TRPA1, TRPM1, TRPM2, TRPM3, TRPM4, TRPM5, TRPM6, TRPM7, TRPM8, TRPC1, TRPC2, TRPC3, TRPC4, TRPC5, TRPC6, TRPC7, TRPV1, TRPV2, TRPV3, TRPV4, TRPV5, TRPV6, and PKD2 (TRPP2) gene expressions in NVAF patients (P < 0.05). However, there was no change in PKD1 (TRPP1) gene expression. This is the first study to provide evidence that elevated gene expressions of TRP channels are associated with the pathogenesis of NVAF.

Gene Expression Profiles for the Identification of Prevalent Atrial Fibrillation

BACKGROUND

Diagnosis of atrial fibrillation (AF) can be difficult, requiring cumbersome investigations. We aimed to determine the association of established whole-blood gene expression scores with prevalent AF and to evaluate their performance for the identification of AF in a SIRS (Steroids in Cardiac Surgery) trial cohort.

METHODS AND RESULTS

Whole-blood, transcriptome-wide gene expression profiling was performed using the Illumina HumanHT-12 Expression BeadChip in 416 participants (65% men) before surgery, including 91 with a diagnosis of AF. An AF gene score (GS) calculated from 7 genes reported to be upregulated in AF and a validated GS for biological age based on 1254 genes related to aging were both independently associated with AF diagnosis before surgery in multivariate logistic regression analyses adjusting for known risk factors (P=0.0006 and P=0.003). Addition of AF and biological age GSs to clinical risk factors led to significant improvement in area under the receiver operating characteristic curve (from 0.77 to 0.80; P=0.03), continuous net reclassification improvement index (P<0.0001), and integrated discrimination improvement index (P=0.0002). When stratifying AF by subtype, AF GS was mainly associated with paroxysmal AF (P=0.003), whereas the biological age GS was mainly associated with permanent AF (P=0.017).

CONCLUSIONS

We validated the existence of a blood gene expression signature for prevalent AF and showed that biological age derived from gene expression is significantly associated with prevalent AF. These findings suggest a potential utility of blood gene expression for the identification of patients with AF, particularly paroxysmal AF. This result could have implications for the prevention and management of cryptogenic stroke.

CLINICAL TRIAL REGISTRATION

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00427388.

Methylome-wide Association Study of Atrial Fibrillation in Framingham Heart Study

Atrial fibrillation (AF) is the most common cardiac arrhythmia, but little is known about the molecular mechanisms associated with AF arrhythmogenesis. DNA methylation is an important epigenetic mechanism that regulates gene expression and downstream biological processes. We hypothesize that DNA methylation might play an important role in the susceptibility to develop AF. A total of 2,639 participants from the Offspring Cohort of Framingham Heart Study were enrolled in the current study. These participants included 183 participants with prevalent AF and 220 with incident AF during up to 9 years follow up. Genome-wide methylation was profiled using the Illumina Infinium HumanMethylation450 BeadChip on blood-derived DNA collected during the eighth examination cycle (2005-2008). Two CpG sites were significantly associated with prevalent AF, and five CpGs were associated with incident AF after correction for multiple testing (FDR < 0.05). Fourteen previously reported genome-wide significant AF-related SNP were each associated with at least one CpG site; the most significant association was rs6490029 at the CUX2 locus and cg10833066 (P = 9.5 × 10^-279). In summary, we performed genome-wide methylation profiling in a community-based cohort and identified seven methylation signatures associated with AF. Our study suggests that DNA methylation might play an important role in AF arrhythmogenesis.

Whole Blood Gene Expression Differentiates between Atrial Fibrillation and Sinus Rhythm after Cardioversion

Background

Treatment to restore sinus rhythm among patients with atrial fibrillation (AF) has limited long-term success rates. Gene expression profiling may provide new insights into AF pathophysiology.

Objective

To identify biomarkers and improve our understanding of AF pathophysiology by comparing whole blood gene expression before and after electrical cardioversion (ECV).

Methods

In 46 patients with persistent AF that underwent ECV, whole blood samples were collected 1–2 hours before and 4 to 6 weeks after successful cardioversion. The paired samples were sent for microarray and plasma biomarker comparison.

Results

Of 13,942 genes tested, expression of SLC25A20 and PDK4 had the strongest associations with AF. Post-cardioversion, SLC25A20 and PDK4 expression decreased by 0.8 (CI 0.7–0.8, p = 2.0x10^-6) and 0.7 (CI 0.6–0.8, p = 3.0x10^-5) fold respectively. Median N-terminal pro B-type natriuretic peptide (NT-proBNP) concentrations decreased from 127.7 pg/mL to 44.9 pg/mL (p = 2.3x10^-13) after cardioversion. AF discrimination models combining NT-proBNP and gene expression (NT-proBNP + SLC25A20 area under the curve = 0.88, NT-proBNP + PDK4 AUC = 0.86) had greater discriminative capacity as compared with NT-proBNP alone (AUC = 0.82). Moreover, a model including NT-proBNP, SLC25A20 and PDK4 significantly improved AF discrimination as compared with other models (AUC = 0.87, Net Reclassification Index >0.56, p<5.8x10^-3). We validated the association between SLC25A20 and PDK4 with AF in an independent sample of 17 patients.

Conclusion

This study demonstrates that SLC25A20, PDK4, and NT-proBNP have incremental utility as biomarkers discriminating AF from sinus rhythm. Elevated SLC25A20 and PDK4 expression during AF indicates an important role for energy metabolism in AF.

Genome-wide association study of new-onset atrial fibrillation after coronary artery bypass grafting surgery

BACKGROUND

Postoperative atrial fibrillation (AF) is a potentially life-threatening complication after coronary artery bypass graft (CABG) surgery. Genetic predisposition may predict risk for developing postoperative AF.

METHODS

Study subjects underwent CABG surgery with cardiopulmonary bypass at Duke University Medical Center. In a discovery cohort of 877 individuals from the Perioperative Genetics and Safety Outcomes Study, we performed a genome-wide association study using a logistic regression model with a covariate adjustment for AF risk index. Single-nucleotide polymorphisms (SNPs) that met a P < 5 × 10(-5) were further tested using a replication dataset of 304 individuals from the CATHeterization GENetics biorepository, followed by meta-analysis. Potential pathways related to postoperative AF were identified through gene enrichment analysis using the top genome-wide association study SNPs (P < 10(-4)).

RESULTS

Nineteen SNPs met the a priori defined discovery threshold for replication, but only 3 met nominal significance (P < .05) in the CATHeterization GENetics group, with only one-rs10504554, in the intronic region in lymphocyte antigen 96 (LY96)-showing the same direction of the effect for postoperative AF (odds ratio [OR] 0.48, 95% CI 0.34-0.68, P = 2.9 × 10(-5) vs OR 0.55, 95% CI 0.31-0.99, P = .046) and strong overall association by meta-analysis (meta-P = 4.0 × 10(-6)). Gene enrichment analysis highlighted the role of LY96 in pathways of biologic relevance to activation and modulation of innate immune responses. Our analysis also showed potential association between LY96 and nuclear factor κ-B interaction and postoperative AF through their relevance to inflammatory signaling pathways.

CONCLUSIONS

In patients undergoing CABG surgery, we found genetic polymorphisms in LY96 associated with decreased risk of postoperative AF.

Assessment of immune status using blood transcriptomics and potential implications for global health

The immune system plays a key role in health maintenance and pathogenesis of a wide range of diseases. Leukocytes that are present in the blood convey valuable information about the status of the immune system. Blood transcriptomics, which consists in profiling blood transcript abundance on genome-wide scales, has gained in popularity over the past several years. Indeed, practicality and simplicity largely makes up for what this approach may lack in terms of cell population-level resolution. An extensive survey of the literature reveals increasingly widespread use across virtually all fields of medicine as well as across a number of different animal species, including model organisms but also animals of economical importance. Dissemination across such a wide range of disciplines holds the promise of adding a new perspective, breadth or context, to the considerable depth afforded by whole genome profiling of blood transcript abundance. Indeed, it is only through such contextualization that a truly global perspective will be gained from the use of systems approaches. Also discussed are opportunities that may arise for the fields of immunology and medicine from using blood transcriptomics as a common denominator for developing interactions and cooperation across fields of research that have traditionally been and largely remain compartmentalized. Finally, an argument is made for building immunology research capacity using blood transcriptomics platforms in low-resource and high-disease burden settings.