- Research article
- Open Access
Long-range regulatory interactions at the 4q25 atrial fibrillation risk locus involve PITX2c and ENPEP
BMC Biologyvolume 13, Article number: 26 (2015)
Recent genome-wide association studies have uncovered genomic loci that underlie an increased risk for atrial fibrillation, the major cardiac arrhythmia in humans. The most significant locus is located in a gene desert at 4q25, approximately 170 kilobases upstream of PITX2, which codes for a transcription factor involved in embryonic left-right asymmetry and cardiac development. However, how this genomic region functionally and structurally relates to PITX2 and atrial fibrillation is unknown.
To characterise its function, we tested genomic fragments from 4q25 for transcriptional activity in a mouse atrial cardiomyocyte cell line and in transgenic mouse embryos, identifying a non-tissue-specific potentiator regulatory element. Chromosome conformation capture revealed that this region physically interacts with the promoter of the cardiac specific isoform of Pitx2. Surprisingly, this regulatory region also interacts with the promoter of the next neighbouring gene, Enpep, which we show to be expressed in regions of the developing mouse heart essential for cardiac electrical activity.
Our data suggest that de-regulation of both PITX2 and ENPEP could contribute to an increased risk of atrial fibrillation in carriers of disease-associated variants, and show the challenges that we face in the functional analysis of genome-wide disease associations.
Genome-wide association studies (GWAS) have exponentially increased our knowledge of the genetic component of human disorders, revealing unsuspected loci that harbour variants linked to an increased risk of disease . However, the majority of GWAS signals fall in non-coding regions of the genome, which has made their functional analysis particularly challenging [2,3]. Even the identification of the genes targeted by disease-associated variants is not straightforward, as mere proximity can result in incorrect identification of the culprit gene .
Atrial fibrillation (AF) is the most common cardiac arrhythmia in humans  although its pathophysiologic basis is still not clearly understood, presenting a challenge for cardiovascular research and therapy. AF is defined as a supraventricular tachyarrhythmia characterised by uncoordinated atrial activation, and is frequently observed as a consequence of various systemic and cardiac disorders (syndromic AF) . However, in 10% to 20% of cases AF is not associated with other cardiovascular disease, and thus is dubbed ‘idiopathic’ or ‘lone’ AF that mostly occurs in patients under 60. The strong association of AF onset with risk factors, such as age, sex, ethnicity, hypertension and other heart diseases , originally suggested it being a non-genetic disorder . Nevertheless, in the last two decades several epidemiological studies pointed to a significant incidence of genetic factors . Furthermore, rare mutations in a dozen genes, mostly encoding ion channel subunits , are associated with AF as part of wider cardiac electrical syndromes.
At least ten loci have been linked to AF by GWAS in large cohorts of non-related patients of distinct ethnic backgrounds [10-12]. The most highly AF associated variants identified in all studies are located on chromosome 4q25 , 170 kilobases (kb) distal to PITX2 and within a 1.5 megabases (Mb) intergenic gene desert. PITX2 encodes an evolutionarily conserved homeodomain transcription factor that is involved in the establishment of left-right asymmetry and cardiovascular development in the vertebrate embryo. In mice and humans the PITX2 gene generates several isoforms. PITX2a and PITX2b are alternative splicing variants produced from a common promoter, whereas PITX2c is the product of an alternative promoter and is the main isoform expressed in the heart .
Mouse Pitx2c is first expressed in the left lateral plate mesoderm of the early embryo, as part of the network regulating the establishment of left-right asymmetry during development , and is then expressed in the left side of the heart at early stages . At later stages, expression follows a dynamic pattern, being present in the left atrium, the myocardium sleeves of the pulmonary veins, the atrio-ventricular cushions or the base of the ventricles [16,17]. Germline homozygous deletion of Pitx2 results in embryonic lethality and numerous cardiac malformations, such as right atrial isomerism and outflow tract defects, varying from double outlet right ventricle or transposition of the great arteries to persistent truncus arteriosus . More recent research showed that Pitx2c is expressed in adult mice and human hearts, predominantly in the left atrium (levels in the right atrium and the ventricles are 100-fold lower), and that its levels of expression decrease in atria of AF patients [13,19]. Furthermore, Pitx2c heterozygous or atrial specific deletion of Pitx2c display molecular and physiological hallmarks of human AF [13,19,20], which also is observed when Pitx2 is deleted in adult mice . Altogether, these data support the hypothesis that PITX2 could play a causal role in the pathogenesis of AF and that its function could be altered by genomic elements located in the vicinity of the single nucleotide polymorphisms (SNP) in 4q25 that correlate with of AF.
In this study, we explored whether the 4q25 region spanning the AF-associated variants identified by GWAS harbours putative regulatory elements that could be acting on neighbouring genes. By tissue culture, in vivo transgenics and analysis of chromatin structure, we have found that this region contains potentiator cis-regulatory elements that interact with the promoters of Pitx2c and, unexpectedly, Enpep, the next gene located downstream of Pitx2. Given the expression of Enpep in the sinoatrial node (SAN) and the co-expression of Pitx2 and Enpep in pro-arrhythmogenic regions of the embryonic heart, such as the sleeves of the pulmonary veins, our data suggest that de-regulation of these genes could underlie increased risk of AF.
Genomic analysis of the 4q25 AF-associated locus
To identify putative cis-regulatory elements located in the 4q25 region we analysed the evolutionary conservation  and distribution of histone modifications associated with active elements  (H3K4me1) in an 85 kb window containing the main SNPs that have been associated with an increased risk of AF by GWAS (Figure 1). This window is centred on the lead rs2200733 SNP  that has repeatedly been identified as the most significant variant associated with AF [10-12,24-26], and spans a region that includes other distal SNPs (rs2634073 and rs17570669) that lie in the proximity of sequences highly conserved between human and mouse (Figure 1A). The region including the majority of AF-associated SNPs in 4q25 is confined to a linkage disequilibrium (LD) block , separated from that containing the coding exons of the gene and from adjacent LD blocks in the 1.5 Mb gene desert located centromeric to PITX2 (Figure 1A). It is also noteworthy that the selected SNPs and the PITX2 gene are contained together in a single topologically associated domain (TAD; Additional file 1A), as defined by Hi-C in three different human cell lines .
We selected seven genomic fragments (AF1-7) for further analysis (Figure 1B). Fragments AF1 and AF2 lie in close proximity to rs2634073 and show high evolutionary conservation. These fragments are included in a region (hs930) tested as part of a large scale screen for tissue specific human enhancers by transgenesis in the mouse embryo  and drive reporter expression in the nervous system and limbs but not in the heart. This region was also tested in transgenic zebrafish, driving expression in similar patterns but again not in the developing heart . AF3 to AF5 are a set of overlapping fragments that include the lead rs2200733 variant and other highly associated SNPs, in a region with high conservation among placental mammals and H3K4me1 marks of active regulatory elements. Finally, AF6 and AF7 map to a region conserved in vertebrates including rs17570669, which has been associated with AF but is independent of rs2200733 .
The 4q25 AF-associated locus contains active regulatory elements
We tested the regulatory activity of these fragments, corresponding to the none-risk haplotype at rs2200733 from a commercial source of human DNA, by linking them to a human minimal beta-globin promoter  and the reporter gene coding for monomeric red fluorescent protein (mRFP) . These constructs were transfected into cultured HL-1 mouse atrial cardiomyocytes , a tissue culture model of the physiological conditions in which PITX2 is active . As a positive control, we used a previously described enhancer from the Nppa gene (encoding atrial natriuretic factor; ANF) that recapitulates its endogenous expression in transgenic mice ; for negative controls we used the reporter constructs lacking any genomic fragment, and containing the pluripotent-specific Oct4 distal enhancer . Of the seven fragments, only AF3, which contains rs2200733 (Figure 1B), showed significant activity compared with negative controls (Figure 2A). While this could seem surprising given that other fragments show a more robust signal for histone modifications associated to regulatory elements (for example AF4, Figure 1B), enhancer prediction based on histone marks is only accurate in a fraction of cases .
Due to the large size of AF3 (over 7 kb), we further analysed its function by generating six overlapping fragments (AF3.1 to AF3.6) of 1 to 2 kb each (Additional file 1B). When tested in HL-1 cells, only AF3.5 and AF3.6 showed activity (Figure 2B). It is noteworthy that AF3.3, which contains rs2200733, is not active in this assay. AF3.5 and AF3.6 overlap in 80 base pairs (bp), so we then tested the activity of this minimal fragment (AF3.5∩3.6) in HL-1 cells as well as versions of AF3.5 and AF3.6 where the overlap was removed (AF3.5Δ, AF3.6Δ). AF3.5∩3.6 showed strong activity in HL-1 cells, and while AF3.5Δ was not active, AF3.6Δ retained activity although at a reduced level (Figure 2B). It is interesting to note that this 80 bp minimal fragment is highly conserved between human and mouse, and that its sequence corresponds to a short interspersed nuclear element of the MIR3 family. It has been shown that repeat sequences can act as enhancers in experimental assays , although their putative function in vivo is still under debate . Therefore, we can conclude that the regulatory activity of AF3 in this assay is located in AF3.6, and that additional activity may be present in the overlapping fragment of AF3.5 and AF3.6.
Next, we assayed the activity of fragments AF3-5 and AF7 in transgenic mouse embryos, using the lacZ gene as a reporter. Again, only AF3 showed enhancer activity (Figure 2C, Additional file 2), confirming the results of the tissue culture assays. AF3 drives reporter expression in facial mesenchyme, limb muscles, and the left gonad, some of which are sites of expression of endogenous Pitx2 [14,38]. Contrary to expectations, AF3 did not drive expression in the developing heart. We reasoned that regulatory elements from PITX2 underlying the association with AF might not be active during development, and instead drive cardiac-specific expression of PITX2 in the adult. We therefore generated transgenic mice and examined reporter expression in the heart at postnatal day 3. Again, we found no expression in cardiac tissues (Additional file 2).
When tested in mouse transgenic assays, both AF3.5 and AF3.6 showed activity (Figure 2D to F). Sites of expression include facial mesenchyme, limb muscles, somites, or pericardium, but we did not observe a reproducible pattern driven by these fragments. We can rule out the possibility that this heterogeneity is due to non-specific reporter activation as a consequence of integration site of the transgene, because of the very low percentage of lacZ positive embryos (out of the total number of transgenics as assessed by genotyping) obtained for genomic fragments tested showing no activity (0% to 5%) as compared to those that do (20% to 45%; Additional file 2).We also tested activity in transgenic embryos of the minimal AF3.5∩3.6 fragment, finding that it was not active (1 weak lacZ+ embryo out of 17 transgenics; Additional file 2).
The 4q25 regulatory elements show non-specific potentiator activity
The above results suggest that these 4q25 elements, while they have regulatory potential, do not confer tissue specificity. To test this hypothesis further, we transfected fragments AF3, AF3.5 and AF3.6 into two cell types unrelated to the cardiac lineage: the mouse teratocarcinoma-derived pluripotent cell line P19 and human embryonic kidney (HEK) cells. We found that all three fragments were active in both cell types, closely matching the degree of activation in HL-1 cells (Additional file 3). As expected, the Oct4-DE was active in P19 but not in HEK cells; in contrast, the Nppa enhancer was not active in P19 but showed activity in HEK cells, as expected given the endogenous expression of NPPA in human kidney . Overall, our results suggest that the regulatory elements detected in 4q25 do not act as cell type-specific enhancers, but rather as accessory elements that can potentiate the activity of tissue-specific enhancers located elsewhere in the locus.
To further prove the putative potentiator activity of 4q25 elements, we assessed the effect of AF3.6 on the activity of a previously identified intronic enhancer from Pitx2, which drives left-sided expression in the embryo . This asymmetric enhancer (ASE) is evolutionarily conserved in sequence and function, but it is noteworthy that the ASE from human PITX2 only drives weak expression in the mice heart compared with its mouse homologue . We generated a chimeric construct containing both human AF3.6 and ASE and compared its activity to that of ASE alone in transgenic mouse embryos at 10.5 (Figure 2G-L). We first observed that when using the chimeric ASE + AF3.6 construct, the variability associated with AF3.6 alone is lost, and all embryos show the characteristic left-sided expression described for the ASE (Figure 2G, H). Importantly, we found that there is no additive effect of both genomic fragments as that observed when placing together different enhancers in the same transgenic construct . In fact, we observed that adding AF3.6 to the ASE apparently increased the levels of reporter expression in the cardiac region (Figure 2G, H), as seen in sections where domains of reporter activity in the right ventricle and in the atrio-ventricular canal are expanded in ASE + AF3.6 embryos compared to ASE (Figure 2I-L). When we examined in detail reporter expression in the developing heart for all transgenic embryos, we found that AF3.6 increases the number of embryos expressing lacZ in the left atrium (two out of five for ASE, as compared to seven out of ten for ASE + AF3.6; Additional file 4). The results of these assays further suggest that 4q25 elements have an accessory role in defining PITX2 expression acting in conjunction with other regulatory elements.
The three-dimensional architecture of the Pitx2 locus identifies promoter-specific long range interactions
Although the above evidence shows that 4q25 includes regulatory elements, there is no direct evidence that this genomic region acts on PITX2 or that, if it does, it shows any specificity regarding the cardiac and non-cardiac isoforms produced from two alternative promoters. To answer these questions, we analysed the three-dimensional organisation of the locus by chromosome conformation capture (3C) [43,44], to address if these regions physically contact the promoters of PITX2. Because we aimed to perform the assays in the physiological context of the heart, we selected the mouse genome region syntenic to human 4q25 and analysed chromatin from the atria and ventricles of adult mouse hearts. Based on the relative order and position of conserved sequence blocks between human and mouse, we could unambiguously map mouse fragments (af1-af7) equivalent to the human AF1-AF7 fragment series tested above (Figure 3A).
We probed the interaction of HindIII restriction fragments containing the Pitx2a,b or Pitx2c promoter with ten fragments spanning over 100 kb of the distal region on mouse chromosome 3 syntenic to the AF-associated 4q25 locus in humans (Figure 3A). Using atria and ventricles from adult mice we observed a clear pattern of long range interactions, with regions interacting specifically with the Pitx2a,b promoter (fragment 2R; Figure 3B), the Pitx2c promoter (fragment 6F; Figure 3C), or with both (fragment 8F; Figure 3B, C). The latter result prompted us to ask if fragment 8F had regulatory activity in HL-1 cells, since it is not included in any of the previously tested fragments. This was not the case, suggesting that this genomic region has other architectural roles in configuring the regulatory landscape of Pitx2. To further examine the specificity of chromatin interaction between the AF-associated region and Pitx2c, we checked the interaction of a fragment containing the promoter of Enpep, the next neighbouring gene distal to Pitx2 in both mouse and humans (Figure 3A). To our surprise, we found a robust interaction between fragment 6F and Enpep (Figure 3D), suggesting that the 4q25 regulatory landscape is partially shared between PITX2 and ENPEP. The specificity of the interactions of fragments 2R, 6F and 8F with Pitx2ab, Pitx2c and Enpep was tested by using a series of control primers located upstream and downstream of the promoters, which showed no interactions (Additional file 5).
We next asked if the interactions we observed showed regional differences in atria, as AF constitutes a disorder of the left atrium and this is the region where PITX2 is prominently expressed . While fragment 2R interacted specifically with the Pitx2ab promoter in both left and right atrium, fragment 6F interacts only in left atrium with Pitx2c and Enpep. On the other hand, fragment 8F shows interaction with Pitx2c in both atria but only in the right atrium with Enpep (Figure 4A). Given the fact that the promoters of Pitx2c and Enpep share many of the interactions tested, we examined if they were physically associated and if this was region-specific. We found a robust and specific promoter-promoter interaction in both atria and in ventricles (Figure 4B), therefore independent of transcription and in line with recent observations on the role of pre-existing promoter-promoter interactions for structuring the genome .
It should be noted that fragment 6F contains the region conserved with human fragment AF3, thereby suggesting that the region with regulatory activity and that contains the lead SNP associated with AF (rs2200733) interacts in a specific manner with the promoters of the cardiac-specific isoform of Pitx2 and the neighbouring gene, Enpep. The 3C analysis of the mouse Pitx2/Enpep locus thus revealed an unexpected complexity of specific and shared chromatin interactions between the regions containing the potentiator elements and the different promoters studied that could be related to their function.
Enpep is expressed in arrhythmogenic sites in the embryonic heart
ENPEP encodes aminopeptidase A, which cleaves angiotensin II to produce angiotensin III as part of the renin-angiotensin system . Therefore ENPEP is involved in the control of blood pressure, and accordingly it is expressed in the renal system and endothelial cells, and knockout mice for Enpep develop hypertension . However, at present there is no report for expression or a role of Enpep in the heart. In light of our results, we examined the expression of Enpep in the E14.5 mouse embryos by in situ hybridization on tissue sections (Figure 5). Enpep is strongly expressed in the endothelial lining of the lungs, but also in a specific and restricted pattern in the developing heart (Figure 5A, D). We compared Enpep expression with that of Pitx2 (Figure 5B, E) and Hcn4 (Figure 5C, F), which encodes a voltage-gated ion channel and at this stage is a marker of most of the cardiac conduction system [48,49]. This analysis showed that Enpep is co-expressed with Pitx2 in the pulmonary veins but not in the myocardium of the left atrium (Figure 5A, B, D, E), and is co-expressed with Hcn4 in the left and right superior venae cavae and in the SAN (Figure 5A, C, D, F). Enpep is thus expressed in the embryonic mouse heart in key components of the cardiac conduction system such as the SAN. Moreover, Enpep is also expressed at the base of the pulmonary veins and the junction of the caval veins, regions prone to initiate ectopic electrical beats, which lead in many cases to the onset of AF .
The advent of GWAS has radically changed our perspective on the genetic analysis of common diseases in humans. On the one hand, a plethora of novel loci linked to increased disease risk have been uncovered, which await further analysis before possible translation to the clinic . On the other hand, the vast majority of risk variants are located in non-coding genomic sequences, pointing to a fundamental role for variation in cis-regulatory elements as the basis of common diseases [50,51]. Understanding the role and function of these genomic elements will be fundamental to making the most of the discoveries of GWAS.
The genomic analysis of AF is a prime example in this regard. All GWAS carried out to date have shown that the major loci for AF lie in an intergenic gene desert in 4q25, located distal to the developmental regulator PITX2 [10-12]. Despite its early role in establishing the left-right patterning of the heart and its prominent expression in the left atrium, no evidence suggested a role for PITX2 in the pathophysiology of AF . Subsequent analysis of loss-of-function mouse models of Pitx2 have confirmed that it plays a pivotal role in regulating different atrial phenotypes by distinguishing electrical from working myocardium in the right and left atria respectively [13,19,20]. However, no reports to date have provided evidence as to how distal variants in 4q25 act on PITX2. Even more surprisingly, a recent report showed that 4q25 variants do not correlate with PITX2 expression in atrial tissue from human patients . This evidences that variants identified by GWAS may have extremely subtle effects, which fall below the threshold of detection of current analytical tools and approaches.
Our analysis of the regulatory structure of the 4q25 locus shows that genomic sequences in close proximity to AF-linked variants can act as transcriptional regulatory elements both in tissue culture and in mouse embryos. Rather unexpectedly, and contrasting with other cases of GWAS-related enhancers in the cardiovascular system [53,54], these elements are not specific to cardiac cell types, either in culture or in vivo. The 4q25 elements show equal activity when transfected into cell types of different origin. Furthermore, in transgenic mouse embryos these elements drive highly variable patterns of reporter expression. These results suggest that the cis-regulatory elements in 4q25 do not act as classical tissue-specific enhancers, but as potentiator elements that would act in cooperation with elements located elsewhere in the locus that dictate tissue restricted expression. In the case of mouse Pitx2, an intronic enhancer (ASE) has been described that drives left-side specific expression in the early embryo and later in the heart, liver and other organs . It is conceivable that precise control of spatial and quantitative expression of PITX2 requires interplay of the ASE and the 4q25 potentiator. In fact, when placed together and tested by transgenics, the variability of 4q25 elements is lost and it can modulate the activity of ASE. It is possible that this potentiator could also modulate the activity of other yet to be identified regulatory elements from the locus.
By analysing the physical interaction between promoters and intergenic sequences of the mouse region syntenic to 4q25, we have found a further degree of complexity in the chromatin structure of the region. First, there is a clear specificity in the interaction of distal elements with the alternative promoters of the different Pitx2 isoforms, despite their being separated by less than 10 kb. In this regard, it is noteworthy that the region containing the potentiator activity we have described interacts specifically with the promoter of the cardiac-specific Pitx2c isoform. Furthermore, this interaction occurs specifically in the left atrium. These results provide additional support to the specific role of the region identified by GWAS in regulating PITX2 in the heart.
More surprising was the fact that this same region physically interacts with the promoter of Enpep, the neighbouring gene located distal to Pitx2. This opens the possibility that Enpep could also be a transcriptional target of the identified cis-regulatory elements. ENPEP, as part of the renin-angiotensin system, has been shown to control blood pressure, and hypertension is a known risk factor for AF . However, 4q25 variants are associated with lone AF, with no co-occurrence of hypertension , and independent variants located in the proximity of ENPEP but not in the 4q25 AF loci are associated with changes in blood pressure . Furthermore, there is no reported correlation between expression of ENPEP and 4q25 variants in the blood or adipose tissue . We can therefore conclude that the possible regulation of ENPEP by the 4q25 potentiator elements would be unrelated to its known role in the control of blood pressure. Our re-evaluation of Enpep expression in the developing mouse heart by in situ hybridization reveals co-expression with Pitx2 in the pulmonary veins, a region with pro-arrhythmogenic potential , and in the SAN of the right atria, a key component of the cardiac conduction system where the electrical impulse is generated. Our preliminary observations suggest that incorrect regulation of ENPEP in these locations could be linked to AF. The precise role of ENPEP in the heart remains to be identified and could offer novel insight into the pathogenesis of AF.
We have shown that novel cis-regulatory elements are located in the region of 4q25 associated with an increased risk for AF. These elements establish complex long-distance interactions with the promoters of both Pitx2c and Enpep, and therefore could regulate the transcription of these genes. A potential limitation of our study is the fact that while we have used human genomic DNA for regulatory assays, the chromatin structure of the Pitx2/Enpep locus and the expression of Enpep in the heart was carried out in mouse. However, the sequence conservation in the regions studied, as well as conserved synteny of the locus and of gene functions strongly suggests that regulatory mechanisms will also be conserved between human and mouse. Overall, our results suggest that de-regulation of either one or both PITX2 and ENPEP might have a causal role in the development of AF. Future work will be needed to identify the causal variants and the upstream regulatory factors that act through the potentiator elements described here.
Our study also highlights the challenges we face in the functional analysis of genetic variation identified by GWAS. Our understanding of the nature and function of non-coding genomic elements is still incomplete, despite the wealth of genome-wide data available through ENCODE and similar projects [58,59]. We are greatly limited by the breadth and specificity of available assays to interrogate the function of a DNA fragment. We can hypothesise that only a fraction of GWAS hits will represent classical tissue-specific enhancers, whose characterisation is feasible with current tools. Many cases will affect other regulatory elements with not such a clear-cut and easily identifiable role in gene transcription, such as potentiators or modulators (as we have identified here), but also silencers, insulators or stabilisers. Novel tools and assays will need to be devised to fully understand the regulatory variation underlying common human disease.
Commercial Clontech (Mountain View, California, USA) human DNA was used for PCR amplification of all the tested genome fragments from chromosome 4q25 (for primers used see Additional file 6). We used the pGem-T Easy Promega (Madison, Wisconsin, USA) vector for initial cloning of the PCR products, followed by digestion with NotI New England BioLabs (Ipswich, Massachusetts, USA) and subsequent cloning in enhancer-detection vectors containing the human minimal beta-globin promoter and either monomeric red fluorescence protein (pβRFP) or lacZ (p1230) reporter genes.
Cell culture and transfections
Mouse HL-1 atrial cardiomyocytes were cultured in Claycomb medium Sigma (St. Louis, Missouri, USA) supplemented with 10% (v/v) inactive (56°C, 30 minutes) fetal bovine serum (FBS) (Sigma), 4 mmol/L L-glutamine (Sigma), 100 μmol/L norepinephrine (Sigma) and 100 U/mL penicillin-streptomycin (Sigma). All seeding supports were previously coated for 24 hours with a solution of gelatin (0.02% w/v, Sigma) and fibronectin (25 μg/mL, Sigma). Mouse P19 embryonic teratocarcinoma cells (a kind gift from Christine Mummery, Leiden University Medical Center, The Netherlands) were cultivated in α-minimal essential medium (α-MEM, Gibco (Grand Island, New York, USA)) containing 10% FBS, 100 U/mL penicillin-streptomycin and 4 mmol/L L-glutamine. HEK293T human embryonic kidney cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Sigma) supplemented with 10% FBS, 4 mmol/L L-glutamine and 100 U/mL penicillin-streptomycin.
One day before transfections, cells were counted and plated at a density of 5 × 105 cells per p12 well (HL-1 cells) or p6 (P19 and HEK293T) with complete growth medium and no antibiotics. Cells were co-transfected with 2 μg of pβRFP vector containing the appropriate 4q25 fragment and 1 μg of pCAGGS-GFP (a kind gift from Joaquín Rodríguez-León, University of Extremadura, Badajoz, Spain) as an internal transfection efficiency control; co-transfections were performed with 6 μL of Lipofectamine 2000 Invitrogen (Waltham, Massachusetts, USA). Cells were transferred to complete medium with antibiotics after five hours. The empty vector pβRFP was used as a negative control.
Forty-eight hours after transfection, cultures were photographed (Zeiss) and fluorescent cells automatically counted (ImageJ) in twelve independent random fields per well (for the transfections of AF1 to AF7 in HL-1 cells; Figure 2A), or were measured by fluorescence activated cell sorting (FACS) (LSRFortessa (BD Biosciences; Franklin Lakes, New Jersey, USA) Flow Cytometer) in all other transfections. Three independent experiments with three technical replicates each were quantified in all cases. Relative regulatory enhancer activity was then calculated as the ratio of red cells (RFP+) to total green (GFP+) control cells, expressed as mean ± standard error of the mean (SEM) and statistically analysed by unpaired Student’s t-test (Prism5), with the significance threshold set at P <0.05.
Transient transgenic mice
p1230-derived constructs were digested with SacII and SalI (New England BioLabs) to remove the plasmid backbone, and the fragment was purified using the Qiagen gel extraction kit. DNA fragments were diluted in microinjection buffer (10 mmol/L Tris–HCl, pH7.4, 0.1 mmol/L ethylenediaminetetraacetic acid (EDTA)) at 5 to 7 ng/μL and injected into zygote pronuclei obtained from crosses of (C57BL/6xCBA/J)F1 mice. Injected zygotes were transferred to CD1 foster mothers, following standard procedures . At the desired stage, mice were euthanised and embryos dissected and stained for β-galactosidase activity . All embryos were genotyped for lacZ by PCR, using primers for Myogenin (Additional file 6) as an internal control for calculating transgenic efficiency and the percentage of embryos expressing lacZ (Additional file 2).
Animal studies were approved by the local ethics committee. All animal procedures conformed to EU Directive 2010/63EU and Recommendation 2007/526/EC regarding the protection of animals used for experimental and other scientific purposes, enforced in Spanish law under Real Decreto 1201/2005.
Chromosome conformation capture (3C) assays
The 3C protocol was performed essentially as described . Hearts from adult (C57BL/6xCBA/J)F1 female mice were dissected into atria and ventricles. After mincing with a scalpel, tissue was mechanically disrupted in 10 volumes of cold PBS, centrifuged at 3,000 g, and the cell supernatants cross-linked with 2% formaldehyde for eight minutes at room temperature. Nuclei were extracted with nuclear extraction buffer and the chromatin was digested with HindIII on a shaking platform at 37°C overnight. The cross-linked and digested chromatin products were ligated with T4 ligase (100 Weiss units) at 15°C for 12 hours in 7 mL 1 × ligation buffer. Sample quality was measured by semi-qPCR of the XPB/Eccr3 locus, as a control of non-tissue-specific chromatin conformation (see Additional file 6 for primer sequences). Only samples with more than 70% amplification efficiency were used as experimental templates. BAC clones (20 μg) containing XPB/Eccr3 (MRC Geneservice, clone 344-C18), Pitx2 (CHORI, clone RP24-215O15), Enpep (CHORI, clone RP24-172B1) and 3:G3 tested region (CHORI, clone RP23-356C23) were treated in parallel, to generate the control templates.
All primers used (Additional file 6) were designed in an approximately 300 kb region of mouse 3:G3 chromosome spanning the syntenic human 4q25 locus in which the GWAS-identified AF-related variants, PITX2 and ENPEP genes are located. Anchor primers were designed within the Pitx2ab, Pitx2c and Enpep promoter sequences (Figure 3). Three technical replicates of three independent experiments were performed for all sets of test-anchor primers for each tissue. Physical interactions among anchor and test primers, in the experimental and control templates, were measured by qPCR (SYBR® Green) and resulting frequencies were calculated and normalised using the XPB/Eccr3 locus as control [61,62]. Statistical analysis, assuming a normal distribution of data, was performed by one-way analysis of variance (ANOVA) test of Student-Newman-Keuls of the significance of differences among biological samples; the significance threshold was set at P <0.05. Error bars represent the SEM for the three biological replicates.
In situ hybridization
In situ hybridization was performed on sections of E14.5 embryos essentially as previously described . A pan-Pitx2 probe was kindly provided by José Luis de la Pompa (CNIC, Madrid, Spain). While this probe recognises all Pitx2 isoforms, only Pitx2c is expressed in the heart . Hcn4 and Enpep dsDNA were amplified by PCR from C57Bl/6 J DNA with primers containing T7 or SP6 RNA polymerase initiation sites (Additional file 6). Sense and anti-sense RNA probes were prepared by PCR using digoxigenin-labelled dsDNA as template Roche (Basel, Switzerland); sense probes were used as negative controls. Embryos used for different probes were processed in parallel in all assays.
chromosome conformation capture
green fluorescent protein
genome-wide association studies
human embryonic kidney
monomeric red fluorescent protein
standard error of the mean
single nucleotide polymorphism
topologically associated domain
Manolio TA. Genomewide association studies and assessment of the risk of disease. N Engl J Med. 2010;363:166–76.
Manolio TA. Bringing genome-wide association findings into clinical use. Nat Rev Genet. 2013;14:549–58.
Paul DS, Soranzo N, Beck S. Functional interpretation of non-coding sequence variation: concepts and challenges. Bioessays. 2014;36:191–9.
Smemo S, Tena JJ, Kim KH, Gamazon ER, Sakabe NJ, Gomez-Marin C, et al. Obesity-associated variants within FTO form long-range functional connections with IRX3. Nature. 2014;507:371–5.
Kirchhof P, Curtis AB, Skanes AC, Gillis AM, Samuel Wann L, John CA. Atrial fibrillation guidelines across the Atlantic: a comparison of the current recommendations of the European Society of Cardiology/European Heart Rhythm Association/European Association of Cardiothoracic Surgeons, the American College of Cardiology Foundation/American Heart Association/Heart Rhythm Society, and the Canadian Cardiovascular Society. Eur Heart J. 2013;34:1471–4.
Fuster V, Ryden LE, Cannom DS, Crijns HJ, Curtis AB, Ellenbogen KA, et al. ACCF/AHA/HRS focused updates incorporated into the ACC/AHA/ESC 2006 guidelines for the management of patients with atrial fibrillation: a report of the American College of Cardiology Foundation/American Heart Association task force on practice guidelines developed in partnership with the European Society of Cardiology and in collaboration with the European Heart Rhythm Association and the Heart Rhythm Society. J Am Coll Cardiol. 2011;2011:e101–98.
Allessie MA, Boyden PA, Camm AJ, Kleber AG, Lab MJ, Legato MJ, et al. Pathophysiology and prevention of atrial fibrillation. Circulation. 2001;103:769–77.
Andalib A, Brugada R, Nattel S. Atrial fibrillation: evidence for genetically determined disease. Curr Opin Cardiol. 2008;23:176–83.
Clauss S, Kaab S. Is Pitx2 growing up? Circ Cardiovasc Genet. 2011;4:105–7.
Ellinor PT, Lunetta KL, Albert CM, Glazer NL, Ritchie MD, Smith AV, et al. Meta-analysis identifies six new susceptibility loci for atrial fibrillation. Nat Genet. 2012;44:670–5.
Gudbjartsson DF, Arnar DO, Helgadottir A, Gretarsdottir S, Holm H, Sigurdsson A, et al. Variants conferring risk of atrial fibrillation on chromosome 4q25. Nature. 2007;448:353–7.
Lubitz SA, Sinner MF, Lunetta KL, Makino S, Pfeufer A, Rahman R, et al. Independent susceptibility markers for atrial fibrillation on chromosome 4q25. Circulation. 2010;122:976–84.
Kirchhof P, Kahr PC, Kaese S, Piccini I, Vokshi I, Scheld HH, et al. PITX2c is expressed in the adult left atrium, and reducing Pitx2c expression promotes atrial fibrillation inducibility and complex changes in gene expression. Circ Cardiovasc Genet. 2011;4:123–33.
Campione M, Steinbeisser H, Schweickert A, Deissler K, van Bebber F, Lowe LA, et al. The homeobox gene Pitx2: mediator of asymmetric left-right signaling in vertebrate heart and gut looping. Development. 1999;126:1225–34.
Campione M, Ros MA, Icardo JM, Piedra E, Christoffels VM, Schweickert A, et al. Pitx2 expression defines a left cardiac lineage of cells: evidence for atrial and ventricular molecular isomerism in the iv/iv mice. Dev Biol. 2001;231:252–64.
Christoffels VM, Smits GJ, Kispert A, Moorman AF. Development of the pacemaker tissues of the heart. Circ Res. 2010;106:240–54.
Furtado MB, Biben C, Shiratori H, Hamada H, Harvey RP. Characterization of Pitx2c expression in the mouse heart using a reporter transgene. Dev Dyn. 2011;240:195–203.
Franco D, Christoffels VM, Campione M. Homeobox transcription factor Pitx2: The rise of an asymmetry gene in cardiogenesis and arrhythmogenesis. Trends Cardiovasc Med. 2014;24:23–31.
Chinchilla A, Daimi H, Lozano-Velasco E, Dominguez JN, Caballero R, Delpon E, et al. PITX2 insufficiency leads to atrial electrical and structural remodeling linked to arrhythmogenesis. Circ Cardiovasc Genet. 2011;4:269–79.
Wang J, Klysik E, Sood S, Johnson RL, Wehrens XH, Martin JF. Pitx2 prevents susceptibility to atrial arrhythmias by inhibiting left-sided pacemaker specification. Proc Natl Acad Sci U S A. 2010;107:9753–8.
Tao Y, Zhang M, Li L, Bai Y, Zhou Y, Moon AM, et al. Pitx2, an atrial fibrillation predisposition gene, directly regulates ion transport and intercalated disc genes. Circ Cardiovasc Genet. 2014;7:23–32.
Wasserman WW, Sandelin A. Applied bioinformatics for the identification of regulatory elements. Nat Rev Genet. 2004;5:276–87.
Heintzman ND, Hon GC, Hawkins RD, Kheradpour P, Stark A, Harp LF, et al. Histone modifications at human enhancers reflect global cell-type-specific gene expression. Nature. 2009;459:108–12.
Benjamin EJ, Rice KM, Arking DE, Pfeufer A, van Noord C, Smith AV, et al. Variants in ZFHX3 are associated with atrial fibrillation in individuals of European ancestry. Nat Genet. 2009;41:879–81.
Ellinor PT, Lunetta KL, Glazer NL, Pfeufer A, Alonso A, Chung MK, et al. Common variants in KCNN3 are associated with lone atrial fibrillation. Nat Genet. 2010;42:240–4.
Gudbjartsson DF, Holm H, Gretarsdottir S, Thorleifsson G, Walters GB, Thorgeirsson G, et al. A sequence variant in ZFHX3 on 16q22 associates with atrial fibrillation and ischemic stroke. Nat Genet. 2009;41:876–8.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.
Visel A, Minovitsky S, Dubchak I, Pennacchio LA. VISTA Enhancer Browser–a database of tissue-specific human enhancers. Nucleic Acids Res. 2007;35:D88–92.
Volkmann BA, Zinkevich NS, Mustonen A, Schilter KF, Bosenko DV, Reis LM, et al. Potential novel mechanism for Axenfeld-Rieger syndrome: deletion of a distant region containing regulatory elements of PITX2. Invest Ophthalmol Vis Sci. 2011;52:1450–9.
Yee SP, Rigby PW. The regulation of myogenin gene expression during the embryonic development of the mouse. Genes Dev. 1993;7:1277–89.
Campbell RE, Tour O, Palmer AE, Steinbach PA, Baird GS, Zacharias DA, et al. A monomeric red fluorescent protein. Proc Natl Acad Sci U S A. 2002;99:7877–82.
Claycomb WC, Lanson Jr NA, Stallworth BS, Egeland DB, Delcarpio JB, Bahinski A, et al. HL-1 cells: a cardiac muscle cell line that contracts and retains phenotypic characteristics of the adult cardiomyocyte. Proc Natl Acad Sci U S A. 1998;95:2979–84.
Horsthuis T, Houweling AC, Habets PE, de Lange FJ, el Azzouzi H, Clout DE, et al. Distinct regulation of developmental and heart disease-induced atrial natriuretic factor expression by two separate distal sequences. Circ Res. 2008;102:849–59.
Yeom YI, Fuhrmann G, Ovitt CE, Brehm A, Ohbo K, Gross M, et al. Germline regulatory element of Oct-4 specific for the totipotent cycle of embryonal cells. Development. 1996;122:881–94.
Kwasnieski JC, Fiore C, Chaudhari HG, Cohen BA. High-throughput functional testing of ENCODE segmentation predictions. Genome Res. 2014;24:1595–602.
Rebollo R, Romanish MT, Mager DL. Transposable elements: an abundant and natural source of regulatory sequences for host genes. Annu Rev Genet. 2012;46:21–42.
de Souza FS, Franchini LF, Rubinstein M. Exaptation of transposable elements into novel cis-regulatory elements: is the evidence always strong? Mol Biol Evol. 2013;30:1239–51.
Piedra ME, Icardo JM, Albajar M, Rodriguez-Rey JC, Ros MA. Pitx2 participates in the late phase of the pathway controlling left-right asymmetry. Cell. 1998;94:319–24.
Annilo T, Kepp K, Laan M. Natural antisense transcript of natriuretic peptide precursor A (NPPA): structural organization and modulation of NPPA expression. BMC Mol Biol. 2009;10:81.
Shiratori H, Sakuma R, Watanabe M, Hashiguchi H, Mochida K, Sakai Y, et al. Two-step regulation of left-right asymmetric expression of Pitx2: initiation by nodal signaling and maintenance by Nkx2. Mol Cell. 2001;7:137–49.
Shiratori H, Yashiro K, Shen MM, Hamada H. Conserved regulation and role of Pitx2 in situs-specific morphogenesis of visceral organs. Development. 2006;133:3015–25.
Visel A, Akiyama JA, Shoukry M, Afzal V, Rubin EM, Pennacchio LA. Functional autonomy of distant-acting human enhancers. Genomics. 2009;93:509–13.
de Wit E, de Laat W. A decade of 3C technologies: insights into nuclear organization. Genes Dev. 2012;26:11–24.
Simonis M, Kooren J, de Laat W. An evaluation of 3C-based methods to capture DNA interactions. Nat Methods. 2007;4:895–901.
Jin F, Li Y, Dixon JR, Selvaraj S, Ye Z, Lee AY, et al. A high-resolution map of the three-dimensional chromatin interactome in human cells. Nature. 2013;503:290–4.
Mizutani S, Ishii M, Hattori A, Nomura S, Numaguchi Y, Tsujimoto M, et al. New insights into the importance of aminopeptidase A in hypertension. Heart Fail Rev. 2008;13:273–84.
Mitsui T, Nomura S, Okada M, Ohno Y, Kobayashi H, Nakashima Y, et al. Hypertension and angiotensin II hypersensitivity in aminopeptidase A-deficient mice. Mol Med. 2003;9:57–62.
Herrmann S, Layh B, Ludwig A. Novel insights into the distribution of cardiac HCN channels: an expression study in the mouse heart. J Mol Cell Cardiol. 2011;51:997–1006.
Mommersteeg MT, Hoogaars WM, Prall OW, de Gier-de VC, Wiese C, Clout DE, et al. Molecular pathway for the localized formation of the sinoatrial node. Circ Res. 2007;100:354–62.
Sakabe NJ, Savic D, Nobrega MA. Transcriptional enhancers in development and disease. Genome Biol. 2012;13:238.
Visel A, Rubin EM, Pennacchio LA. Genomic views of distant-acting enhancers. Nature. 2009;461:199–205.
Gore-Panter SR, Hsu J, Hanna P, Gillinov AM, Pettersson G, Newton DW, et al. Atrial Fibrillation associated chromosome 4q25 variants are not associated with PITX2c expression in human adult left atrial appendages. PLoS One. 2014;9:e86245.
Arnolds DE, Liu F, Fahrenbach JP, Kim GH, Schillinger KJ, Smemo S, et al. TBX5 drives Scn5a expression to regulate cardiac conduction system function. J Clin Invest. 2012;122:2509–18.
van den Boogaard M, Wong LY, Tessadori F, Bakker ML, Dreizehnter LK, Wakker V, et al. Genetic variation in T-box binding element functionally affects SCN5A/SCN10A enhancer. J Clin Invest. 2012;122:2519–30.
Healey JS, Connolly SJ. Atrial fibrillation: hypertension as a causative agent, risk factor for complications, and potential therapeutic target. Am J Cardiol. 2003;91:9G–14G.
Kato N, Takeuchi F, Tabara Y, Kelly TN, Go MJ, Sim X, et al. Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in east Asians. Nat Genet. 2011;43:531–8.
Mommersteeg MT, Brown NA, Prall OW, de Gier-de VC, Harvey RP, Moorman AF, et al. Pitx2c and Nkx2-5 are required for the formation and identity of the pulmonary myocardium. Circ Res. 2007;101:902–9.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.
Yue F, Cheng Y, Breschi A, Vierstra J, Wu W, Ryba T, et al. A comparative encyclopedia of DNA elements in the mouse genome. Nature. 2014;515:355–64.
Nagy A, Gertsensten M, Vintersten K, Behringer R. Manipulating the mouse embryo: a laboratory manual. 3rd ed. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press; 2003.
Hagege H, Klous P, Braem C, Splinter E, Dekker J, Cathala G, et al. Quantitative analysis of chromosome conformation capture assays (3C-qPCR). Nat Protoc. 2007;2:1722–33.
Palstra RJ, Tolhuis B, Splinter E, Nijmeijer R, Grosveld F, de Laat W. The beta-globin nuclear compartment in development and erythroid differentiation. Nat Genet. 2003;35:190–4.
Acloque H, Wilkinson DG, Nieto MA. In situ hybridization analysis of chick embryos in whole-mount and tissue sections. Methods Cell Biol. 2008;87:169–85.
Schweickert A, Campione M, Steinbeisser H, Blum M. Pitx2 isoforms: involvement of Pitx2c but not Pitx2a or Pitx2b in vertebrate left-right asymmetry. Mech Dev. 2000;90:41–51.
We thank Miguel Torres and members of the Manzanares lab for support and comments; Christine Mummery, José Luis de la Pompa and Joaquín Rodríguez-León for reagents; the CNIC Transgenic Unit for generation of embryos; Stuart Pocock for statistical advice; and Simon Bartlett for English editing. This study was funded by the CNIC Translational Grant Programme (CNIC-08-2009 to MM and DF), the Spanish Ministerio de Economia y Competitividad (grants BFU2011-23083 to MM, BFU2013-41322-P to JLGS, BFU2012-38111 to AA, and CSD2007-00008 to JLGS and MM), the Comunidad Autónoma de Madrid (grant CELLDD-CM to MM), and the Andalusian Government (grant BIO-396 to JLGS). The CNIC is supported by the Spanish Ministerio de Economia y Competitividad and the Pro-CNIC Foundation.
The authors declare that they have no competing interests.
LAA, JLGS, DF and MM designed the study. LAA performed the majority of experiments with contributions from MEA and IR (transgenics); CBC (histology and in situs); CA, ELJ and AA (cell transfections); and AFM and JLGS (3C). LAA and MM wrote the manuscript with input from AA, JLGS and DF. All authors read and approved the final manuscript.
Genomic analysis of 4q25. (A) TAD structure of the 4q25 genomic region (hg19; chr4:110,940,551-113,100,551). Hi-C analysis in three different human cell lines (IMR90, lung fibroblasts; hES, embryonic stem cells; GM12878, lymphoblastoid cells) identifies stable TADs in the gene desert surrounding PITX2, one of which includes both the gene promoter and the AF-associated SNPs (highlighted in green). Hi-C data were obtained from http://yuelab.org/hi-c/  (B) Overlapping fragments from AF3. Detailed view of the genomic landscape and evolutionary conservation of the region surrounding fragment AF3 and of the sub-fragments (AF3.1 to AF3.6) used in this study (hg19; chr4:111,706,648-111,715,384). Legend as in Figure 1.
Results of transgenic experiments.
4q25 regulatory elements do not show cell-type specificity. Activity of 4q25 regulatory elements (AF3, AF3.5 and AF3.6) in mouse HL-1 cardiomyocytes (dark grey), compared to mouse P19 teratocarcinoma (light grey) and human HEK293T embryonic kidney (black) cells. Nppa proximal (NppaPE) and Oct4 distal (Oct4DE) enhancers were used as controls of cell type specificity and the empty pβRFP as a control of basal activity. Data are expressed as mean ± SEM. Statistical significance versus empty pβRFP was calculated with the unpaired Student’s t-test. *P <0.05, **P <0.01 and ***P <0.001.
Assessment of reporter activity in ASE and ASE+AF3.6 transgenic embryos.
Specificity of 3C interactions with Pitx2 and Enpep promoters. (A) Schematic representation of the interacting regions (2R, 6F and 8F; Figure 3) and of the Pitx2 and Enpep genes showing the location of promoter specific anchor primers (Pitx2ab, Pitx2c, Enpep) and control anchor primers located upstream of Pitx2ab (c1), in between Pitx2ab and Pitx2c (c2), downstream of Pitx2c (c3), and upstream (c4) or downstream (c5) of Enpep. (B-D) Normalised 3C interactions, expressed as relative crosslinking frequencies (y-axis), between 2R (B), 6F (C) and 8F (D) fragments and controls (c1-c5) and promoters (Pitx2ab, Pitx2c and Enpep), in atria (red) and ventricles (blue). In each graph, the highest crosslinking frequency values were set to 1. Statistical significance was assessed following one-way ANOVA test of Student-Newman-Keuls. **P <0.01, ***P <0.001. Error bars represent ± SEM.
Primers used in this study.