Genetic loci associated with heart rate variability and their effects on cardiac disease risk

Publication guidelines for human heart rate and heart rate variability studies in psychophysiology-Part 1: Physiological underpinnings and foundations of measurement

This Committee Report provides methodological, interpretive, and reporting guidance for researchers who use measures of heart rate (HR) and heart rate variability (HRV) in psychophysiological research. We provide brief summaries of best practices in measuring HR and HRV via electrocardiographic and photoplethysmographic signals in laboratory, field (ambulatory), and brain-imaging contexts to address research questions incorporating measures of HR and HRV. The Report emphasizes evidence for the strengths and weaknesses of different recording and derivation methods for measures of HR and HRV. Along with this guidance, the Report reviews what is known about the origin of the heartbeat and its neural control, including factors that produce and influence HRV metrics. The Report concludes with checklists to guide authors in study design and analysis considerations, as well as guidance on the reporting of key methodological details and characteristics of the samples under study. It is expected that rigorous and transparent recording and reporting of HR and HRV measures will strengthen inferences across the many applications of these metrics in psychophysiology. The prior Committee Reports on HR and HRV are several decades old. Since their appearance, technologies for human cardiac and vascular monitoring in laboratory and daily life (i.e., ambulatory) contexts have greatly expanded. This Committee Report was prepared for the Society for Psychophysiological Research to provide updated methodological and interpretive guidance, as well as to summarize best practices for reporting HR and HRV studies in humans.

Genetic and polygenic investigation of heart rate variability to identify biomarkers associated with Anxiety disorders

Given that anxiety disorders (AD) are associated with reduced vagally-mediated heart rate variability (HRV), genetic variants related to HRV may provide insight into anxiety etiology. This study used polygenic risk scores (PRS) to explore the genetic overlap between AD and HRV, and investigated whether HRV-related polymorphisms influence anxiety risk. Resting vagally-mediated HRV was measured using a wearable device in 188 European individuals (AD=101, healthy controls=87). AD PRS was tested for association with resting HRV, and HRV PRS for association with AD. We also investigated 15 significant hits from an HRV genome-wide association study (GWAS) for association with resting HRV and AD and if this association is mediated through resting HRV. The AD PRS and HRV PRS showed nominally significant associations with resting HRV and anxiety disorders, respectively. HRV GWAS variants associated with resting HRV were rs12980262 (NDUFA11), rs2680344 (HCN4), rs4262 and rs180238 (GNG11), and rs10842383 (LINC00477). Mediation analyses revealed that NDUFA11 rs12980262 A-carriers and GNG11 rs180238 and rs4262 C-carriers had higher anxiety risk through lower HRV. This study supports an anxiety-HRV genetic relationship, with HRV-related genetic variants translating to AD. This study encourages exploration of HRV genetics to understand mechanisms and identify novel treatment targets for anxiety.

Exploring the impact of electrocardiographic parameters on the risk of common arrhythmias: a two-sample Mendelian randomization study

Background

Observational studies have shown that heart rate (HR), heart rate variability (HRV), P-wave terminal force, P-wave duration, T-wave amplitude and PR interval are associated with risk factors for atrial fibrillation (AF) or bradycardia. Arrhythmias are associated with many causes of hospitalization. However, observational studies are susceptible to confounding factors that have not yet been identified. The objective of this study was to clarify the causal relationships by Mendelian randomization analysis.

Methods

We conducted a two-sample and multivariate Mendelian randomization (MVMR) analysis using genome-wide association study (GWAS) data from a European population to assess the total and direct causal effects of HR, three HRV traits, P-wave terminal force, P-wave duration, T-wave top amplitude in five-lead modes, and the PR interval on the risk of AF (N=191,205), bradycardia (N=463,010), and supraventricular tachycardia (SVT) (N=463,010).

Results

The results of the univariate MR analysis revealed the following significant causal effects: the higher the genetically predicted PR interval, the lower the risk of AF; the higher the HR and T-wave top amplitude (aVR leads and V3 + V4 + aVL leads), the lower the risk of bradycardia; and the higher HR and the lower PR interval, the higher the risk of SVT. The multivariate MR results indicated that the HRV_standard deviation of the normal-to-normal (SDNN) interval had an independent causal effect on the risk of AF [odds ratio (OR): 0.515; 95% confidence interval (CI): 0.278-0.954; P=0.03], and the T-wave top amplitude in the aVR leads (OR: 0.998; 95% CI: 0.996-0.999; P<0.001) and the HRV_SDNN (OR: 0.988; 95% CI: 0.976-1.000; P=0.045) had independent causal effects on the risk of bradycardia.

Conclusions

The HRV_SDNN had an independent causal effect on AF, while the HRV_SDNN and T-wave top amplitude in the aVR leads had independent causal effects on bradycardia, which suggests that some of the electrocardiographic parameters have preventive effects on the incidence of AF and bradycardia.

Heart Rate Variability, Deceleration Capacity of Heart Rate, and Death: A Veteran Twins Study

BACKGROUND

Autonomic function can be measured noninvasively using heart rate variability (HRV), which indexes overall sympathovagal balance. Deceleration capacity (DC) of heart rate is a more specific metric of vagal modulation. Higher values of these measures have been associated with reduced mortality risk primarily in patients with cardiovascular disease, but their significance in community samples is less clear.

METHODS AND RESULTS

This prospective twin study followed 501 members from the VET (Vietnam Era Twin) registry. At baseline, frequency domain HRV and DC were measured from 24-hour Holter ECGs. During an average 12-year follow-up, all-cause death was assessed via the National Death Index. Multivariable Cox frailty models with random effect for twin pair were used to examine the hazard ratios of death per 1-SD increase in log-transformed autonomic metrics. Both in the overall sample and comparing twins within pairs, higher values of low-frequency HRV and DC were significantly associated with lower hazards of all-cause death. In within-pair analysis, after adjusting for baseline factors, there was a 22% and 27% lower hazard of death per 1-SD increment in low-frequency HRV and DC, respectively. Higher low-frequency HRV and DC, measured during both daytime and nighttime, were associated with decreased hazard of death, but daytime measures showed numerically stronger associations. Results did not substantially vary by zygosity.

CONCLUSIONS

Autonomic inflexibility, and especially vagal withdrawal, are important mechanistic pathways of general mortality risk, independent of familial and genetic factors.

The relationship between emotional disorders and heart rate variability: A Mendelian randomization study

OBJECTIVE

Previous studies have shown that emotional disorders are negatively associated with heart rate variability (HRV), but the potential causal relationship between genetic susceptibility to emotional disorders and HRV remains unclear. We aimed to perform a Mendelian randomization (MR) study to investigate the potential association between emotional disorders and HRV.

METHODS

The data used for this study were obtained from publicly available genome-wide association study datasets. Five models, including the inverse variance weighted model (IVW), the weighted median estimation model (WME), the weighted model-based method (WM), the simple model (SM) and the MR-Egger regression model (MER), were utilized for MR. The leave-one-out sensitivity test, MR pleiotropy residual sum and outlier test (MR-PRESSO) and Cochran's Q test were used to confirm heterogeneity and pleiotropy.

RESULTS

MR analysis revealed that genetic susceptibility to broad depression was negatively correlated with HRV (pvRSA/HF) (OR = 0.380, 95% CI 0.146-0.992; p = 0.048). However, genetic susceptibility to irritability was positively correlated with HRV (pvRSA/HF, SDNN) (OR = 2.017, 95% CI 1.152-3.534, p = 0.008) (OR = 1.154, 95% CI 1.000-1.331, p = 0.044). Genetic susceptibility to anxiety was positively correlated with HRV (RMSSD) (OR = 2.106, 95% CI 1.032-4.299; p = 0.041). No significant directional pleiotropy or heterogeneity was detected. The accuracy and robustness of these findings were confirmed through a sensitivity analysis.

CONCLUSIONS

Our MR study provides genetic support for the causal effects of broad depression, irritable mood, and anxiety on HRV.

Effects of endocrine disrupting chemicals and their interactions with genetic risk scores on cardiometabolic traits

Ubiquitous non-persistent endocrine disrupting chemicals (EDCs) have inconsistent associations with cardiometabolic traits. Additionally, large-scale genome-wide association studies (GWASs) have yielded many genetic risk variants for cardiometabolic traits and diseases. This study aimed to investigate the associations between a wide range of EDC exposures (parabens, bisphenols, and phthalates) and 14 cardiometabolic traits and whether these are moderated by their respective genetic risk scores (GRSs). Data were from 1074 participants aged 18 years or older of the Lifelines Cohort Study, a large population-based biobank. GRSs for 14 cardiometabolic traits were calculated based on genome-wide significant common variants from recent GWASs. The concentrations of 15 EDCs in 24-hour urine were measured by isotope dilution liquid chromatography tandem mass spectrometry technology. The main effects of trait-specific GRSs and each of the EDC exposures and their interaction effects on the 14 cardiometabolic traits were examined in multiple linear regression. The present study confirmed significant main effects for all GRSs on their corresponding cardiometabolic trait. Regarding the main effects of EDC exposures, 26 out of 280 EDC-trait tests were significant with explained variances ranging from 0.43 % (MMP- estimated glomerular filtration rate (eGFR)) to 2.37 % (PrP-waist-hip ratio adjusted body mass index (WHRadjBMI)). We confirmed the association of MiBP and MBzP with WHRadjBMI and body mass index (BMI), and showed that parabens, bisphenol F, and many other phthalate metabolites significantly contributed to the variance of WHRadjBMI, BMI, high-density lipoprotein (HDL), eGFR, fasting glucose (FG), and diastolic blood pressure (DBP). Only one association between BMI and bisphenol F was nominally significantly moderated by the GRS explaining 0.36 % of the variance. However, it did not survive multiple testing correction. We showed that non-persistent EDC exposures exerted effects on BMI, WHRadjBMI, HDL, eGFR, FG, and DBP. However no evidence for a modulating role of GRSs was found.

RGS6 drives cardiomyocyte death following nucleolar stress by suppressing Nucleolin/miRNA-21

BACKGROUND

Prior evidence demonstrated that Regulator of G protein Signaling 6 (RGS6) translocates to the nucleolus in response to cytotoxic stress though the functional significance of this phenomenon remains unknown.

METHODS

Utilizing in vivo gene manipulations in mice, primary murine cardiac cells, human cell lines and human patient samples we dissect the participation of a RGS6-nucleolin complex in chemotherapy-dependent cardiotoxicity.

RESULTS

Here we demonstrate that RGS6 binds to a key nucleolar protein, Nucleolin, and controls its expression and activity in cardiomyocytes. In the human myocyte AC-16 cell line, induced pluripotent stem cell derived cardiomyocytes, primary murine cardiomyocytes, and the intact murine myocardium tuning RGS6 levels via overexpression or knockdown resulted in diametrically opposed impacts on Nucleolin mRNA, protein, and phosphorylation.RGS6 depletion provided marked protection against nucleolar stress-mediated cell death in vitro, and, conversely, RGS6 overexpression suppressed ribosomal RNA production, a key output of the nucleolus, and triggered death of myocytes. Importantly, overexpression of either Nucleolin or Nucleolin effector miRNA-21 counteracted the pro-apoptotic effects of RGS6. In both human and murine heart tissue, exposure to the genotoxic stressor doxorubicin was associated with an increase in the ratio of RGS6/Nucleolin. Preventing RGS6 induction via introduction of RGS6-directed shRNA via intracardiac injection proved cardioprotective in mice and was accompanied by restored Nucleolin/miRNA-21 expression, decreased nucleolar stress, and decreased expression of pro-apoptotic, hypertrophy, and oxidative stress markers in heart.

CONCLUSION

Together, these data implicate RGS6 as a driver of nucleolar stress-dependent cell death in cardiomyocytes via its ability to modulate Nucleolin. This work represents the first demonstration of a functional role for an RGS protein in the nucleolus and identifies the RGS6/Nucleolin interaction as a possible new therapeutic target in the prevention of cardiotoxicity.

Sex-specific genetic architecture of late-life memory performance

BACKGROUND

Women demonstrate a memory advantage when cognitively healthy yet lose this advantage to men in Alzheimer's disease. However, the genetic underpinnings of this sex difference in memory performance remain unclear.

METHODS

We conducted the largest sex-aware genetic study on late-life memory to date (Nmales  = 11,942; Nfemales  = 15,641). Leveraging harmonized memory composite scores from four cohorts of cognitive aging and AD, we performed sex-stratified and sex-interaction genome-wide association studies in 24,216 non-Hispanic White and 3367 non-Hispanic Black participants.

RESULTS

We identified three sex-specific loci (rs67099044-CBLN2, rs719070-SCHIP1/IQCJ-SCHIP), including an X-chromosome locus (rs5935633-EGL6/TCEANC/OFD1), that associated with memory. Additionally, we identified heparan sulfate signaling as a sex-specific pathway and found sex-specific genetic correlations between memory and cardiovascular, immune, and education traits.

DISCUSSION

This study showed memory is highly and comparably heritable across sexes, as well as highlighted novel sex-specific genes, pathways, and genetic correlations that related to late-life memory.

HIGHLIGHTS

Demonstrated the heritable component of late-life memory is similar across sexes. Identified two genetic loci with a sex-interaction with baseline memory. Identified an X-chromosome locus associated with memory decline in females. Highlighted sex-specific candidate genes and pathways associated with memory. Revealed sex-specific shared genetic architecture between memory and complex traits.

Heart rate variability: Evaluating a potential biomarker of anxiety disorders

Establishing quantifiable biological markers associated with anxiety will increase the objectivity of phenotyping and enhance genetic research of anxiety disorders. Heart rate variability (HRV) is a physiological measure reflecting the dynamic relationship between the sympathetic and parasympathetic nervous systems, and is a promising target for further investigation. This review summarizes evidence evaluating HRV as a potential physiological biomarker of anxiety disorders by highlighting literature related to anxiety and HRV combined with investigations of endophenotypes, neuroimaging, treatment response, and genetics. Deficient HRV shows promise as an endophenotype of pathological anxiety and may serve as a noninvasive index of prefrontal cortical control over the amygdala, and potentially aid with treatment outcome prediction. We propose that the genetics of HRV can be used to enhance the understanding of the genetics of pathological anxiety for etiological investigations and treatment prediction. Given the anxiety-HRV link, strategies are offered to advance genetic analytical approaches, including the use of polygenic methods, wearable devices, and pharmacogenetic study designs. Overall, HRV shows promising support as a physiological biomarker of pathological anxiety, potentially in a transdiagnostic manner, with the heart-brain entwinement providing a novel approach to advance anxiety treatment development.

Obstructive sleep apnea heterogeneity and autonomic function: a role for heart rate variability in therapy selection and efficacy monitoring

Obstructive sleep apnea is a highly prevalent sleep-related breathing disorder, resulting in a disturbed breathing pattern, changes in blood gases, abnormal autonomic regulation, metabolic fluctuation, poor neurocognitive performance, and increased cardiovascular risk. With broad inter-individual differences recognised in risk factors, clinical symptoms, gene expression, physiological characteristics, and health outcomes, various obstructive sleep apnea subtypes have been identified. Therapeutic efficacy and its impact on outcomes, particularly for cardiovascular consequences, may also vary depending on these features in obstructive sleep apnea. A number of interventions such as positive airway pressure therapies, oral appliance, surgical treatment, and pharmaceutical options are available in clinical practice. Selecting an effective obstructive sleep apnea treatment and therapy is a challenging medical decision due to obstructive sleep apnea heterogeneity and numerous treatment modalities. Thus, an objective marker for clinical evaluation is warranted to estimate the treatment response in patients with obstructive sleep apnea. Currently, while the Apnea-Hypopnea Index is used for severity assessment of obstructive sleep apnea and still considered a major guide to diagnosis and managements of obstructive sleep apnea, the Apnea-Hypopnea Index is not a robust marker of symptoms, function, or outcome improvement. Abnormal cardiac autonomic modulation can provide additional insight to better understand obstructive sleep apnea phenotyping. Heart rate variability is a reliable neurocardiac tool to assess altered autonomic function and can also provide cardiovascular information in obstructive sleep apnea. Beyond the Apnea-Hypopnea Index, this review aims to discuss the role of heart rate variability as an indicator and predictor of therapeutic efficacy to different modalities in order to optimise tailored treatment for obstructive sleep apnea.

Resting heart rate causally affects the brain cortical structure: Mendelian randomization study

Resting heart rate (RHR) has been linked to impaired cortical structure in observational studies. However, the extent to which this association is potentially causal has not been determined. Using genetic data, this study aimed to reveal the causal effect of RHR on brain cortical structure. A Two-Sample Mendelian randomization (MR) analysis was conducted. Sensitivity analyses, weighted median, MR Pleiotropy residual sum and outlier, and MR-Egger regression were conducted to evaluate heterogeneity and pleiotropy. A causal relationship between RHR and cortical structures was identified by MR analysis. On the global scale, elevated RHR was found to decrease global surface area (SA; P < 0.0125). On a regional scale, the elevated RHR significantly decreased the SA of pars triangularis without global weighted (P = 1.58 × 10-4) and the thickness (TH) of the paracentral with global weighted (P = 3.56 × 10-5), whereas it increased the TH of banks of the superior temporal sulcus in the presence of global weighted (P = 1.04 × 10-4). MR study provided evidence that RHR might be causally linked to brain cortical structure, which offers a different way to understand the heart-brain axis theory.

A genome-wide association study of 24-hour urinary excretion of endocrine disrupting chemicals

Ubiquitous exposure to environmental endocrine disrupting chemicals (EDCs) instigates a major public health problem, but much remains unknown on the inter-individual differences in metabolism and excretion of EDCs. To examine this we performed a two-stage genome-wide association study (GWAS) for 24-hour urinary excretions of four parabens, two bisphenols, and nine phthalate metabolites. Results showed five genome-wide significant (p-value < 5x10-8) and replicated single nucleotide polymorphisms (SNPs) representing four independent signals that associated with mono-(2-ethyl-5-carboxypentyl) phthalate (MECPP) and mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP). Three of the four signals were located on chromosome 10 in a locus harboring the cytochrome P450 (CYP) genes CYP2C9, CYP2C58P, and CYP2C19 (rs117529685, pMECPP = 5.38x10-25; rs117033379, pMECPP = 1.96x10-19; rs4918798, pMECPP = 4.01x10-71; rs7895726, pMEHHP = 1.37x10-15, r2 with rs4918798 = 0.93). The other signal was on chromosome 6 close to the solute carrier (SLC) genes SLC17A1, SLC17A3, SLC17A4, and SCGN (rs1359232, pMECPP = 7.6x10-16). These four SNPs explained a substantial part (8.3 % - 9.2 %) of the variance in MECPP in the replication cohort. Bioinformatics analyses supported a likely causal role of CYP2C9 and SLC17A1 in metabolism and excretion of MECPP and MEHHP. Our results provide biological insights into mechanisms of phthalate metabolism and excretion with a likely causal role for CYP2C9 and SLC17A1.

Heart rate variability and cardiovascular diseases: A Mendelian randomization study

BACKGROUND

The causal relationship between heart rate variability and cardiovascular diseases and the associated events is still unclear, and the conclusions of current studies are inconsistent. We aimed to explore the relationship between heart rate variability and cardiovascular diseases and the associated events with the Mendelian randomization study.

METHODS

We selected normal-to-normal inter-beat intervals (SDNN), root mean square of the successive differences of inter-beat intervals (RMSSD) and peak-valley respiratory sinus arrhythmia or high-frequency power (pvRSA/HF) as the three sets of instrumental variables for heart rate variability. The outcome for cardiovascular diseases included essential hypertension, heart failure, angina pectoris, myocardial infarction, nonischemic cardiomyopathy and arrhythmia. Cardiac arrest, cardiac death and major coronary heart disease event were defined as the related events of cardiovascular diseases. The data for exposures and outcomes were derived from publicly available genome-wide association studies. Inverse variance weighted was used for the main causal estimation. Analyses of heterogeneity and pleiotropy were conducted using the Cochran Q test of Inverse variance weighted and MR-Egger, leave-one-out analysis, and MR-Pleiotropy Residual Sum and Outlier methods.

RESULTS

The Inverse variance weighted method indicated that genetically predicted pvRSA/HF was associated with the increased risk of cardiac arrest (odds ratio 2.02, 95% confidence interval 1.25-3.28, p = .004). The results were free of heterogeneity and pleiotropy. There were no outliers and the leave-one-out analysis proved that the results were reliable.

CONCLUSIONS

This study provides genetic evidence that pvRSA/HF is causally related to cardiac arrest.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Integrated omics analysis of coronary artery calcifications and myocardial infarction: the Framingham Heart Study

Gene function can be described using various measures. We integrated association studies of three types of omics data to provide insights into the pathophysiology of subclinical coronary disease and myocardial infarction (MI). Using multivariable regression models, we associated: (1) single nucleotide polymorphism, (2) DNA methylation, and (3) gene expression with coronary artery calcification (CAC) scores and MI. Among 3106 participants of the Framingham Heart Study, 65 (2.1%) had prevalent MI and 60 (1.9%) had incident MI, median CAC value was 67.8 [IQR 10.8, 274.9], and 1403 (45.2%) had CAC scores > 0 (prevalent CAC). Prevalent CAC was associated with AHRR (linked to smoking) and EXOC3 (affecting platelet function and promoting hemostasis). CAC score was associated with VWA1 (extracellular matrix protein associated with cartilage structure in endomysium). For prevalent MI we identified FYTTD1 (down-regulated in familial hypercholesterolemia) and PINK1 (linked to cardiac tissue homeostasis and ischemia-reperfusion injury). Incident MI was associated with IRX3 (enhancing browning of white adipose tissue) and STXBP3 (controlling trafficking of glucose transporter type 4 to plasma). Using an integrative trans-omics approach, we identified both putatively novel and known candidate genes associated with CAC and MI. Replication of findings is warranted.

Polygenic Risk Score for Cardiovascular Diseases in Artificial Intelligence Paradigm: A Review

Cardiovascular disease (CVD) related mortality and morbidity heavily strain society. The relationship between external risk factors and our genetics have not been well established. It is widely acknowledged that environmental influence and individual behaviours play a significant role in CVD vulnerability, leading to the development of polygenic risk scores (PRS). We employed the PRISMA search method to locate pertinent research and literature to extensively review artificial intelligence (AI)-based PRS models for CVD risk prediction. Furthermore, we analyzed and compared conventional vs. AI-based solutions for PRS. We summarized the recent advances in our understanding of the use of AI-based PRS for risk prediction of CVD. Our study proposes three hypotheses: i) Multiple genetic variations and risk factors can be incorporated into AI-based PRS to improve the accuracy of CVD risk predicting. ii) AI-based PRS for CVD circumvents the drawbacks of conventional PRS calculators by incorporating a larger variety of genetic and non-genetic components, allowing for more precise and individualised risk estimations. iii) Using AI approaches, it is possible to significantly reduce the dimensionality of huge genomic datasets, resulting in more accurate and effective disease risk prediction models. Our study highlighted that the AI-PRS model outperformed traditional PRS calculators in predicting CVD risk. Furthermore, using AI-based methods to calculate PRS may increase the precision of risk predictions for CVD and have significant ramifications for individualized prevention and treatment plans.

Long-term association of ultra-short heart rate variability with cardiovascular events

Heart rate variability (HRV) is a cardiac autonomic marker with predictive value in cardiac patients. Ultra-short HRV (usHRV) can be measured at scale using standard and wearable ECGs, but its association with cardiovascular events in the general population is undetermined. We aimed to validate usHRV measured using ≤ 15-s ECGs (using RMSSD, SDSD and PHF indices) and investigate its association with atrial fibrillation, major adverse cardiac events, stroke and mortality in individuals without cardiovascular disease. In the National Survey for Health and Development (n = 1337 participants), agreement between 15-s and 6-min HRV, assessed with correlation analysis and Bland-Altman plots, was very good for RMSSD and SDSD and good for PHF. In the UK Biobank (n = 51,628 participants, 64% male, median age 58), after a median follow-up of 11.5 (11.4-11.7) years, incidence of outcomes ranged between 1.7% and 4.3%. Non-linear Cox regression analysis showed that reduced usHRV from 15-, 10- and 5-s ECGs was associated with all outcomes. Individuals with low usHRV (< 20th percentile) had hazard ratios for outcomes between 1.16 and 1.29, p < 0.05, with respect to the reference group. In conclusion, usHRV from ≤ 15-s ECGs correlates with standard short-term HRV and predicts increased risk of cardiovascular events in a large population-representative cohort.

Phenotypic but not genetically predicted heart rate variability associated with all-cause mortality

Low heart rate variability (HRV) has been widely reported as a predictor for increased mortality. However, the molecular mechanisms are poorly understood. Therefore, this study aimed to identify novel genetic loci associated with HRV and assess the association of phenotypic HRV and genetically predicted HRV with mortality. In a GWAS of 46,075 European ancestry individuals from UK biobank, we identified 17 independent genome-wide significant genetic variants in 16 loci associated with HRV traits. Notably, eight of these loci (RNF220, GNB4, LINCR-002, KLHL3/HNRNPA0, CHRM2, KCNJ5, MED13L, and C160rf72) have not been reported previously. In a prospective phenotypic relationship between HRV and mortality during a median follow-up of seven years, individuals with lower HRV had higher risk of dying from any cause. Genetically predicted HRV, as determined by the genetic risk scores, was not associated with mortality. To the best of our knowledge, the findings provide novel biological insights into the mechanisms underlying HRV. These results also underline the role of the cardiac autonomic nervous system, as indexed by HRV, in predicting mortality.

Heart Rate Variability and Incident Type 2 Diabetes in General Population

CONTEXT

Hyperglycemia and autonomic dysfunction are bidirectionally related.

OBJECTIVE

We investigated the association of longitudinal evolution of heart rate variability (HRV) with incident type 2 diabetes (T2D) among the general population.

METHODS

We included 7630 participants (mean age 63.7 years, 58% women) from the population-based Rotterdam Study who had no history of T2D and atrial fibrillation at baseline and had repeated HRV assessments at baseline and during follow-up. We used joint models to assess the association between longitudinal evolution of heart rate and different HRV metrics (including the heart rate-corrected SD of the normal-to-normal RR intervals [SDNNc], and root mean square of successive RR-interval differences [RMSSDc]) with incident T2D. Models were adjusted for cardiovascular risk factors. Bidirectional Mendelian randomization (MR) using summary-level data was also performed.

RESULTS

During a median follow-up of 8.6 years, 871 individuals developed incident T2D. One SD increase in heart rate (hazard ratio [HR] 1.20; 95% CI, 1.09-1.33), and log(RMSSDc) (HR 1.16; 95% CI, 1.01-1.33) were independently associated with incident T2D. The HRs were 1.54 (95% CI, 1.08-2.06) for participants younger than 62 years and 1.15 (95% CI, 1.01-1.31) for those older than 62 years for heart rate (P for interaction <.001). Results from bidirectional MR analyses suggested that HRV and T2D were not significantly related to each other.

CONCLUSION

Autonomic dysfunction precedes development of T2D, especially among younger individuals, while MR analysis suggests no causal relationship. More studies are needed to further validate our findings.

RGS7 balances acetylation/de-acetylation of p65 to control chemotherapy-dependent cardiac inflammation

Cardiotoxicity remains a major limitation in the clinical utility of anthracycline chemotherapeutics. Regulator of G-protein Signaling 7 (RGS7) and inflammatory markers are up-regulated in the hearts of patients with a history of chemotherapy particularly those with reduced left-ventricular function. RGS7 knockdown in either the murine myocardium or isolated murine ventricular cardiac myocytes (VCM) or cultured human VCM provided marked protection against doxorubicin-dependent oxidative stress, NF-κB activation, inflammatory cytokine production, and cell death. In exploring possible mechanisms causally linking RGS7 to pro-inflammatory signaling cascades, we found that RGS7 forms a complex with acetylase Tip60 and deacetylase sirtuin 1 (SIRT1) and controls the acetylation status of the p65 subunit of NF-κB. In VCM, the detrimental impact of RGS7 could be mitigated by inhibiting Tip60 or activating SIRT1, indicating that the ability of RGS7 to modulate cellular acetylation capacity is critical for its pro-inflammatory actions. Further, RGS7-driven, Tip60/SIRT1-dependent cytokines released from ventricular cardiac myocytes and transplanted onto cardiac fibroblasts increased oxidative stress, markers of transdifferentiation, and activity of extracellular matrix remodelers emphasizing the importance of the RGS7-Tip60-SIRT1 complex in paracrine signaling in the myocardium. Importantly, while RGS7 overexpression in heart resulted in sterile inflammation, fibrotic remodeling, and compromised left-ventricular function, activation of SIRT1 counteracted the detrimental impact of RGS7 in heart confirming that RGS7 increases acetylation of SIRT1 substrates and thereby drives cardiac dysfunction. Together, our data identify RGS7 as an amplifier of inflammatory signaling in heart and possible therapeutic target in chemotherapeutic drug-induced cardiotoxicity.

Integrating Genome-wide information and Wearable Device Data to Explore the Link of Anxiety and Antidepressants with Heart Rate Variability

Background

Anxiety disorders are associated with decreased heart rate variability (HRV), but the underlying mechanisms remain elusive.

Methods

We selected individuals with whole-genome sequencing, Fitbit, and electronic health record data (N=920; 61,333 data points) from the All of Us Research Program. Anxiety PRS were derived with PRS-CS after meta-analyzing anxiety genome-wide association studies from three major cohorts-UK Biobank, FinnGen, and the Million Veterans Program (N Total =364,550). The standard deviation of average RR intervals (SDANN) was calculated using five-minute average RR intervals over full 24-hour heart rate measurements. Antidepressant exposure was defined as an active antidepressant prescription at the time of the HRV measurement in the EHR. The associations of daily SDANN measurements with the anxiety PRS, antidepressant classes, and antidepressant substances were tested. Participants with lifetime diagnoses of cardiovascular disorders, diabetes mellitus, and major depression were excluded in sensitivity analyses. One-sample Mendelian randomization (MR) was employed to assess potential causal effect of anxiety on SDANN.

Results

Anxiety PRS was independently associated with reduced SDANN (beta=-0.08; p=0.003). Of the eight antidepressant medications and four classes tested, venlafaxine (beta=-0.12, p=0.002) and bupropion (beta=-0.071, p=0.01), tricyclic antidepressants (beta=-0.177, p=0.0008), selective serotonin reuptake inhibitors (beta=-0.069; p=0.0008) and serotonin and norepinephrine reuptake inhibitors (beta=-0.16; p=2×10 -6 ) were associated with decreased SDANN. One-sample MR indicated an inverse effect of anxiety on SDANN (beta=-2.22, p=0.03).

Conclusions

Anxiety and antidepressants are independently associated with decreased HRV, and anxiety appears to exert a causal effect on HRV. Our observational findings provide novel insights into the impact of anxiety on HRV.

Genome-wide association studies: utility and limitations for research in physiology

Physiological systems are subject to interindividual variation encoded by genetics. Genome-wide association studies (GWAS) operate by surveying thousands of genetic variants from a substantial number of individuals and assessing their association to a trait of interest, be it a physiological variable, a molecular phenotype (e.g. gene expression), or even a disease or condition. Through a myriad of methods, GWAS downstream analyses then explore the functional consequences of each variant and attempt to ascertain a causal relationship to the phenotype of interest, as well as to delve into its links to other traits. This type of investigation allows mechanistic insights into physiological functions, pathological disturbances and shared biological processes between traits (i.e. pleiotropy). An exciting example is the discovery of a new thyroid hormone transporter (SLC17A4) and hormone metabolising enzyme (AADAT) from a GWAS on free thyroxine levels. Therefore, GWAS have substantially contributed with insights into physiology and have been shown to be useful in unveiling the genetic control underlying complex traits and pathological conditions; they will continue to do so with global collaborations and advances in genotyping technology. Finally, the increasing number of trans-ancestry GWAS and initiatives to include ancestry diversity in genomics will boost the power for discoveries, making them also applicable to non-European populations.

Multi-Omics Data Analysis Identifies Prognostic Biomarkers across Cancers

Combining omics data from different layers using integrative methods provides a better understanding of the biology of a complex disease such as cancer. The discovery of biomarkers related to cancer development or prognosis helps to find more effective treatment options. This study integrates multi-omics data of different cancer types with a network-based approach to explore common gene modules among different tumors by running community detection methods on the integrated network. The common modules were evaluated by several biological metrics adapted to cancer. Then, a new prognostic scoring method was developed by weighting mRNA expression, methylation, and mutation status of genes. The survival analysis pointed out statistically significant results for GNG11, CBX2, CDKN3, ARHGEF10, CLN8, SEC61G and PTDSS1 genes. The literature search reveals that the identified biomarkers are associated with the same or different types of cancers. Our method does not only identify known cancer-specific biomarker genes, but also proposes new potential biomarkers. Thus, this study provides a rationale for identifying new gene targets and expanding treatment options across cancer types.

Resting heart rate and antisocial behaviour: a Mendelian randomisation study

Observational studies frequently report phenotypic associations between low resting heart rate (RHR) and higher levels of antisocial behaviour (ASB), although it remains unclear whether this relationship reflects causality. To triangulate evidence, we conducted two-sample univariable Mendelian randomisation (MR), multivariable MR and linkage disequilibrium score regression (LDSC) analyses. Genetic data were accessed from published genome-wide association studies (GWAS) for RHR (n = 458,835) and ASB (n = 85,359) for the univariable analyses, along with a third GWAS for heart rate variability (HRV; n = 53,174) for all other analyses. Genome-wide significant (p < 5 × 10-8) single-nucleotide polymorphisms associated with RHR (n = 278) were selected as instrumental variables and the outcome was a composite measure of ASB. No causal association was observed between RHR and ASB (BIVW =  - 0.0004, p = 0.841). The multivariable MR analyses including RHR and HRV also suggested no causal associations (BIVW = 0.016, p = 0.914) and no genetic correlations between the heart rate measures and ASB were observed using LDSC (rg = 0.057, p = 0.169). Sensitivity analyses suggested that our results are not likely to be affected by heterogeneity, pleiotropic effects, or reverse causation. These findings suggest that individual differences in autonomic nervous system functioning indexed by RHR are not likely to directly contribute to the development of ASB. Therefore, previously observed associations between RHR and ASB may arise from confounding, reverse causation, and/or additional study characteristics. Further causally informative longitudinal research is required to confirm our findings, and caution should be applied when using measures of RHR in interventions targeting ASB.

Associations of schizophrenia with arrhythmic disorders and electrocardiogram traits: an in-depth genetic exploration of population samples

Background

An important contributor to the decreased life expectancy of individuals with schizophrenia is sudden cardiac death. While arrhythmic disorders play an important role in this, the nature of the relation between schizophrenia and arrhythmia is not fully understood.

Methods

We leveraged summary-level data of large-scale genome-wide association studies of schizophrenia (53,386 cases 77,258 controls), arrhythmic disorders (atrial fibrillation, 55,114 cases 482,295 controls; Brugada syndrome, 2,820 cases 10,001 controls) and electrocardiogram traits (heart rate (variability), PR interval, QT interval, JT interval, and QRS duration, n=46,952-293,051). First, we examined shared genetic liability by assessing global and local genetic correlations and conducting functional annotation. Next, we explored bidirectional causal relations between schizophrenia and arrhythmic disorders and electrocardiogram traits using Mendelian randomization.

Outcomes

There was no evidence for global genetic correlations, except between schizophrenia and Brugada (rg=0·14, p=4·0E-04). In contrast, strong positive and negative local genetic correlations between schizophrenia and all cardiac traits were found across the genome. In the strongest associated regions, genes related to immune system and viral response mechanisms were overrepresented. Mendelian randomization indicated a causal, increasing effect of liability to schizophrenia on Brugada syndrome (OR=1·15, p=0·009) and heart rate during activity (beta=0·25, p=0·015).

Interpretation

While there was little evidence for global genetic correlations, specific genomic regions and biological pathways important for both schizophrenia and arrhythmic disorders and electrocardiogram traits emerged. The putative causal effect of liability to schizophrenia on Brugada warrants increased cardiac monitoring and potentially early medical intervention in patients with schizophrenia.

Funding

European Research Council Starting Grant.

Phenome-wide Mendelian randomization study evaluating the association of circulating vitamin D with complex diseases

Background

Circulating vitamin D has been associated with multiple clinical diseases in observational studies, but the association was inconsistent due to the presence of confounders. We conducted a bidirectional Mendelian randomization (MR) study to explore the healthy atlas of vitamin D in many clinical traits and evaluate their causal association.

Methods

Based on a large-scale genome-wide association study (GWAS), the single nucleotide polymorphism (SNPs) instruments of circulating 25-hydroxyvitamin D (25OHD) from 443,734 Europeans and the corresponding effects of 10 clinical diseases and 42 clinical traits in the European population were recruited to conduct a bidirectional two-sample Mendelian randomization study. Under the network of Mendelian randomization analysis, inverse-variance weighting (IVW), weighted median, weighted mode, and Mendelian randomization (MR)–Egger regression were performed to explore the causal effects and pleiotropy. Mendelian randomization pleiotropy RESidual Sum and Outlier (MR-PRESSO) was conducted to uncover and exclude pleiotropic SNPs.

Results

The results revealed that genetically decreased vitamin D was inversely related to the estimated BMD (β = −0.029 g/cm2, p = 0.027), TC (β = −0.269 mmol/L, p = 0.006), TG (β = −0.208 mmol/L, p = 0.002), and pulse pressure (β = −0.241 mmHg, p = 0.043), while positively associated with lymphocyte count (β = 0.037%, p = 0.015). The results did not reveal any causal association of vitamin D with clinical diseases. On the contrary, genetically protected CKD was significantly associated with increased vitamin D (β = 0.056, p = 2.361 × 10−26).

Conclusion

The putative causal effects of circulating vitamin D on estimated bone mass, plasma triglyceride, and total cholesterol were uncovered, but not on clinical diseases. Vitamin D may be linked to clinical disease by affecting health-related metabolic markers.

Implication of heart rate variability on cerebral small vessel disease: A potential therapeutic target

OBJECTIVE

This study aimed to investigate the relationships of heart rate variability (HRV) with the presence, severity, and individual neuroimaging markers of cerebral small vessel disease (CSVD).

METHOD

A total of 4676 participants from the Third China National Stroke Registry (CNSR-III) study were included in this cross-sectional analysis. CSVD and its markers, including white matter hyperintensity (WMH), lacunes, enlarged perivascular spaces (EPVS), cerebral microbleeds (CMBs), and brain atrophy (BA), were evaluated. Two common HRV parameters, including the square root of the mean of the sum of the squares of differences between adjacent N-N intervals (RMSSD) and the standard deviation of all N-N intervals (SDNN), were used to evaluate the function of the autonomic nervous system (ANS). Binary or ordinal logistic regression analyses were performed to investigate the association between HRV and CSVD. In addition, two-sample mendelian randomization (MR) analyses were performed to investigate the causality of HRV with CSVD.

RESULTS

RMSSD was significantly associated with total burden of CSVD (Wardlaw's scale, common odds ratio [cOR] 0.80, 95% confidence interval [CI] 0.67-0.96, p = 0.02; Rothwell's scale, cOR 0.75, 95% CI 0.60-0.93, p = 0.008) and the presence of CSVD (Rothwell, OR 0.75, 95% CI 0.60-0.93, p = 0.008). However, no significant associations between SDNN and the presence or total burden of CSVD were observed. Moreover, RMSSD was related to WMH burden (OR 0.80, 95% CI 0.66-0.96, p = 0.02), modified WMH burden (cOR 0.82, 95% CI 0.69-0.97, p = 0.02), and Deep-WMH (OR 0.75, 95% CI 0.62-0.91, p = 0.003), while SDNN was related to Deep-WMH (OR 0.80, 95% CI 0.66-0.96, p = 0.02) and BA (cOR 0.80, 95% CI 0.68-0.95, p = 0.009). Furthermore, adding HRV to the conventional model based on vascualr risk factors enhanced the predictive performance for CSVD, as validated by the integrated discrimination index (p < 0.05). In addition, no causality between HRV and CSVD was observed in two-sample MR analyses.

CONCLUSION

Decreased HRV may be a potential risk factor of CSVD, implying the possible role of the ANS in the pathogenesis of CSVD.

A comparative evaluation of resting state proxies of sympathetic and parasympathetic nervous system activity in adolescent major depression

Major depressive disorder (MDD) is associated with autonomic nervous system (ANS) dysfunction, characterized by decreased parasympathetic (PNS) and increased sympathetic (SNS) activity. Although findings on reduced PNS activity in adult MDD have been replicated in adolescents, comprehensive studies assessing PNS and SNS proxies in underage patients with MDD are scarce. Proxies of resting PNS (heart rate variability (HRV) and SNS activity (skin conductance response [SCR] and salivary alpha amylase [sAA], as well as mixed activity (heart rate [HR]) were collected in adolescents with MDD (n = 29) and non-depressed controls (n = 29). Primary analyses addressed differences between groups and correlations with depression severity. Patients with MDD showed significantly decreased HRV (g = - 0.87; 95% CI [- 1.39; - 0.35]) and increased HR (g = 0.66; 95% CI [0.14; 1.18]). Proxies of pure SNS activity showed no significant differences between groups. HR (positive) and HRV (negative) were significantly correlated with self- and clinician-rated depression severity. Alterations of ANS activity are evident in adolescent MDD, but characterized by decreased PNS activity only. We found no evidence for altered SNS activity. Findings suggest that ANS dysfunction early in the course of MDD might be predominantly driven by decreased PNS activity.

Cross-ancestry genome-wide analysis of atrial fibrillation unveils disease biology and enables cardioembolic risk prediction

Atrial fibrillation (AF) is a common cardiac arrhythmia resulting in increased risk of stroke. Despite highly heritable etiology, our understanding of the genetic architecture of AF remains incomplete. Here we performed a genome-wide association study in the Japanese population comprising 9,826 cases among 150,272 individuals and identified East Asian-specific rare variants associated with AF. A cross-ancestry meta-analysis of >1 million individuals, including 77,690 cases, identified 35 new susceptibility loci. Transcriptome-wide association analysis identified IL6R as a putative causal gene, suggesting the involvement of immune responses. Integrative analysis with ChIP-seq data and functional assessment using human induced pluripotent stem cell-derived cardiomyocytes demonstrated ERRg as having a key role in the transcriptional regulation of AF-associated genes. A polygenic risk score derived from the cross-ancestry meta-analysis predicted increased risks of cardiovascular and stroke mortalities and segregated individuals with cardioembolic stroke in undiagnosed AF patients. Our results provide new biological and clinical insights into AF genetics and suggest their potential for clinical applications.

A RGS7-CaMKII complex drives myocyte-intrinsic and myocyte-extrinsic mechanisms of chemotherapy-induced cardiotoxicity

Dose-limiting cardiotoxicity remains a major limitation in the clinical use of cancer chemotherapeutics. Here, we describe a role for Regulator of G protein Signaling 7 (RGS7) in chemotherapy-dependent heart damage, the demonstration for a functional role of RGS7 outside of the nervous system and retina. Though expressed at low levels basally, we observed robust up-regulation of RGS7 in the human and murine myocardium following chemotherapy exposure. In ventricular cardiomyocytes (VCM), RGS7 forms a complex with Ca²⁺/calmodulin-dependent protein kinase (CaMKII) supported by key residues (K412 and P391) in the RGS domain of RGS7. In VCM treated with chemotherapeutic drugs, RGS7 facilitates CaMKII oxidation and phosphorylation and CaMKII-dependent oxidative stress, mitochondrial dysfunction, and apoptosis. Cardiac-specific RGS7 knockdown protected the heart against chemotherapy-dependent oxidative stress, fibrosis, and myocyte loss and improved left ventricular function in mice treated with doxorubicin. Conversely, RGS7 overexpression induced fibrosis, reactive oxygen species generation, and cell death in the murine myocardium that were mitigated following CaMKII inhibition. RGS7 also drives production and release of the cardiokine neuregulin-1, which facilitates paracrine communication between VCM and neighboring vascular endothelial cells (EC), a maladaptive mechanism contributing to VCM dysfunction in the failing heart. Importantly, while RGS7 was both necessary and sufficient to facilitate chemotherapy-dependent cytotoxicity in VCM, RGS7 is dispensable for the cancer-killing actions of these same drugs. These selective myocyte-intrinsic and myocyte-extrinsic actions of RGS7 in heart identify RGS7 as an attractive therapeutic target in the mitigation of chemotherapy-driven cardiotoxicity.

The Relevance of GIRK Channels in Heart Function

Among the large number of potassium-channel families implicated in the control of neuronal excitability, G-protein-gated inwardly rectifying potassium channels (GIRK/Kir3) have been found to be a main factor in heart control. These channels are activated following the modulation of G-protein-coupled receptors and, although they have been implicated in different neurological diseases in both human and animal studies of the central nervous system, the therapeutic potential of different subtypes of these channel families in cardiac conditions has remained untapped. As they have emerged as a promising potential tool to treat a variety of conditions that disrupt neuronal homeostasis, many studies have started to focus on these channels as mediators of cardiac dynamics, thus leading to research into their implication in cardiovascular conditions. Our aim is to review the latest advances in GIRK modulation in the heart and their role in the cardiovascular system.

Consequences of somatic mutations of GIRK1 detected in primary malign tumors on expression and function of G-protein activated, inwardly rectifying, K+ channels

A search in the GDC Data Portal revealed 304 documented somatic mutations of the KCNJ3 gene in primary tumors (out of 10.202 cases). Most affected tumor types were carcinomas from uterus, skin and lung, while breast cancer exerted the lowest number of somatic mutations. We focused our research on 15 missense mutations within the region between TM1 and TM2, comprising the pore helix and ion selectivity signature. Expression was measured by confocal laser scan microscopy of eGFP tagged GIRK1 subunits, expressed with and without GIRK4 in oocytes of Xenopus laevis. GIRK ion currents were activated via coexpressed m2Rs and measured by the Two Electrode Voltage Clamp technique. Magnitude of the total GIRK current, as well as the fraction of current inducible by the agonist, were measured. Ion selectivity was gauged by assessment of the PNa+/PK+ ratio, calculated by the GIRK current reversal potential in extracellular media at different Na+ and K+ concentrations. None of the tested mutations was able to form functional GIRK1 homooligomeric ion channels. One of the mutations, G145A, which locates directly to the ion selectivity signature, exerted an increased PNa+/PK+ ratio. Generally, the missense mutations studied can be categorized into three groups: (i) normal/reduced expression accompanied by reduced/absent function (S132Y, F136L, E139K, G145A, R149Q, R149P, G178D, S185Y, Q186R), (ii) normal/increased expression as well as increased function (E140M, A142T, M184I) and (iii) miniscule expression but increased function relative to expression levels (I151N, G158S). We conclude, that gain of function mutations, identical or similar to categories (ii) and (iii), may potentially be involved in genesis and progression of malignancies in tissues that exert a high rate of occurrence of somatic mutations of KCNJ3.

Heart rate variability and atrial fibrillation in the general population: a longitudinal and Mendelian randomization study

BACKGROUND

Sex differences and causality of the association between heart rate variability (HRV) and atrial fibrillation (AF) in the general population remain unclear.

METHODS

12,334 participants free of AF from the population-based Rotterdam Study were included. Measures of HRV including the standard deviation of normal RR intervals (SDNN), SDNN corrected for heart rate (SDNNc), RR interval differences (RMSSD), RMSSD corrected for heart rate (RMSSDc), and heart rate were assessed at baseline and follow-up examinations. Joint models, adjusted for cardiovascular risk factors, were used to determine the association between longitudinal measures of HRV with new-onset AF. Genetic variants for HRV were used as instrumental variables in a Mendelian randomization (MR) analysis using genome-wide association studies (GWAS) summary-level data.

RESULTS

During a median follow-up of 9.4 years, 1302 incident AF cases occurred among 12,334 participants (mean age 64.8 years, 58.3% women). In joint models, higher SDNN (fully-adjusted hazard ratio (HR), 95% confidence interval (CI) 1.24, 1.04-1.47, p = 0.0213), and higher RMSSD (fully-adjusted HR, 95% CI 1.33, 1.13-1.54, p = 0.0010) were significantly associated with new-onset AF. Sex-stratified analyses showed that the associations were mostly prominent among women. In MR analyses, a genetically determined increase in SDNN (odds ratio (OR), 95% CI 1.60, 1.27-2.02, p = 8.36 × 10^-05), and RMSSD (OR, 95% CI 1.56, 1.31-1.86, p = 6.32 × 10^-07) were significantly associated with an increased odds of AF.

CONCLUSION

Longitudinal measures of uncorrected HRV were significantly associated with new-onset AF, especially among women. MR analyses supported the causal relationship between uncorrected measures of HRV with AF. Our findings indicate that measures to modulate HRV might prevent AF in the general population, in particular in women. AF; atrial fibrillation, GWAS; genome-wide association study, IVW; inverse variance weighted, MR; Mendelian randomization, MR-PRESSO; MR-egger and mendelian randomization pleiotropy residual sum and outlier, RMSSD; root mean square of successive RR interval differences, RMSSDc; root mean square of successive RR interval differences corrected for heart rate, SDNN; standard deviation of normal to normal RR intervals, SDNNc; standard deviation of normal to normal RR intervals corrected for heart rate, WME; weighted median estimator. ^aRotterdam Study n=12,334 ^bHRV GWAS n=53,174 ^cAF GWAS n=1,030,836.

A Polygenic Risk Score Based on a Cardioembolic Stroke Multitrait Analysis Improves a Clinical Prediction Model for This Stroke Subtype

Background

Occult atrial fibrillation (AF) is one of the major causes of embolic stroke of undetermined source (ESUS). Knowing the underlying etiology of an ESUS will reduce stroke recurrence and/or unnecessary use of anticoagulants. Understanding cardioembolic strokes (CES), whose main cause is AF, will provide tools to select patients who would benefit from anticoagulants among those with ESUS or AF. We aimed to discover novel loci associated with CES and create a polygenetic risk score (PRS) for a more efficient CES risk stratification.

Methods

Multitrait analysis of GWAS (MTAG) was performed with MEGASTROKE-CES cohort (n = 362,661) and AF cohort (n = 1,030,836). We considered significant variants and replicated those variants with MTAG p-value < 5 × 10⁻⁸ influencing both traits (GWAS-pairwise) with a p-value < 0.05 in the original GWAS and in an independent cohort (n = 9,105). The PRS was created with PRSice-2 and evaluated in the independent cohort.

Results

We found and replicated eleven loci associated with CES. Eight were novel loci. Seven of them had been previously associated with AF, namely, CAV1, ESR2, GORAB, IGF1R, NEURL1, WIPF1, and ZEB2. KIAA1755 locus had never been associated with CES/AF, leading its index variant to a missense change (R1045W). The PRS generated has been significantly associated with CES improving discrimination and patient reclassification of a model with age, sex, and hypertension.

Conclusion

The loci found significantly associated with CES in the MTAG, together with the creation of a PRS that improves the predictive clinical models of CES, might help guide future clinical trials of anticoagulant therapy in patients with ESUS or AF.

Cardiovascular Health of Youth During Gender-Affirming Testosterone Treatment: A Review

PURPOSE

Many birth-assigned female/transgender male and nonbinary people (identified as masculine spectrum here) begin gender-affirming testosterone therapy by the age of 24 years. Few data inform assessment of cardiovascular health of masculine spectrum youth as a specific subgroup of the 1.5 million transgender people in the United States. The purpose of this review is to help youth-serving practitioners consider, understand, and evaluate cardiovascular health in adolescent and young adult masculine spectrum patients receiving gender-affirming testosterone treatment.

METHODS

This is a narrative review intended to synthesize a broad body of clinical and research literature.

RESULTS

Common cardiovascular health changes associated with testosterone include increased red blood cell mass and likely insignificant changes in high-density lipoprotein and low-density lipoprotein levels. Changes in heart mass, heart electrophysiology, and endothelial reactivity are likely, based on extrapolation of data from adults. Testosterone may have indirect effects on cardiovascular health through influences on depression, anxiety, stress, and anorexia nervosa as well as on behaviors such as tobacco use.

CONCLUSIONS

Testosterone contributes importantly to the cardiovascular health and well-being of masculine spectrum gender-diverse youth. We need to do a better job of supporting these young people with data on cardiovascular health over the life span.

Detecting Genetic Ancestry and Adaptation in the Taiwanese Han People

The Taiwanese people are composed of diverse indigenous populations and the Taiwanese Han. About 95% of the Taiwanese identify themselves as Taiwanese Han, but this may not be a homogeneous population because they migrated to the island from various regions of continental East Asia over a period of 400 years. Little is known about the underlying patterns of genetic ancestry, population admixture, and evolutionary adaptation in the Taiwanese Han people. Here, we analyzed the whole-genome single-nucleotide polymorphism genotyping data from 14,401 individuals of Taiwanese Han collected by the Taiwan Biobank and the whole-genome sequencing data for a subset of 772 people. We detected four major genetic ancestries with distinct geographic distributions (i.e., Northern, Southeastern, Japonic, and Island Southeast Asian ancestries) and signatures of population mixture contributing to the genomes of Taiwanese Han. We further scanned for signatures of positive natural selection that caused unusually long-range haplotypes and elevations of hitchhiked variants. As a result, we identified 16 candidate loci in which selection signals can be unambiguously localized at five single genes: CTNNA2, LRP1B, CSNK1G3, ASTN2, and NEO1. Statistical associations were examined in 16 metabolic-related traits to further elucidate the functional effects of each candidate gene. All five genes appear to have pleiotropic connections to various types of disease susceptibility and significant associations with at least one metabolic-related trait. Together, our results provide critical insights for understanding the evolutionary history and adaption of the Taiwanese Han population.

Key Characteristics of Cardiovascular Toxicants

BACKGROUND

The concept of chemical agents having properties that confer potential hazard called key characteristics (KCs) was first developed to identify carcinogenic hazards. Identification of KCs of cardiovascular (CV) toxicants could facilitate the systematic assessment of CV hazards and understanding of assay and data gaps associated with current approaches.

OBJECTIVES

We sought to develop a consensus-based synthesis of scientific evidence on the KCs of chemical and nonchemical agents known to cause CV toxicity along with methods to measure them.

METHODS

An expert working group was convened to discuss mechanisms associated with CV toxicity.

RESULTS

The group identified 12 KCs of CV toxicants, defined as exogenous agents that adversely interfere with function of the CV system. The KCs were organized into those primarily affecting cardiac tissue (numbers 1-4 below), the vascular system (5-7), or both (8-12), as follows: 1) impairs regulation of cardiac excitability, 2) impairs cardiac contractility and relaxation, 3) induces cardiomyocyte injury and death, 4) induces proliferation of valve stroma, 5) impacts endothelial and vascular function, 6) alters hemostasis, 7) causes dyslipidemia, 8) impairs mitochondrial function, 9) modifies autonomic nervous system activity, 10) induces oxidative stress, 11) causes inflammation, and 12) alters hormone signaling.

DISCUSSION

These 12 KCs can be used to help identify pharmaceuticals and environmental pollutants as CV toxicants, as well as to better understand the mechanistic underpinnings of their toxicity. For example, evidence exists that fine particulate matter [PM ≤2.5μm in aerodynamic diameter (PM2.5)] air pollution, arsenic, anthracycline drugs, and other exogenous chemicals possess one or more of the described KCs. In conclusion, the KCs could be used to identify potential CV toxicants and to define a set of test methods to evaluate CV toxicity in a more comprehensive and standardized manner than current approaches. https://doi.org/10.1289/EHP9321.

Rare Coding Variants Associated With Electrocardiographic Intervals Identify Monogenic Arrhythmia Susceptibility Genes: A Multi-Ancestry Analysis

BACKGROUND

Alterations in electrocardiographic (ECG) intervals are well-known markers for arrhythmia and sudden cardiac death (SCD) risk. While the genetics of arrhythmia syndromes have been studied, relations between electrocardiographic intervals and rare genetic variation at a population level are poorly understood.

METHODS

Using a discovery sample of 29 000 individuals with whole-genome sequencing from Trans-Omics in Precision Medicine and replication in nearly 100 000 with whole-exome sequencing from the UK Biobank and MyCode, we examined associations between low-frequency and rare coding variants with 5 routinely measured electrocardiographic traits (RR, P-wave, PR, and QRS intervals and corrected QT interval).

RESULTS

We found that rare variants associated with population-based electrocardiographic intervals identify established monogenic SCD genes (KCNQ1, KCNH2, and SCN5A), a controversial monogenic SCD gene (KCNE1), and novel genes (PAM and MFGE8) involved in cardiac conduction. Loss-of-function and pathogenic SCN5A variants, carried by 0.1% of individuals, were associated with a nearly 6-fold increased odds of the first-degree atrioventricular block (P=8.4×10^-5). Similar variants in KCNQ1 and KCNH2 (0.2% of individuals) were associated with a 23-fold increased odds of marked corrected QT interval prolongation (P=4×10^-25), a marker of SCD risk. Incomplete penetrance of such deleterious variation was common as over 70% of carriers had normal electrocardiographic intervals.

CONCLUSIONS

Our findings indicate that large-scale high-depth sequence data and electrocardiographic analysis identifies monogenic arrhythmia susceptibility genes and rare variants with large effects. Known pathogenic variation in conventional arrhythmia and SCD genes exhibited incomplete penetrance and accounted for only a small fraction of marked electrocardiographic interval prolongation.

A two-sample Mendelian randomization analysis of heart rate variability and cerebral small vessel disease

Cerebral small vessel disease (cSVD) is correlated with a high risk of stroke and cognitive impairment. Previous studies between heart rate variability (HRV) and cSVD revealed paradoxical results. The authors aimed to investigate the relationship between HRV and cSVD using Mendelian randomization analysis. Genetic instruments for HRV were obtained from previous genome‐wide association studies. They applied inverse variance‐weighted analysis, weighted median analysis, simple median analysis, and Mendelian randomization–Egger regression to evaluate the associations of HRV with white matter hyperintensity (WMH) and small vessel stroke (SVS) in the UK Biobank neuroimaging dataset and the MEGASTROKE genome‐wide association study dataset. Two genetically predicted traits of HRV (the root mean square of the successive differences of inter beat intervals [RMSSD] and the peak‐valley respiratory sinus arrhythmia or high frequency power [pvRSA/HF]) were suggestively associated with WMH (β 0.26, 95% confidence interval [CI] 0.04–0.49, p = .02; β 0.14, 95% CI 0.02–0.27, p = .03, respectively). Genetically predicted traits of HRV were not significantly associated with SVS. This study provides genetic support for a suggestive causal effect of HRV (RMSSD, pvRSA/HF) on WMH but not SVS.

Genetic associations for two biological age measures point to distinct aging phenotypes

Biological age measures outperform chronological age in predicting various aging outcomes, yet little is known regarding genetic predisposition. We performed genome-wide association scans of two age-adjusted biological age measures (PhenoAgeAcceleration and BioAgeAcceleration), estimated from clinical biochemistry markers (Levine et al., 2018; Levine, 2013) in European-descent participants from UK Biobank. The strongest signals were found in the APOE gene, tagged by the two major protein-coding SNPs, PhenoAgeAccel-rs429358 (APOE e4 determinant) (p = 1.50 × 10-72 ); BioAgeAccel-rs7412 (APOE e2 determinant) (p = 3.16 × 10-60 ). Interestingly, we observed inverse APOE e2 and e4 associations and unique pathway enrichments when comparing the two biological age measures. Genes associated with BioAgeAccel were enriched in lipid related pathways, while genes associated with PhenoAgeAccel showed enrichment for immune system, cell function, and carbohydrate homeostasis pathways, suggesting the two measures capture different aging domains. Our study reaffirms that aging patterns are heterogeneous across individuals, and the manner in which a person ages may be partly attributed to genetic predisposition.

Application of Nanotechnology for Sensitive Detection of Low-Abundance Single-Nucleotide Variations in Genomic DNA: A Review

Single-nucleotide polymorphisms (SNPs) are the simplest and most common type of DNA variations in the human genome. This class of attractive genetic markers, along with point mutations, have been associated with the risk of developing a wide range of diseases, including cancer, cardiovascular diseases, autoimmune diseases, and neurodegenerative diseases. Several existing methods to detect SNPs and mutations in body fluids have faced limitations. Therefore, there is a need to focus on developing noninvasive future polymerase chain reaction (PCR)-free tools to detect low-abundant SNPs in such specimens. The detection of small concentrations of SNPs in the presence of a large background of wild-type genes is the biggest hurdle. Hence, the screening and detection of SNPs need efficient and straightforward strategies. Suitable amplification methods are being explored to avoid high-throughput settings and laborious efforts. Therefore, currently, DNA sensing methods are being explored for the ultrasensitive detection of SNPs based on the concept of nanotechnology. Owing to their small size and improved surface area, nanomaterials hold the extensive capacity to be used as biosensors in the genotyping and highly sensitive recognition of single-base mismatch in the presence of incomparable wild-type DNA fragments. Different nanomaterials have been combined with imaging and sensing techniques and amplification methods to facilitate the less time-consuming and easy detection of SNPs in different diseases. This review aims to highlight some of the most recent findings on the aspects of nanotechnology-based SNP sensing methods used for the specific and ultrasensitive detection of low-concentration SNPs and rare mutations.

Potential functional variants of KIAA genes are associated with breast cancer risk in a case control study

Background

KIAA genes identified in the Kazusa cDNA-sequencing project may play important roles in biological processes and are involved in carcinogenesis of many cancers. Genetic variants of KIAA genes are implicated in the abnormal expression of these genes and are linked to susceptibility of several human complex diseases.

Methods

The differentially expressed KIAA genes were screened and identified in The Cancer Genome Atlas (TCGA) database of breast cancer. A total of 48 variants located in the 28 KIAA genes were selected to investigate the associations between polymorphism and breast cancer in 1,032 cases and 1,063 cancer-free controls in a Chinese population.

Results

Two coding variants, which included a SNP rs2306369 in KIAA1109 and a SNP rs1205434 in KIAA1755, were identified to be associated with the incidences of breast cancer. Logistic regression analysis showed that the SNP rs2306369 G allele was associated with a decreased risk of breast cancer (additive model: OR =0.81, 95% CI: 0.66-0.99, P=0.038), whereas the SNP rs1205434 A allele was involved with a higher risk of breast cancer (additive model: OR =1.19, 95% CI: 1.02-1.38, P= 0.025). Further stratified analysis revealed that the SNP rs1205434 showed a significant difference for age at menarche strata (heterogeneity test P=0.009). Multiplicative interaction analysis indicated that there was positive multiplicative interaction between the SNP rs1205434 and menarche age (OR =1.09, 95% CI: 1.01-1.17, P=0.036). Additionally, expression quantitative trait loci analysis revealed that the SNP rs1205434 A allele could decrease the KIAA1755 expression in the Genotype-Tissue Expression (GTEx) database (P=0.002). The Kaplan-Meier plotter showed that breast cancer patients with high KIAA1755 expression have significantly better outcomes than those with low levels of expression (HR =0.84, 95% CI: 0.72-0.99, P=0.033).

Conclusions

The results indicate that the genetic variants (rs2306369 and rs1205434) in the coding region of KIAA1109 and KIAA1755 respectively may affect Chinese females' breast cancer susceptibility and act as potential predictive biomarkers for breast cancer.

Two-Year Responses of Heart Rate and Heart Rate Variability to First Occupational Lead Exposure

[Figure: see text].

Functional and Adaptive Significance of Promoter Mutations That Affect Divergent Myocardial Expressions of TRIM72 in Primates

Cis-regulatory elements play important roles in tissue-specific gene expression and in the evolution of various phenotypes, and mutations in promoters and enhancers may be responsible for adaptations of species to environments. TRIM72 is a highly conserved protein that is involved in energy metabolism. Its expression in the heart varies considerably in primates, with high levels of expression in Old World monkeys and near absence in hominids. Here, we combine phylogenetic hypothesis testing and experimentation to demonstrate that mutations in promoter are responsible for the differences among primate species in the heart-specific expression of TRIM72. Maximum likelihood estimates of lineage-specific substitution rates under local-clock models show that relative to the evolutionary rate of introns, the rate of promoter was accelerated by 78% in the common ancestor of Old World monkeys, suggesting a role for positive selection in the evolution of the TRIM72 promoter, possibly driven by selective pressure due to changes in cardiac physiology after species divergence. We demonstrate that mutations in the TRIM72 promoter account for the differential myocardial TRIM72 expression of the human and the rhesus macaque. Furthermore, changes in TRIM72 expression alter the expression of genes involved in oxidative phosphorylation, which in turn affects mitochondrial respiration and cardiac energy capacity. On a broader timescale, phylogenetic regression analyses of data from 29 mammalian species show that mammals with high cardiac expression of TRIM72 have high heart rate, suggesting that the expression changes of TRIM72 may be related to differences in the heart physiology of those species.

An integrated approach to understand biological stress system dysregulation across depressive and anxiety disorders

BACKGROUND

Affective disorders involve dysregulation of major biological stress systems (hypothalamic-pituitary-adrenal (HPA)-axis, immune system, autonomic nervous system (ANS)). Suchdysregulationshave rarely beensimultaneously examined across different stress systems.

METHODS

In the Netherlands Study of Depression and Anxiety (n=2789), we investigated whether current or remitted depressive and/or anxiety disorders (based on the CIDI semi-structured interview), including specific symptom profiles, were associated with separate markers and cumulative indexes of the HPA-axis (cortisol awakening response, evening cortisol, dexamethasone suppression test cortisol), immune system (C-reactive protein, interleukin-6, tumor necrosis factor-α), and ANS (heart rate, respiratory sinus arrhythmia, pre-ejection period).

RESULTS

Depressive andanxiety disorderswere significantlyassociated with changes in three biological stress systemsincluding HPA-axis hyperactivity, increased inflammatory activity, and a higher ANS tone, particularly for integrative and cumulative indexes of these stress systems (pFDR <.05) vs. controls. The strongest associations were seen with current disorders andcumulative indexes of the HPA-axis (β=.124, pFDR=.001), the immune system (β =.057, pFDR=.032), and total cumulative index across stress systems (β=.102, pFDR=.004). Atypical, energy-related depression severity was linked to immune system markers (pFDR<0.001), melancholic depression severity to HPA-axis markers (pFDR=.032), and anxiety arousal severity to both HPA-axis and immune system markers (pFDR<0.05). Findings were partially explained by poorer lifestyle, more chronic diseases,or (especially for ANS-function) antidepressant use.

LIMITATIONS

Cross-sectional analyses limit examination of temporal associations.

CONCLUSION

Patients withdepressive and anxiety disorders showed consistent dysregulation across biological stress systems, particularly for current episodes.To understand stress system functionality in affective disorders, an integrated approach capturing cumulative stress indices within and across biological stress systems is important.

Heart rate, intelligence in adolescence, and Parkinson's disease later in life

To investigate whether physical and cognitive fitness measured in late adolescence was associated with future risk of Parkinson’s disease (PD). The cohort included 1,259,485 Swedish men with physical fitness, body mass index (BMI), resting heart rate (RHR), blood pressure, intelligence quotient (IQ), and stress resilience measured at the age of 17–20 in relation to conscription. Incident cases of PD were ascertained from the Swedish Patient Register. Hazard ratios were estimated from Cox models, after controlling for multiple confounders. We further performed Mendelian randomization (MR) analyses to assess the causality of the associations, using GWAS summary statistics with > 800,000 individuals. During follow-up, we identified 1,034 cases of PD (mean age at diagnosis = 53). Men with an RHR > 100 beats per minute had a higher risk of PD compared to men with an RHR of 60–100 beats per minute (HR = 1.47; 95% CI = 1.08–1.99). Men with IQ above the highest tertile had a higher risk of PD compared to men with an IQ below the lowest tertile (HR = 1.46; 95% CI = 1.19–1.79). We found no association for physical fitness, BMI, blood pressure, or stress resilience. A causal relationship was suggested by the MR analysis between IQ and PD, but not between RHR and PD. RHR and IQ in late adolescence were associated with a higher risk of PD diagnosed at relatively young age. The association of IQ with PD is likely causal, whereas the association of RHR with PD suggests that altered cardiac autonomic function might start before 20 years of age in PD.

Coronary artery disease: differential expression of ceRNAs and interaction analyses

Background

Previous studies have demonstrated associations between cardiovascular disease and the expression of various messenger RNAs (mRNAs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs). This study aimed to investigate the differential expression of mRNAs, lncRNAs, and miRNAs between tissues from patients with coronary artery disease (CAD) and healthy controls, and to determine the interactions between these molecules in CAD.

Methods

We investigated the differential expression of competitive endogenous RNAs (ceRNAs) between patients with CAD and healthy controls by collecting data from Gene Expression Omnibus (GEO) microarrays. We also investigated the biological function of these differentially expressed ceRNAs by performing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. We then created a protein-protein interaction (PPI) network to identify the hub genes. Biosystems and literature searches were also carried out to identify relevant signaling pathways and the potential function of the differentially expressed ceRNAs.

Results

We identified 456 expression profiles for miRNAs, 16,325 mRNA expression profiles, and 2,869 lncRNA expression profiles. With regards to connectivity, GO and KEGG analyses (count ≥9) identified the top 11 PPI network nodes in rank order. We also identified the top 15 significant nodes for the ceRNAs identified according to degree centrality (DC) (P<0.05). Collectively, our analyses confirmed that the differential expression of certain ceRNAs, and their respective signaling pathways were associated with CAD.

Conclusions

Data arising from 11 GO and KEGG pathways, the top 15 PPI network nodes with the best connectivity rank, and the top 15 ceRNA network nodes, as determined by DC rank in CAD population, indicated that the differential expression of these ceRNAs plays a key role in the CAD. Our findings highlight new molecular mechanisms for CAD and provide new options for the development of therapeutic targets.

Effects of short-term aerobic training versus CPAP therapy on heart rate variability in moderate to severe OSA patients

We compared the effects of 2-month CPAP or exercise training (ET) therapies on the autonomic balance in moderate to severe obstructive sleep apnea (OSA) through heart rate variability (HRV) analysis. Thirty-nine OSA patients were divided into CPAP (n = 18) and ET (n = 21) groups, being further split into hypertensive and non-hypertensive subgroups. All patients were submitted to continuous ECG recordings for HRV analysis. Hemodynamic parameters were recorded by oscillometry. Excessive daytime sleepiness and sleep quality were assessed through the Epworth Sleepiness Scale and the Pittsburgh questionnaire, respectively. ET decreased systolic arterial pressure in hypertensive and non-hypertensive participants when compared to baseline values, whereas diastolic arterial pressure was decreased only in non-hypertensive ones. CPAP had no effect over hemodynamic parameters in either subgroup. ET significantly increased the HRV parameters SDNN and pNN50 in non-hypertensive participants, while reducing the LF/HF ratio in both subgroups. CPAP significantly decreased SDNN in both subgroups. ET significantly decreased excessive daytime sleepiness in both subgroups, but did not affect sleep quality. CPAP significantly improved sleep quality in both subgroups, although global scores were still those of poor sleepers, while excessive daytime sleepiness was normalized only in hypertensive patients. In conclusion, while short-term ET modulated different HRV parameters, leading to a predominant vagal tone in the cardiac sympathovagal balance and decreasing blood pressure in moderate to severe OSA, short-term CPAP had next to no effect in these parameters. We believe ET should be considered as an adjunct interventional strategy in the conservative management of hypertensive or non-hypertensive OSA patients.

GWASinspector: comprehensive quality control of genome-wide association study results

Summary

Quality control (QC) of genome wide association study (GWAS) result files has become increasingly difficult due to advances in genomic technology. The main challenges include continuous increases in the number of polymorphic genetic variants contained in recent GWASs and reference panels, the rising number of cohorts participating in a GWAS consortium, and inclusion of new variant types. Here, we present GWASinspector, a flexible R package for comprehensive QC of GWAS results. This package is compatible with recent imputation reference panels, handles insertion/deletion and multi-allelic variants, provides extensive QC reports and efficiently processes big data files. Reference panels covering three human genome builds (NCBI36, GRCh37 and GRCh38) are available. GWASinspector has a user friendly design and allows easy set-up of the QC pipeline through a configuration file. In addition to checking and reporting on individual files, it can be used in preparation of a meta-analysis by testing for systemic differences between studies and generating cleaned, harmonized GWAS files. Comparison with existing GWAS QC tools shows that the main advantages of GWASinspector are its ability to more effectively deal with insertion/deletion and multi-allelic variants and its relatively low memory use.

Availability and implementation

Our package is available at The Comprehensive R Archive Network (CRAN): https://CRAN.R-project.org/package=GWASinspector. Reference datasets and a detailed tutorial can be found at the package website at http://gwasinspector.com/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Angry in America: Psychophysiological Responses to Unfair Treatment

Background

African Americans have the highest rates of hypertension-related disease of any ethnic group in the USA. Importantly, racism and discrimination have been linked to these higher rates of morbidity and mortality. Discrimination is deleterious not only to those that are the recipients of this unfair treatment but also to the partners and family members of those affected as well to those that perpetrate this bias.

Purpose

In this paper, we identify a unique pattern of physiological response to unfair treatment, we have called the “cardiovascular conundrum.” This pattern is characterized by greater heart rate variability and greater total peripheral resistance in African Americans compared to their European American counterparts.

Methods and Results

We review the evidence supporting the existence of this pattern and propose several physiological and psychological factors that might underpin it. We also propose a number of factors that might help to mitigate the deleterious effects associated with it.

Conclusions

Whereas the context of the current review is on Black/White disparities the framework we propose may be relevant to others exposed to unfair treatment. Ultimately, the systemic factors that perpetuate these inequalities will require that we first acknowledge and then face the challenges they present if we are to address the wealth and health disparities in our country.