Decreased activity of RCAN1.4 is a potential risk factor for congenital heart disease in a Han Chinese population

Protein & Cell, Mar 2018

Liangping Cheng, Peiqiang Li, He Wang, Xueyan Yang, Huiming Zhou, Wufan Tao, Jie Tian, Hongyan Wang

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Decreased activity of RCAN1.4 is a potential risk factor for congenital heart disease in a Han Chinese population

Decreased activity of RCAN1.4 is a potential risk factor for congenital heart disease in a Han Chinese population Dear Editor - Congenital heart disease (CHD) is a major birth defect worldwide. However, the aetiology of CHD remains unclear and the detailed molecular mechanisms underlying CHD pathogenesis have not been fully understood. The human regulator of calcineurin (RCAN1) gene is highly expressed in human hearts and brains (Fuentes et al., 1995) and composed of seven exons that are alternatively spliced and/or transcribed by differential promoters to produce different isoforms (Fig. 1A) (Yang et al., 2000) . RCAN1.1 (NM_004414) and RCAN1.4 (NM_203418) are the two major isoforms differentially expressed in many tissues and cells. RCAN1.1 is constitutively expressed while the expression of RCAN1.4 is induced by diverse stimuli (Yang et al., 2000; Lange et al., 2004) . RCAN1.4 negatively regulates NFAT-mediated transcription (Lange and Yutzey, 2006; Qin et al., 2006) , which is an important regulator for early heart development (de la Pompa et al., 1998) . Although SNPs rs149048873, rs193289374 and rs143081213 in the RCAN1.1 promoter were not associated with sporadic CHD (Guo et al., 2015; Li et al., 2015) , the function of RCAN1.4 in CHD remains elusive. To investigate if the RCAN1.4 plays a role in heart development, we evaluated the expression of the RCAN1.1 and RCAN1.4 proteins in E10.5 mouse embryonic hearts and human hearts. The results revealed that both the RCAN1.1 and RCAN1.4 were expressed in mouse embryonic hearts and human hearts (Fig. 1B). In addition, the RCAN1.4 protein was easily detected in embryonic mouse hearts at E10.5 but not in embryonic mouse heads at E10.5, while RCAN1.1 protein levels were similar between embryonic mouse hearts and heads (Figs. 1B and S1A). These results demonstrate that not only RCAN1.4 is expressed in hearts, but also its expression is differentially regulated, at least during mouse heart development, suggesting that the expression level of RCAN1.4 might contribute to human CHD. To evaluate the association of RCAN1 gene with human CHD, we performed an association study of variants in RCAN1 with CHD. The −2 kb of 5′-upstream regulatory regions, coding regions and 3′UTRs for RCAN1.1 and RCAN1.4 transcripts were first sequenced by target-capture sequencing in 412 blood samples of sporadic CHD patients and 213 matched controls from Shandong populations (Table S1). Eight SNPs with a minor allele frequency (MAF) > 0.05 were identified within the −2 kb of 5′ regulatory region of RCAN1.1 and RCAN1.4 transcripts. Three SNPs (rs765610, rs1248192 and rs36012405) and the other five (rs10550296, rs8133540, rs2243890, rs4816491 and rs2300385) were positioned in 5′ regulatory regions of RCAN1.1 and RCAN1.4 transcripts, respectively (Fig. 1A). We also identified one rare mutation (rs149533318, RCAN1.4E25K) in the coding region of RCAN1.4. Among them, only the genotype distribution of SNP rs2243890 A>G (HGVS Name: NM_203418:c.-1712 A>G) located in RCAN1.4 promoter was normally different between the CHDs and control (Table S2). We observed the higher frequencies of G allele in CHD (MAF = 0.09) than control (MAF = 0.05) (OR = 1.98, 95% CI = 1.20–3.27, P = 0.007 in additive model; OR = 2.11, 95% CI = 1.26–3.53, P = 0.004 in dominant model) (Table S2). To further confirm the validity of the association of rs2243890 A>G with CHD, we not only verified the status of rs2243890 A>G by re-sequencing a 0.5-kb genomic sequence around the SNP in the previously analysed 412 sporadic CHD patients and 213 matched controls, but also examined the same sequence in 458 newly collected sporadic CHD patients and 1,107 healthy controls from Shanghai by Sanger sequencing. Consistent results were obtained for this SNP (OR = 1.57, 95% CI = 1.19–2.07, P = 0.002) (Table 1). Subsequently, we combined the data from both cohorts for combined analysis and the result showed that the distribution of the minor G allele of rs2243890 was significantly different between the CHD and control groups. The CHD risk was significantly increased by 1.66-fold due to the G allele (OR = 1.66, 95% CI= 1.30–2.11, P = 4.68E−05) (Table 1). The frequencies of all genotypes among the control subjects were in accordance with the Hardy-Weinberg expectation (P > 0.05). The statistical power of the study is 0.924 at the 0.05 significance level, which is calculated with the power and sample size RCAN1.1 Promoter rs36012405 rs12481927 rs765610 A T C T T C A G C A/G A A A G T G A A A C C A H 220 *** 286 360 ** *** 368 ns 390 320 Normal Abnormal l a Norm Abnorm l a E1 E2 E3 E4 E5 E6 E7 Rcan1.1 Rcan1.4 β-Actin E 0 -1712A/A -1712A/G (N = 28) (N = 22) 0 -1885C/C -1885G/G (N = 39) (N = 11) Labeled A probe Labeled G probe L l/ o m μ 0 0 5 b calculation (version 3.1.2, wiki/Main/PowerSampleSize). A further stratified analysis of rs2243890 A>G was performed based on different subtypes of CHDs. The most significant differences in CHD vs. control were observed in an additive genetic model in septal defects, which had a 1.51-fold higher risk in sporadic CHD patients (OR = 1.51, 95% CI = 1.11–2.07, P = 0.01) and conotruncal defects, which had a 2.18-fold higher (OR = 2.18, 95% CI = 1.59– 2.99, P = 1.33E−06) risk in sporadic CHD patients (Table S3). There were no significant associations observed between rs2243890 A>G and other CHD subtypes, including right ventricular outflow tract obstruction (RVOTO) (P = 0.12), left ventricular outflow tract obstruction (LVOTO) (P = 0.34) and patent ductus arteriosus (PDA) (P = 0.41) (Table S3). Since the minor G allele of rs2243890 (at −1,712 bp) was located in the promoter region for RCAN1.4 transcript, we hypothesized that the G allele may contribute to the risk of CHD by influencing RCAN1.4 transcription. To test our hypothesis, we evaluated the RCAN1.4 mRNA level in human heart tissue samples from sporadic CHD patients by qRT-PCR. We found that RCAN1.4 mRNA level in human heart samples with heterozygous A/G alleles was almost threefold lower than that of samples with homozygous A/A alleles (Fig. 1C), while as a control, the difference in RCAN1.4 mRNA level between the heart samples carrying either heterozygous C/G alleles or homozygous G/G alleles of rs2300385 (at −1,885 bp) was not significant (Fig. 1C). These results suggest that the minor G allele of rs2243890 in RCAN1 may affect CHD by decreasing rather than enhancing the transcription of RCAN1.4. To further confirm the effect of rs2243890 A>G, we evaluated the RCAN1.4 promoter activity of the A and G alleles of rs2243890 using a dual-luciferase reporter assay in HEK293, H9C2 and C2C12 cells. The results showed that the transcriptional activity of the G allele was significantly lower than that of the A allele in all three cell lines (Fig. 1D), demonstrating that the G allele of rs2243890 attenuates RCAN1.4 transcription. Taken together, these results suggest that the minor G allele of rs2243890 in RCAN1 may contribute to CHD by inhibiting rather than enhancing RCAN1.4 transcription. As the minor G allele of rs2243890 decreased RCAN1.4 promoter activity (Fig. 1D), we questioned whether these changes might be due to its different binding affinity for some transcription factors. Using a major probe and a minor probe containing the A allele and G allele, respectively (Table S5), we performed an EMSA to evaluate the binding affinity of these alleles for unknown transcription factors in 293T cells. The results showed that the G allele displayed a higher binding affinity for some nuclear protein(s) than the A allele did (Fig. 1E). As the minor G allele decreased RCAN1.4 promoter activity compared with the major A allele, we inferred that the specific binding protein(s) probably function as transcriptional repressors. l l i t a AA/AG/GG. b OR (95% confidence interval (CI)) and P values were derived from logistic regression analysis in genetic models with adjustment for sex. c P for the HWE test in the controls. d Pooled from fixed effects model (Pheterogeneity > 0.05). Study Shandong samples Shanghai samples Combinedanalysis By Sanger sequencing we also confirmed the CHDspecific rare mutation RCAN1.4E25K identified from targetcapture sequencing (Fig. 1F). The MAF for RCAN1.4E25K in the ExAC database was as rare as 3.35E−05 (4/119506). To assess whether this RCAN1.4E25K mutant plays a role in CHD development, we performed an NFAT promoter-mediated luciferase reporter assay. Results showed that mutant RCAN1.4E25K significantly impaired suppression of NFATmediated transcription (Figs. 1G and S1B), suggesting that RCAN1.4E25K may be a loss-of-function mutant. The E25 of human RCAN1.4 protein is conserved from the zebrafish to human (Fig. S1C). About 70% of zebrafishes injected with antisense morpholino-modified oligonucleotides (MO) for zebrafish rcan1a-4 (a zebrafish homolog of human RCAN1.4) (Alghanem et al., 2017) displayed an enlarged pericardium (a CHD phenotype) and this rcan1a-4 MO-mediated CHD phenotype could be partially rescued by co-injected human RCAN1.4 mRNA in a dosage-dependent manner (Fig. 1H). However, human RCAN1.4E25K mRNA failed to rescue rcan1a-4 MO-mediated CHD phenotype, suggesting again that RCAN1.4E25K is a loss-of-function mutation in vivo. In this study, the SNP rs2243890 in RCAN1.4 promoter region was found to be significantly associated with CHD and further validated to be significantly associated with septal and conotruncal defects of CHD by analysing blood samples from 870 sporadic CHD patients and 1,320 normal controls from a Han Chinese population. Our results, therefore, adds one more piece of evidence supporting the conclusion that genetic associations with CHD have a high phenotypic specificity (Soemedi et al., 2012; Cordell et al., 2013) and also suggest a population specificity for genetic association of SNP rs2243890 with CHD because reported GWAS studies did not find the association of SNP rs2243890 with CHD in Caucasian (Cordell et al., 2013). Consistent with the results from two previous genetic association studies (Guo et al., 2015; Li et al., 2015) , we confirmed that three other SNPs (rs765610, rs12481297 and rs36012405) in the RCAN1.1 promoter region were also not associated with CHD. Our Western blot analyses further showed that the RCAN1.4 protein was expressed in both mouse and human embryonic hearts and the RCAN1.4, but not RCAN1.1, protein level was differentially up-regulated in developing mouse embryonic hearts (Fig. 1B). These results support the notion that RCAN1.4, but not RCAN1.1, may contribute to CHD pathogenesis in the Han Chinese population. RCAN1, also known as Down syndrome candidate region 1 (DSCR1), is located in syntenic regions of human chromosome 21. Approximately 40%–60% of DS patients are accompanied by CHD (Stoll et al., 1998; Vis et al., 2009; Elmagrpy et al., 2011) . RCAN1 triplication was thought to be a potential cause for CHD in DS (Lange et al., 2004) . However, Lyle et al. identified eight new DS cases of partial monosomy 21. One patient carried only a 1.48 Mb deletion containing eight genes, including RCAN1, and displayed a relatively severe DS phenotype, including cardiac anomaly (Lyle et al., 2009) . Now, our study revealed that rs2243890 A>G in the RCAN1.4 promoter increases the risk of CHD by decreasing RCAN1.4 transcription. We not only demonstrated that the minor G allele of rs2243890 significantly decreased RCAN1.4 transcription but also showed that the rare mutation RCAN1.4E25K identified in a VSD patient almost completely abolished RCAN1.4 activity in both in vivo zebrafish model as well as in an in vitro luciferase reporter assay. Altogether, these studies suggest that the expression or activity of RCAN1, rather than the RCAN1 copy number, may contribute to DS-associated CHD. Our study may provide a new understanding of the role of the RCAN1 gene in human CHD. FOOTNOTES This work was supported by the grants from National Key Research and Development Program of China (2016YFC1000500) to Hongyan Wang and Wufan Tao, the National Natural Science Foundation of China (Grant Nos. 31671330 and 31471380) to Wufan Tao and (81472050) to Xueyan Yang, the National Natural Science Foundation of China (Grant Nos. 81430005 and 31521003) and the Commission for Science and Technology of Shanghai Municipality (13JC1407600) to Hongyan Wang. Liangping Cheng, Peiqiang Li, He Wang, Xueyan Yang, Huiming Zhou, Wufan Tao, Jie Tian and Hongyan Wang declare that they have no conflict of interest. All procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation (institutional and national) and with the Helsinki Declaration of 1975, as revised in 2000 (5). Informed consent was obtained from all patients for being included in the study. Additional informed consent was obtained from all patients for which identifying information is included in this article. All institutional and national guidelines for the care and use of laboratory animals were followed. Liangping Cheng1,2 , Peiqiang Li3,4, He Wang4, Xueyan Yang4, Huiming Zhou1, Wufan Tao2,5&, Jie Tian1&, Hongyan Wang2,4& 1 Key Laboratory of Pediatrics in Chongqing, Chongqing International Science and Technology Cooperation Center for Child Development and Disorders, Department of Cardiovascular Medicine, Children’s Hospital of Chongqing Medical University, Chongqing 400014, China 2 Obstetrics Gynecology Hospital, The Institute of Reproduction and Developmental Biology, Fudan University, Shanghai 200011, China 3 Institute of Genetics, School of Basic Medical Sciences, Lanzhou University, Lanzhou 730000, China 4 The State Key Laboratory of Genetic Engineering, MOE Key Laboratory of Contemporary Anthropology, Collaborative Innovation Center of Genetics and Development, School of Life Science, Fudan University, Shanghai 200433, China 5 Institute of Developmental Biology and Molecular Medicine, Fudan University, Shanghai 200433, China & Correspondence: (W. Tao), (J. Tian), (H. Wang) OPEN ACCESS This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. l l i t Electronic supplementary material The online version of this article ( contains supplementary material, which is available to authorized users. Alghanem AF , Wilkinson EL , Emmett MS , Aljasir MA , Holmes K , Rothermel BA , Simms VA , Heath VL , Cross MJ ( 2017 ) RCAN1.4 regulates VEGFR-2 internalisation, cell polarity and migration in human microvascular endothelial cells . Angiogenesis 20 : 341 - 358 Cordell HJ , Bentham J , Topf A , Zelenika D , Heath S , Mamasoula C , Cosgrove C , Blue G , Granados-Riveron J , Setchfield K et al ( 2013 ) Genome-wide association study of multiple congenital heart disease phenotypes identifies a susceptibility locus for atrial septal defect at chromosome 4p16 . Nat Genet 45 : 822 - 824 de la Pompa JL , Timmerman LA , Takimoto H , Yoshida H , Elia AJ , Samper E , Potter J , Wakeham A , Marengere L , Langille BL et al ( 1998 ) Role of the NF-ATc transcription factor in morphogenesis of cardiac valves and septum . Nature 392 : 182 - 186 Elmagrpy Z , Rayani A , Shah A , Habas E , Aburawi EH ( 2011 ) Down syndrome and congenital heart disease: why the regional difference as observed in the Libyan experience? Cardiovasc J Afr 22 : 306 - 309 Fuentes JJ , Pritchard MA , Planas AM , Bosch A , Ferrer I , Estivill X ( 1995 ) A new human gene from the down syndrome critical region encodes a proline-rich protein highly expressed in fetal brain and heart . Hum Mol Genet 4 : 1935 - 1944 Guo RY , Li XF , Bai S , Guo J , Ding N , Li ZZ ( 2015 ) Association between DSCR1 variations and congenital heart disease susceptibility . Med Sci Monit 21 : 3536 - 3539 Lange AW , Molkentin JD , Yutzey KE ( 2004 ) DSCR1 gene expression is dependent on NFATc1 during cardiac valve formation and colocalizes with anomalous organ development in trisomy 16 mice . Dev Biol 266 : 346 - 360 Lange AW , Yutzey KE ( 2006 ) NFATc1 expression in the developing heart valves is responsive to the RANKL pathway and is required for endocardial expression of cathepsin K. Dev Biol 292 : 407 - 417 Li X , Wang G , An Y , Li H , Li Y , Wu C ( 2015 ) Association between sequence variations in RCAN1 promoter and the risk of sporadic congenital heart disease in a Chinese population . Pediatr Cardiol 36 : 1393 - 1399 Lyle R , Bena F , Gagos S , Gehrig C , Lopez G , Schinzel A , Lespinasse J , Bottani A , Dahoun S , Taine L et al ( 2009 ) Genotype-phenotype correlations in Down syndrome identified by array CGH in 30 cases of partial trisomy and partial monosomy chromosome 21 . Eur J Hum Genet 17 : 454 - 466 Qin L , Zhao D , Liu X , Nagy JA , Hoang MV , Brown LF , Dvorak HF , Zeng H ( 2006 ) Down syndrome candidate region 1 isoform 1 mediates angiogenesis through the calcineurin-NFAT pathway . Mol Cancer Res 4 : 811 - 820 Soemedi R , Topf A , Wilson IJ , Darlay R , Rahman T , Glen E , Hall D , Huang N , Bentham J , Bhattacharya S et al ( 2012 ) Phenotypespecific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls . Hum Mol Genet 21 : 1513 - 1520 Stoll C , Alembik Y , Dott B , Roth MP ( 1998 ) Study of Down syndrome in 238,942 consecutive births . Ann Genet 41 : 44 - 51 Vis JC , Duffels MG , Winter MM , Weijerman ME , Cobben JM , Huisman SA , Mulder BJ ( 2009 ) Down syndrome: a cardiovascular perspective . J Intellect Disabil Res 53 : 419 - 425 Yang J , Rothermel B , Vega RB , Frey N , McKinsey TA , Olson EN , Bassel-Duby R , Williams RS ( 2000 ) Independent signals control expression of the calcineurin inhibitory proteins MCIP1 and MCIP2 in striated muscles . Circ Res 87 : E61 - E68

This is a preview of a remote PDF:

Liangping Cheng, Peiqiang Li, He Wang, Xueyan Yang, Huiming Zhou, Wufan Tao, Jie Tian, Hongyan Wang. Decreased activity of RCAN1.4 is a potential risk factor for congenital heart disease in a Han Chinese population, Protein & Cell, 2018, 1-6, DOI: 10.1007/s13238-018-0525-8