The potential causal association of systemic lupus erythematosus with congestive heart failure in the East Asian population: A two-sample mendelian randomization study.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that increases the risk of cardiovascular disease, but the causal relationship has remained unknown in the East Asian population. We aim to determine the causal relationship between SLE and congestive heart failure (CHF) in the East Asian population. METHODS: We conducted a two-sample Mendelian randomization MR study to examine the potential causal association of SLE with CHF, using the East Asian genome-wide association study datasets for SLE (4222 cases, 8431 controls) and CHF (9413 cases, 203,040 controls). RESULTS: The MR analysis showed that genetic susceptibility to SLE was associated with the increased risk of CHF (odds ratio [OR] 1.032, 95% confidence interval [CI] 1.004-1.061, p = .023). After excluding the human leukocyte antigen (HLA) regions, SLE was also associated with a higher risk of CHF (OR 1.034, 95% CI 1.007-1.062, p = .015). The multiple MR sensitivity analyses confirmed that this potential causal association was unlikely to be biased by horizontal pleiotropy. CONCLUSIONS: The findings of this MR study suggest that SLE potentially increases the risk of CHF in the East Asian population. Genetic predisposition to SLE may play a significant role in developing CHF in the East Asian population.

Effects of hydroxychloroquine on atrial electrophysiology in in silico wild-type and PITX2+/- atrial cardiomyocytes.

BACKGROUND: Hydroxychloroquine (HCQ) is commonly used in the treatment of autoimmune diseases and increases the risk of QT interval prolongation. However, it is unclear how HCQ affects atrial electrophysiology and the risk of atrial fibrillation (AF). METHODS: We quantitatively examined the potential atrial arrhythmogenic effects of HCQ on AF using a computational model of human atrial cardiomyocytes. We measured atrial electrophysiological markers after systematically varying HCQ concentrations. RESULTS: The HCQ concentrations were positively correlated with the action potential duration (APD), resting membrane potential, refractory period, APD alternans threshold, and calcium transient alternans threshold (p < 0.05). By contrast, HCQ concentrations were inversely correlated with the maximum upstroke velocity and calcium transient amplitude (p < 0.05). When the therapeutic concentration (Cmax) of HCQ was applied, HCQ increased APD90 by 1.4% in normal sinus rhythm, 1.8% in wild-type AF, and 2.6% in paired-like homeodomain transcription factor 2 (PITX2)+/- AF, but did not affect the alternans thresholds. The overall in silico results suggest no significant atrial arrhythmogenic effects of HCQ at Cmax, instead implying a potential antiarrhythmic role of low-dose HCQ in AF. However, at an HCQ concentration of fourfold Cmax, a rapid pacing rate of 4 Hz induced prominent APD alternans, particularly in the PITX2+/- AF model. CONCLUSION: Our in silico analysis suggests a potential antiarrhythmic role of low-dose HCQ in AF. Concomitant PITX2 mutations and high-dose HCQ treatments may increase the risk of AF, and this potential genotype/dose-dependent arrhythmogenic effect of HCQ should be investigated further.

Impact of noise on the instability of spiral waves in stochastic 2D mathematical models of human atrial fibrillation.

Sustained spiral waves, also known as rotors, are pivotal mechanisms in persistent atrial fibrillation (AF). Stochasticity is inevitable in nonlinear biological systems such as the heart; however, it is unclear how noise affects the instability of spiral waves in human AF. This study presents a stochastic two-dimensional mathematical model of human AF and explores how Gaussian white noise affects the instability of spiral waves. In homogeneous tissue models, Gaussian white noise may lead to spiral-wave meandering and wavefront break-up. As the noise intensity increases, the spatial dispersion of phase singularity (PS) points increases. This finding indicates the potential AF-protective effects of cardiac system stochasticity by destabilizing the rotors. By contrast, Gaussian white noise is unlikely to affect the spiral-wave instability in the presence of localized scar or fibrosis regions. The PS points are located at the boundary or inside the scar/fibrosis regions. Localized scar or fibrosis may play a pivotal role in stabilizing spiral waves regardless of the presence of noise. This study suggests that fibrosis and scars are essential for stabilizing the rotors in stochastic mathematical models of AF. Further patient-derived realistic modeling studies are required to confirm the role of scar/fibrosis in AF pathophysiology.

Interpretable machine learning of action potential duration restitution kinetics in single-cell models of atrial cardiomyocytes.

Action potential duration (APD) restitution curve and its maximal slope (Smax) reflect single cell-level dynamic instability for inducing chaotic heart rhythms. However, conventional parameter sensitivity analysis often fails to describe nonlinear relationships between ion channel parameters and electrophysiological phenotypes, such as Smax. We explored the parameter-phenotype mapping in a population of 5000 single-cell atrial cell models through interpretable machine learning (ML) approaches. Parameter sensitivity analyses could explain the linear relationships between parameters and electrophysiological phenotypes, including APD90, resting membrane potential, Vmax, refractory period, and APD/calcium alternans threshold, but not for Smax. However, neural network models had better prediction performance for Smax. To interpret the ML model, we evaluated the parameter importance at the global and local levels by computing the permutation feature importance and the local interpretable model-agnostic explanations (LIME) values, respectively. Increases in ICaL, INCX, and IKr, and decreases in IK1, Ib,Cl, IKur, ISERCA, and Ito are correlated with higher Smax values. The LIME algorithm determined that INaK plays a significant role in determining Smax as well as Ito and IKur. The atrial cardiomyocyte population was hierarchically clustered into three distinct groups based on the LIME values and the single-cell simulation confirmed that perturbations in INaK resulted in different behaviors of APD restitution curves in three clusters. Our combined top-down interpretable ML and bottom-up mechanistic simulation approaches uncovered the role of INaK in heterogeneous behaviors of Smax in the atrial cardiomyocyte population.

Network determinants of cardiovascular calcification and repositioned drug treatments.

Ectopic cardiovascular calcification is a highly prevalent pathology for which there are no effective novel or repurposed pharmacotherapeutics to prevent disease progression. We created a human calcification endophenotype module (ie, the "calcificasome") by mapping vascular calcification genes (proteins) to the human vascular smooth muscle-specific protein-protein interactome (218 nodes and 632 edges, P < 10-5 ). Network proximity analysis was used to demonstrate that the calcificasome overlapped significantly with endophenotype modules governing inflammation, thrombosis, and fibrosis in the human interactome (P < 0.001). A network-based drug repurposing analysis further revealed that everolimus, temsirolimus, and pomalidomide are predicted to target the calcificasome. The efficacy of these agents in limiting calcification was confirmed experimentally by treating human coronary artery smooth muscle cells in an in vitro calcification assay. Each of the drugs affected expression or activity of their predicted target in the network, and decreased calcification significantly (P < 0.009). An integrated network analytical approach identified novel mediators of ectopic cardiovascular calcification and biologically plausible candidate drugs that could be repurposed to target calcification. This methodological framework for drug repurposing has broad applicability to other diseases.

Persistent homology analysis of type 2 diabetes genome-wide association studies in protein-protein interaction networks.

Genome-wide association studies (GWAS) involving increasing sample sizes have identified hundreds of genetic variants associated with complex diseases, such as type 2 diabetes (T2D); however, it is unclear how GWAS hits form unique topological structures in protein-protein interaction (PPI) networks. Using persistent homology, this study explores the evolution and persistence of the topological features of T2D GWAS hits in the PPI network with increasing p-value thresholds. We define an n-dimensional persistent disease module as a higher-order generalization of the largest connected component (LCC). The 0-dimensional persistent T2D disease module is the LCC of the T2D GWAS hits, which is significantly detected in the PPI network (196 nodes and 235 edges, P<0.05). In the 1-dimensional homology group analysis, all 18 1-dimensional holes (loops) of the T2D GWAS hits persist over all p-value thresholds. The 1-dimensional persistent T2D disease module comprising these 18 persistent 1-dimensional holes is significantly larger than that expected by chance (59 nodes and 83 edges, P<0.001), indicating a significant topological structure in the PPI network. Our computational topology framework potentially possesses broad applicability to other complex phenotypes in identifying topological features that play an important role in disease pathobiology.