A phenotype-based forward genetic screen identifies Dnajb6 as a sick sinus syndrome gene
Abstract
Previously we showed the generation of a protein trap library made with the gene-break transposon (GBT) in zebrafish (Danio rerio) that could be used to facilitate novel functional genome annotation towards understanding molecular underpinnings of human diseases (Ichino et al, 2020). Here, we report a significant application of this library for discovering essential genes for heart rhythm disorders such as sick sinus syndrome (SSS). SSS is a group of heart rhythm disorders caused by malfunction of the sinus node, the heart’s primary pacemaker. Partially owing to its aging-associated phenotypic manifestation and low expressivity, molecular mechanisms of SSS remain difficult to decipher. From 609 GBT lines screened, we generated a collection of 35 zebrafish insertional cardiac (ZIC) mutants in which each mutant traps a gene with cardiac expression. We further employed electrocardiographic measurements to screen these 35 ZIC lines and identified three GBT mutants with SSS-like phenotypes. More detailed functional studies on one of the arrhythmogenic mutants, GBT411, in both zebrafish and mouse models unveiled Dnajb6 as a novel SSS causative gene with a unique expression pattern within the subpopulation of sinus node pacemaker cells that partially overlaps with the expression of hyperpolarization activated cyclic nucleotide gated channel 4 (HCN4), supporting heterogeneity of the cardiac pacemaker cells.
Editor's evaluation
This study presents a valuable discovery of a gene important for the function of the cardiac pacemaker. The evidence is convincing as mutation in this gene causes sick sinus syndrome (SSS) in both zebrafish and mice, and potentially in humans. This manuscript is of interest to scientists in the field of cardiology, particular cardiac electrophysiology and arrhythmia.
https://doi.org/10.7554/eLife.77327.sa0Introduction
Cardiac arrhythmia affects >2% of individuals in community-dwelling adults (Khurshid et al., 2018). Sick sinus syndrome (SSS), also known as sinus node dysfunction or sinoatrial node (SAN) disease, is a group of heart rhythm disorders affecting cardiac impulse formation and/or propagation from the SAN, the heart’s primary pacemaker. SSS manifests a spectrum of presentations such as sinus pause or arrest (SA), bradycardia, sinoatrial exit block, or tachy-brady syndrome accompanied by atrial fibrillation (AF) (Semelka et al., 2013; De Ponti et al., 2018). In addition, 20% to 60% SSS patients show abnormal response to autonomic stresses (Dakkak and Doukky, 2020). SSS occurs most commonly in elderly, with an estimated prevalence of 1 case per 600 adults over age 65 (Dobrzynski et al., 2007). Symptomatic SSS can lead to inadequate blood supply to the heart and body and contribute significantly to life-threatening problems such as heart failure and cardiac arrest. While SSS is the most common indication for pacemaker implantation worldwide (Mond and Proclemer, 2011), the mechanisms of SSS remain poorly understood, making it difficult to stratify SSS risk in vulnerable cohorts of patients and development of effective pharmacologic therapy for pacemaker abnormalities.
To develop mechanism-based diagnostic and therapeutic strategies for SSS, it is desirable to discover genes that are expressed in the SAN and may contribute to SSS. Unfortunately, very limited number of SSS genes and related animal models are currently available. While mutations in the cardiac sodium channel α-subunit encoding gene (SCN5A) (Nof et al., 2007; Tan et al., 2001) and hyperpolarization-activated cyclic nucleotide-aged channel encoding gene (HCN4) Schulze-Bahr et al., 2003; Verkerk and Wilders, 2015 have been found to cause SSS, only a few other genes affecting the structure and/or function of the SAN were identified to increase the risk of developing SSS (Anderson and Benson, 2010; Holm et al., 2011). Classic human genetic linkage analysis-based approach has played important roles in gene discovery, but it is largely limited by the availability of suitable pedigree, especially in this age-dependent disease (Zhu et al., 2018). More recently, the genome-wide association studies (GWASs) have been used to identify novel genetic susceptibility factors associated with SSS (Holm et al., 2011; Monfredi and Boyett, 2015). However, owing to its statistic and associative nature, it has been difficult to confidently establish genotype-phenotype relationships for the vast amount of variants (Lin and Musunuru, 2018; Tam et al., 2019). Alternative approaches for effective identification of essential genes for SSS are thus needed.
Phenotype-based forward genetic screen in model organisms is a powerful strategy for deciphering genetic basis of a biological process. Without any a prior assumption, new genes can be identified that shed light on key signaling pathways. However, this approach is difficult to carry out in adult vertebrates, because of significantly increased burden of colony management efforts (Kamp et al., 2010; Shen et al., 2005). To address this bottleneck, zebrafish, a vertebrate with higher throughput than rodents, has been explored to study cardiac diseases (Gut et al., 2017). Despite its small body size, a zebrafish heart has conserved myocardium, endocardium, and epicardium as found in human, and adult zebrafish shows strikingly similar cardiac physiology to humans (Bakkers, 2011). Its heart rate is around 100 beats per minute (bpm), which is much comparable to that in human than in rodents. Adult zebrafish models for human cardiac diseases such as cardiomyopathies have been successfully generated (Ding et al., 2020a). Besides N-ethyl-N-nitrosourea (ENU)-based mutagenesis screens that have been conducted to identify embryonic recessive mutants, insertional mutagens such as those based on viruses and/or transposons have been developed to further increase the throughput of the screen, opening doors to screening genes affecting adult phenotypes (Amsterdam et al., 1999; Wang et al., 2007). Our team recently reported a gene-breaking transposon (GBT)-based gene-trap system in zebrafish which enables to disrupt gene function reversibly at high efficiency (>99% at the RNA level) (Clark et al., 2011). Approximately 1,200 GBT lines have been generated, laying a foundation for adult phenotype-based forward genetic screens (Ichino et al., 2020).Because the expression pattern of the affected genes in each GBT line is reported by a fluorescence reporter, we enriched GBT lines with cardiac expression and generated a zebrafish insertional cardiac (ZIC) mutant collection (Ding et al., 2013). Through stressing the ZIC collection with doxorubicin, an anti-cancer drug, we demonstrated that novel genetic factors of doxorubicin-induced cardiomyopathy (DIC), such as Dnaj (Hsp40) homology, subfamily B, member 6b (dnajb6b), sorbin and SH3 domain-containing 2b (sorbs2b) and retinoid x receptor alpha a (rxraa), could be successfully identified (Ding et al., 2016; Ding et al., 2020b; Ma et al., 2020). Follow up studies on these hits confirmed their identity as important cardiomyopathy genes.
Encouraged by our success in identifying new genetic factors for DIC, we reasoned that genes for rhythm disorders could be similarly identified by directly screening adult ZIC lines using echocardiographic measurement. We had recently optimized a commercially available ECG system to define SA episodes in an adult zebrafish, and the baseline frequency of aging-associated SSS in wild-type (WT) adult zebrafish (Yan et al., 2020). Here, we reported a pilot screen of our ZIC collection using this ECG platform and the resultant discovery of three positive hits, followed by comprehensive expressional and functional analysis of dnajb6b gene that is linked to one of the hits. Together, our data prove the feasibility of a phenotype-based screening strategy in adult zebrafish for discovering new rhythm genes.
Results
Identification of 35 zebrafish insertional cardiac (ZIC) mutants
We recently reported the generation of more than 1200 zebrafish mutant strains using the gene-break transposon (GBT) vector (Ichino et al., 2020). The tagged gene in each GBT mutant is typically disrupted with 99% knockdown efficiency and its expression pattern is reported by a monomeric red fluorescent protein (mRFP) reporter Ichino et al., 2020. We screened 609 GBT lines based on their mRFP expression and identified 44 mutants with either the embryonic or adult heart expression Ding et al., 2016. Then, we outcrossed these 44 lines, aided by Southern blotting to identify offsprings with a lower copy number of insertions, Ding et al., 2013 and identified 35 mutants with a single copy of the GBT insertion after 2–4 generations of outcross (Table 1; Ding et al., 2016). Using a combination of inverse PCR and/or 5’- and 3’-RACE PCR cloning approaches, we mapped the genetic loci of GBT inserts in these 35 mutants (Table 1; Ding et al., 2013). Most of the affected genes have human orthologs with a corresponding Online Mendelian Inheritance in Man (OMIM) number. Because each GBT line contains a single GBT insertion that traps a gene with cardiac expression, these 35 GBT lines were termed as zebrafish insertional cardiac (ZIC) mutants.
An ECG screen of 35 ZIC lines identified three mutants with increased incidence of SA and/or AV block episodes
Because each ZIC mutant disrupts a gene with cardiac expression, we enquired whether an ECG screening can be conducted to identify genetic lesions that result in arrhythmia. Since aging is a strong risk factor for heart rhythm disorders, we initially carried our screen in 35 aged ZIC fish lines generated from incrosses to facilitate the manifestation of cardiac rhythm abnormalities (Table 1). Because these fish are offsprings of incrosses and have been preselected based on the mRFP tag, their genotypes consist of both heterozygous and homozygous for the affected genes. As reported recently, in WT fish aged around 2 years old, we noted baseline SA episodes in about 1 out of 20 fish (5%) fish Yan et al., 2020. By contrast, among the 35 ZIC mutants with mixed heterozygous and homozygous genotypes, we noted an increased incidence of SA in 3 lines, including 3 out of 13 GBT103 mutant fish at 1.5 years old, 4 out of 10 GBT410 mutant fish at 2 years old, and 3 out of 8 GBT411 mutant fish at 2 years old (Figure 1A). In addition to SA, we also noted incidence of atrioventricular block (AVB) in 4 different GBT103 mutant animals at 1.5 years of age. Because the increased incidence of SA and/or AVB are hallmarks of SSS, these three lines were thus identified as three candidate SSS-like mutants.
To confirm the linkage between genetic lesions and the SSS-like phenotypes, we focused on homozygous animals for further validation. Because the precise insertional positions for all the 35 ZIC lines have been mapped, all these three homozygous ZIC mutants were easily identified by genotyping PCR Table 1, Figure 1B; Clark et al., 2011; Ichino et al., 2020; Ding et al., 2013. In contrast to 5% WT fish whereby SA episodes can be detected, significantly increased SA incidence was noted in all three homozygous mutants at 16 months of age, with an incidence of 57.1% in the GBT103/cyth3a, 44.4% in the GBT410/vapal, and 40% in the GBT411/dnajb6b homozygous mutants, respectively (Table 2). There was one animal manifesting AVB in the GBT103/cyth3a and GBT411/dnajb6b homozygous mutants, respectively. In addition, we also noted a reduced heart rate, another SSS phenotypic trait in the GBT411/dnajb6b, but not the other two GBT homozygous mutants (Table 2).
To seek additional evidence supporting our screening strategy, we decided to focus on the GBT411/dnajb6b mutant that is also characterized with significantly reduced heart rate phenotype. Detailed analysis of ECG indices showed increased RR interval in the GBT411/dnajb6b homozygous mutants (GBT411-/-) (Supplementary file 2), which is consistent with reduced heart rate. No obvious abnormality on other ECG indices was detected. Because arrhythmic mutants often manifest an aberrant response to extrinsic regulation of the heart rate, we examined responses of GBT411-/- to autonomic stimuli by stressing them with three compounds, including isoproterenol, a β-adrenoreceptor agonist for sympathetic nervous system; atropine, an anticholinergic inhibitor; and carbachol, a cholinergic agonist for parasympathetic nervous system. After administrating these drugs to the GBT411-/- fish at 1 year old via intraperitoneal (IP) injection, we noted aberrant heart rate response to both atropine and carbachol, while its response to isoproterenol appeared to be similar to that in WT control animals (Figure 1—figure supplement 1).
Next, we stressed the GBT411-/- fish with verapamil, an L-type Ca2+ channel antagonists, to stress out cardiac pacemaking and unmask SSS phenotype. Indeed, SA incidence was significantly increased in the GBT411-/- fish at 10 months of age (Supplementary file 1). Similarly, the heart rate was significantly reduced in the GBT411-/- fish compared to WT controls. Together, these data provided additional evidence to support GBT411/dnajb6b as an arrhythmia mutant.
Dnajb6b and its mouse ortholog exhibit unique expression patterns in the cardiac conduction system
Dnajb6b was previously identified as a cardiomyopathy-associated gene, disruption of which led to abnormal cardiac remodeling in zebrafish, Ding et al., 2016 raising concerns on whether the arrhythmic phenotype in the GBT411/dnajb6b mutant is a primary defect in the cardiac conduction system or a consequence of cardiac remodeling in cardiomyocytes. To address this concern, we firstly defined the expression of the Dnajb6b protein in the zebrafish heart. Our previous characterization of the mRFP reporter in the GBT411/dnajb6b fish revealed expression of Dnajb6b protein in both the embryonic and the adult hearts Ding et al., 2013; Ding et al., 2016. To enquire its expression in the cardiac conduction system (CCS), we crossed the GBT411/dnajb6b line into the sqET33-mi59B transgenic line in which EGFP labels the zebrafish SAN and atrio-ventricular canal (AVC) cells (Poon et al., 2016). Co-localization analysis demonstrated that the mRFP positive, Dnajb6b-expressing cells partially overlap with the EGFP signal labeling both AVC and SAN cells at the base of atrium in the embryonic heart at 3 days post-fertilization (Figure 2A-C). In the GBT411/dnajb6b heterozygous (GBT411+/-) adult hearts crossed with the sqET33-mi59B line, EGFP signal labeling the AVC and SAN cells were consistently detected in all animals (Figure 2D-F; Poon et al., 2016). However, in the GBT411/dnajb6b homozygous (GBT411-/-) adult hearts crossed with the sqET33-mi59B line, EGFP signal in the AVC region appeared to be more diffused compared to that in GBT411+/-, while no EGFP-positive SAN cells was detected in 3 out of 9 fish hearts examined (Figure 2G and H). Together, these results underscored the expression of Dnajb6b in the CCS, and disruption of dnajb6b in the GBT411-/- mutant altered the CCS expression which might contribute to the observed SSS-like phenotypes. It should be noted that the Dnajb6b-mRFP-positive expression patterns overlap with but extend beyond the sqET33-mi59B EGFP-positive expression patterns in both embryonic and adult fish hearts (Figure 2).
To seek additional evidence supporting expression and function of dnajb6b in the CCS, we turned to the mouse model. The mouse DNAJB6 protein can be detected in all four cardiac chambers in a sectioned mouse heart tissue (Figure 3—figure supplement 1). Interestingly, we found a highly enriched expression of DNAJB6 in the SAN region, as defined by the expression of HCN4 channels which are responsible for the generation of hyperpolarization-activated pacemaker ‘funny’ current in pacemaker cells (Figure 3A). However, at higher magnification images, only a proportion of DNAJB6-positive cells showed colocalization with the HCN4-positive cells (arrows for colocalized cells vs. stars for non-colocalized cells in Figure 3B). In addition, co-localization of DNAJB6 with TBX3, a transcription factor that specifies the formation of the SAN cells, was noted (Figure 3C and D). More interestingly, a negative correlation between DNAJB6 and TBX3 signal intensity was appreciated: cells with strong DNAJB6 expression tend to overlap with cells that show weak TBX3 signal, while cells with weak DNAJB6 expression tend to overlap with the cells with strong TBX3 signal (Figure 3C-E). Furthermore, the overall TBX3 signal in the SAN tissue of DNAJB6 heterozygous knock out (KO) mouse (see below) was significantly increased compared to that in WT control (Figure 3F). Together, these results uncovered a unique expression of DNAJB6 in the SAN region which might contribute to SSS development; however, its unique expression patterns also underscored heterogeneity of pacemaker cells within the SAN (Boyett et al., 2000; Liang et al., 2021).
The Dnajb6+/- mice manifest features of SSS when there is no sign of cardiomyopathy
To test the conservation of the cardiac arrhythmic functions of Dnajb6b suggested from zebrafish, we obtained a global Dnajb6 KO mouse line. The mutant harbors a deletion of 36,843 bp nucleotides spanning from the first intron to the last intron of Dnajb6 gene located in the Chromosome 5, which was created by the insertion of the Velocigene ZEN-Ub1 cassette and subsequent LoxP excision using Cre (Figure 4A). Genotyping PCR using a combination of the Dnajb6 gene-specific and the Zen-Ubi cassette-specific primers was carried out to identify both Dnajb6 heterozygous (Dnajb6+/-) and homozygous (Dnajb6-/-) KO mice (Figure 4B). At the protein level, both the DNAJB6 short (S) and long (L) isoforms were reduced by ~50% in Dnajb6+/- mouse embryonic hearts at E12.5 stage (+/-), and near completely depleted in Dnajb6-/- mutant hearts. Consistent with a previous report, Hunter et al., 1999 Dnajb6-/- KO mice were embryonic lethal, died in the uterus at about E13.5 stage, likely due to the placental defects (data not shown). The Dnajb6 +/- mice were able to grow to adulthood without visually noticeable phenotypes until at least 1 year of age. Cardiac mechanical function remained normal, as indicated by indistinguishable cardiac echocardiography indices from those of WT siblings at the same age (Table 3). No abnormal myocardial structural morphology was detected in the left ventricle (LV) of Dnajb6+/- mice (Figure 4—figure supplement 1). However, increased frequency of SA and AVB episodes, as well as bradycardia phenotype, were noted in Dnajb6+/- mice at 6 months old (Figure 4D and E, and Table 4). Other ECG indices such as PR interval, QRS duration and QT interval remained comparable to WT control (Supplementary file 3). Similar to the GBT411/dnajb6b mutant in zebrafish, the Dnajb6+/- mice exhibited an impaired response to autonomic stimuli including isoproterenol and carbachol (Figure 4E). Together, these studies suggest that Dnajb6+/- mice manifest SSS phenotype without structural/functional remodeling of the heart.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/77327/elife-77327-fig4-data1-v2.zip
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/77327/elife-77327-fig4-data2-v2.zip
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/77327/elife-77327-fig4-data3-v2.zip
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/77327/elife-77327-fig4-data4-v2.zip
Ex vivo evidences of SAN dysfunction in the Dnajb6+/- mice
To further prove SAN dysfunction in Dnajb6+/- mice, we performed electrophysiological assessment of SAN pacemaker function by high-resolution fluorescent optical mapping of action potentials from isolated mouse atria at 1 year of age. We firstly analyzed the distribution of the leading pacemaker location site in Dnajb6+/- mice compared to WT control. In WT mice, leading pacemakers were mostly located within the anatomically and functionally defined SAN region Figure 5A and B; Gut et al., 2017; Glukhov et al., 2010; Liu et al., 2007; Verheijck et al., 2001. In contrast, significant increase in the number of leading pacemakers located outside of the SAN, including the subsidiary atrial pacemakers and inter-atrial septum pacemakers, was observed in Dnajb6+/- mice (p=0.039 vs. WT mice). In addition, in Dnajb6+/- mice, we also found a highly irregular heart rate, accompanied by the presence of multiple competing pacemakers and a beat-to-beat migration of the leading pacemaker between various sites which included SAN, right atrial ectopic (subsidiary) pacemakers, and inter-atrial septum (Figure 5C and D). Similar to the results from the in vivo studies, bradycardia phenotype was consistently detected in the isolated atrial preparations as well (Figure 5E). Optical mapping on isolated atrial preparations further revealed different responses of heart rate during isoproterenol, atropine, and carbachol stimulations in Dnajb6+/- mice. Significantly increased cycle length (CL) variations were also observed at baseline and upon carbachol stimulation (Figure 5F).
Furthermore, in Dnajb6+/- mice, we found significant prolongation of the SAN recovery time corrected to beating rate (cSANRT) measured both at baseline and under autonomic stresses, including stimulation by isoproterenol, carbachol, and atropine (Figure 5—figure supplement 1), confirming the presence of SAN dysfunction in Dnajb6+/- mice. Optical mapping also showed that, unlike WT, the first spontaneous post-pacing atrial beats during SANRT measurements in Dnajb6+/- mice were originated from ectopic locations outside of the SAN (Figure 5—figure supplement 1A, B), further supporting a suppressed SAN function. Histological evaluation of fibrosis tissue content in Dnajb6+/- mouse atria, however, did not reveal any significant difference compared to WT mice, for both atria and SAN.
Computational analysis of the cellular mechanisms underlying the SSS phenotype
To determine the potential cellular/ionic mechanisms underlying the observed SSS phenotype in the Dnajb6+/- mice , we utilized a population-based computational modeling approach. We used our previously published model of the mouse SAN myocyte to generate a population of 10,000 model variants by randomly varying selected model parameters (Figure 6A and B; Morotti et al., 2021). In each variant, we simulated both sympathetic and parasympathetic stimulations and recorded baseline heart rate and heart rate responses to autonomic stimuli. Simulations of both our original model and population predicted an increase in heart rate with isoproterenol and heart rate slowing with carbachol. Nevertheless, the cell-cell variability in heart rate response allowed identifying two subpopulations of model variants, whereby several models (n=438) displayed a slower firing rate at baseline, an increased response to isoproterenol, and a decreased response to carbachol administration (Figure 6C and D), thus recapitulating the Dnajb6+/- mice as measured in our ex vivo functional experiments (Figure 5E). The WT subpopulation comprised of the remaining n=6995 models. To reveal the ionic processes that are associated with the observed electrophysiological differences in Dnajb6+/- vs. WT, we then compared the parameter values ( the randomly applied scaling factors) in the two model subpopulations and found significant differences in several model parameters (Figure 6D–F). The analysis revealed a significant decrease in the maximal conductance of the fast (Nav1.5) Na+ current, the L-type Ca2+ current (ICa,L), the transient outward, sustained, and acetylcholine-activated K+ currents, the background Na+ and Ca2+ currents, as well as the ryanodine receptor maximal release flux of the Dnajb6+/- vs. WT model variants. We also found a significant increase in the Na+/Ca2+-exchanger maximal transport rate, and conductance of the T-type Ca2+ current and the slowly-activating delayed rectifier K+ current.
Transcriptome analysis of the Dnajb6+/- mutant hearts identifies altered genes encoding ion channels and proteins in the Wnt/beta-catenin pathway
To further seek molecular mechanisms underlying the SSS phenotypes observed in Dnajb6+/- mice, we performed whole transcriptome RNA-sequencing experiments using right atrial tissues isolated from Dnajb6+/- mice WT mice at 1 year of age. Transcriptomes of biological replicates for Dnajb6+/- mice did form a cluster that differs from the cluster for WT control samples, as indicated by principal component analysis (PCA) (Figure 7—figure supplement 1A). Based on a cut-off of adjusted p-value <0.05 and≥2 folds change, 107 differentially expressed (DE) genes were identified, among which 37 genes were upregulated and 70 genes were downregulated in the Dnajb6+/- mice compared with WT controls (Figure 7—figure supplement 1B, C). Through Ingenuity pathway analysis (IPA), several diverse signaling pathways were identified to be altered in the Dnajb6+/- mice (Figure 7—figure supplement 1D). Among these 107 differentially expressed genes, we noted calcium handling related protein-encoding genes like Slc24a2 and Cdh20, ion channel-encoding genes including Slc9a3r1, Kcnh7, Fxyd5, and Gjb5 (Figure 7A), as well as 4 Wnt pathway related genes (Figure 7B). We then performed quantitative RT-PCR analysis and experimentally confirmed dysregulation of these genes in the Dnajb6+/- mice (Figure 7C). The data on calcium handling and ion channel-encoding genes are in line with the SAN dysfunction phenotype observed in the Dnajb6+/- mice. Because Wnt signaling has been shown to direct pacemaker cell specification during SAN morphogenesis, Liang et al., 2020; Ren et al., 2019 the identification of 4 Wnt pathway related genes suggested that this SAN developmental pathway might play an important role in the observed SAN dysfunction phenotype in the adult Dnajb6+/- mice. Future studies are warranted to test this possibility.
Identification of DNAJB6 sequence variants associated with human SSS patients
To investigate the potential role of DNAJB6 in human SSS, we queried a sequence variant dataset derived from a genome-wide association study (GWAS) of 6,469 SSS cases and 1,000,187 controls Thorolfsdottir et al., 2021. Out of 313 variants with minor allele frequency ≥1% in DNAJB6, four variants showed nominal association (p<0.05, Supplementary file 4). Although none of the four variants survived Bonferroni correction, it’s interesting that two were located in untranslated regions. The most significant variant was observed for rs754941044 (p=0.0193), which was predicted as a splice acceptor variant by the Ensembl variant effect predictor McLaren et al., 2016. Thus, it is likely this variant has a significant impact on DNAJB6 gene function.
Discussion
GBT lines enable a phenotype-based screening approach for discovering new SSS genes
This work is based on recent establishment of a GBT protein trap-based insertional mutagenesis screening strategy and the generation of a collection of 1,200 zebrafish mutant strains Ichino et al., 2020. Here, we demonstrated the feasibility of screening these GBT lines for discovering new genetic factors for SSS, an aging-associated human disease. To overcome the challenge of colony management efforts that is associated with an adult screen, we leveraged the following unique advantages of the GBT vectors and zebrafish models. First, the knockdown efficiency for the tagged gene in each GBT homozygous mutant is consistently high, which is typically >99%, which ensued the success of an adult screen. Second, because of a fluorescence tag, heterozygous GBT fish can be easily identified under a fluorescent microscope without the need of genotyping. As a consequence, a cardiac expression-based enrichment strategy can be used to identify ZIC lines. Instead of screening 609 GBT lines, only 35 ZIC lines need to be screened, which significantly reduced the workload. We acknowledge that some genes with extremely weak cardiac expression might be missed; however, this is not a concern during the early phase of a genome-wide screen. Third, it is economically feasible to house hundreds of mutant fish lines with different genetic lesions to 1–3 years old. Finally, we optimized an ECG technology, defined the baseline SSS in WT fish, and implemented heat-stress to zebrafish at old ages, which shall increase the SSS phenotypic expressivity.
While the forward genetic screening approach has been successfully utilized to pinpoint genetic basis of cardiogenesis in embryonic fish and doxorubicin-induced cardiomyopathy (DIC) in adult zebrafish, Amsterdam et al., 1999; Ding et al., 2016 this study extended this powerful genetic approach to adult zebrafish for discovering genetic factors associated with rhythm disorders. Given very little knowledge of molecular underpinnings of SSS, the development of this novel approach is significant. Human genetics approach has been difficult, partially owing to the aging associated nature - SSS-like phenotypes at its early stage are often missed, because SA episodes cannot be detected if the ECG measurement only covers a short time window. It takes years in patients to develop from asymptotic to onset of SSS symptoms. Moreover, human genetic studies of SSS are typically confounded by complicated environmental factors, which are minimalized in our zebrafish forward genetic approach - each ZIC mutant is maintained in a well-controlled living environment, and the only difference among different ZIC lines is a single genetic deficiency.
DNAJB6 is a new SSS gene with a unique expression in SAN
The human DNAJB6 gene encodes a molecular chaperone protein of the heat shock protein 40 (Hsp40) family. DNAJB6 has been previously linked to neurodegenerative diseases via its function in protein folding and the clearance of polyglutamine stretches (polyQ), Gillis et al., 2013; Hageman et al., 2010 and to muscular dystrophy via its protein-protein interaction with Bag3 in the sarcomere Sarparanta et al., 2012. Our previous forward genetic screen in adult zebrafish identified GBT411/dnajb6b as a deleterious modifier for DIC Ding et al., 2016. Here, we provided several lines of evidence in both fish and mouse models, suggesting new functions of Dnajb6 as a genetic factor for arrhythmia/SSS. First, GBT411/dnajb6b is one of three ZIC lines with SSS-like phenotypes that were identified from a screen of 607 GBT lines that is independent of the previous DIC screen. Second, in zebrafish, the increased incidence of SA episodes and reduced heart rate, two main features of SSS, were detected in as early as 10-month-old GBT411/dnajb6b homozygous fish. Similarly, bradycardia and SA episodes were noted in Dnajb6+/- KO mice at 6 months old, when the structural remodeling in both left ventricular and atrial myocardium have not occurred yet, and echocardiography indices remained indistinguishable from their age-matched siblings. Depletion of Dnajb6 in mice manifests severer phenotypes than in zebrafish, probably because mouse has only one DNAJB6 homologue, while zebrafish has two DNAJB6 homologues, dnajb6b and dnajb6a. Third, consistent with loss-of-function studies, DNAJB6 expression was detected in the SAN of both zebrafish and mice. Importantly, DNAJB6 is highly enriched in the SAN region of the mouse comparing to the surrounding atrial tissue. Fourth, transcriptome analysis of Dnajb6+/- mice uncovered altered expression of genes involved in calcium handing, ion channels, and Wnt signaling pathway, which have been linked to the formation/function of the SAN during development. Thus, our data from mice strongly suggested that the observed SSS is not a consequence of Dnajb6 cardiomyopathy, Instead, the irregular heartbeat is most likely a direct consequence of Dnajb6 depletion in pacemaker cells, subsequently contributing to the pathogenesis of cardiomyopathy that occurs later. To ultimately confirm this hypothesis and to discern functions of Dnajb6 in SAN pacemakers from working cardiomyocytes, a tissue-specific KO line for Dnajb6 needs to be generated and studied. Prompted by our preliminary success in identifying potentially significant sequence variances for DNAJB6 from human SSS patients, future human genetic studies are warranted to search more sequence variants and to confirm their pathogenicity, which are required to firmly establish DNAJB6 as a new SSS causative gene in human.
Detailed examination of DNAJB6 expression in the SAN uncovered unique expression patterns. While the expression of DNAJB6 is detected in the SAN (Figure 3A), we found a partial co-expression with one of the main pacemaker protein HCN4: DNAJB6-positive cells overlap only with a portion of the HCN4-positive cells (Figure 3B). The unique expression pattern of DNAJB6 is also underscored by a negative correlation between the expression level of DNAJB6 and TBX3 (Figure 3C and D), as TBX3 is one of the main transcriptional regulators to define pacemaker cell specificity Hoogaars et al., 2007; Mohan et al., 2020 While these results may sound surprising, studies on isolated SAN cells reported dramatic variability in the density of HCN4-formed ‘funny’ current If Honjo et al., 1996; Mangoni and Nargeot, 2001; Monfredi et al., 2018; Wilders et al., 1996. In spontaneously beating cardiomyocytes isolated from the rabbit SAN, Monfredi et al. showed that If density can range from 0 to ~50 pA/pF and some the spontaneously beating SAN cells may have little to zero If Monfredi et al., 2018. The authors further observed SAN cells with lower If current densities, demonstrating a significantly greater sensitivity to inhibition of Ca2+ clock component of the SAN pacemaking machinery by cyclopiazonic acid, a moderate disruptor of Ca2+ cycling, in terms of beating rate slowing. The authors also noted that a relatively large cell population (21 of 90 cells) stopped beating when the sarcoplasmic reticulum pumping rate decreased in the presence of cyclopiazonic acid, despite a relatively high If density. Together with other studies, Boyett et al., 2000; Kim et al., 2018 these results may indicate a significant functional heterogeneity of pacemaker cells within the SAN in terms of their spontaneous beating rate, ion channel and calcium handling protein expression repertoire, and molecular mechanisms of their pacemaker activities. The latter was recently linked to the balance between the voltage and calcium components of the coupled-clock pacemaker system describing mechanisms of SAN automaticity Lakatta et al., 2010. As summarized in details in our recent review, Lang and Glukhov, 2021 it was suggested that pacemaker cells, which primary rely on the Ca2+ clock, are more sensitive to the autonomic modulation through cAMP-mediated regulation of phosphorylation of Ca2+ handling proteins Kim et al., 2018. This is in line with our findings indicating that DNAJB6 is mainly expressed in SAN cells with low HCN4 density (Figure 3B) and that Dnajb6 knock-out activates the expression of transcription factor TBX3 and affects calcium homeostasis genes (Figure 7) and leads to abnormal autonomic regulation of the SAN (Figure 4E and Figure 5).
Potential mechanisms underlying the role of Dnajb6 in SSS
Besides uncovering a crucial role of DNAJB6 in SAN automaticity and autonomic regulation and specification of SAN pacemaking, our studies raised several hypotheses on the underlying cellular and molecular mechanisms. Our model-based analysis provided pilot screening of potential cellular/ionic targets that could contribute to the observed SSS phenotype in DNAJB6+/- mice. Direct testing of these mechanisms would require a substantial amount of single SAN cell patch clamp and confocal microscopy experiments that can be further pursued in a follow-up study. Importantly, these new in silico experiments add another conceptual level to a phenotype-based high-throughput screening approach introduced in the current study to identify genetic factors associated with SAN dysfunction.
In addition to SSS phenotype, we observed enhanced ectopic activity in the Dnajb6+/- mice that was associated with subsidiary atrial pacemakers (Figure 5A and B). Based on the diffused AV canal signal and SAN signal loss in the GBT411/dnajb6b homozygous mutant fish hearts (Figure 2D), as well as a negatively correlated expression of DNAJB6 with TBX3 in the mouse SAN tissues (Figure 3C and D), we speculate that DNAJB6 might act as a suppressor of TBX3 transcription factor to define SAN cell specification. This potential mechanism is also supported by the observation in mice that loss-of-function of Dnaj6b results in conduction system defects and ectopic pacemaker activity. Since TBX3 suppresses chamber myocardial differentiation, Bakker et al., 2008 upregulation of TBX3 may thus contribute to enhanced atrial ectopic activity observed in Dnajb6+/-mic. Furthermore, TBX3 has been recently identified as component of the Wnt/β-catenin-dependent transcriptional complex, Zimmerli et al., 2020 which is significantly affected in Dnajb6+/- mice (Figure 6B). This further indicates a possible role of TBX3 in both SAN and atrial remodeling.
In the human hearts, all the observed ectopic pacemakers were located within the region of extensive distributed system of atrial pacemakers (i.e. atrial pacemaker complex), which includes but extends well beyond an anatomically defined SAN Boineau et al., 1988. Under physiological conditions, spontaneous activity of subsidiary pacemakers is overdrive suppressed by the SAN. However, when SAN function is diminished (i.e. during SSS), subsidiary pacemakers can produce escape beats leading to pacemaker irregularities and significant heart rate lability. Though the subsidiary pacemakers can provide a relatively regular rhythm, they are characterized by a slower resting heart, slower exertional heart rates, a prolonged post-pacing recovery time (a parameter similar to SAN recovery time but for non-SAN pacemakers), and an increased beat-to-beat heart rate variability Morris et al., 2013. Furthermore, the electrical activity of this subsidiary pacemakers is more akin to that of the SAN than to the surrounding atrial muscle; the subsidiary pacemaker action potential exhibits prominent diastolic depolarization and a significantly lower maximum diastolic potential, take-off potential, overshoot, rate of rise, and amplitude than typical atrial muscle Rozanski et al., 1984. Finally, while being bradycardic in general, subsidiary atrial pacemakers can contribute to the development of atrial tachycardia Kistler et al., 2006. Therefore, we hypothesize that TBX3 overexpression observed in Dnajb6+/- mice, could further facilitate pacemaker activity in cells within the extended atrial pacemaker complex and, maybe, promote atrial arrhythmogenesis in the setting of profound structural remodeling.
A phenotype-based screening approach would facilitate the elucidation of molecular basis of SSS
Besides dnajb6b, our pilot forward genetic screen also suggested two additional candidate SSS genes like cyth3a and vapal, pending more experimental evidence to confirm their function. This forward genetic screening approach is scalable to the genome, which would generate a comprehensive list of candidate genes for SSS. Because there are at least three major cell types in the SAN region, including pacemaker cells in SAN that generate rhythm, paranodal areas and transition cells in the atrium that transmit the signal from pacemaker cells to govern coordinated contraction of the heart from atrium and then to the ventricle, Li et al., 2017; Li et al., 2020 newly identified SSS genes could be categorized into different groups based on their expression pattern and phenotypic traits. We anticipate that systematic studies of these candidate genes identified from zebrafish will significantly advance our understanding of pathophysiology of SSS.
Materials and methods
Animals
All experiments were conducted in accordance with the Guidelines for the Care and Use of Laboratory Animals published by the US National Institutes of Health (publication No. 85–23, revised 1996). All animal procedures and protocols used in these studies (for zebrafish, #: A00005409-20; for mouse, #: A00003511-20 and M005490-R02) have been approved by the Mayo Clinic Institutional Animal Care and Use Committee (Permit number: D16-00187) and by the Animal Care and Use Committee of University of Wisconsin-Madison (Permit number: D16-00239). The zebrafish (Danio rerio) WIK line was maintained under a 14 hr light/10 hr dark cycle at 28.5 °C. All GBT lines were generated previously Clark et al., 2011; Ichino et al., 2020; Ding et al., 2013. The Dnajb6 knockout (KO) mice, originally named Dnajb6tm1.1(KOMP)Vlcg, were generated from the Jackson Laboratory (Original catalog #018623). Briefly, the insertion of Velocigene cassette ZEN-Ub1 created a deletion sized 36,843 bp nucleotides spanning from the first to the last intron of the Dnajb6 gene at the Chromosome 5 (Genome Build37) of the C57BL/6 N mice. The mouse was subsequently bred to a ubiquitous Cre deletion mouse line for recombination of the LoxP sites that excised the neomycin selection cassette. The following genotyping PCR primers for the Dnajb6 mutant mice were used: mutant primer F2, 5’-AAACTGCGCACTGTACCACC-3’ and mutant primer R2, 5’-CGGTCGCTACCATTACCAGT-3’ for detecting the mutant allele (predicted size of 700 bp); and WT primer F1, 5’-TACTCCAGCCCCACTCTTACTC-3’ and WT primer R1, 5’- ACTGCCCATCTTCTTCAACTTC-3’ for detecting the WT allele (predicted size of 300 bp).
Enrichment and cloning of 35 ZIC mutants
Request a detailed protocolZebrafish cardiac insertional (ZIC) mutants were identified and collected based on the mRFP expression in the embryonic heart from 2 to 4 days post-fertilization (dpf) and/or in the dissected adult heart at 6 months to 1 year of age. All ZIC lines, each with a single copy of GBT insertion, were obtained after 2–4 generations of outcrosses, guided by Southern blotting using the GFP probe primed to the GBT vector Ding et al., 2013. A combination of three different methods including Inverse PCR, 5’-RACE and/or 3’-RACE were employed to clone the GBT transposon integration sites accordingly to previously published protocols Clark et al., 2011; Ichino et al., 2020; Ding et al., 2013. A combination of gene-specific primers flanking the GBT integration site coupled with GBT vector-specific primer were used for genotyping PCR to identify homozygous mutants for the three candidate ZIC including GBT103/cyth3a, GBT410/vapal, and GBT411/dnajb6b lines using genomic DNA isolated from tail fins Clark et al., 2011; Ichino et al., 2020; Ding et al., 2013.
Zebrafish electrocardiogram (ECG)
Request a detailed protocolMicrosurgery was operated under a dissection microscope to remove the silvery epithelial layer of the hypodermis one week before fish were subjected to the ECG Yan et al., 2020. Fish were initially acclimated for 1 hr after transferred from the circulating fish facility to the laboratory bench, followed by anesthesia in the solution of pH 7.0 adjusted tricaine (MS-222, Sigma) at the concentration of 0.02% dissolved in E3 medium (containing 5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, and 0.33 mM MgSO4) for 6 min. Two minutes of ECG recording were then obtained with the ECG recording system, according to the instructions (ZS-200, iWorx Systems, Inc) and a recently published protocol Yan et al., 2020. Initial ECG screens of ZIC heterozygous mutants were performed at 32 °C using a temperature-controlled chamber set-up, made by covering the ECG recording system with a foam box. 6 to 25 fish per ZIC line were initially analyzed, depending on the fish availability. The ECG machine was held on top of a heating plate controlled by a heating machine. The subsequent ECG validation in the homozygous mutants was performed at room temperature (25 °C). To analyze the ECG recording, ECG signals were amplified and filtered at 0.5 Hz high pass and 200 Hz low pass. ECG variables, including heart rate, PR interval, QRS duration, QT interval and R-wave amplitude, and PP and RR intervals were calculated using an in-house Matlab code Lenning et al., 2017. A SA episode was defined in zebrafish when the PP interval is more than 1.5 s.
Mouse ECG and echocardiography
Request a detailed protocolMouse echocardiography and ECG measurements were performed according to a previously published protocol with modifications Ding et al., 2016; Ding et al., 2020b. For ECG, mice were anesthetized with isoflurane (0.5%–1.0% v/v) via a nose cone. Mice were placed on an ECG-heater board with 4 paws on individual electrodes. The ECG-heater board maintained the body temperature at 37 °C. The ECG signal was amplified through an amplifier (Axon CNS digital 1440 A) and recorded using > Chart 5 software. For each mouse, 10 min of ECG signal were recorded. Series of ECG parameters, including heart rate, PR interval, QRS duration, QT interval and RR interval were calculated by an in-house Matlab code Lenning et al., 2017. For echocardiography, mice were anesthetized under light isoflurane (0.5%–1.0% v/v) administered via a nose cone. Echocardiography gel was placed on the shaved chest, and the mouse heart was imaged with a 13-MHz probe using two-dimensional echocardiography (GE Healthcare). All measurements were made by an independent operator to whom the study groups were masked.
Administration of autonomic response drugs
Request a detailed protocolFor zebrafish, 0.6 µg/g isoproterenol (Millipore Sigma, Cat# 1351005), 4 µg/g atropine (Millipore Sigma, Cat# A0132), and 0.3 µg/g carbachol (Millipore Sigma, cat# C4382) were administrated via intraperitoneal injection. For in vivo mouse studies, 0.2 mg/kg isoproterenol, 1 mg/kg atropine, and 0.3 mg/kg carbachol was injected intraperitoneally. For ex vivo mouse atrial studies, 100 nM isoproterenol, 2 µM atropine, and 300 nM carbachol was administered via superfusion for 10–20 min.
Antibody immunostaining
Request a detailed protocolHeart samples harvested from mouse SAN tissues were embedded in a tissue freezing medium, followed by sectioning at 10 μm using a cryostat (Leica CM3050 S). The slides were subjected to immunostaining using a previously described protocol Sun et al., 2009. The following antibodies were used: anti-HCN4 (Millipore, Cat#: AB5805; Novus biologicals, Cat#: NB100-74439) at 1:200, anti-DNAJB6 (Novus, Cat#: H00010049-M0; Santa Cruz Biotechnology Inc, Cat#: sc-104204) at 1:200, anti-TBX3 (abcam, Cat#: ab99302). All images were captured either using a Zeiss Axioplan II microscope equipped with ApoTome and AxioVision software (Carl Zeiss Microscopy) or a Zeiss LSM 780 confocal microscope. Signal intensity from DNAJB6 and TBX3 antibodies immunostaining was quantified using Zen 2.3 Pro software.
Western blotting
Request a detailed protocolFor Western blotting, mice embryonic hearts at E12.5 stage were dissected after genotyping PCR using genomic DNA isolated from tail and transferred to RIPA lysis buffer supplemented with complete protease inhibitor cocktail for homogenization. About 1 µg resultant protein lysates were subject to western blotting using a standard protocol. The following primary antibodies were used: anti-Gapdh (1:4000, Santa Cruz Biotechnology Inc, Cat#: sc-25778); anti-DNAJB6 (1:6000, abcam, Cat#: 198995).
Isolated mouse atrial preparations
Request a detailed protocolThe mouse atrial preparation was performed as previously described Glukhov et al., 2015. After the mice were anesthetized with isoflurane, a mid-sternal incision was applied. The heart was then removed and cannulated to a custom made 21-gauge cannula. The heart was then perfused and superfused with oxygenated (95% O2, 5% CO2), 37 °C modified Tyrode solution (in mM: 128.2 NaCl, 4.7 KCl, 1.19 NaH2PO4, 1.05 MgCl2, 1.3 CaCl2, 20.0 NaHCO3, and 11.1 glucose; pH = 7.35 ± 0.05). Lung, thymus, and fat tissue was then removed. Perfusion was maintained under constant aortic pressure of 60–80 mmHg. After 10 min stabilization, the ventricles were dissected. The atrial were cut open as previously described Lang and Glukhov, 2016. The medial limb of the crista terminalis was cut to open right atrial appendage. The preparation was superfused with Tyrode solution at a constant rate of ~15 ml/min.
Optical mapping
Request a detailed protocolHigh spatial (100x100 pixels, 60±10 μm per pixel) and temporal (1,000–3,000 frames/sec) resolution optical mapping of electrical activity was applied on the isolated mouse atrial preparations as previously described Lang and Glukhov, 2016; Lang et al., 2011. The isolated mouse atrial preparations were coronary and surface stained with voltage-sensitive dye RH-237 (1.25 mg/ml in dimethyl sulfoxide ThermoFisher Scientific, USA). Blebbistatin (10 μM, Tocris Bioscience, USA) was then applied to reduce the motion artifact. A 150 W halogen lamp (MHAB-150W, Moritex USA Inc, CA, USA) with band pass filter (530/40 nm) was used as excitation light source. The fluorescent light emitted from the preparation was recorded by a MiCAM Ultima-L camera (SciMedia, CA, USA) after a long-pass filter (>650 nm). The acquired fluorescent signal was digitized, amplified, and visualized using custom software (SciMedia, CA, USA). After 20–30 min stabilization, activation map was collected during baseline spontaneous rhythm. To estimate the pacemaker location and a possible pacemaker shift during autonomic stimulation, the originations of action potentials were plotted with orthogonal axes crossing at the inferior vena cava. The superior to inferior direction is along the ordinate. The lateral to media direction is along the abscissa. SAN recovery time (SANRT) was measured as the time-period between the last S1S1 pacing (10 Hz) beat and the first spontaneous beat. Corrected SANRT (SANRTc) was calculated as the difference between the SANRT and the resting cycle length measure before the SANRT pacing protocol. After baseline measurement, 100 nM isoproterenol was applied. Recordings were collected after 10 min which allows the stimulation to reach steady-state effect. Complete washout was then performed which is characterized by the recovery of the heart rate back to baseline values. Additional staining and blebbistatin was applied as needed. 300 nM carbachol then was applied. 2 µM atropine was used after protocols completed during carbachol stimulation.
RNA-seq data collection and analysis
Request a detailed protocolTotal RNA was extracted from the right atrium (RA) tissues of 1-year-old Dnajb6+/- heterozygous mutant hearts and WT sibling controls. Six total samples ( +/- biological replicates for each genotype) were submitted for RNA sequencing (Azenta Life Science, NJ). Genes were considered to be differentially expressed between the two groups if they exhibited a greater than twofold change and an FDR of less than 0.05 according to the DESeq approach Love et al., 2014. Unsupervised hierarchical clustering was performed with Pearson correlation and scaled based on the fragments per kilobase of transcript per million mapped reads (FPKM) value using the pheatmap R package (https://github.com/raivokolde/pheatmap; Kolde et al., 2018 R Development Core Team, 2022). The gene lists of interest were annotated by IPA (QIAGEN) (http://www.ingenuity.com/). We queried the IPA with the gene list of interest to map and generate putative biological processes/functions, networks, and pathways based on the manually curated knowledge database of molecular interactions extracted from the public literature. The enriched pathways and gene networks were generated using both direct and indirect relationships/connectivity. These pathways and networks were ranked by their enrichment score, which measures the probability that the genes were included in a network by chance.
Quantitative reverse transcription (RT) PCR
Request a detailed protocolTotal RNA was extracted from ~2 mg of right atrium (RA) tissues of 1-year-old Dnajb6+/- heterozygous mutant hearts and WT sibling controls using Trizol reagent (ThermoFisher Scientific) following the manufacturer’s instruction. Approximately+/- µg total RNA was used for reverse transcription (RT) and cDNA synthesis using Superscript III First-Strand Synthesis System (ThermoFisher Scientific). Real-time quantitative RT-PCR was run in 96-well optical plates (ThermoFisher Scientific) using an Applied Biosystem VAii 7 System (ThermoFisher Scientific). Gene expression levels were normalized using the expression level of glyceraldehyde 3-phosphate dehydrogenase (gapdh) by –ΔΔCt (cycle threshold) values. All quantitative RT-PCR primer sequences were listed in Supplementary file 5.
Histology
Request a detailed protocolFor H&E and Masson’s trichrome staining of Dnajb6+/- mice left ventricle, mouse hearts were dissected and harvested after mice were euthanized by administration of high-dose (5%) isoflurane anesthesia and after ventilation was ceased. Dissected mice hearts were immediately fixed in 4% PBS buffered formaldehyde overnight at 4 °C and sent to the Mayo Clinic Histology Core Laboratory for sample processing and H&E staining.
For transmission electron microscopy (TEM) analysis, the left ventricle apexes of dissected hearts from either zebrafish or mice were fixed immediately in Trump’s solution (4% paraformaldehyde and 1% glutaraldehyde in 0.1 M phosphate buffer [pH 7.2]) at room temperature for 1 hr, followed by overnight incubation at 4 °C. Fixed samples were subsequently processed and imaged at the Mayo Clinic Electron Microscopy Core Facility using a Philips CM10 transmission electron microscope.
To quantify the amount of fibrosis in Dnajb6+/- mouse atria, the isolated atrial preparations were fixed after optical mapping experiments overnight in +/-% paraformaldehyde buffered with 0.1 M sodium phosphate, pH 7.4; and then paraffin embedded. The preparations were sectioned parallel to the epicardial surface at 3–5 μm thickness. Tissue sections were mounted on Superfrost Plus glass slides (Fisher Scientific, Pittsburgh, PA) and maintained at room temperature until use. Sections were stained for histology with Masson’s trichrome (International Medical Equipment, San Marcos, CA, USA) and Picrosirius Red (International Medical Equipment, San Marcos, CA, USA). The density of fibrosis was estimated as a ratio of cardiac tissue to connective tissue measured at different transmural layers, quantified using ImageJ (National Institutes of Health) as previously described Glukhov et al., 2015.
Computational modeling
Request a detailed protocolTo investigate the cellular mechanisms underlying the SSS phenotype, we used our model of the mouse SAN myocyte, Morotti et al., 2021 based on the original model, Kharche et al., 2011 and including the formulation of the acetylcholine-activated K+ current developed by Arbel-Ganon et al. to simulate carbachol administration Arbel-Ganon et al., 2020. Functional effects of isoproterenol administration on ion channels and transporters (listed in Supplementary file 6) were simulated as in the parent model, Kharche et al., 2011 wherein properties of isoproterenol-dependent modulation of voltage-gated Ca2+ currents and funny current If were updated to reflect experimental observations in mice Larson et al., 2013; Peters et al., 2021. Using an established approach, Sobie, 2009 we randomly varying selected model parameters describing maximum ion channel conductances and ion transport rates (defined in Supplementary file 6) to generate a population of 10,000 model variants. For each variant, the baseline value of each parameter was independently varied with a log-normal distribution (σ=0.26). We assessed the steady-state firing rate in each model in the population at baseline, and upon stimulation with either isoproterenol or carbachol. Model variants showing non-physiological behavior (e.g., lack of firing activity) at baseline or in response to autonomic stimulation were discarded from subsequent analysis. We separated the population of models in two subpopulations mimicking the WT and Dnajb6+/- mice phenotypes. Namely, we extracted the model variants that recapitulate changes observed in Dnajb6+/- vs. WT mice, including a slower firing rate at baseline, an increased response to isoproterenol, and a diminished response to carbachol administration . We analyzed the parameter value differences in these two subgroups to reveal several ionic processes that are significantly correlated with the observed electrophysiological changes. The nonparametric two-sided Wilcoxon rank sum test was used to compare the two groups and p value less than 0.05 was considered statistically significant. All the codes used to perform in silico simulations were generated in MATLAB (MathWorks, Natick, MA, USA) and are freely available for download at http://elegrandi.wixsite.com/grandilab/downloads and https://github.com/drgrandilab/Ding-et-al-2022-mouse-sinoatrial-model; Grandi Lab, 2022 (copy archived at swh:1:rev:9ffd9fee426ef9e4b26826b4b8700a93821ba9ab).
Statistics
No sample sizes were calculated before performing the experiments. No animals were excluded for analysis. Unpaired two-tailed student’s t-test was used to compare two groups. One-way Analysis of Variance (ANOVA) or Kruskal-Wallis test followed by post hoc Tukey’s test was used for comparing three and more groups. Chi-square test was used for rate comparison. p Values less than 0.05 was considered statistically significant. For dot plot graphs, values are displayed as mean ± standard deviation (SD). Sample size (N) represents animal number, otherwise specifically designated as biological or technical replicates. All statistical analyses were conducted with the Graphpad Prism 7 and/or R Statistical Software Version 3.6.1.
Availability of the materials and resources
Request a detailed protocolAll reagents are available upon reasonable request. Zebrafish GBT mutant lines are available either from the Zebrafish International Recourse Center (ZIRC, http://zebrafish.org) or the Mayo Clinic Zebrafish Facility, respectively. Both RNAseq raw data and processed data are deposited to GEO (Access number GSE195953) associated with the token: kvmhesayzryljop. The code of our computational model of the mouse SAN myocyte is freely available for download at http://elegrandi.wixsite.com/grandilab/downloads andhttps://github.com/drgrandilab/Ding-et-al-2022-mouse-sinoatrial-model.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Source data files have been provided for Figure 4 and Figure 7.
References
-
A large-scale insertional mutagenesis screen in zebrafishGenes & Development 13:2713–2724.https://doi.org/10.1101/gad.13.20.2713
-
Genetics of sick sinus syndromeCardiac Electrophysiology Clinics 2:499–507.https://doi.org/10.1016/j.ccep.2010.09.001
-
Distinct mechanisms mediate pacemaker dysfunction associated with catecholaminergic polymorphic ventricular tachycardia mutations: insights from computational modelingJournal of Molecular and Cellular Cardiology 143:85–95.https://doi.org/10.1016/j.yjmcc.2020.04.017
-
Zebrafish as a model to study cardiac development and human cardiac diseaseCardiovascular Research 91:279–288.https://doi.org/10.1093/cvr/cvr098
-
The sinoatrial node, a heterogeneous pacemaker structureCardiovascular Research 47:658–687.https://doi.org/10.1016/s0008-6363(00)00135-8
-
Sick sinus syndromeCardiac Electrophysiology Clinics 10:183–195.https://doi.org/10.1016/j.ccep.2018.02.002
-
Modeling inherited cardiomyopathies in adult zebrafish for precision medicineFrontiers in Physiology 11:599244.https://doi.org/10.3389/fphys.2020.599244
-
Knockout of SORBS2 protein disrupts the structural integrity of intercalated disc and manifests features of arrhythmogenic cardiomyopathyJournal of the American Heart Association 9:e017055.https://doi.org/10.1161/JAHA.119.017055
-
The DNAJB6 and DNAJB8 protein chaperones prevent intracellular aggregation of polyglutamine peptidesThe Journal of Biological Chemistry 288:17225–17237.https://doi.org/10.1074/jbc.M112.421685
-
Functional anatomy of the murine sinus node: high-resolution optical mapping of ankyrin-B heterozygous miceAmerican Journal of Physiology. Heart and Circulatory Physiology 299:H482–H491.https://doi.org/10.1152/ajpheart.00756.2009
-
Little fish, big data: zebrafish as a model for cardiovascular and metabolic diseasePhysiological Reviews 97:889–938.https://doi.org/10.1152/physrev.00038.2016
-
A rare variant in MYH6 is associated with high risk of sick sinus syndromeNature Genetics 43:316–320.https://doi.org/10.1038/ng.781
-
Correlation between electrical activity and the size of rabbit sino-atrial node cellsThe Journal of Physiology 496 (Pt 3):795–808.https://doi.org/10.1113/jphysiol.1996.sp021728
-
TBX3 controls the sinoatrial node gene program and imposes pacemaker function on the atriaGenes & Development 21:1098–1112.https://doi.org/10.1101/gad.416007
-
Genome-Wide identification of mouse congenital heart disease lociHuman Molecular Genetics 19:3105–3113.https://doi.org/10.1093/hmg/ddq211
-
A mathematical model of action potentials of mouse sinoatrial node cells with molecular basesAmerican Journal of Physiology. Heart and Circulatory Physiology 301:H945–H963.https://doi.org/10.1152/ajpheart.00143.2010
-
Frequency of cardiac rhythm abnormalities in a half million adultsCirculation. Arrhythmia and Electrophysiology 11:e006273.https://doi.org/10.1161/CIRCEP.118.006273
-
P-wave morphology in focal atrial tachycardia: development of an algorithm to predict the anatomic site of originJournal of the American College of Cardiology 48:1010–1017.https://doi.org/10.1016/j.jacc.2006.03.058
-
Optical mapping of action potentials and calcium transients in the mouse heartJournal of Visualized Experiments 2011:3275.https://doi.org/10.3791/3275
-
High-resolution optical mapping of the mouse sino-atrial nodeJournal of Visualized Experiments 1:54773.https://doi.org/10.3791/54773
-
Cellular and molecular mechanisms of functional hierarchy of pacemaker clusters in the sinoatrial node: new insights into sick sinus syndromeJournal of Cardiovascular Development and Disease 8:43.https://doi.org/10.3390/jcdd8040043
-
Redundant and diverse intranodal pacemakers and conduction pathways protect the human sinoatrial node from failureScience Translational Medicine 9:eaam5607.https://doi.org/10.1126/scitranslmed.aam5607
-
From genotype to phenotype: a primer on the functional follow-up of genome-wide association studies in cardiovascular diseaseCirculation. Genomic and Precision Medicine 11:e001946.
-
Organisation of the mouse sinoatrial node: structure and expression of HCN channelsCardiovascular Research 73:729–738.https://doi.org/10.1016/j.cardiores.2006.11.016
-
The Ensembl variant effect predictorGenome Biology 17:122.https://doi.org/10.1186/s13059-016-0974-4
-
The 11th world survey of cardiac pacing and implantable cardioverter-defibrillators: calendar year 2009 -- a world Society of arrhythm’a's projectPacing and Clinical Electrophysiology 34:1013–1027.https://doi.org/10.1111/j.1540-8159.2011.03150.x
-
Sick sinus syndrome and atrial fibrillation in older persons-a view from the sinoatrial nodal myocyteJournal of Molecular and Cellular Cardiology 83:88–100.https://doi.org/10.1016/j.yjmcc.2015.02.003
-
Electrophysiological heterogeneity of pacemaker cells in the rabbit intercaval region, including the SA node: insights from recording multiple ion currents in each cellAmerican Journal of Physiology. Heart and Circulatory Physiology 314:H403–H414.https://doi.org/10.1152/ajpheart.00253.2016
-
Intracellular na+ modulates pacemaking activity in murine sinoatrial node myocytes: an in silico analysisInternational Journal of Molecular Sciences 22:5645.https://doi.org/10.3390/ijms22115645
-
Development of the cardiac conduction system in zebrafishGene Expression Patterns 21:89–96.https://doi.org/10.1016/j.gep.2016.08.003
-
SoftwareR: A language and environment for statistical computing, version 2.6.2R Foundation for Statistical Computing, Vienna, Austria.
-
Alterations in subsidiary pacemaker function after prolonged subsidiary pacemaker dominance in the canine right atriumJournal of the American College of Cardiology 4:535–542.https://doi.org/10.1016/s0735-1097(84)80098-4
-
Pacemaker channel dysfunction in a patient with sinus node diseaseThe Journal of Clinical Investigation 111:1537–1545.https://doi.org/10.1172/JCI16387
-
Benefits and limitations of genome-wide association studiesNature Reviews. Genetics 20:467–484.https://doi.org/10.1038/s41576-019-0127-1
-
Genetic insight into sick sinus syndromeEuropean Heart Journal 42:1959–1971.https://doi.org/10.1093/eurheartj/ehaa1108
-
Pacemaker activity of the human sinoatrial node: an update on the effects of mutations in HCN4 on the hyperpolarization-activated currentInternational Journal of Molecular Sciences 16:3071–3094.https://doi.org/10.3390/ijms16023071
-
Model clamp and its application to synchronization of rabbit sinoatrial node cellsThe American Journal of Physiology 271:H2168–H2182.https://doi.org/10.1152/ajpheart.1996.271.5.H2168
-
Genetic analysis of sick sinus syndrome in a family harboring compound CACNA1C and TTN mutationsMolecular Medicine Reports 17:7073–7080.https://doi.org/10.3892/mmr.2018.8773
Article and author information
Author details
Funding
Mayo Foundation for Medical Education and Research
- Xiaolei Xu
Science and Technology Innovation Action Plan of Shanghai (201409005600)
- Yigang Li
National Institute of Health (GM063904)
- Stephen C Ekker
American Heart Association (16SDG29120011)
- Alexey V Glukhov
National Institute of Health (R00HL138160)
- Stefano Morotti
National Institute of Health (R01HL131517)
- Eleonora Grandi
American Heart Association (15SDG24910015)
- Eleonora Grandi
American Heart Association (20POST35120462)
- Haibo Ni
National Institute of Health (NIH R01HL141214)
- Alexey V Glukhov
Wisconsin Partnership Program 4140
- Alexey V Glukhov
American Heart Association (17POST33370089)
- Di Lang
American Heart Association (846898)
- Di Lang
National Institute of Health (1OT2OD026580-01)
- Eleonora Grandi
National Institute of Health (P01HL141084)
- Eleonora Grandi
National Institutes of Health (R44OD024874)
- Tai Le
National Institutes of Health (R44OD024874)
- Hung Cao
National Institutes of Health (R01HL142704)
- Jun Wang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Beninio Gores and Kashia Stragey for managing zebrafish facility and Ronald H May for murine echocardiography. This work was supported in part by Mayo Foundation to XX, grant from Science and Technology Innovation Action Plan of Shanghai, experimental animal research project 201409005600 to YGL, NIH GM063904 to SCE, NIH R01HL141214, American Heart Association 16SDG29120011, and the Wisconsin Partnership Program 4140 to AVG, American Heart Association 17POST33370089 and American Heart Association Career Development Award 846898 to DL, NIH R00HL138160 to SM, NIH R01HL131517, NIH P01HL141084 NIH Stimulating Peripheral Activity to Relieve Conditions Grant 1OT2OD026580-01, UC Davis School of Medicine Dean’s Fellow Award, and American Heart Association Scientist Development Award 15SDG24910015 to EG, and American Heart Association Postdoctoral Fellowship 20POST35120462 to HN, NIH R44OD024874 to TL and HC, NIH R01HL142704 to JW.
Ethics
All experiments were conducted in accordance with the Guidelines for the Care and Use of Laboratory Animals published by the US National Institutes of Health (publication No. 85-23, revised 1996). All animal procedures and protocols used in these studies (for zebrafish, #: A00005409-20; for mouse, #: A00003511-20 and M005490-R02) have been approved by the Mayo Clinic Institutional Animal Care and Use Committee (Permit number: D16-00187) and by the Animal Care and Use Committee of University of Wisconsin-Madison (Permit number: D16-00239).
Copyright
© 2022, Ding et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 977
- views
-
- 192
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Genetics and Genomics
Calcineurin is a highly conserved calcium/calmodulin-dependent serine/threonine protein phosphatase with diverse functions. Inhibition of calcineurin is known to enhance the lifespan of Caenorhabditis elegans through multiple signaling pathways. Aiming to study the role of calcineurin in regulating innate immunity, we discover that calcineurin is required for the rhythmic defecation motor program (DMP) in C. elegans. Calcineurin inhibition leads to defects in the DMP, resulting in intestinal bloating, rapid colonization of the gut by bacteria, and increased susceptibility to bacterial infection. We demonstrate that intestinal bloating caused by calcineurin inhibition mimics the effects of calorie restriction, resulting in enhanced lifespan. The TFEB ortholog, HLH-30, is required for lifespan extension mediated by calcineurin inhibition. Finally, we show that the nuclear hormone receptor, NHR-8, is upregulated by calcineurin inhibition and is necessary for the increased lifespan. Our studies uncover a role for calcineurin in the C. elegans DMP and provide a new mechanism for calcineurin inhibition-mediated longevity extension.
-
- Evolutionary Biology
- Genetics and Genomics
The evolutionary origins of Bilateria remain enigmatic. One of the more enduring proposals highlights similarities between a cnidarian-like planula larva and simple acoel-like flatworms. This idea is based in part on the view of the Xenacoelomorpha as an outgroup to all other bilaterians which are themselves designated the Nephrozoa (protostomes and deuterostomes). Genome data can provide important comparative data and help understand the evolution and biology of enigmatic species better. Here, we assemble and analyze the genome of the simple, marine xenacoelomorph Xenoturbella bocki, a key species for our understanding of early bilaterian evolution. Our highly contiguous genome assembly of X. bocki has a size of ~111 Mbp in 18 chromosome-like scaffolds, with repeat content and intron, exon, and intergenic space comparable to other bilaterian invertebrates. We find X. bocki to have a similar number of genes to other bilaterians and to have retained ancestral metazoan synteny. Key bilaterian signaling pathways are also largely complete and most bilaterian miRNAs are present. Overall, we conclude that X. bocki has a complex genome typical of bilaterians, which does not reflect the apparent simplicity of its body plan that has been so important to proposals that the Xenacoelomorpha are the simple sister group of the rest of the Bilateria.