GWAS of the electrocardiographic QT interval in Hispanics/Latinos generalizes previously identified loci and identifies population-specific signals

Artificial intelligence-enhanced electrocardiography derived body mass index as a predictor of future cardiometabolic disease

The electrocardiogram (ECG) can capture obesity-related cardiac changes. Artificial intelligence-enhanced ECG (AI-ECG) can identify subclinical disease. We trained an AI-ECG model to predict body mass index (BMI) from the ECG alone. Developed from 512,950 12-lead ECGs from the Beth Israel Deaconess Medical Center (BIDMC), a secondary care cohort, and validated on UK Biobank (UKB) (n = 42,386), the model achieved a Pearson correlation coefficient (r) of 0.65 and 0.62, and an R2 of 0.43 and 0.39 in the BIDMC cohort and UK Biobank, respectively for AI-ECG BMI vs. measured BMI. We found delta-BMI, the difference between measured BMI and AI-ECG-predicted BMI (AI-ECG-BMI), to be a biomarker of cardiometabolic health. The top tertile of delta-BMI showed increased risk of future cardiometabolic disease (BIDMC: HR 1.15, p < 0.001; UKB: HR 1.58, p < 0.001) and diabetes mellitus (BIDMC: HR 1.25, p < 0.001; UKB: HR 2.28, p < 0.001) after adjusting for covariates including measured BMI. Significant enhancements in model fit, reclassification and improvements in discriminatory power were observed with the inclusion of delta-BMI in both cohorts. Phenotypic profiling highlighted associations between delta-BMI and cardiometabolic diseases, anthropometric measures of truncal obesity, and pericardial fat mass. Metabolic and proteomic profiling associates delta-BMI positively with valine, lipids in small HDL, syntaxin-3, and carnosine dipeptidase 1, and inversely with glutamine, glycine, colipase, and adiponectin. A genome-wide association study revealed associations with regulators of cardiovascular/metabolic traits, including SCN10A, SCN5A, EXOG and RXRG. In summary, our AI-ECG-BMI model accurately predicts BMI and introduces delta-BMI as a non-invasive biomarker for cardiometabolic risk stratification.

Arrhythmias and ion channelopathies causing sudden cardiac death in Hispanic/Latino and Indigenous populations

The limited literature and increasing interest in studies on cardiac electrophysiology, explicitly focusing on cardiac ion channelopathies and sudden cardiac death in diverse populations, has prompted a comprehensive examination of existing research. Our review specifically targets Hispanic/Latino and Indigenous populations, which are often underrepresented in healthcare studies. This review encompasses investigations into genetic variants, epidemiology, etiologies, and clinical risk factors associated with arrhythmias in these demographic groups. The review explores the Hispanic paradox, a phenomenon linking healthcare outcomes to socioeconomic factors within Hispanic communities in the United States. Furthermore, it discusses studies exemplifying this observation in the context of arrhythmias and ion channelopathies in Hispanic populations. Current research also sheds light on disparities in overall healthcare quality in Indigenous populations. The available yet limited literature underscores the pressing need for more extensive and comprehensive research on cardiac ion channelopathies in Hispanic/Latino and Indigenous populations. Specifically, additional studies are essential to fully characterize pathogenic genetic variants, identify population-specific risk factors, and address health disparities to enhance the detection, prevention, and management of arrhythmias and sudden cardiac death in these demographic groups.

POPDC1 Variants Cause Atrioventricular Node Dysfunction and Arrhythmogenic Changes in Cardiac Electrophysiology and Intracellular Calcium Handling in Zebrafish

Popeye domain-containing (POPDC) proteins selectively bind cAMP and mediate cellular responses to sympathetic nervous system (SNS) stimulation. The first discovered human genetic variant (POPDC1S201F) is associated with atrioventricular (AV) block, which is exacerbated by increased SNS activity. Zebrafish carrying the homologous mutation (popdc1S191F) display a similar phenotype to humans. To investigate the impact of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling, homozygous popdc1S191F and popdc1 knock-out (popdc1KO) zebrafish larvae and adult isolated popdc1S191F hearts were studied by functional fluorescent analysis. It was found that in popdc1S191F and popdc1KO larvae, heart rate (HR), AV delay, action potential (AP) and calcium transient (CaT) upstroke speed, and AP duration were less than in wild-type larvae, whereas CaT duration was greater. SNS stress by β-adrenergic receptor stimulation with isoproterenol increased HR, lengthened AV delay, slowed AP and CaT upstroke speed, and shortened AP and CaT duration, yet did not result in arrhythmias. In adult popdc1S191F zebrafish hearts, there was a higher incidence of AV block, slower AP upstroke speed, and longer AP duration compared to wild-type hearts, with no differences in CaT. SNS stress increased AV delay and led to further AV block in popdc1S191F hearts while decreasing AP and CaT duration. Overall, we have revealed that arrhythmogenic effects of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling in zebrafish are varied, but already present in early development, and that AV node dysfunction may underlie SNS-induced arrhythmogenesis associated with popdc1 mutation in adults.

Genetic variants and effect modifiers of QT interval prolongation in patients with sickle cell disease

BACKGROUND

Sickle cell disease (SCD) is a common inherited blood disorder among African Americans (AA), with premature mortality which has been associated with prolongation of the heart rate-corrected QT interval (QTc), a known risk factor for sudden cardiac death. Although numerous genetic variants have been identified as contributors to QT interval prolongation in the general population, their impact on SCD patients remains unclear. This study used an unweighted polygenic risk score (PRS) to validate the previously identified associations between SNPs and QTc interval in SCD patients, and to explore possible interactions with other factors that prolong QTc interval in AA individuals with SCD.

METHODS

In SCD patients, candidate genetic variants associated with the QTc interval were genotyped. To identify any risk SNPs that may be correlated with QTc interval prolongation, linear regression was employed, and an unweighted PRS was subsequently constructed. The effect of PRS on the QTc interval was evaluated using linear regression, while stratification analysis was used to assess the influence of serum alanine transaminase (ALT), a biomarker for liver disease, on the PRS effect. We also evaluated the PRS with the two subcomponents of QTc, the QRS and JTc intervals.

RESULTS

Out of 26 candidate SNPs, five risk SNPs were identified for QTc duration under the recessive model. For every unit increase in PRS, the QTc interval prolonged by 4.0 ms (95% CI: [2.0, 6.1]; p-value: <0.001) in the additive model and 9.4 ms in the recessive model (95% CI: [4.6, 14.1]; p-value: <0.001). Serum ALT showed a modification effect on PRS-QTc prolongation under the recessive model. In the normal ALT group, each PRS unit increased QTc interval by 11.7 ms (95% CI: [6.3, 17.1]; p-value: 2.60E-5), whereas this effect was not observed in the elevated ALT group (0.9 ms; 95% CI: [-7.0, 8.8]; p-value: 0.823).

CONCLUSION

Several candidate genetic variants are associated with QTc interval prolongation in SCD patients, and serum ALT acts as a modifying factor. The association of a CPS1 gene variant in both QTc and JTc duration adds to NOS1AP as evidence of involvement of the urea cycle and nitric oxide metabolism in cardiac repolarization in SCD. Larger replication studies are needed to confirm these findings and elucidate the underlying mechanisms.

Cardiac Repolarization in Health and Disease

Abnormal cardiac repolarization is at the basis of life-threatening arrhythmias in various congenital and acquired cardiac diseases. Dysfunction of ion channels involved in repolarization at the cellular level are often the underlying cause of the repolarization abnormality. The expression pattern of the gene encoding the affected ion channel dictates its impact on the shape of the T-wave and duration of the QT interval, thereby setting the stage for both the occurrence of the trigger and the substrate for maintenance of the arrhythmia. Here we discuss how research into the genetic and electrophysiological basis of repolarization has provided us with insights into cardiac repolarization in health and disease and how this in turn may provide the basis for future improved patient-specific management.

A large genome-wide association study of QT interval length utilizing electronic health records

QT interval length is an important risk factor for adverse cardiovascular outcomes; however, the genetic architecture of QT interval remains incompletely understood. We conducted a genome-wide association study of 76,995 ancestrally diverse Kaiser Permanente Northern California members enrolled in the Genetic Epidemiology Research on Adult Health and Aging cohort using 448,517 longitudinal QT interval measurements, uncovering 9 novel variants, most replicating in 40,537 individuals in the UK Biobank and Population Architecture using Genomics and Epidemiology studies. A meta-analysis of all 3 cohorts (n = 117,532) uncovered an additional 19 novel variants. Conditional analysis identified 15 additional variants, 3 of which were novel. Little, if any, difference was seen when adjusting for putative QT interval lengthening medications genome-wide. Using multiple measurements in Genetic Epidemiology Research on Adult Health and Aging increased variance explained by 163%, and we show that the ≈6 measurements in Genetic Epidemiology Research on Adult Health and Aging was equivalent to a 2.4× increase in sample size of a design with a single measurement. The array heritability was estimated at ≈17%, approximately half of our estimate of 36% from family correlations. Heritability enrichment was estimated highest and most significant in cardiovascular tissue (enrichment 7.2, 95% CI = 5.7-8.7, P = 2.1e-10), and many of the novel variants included expression quantitative trait loci in heart and other relevant tissues. Comparing our results to other cardiac function traits, it appears that QT interval has a multifactorial genetic etiology.

Complex interactions between p.His558Arg and linked variants in the sodium voltage-gated channel alpha subunit 5 (Na V 1.5)

Common genetic polymorphisms may modify the phenotypic outcome when co-occurring with a disease-causing variant, and therefore understanding their modulating role in health and disease is of great importance. The polymorphic p.His558Arg variant of the sodium voltage-gated channel alpha subunit 5 (Na _V 1.5) encoded by the SCN5A gene is a case in point, as several studies have shown it can modify the clinical phenotype in a number of cardiac diseases. To evaluate the genetic backgrounds associated with this modulating effect, we reanalysed previous electrophysiological findings regarding the p.His558Arg variant and further assessed its patterns of genetic diversity in human populations. The Na _V 1.5 p.His558Arg variant was found to be in linkage disequilibrium with six other polymorphic variants that previously were also associated with cardiac traits in GWAS analyses. On account of this, incongruent reports that Arg558 allele can compensate, aggravate or have no effect on Na _V 1.5, likely might have arose due to a role of p.His558Arg depending on the additional linked variants. Altogether, these results indicate a major influence of the epistatic interactions between SCN5A variants, revealing also that phenotypic severity may depend on the polymorphic background associated to each individual genome.

Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.

The Role of POPDC Proteins in Cardiac Pacemaking and Conduction

The Popeye domain-containing (POPDC) gene family, consisting of Popdc1 (also known as Bves), Popdc2, and Popdc3, encodes transmembrane proteins abundantly expressed in striated muscle. POPDC proteins have recently been identified as cAMP effector proteins and have been proposed to be part of the protein network involved in cAMP signaling. However, their exact biochemical activity is presently poorly understood. Loss-of-function mutations in animal models causes abnormalities in skeletal muscle regeneration, conduction, and heart rate adaptation after stress. Likewise, patients carrying missense or nonsense mutations in POPDC genes have been associated with cardiac arrhythmias and limb-girdle muscular dystrophy. In this review, we introduce the POPDC protein family, and describe their structure function, and role in cAMP signaling. Furthermore, the pathological phenotypes observed in zebrafish and mouse models and the clinical and molecular pathologies in patients carrying POPDC mutations are described.

Fish-Ing for Enhancers in the Heart

Precise control of gene expression is crucial to ensure proper development and biological functioning of an organism. Enhancers are non-coding DNA elements which play an essential role in regulating gene expression. They contain specific sequence motifs serving as binding sites for transcription factors which interact with the basal transcription machinery at their target genes. Heart development is regulated by intricate gene regulatory network ensuring precise spatiotemporal gene expression program. Mutations affecting enhancers have been shown to result in devastating forms of congenital heart defect. Therefore, identifying enhancers implicated in heart biology and understanding their mechanism is key to improve diagnosis and therapeutic options. Despite their crucial role, enhancers are poorly studied, mainly due to a lack of reliable way to identify them and determine their function. Nevertheless, recent technological advances have allowed rapid progress in enhancer discovery. Model organisms such as the zebrafish have contributed significant insights into the genetics of heart development through enabling functional analyses of genes and their regulatory elements in vivo. Here, we summarize the current state of knowledge on heart enhancers gained through studies in model organisms, discuss various approaches to discover and study their function, and finally suggest methods that could further advance research in this field.

Semi-Lethal Primary Ciliary Dyskinesia in Rats Lacking the Nme7 Gene

NME7 (non-metastatic cells 7, nucleoside diphosphate kinase 7) is a member of a gene family with a profound effect on health/disease status. NME7 is an established member of the ciliome and contributes to the regulation of the microtubule-organizing center. We aimed to create a rat model to further investigate the phenotypic consequences of Nme7 gene deletion. The CRISPR/Cas9 nuclease system was used for the generation of Sprague Dawley Nme7 knock-out rats targeting the exon 4 of the Nme7 gene. We found the homozygous Nme7 gene deletion to be semi-lethal, as the majority of SD^Nme7−/− pups died prior to weaning. The most prominent phenotypes in surviving SD^Nme7−/− animals were hydrocephalus, situs inversus totalis, postnatal growth retardation, and sterility of both sexes. Thinning of the neocortex was histologically evident at 13.5 day of gestation, dilation of all ventricles was detected at birth, and an external sign of hydrocephalus, i.e., doming of the skull, was usually apparent at 2 weeks of age. Heterozygous SD^Nme7+/− rats developed normally; we did not detect any symptoms of primary ciliary dyskinesia. The transcriptomic profile of liver and lungs corroborated the histological findings, revealing defects in cell function and viability. In summary, the knock-out of the rat Nme7 gene resulted in a range of conditions consistent with the presentation of primary ciliary dyskinesia, supporting the previously implicated role of the centrosomally located Nme7 gene in ciliogenesis and control of ciliary transport.

A genome-wide association and polygenic risk score study on abnormal electrocardiogram in a Chinese population

Electrocardiography is a common and widely-performed medical examination based on the measurement and evaluation of electrocardiogram (ECG) to assess the up-to-date cardiac rhythms and thus suggest the health conditions of cardiovascular system and on a larger level the individual’s wellness. Abnormal ECG assessment from the detection of abnormal heart rhythms may have clinical implications including blood clots in formation, ongoing heart attack, coronary artery blockage, etc. Past genetic-phenotypic research focused primarily on the physical parameters of ECG but not the medical evaluation. To unbiasedly uncover the underlying links of genetic variants with normal vs. abnormal ECG assessment, a genome-wide association study (GWAS) is carried out in a 1006-participant cohort of Chinese population effectively genotyped for 243487 single nucleotide polymorphisms (SNPs). Both age and sex are influential factors, and six novel SNPs are identified for potential association with abnormal ECG. With the selected SNPs, a polygenic risk score (PRS) differentiates the case–control subgroups, and correlates well with increased risk of abnormal ECG. The findings are reproduced in an independent validation cohort. The derived PRS may function as a potential biomarker for prospectively screening the high-risk subgroup of heart issues in the Chinese population.

Downregulation of M Current Is Coupled to Membrane Excitability in Sympathetic Neurons Before the Onset of Hypertension

Supplemental Digital Content is available in the text.

Neurohumoral activation is an early hallmark of cardiovascular disease and contributes to the etiology of the pathophysiology. Stellectomy has reemerged as a positive therapeutic intervention to modify the progression of dysautonomia, although the biophysical properties underpinning abnormal activity of this ganglia are not fully understood in the initial stages of the disease. We investigated whether stellate ganglia neurons from prehypertensive SHRs (spontaneously hypertensive rats) are hyperactive and describe their electrophysiological phenotype guided by single-cell RNA sequencing, molecular biology, and perforated patch clamp to uncover the mechanism of abnormal excitability. We demonstrate the contribution of a plethora of ion channels, in particular inhibition of M current to stellate ganglia neuronal firing, and confirm the conservation of expression of key ion channel transcripts in human stellate ganglia. We show that hyperexcitability was curbed by M-current activators, nonselective sodium current blockers, or inhibition of Na_v1.1-1.3, Na_v1.6, or I_NaP. We conclude that reduced activity of M current contributes significantly to abnormal firing of stellate neurons, which, in part, contributes to the hyperexcitability from rats that have a predisposition to hypertension. Targeting these channels could provide a therapeutic opportunity to minimize the consequences of excessive sympathetic activation.

Identification and Characterization of a Transcribed Distal Enhancer Involved in Cardiac Kcnh2 Regulation

The human ether-a-go-go-related gene KCNH2 encodes the voltage-gated potassium channel underlying I_Kr, a current critical for the repolarization phase of the cardiac action potential. Mutations in KCNH2 that cause a reduction of the repolarizing current can result in cardiac arrhythmias associated with long-QT syndrome. Here, we investigate the regulation of KCNH2 and identify multiple active enhancers. A transcribed enhancer ∼85 kbp downstream of Kcnh2 physically contacts the promoters of two Kcnh2 isoforms in a cardiac-specific manner in vivo. Knockdown of its ncRNA transcript results in reduced expression of Kcnh2b and two neighboring mRNAs, Nos3 and Abcb8, in vitro. Genomic deletion of the enhancer, including the ncRNA transcription start site, from the mouse genome causes a modest downregulation of both Kcnh2a and Kcnh2b in the ventricles. These findings establish that the regulation of Kcnh2a and Kcnh2b is governed by a complex regulatory landscape that involves multiple partially redundantly acting enhancers.

Population-based toxicity screening in human induced pluripotent stem cell-derived cardiomyocytes

The potential for cardiotoxicity is carefully evaluated for pharmaceuticals, as it is a major safety liability. However, environmental chemicals are seldom tested for their cardiotoxic potential. Moreover, there is a large variability in both baseline and drug-induced cardiovascular risk in humans, but data are lacking on the degree to which susceptibility to chemically-induced cardiotoxicity may also vary. Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes have become an important in vitro model for drug screening. Thus, we hypothesized that a population-based model of iPSC-derived cardiomyocytes from a diverse set of individuals can be used to assess potential hazard and inter-individual variability in chemical effects on these cells. We conducted concentration-response screening of 134 chemicals (pharmaceuticals, industrial and environmental chemicals and food constituents) in iPSC-derived cardiomyocytes from 43 individuals, comprising both sexes and diverse ancestry. We measured kinetic calcium flux and conducted high-content imaging following chemical exposure, and utilized a panel of functional and cytotoxicity parameters in concentration-response for each chemical and donor. We show reproducible inter-individual variability in both baseline and chemical-induced effects on iPSC-derived cardiomyocytes. Further, chemical-specific variability in potency and degree of population variability were quantified. This study shows the feasibility of using an organotypic population-based human in vitro model to quantitatively assess chemicals for which little cardiotoxicity information is available. Ultimately, these results advance in vitro toxicity testing methodologies by providing an innovative tool for population-based cardiotoxicity screening, contributing to the paradigm shift from traditional animal models of toxicity to in vitro toxicity testing methods.

GWAS of QRS duration identifies new loci specific to Hispanic/Latino populations

BACKGROUND

The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

METHODS

We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

RESULTS

We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

CONCLUSIONS

Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

A Multi-Trait Approach Identified Genetic Variants Including a Rare Mutation in RGS3 with Impact on Abnormalities of Cardiac Structure/Function

Heart failure is a major cause for premature death. Given the heterogeneity of the heart failure syndrome, identifying genetic determinants of cardiac function and structure may provide greater insights into heart failure. Despite progress in understanding the genetic basis of heart failure through genome wide association studies, the heritability of heart failure is not well understood. Gaining further insights into mechanisms that contribute to heart failure requires systematic approaches that go beyond single trait analysis. We integrated a Bayesian multi-trait approach and a Bayesian networks for the analysis of 10 correlated traits of cardiac structure and function measured across 3387 individuals with whole exome sequence data. While using single-trait based approaches did not find any significant genetic variant, applying the integrative Bayesian multi-trait approach, we identified 3 novel variants located in genes, RGS3, CHD3, and MRPL38 with significant impact on the cardiac traits such as left ventricular volume index, parasternal long axis interventricular septum thickness, and mean left ventricular wall thickness. Among these, the rare variant NC_000009.11:g.116346115C > A (rs144636307) in RGS3 showed pleiotropic effect on left ventricular mass index, left ventricular volume index and maximal left atrial anterior-posterior diameter while RGS3 can inhibit TGF-beta signaling associated with left ventricle dilation and systolic dysfunction.

Expression profiling of Nme7 interactome in experimental models of metabolic syndrome

Nucleoside diphosphate kinase 7, non-metastatic cells 7 (NME7) is an acknowledged member of ciliome and is involved in the biogenesis or function of cilia. As obesity and diabetes are common in several ciliopathies, we aimed to analyze changes of gene expression within Nme7 interactome in genetically designed rat models of metabolic syndrome. We assessed the liver transcriptome by Affymetrix microarrays in adult males of 14 PXO recombinant inbred rat strains and their two progenitor strains, SHR-Lx and BXH2. In the strains with the lowest expression of Nme7, we have identified significant enrichment of transcripts belonging to Nme7 interactome. In the subsequent network analysis, we have identified three major upstream regulators - Hnf4a, Ppara and Nr1h4 and liver steatosis (p=0.0001) and liver necrosis/cell death (apoptosis of liver cells, p=0.0003) among the most enriched Tox categories. The mechanistic network reaching the top score showed substantial overlap with Assembly of non-motile cilium and Glucose metabolism disorder gene lists. In summary, we show in a genetic model of metabolic syndrome that rat strains with the lowest expression of Nme7 present gene expression shifts of Nme7 interactome that are perturbing networks relevant for carbohydrate and lipid metabolism as well as ciliogenesis.