Genetic study links components of the autonomous nervous system to heart-rate profile during exercise

Fine mapping of candidate effector genes for heart rate

An elevated resting heart rate (RHR) is associated with increased cardiovascular mortality. Genome-wide association studies (GWAS) have identified > 350 loci. Uniquely, in this study we applied genetic fine-mapping leveraging tissue specific chromatin segmentation and colocalization analyses to identify causal variants and candidate effector genes for RHR. We used RHR GWAS summary statistics from 388,237 individuals of European ancestry from UK Biobank and performed fine mapping using publicly available genomic annotation datasets. High-confidence causal variants (accounting for > 75% posterior probability) were identified, and we collated candidate effector genes using a multi-omics approach that combined evidence from colocalisation with molecular quantitative trait loci (QTLs), and long-range chromatin interaction analyses. Finally, we performed druggability analyses to investigate drug repurposing opportunities. The fine mapping pipeline indicated 442 distinct RHR signals. For 90 signals, a single variant was identified as a high-confidence causal variant, of which 22 were annotated as missense. In trait-relevant tissues, 39 signals colocalised with cis-expression QTLs (eQTLs), 3 with cis-protein QTLs (pQTLs), and 75 had promoter interactions via Hi-C. In total, 262 candidate genes were highlighted (79% had promoter interactions, 15% had a colocalised eQTL, 8% had a missense variant and 1% had a colocalised pQTL), and, for the first time, enrichment in nervous system pathways. Druggability analyses highlighted ACHE, CALCRL, MYT1 and TDP1 as potential targets. Our genetic fine-mapping pipeline prioritised 262 candidate genes for RHR that warrant further investigation in functional studies, and we provide potential therapeutic targets to reduce RHR and cardiovascular mortality.

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Causal relationship between sex hormones and cutaneous melanoma: a two-sample Mendelian randomized study

OBJECTIVE

This study aimed to elucidate the genetic aspects of the relationship between sex hormones and cutaneous melanoma risk, providing valuable insights into this complex association.

METHODS

In this study, we used estradiol, bioavailable testosterone, sex hormone-binding globulin, and total testosterone as the exposure and melanoma as the outcome for two-sample Mendelian randomization analysis. In this study, a random-effects inverse-variance weighting (IVW) model was used as the main analysis model, and the corresponding weighted median, simple mode, weighted mode, and Mendelian randomization‒Egger methods were used as supplementary methods. We assessed both heterogeneity and horizontal pleiotropy in our study, scrutinizing whether the analysis results were affected by any individual single nucleotide polymorphism.

RESULTS

The random-effects IVW method indicated that estradiol [odds ratio (OR), 1.000; 95% confidence interval (CI), 0.998-1.003; P = 0.658], bioavailable testosterone (OR = 1.001, 95% CI, 0.999-1.003; P = 0.294), sex hormone-binding globulin (IVW: OR, 1.000; 95% CI, 0.998-1.003; P = 0.658), and total testosterone (IVW: OR, 1.002; 95% CI, 0.999-1.005; P = 0.135) were not genetically linked to cutaneous melanoma. No analyses exhibited heterogeneity, horizontal pleiotropy, or deviations.

CONCLUSION

We were unable to find genetic evidence for a causal relationship between sex hormones and the occurrence of cutaneous melanoma in this study. These results are limited by sample size and population, so the causal relationship between sex hormones and cutaneous melanoma needs to be further studied.

Genetic causal relationship between gut microbiota and basal cell carcinoma: A two-sample mendelian randomization study

OBJECTIVE

Research has previously established connections between the intestinal microbiome and the progression of some cancers. However, there is a noticeable gap in the literature in regard to using Mendelian randomisation (MR) to delve into potential causal relationships between the gut microbiota (GM) and basal cell carcinoma (BCC). Therefore, the purpose of our study was to use MR to explore the causal relationship between four kinds of GM (Bacteroides, Streptococcus, Proteobacteria and Lachnospiraceae) and BCC.

METHODS

We used genome-wide association study (GWAS) data and MR to explore the causal relationship between four kinds of GM and BCC. This study primarily employed the random effect inverse variance weighted (IVW) model for analysis, as complemented by additional methods including the simple mode, weighted median, weighted mode and MR‒Egger methods. We used heterogeneity and horizontal multiplicity to judge the reliability of each analysis. MR-PRESSO was mainly used to detect and correct outliers.

RESULTS

The random-effects IVW results showed that Bacteroides (OR = 0.936, 95% CI = 0.787-1.113, p = 0.455), Streptococcus (OR = 0.974, 95% CI = 0.875-1.083, p = 0.629), Proteobacteria (OR = 1.113, 95% CI = 0.977-1.267, p = 0.106) and Lachnospiraceae (OR = 1.027, 95% CI = 0.899-1.173, p = 0.688) had no genetic causal relationship with BCC. All analyses revealed no horizontal pleiotropy, heterogeneity or outliers.

CONCLUSION

We found that Bacteroides, Streptococcus, Proteobacteria and Lachnospiraceae do not increase the incidence of BCC at the genetic level, which provides new insight for the study of GM and BCC.

Epistasis regulates genetic control of cardiac hypertrophy

The combinatorial effect of genetic variants is often assumed to be additive. Although genetic variation can clearly interact non-additively, methods to uncover epistatic relationships remain in their infancy. We develop low-signal signed iterative random forests to elucidate the complex genetic architecture of cardiac hypertrophy. We derive deep learning-based estimates of left ventricular mass from the cardiac MRI scans of 29,661 individuals enrolled in the UK Biobank. We report epistatic genetic variation including variants close to CCDC141 , IGF1R , TTN , and TNKS. Several loci where variants were deemed insignificant in univariate genome-wide association analyses are identified. Functional genomic and integrative enrichment analyses reveal a complex gene regulatory network in which genes mapped from these loci share biological processes and myogenic regulatory factors. Through a network analysis of transcriptomic data from 313 explanted human hearts, we found strong gene co-expression correlations between these statistical epistasis contributors in healthy hearts and a significant connectivity decrease in failing hearts. We assess causality of epistatic effects via RNA silencing of gene-gene interactions in human induced pluripotent stem cell-derived cardiomyocytes. Finally, single-cell morphology analysis using a novel high-throughput microfluidic system shows that cardiomyocyte hypertrophy is non-additively modifiable by specific pairwise interactions between CCDC141 and both TTN and IGF1R . Our results expand the scope of genetic regulation of cardiac structure to epistasis.

Mitral valve transcriptome analysis in thirty-four age-matched Cavalier King Charles Spaniels with or without congestive heart failure caused by myxomatous mitral valve disease

We here report the results of a mitral valve transcriptome study designed to identify genes and molecular pathways involved in development of congestive heart failure (CHF) following myxomatous mitral valve disease (MMVD) in dogs. The study is focused on a cohort of elderly age-matched dogs (n = 34, age ~ 10 years) from a single breed-Cavalier King Charles Spaniels (CKCS)-with a high incidence of MMVD. The cohort comprises 19 dogs (10♀, 9♂) without MMVD-associated CHF, and 15 dogs (6♀, 9♂) with CHF caused by MMVD; i.e., we compare gene expression in breed and age-matched groups of dogs, which only differ with respect to CHF status. We identify 56 genes, which are differentially expressed between the two groups. In this list of genes, we confirm an enrichment of genes related to the TNFβ-signaling pathway, extracellular matrix organization, vascular development, and endothelium damage, which also have been identified in previous studies. However, the genes with the greatest difference in expression between the two groups are CNTN3 and MYH1. Both genes encode proteins, which are predicted to have an effect on the contractile activity of myocardial cells, which in turn may have an effect on valvular performance and hemodynamics across the mitral valve. This may result in shear forces with impact on MMVD progression.

Genetic architecture of cardiac dynamic flow volumes

Cardiac blood flow is a critical determinant of human health. However, the definition of its genetic architecture is limited by the technical challenge of capturing dynamic flow volumes from cardiac imaging at scale. We present DeepFlow, a deep-learning system to extract cardiac flow and volumes from phase-contrast cardiac magnetic resonance imaging. A mixed-linear model applied to 37,653 individuals from the UK Biobank reveals genome-wide significant associations across cardiac dynamic flow volumes spanning from aortic forward velocity to aortic regurgitation fraction. Mendelian randomization reveals a causal role for aortic root size in aortic valve regurgitation. Among the most significant contributing variants, localizing genes (near ELN, PRDM6 and ADAMTS7) are implicated in connective tissue and blood pressure pathways. Here we show that DeepFlow cardiac flow phenotyping at scale, combined with genotyping data, reinforces the contribution of connective tissue genes, blood pressure and root size to aortic valve function.

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.

Causal relationship between psychiatric traits and temporomandibular disorders: a bidirectional two-sample Mendelian randomization study

OBJECTIVES

This study was to investigate the causal relationship between temporomandibular disorders (TMD) and psychiatric disorders by Mendelian randomization (MR) analysis.

MATERIALS AND METHODS

A two-sample bidirectional MR analysis was adopted to systematically explore the causal relationship between TMD and eight psychiatric traits, including anxiety disorder (AD), panic disorder (PD), major depressive disorder (MDD), neuroticism, attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BIP), and schizophrenia (SCZ). Inverse variance weighted (IVW), weighted median, and MR-Egger regression were used in my study. Furthermore, we also performed three sensitivity analyses to illustrate the reliability of the analysis.

RESULTS

Two psychiatric traits have risk effects on TMD: PD (OR = 1.118, 95% CI: 1.047-1.194, P = 8.161 × 10-4, MDD (OR = 1.961, 95% CI: 1.450-2.653, P = 1.230 × 10-5). Despite not surpassing the strict Bonferroni correction applied (P > 0.00625), we could think that there was a suggestive causal effect of neuroticism and SCZ increasing the risk of TMD. On the reverse MR analysis, we found no significant evidence of causal effects of TMD on these psychiatric traits. Except for heterogeneity in the causal analysis for SCZ on TMD, no heterogeneity and horizontal pleiotropy were detected in the other analyses.

CONCLUSIONS

Our two-sample MR study has provided further evidence of PD and MDD being related to a higher risk of TMD.

CLINICAL RELEVANCE

These findings highlight the importance of closely monitoring mental traits during future TMD treatments to prevent an increased risk of TMD.

Heart rate recovery after orthostatic challenge and cardiopulmonary exercise testing in older individuals: prospective multicentre observational cohort study

Background

Impaired vagal function in older individuals, quantified by the 'gold standard' delayed heart rate recovery after maximal exercise (HRRexercise), is an independent predictor of cardiorespiratory capacity and mortality (particularly when HRR ≤12 beats min-1). Heart rate also often declines after orthostatic challenge (HRRorthostatic), but the mechanism remains unclear. We tested whether HRRorthostatic reflects similar vagal autonomic characteristics as HRRexercise.

Methods

Prospective multicentre cohort study of subjects scheduled for cardiopulmonary exercise testing (CPET) as part of routine care. Before undergoing CPET, heart rate was measured with participants seated for 3 min, before standing for 3 min (HRRorthostatic). HRRexercise 1 min after the end of CPET was recorded. The primary outcome was the correlation between mean heart rate change every 10 s for 1 min after peak heart rate was attained on standing and after exercise for each participant. Secondary outcomes were HRRorthostatic and peak VO2 compared between individuals with HRRexercise <12 beats min-1.

Results

A total of 87 participants (mean age: 64 yr [95%CI: 61-66]; 48 (55%) females) completed both tests. Mean heart rate change every 10 s for 1 min after peak heart rate after standing and exercise was significantly correlated (R2=0.81; P<0.0001). HRRorthostatic was unchanged in individuals with HRRexercise ≤12 beats min-1 (n=27), but was lower when HRRexercise >12 beats min-1 (n=60; mean difference: 3 beats min-1 [95% confidence interval 1-5 beats min-1]; P<0.0001). Slower HRRorthostatic was associated with lower peak VO2 (mean difference: 3.7 ml kg-1 min-1 [95% confidence interval 0.7-6.8 ml kg-1 min-1]; P=0.039).

Conclusion

Prognostically significant heart rate recovery after exhaustive exercise is characterised by quantitative differences in heart rate recovery after orthostatic challenge. These data suggest that orthostatic challenge is a valid, simple test indicating vagal impairment.

Clinical trial registration

researchregistry6550.

Analysis of Zinc and Stromal Immunity in Disuse Osteoporosis: Mendelian Randomization and Transcriptomic Analysis

OBJECTIVE

Disuse osteoporosis is known to be primarily caused by a lack of exercise. However, the causal relationships between zinc and immunity and disuse osteoporosis remain unknown. This study investigated these relationships and their potential mechanisms.

METHODS

This study was an integrative study combining genome-wide association studies and transcriptomics. Two-sample Mendelian randomization analysis (MR) was used to analyze the causal relationships between exposures (zinc, immunity, physical activity) and the outcome (osteoporosis) with the aid of single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs). Four models, MR-Egger, inverse variance weighted, weighted median and MR-Pleiotrophy RESidual Sum and Outlier (MRPRESSO), were used to calculate odds ratio values. Sensitivity and heterogeneity analyses were also performed using MRPRESSO and MR-Egger methods. The mRNA transcriptomic analysis was subsequently conducted. Zinc metabolism scores were acquired through single-sample Gene Set Enrichment Analysis algorithms. Stromal scores were obtained using the R Package "estimate" algorithms. Important Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathways were also derived through gene set variation analysis. Cytoscape software helped construct the transcription factor (TF)-mRNA-microRNA (miRNA) network. Virtual screening and molecular docking were performed. Polymerase chain reaction validation was also carried out in vivo.

RESULTS

Causal relationships were demonstrated between zinc and exercise (95% confidence interval [CI] = 1.30-2.95, p = 0.001), exercise and immunity (95% CI = 0.36-0.80, p = 0.002), exercise and osteoporosis (95% CI = 0.97-0.99, p = 0.0007), and immunity disorder and osteoporosis (95% CI = 1.30-2.03, p = 0.00002). One hundred and seventy-nine mRNAs in important modules were screened. Combining the differential expressional genes (DEGs) and the Boruta selection, six DEGs were screened (AHNAK, CSF2, ADAMTS12, SRA1, RUNX2, and SLC39A14). TF HOXC10 and miRNA hsa-miR-204 were predicted. Then, the TF-mRNA-miRNA network was successfully constructed. RUNX2 and SLC39A14 were identified as hub mRNAs in the TF-mRNA-miRNA network. Eventually, the novel small drug C6O4NH5 was designed according to the pharmacophore structure of SLC39A14. The docking energy for the novel drug was -5.83 kcal/mol. SLC39A14 and RUNX2 were downregulated (of statistical significance p-value < 0.05) in our animal experiment.

CONCLUSION

This study revealed that zinc had a protective causal relationship with disuse osteoporosis by promoting exercise and immunity. SLC39A14 and RUNX2 mRNA participated in this zinc-related mechanism.

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.

Genetic insights into resting heart rate and its role in cardiovascular disease

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Causal relationship between osteoarthritis with atrial fibrillation and coronary atherosclerosis: a bidirectional Mendelian randomization study of European ancestry

Background

Osteoarthritis (OA) is a degenerative disease with high prevalence. Some observational studies have shown that patients with osteoarthritis often have co-existing cardiovascular diseases (CVD) such as atrial fibrillation (AF) and coronary atherosclerosis (CA). However, there is still a lack of stronger evidence confirming the association between osteoarthritis and cardiovascular disease. In this study, we used a bidirectional two-sample Mendelian randomization study to investigate the relationship between OA with AF and CA.

Methods

OA data from the UK Biobank and arcOGEN (Arthritis Research UK Osteoarthritis Genetics, a study that aimed to find genetic determinants of osteoarthritis and elucidate the genetic architecture of the disease) integration were selected for the study (n = 417,596), AF data were obtained from six studies (n = 1,030,836), and coronary atherosclerosis data were derived from the FinnGen (n = 218,792). MR analysis was performed primarily using the Inverse variance weighted (IVW) method, with MR Egger, weighted median, simple mode, weighted mode as supplements, sensitivity analysis was performed using Cochran Q statistic, and leave-one-out analysis.

Results

We found that OA and AF were positively associated [IVW: OR (95% CI): 1.11 (1.04, 1.19), P = 0.002], while OA and CA were negatively associated [IVW: OR (95% CI): 0.88 (0.79, 0.98), P = 0.02]. In the reverse MR analysis, no effect of AF on OA was found [IVW: OR (95% CI): 1.00 (0.97, 1.03), P = 0.84], meanwhile, CA and OA were found to be associated negatively [IVW: OR (95% CI): 0.95 (0.92, 0.99), P = 0.01]. No violations of MR assumptions were found in the sensitivity analysis.

Conclusion

This research confirms that OA is a risk factor for AF, and there is a mutual protective factor between OA and CA. However, further studies are still necessary to elucidate the underlying mechanisms.

Spatially resolved multiomics of human cardiac niches

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

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.

Observational and genetic evidence support a relationship between cardiac autonomic function and blood pressure

Background

It is unclear how cardiac autonomic function, as indicated by heart rate (HR), heart rate variability (HRV), HR increase during exercise, and HR recovery after exercise, is related to blood pressure (BP). We aimed to examine the observational and genetic evidence for a potential causal effect of these HR(V) traits on BP.

Methods

We performed multivariable adjusted linear regression using Lifelines and UK Biobank cohorts to investigate the relationship between HR(V) traits and BP. Linkage disequilibrium score regression was conducted to examine genetic correlations. We used two-sample Mendelian randomization (2SMR) to examine potential causal relations between HR(V) traits and BP.

Results

Observational analyses showed negative associations of all HR(V) traits with BP, except for HR, which was positively associated. Genetic correlations were directionally consistent with the observational associations, but most significant genetic correlations between HR(V) traits and BP were limited to diastolic blood pressure (DBP). 2SMR analyses suggested a potentially causal relationship between HR(V) traits and DBP but not systolic blood pressure (SBP). No reverse effect of BP on HR(V) traits was found. One standard deviation (SD) unit increase in HR was associated with a 1.82 mmHg elevation of DBP. In contrast, one ln(ms) unit increase of the root mean square of the successive differences (RMSSD) and corrected RMSSD (RMSSDc), decreased DBP by 1.79 and 1.83 mmHg, respectively. For HR increase and HR recovery at 50 s, every additional SD increase was associated with a lower DBP by 2.05 and 1.47 mmHg, respectively. Results of secondary analyses with pulse pressure as outcome were inconsistent between observational and 2SMR analyses, as well as between HR(V) traits, and therefore inconclusive.

Conclusion

Both observational and genetic evidence show strong associations between indices of cardiac autonomic function and DBP, suggesting that a larger relative contribution of the sympathetic versus the parasympathetic nervous system to cardiac function may cause elevated DBP.

Multi-omic molecular profiling and network biology for precision anaesthesiology: a narrative review

Technological advancement, data democratisation, and decreasing costs have led to a revolution in molecular biology in which the entire set of DNA, RNA, proteins, and various other molecules - the 'multi-omic' profile - can be measured in humans. Sequencing 1 million bases of human DNA now costs US$0.01, and emerging technologies soon promise to reduce the cost of sequencing the whole genome to US$100. These trends have made it feasible to sample the multi-omic profile of millions of people, much of which is publicly available for medical research. Can anaesthesiologists use these data to improve patient care? This narrative review brings together a rapidly growing literature in multi-omic profiling across numerous fields that points to the future of precision anaesthesiology. Here, we discuss how DNA, RNA, proteins, and other molecules interact in molecular networks that can be used for preoperative risk stratification, intraoperative optimisation, and postoperative monitoring. This literature provides evidence for four fundamental insights: (1) Clinically similar patients have different molecular profiles and, as a consequence, different outcomes. (2) Vast, publicly available, and rapidly growing molecular datasets have been generated in chronic disease patients and can be repurposed to estimate perioperative risk. (3) Multi-omic networks are altered in the perioperative period and influence postoperative outcomes. (4) Multi-omic networks can serve as empirical, molecular measurements of a successful postoperative course. With this burgeoning universe of molecular data, the anaesthesiologist-of-the-future will tailor their clinical management to an individual's multi-omic profile to optimise postoperative outcomes and long-term health.

Genetic variants associated with syncope implicate neural and autonomic processes

AIMS

Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications.

METHODS AND RESULTS

This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders.

CONCLUSION

The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

A Mendelian randomization study to assess the genetic liability of gastroesophageal reflux disease for cardiovascular diseases and risk factors

Observational studies have reported that gastroesophageal reflux disease (GERD) is a risk factor for cardiovascular diseases (CVD); however, the causal inferences between them remain unknown. We conducted a Mendelian randomization (MR) study to estimate the causal associations between GERD and 10 CVD outcomes, as well as 14 cardiovascular risk factors. We used summary statistics from genome-wide association studies for GERD and the FinnGen consortium for CVD. We further investigated whether GERD correlated with cardiovascular risk factors and performed multivariable MR and mediation analyses to estimate the mediating effects of these risk factors on GERD-CVD progression. Sensitivity analyses and replication analyses were also performed. Our results indicated that GERD was positively associated with seven CVD outcomes with odds ratios of 1.26 [95% confidence interval (CI), 1.15, 1.37] for coronary artery disease, 1.41 (95% CI, 1.28, 1.57) for myocardial infarction, 1.34 (95% CI, 1.19, 1.51) for atrial fibrillation, 1.34 (95% CI, 1.21, 1.50) for heart failure, 1.30 (95% CI, 1.18, 1.43) for any stroke, 1.19 (95% CI, 1.06, 1.34) for ischemic stroke and 1.29 (95% CI, 1.16, 1.44) for venous thromboembolism. Furthermore, GERD was associated with nine cardiovascular risk factors and major depressive disorder demonstrated significant mediation effects on the causal pathway linking GERD and any stroke. This study demonstrates that GERD is associated with seven CVD outcomes and nine cardiovascular risk factors. Importantly, GERD treatment may help prevent common CVD events.

Genome-Wide Association Analysis of Over 170,000 Individuals from the UK Biobank Identifies Seven Loci Associated with Dietary Approaches to Stop Hypertension (DASH) Diet

Diet is a modifiable risk factor for common chronic diseases and mental health disorders, and its effects are under partial genetic control. To estimate the impact of diet on individual health, most epidemiological and genetic studies have focused on individual aspects of dietary intake. However, analysing individual food groups in isolation does not capture the complexity of the whole diet pattern. Dietary indices enable a holistic estimation of diet and account for the intercorrelations between food and nutrients. In this study we performed the first ever genome-wide association study (GWA) including 173,701 individuals from the UK Biobank to identify genetic variants associated with the Dietary Approaches to Stop Hypertension (DASH) diet. DASH was calculated using the 24 h-recall questionnaire collected by UK Biobank. The GWA was performed using a linear mixed model implemented in BOLT-LMM. We identified seven independent single-nucleotide polymorphisms (SNPs) associated with DASH. Significant genetic correlations were observed between DASH and several educational traits with a significant enrichment for genes involved in the AMP-dependent protein kinase (AMPK) activation that controls the appetite by regulating the signalling in the hypothalamus. The colocalization analysis implicates genes involved in body mass index (BMI)/obesity and neuroticism (ARPP21, RP11-62H7.2, MFHAS1, RHEBL1). The Mendelian randomisation analysis suggested that increased DASH score, which reflect a healthy diet style, is causal of lower glucose, and insulin levels. These findings further our knowledge of the pathways underlying the relationship between diet and health outcomes. They may have significant implications for global public health and provide future dietary recommendations for the prevention of common chronic diseases.

Integrated Analysis of the microRNA–mRNA Network Predicts Potential Regulators of Atrial Fibrillation in Humans

Atrial fibrillation (AF) is a form of sustained cardiac arrhythmia and microRNAs (miRs) play crucial roles in the pathophysiology of AF. To identify novel miR–mRNA pairs, we performed RNA-seq from atrial biopsies of persistent AF patients and non-AF patients with normal sinus rhythm (SR). Differentially expressed miRs (11 down and 9 up) and mRNAs (95 up and 82 down) were identified and hierarchically clustered in a heat map. Subsequently, GO, KEGG, and GSEA analyses were run to identify deregulated pathways. Then, miR targets were predicted in the miRDB database, and a regulatory network of negatively correlated miR–mRNA pairs was constructed using Cytoscape. To select potential candidate genes from GSEA analysis, the top-50 enriched genes in GSEA were overlaid with predicted targets of differentially deregulated miRs. Further, the protein–protein interaction (PPI) network of enriched genes in GSEA was constructed, and subsequently, GO and canonical pathway analyses were run for genes in the PPI network. Our analyses showed that TNF-α, p53, EMT, and SYDECAN1 signaling were among the highly affected pathways in AF samples. SDC-1 (SYNDECAN-1) was the top-enriched gene in p53, EMT, and SYDECAN1 signaling. Consistently, SDC-1 mRNA and protein levels were significantly higher in atrial samples of AF patients. Among negatively correlated miRs, miR-302b-3p was experimentally validated to suppress SDC-1 transcript levels. Overall, our results suggested that the miR-302b-3p/SDC-1 axis may be involved in the pathogenesis of AF.

Deep learning on resting electrocardiogram to identify impaired heart rate recovery

Background and Objective

Postexercise heart rate recovery (HRR) is an important indicator of cardiac autonomic function and abnormal HRR is associated with adverse outcomes. We hypothesized that deep learning on resting electrocardiogram (ECG) tracings may identify individuals with impaired HRR.

Methods

We trained a deep learning model (convolutional neural network) to infer HRR based on resting ECG waveforms (HRR_pred) among UK Biobank participants who had undergone exercise testing. We examined the association of HRR_pred with incident cardiovascular disease using Cox models, and investigated the genetic architecture of HRR_pred in genome-wide association analysis.

Results

Among 56,793 individuals (mean age 57 years, 51% women), the HRR_pred model was moderately correlated with actual HRR (r = 0.48, 95% confidence interval [CI] 0.47-0.48). Over a median follow-up of 10 years, we observed 2060 incident diabetes mellitus (DM) events, 862 heart failure events, and 2065 deaths. Higher HRR_pred was associated with lower risk of DM (hazard ratio [HR] 0.79 per 1 standard deviation change, 95% CI 0.76-0.83), heart failure (HR 0.89, 95% CI 0.83-0.95), and death (HR 0.83, 95% CI 0.79-0.86). After accounting for resting heart rate, the association of HRR_pred with incident DM and all-cause mortality were similar. Genetic determinants of HRR_pred included known heart rate, cardiac conduction system, cardiomyopathy, and metabolic trait loci.

Conclusion

Deep learning-derived estimates of HRR using resting ECG independently associated with future clinical outcomes, including new-onset DM and all-cause mortality. Inferring postexercise heart rate response from a resting ECG may have potential clinical implications and impact on preventive strategies warrants future study.

The Omnipresence of Autonomic Modulation in Health and Disease

The Autonomic Nervous System (ANS) is a critical part of the homeostatic machinery with both central and peripheral components. However, little is known about the integration of these components and their joint role in the maintenance of health and in allostatic derailments leading to somatic and/or neuropsychiatric (co)morbidity. Based on a comprehensive literature search on the ANS neuroanatomy we dissect the complex integration of the ANS: (1) First we summarize Stress and Homeostatic Equilibrium - elucidating the responsivity of the ANS to stressors; (2) Second we describe the overall process of how the ANS is involved in Adaptation and Maladaptation to Stress; (3) In the third section the ANS is hierarchically partitioned into the peripheral/spinal, brainstem, subcortical and cortical components of the nervous system. We utilize this anatomical basis to define a model of autonomic integration. (4) Finally, we deploy the model to describe human ANS involvement in (a) Hypofunctional and (b) Hyperfunctional states providing examples in the healthy state and in clinical conditions.

Causal relationship between osteoporosis and osteoarthritis: A two-sample Mendelian randomized study

INTRODUCTION

At present, clinical studies have confirmed that osteoporosis (OP) has an inverse relationship with osteoarthritis (OA), but it has not been proven from the point of view of genetics, so our study hopes to clarify the potential effect of OP on OA at the level of gene prediction through two-sample Mendelian randomization (MR) analysis.

METHODS

A two-sample MR was adopted to research the causal relationship of OP with OA (including total OA, knee OA and hip OA). All data come from a public shared database. Such traditional methods as simple and weighted models, inverse variance weighted, weighted median, and Mendelian Randomization (MR-Egger) regression were employed to assess the causal effect of OP on OA. We used the Pleiotrophy RESidual Sum and Outlier (MR-PRESSO) method and MR-Egger method to study sensitivity. The leave-one-out test is used to determine the influence of outliers. The heterogeneity was calculated by using Cochran Q statistics and MR-Egger regression in the inverse variance-weighted (IVW) method. P > 0.05 indicates that there is a large heterogeneity. MR-Robust Adjustment Profile Score (RAPS) is stable to both systematic and specific multiplicity, so we used MR-RAPS as a supplementary method to verify the results of IVW.

RESULTS

According to the results of IVW, we found that there was a causal relationship between OP and total OA, and OP reduced the incidence of total OA (beta=-0.285, OR=0.751, P value< 0.016). The MR estimation of the causal effect of OP on knee OA suggested that the genetic prediction of OP was negatively correlated with knee osteoarthritis (KOA) (IVW: beta=-6.11, OR=0.002, P value< 0.016). The IVW results suggested that OP was causally related to hip OA, and OP had a protective effect on hip OA (beta=-5.48, OR=4.15e-3, P value= 3.99e-3). Except for heterogeneity in the analysis of OP and knee OA, there was no horizontal pleiotropy or heterogeneity in the other analyses.

CONCLUSION

We explored the causal relationship between OP and OA through a two-sample MR analysis and found that OP can reduce the incidence of OA (including knee OA and hip OA).

The genetic case for cardiorespiratory fitness as a clinical vital sign and the routine prescription of physical activity in healthcare

BACKGROUND

Cardiorespiratory fitness (CRF) and physical activity (PA) are well-established predictors of morbidity and all-cause mortality. However, CRF is not routinely measured and PA not routinely prescribed as part of standard healthcare. The American Heart Association (AHA) recently presented a scientific case for the inclusion of CRF as a clinical vital sign based on epidemiological and clinical observation. Here, we leverage genetic data in the UK Biobank (UKB) to strengthen the case for CRF as a vital sign and make a case for the prescription of PA.

METHODS

We derived two CRF measures from the heart rate data collected during a submaximal cycle ramp test: CRF-vo2max, an estimate of the participants' maximum volume of oxygen uptake, per kilogram of body weight, per minute; and CRF-slope, an estimate of the rate of increase of heart rate during exercise. Average PA over a 7-day period was derived from a wrist-worn activity tracker. After quality control, 70,783 participants had data on the two derived CRF measures, and 89,683 had PA data. We performed genome-wide association study (GWAS) analyses by sex, and post-GWAS techniques to understand genetic architecture of the traits and prioritise functional genes for follow-up.

RESULTS

We found strong evidence that genetic variants associated with CRF and PA influenced genetic expression in a relatively small set of genes in the heart, artery, lung, skeletal muscle and adipose tissue. These functionally relevant genes were enriched among genes known to be associated with coronary artery disease (CAD), type 2 diabetes (T2D) and Alzheimer's disease (three of the top 10 causes of death in high-income countries) as well as Parkinson's disease, pulmonary fibrosis, and blood pressure, heart rate, and respiratory phenotypes. Genetic variation associated with lower CRF and PA was also correlated with several disease risk factors (including greater body mass index, body fat and multiple obesity phenotypes); a typical T2D profile (including higher insulin resistance, higher fasting glucose, impaired beta-cell function, hyperglycaemia, hypertriglyceridemia); increased risk for CAD and T2D; and a shorter lifespan.

CONCLUSIONS

Genetics supports three decades of evidence for the inclusion of CRF as a clinical vital sign. Given the genetic, clinical and epidemiological evidence linking CRF and PA to increased morbidity and mortality, regular measurement of CRF as a marker of health and routine prescription of PA could be a prudent strategy to support public health.

The many faces of the E3 ubiquitin ligase, RNF220, in neural development and beyond

Ubiquitin modification plays important roles in many cellular processes that are fundamental for vertebrate embryo development, such as cell division, differentiation, and migration. Aberrant function or deregulation of ubiquitination enzymes can cause developmental disorders, cancer progression, and neurodegenerative diseases in humans. RING finger protein 220 (RNF220) is an evolutionarily conserved RING-type ubiquitin E3 ligase. Recent studies have revealed the roles and mechanisms of RNF220 and its partner protein, zinc finger C4H2-type containing protein (ZC4H2), in embryonic development and human diseases. Using mouse and zebrafish models, it has been shown that RNF220 regulates sonic hedgehog (Shh) signaling via Gli and embryonic ectoderm development (EED), a polycomb repressive complex 2 (PRC2) component, during ventral neural patterning and cerebellum development. In addition, RNF220 also regulates the development and functions of central noradrenergic and motor neurons in mice. By stabilizing β-catenin and signal transducer and activator of transcription 1 (STAT1), RNF220 is also involved in Wnt and interferon (IFN)-STAT1 signaling and thus the regulation of tumorigenesis and immune response, respectively. In humans, both RNF220 and ZC4H2 mutations have been reported to be associated with diseases accompanied by complicated neural defects. In this review, we summarize the current knowledge of RNF220 with special emphasis on its roles and mechanisms of action in signal transduction, vertebrate neural development, and related human disorders.

The genetics of human performance

Human physiology is likely to have been selected for endurance physical activity. However, modern humans have become largely sedentary, with physical activity becoming a leisure-time pursuit for most. Whereas inactivity is a strong risk factor for disease, regular physical activity reduces the risk of chronic disease and mortality. Although substantial epidemiological evidence supports the beneficial effects of exercise, comparatively little is known about the molecular mechanisms through which these effects operate. Genetic and genomic analyses have identified genetic variation associated with human performance and, together with recent proteomic, metabolomic and multi-omic analyses, are beginning to elucidate the molecular genetic mechanisms underlying the beneficial effects of physical activity on human health.

Genetically Predicted Insomnia in Relation to 14 Cardiovascular Conditions and 17 Cardiometabolic Risk Factors: A Mendelian Randomization Study

Background

This Mendelian randomization study aims to investigate causal associations between genetically predicted insomnia and 14 cardiovascular diseases (CVDs) as well as the potential mediator role of 17 cardiometabolic risk factors.

Methods and Results

Using genetic association estimates from large genome‐wide association studies and UK Biobank, we performed a 2‐sample Mendelian randomization analysis to estimate the associations of insomnia with 14 CVD conditions in the primary analysis. Then mediation analysis was conducted to explore the potential mediator role of 17 cardiometabolic risk factors using a network Mendelian randomization design. After correcting for multiple testing, genetically predicted insomnia was consistent significantly positively associated with 9 of 14 CVDs, those odds ratios ranged from 1.13 (95% CI, 1.08–1.18) for atrial fibrillation to 1.24 (95% CI, 1.16–1.32) for heart failure. Moreover, genetically predicted insomnia was consistently associated with higher body mass index, triglycerides, and lower high‐density lipoprotein cholesterol, each of which may act as a mediator in the causal pathway from insomnia to several CVD outcomes. Additionally, we found very little evidence to support a causal link between insomnia with abdominal aortic aneurysm, thoracic aortic aneurysm, total cholesterol, low‐density lipoprotein cholesterol, glycemic traits, renal function, and heart rate increase during exercise. Finally, we found no evidence of causal associations of genetically predicted body mass index, high‐density lipoprotein cholesterol, or triglycerides on insomnia.

Conclusions

This study provides evidence that insomnia is associated with 9 of 14 CVD outcomes, some of which may be partially mediated by 1 or more of higher body mass index, triglycerides, and lower high‐density lipoprotein cholesterol.

Analysing electrocardiographic traits and predicting cardiac risk in UK biobank

The electrocardiogram (ECG) is a commonly used clinical tool that reflects cardiac excitability and disease. Many parameters are can be measured and with the improvement of methodology can now be quantified in an automated fashion, with accuracy and at scale. Furthermore, these measurements can be heritable and thus genome wide association studies inform the underpinning biological mechanisms. In this review we describe how we have used the resources in UK Biobank to undertake such work. In particular, we focus on a substudy uniquely describing the response to exercise performed at scale with accompanying genetic information.

Genetic insight into sick sinus syndrome

Aims

The aim of this study was to use human genetics to investigate the pathogenesis of sick sinus syndrome (SSS) and the role of risk factors in its development.

Methods and results

We performed a genome-wide association study of 6469 SSS cases and 1 000 187 controls from deCODE genetics, the Copenhagen Hospital Biobank, UK Biobank, and the HUNT study. Variants at six loci associated with SSS, a reported missense variant in MYH6, known atrial fibrillation (AF)/electrocardiogram variants at PITX2, ZFHX3, TTN/CCDC141, and SCN10A and a low-frequency (MAF = 1.1–1.8%) missense variant, p.Gly62Cys in KRT8 encoding the intermediate filament protein keratin 8. A full genotypic model best described the p.Gly62Cys association (P = 1.6 × 10⁻²⁰), with an odds ratio (OR) of 1.44 for heterozygotes and a disproportionally large OR of 13.99 for homozygotes. All the SSS variants increased the risk of pacemaker implantation. Their association with AF varied and p.Gly62Cys was the only variant not associating with any other arrhythmia or cardiovascular disease. We tested 17 exposure phenotypes in polygenic score (PGS) and Mendelian randomization analyses. Only two associated with the risk of SSS in Mendelian randomization, AF, and lower heart rate, suggesting causality. Powerful PGS analyses provided convincing evidence against causal associations for body mass index, cholesterol, triglycerides, and type 2 diabetes (P > 0.05).

Conclusion

We report the associations of variants at six loci with SSS, including a missense variant in KRT8 that confers high risk in homozygotes and points to a mechanism specific to SSS development. Mendelian randomization supports a causal role for AF in the development of SSS.

No Clinically Relevant Effect of Heart Rate Increase and Heart Rate Recovery During Exercise on Cardiovascular Disease: A Mendelian Randomization Analysis

BACKGROUND

Reduced heart rate (HR) increase (HRI), recovery (HRR), and higher resting HR are associated with cardiovascular (CV) disease, but causal inferences have not been deduced. We investigated causal effects of HRI, HRR, and resting HR on CV risk, all-cause mortality (ACM), atrial fibrillation (AF), coronary artery disease (CAD), and ischemic stroke (IS) using Mendelian Randomization.

METHODS

11 variants for HRI, 11 for HRR, and two sets of 46 and 414 variants for resting HR were obtained from four genome-wide association studies (GWASs) on UK Biobank. We performed a lookup on GWASs for CV risk and ACM in UK Biobank (N = 375,367, 5.4% cases and N = 393,165, 4.4% cases, respectively). For CAD, AF, and IS, we used publicly available summary statistics. We used a random-effects inverse-variance weighted (IVW) method and sensitivity analyses to estimate causality.

RESULTS

IVW showed a nominally significant effect of HRI on CV events (odds ratio [OR] = 1.0012, P = 4.11 × 10-2) and on CAD and AF. Regarding HRR, IVW was not significant for any outcome. The IVW method indicated statistically significant associations of resting HR with AF (OR = 0.9825, P = 9.8 × 10-6), supported by all sensitivity analyses, and a nominally significant association with IS (OR = 0.9926, P = 9.82 × 10-3).

CONCLUSION

Our findings suggest no strong evidence of an association between HRI and HRR and any outcome and confirm prior work reporting a highly significant effect of resting HR on AF. Future research is required to explore HRI and HRR associations further using more powerful predictors, when available.

Resting Heartbeat Complexity Predicts All-Cause and Cardiorespiratory Mortality in Middle- to Older-Aged Adults From the UK Biobank

Background

Spontaneous heart rate fluctuations contain rich information related to health and illness in terms of physiological complexity, an accepted indicator of plasticity and adaptability. However, it is challenging to make inferences on complexity from shorter, more practical epochs of data. Distribution entropy (DistEn) is a recently introduced complexity measure that is designed specifically for shorter duration heartbeat recordings. We hypothesized that reduced DistEn predicted increased mortality in a large population cohort.

Method and Results

The prognostic value of DistEn was examined in 7631 middle‐older–aged UK Biobank participants who had 2‐minute resting ECGs conducted (mean age, 59.5 years; 60.4% women). During a median follow‐up period of 7.8 years, 451 (5.9%) participants died. In Cox proportional hazards models with adjustment for demographics, lifestyle factors, physical activity, cardiovascular risks, and comorbidities, for each 1‐SD decrease in DistEn, the risk increased by 36%, 56%, and 73% for all‐cause, cardiovascular, and respiratory disease–related mortality, respectively. These effect sizes were equivalent to the risk of death from being >5 years older, having been a former smoker, or having diabetes mellitus. Lower DistEn was most predictive of death in those <55 years with a prior myocardial infarction, representing an additional 56% risk for mortality compared with older participants without prior myocardial infarction. These observations remained after controlling for traditional mortality predictors, resting heart rate, and heart rate variability.

Conclusions

Resting heartbeat complexity from short, resting ECGs was independently associated with mortality in middle‐ to older‐aged adults. These risks appear most pronounced in middle‐aged participants with prior MI, and may uniquely contribute to mortality risk screening.

Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer

RATIONALE

The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined.

OBJECTIVE

Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation.

METHODS AND RESULTS

We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell-derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell-specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos.

CONCLUSIONS

We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2, TBX3, and ISL1, and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3.

Adaptive selection drives TRPP3 loss-of-function in an Ethiopian population

TRPP3 (also called PKD2L1) is a nonselective, cation-permeable channel activated by multiple stimuli, including extracellular pH changes. TRPP3 had been considered a candidate for sour sensor in humans, due to its high expression in a subset of tongue receptor cells detecting sour, along with its membership to the TRP channel family known to function as sensory receptors. Here, we describe the functional consequences of two non-synonymous genetic variants (R278Q and R378W) found to be under strong positive selection in an Ethiopian population, the Gumuz. Electrophysiological studies and 3D modelling reveal TRPP3 loss-of-functions produced by both substitutions. R278Q impairs TRPP3 activation after alkalinisation by mislocation of H⁺ binding residues at the extracellular polycystin mucolipin domain. R378W dramatically reduces channel activity by altering conformation of the voltage sensor domain and hampering channel transition from closed to open state. Sour sensitivity tests in R278Q/R378W carriers argue against both any involvement of TRPP3 in sour detection and the role of such physiological process in the reported evolutionary positive selection past event.

The Genetic Makeup of the Electrocardiogram

The electrocardiogram (ECG) is one of the most useful non-invasive diagnostic tests for a wide array of cardiac disorders. Traditional approaches to analyzing ECGs focus on individual segments. Here, we performed comprehensive deep phenotyping of 77,190 ECGs in the UK Biobank across the complete cycle of cardiac conduction, resulting in 500 spatial-temporal datapoints, across 10 million genetic variants. In addition to characterizing polygenic risk scores for the traditional ECG segments, we identified over 300 genetic loci that are statistically associated with the high-dimensional representation of the ECG. We established the genetic ECG signature for dilated cardiomyopathy, associated the BAG3, HSPB7/CLCNKA, PRKCA, TMEM43, and OBSCN loci with disease risk and confirmed this association in an independent cohort. In total, our work demonstrates that a high-dimensional analysis of the entire ECG provides unique opportunities for studying cardiac biology and disease and furthering drug development. A record of this paper's transparent peer review process is included in the Supplemental Information.

Biological insights from multi-omic analysis of 31 genomic risk loci for adult hearing difficulty

Age-related hearing impairment (ARHI), one of the most common medical conditions, is strongly heritable, yet its genetic causes remain largely unknown. We conducted a meta-analysis of GWAS summary statistics from multiple hearing-related traits in the UK Biobank (n = up to 330,759) and identified 31 genome-wide significant risk loci for self-reported hearing difficulty (p < 5x10-8), of which eight have not been reported previously in the peer-reviewed literature. We investigated the regulatory and cell specific expression for these loci by generating mRNA-seq, ATAC-seq, and single-cell RNA-seq from cells in the mouse cochlea. Risk-associated genes were most strongly enriched for expression in cochlear epithelial cells, as well as for genes related to sensory perception and known Mendelian deafness genes, supporting their relevance to auditory function. Regions of the human genome homologous to open chromatin in epithelial cells from the mouse were strongly enriched for heritable risk for hearing difficulty, even after adjusting for baseline effects of evolutionary conservation and cell-type non-specific regulatory regions. Epigenomic and statistical fine-mapping most strongly supported 50 putative risk genes. Of these, 39 were expressed robustly in mouse cochlea and 16 were enriched specifically in sensory hair cells. These results reveal new risk loci and risk genes for hearing difficulty and suggest an important role for altered gene regulation in the cochlear sensory epithelium.

Thyroid Function Affects the Risk of Stroke via Atrial Fibrillation: A Mendelian Randomization Study

CONTEXT

Observational studies suggest that variations in normal range thyroid function are associated with cardiovascular diseases. However, it remains to be determined whether these associations are causal or not.

OBJECTIVE

To test whether genetically determined variation in normal range thyroid function is causally associated with the risk of stroke and coronary artery disease (CAD) and investigate via which pathways these relations may be mediated.

DESIGN, SETTING, AND PARTICIPANTS

Mendelian randomization analyses for stroke and CAD using genetic instruments associated with normal range thyrotropin (TSH) and free thyroxine levels or Hashimoto's thyroiditis and Graves' disease. The potential mediating role of known stroke and CAD risk factors was examined. Publicly available summary statistics data were used.

MAIN OUTCOME MEASURES

Stroke or CAD risk per genetically predicted increase in TSH or FT4 levels.

RESULTS

A 1 standard deviation increase in TSH was associated with a 5% decrease in the risk of stroke (odds ratio [OR], 0.95; 95% confidence interval [CI], 0.91-0.99; P = 0.008). Multivariable MR analyses indicated that this effect is mainly mediated via atrial fibrillation. MR analyses did not show a causal association between normal range thyroid function and CAD. Secondary analyses showed a causal relationship between Hashimoto's thyroiditis and a 7% increased risk of CAD (OR, 1.07; 95% CI, 1.01-1.13; P = 0.026), which was mainly mediated via body mass index.

CONCLUSION

These results provide important new insights into the causal relationships and mediating pathways between thyroid function, stroke, and CAD. We identify variation in normal range thyroid function and Hashimoto's thyroiditis as risk factors for stroke and CAD, respectively.

Translating GWAS-identified loci for cardiac rhythm and rate using an in vivo image- and CRISPR/Cas9-based approach

A meta-analysis of genome-wide association studies (GWAS) identified eight loci that are associated with heart rate variability (HRV), but candidate genes in these loci remain uncharacterized. We developed an image- and CRISPR/Cas9-based pipeline to systematically characterize candidate genes for HRV in live zebrafish embryos. Nine zebrafish orthologues of six human candidate genes were targeted simultaneously in eggs from fish that transgenically express GFP on smooth muscle cells (Tg[acta2:GFP]), to visualize the beating heart. An automated analysis of repeated 30 s recordings of beating atria in 381 live, intact zebrafish embryos at 2 and 5 days post-fertilization highlighted genes that influence HRV (hcn4 and si:dkey-65j6.2 [KIAA1755]); heart rate (rgs6 and hcn4); and the risk of sinoatrial pauses and arrests (hcn4). Exposure to 10 or 25 µM ivabradine—an open channel blocker of HCNs—for 24 h resulted in a dose-dependent higher HRV and lower heart rate at 5 days post-fertilization. Hence, our screen confirmed the role of established genes for heart rate and rhythm (RGS6 and HCN4); showed that ivabradine reduces heart rate and increases HRV in zebrafish embryos, as it does in humans; and highlighted a novel gene that plays a role in HRV (KIAA1755).

Trait-associated noncoding variant regions affect TBX3 regulation and cardiac conduction

Genome-wide association studies have implicated common genomic variants in the gene desert upstream of TBX3 in cardiac conduction velocity. Whether these noncoding variants affect expression of TBX3 or neighboring genes and how they affect cardiac conduction is not understood. Here, we use high-throughput STARR-seq to test the entire 1.3 Mb human and mouse TBX3 locus, including two cardiac conduction-associated variant regions, for regulatory function. We identified multiple accessible and functional regulatory DNA elements that harbor variants affecting their activity. Both variant regions drove gene expression in the cardiac conduction tissue in transgenic reporter mice. Genomic deletion from the mouse genome of one of the regions caused increased cardiac expression of only Tbx3, PR interval shortening and increased QRS duration. Combined, our findings address the mechanistic link between trait-associated variants in the gene desert, TBX3 regulation and cardiac conduction.

Genome-wide association studies and Mendelian randomization analyses for leisure sedentary behaviours

Leisure sedentary behaviours are associated with increased risk of cardiovascular disease, but whether this relationship is causal is unknown. The aim of this study is to identify genetic determinants associated with leisure sedentary behaviours and to estimate the potential causal effect on coronary artery disease (CAD). Genome wide association analyses of leisure television watching, leisure computer use and driving behaviour in the UK Biobank identify 145, 36 and 4 genetic loci (P < 1×10-8), respectively. High genetic correlations are observed between sedentary behaviours and neurological traits, including education and body mass index (BMI). Two-sample Mendelian randomization (MR) analysis estimates a causal effect between 1.5 hour increase in television watching and CAD (OR 1.44, 95%CI 1.25-1.66, P = 5.63 × 10-07), that is partially independent of education and BMI in multivariable MR analyses. This study finds independent observational and genetic support for the hypothesis that increased sedentary behaviour by leisure television watching is a risk factor for CAD.

White Blood Cells and Blood Pressure: A Mendelian Randomization Study

BACKGROUND

High blood pressure (BP) is a risk factor for cardiovascular morbidity and mortality. While BP is regulated by the function of kidney, vasculature, and sympathetic nervous system, recent experimental data suggest that immune cells may play a role in hypertension.

METHODS

We studied the relationship between major white blood cell types and blood pressure in the UK Biobank population and used Mendelian randomization (MR) analyses using the ≈750 000 UK-Biobank/International Consortium of Blood Pressure-Genome-Wide Association Studies to examine which leukocyte populations may be causally linked to BP.

RESULTS

A positive association between quintiles of lymphocyte, monocyte, and neutrophil counts, and increased systolic BP, diastolic BP, and pulse pressure was observed (eg, adjusted systolic BP mean±SE for 1st versus 5th quintile respectively: 140.13±0.08 versus 141.62±0.07 mm Hg for lymphocyte, 139.51±0.08 versus 141.84±0.07 mm Hg for monocyte, and 137.96±0.08 versus 142.71±0.07 mm Hg for neutrophil counts; all P<10^-50). Using 121 single nucleotide polymorphisms in MR, implemented through the inverse-variance weighted approach, we identified a potential causal relationship of lymphocyte count with systolic BP and diastolic BP (causal estimates: 0.69 [95% CI, 0.19-1.20] and 0.56 [95% CI, 0.23-0.90] of mm Hg per 1 SD genetically elevated lymphocyte count, respectively), which was directionally concordant to the observational findings. These inverse-variance weighted estimates were consistent with other robust MR methods. The exclusion of rs3184504 SNP in the SH2B3 locus attenuated the magnitude of the signal in some of the MR analyses. MR in the reverse direction found evidence of positive effects of BP indices on counts of monocytes, neutrophils, and eosinophils but not lymphocytes or basophils. Subsequent MR testing of lymphocyte count in the context of genetic correlation with renal function or resting and postexercise heart rate demonstrated a positive association of lymphocyte count with urine albumin-to-creatinine ratio.

CONCLUSIONS

Observational and genetic analyses demonstrate a concordant, positive and potentially causal relationship of lymphocyte count with systolic BP and diastolic BP.

Genetics and the heart rate response to exercise

The acute heart rate response to exercise, i.e., heart rate increase during and heart rate recovery after exercise, has often been associated with all-cause and cardiovascular mortality. The long-term response of heart rate to exercise results in favourable changes in chronotropic function, including decreased resting and submaximal heart rate as well as increased heart rate recovery. Both the acute and long-term heart rate response to exercise have been shown to be heritable. Advances in genetic analysis enable researchers to investigate this hereditary component to gain insights in possible molecular mechanisms underlying interindividual differences in the heart rate response to exercise. In this review, we comprehensively searched candidate gene, linkage, and genome-wide association studies that investigated the heart rate response to exercise. A total of ten genes were associated with the acute heart rate response to exercise in candidate gene studies. Only one gene (CHRM2), related to heart rate recovery, was replicated in recent genome-wide association studies (GWASs). Additional 17 candidate causal genes were identified for heart rate increase and 26 for heart rate recovery in these GWASs. Nine of these genes were associated with both acute increase and recovery of the heart rate during exercise. These genes can be broadly categorized into four categories: (1) development of the nervous system (CCDC141, PAX2, SOX5, and CAV2); (2) prolongation of neuronal life span (SYT10); (3) cardiac development (RNF220 and MCTP2); (4) cardiac rhythm (SCN10A and RGS6). Additional 10 genes were linked to long-term modification of the heart rate response to exercise, nine with heart rate increase and one with heart rate recovery. Follow-up will be essential to get functional insights in how candidate causal genes affect the heart rate response to exercise. Future work will be required to translate these findings to preventive and therapeutic applications.

The hypertension advantage and natural selection: Since type 2 diabetes associates with co-morbidities and premature death, why have the genetic variants remained in the human genome?

Type 2 diabetes is a major public health crisis around the world. It is estimated that more than 300 million people worldwide have type 2 diabetes. Furthermore, the World Health Organization estimates that deaths from the complications of diabetes will increase by two thirds between 2008 and 2030. Since type 2 diabetes is a major public health crisis, why have the genetic variants for diabetes not been removed from the genome by natural selection? We hypothesize that insulin resistance, a predisposition to type 2 diabetes, and the associated elevation in sympathetic nervous system activity and arterial blood pressure provided an advantage to humans who lived as hunter-gatherers. Specifically, sympathetic hyperactivity stimulates the renin-angiotensin aldosterone system, promotes sodium reabsorption, and increases blood volume, heart rate, stroke volume and peripheral vascular resistance, thus inducing hypertension. The hypertension in turn provides a hemodynamic advantage for hunter-gatherers. Specifically, sympathetic hyperactivity and increased blood pressure increases blood flow delivery to working muscles by increasing cardiac output and shunting blood from non-active tissue. This natural selection for hypertension occurred during the time in human evolutionary history when the lifespan of most individuals was probably 30-40 years, and morbidity and mortality from cardiovascular disorders was limited. Thus, the selection pressure for elevation in sympathetic nervous system activity and blood pressure provided an advantage for hunting and gathering that would be greater than the selection pressure exerted by the manifestations of cardiovascular disease in aged individuals.

Genetic and epigenetic sex-specific adaptations to endurance exercise

In recent years, the interest in personalised interventions such as medicine, nutrition, and exercise is rapidly rising to maximize health outcomes and ensure the most appropriate treatments. Exercising regularly is recommended for both healthy and diseased populations to improve health. However, there are sex-specific adaptations to exercise that often are not taken into consideration. While endurance exercise training alters the human skeletal muscle epigenome and subsequent gene expression, it is still unknown whether it does so differently in men and women, potentially leading to sex-specific physiological adaptations. Elucidating sex differences in genetics, epigenetics, gene regulation and expression in response to exercise will have great health implications, as it may enable gene targets in future clinical interventions and may better individualised interventions. This review will cover this topic and highlight the recent findings of sex-specific genetic, epigenetic, and gene expression studies, address the gaps in the field, and offer recommendations for future research.

Cardiac Vagus and Exercise

Lower resting heart rate and high autonomic vagal activity are strongly associated with superior exercise capacity, maintenance of which is essential for general well-being and healthy aging. Recent evidence obtained in experimental studies using the latest advances in molecular neuroscience, combined with human exercise physiology, physiological modeling, and genomic data suggest that the strength of cardiac vagal activity causally determines our ability to exercise.

Cardiac vagal dysfunction and myocardial injury after non-cardiac surgery: a planned secondary analysis of the measurement of Exercise Tolerance before surgery study

Background

The aetiology of perioperative myocardial injury is poorly understood and not clearly linked to pre-existing cardiovascular disease. We hypothesised that loss of cardioprotective vagal tone [defined by impaired heart rate recovery ≤12 beats min⁻¹ (HRR ≤12) 1 min after cessation of preoperative cardiopulmonary exercise testing] was associated with perioperative myocardial injury.

Methods

We conducted a pre-defined, secondary analysis of a multi-centre prospective cohort study of preoperative cardiopulmonary exercise testing. Participants were aged ≥40 yr undergoing non-cardiac surgery. The exposure was impaired HRR (HRR≤12). The primary outcome was postoperative myocardial injury, defined by serum troponin concentration within 72 h after surgery. The analysis accounted for established markers of cardiac risk [Revised Cardiac Risk Index (RCRI), N-terminal pro-brain natriuretic peptide (NT pro-BNP)].

Results

A total of 1326 participants were included [mean age (standard deviation), 64 (10) yr], of whom 816 (61.5%) were male. HRR≤12 occurred in 548 patients (41.3%). Myocardial injury was more frequent amongst patients with HRR≤12 [85/548 (15.5%) vs HRR>12: 83/778 (10.7%); odds ratio (OR), 1.50 (1.08–2.08); P=0.016, adjusted for RCRI). HRR declined progressively in patients with increasing numbers of RCRI factors. Patients with ≥3 RCRI factors were more likely to have HRR≤12 [26/36 (72.2%) vs 0 factors: 167/419 (39.9%); OR, 3.92 (1.84–8.34); P<0.001]. NT pro-BNP greater than a standard prognostic threshold (>300 pg ml⁻¹) was more frequent in patients with HRR≤12 [96/529 (18.1%) vs HRR>12 59/745 (7.9%); OR, 2.58 (1.82–3.64); P<0.001].

Conclusions

Impaired HRR is associated with an increased risk of perioperative cardiac injury. These data suggest a mechanistic role for cardiac vagal dysfunction in promoting perioperative myocardial injury.

Examining the current standards for genetic discovery and replication in the era of mega-biobanks

With the recent deluge of mega-biobank data, it is time to revisit what constitutes “replication” for genome-wide association studies. Many replication samples are unavailable or underpowered, therefore alternatives beyond strict statistical replication are needed until the required resources become available.

Genome-wide association studies (GWAS) have become a mainstay in genetics research to understand genotype-phenotype relationships. Following the second release of UK Biobank data and the flood of publications using these data, here the author revisits the standards for discovery, replication and follow-up in GWAS today.

Thirty loci identified for heart rate response to exercise and recovery implicate autonomic nervous system

Impaired capacity to increase heart rate (HR) during exercise (ΔHR^ex), and a reduced rate of recovery post-exercise (ΔHR^rec) are associated with higher cardiovascular mortality rates. Currently, the genetic basis of both phenotypes remains to be elucidated. We conduct genome-wide association studies (GWASs) for ΔHR^ex and ΔHR^rec in ~40,000 individuals, followed by replication in ~27,000 independent samples, all from UK Biobank. Six and seven single-nucleotide polymorphisms for ΔHR^ex and ΔHR^rec, respectively, formally replicate. In a full data set GWAS, eight further loci for ΔHR^ex and nine for ΔHR^rec are genome-wide significant (P ≤ 5 × 10⁻⁸). In total, 30 loci are discovered, 8 being common across traits. Processes of neural development and modulation of adrenergic activity by the autonomic nervous system are enriched in these results. Our findings reinforce current understanding of HR response to exercise and recovery and could guide future studies evaluating its contribution to cardiovascular risk prediction.

Genome-wide association studies have identified multiple loci for resting heart rate (HR) but the genetic factors associated with HR increase during and HR recovery after exercise are less well studied. Here, the authors examine both traits in a two-stage GWAS design in up to 67,257 individuals from UK Biobank.