Pulmonary veins (PVs) are the major origin of atrial fibrillation. Recently, we recorded hyperpolarization-activated Cl⁻ current (I_{Cl, h}) in rat PV cardiomyocytes. Unlike the well-known chloride channel protein 2 (CLCN2) current, the activation curve of I_{Cl, h} was hyperpolarized as the Cl⁻ ion concentration ([Cl⁻]_i) increased. This current could account for spontaneous activity in PV cardiomyocytes linked to atrial fibrillation. In this study, we aimed to identify the channel underlying I_{Cl, h}. Using RT-PCR amplification specific for Clcn2 or its homologs, a chloride channel was cloned from rat PV and detected in rat PV cardiomyocytes using immunocytochemistry. The gene sequence and electrophysiological functions of the protein were identical to those previously reported for Clcn2, with protein activity observed as a hyperpolarization-activated current by the patch-clamp method. However, the [Cl⁻]_i dependence of activation was entirely different from the observed I_{Cl, h} of PV cardiomyocytes; the activation curve of the Clcn2-transfected cells shifted toward positive potential with increased [Cl⁻]_i, whereas the I_{Cl, h} of PV and left ventricular cardiomyocytes showed a leftward shift. Therefore, we used MS to explore the possibility of additional proteins interacting with CLCN2 and identified an individual 71-kDa protein, HSPA8, that was strongly expressed in rat PV cardiomyocytes. With co-expression of HSPA8 in HEK293 and PC12 cells, the CLCN2 current showed voltage-dependent activation and shifted to negative potential with increasing [Cl⁻]_i. Molecular docking simulations further support an interaction between CLCN2 and HSPA8. These findings suggest that CLCN2 in rat heart contains HSPA8 as a unique accessory protein.

The heart is a significant organ in mammalian life, and the heartbeat mechanism has been an essential focus of science. However, few studies have focused on species differences. Accordingly, challenges remain in studying genes that have universal functions across species and genes that determine species differences. Here, we analyzed transcriptome data in mouse, rat, and human atria, ventricles, and sinoatrial nodes (SA) obtained from different platforms and compared them by calculating specificity measure (SPM) values in consideration of species differences. Among the three heart regions, the species differences in SA were the greatest, and we searched for genes that determined the essential characteristics of SA, which was SHOX2 in our criteria. The SPM value of SHOX2 was prominently high across species. Similarly, by calculating SPM values, we identified 3 atrial-specific, 11 ventricular-specific, and 17 SA-specific markers. Ontology analysis identified 70 cardiac region- and species-specific ontologies. These results suggest that reanalyzing existing data by calculating SPM values may identify novel tissue-specific genes and species-dependent gene expression. This study identified the importance of SHOX2 as an SA-specific transcription factor, a novel cardiac regional marker, and species-dependent ontologies.

Atrial fibrillation (AF) is the most frequent persistent arrhythmia. Many genes have been reported as a genetic background for AF. However, most transcriptome analyses of AF are limited to the atrial samples and have not been evaluated by multiple cardiac regions. In this study, we analyzed the expression levels of protein-coding and long noncoding RNAs (lncRNAs) in six cardiac regions by RNA-seq. Samples were donated from six subjects with or without persistent AF for left atria, left atrial appendages, right atria, sinoatrial nodes, left ventricles, right ventricles, and pulmonary veins (PVs), and additional four right atrial appendages samples were collected from patients undergoing mitral valve replacement. In total, 23 AF samples were compared to 23 non-AF samples. Surprisingly, the most influenced heart region in gene expression by AF was the PV, not the atria. The ion channel-related gene set was significantly enriched upon analysis of these significant genes. In addition, some significant genes are cancer-related lncRNAs in PV in AF. A co-expression network analysis could detect the functional gene clusters. In particular, the cancer-related lncRNA, such as SAMMSON and FOXCUT, belong to the gene network with the cancer-related transcription factor FOXC1. Thus, they may also play an aggravating role in the pathogenesis of AF, similar to carcinogenesis. In the least, this study suggests that (1) RNA alteration is most intense in PVs and (2) post-transcriptional gene regulation by lncRNA may contribute to the progression of AF. Through the screening analysis across the six cardiac regions, the possibility that the PV region can play a role other than paroxysmal triggering in the pathogenesis of AF was demonstrated for the first time. Future research with an increase in the number of PV samples will lead to a novel understanding of the pathophysiology of AF.

Extracellular signal-regulated kinase 5 (ERK5) regulates diverse physiological responses such as proliferation, differentiation, and gene expression. Previously, we demonstrated that ERK5 is essential for neurite outgrowth and catecholamine biosynthesis in PC12 cells and sympathetic neurons. However, it remains unclear how ERK5 regulates the activity of ion channels, which are important for membrane excitability. Thus, we examined the effect of ERK5 on the ion channel activity in the PC12 cells that overexpress both ERK5 and the constitutively active MEK5 mutant. The gene and protein expression levels of voltage-dependent Ca²⁺ and K⁺ channels were determined by RT-qPCR or Western blotting. The A-type K⁺ current was recorded using the whole-cell patch clamp method. In these ERK5-activated cells, the gene expression levels of voltage-dependent L- and P/Q-type Ca²⁺ channels did not alter, but the N-type Ca²⁺ channel was slightly reduced. In contrast, those of K_v4.2 and K_v4.3, which are components of the A-type current, were significantly enhanced. Unexpectedly, the protein levels of K_v4.2 were not elevated by ERK5 activation, but the phosphorylation levels were increased by ERK5 activation. By electrophysiological analysis, the inactivation time constant of the A-type current was prolonged by ERK5 activation, without changes in the peak current. Taken together, ERK5 inhibits an inactivation of the A-type current by phosphorylation of K_v4.2, which may contribute to the neuronal differentiation process.

Ectopic excitability in pulmonary veins (PVs) is the major cause of atrial fibrillation. We previously reported that the inositol trisphosphate receptor in rat PV cardiomyocytes cooperates with the Na⁺-Ca²⁺ exchanger to provoke ectopic automaticity in response to norepinephrine. Here, we focused on adenylyl cyclase (AC) as another effector of norepinephrine stimulation. RT-PCR, immunohistochemistry, and Western blotting revealed that the abundant expression of Ca²⁺-stimulable AC3 was restricted to the supraventricular area, including the PVs. All the other AC isotypes hardly displayed any region-specific expressions. Immunostaining of isolated cardiomyocytes showed an enriched expression of AC3 along the t-tubules in PV myocytes. The cAMP-dependent response of L-type Ca²⁺ currents in the PV and LA cells is strengthened by the 0.1 mM intracellular Ca²⁺ condition, unlike in the ventricular cells. The norepinephrine-induced automaticity of PV cardiomyocytes was reversibly suppressed by 100 µM SQ22536, an adenine-like AC inhibitor. These findings suggest that the specific expression of AC3 along t-tubules may contribute to arrhythmogenic automaticity in rat PV cardiomyocytes.

Cardiomyocytes and myocardial sleeves dissociated from pulmonary veins (PVs) potentially generate ectopic automaticity in response to noradrenaline (NA), and thereby trigger atrial fibrillation. We developed a mathematical model of rat PV cardiomyocytes (PVC) based on experimental data that incorporates the microscopic framework of the local control theory of Ca²⁺ release from the sarcoplasmic reticulum (SR), which can generate rhythmic Ca²⁺ release (limit cycle revealed by the bifurcation analysis) when total Ca²⁺ within the cell increased. Ca²⁺ overload in SR increased resting Ca²⁺ efflux through the type II inositol 1,4,5-trisphosphate (IP₃) receptors (InsP₃R) as well as ryanodine receptors (RyRs), which finally triggered massive Ca²⁺ release through activation of RyRs via local Ca²⁺ accumulation in the vicinity of RyRs. The new PVC model exhibited a resting potential of −68 mV. Under NA effects, repetitive Ca²⁺ release from SR triggered spontaneous action potentials (APs) by evoking transient depolarizations (TDs) through Na⁺/Ca²⁺ exchanger (AP_TDs). Marked and variable latencies initiating AP_TDs could be explained by the time courses of the α1- and β1-adrenergic influence on the regulation of intracellular Ca²⁺ content and random occurrences of spontaneous TD activating the first AP_TD. Positive and negative feedback relations were clarified under AP_TD generation.