2.1 cftr mutation leads to defects in cardiogenesis beginning at the stage where cardiac progenitor cells emerge
To determine whether CFTR plays an important role in cardiac development during zebrafish embryogenesis, we first performed morphological phenotypic analysis on two existing cftr mutant zebrafish lines, cftrscu101 [18] (http://zfin.org/action/feature/view/ZDB-ALT-181127-1) and cftrscu102 [19] (http://zfin.org/action/feature/view/ZDB-ALT-190307-1). cftrscu101 carried a 9-bp deletion, which was predicted to delete the start codon of Cftr; cftrscu102 carried a 2-bp deletion in Exon 6, causing frameshifts leading to premature stop codons at 219 AA. This revealed that embryos from the two cftr mutant lines similarly develop pericardial edema by 60 hpf (Fig. 1A), accompanied by a lower heart rate (Fig. 1B) and a dramatic reduction in systolic function (Fig. 1C) as measured by ventricular fractional shortening (cftrscu102data not shown).
The cftrscu101 mutant line was used for the following analyses. To further determine the effect of cftr mutation on heart development, we investigated the expression levels of several common cardiac markers by ISH in cftrscu101 embryos. At the 8-somite stage (13 hpf), when cardiogenic differentiation is initiated, the expression of the cardiac mesoderm marker nkx2.5 was dramatically affected (Fig. 1D). Staining for the two populations of nkx2.5+ cardiac precursors, which lie adjacent to the midline, was not only reduced in intensity (Fig. 1E), but their distance from one another was also increased (Fig. 1F). These data suggest that the specification and migration of cardiac mesoderm is greatly impaired in cftr mutants.
To further investigate the effect of cftr mutation on cardiac morphogenesis in zebrafish, cardiomyocyte marker myl7 (formerly indicated as cmlc2) expression was examined by ISH. In WT embryos, the heart tube had formed at 24 hpf, that is, cardiac progenitors were specified correctly and extended from the midline to the region under the left eye to form the linear heart tube. However, consistent with a previously described role for cftr in the establishment of embryonic laterality [20], cftrscu101 had heart tubes located at the midline and to the right, resulting in the formation of a defective heart tube by 48 hpf (Fig. 1G). Because of the variation existence of defective heart phenotype among larvae, the angle between the longitudinal axes of the ventricle and atrium was measured to demonstrate the looping defect. Results indicated that the statistical looping angle (average 3°) of cftrscu101 was decreased significantly when compared to WT embryos (average 21°) (Fig. 1H). Furthermore, cftrscu101 embryos also showed markedly dilated atria (Fig. 1I).
We also assessed the expression of the atrioventricular boundary marker genes bmp4 (Fig. 1J and K) and has2 (Fig. 1L and M) in the cftr mutants, and found that the expression of these two marker genes was also dramatically weakened in cftrscu101 embryos. Taken together, these results suggest that cftr plays an essential role in early heart morphogenesis, particularly at the stage of cardiogenic differentiation.
2.2 The blastula stage is the key temporal window of cftr function on cardiogenesis
To further characterize the role of cftr ion channel function during cardiogenesis, we investigated the effect of the CFTR inhibitor CFTRinh-172 on WT embryos at different time points. CFTRinh-172 is a potent and specific inhibitor of the CFTR channel that was identified by high-throughput screening [21, 22]. During the blastula stage (5 hpf), cardiac progenitor cells are located bilaterally in the lateral marginal zone [16]. We first treated WT embryos at 3–5 hpf with 10 µm CFTRinh-172, a dose widely used to inhibit cftr ion channel activity [23]. Interestingly, the embryos treated by CFTRinh-172 at this time point showed cardiac dysplasia similar to cftrscu101 (Fig. 2A-H).
Next, we treated WT embryos at the gastrula stage, when cardiac progenitor cells move dorsally towards the mid-line before migrating to the anterior lateral plate mesoderm (ALPM) [16]. However, expression and localization of the cardiac marker nkx2.5 did not show significant abnormality when embryos were treated at either the early (treatment at 6–8 hpf) and late (treatment at 8–10 hpf) gastrula stages (Supplemental Fig. 1A-F). In brief, these results indicate that the blastula stage is the stage at which cftr affects cardiogenesis.
2.3 RNA-seq analysis reveals significant changes in the expression levels of factors essential for cardiogenesis in cftrscu101 embryos at the blastula stage
To uncover the molecular mechanisms of cftr regulation during cardiogenesis, and to quantify the dynamic changes across the entire transcriptome, we performed RNA-seq analysis in WT and cftr mutant zebrafish embryos at the blastula stage (5 hpf). Results showed that cftr mutants differentially expressed 1334 transcripts (412 up-regulated and 922 down-regulated) (Supplemental Table 1). To characterize the function of these altered mRNAs, KEGG pathway-based classification analysis on the ontology of cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems was performed, and revealed a widely different distribution in cftrscu101 mutant embryos compared to WT (Fig. 3A and Supplemental Table 2).
Of note, human disease-based enrichment results found that transcripts for genes involved in cardiovascular diseases, including dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM) (Fig. 3B and C, Supplemental Tables 3 and 4), were significantly affected in cftrscu101 mutants. Then, we selected 10 representative genes with significant expression change for further identification, quantitative PCR results showed that all detected genes demonstrated similar expression pattern in embryos (Fig. 3D). In conclusion, gene expression analysis comparing WT and cftr mutant embryos indicated that cftr deficiency induces significant changes in transcript levels, resulting in the altered expression of many mRNAs critical to cardiovascular diseases.
2.4 CFTR I556V mutation in human dilated cardiomyopathy
Due to the severe cardiomyopathy and alteration of cardiovascular disease-related mRNAs resulting from the loss of cftr in zebrafish, as well as to provide more convincing evidence for the role of CFTR in myocardial development, we considered that whether CFTR was associated with human heart trouble. So, the potential involvement of CFTR in a Chinese population of dilated cardiomyopathy was evaluated.
Fortunately, we enrolled 1 propositus (see Supplemental Table 5 for clinical information and Supplemental Fig. 2 for sequencing data) with family history of diagnosed idiopathic dilated cardiomyopathy for CFTR mutation screening by whole-exome sequencing. Interestingly, we identified a heterozygous missense mutation (A→G at position 1666, resulting in the amino acid change I556V) in this propositus. Furthermore, this CFTR I556V mutation was not identified in over gender/age/geography/ethnicity matched 600 healthy controls, demonstrating a highly significant association between the CFTR mutation and dilated cardiomyopathy through exact Fisher’s test analysis (P < 0.001).
2.5 CFTR I556V mutation fails to rescue the cardiac dysplasia in zebrafish cftr mutant model
The CFTR I556V mutation was initially reported in a French male who presented with asthma-like bronchopathy and chronic diarrhea [24, 25]. I556V and a second mutation, F508del, are the most common types of CFTR mutations in Chinese patients with CF and congenital bilateral absence of the vas deferens (CBAVD) [26, 27]. Importantly, compared with WT, open probability (Po) of Cl− channel function was significantly reduced in the presence of I556V (by 34%) [28]. However, the association of the CFTR I556V mutation with cardiac development and dilated cardiomyopathy has not yet been investigated.
To investigate the function of CFTR I556V mutation in the heart, we sought to restore the cardiac dysplasia in cftr mutants by by injecting human CFTR I556V mRNAs into mutant embryos. CFTR wild-type or channel-defective human CFTR (∆F508 and G551D) mRNAs were injected in parallel for comparison. The CFTR ΔF508 mutation results in protein misfolding and subsequent protein retention in the endoplasmic reticulum (ER), with a negligible amount of mutant CFTR reaching the plasma membrane, leading to impaired ion channel function [29], while CFTR G551D is a well-known mutation causing gating defect in CFTR channel function [30, 31].
As with the channel-defective ΔF508 and G551D mutants, CFTR I556V failed to restore the cardiac dysplasia in cftr mutants, as evidenced by altered nkx2.5 expression at 13 hpf (Fig. 4A-C), myl7 expression at 48 hpf (Fig. 5A-C), heart rate (Fig. 6A) and systolic function at 60 hpf (Fig. 6B). In contrast, WT CFTR mRNA markedly restored the cardiogenesis in offspring embryos derived from the mutant line. Interestingly, a ΔPDZ mutant, which lacks the protein interaction module but has intact channel function, also rescued the cardiogenic defects. Taken together, these results, along with those demonstrating impaired cardiogenesis in the presence of CFTRinh-172, indicate that cftr ion channel activity is essential to cardiac development.