Human Molecular Genetics Advance Access originally published online on May 25, 2005
Human Molecular Genetics 2005 14(14):1965-1975; doi:10.1093/hmg/ddi202
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Published by Oxford University Press 2005
Functional dissection of sequence-specific NKX2-5 DNA binding domain mutations associated with human heart septation defects using a yeast-based system



1Chromosome Stability Section, Laboratory of Molecular Genetics, National Institute of Environmental Health Sciences, NIH, PO Box 12233, Research Triangle Park, NC 27709, USA and 2Department of Drug Research and Medical Biotechnology, Fraunhofer Institute of Toxicology and Experimental Medicine, Nikolai-Fuchs-Strasse 1, D-30625 Hannover, Germany
* To whom correspondence should be addressed. Email: resnick{at}niehs.nih.gov and Juergen.Borlak{at}item.fraunhofer.de
Received March 18, 2005; Revised May 10, 2005; Accepted May 20, 2005
| ABSTRACT |
|---|
|
|
|---|
Human heart development requires an orderly coordination of transcriptional programs, with the homeodomain protein NKX2-5 being one of the key transcription factors required for the differentiation of mesodermal progenitor cells. Indeed, lack of Nkx2-5 in mice arrests heart development prior to looping, resulting in embryonic lethality. There are 28 germline NKX2-5 mutations identified in humans that are associated with congenital heart disease, and we recently reported multiple somatic mutations in patients with complex cardiac malformations. To address the functional consequences of single and multiple mutations of NKX2-5, we developed a functional assay in the budding yeast Saccharomyces cerevisiae, which could determine transactivation capacity and specificity of expressed NKX2-5 alleles towards targeted response element (RE) sequences. We focused on mutants of the third helix, which provides DNA binding specificity, and characterized mutations that were highly associated with either ventricular (VSD) or atrioventricular (AVSD) septal defects. Individual mutants exhibited partial to complete loss of function and differences in transactivation capacity between the various REs. The mutants also exhibited gene dosage rather than dominant effects on transcription. Surprisingly, all AVSD patients (22/23) had a single K183E mutation in the DNA binding domain, which resulted in transcriptional inactivation. None of the VSD patients had this mutation; yet 14/29 had at least one mutation in the third helix leading to either inactivation or reduction of NKX2-5 transactivation. Therefore, mutations of somatic origin in the binding domains of NKX2-5 were associated specifically with AVSD or VSD and resulted in loss of protein function.
| INTRODUCTION |
|---|
|
|
|---|
The homeodomain (HD) protein Nkx2-5 is an evolutionarily conserved transcription factor required for the organogenesis of the heart (1
Several transcription factors cooperate in the regulation of heart development. Results from whole animals and cellular models suggest complex functional interplays and the necessity for proper gene dosage of these regulators (reviewed in 4
6
). Deletion of both Nkx2-5 alleles in the mouse results in lethality due to impaired cardiac looping (7
,8
), whereas hemizygosity can result in hypoplasia of the conduction system with variable penetrance, leading to a mild first-degree atrioventricular (AV) block (2
,9
). A ventricular myocyte-restricted knock-out mouse model further established a role for Nkx2-5 in the formation, maturation and maintenance of the conduction system (10
). Although no structural heart defects were observed during embryonic development, all adult mice exhibited progressive cardiomyopathy with complete heart block and trabecular muscle overgrowth.
In humans, 28 different heterozygous germline NKX2-5 mutations have been identified in patients suffering from congenital heart disease (CHD) (reviewed in 11
). Although most are in sequences affecting the HD, there is no clear genotypephenotype correlation. Similar to what is observed in the heterozygous Nkx2-5 mouse model, progressive AV conduction delay is the common phenotype, but a near complete penetrance of conduction abnormalities (12
) and a diversity of cardiac malformations are observed. Cardiac malformations of patients with NKX2-5 germline mutations were mainly ASD with or without AV block, although there were reports of patient(s) with tetralogy of Fallot, ventricular septal defect (VSD), double-outlet right ventricle, interrupted aortic arch, truncus arteriosus, L-transposition of the great arteries, hypoplastic left heart syndrome and coarctation of the aorta.
The phenotypic discrepancies between CHD associated with NKX2-5 germline mutations in humans and in the haploinsufficient Nkx2-5 mouse models have been ascribed to the dominant negative impact of mutant Nkx2-5 proteins towards wild-type Nkx2-5 or other cardiac transcription factors that act in cooperation and physically associate with Nkx2-5, such as Gata4 and Tbx5 (5
). However, the spectrum of germline human NKX2-5 mutations is diverse in terms of mutation type, position of the affected amino acids and predicted impact of the mutations on proteinprotein interactions (11
,13
,14
). Indeed, there appears to be no common mutation identified in familial or sporadic cases of CHD. Although evidence of dominant negative properties of NKX2-5 mutations has been presented on the basis of transcription assays in model systems (11
,14
), the in vivo relevance of these observations remains to be established. In fact, the transcriptional targets of NKX2-5 that contribute to the heart chamber maturation and the development of the conduction system are largely unknown (5
), and the results in animal models showed weak correspondence between alterations of gene expression in vivo and identification of direct target genes obtained using transfected cell lines, with over-expression of NKX2-5 alone or of other regulators of its transcriptional network (15
).
NKX2-5 genetic heterogeneity and mosaicism in the diseased tissue could be a source of the phenotypic variability that is associated with human CHD, but is not recapitulated in animal models. This hypothesis is based on recent results from this (J. Borlak) laboratory of NKX2-5 mutations identified in diseased tissues from explanted formalin-fixed hearts of patients affected by complex cardiac malformations that included ASD, VSD and AVSD (16
,17
). Particularly interesting was the observation that non-synonymous multiple NKX2-5 mutations were identified in almost all diseased tissue samples. However, there were no mutations in unaffected heart sections of the same patients or lymphocytic DNA from unrelated individuals with CHD. Although mutations in the genes encoding the cardiac transcription factors TBX5 and GATA4 were also found in a fraction of these samples (18
,19
), there did not appear to be a general mutagenesis for other genes coding for cardiac transcription factors such as MEF2C (unpublished data). Indeed, MEF2C maps to chromosome 5 (5q14) close to NKX2-5 (5q34) but did not carry mutations in the sequences encoding its DNA binding domains in 68 hearts and 12 lymphocytic DNA samples of patients with CHD, as well as in 50 lymphocytic DNA samples of normal individuals. Thus, the results were indicative of NKX2-5 mutations of somatic origin, with multiple haplotypes and mosaicism generated through an unknown mechanism of localized hyper-mutagenesis. Although not nearly as extensive as for NKX2-5, multiple somatic mutations were also reported for the HPRT gene in kidney tubular epithelial cells (20
). We have identified an additional 16 cases in which germline NKX2-5 mutations were detected in the affected heart tissues in conjunction with other mutations (17
). Taken together, these results suggest that the wide phenotypic spectrum of congenital heart disease may be due, at least in part, to somatic mutations in transcription factor genes that control heart development in progenitor cells.
To address the functional significance of individual and multiple somatic NKX2-5 mutants in CHD, particularly VSD and AVSD, we created a functional assay for this homeobox transcription factor in a model organism, the budding yeast Saccharomyces cerevisiae. We focused mainly on newly observed somatic NKX2-5 gene mutations affecting amino acids in the third helix of the HD, a region that provides for specific binding to the DNA target response elements (NK-REs). Similar to a system that we had developed to examine the tumor suppressor p53 (21
), the yeast-based assay measures transactivation capacity under conditions of tight, rheostatable regulation of NKX2-5 expression and utilizes a set of reporter constructs under isogenic conditions. Different defined NK-REs derived from the promoters of known NKX2-5 target genes are placed upstream of the reporters to stimulate transcription in an otherwise constant chromatin environment. The tight regulation of expression of the NKX2-5 transcription factor also renders the functional assay sensitive to gene dosage, thus providing for analysis of the dominant negative potential of NKX2-5 mutants when heterozygous with wild-type NKX2-5.
The newly developed functional assay in yeast provides a robust tool for the in vivo assessment of the transactivation capacity of wild-type NKX2-5 and of single or multiple mutants towards a variety of REs. Nearly all AVSD patients had a single mutation which was found to be transcriptionally inactive. In contrast, none of the VSD patients had this mutation. However, there was always at least one mutation in this region that led to inactivation of NKX2-5 transactivation capacity, with some individual mutations exhibiting partial function. Loss-of-function NKX2-5 mutants did not exhibit dominance when co-expressed at equal levels with the wild-type allele, but instead resulted in a gene-dosage reduction in function.
| RESULTS |
|---|
|
|
|---|
NKX2-5 can act as a sequence-specific transcription factor in yeast
The model organism S. cerevisiae has provided opportunities to examine the transactivation capacity of wild-type and mutant forms of various mammalian transcription factors, often in the absence of a yeast homolog (reviewed in 22
Our recent development of a system for in vivo site-directed mutagenesis by oligonucleotides (25
) has greatly improved the ability to modify yeast promoters for transcription studies (21
,26
). By applying this methodology to the analysis of the transactivation capacity of human p53 towards many different REs, we have also recently developed conditions that result in tight regulation of expression of the transcription factor along with rheostatable induction of the GAL1 promoter using variable levels of galactose in the culture medium (21
). In a separate study demonstrating the utility of our approach, the combination of tight regulation of p53 levels and the transactivation analysis towards many REs greatly improved the analysis and classification of p53 mutations associated with cancer (26
).
Drawing upon the experience with p53, we constructed expression vectors pTSGNKX2-5 and pUSGNKX2-5, which are low copy centromeric plasmids, to examine the human homoeodomain NKX2-5 protein in yeast, as described in Materials and Methods. In addition to the HD, which provides for DNA binding, two other domains are conserved in the NK2 class of homeoproteins to which NKX2-5 belongs. The specific domain (NK2-SD) may mediate proteinprotein interactions and appears to have an inhibitory function on sequence-specific transactivation in model systems (27
). The TN domain contains a motif weakly conserved in other transcription factors and may contribute to transactivation (27
). Thus, to establish a system based in yeast to address homoeodomain transcription factors, four different NKX2-5 constructs were tested. They included the wild-type full-length protein, a C-terminal truncation that lacked the NK2-SD and might result in higher activity and also those same two constructs fused to an acidic TAD derived from p53 (Fig. 1A). Given the absence of a well-defined NKX2-5 activation domain and previous reports of weak activity in yeast of some classes of TADs (22
), constructs with the p53 chimeric TAD were developed with the view that they could provide transcription if sequence-specific binding were obtained.
|
For the initial functional analyses of NKX2-5 towards an RE, we chose the distal NK-RE derived from the atrial natriuretic factor (ANF) promoter that has been used in transactivation assays in mammalian cells (28
The ADE2 gene was chosen as the NKX2-5 reporter because it provides for a simple visual transactivation assay based on change in color of the yeast colonies from red to pink or white, depending on level of transactivation (21
). We also used a quantitative assay based on the firefly luciferase reporter, as previously described (26
). The expression levels of NKX2-5 proteins were adjusted in the functional assays by varying the amounts of galactose inducer (21
).
As shown in Figure 1B and C, wild-type NKX2-5 can greatly stimulate transcription in yeast. This activity required a wild-type NKX2-5 HD and the presence of an NKX-RE in the promoter. This is the first report of a human HD protein functioning in transcription in yeast. Transcription also required the acidic TAD provided in the chimeric construct. The C-terminal truncation that lacks the NK2-SD exhibited up to 3-fold higher activity when compared with the full-length construct (Fig. 1C), suggesting the conservation in the yeast system of the inhibitory effect seen in mammalian cells (29
). At high expression, SD-truncated NKX2-5 showed low levels of transactivation even in the absence of the exogenous activation domain, whereas the chimeric construct is at least 50 times more active.
Our results show that the HD-containing protein NKX2-5 can act as a sequence-specific transcription factor in yeast but it is intrinsically very weak. This observation suggests that it would be possible to use the system to screen for additional transcription factors or cofactors that can physically associate and cooperate with NKX2-5 function, as in the case of the cardiac transcription factors GATA4 and TBX5 in human cells (30
,31
). As our aim was to study the impact of HD NKX2-5 mutations on sequence-specific transactivation, in the subsequent analyses, we took advantage of the robust activity of the chimeric full-length NKX2-5 with an exogenous acidic TAD and also examined transactivation at several different REs.
Impact of RE sequence variation on NKX2-5 transactivation capacity
We explored the effect that variation in NK-RE sequence and structure can have on NKX2-5 transactivation capacity. Besides allowing comparison of results between an in vivo functional assay and an in vitro DNA binding measurements with naked DNA (29
), this analysis was expected to generate a set of reporter strains for the subsequent analysis of NKX2-5 mutations.
The distal ANF-RE, which provides for a dimer-binding site, showed relatively weak activity. Two copies of this RE could mediate only low levels of transactivation as evidenced by the appearance of pink colonies in the ADE2-based assay at the high NKX2-5 expression obtained with 0.5% galactose (data not shown). Increasing the number of ANF-REs resulted in cooperative interaction. Seven copies of the RE (i.e. seven adjacent NKX2-5 dimer-binding sites) were nearly 30 times more active (pink colonies at 0.004% galactose, Fig. 1B). Each monomer-binding site in the distal ANF-RE has one mismatch from the NK-RE consensus TYAAGTG (27
) (Table 1). To explore the impact of single nucleotide changes in the NK-RE on NKX2-5 functionality, we tested six sequence variants that are present in proposed NKX2-5 target promoters. On the basis of the number of distal ANF-REs required for transactivation, each sequence motif was examined by cloning four adjacent monomer-binding sites into the ADE2 reporter construct in direct and in inverse orientation. None of the sequences was capable of mediating NKX2-5 transactivation in direct orientation (Table 1). Instead, stimulation of transcription was detected with two sequence motifs placed in inverse orientation, thus forming two palindromic dimer-binding sites. The endothelin converting enzyme (ECE) consensus sequence was 14 times stronger than the distal ANF-RE of the same structure. There was no transactivation with nucleotide changes at the following position of the consensus sequence: first (T>G, pANF), sixth (T>A, dANF-b or T>G, pDIO2) and seventh (G>T, dDIO2). The C>G mismatch (dANF-a) at the second position retained
30% function. A combination of a weak and an inactive sequence motif explains the very low activity (
7% of the consensus) of the distal ANF-a+b RE. (Note: GTG is an E-box half-site, frequently contacted by bHLH TFs. Most likely the thymine is contacted by the NKX2-5 HD.)
|
Although based on a relatively small sampling of NK-RE sequence variants, our results suggest that the REs that are active in our in vivo functional assay are only a subset of the sequences than can bind the NKX2-5 protein in vitro (29
Functional consequences of novel somatic NKX2-5 mutations associated with cardiac malformations
Although not identical in structure to the NK-REs found in the human promoters, palindromic ECE and ANF-REs provide opportunities with the yeast system to analyze the functional consequences of NKX2-5 DNA binding mutations. On the basis of our previous studies of mutant forms of the p53 transcription factor, where the impact of mutations on transactivation capacity differed between RE sequences, we chose to measure transactivation by NKX2-5 mutant proteins using various RE sequences and number of dimer-binding sites.
For this study, we focused on NKX2-5 gene mutations of both somatic and germline origin that were detected in VSD and AVSD heart malformations. The mutations resulted in amino acid changes at the third helix of the HD which mediates specific interactions with the major groove of the DNA at the NK-RE site. As molecular analysis of the NKX2-5 gene from diseased tissues revealed the presence of multiple mutants and different haplotypes (17
), we were also interested in determining the functional consequences of individual and multiple mutations within a given haplotype. This information might provide clues to possible selection processes that give rise to the disease. Six different haplotypes from three VSD patients were analyzed (Fig. 2 and Table 2). We also investigated the mutant K183E (corresponding to amino acid 46 of the HD) that was highly associated with AVSD (observed in 22 out of 23 patients) and the stop codon mutant Q187X (corresponding to amino acid 50 of the HD). Except for the T178M mutation, which has also been found in the germline (32
,33
), none of the other HD mutants given in Table 2 have been analyzed in a functional assay. Further, we analyzed A219V which is located in the NK2-SD and is also a germline mutation (34
).
|
|
The single and multiple mutants were constructed in the pTSG-TAD::NKX2-5 vector backbone (see Materials and Methods) and the functional analysis was conducted using the quantitative luciferase reporter at moderate levels of NKX2-5 allele expression (0.12% galactose). The normalized average light units for four mutants and the controls are presented in Figure 2. The frequently appearing AVSD-associated K183E amino acid change resulted in complete loss of transactivation capacity, whereas the double mutant K192T::K194R found in one haplotype from VSD patient E04 (17
10% residual function. When analyzed separately, both mutations showed a milder effect, particularly in the case of K194R, which can exhibit nearly 75% functionality depending on the RE tested. Thus, the two mutations can contribute to the phenotype of the disease-associated haplotype. For the case of ANF (14 dimers), there was a clear synergy in the reduction of transactivation capacity by the double mutant when compared with the single mutations, demonstrating that multiple mutations can act in concert to reduce NKX2-5 function.
In all, we tested 20 different NKX2-5 mutants with four different sets of REs. The average transactivation capacities relative to wild-type NKX2-5 are summarized in Table 2. Disease-associated haplotypes are indicated by the patient identifier in the first column and the disease name in the third column. Overall, disease-associated somatic NKX2-5 mutants were non-functional or retained little activity as haplotypes, generally below the effect obtained with hemizygosity (discussed subsequently). Almost all the individual mutations exhibited some degree of functional alteration, particularly when tested with the strain containing the most sensitive reporter (two NKX dimer-binding sites of the ECE consensus sequence ECE 2d). A mild effect was seen with the frequent mutation A219V that affects an amino acid in the NK2-SD regulatory domain and is not directly involved in sequence-specific DNA binding. The negative phenotype of the triple mutant haplotypes from patient D03 appears to be due to the T178M mutation, an amino acid change that causes a loss-of-function. This germline mutation, also found in familial CHD patients, was shown to have reduced in vitro DNA binding resulting in loss of transactivation in human cells (14
).
It is important to note that the impact of various mutations and combination of mutations can be dependent on the particular RE used to address functionality. This observation is illustrated in Figure 3, where the relative transactivation capacities of NKX2-5 mutations exhibiting partial function are presented in the form of a line graph to emphasize the different patterns of response with the various REs. The overall pattern of activity was similar for the ECE (2d)- and ECE (4d)-REs, whereas responses were more variable and dependent on specific NKX2-5 mutations with the ANF-REs, particularly with the high activity ANF (14d), which is generally less sensitive to NKX2-5 functional mutants. For example, K192R and A219V, either together or individually, have only a modest decrease in transactivation for ANF (14d), but the combined mutations have a strong impact when measured against the ECE-REs. In addition, some mutants (e.g. K192T, A219V) led to changes in spectrum of activation, not simply a general reduction, with ANF (14d) being particularly affected. The pattern of transactivation responses also shows how different amino acid changes at the same position can lead to significant differences in transactivation, as found for K192R versus K192T. Overall, these results are reminiscent of those previously reported for mutations in the tumor suppressor p53 and suggest that mutations in the sequence binding domain can have markedly different effects, depending on the RE being transactivated and the specific mutation. This is particularly important for genes such as NKX2-5 that can regulate many genes from various REs. As for the case of p53, we suggest that NKX2-5 may also be a master gene of diversity (26
).
|
Addressing dominance of wild-type and mutant NKX2-5
The yeast-based functional assay can be easily adapted to the analysis of dominance of mutant transcription factor alleles when heterozygous with the wild-type. The tight regulation of protein expression obtained with the inducible GAL1 promoter combined with the availability of reporters with different responsiveness also provides the opportunity to develop functional assays that are sensitive to gene dosage. We constructed vector pUSG-TAD::NKX2-5 that expresses wild-type NKX2-5 fused to the acidic TAD and contains the URA3 selectable marker (see Materials and Methods). This plasmid can then be selected along with the original pTSG-TAD::NKX2-5 vectors (TRP1 marker) that contain either wild-type NKX2-5 or mutants. We chose reporter strain with the ECE (2d)-RE because, overall, it showed the highest discrimination between NKX2-5 mutants.
The results of the functional assay using moderate co-expression of two NKX2-5 alleles are presented in Table 3. The induction of the luciferase reporter was clearly sensitive to the copy number of wild-type NKX2-5. As expected, expression from the two vectors was nearly identical as they only differ in the selection marker. Selection of both clones resulted in nearly double the amount of transactivation when the two wild-type NKX2-5 proteins were expressed. There was no apparent dominance by the mutants because co-expression of loss-of-function mutants (from pTSG-TAD vectors) with wild-type NKX2-5 (from pUSG-TAD::NKX2-5) was functionally equivalent to the expression of the wild-type allele only, and a partial function allele showed an additive effect with the wild-type. The lack of dominance for these particular VSD-associated haplotypes may also result from an inability of the mutant proteins to form heterodimers with the wild-type protein (28
).
|
| DISCUSSION |
|---|
|
|
|---|
Nkx2-5 is a transcription factor that is required during cardiac development, and mutations in the human NKX2-5 gene have been associated with CHD. Mutations of the NKX2-5 gene are typically identified in lymphocytic DNAs and most are familial. Among a total of 28 reported NKX2-5 germline mutations, 10 are located in the HD. Of these, eight are missense (R142C, L171P, T178M, Q187H, N188K, R189G, R190H, Y191C) and two are nonsense (Q149X, Q170X) mutations. Five are located in the third helix (Q187H, N188K, R189G, R190H, Y191C), two are in the second helix (Q170X, L171P), one is in the first helix (Q149X), one is in the loop between the second and third helix (T178M) and the remaining mutation is in the NLS region (R142C). The appearance of NKX2-5 mutations in the lymphocytic DNA of patients with non-familial heart defects is uncommon. Among 608 patients, only 18 (3%) revealed an NKX2-5 mutation in their lymphocytic DNA, with 12 distinct mutations (35
Previous functional analysis of the eight missense HD mutations showed localization of all mutant proteins to the nuclei, but some exhibited anomalous nuclear distribution (e.g. T178M) (11
). All mutants had exhibited markedly decreased DNA binding and reduced transcriptional activation, whereas interaction with transcriptional partners varied. We have undertaken functional studies on mutations affecting the HD of NKX2-5. Unlike the situation for single mutations being identified in lymphocytic DNA, the NKX2-5 gene from diseased cardiac tissues of patients with complex cardiac malformations typically contained multiple mutations (17
). Neither have the reasons for the multiple mutations been established nor has the role of the individual mutations in the appearance of the disease been directly addressed. It remains to be determined whether heart disease associated with the single germinal mutations is due only to the single mutations or, alternatively, whether additional mutations arise in the heart that amplify the effect of the germinal mutations. Further, mutations in other transcription factor genes including GATA4, TBX5, FOG2 and CFC1 have been associated with CHD in humans, and there is emerging evidence from mutational studies for additional transcription factors to have a role in CHD (37
). Recently, three individuals with common arterial trunk (CAT) were found to be carrying a homozygous c.451T>C (F151L) mutation in the HD of NKX2-6 (38
). Similar to other germline mutations in cardiac-specific transcription factor genes, F151L was familial and not detected in 10 sporadic cases of CAT. Therefore, studies in which there is only an examination of lymphocytic DNA may not reveal the molecular basis of CHD.
The assay we developed, which is the first to examine HD transcription factors in yeast, provides a convenient and powerful system for the functional characterization of single and multiple mutations to elucidate mechanisms of CHD. Noteworthy among the mutations examined with the yeast system is K183E, which was detected in 22 of 23 patients exhibiting AVSD, in seven of 16 with ASD, but was not detected among 29 VSD patients. The single K183E mutation led to complete loss of transactivation activity. In constrast to AVSD, where K183E was the only mutation in the third helix and specific for AVSD, there were four mutations in patients with VSD (Q187X, K192R, K192T, K194R). Two patients had combined K192R and K194R, which were located in the same allele. In VSDs, 14/29 had at least one mutation in the third helix, which led to inactivaction or reduction of transactivation.
Dissection of the triple mutant haplotypes revealed that the individual K192R, K192T and K194R reduced transactivation to a level that ranged from 30 to
80% of wild-type; some pair wise combinations of mutations had a much more dramatic effect (i.e. K192T, K194R). Also, there were examples where changes in transactivation differed greatly between REs, as found for ECE (2d) versus ANF (14d) with the K192R, A219V mutant protein. Although the transition from these results with a model system to clinical manifestations must be done with reservation, these observations with various single and multiple mutants and with different REs suggest that the various mutants might give rise to a variety of cardiac phenotypes.
In this regard, it is interesting that T178M and A219V are common mutations in three haplotypes from the same patient that also contain a third distinct mutation, suggesting that they originated prior to the third mutation (Table 1 and patient D03 in Table 2). T178M is a loss-of-function mutant in our assay, as well as in mammalian cell-based studies (13
,14
). A219V showed instead a mild reduction in function for all REs examined. T178 is located in the loop between the second and third helix of the HD, whereas A219 is located in the NK2-SD domain. Both T178M and A219V changes have been identified as germline CHD-associated mutations (33
,34
). Interestingly, there is preliminary evidence for mutations within NK2-SD contributing to CHD (http://cardiogenomics.med.harvard.edu/groups/proj1/pages/nk2-sd_home.html). On the basis of our molecular identification of mutations in disease-associated tissues (17
) along with functional analyses of mutations, we propose that a germinal A219V mutation is a risk factor that when combined with other somatic mutations can increase the likelihood of CHD.
As the double mutant T178M/A219V is loss-of-function, a disease selection hypothesis would imply that the subsequent mutation confers additional negative impact on NKX2-5 functions. Given the mosaic nature of the diseased tissues, the presence of functional NKX2-5 proteins in individual cells cannot be excluded (17
). Hence, the third mutations in the haplotypes might confer a dominant negative phenotype to the already transactivation defective mutant protein when engaged in dimers with functional NKX2-5 proteins. However, we were unable to detect dominance from the triple mutant protein in the yeast-based assay (Table 3). Alternatively, the multiple mutations may disrupt functions of NKX2-5 separate from the sequence-specific DNA binding. Along this line, it was shown that the NKX2-5 region between amino acids 182 and 199 is responsible for the physical interaction with the cardiac transcription factor GATA4 and that the K194I mutants abolished the interaction in vitro (28
). Physical interaction between GATA4 and Nkx2-5 can result in transcriptional synergy (31
,39
41
). Recently, three germline GATA4 gene mutations have been associated with CHD (30
,42
): missense G296S, immediately after the C-terminal finger; a frameshift E359del and another frameshift mutation leading to a premature stop codon at amino acid 403. In addition, 12 somatic GATA4 mutations affecting the zinc fingers were detected in our patient cohort (19
). Specifically, mutant C292R, which affects a zinc coordinating cysteine in the C-terminal zinc finger, was detected in 19 of 29 patients with VSD, suggesting that targeting the functional synergy between GATA4 and NKX2-5 may result in septation defects. Taken together, these results suggest that disrupting the interaction between NKX2-5 and GATA4 can contribute to the clinical manifestations of CHD.
In conclusion, the yeast-based functional assay provides opportunities to examine the consequences of single and multiple mutations in the master regulatory NKX2-5 protein that is essential in normal heart development. Overall, the functional dissection of individual mutations contained in different haplotypes from the same patient suggested the selection in the disease tissues of NKX2-5 variants impaired both in sequence-specific DNA binding and in physical interaction with other cardiac transcription factors. The approach that we have taken can be adapted to the direct cloning of separated NKX2-5 alleles from genomic DNA isolated from heart tissue, especially because the gene comprises only two large exons. Similar approaches have been developed for the detection of p53 mutations in tumor tissue (43
,44
). This would facilitate the identification of mutations with reduced functionality in mixed samples. The system could also be extended to the analysis of other transcription factors involved in heart development whose function is altered in CHD, such as GATA4 and TBX5, and their functional interplay with NKX2-5 in the activation of the transcriptional network of heart development.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Cardiac malformations, genomic DNA isolation and detection of NKX2-5 mutations
The characterization of malformed hearts, genomic DNA isolation and mutation detection in NKX2-5 gene have been described previously (16
Genomic DNA was isolated with NucleoSpin Tissue Kit (Macherey-Nagel, Düren Germany). NKX2-5 fragments were amplified by PCR reaction and directly sequenced with specific primers. Mutations were confirmed by PCRRFLP assay, when possible, or by cloning of the PCR fragments into TOPO TA Cloning Kit for re-sequencing (Invitrogen, Karlsruhe, Germany).
Cloning of wild-type NKX2-5 into yeast expression vectors
Total RNA was isolated from normal heart tissue with RNAeasy Kit and RTPCR was carried out with Omniscript RT Kit, according to the manufacturer's procedure (Qiagen, Hilden, Germany). The coding region of NKX2-5 was amplified using the same PCR conditions described for formalin-fixed samples (16
) and PCR primers with incorporated BamHI and XhoI sites. PCR fragments were digested with BamHI and XhoI and cloned into the pTSG and the pTSG-TAD vectors using standard procedure. pTSG is a centromeric yeast expression vector based on the inducible GAL1 promoter (Fig. 1A) (21
). pTSG-TAD is identical to pTSG but it contains the TAD derived from the human transcription factor p53 downstream GAL1 and allows the construction of A-terminal fusions with full-length NKX2-5. Cloning of the NKX2-5-
-SD was performed using a similar approach. DNA from the resulting plasmids pTSG-NKX2-5, pTSG-TAD::NKX2-5, pTSG-NKX2-5-
-SD and pTSG-TAD::NKX2-5-
-SD were isolated using QIAprep Spin Miniprep Kit (Qiagen) and inserts were analyzed by restriction digestion and confirmed by DNA sequencing. To allow simultaneous selection of two NKX2-5 expression vectors for the dominance assay, plasmid pUSG-TAD::NKX2-5 was constructed by replacing the TRP1 selection marker in pTSG-TAD::NKX2-5 with URA3 using standard cloning procedures.
Cloning of wild-type NKX2-5 into yeast expression vector
Total RNA was isolated from normal heart tissue with RNAeasy Kit and RTPCR was carried out with Omniscript RT Kit, according to the manufacturer's procedure (Qiagen). The coding region of NKX2-5 was amplified using the same PCR conditions described for formalin-fixed samples (16
) and PCR primers with incorporated BamHI and XhoI sites. PCR fragments were digested with BamHI and XhoI and cloned into the pTSG and the pTSG-TAD vectors using standard procedure. pTSG is a centromeric yeast expression vector based on the inducible GAL1 promoter (Fig. 1A) (21
). pTSG-TAD is identical to pTSG but it contains the TAD derived from the human transcription factor p53 downstream GAL1 and allows the construction of A-terminal fusions with full-length NKX2-5. Cloning of the NKX2-5
-SD construct, with deletion of the C-terminal domain starting from amino acid 208, was performed using a similar approach. DNA from the resulting plasmids pTSG-NKX2-5, pTSG-TAD::NKX2-5, pTSG-NKX2-5-
-SD and pTSG-TAD::NKX2-5-
-SD were isolated using QIAprep Spin Miniprep Kit (Qiagen) and inserts were analyzed by restriction digestion and confirmed by DNA sequencing. To allow simultaneous selection of two NKX2-5 expression vectors for the dominance assay, plasmid pUSG-TAD::NKX2-5 was constructed by replacing the TRP1 selection marker in pTSG-TAD::NKX2-5 with URA3 using standard cloning procedures.
Preparation of yeast expression vectors with different NKX2-5 mutations
Multiple NKX2-5 mutations (e.g. triplets) were obtained from TOPO-TA cloned PCR fragments as part of the verification procedure for NKX2-5 variants from fixed heart samples (previously mentioned). Multiple mutations were brought into the wild-type NKX2-5 plasmids by standard digestion (KpnI/BsrI) and ligation protocols. Sequential dissections of triple NKX2-5 mutants to obtain double and single mutations were carried out with the GeneTailorTM Site-Directed (SD) Mutagenesis System (Invitrogen). To facilitate amplification, NKX2-5 multiple mutants cloned in the small pBluescript plasmid (Stratagene, 3 kb) were used as templates for GeneTailor SD-mutagenesis. Q-solution (Qiagen) was also added during amplification step. Clones containing only the desired NKX2-5 mutation (s) were identified by DNA sequencing.
Phenotypic and quantitative yeast-based NKX2-5 functional assays
The construction of NKX2-5 reporter strains was performed using the master reporter strains yAFM-ICORE (ADE2 color reporter) and yLFM-ICORE (luciferase reporter), as previously described (26
). Briefly, the ICORE cassette contains the counter-selectable KLURA3, the reporter KanMX4 and the I-SceI endonuclease gene under control of the GAL1 promoter and its unique 18 nt recognition site. Induction of I-SceI on galactose leads to the generation of a single DNA double-strand break at the ICORE site, resulting in a dramatic increase in homologous DNA targeting by oligonucleotides along with coincident loss of the ICORE (25
). The cassette is integrated upstream of a minimal cyc1 promoter which was placed to regulate the expression of the reporters at the ADE2 locus on chromosome XV. To generate NKX2-5 responsive strains, the ICORE was replaced with oligonucleotides containing the desired NK-REs, according to the delitto-perfetto protocol (25
). NK-REs were derived from the promoters of proposed and validated NKX2-5 targets. The sequence and the structure of the NK-REs examined in this study are reported in Table 1. The ANF-, ECE- and DIO2-REs are derived from the promoters of the ANF, ECE and isoform 1b and the iodothyronine deiodinase genes, respectively. Details on strain construction are available upon request. The ADE2-based visual transactivation assay was performed on plates with low adenine (5 mg/l) to allow color detection (21
). Transactivation of ADE2 results in a change of the colony color from red to pink or white. Differences in transactivation capacities with the different REs are estimated based on the amount of NKX2-5 required for transactivation, which is dependent on the amount of galactose inducer, and is based on the approach we developed to analyze p53 transactivation capacity (21
). Luciferase assays were performed using protein extracts obtained from the mechanical lysis of yeast cells grown in liquid cultures, as previously described (26
). The averages for induced light units and the standard deviations of at least three biological replicates are presented in Figures 1C and 2 and reported in Tables 2 and 3.
| ACKNOWLEDGEMENTS |
|---|
We thank Annika Roskowetz and Andreas Hiemisch for technical support, and we appreciate the financial support to J.B. (grant no. 25A.5-76251-99-3/00) of the Lower Saxony Ministry of Science and Culture, Germany.
Conflict of Interest statement. None declared.
| FOOTNOTES |
|---|
Present address: Laboratory of Experimental Oncology B, Department of Translational Oncology, National Cancer Research Institute, IST, Largo R. Benzi X, 16132 Genoa, Italy.
These two authors contributed equally as Senior Authors. ![]()
| REFERENCES |
|---|
|
|
|---|
- Harvey, R.P., Lai, D., Elliott, D., Biben, C., Solloway, M., Prall, O., Stennard, F., Schindeler, A., Groves, N., Lavulo, L. et al. (2002) Homeodomain factor Nkx2-5 in heart development and disease. Cold Spring Harb. Symp. Quant. Biol., 67, 107114.[CrossRef][Web of Science][Medline]
- Jay, P.Y., Harris, B.S., Maguire, C.T., Buerger, A., Wakimoto, H., Tanaka, M., Kupershmidt, S., Roden, D.M., Schultheiss, T.M., O'Brien, T.X. et al. (2004) Nkx2-5 mutation causes anatomic hypoplasia of the cardiac conduction system. J. Clin. Invest., 113, 11301137.[CrossRef][Web of Science][Medline]
-
Watada, H., Mirmira, R.G., Kalamaras, J. and German, M.S. (2000) Intramolecular control of transcriptional activity by the NK2-specific domain in NK-2 homeodomain proteins. Proc. Natl Acad. Sci. USA, 97, 94439448.
[Abstract/Free Full Text] - Srivastava, D. and Olson, E.N. (2000) A genetic blueprint for cardiac development. Nature, 407, 221226.[CrossRef][Medline]
-
Bruneau, B.G. (2002) Transcriptional regulation of vertebrate cardiac morphogenesis. Circ. Res., 90, 509519.
[Abstract/Free Full Text] - Harvey, R.P. (2002) Patterning the vertebrate heart. Nat. Rev. Genet., 3, 544556.[CrossRef][Web of Science][Medline]
-
Lyons, I., Parsons, L.M., Hartley, L., Li, R., Andrews, J.E., Robb, L. and Harvey, R.P. (1995) Myogenic and morphogenetic defects in the heart tubes of murine embryos lacking the homeo box gene Nkx2-5. Genes Dev., 9, 16541666.
[Abstract/Free Full Text] -
Tanaka, M., Schinke, M., Liao, H.S., Yamasaki, N. and Izumo, S. (2001) Nkx2.5 and Nkx2.6, homologs of Drosophila tinman, are required for development of the pharynx. Mol. Cell. Biol., 21, 43914398.
[Abstract/Free Full Text] - Jay, P.Y., Harris, B.S., Buerger, A., Rozhitskaya, O., Maguire, C.T., Barbosky, L.A., McCusty, E., Berul, C.I., O'Brien, T.X., Gourdie, R.G. and Izumo, S. (2004) Function follows form: Cardiac conduction system defects in Nkx2-5 mutation. Anat. Rec., 280A, 966972.
- Pashmforoush, M., Lu, J.T., Chen, H., Amand, T.S., Kondo, R., Pradervand, S., Evans, S.M., Clark, B., Feramisco, J.R., Giles, W. et al. (2004) Nkx2-5 pathways and congenital heart disease; loss of ventricular myocyte lineage specification leads to progressive cardiomyopathy and complete heart block. Cell, 117, 373386.[CrossRef][Web of Science][Medline]
-
Kasahara, H. and Benson, D.W. (2004) Biochemical analyses of eight NKX2.5 homeodomain missense mutations causing atrioventricular block and cardiac anomalies. Cardiovasc. Res., 64, 4051.
[Abstract/Free Full Text] - Benson, D.W., Silberbach, G.M., Kavanaugh-McHugh, A., Cottrill, C., Zhang, Y., Riggs, S., Smalls, O., Johnson, M.C., Watson, M.S., Seidman, J.G. et al. (1999) Mutations in the cardiac transcription factor NKX2.5 affect diverse cardiac developmental pathways. J. Clin. Invest., 104, 15671573.[Web of Science][Medline]
-
Zhu, W., Shiojima, I., Hiroi, Y., Zou, Y., Akazawa, H., Mizukami, M., Toko, H., Yazaki, Y., Nagai, R. and Komuro, I. (2000) Functional analyses of three Csx/Nkx-2.5 mutations that cause human congenital heart disease. J. Biol. Chem., 275, 3529135296.
[Abstract/Free Full Text] - Kasahara, H., Lee, B., Schott, J.J., Benson, D.W., Seidman, J.G., Seidman, C.E. and Izumo, S. (2000) Loss of function and inhibitory effects of human CSX/NKX2.5 homeoprotein mutations associated with congenital heart disease. J. Clin. Invest., 106, 299308.[Web of Science][Medline]
- Kasahara, H., Wakimoto, H., Liu, M., Maguire, C.T., Converso, K.L., Shioi, T., Huang, W.Y., Manning, W.J., Paul, D., Lawitts, J., Berul, C.I. and Izumo, S. (2001) Progressive atrioventricular conduction defects and heart failure in mice expressing a mutant Csx/Nkx2.5 homeoprotein. J. Clin. Invest., 108, 189201.[CrossRef][Web of Science][Medline]
-
Reamon-Buettner, S.M., Hecker, H., Spanel-Borowski, K., Craatz, S., Kuenzel, E., and Borlak, J. (2004) Novel NKX2-5 mutations in diseased heart tissues of patients with cardiac malformations. Am. J. Pathol., 164, 21172125.
[Abstract/Free Full Text] -
Reamon-Buettner, S.M. and Borlak, J. (2004) Somatic NKX2-5 mutations as a novel mechanism of disease in complex congenital heart disease. J. Med. Genet., 41, 684690.
[Abstract/Free Full Text] - Reamon-Buettner, S.M. and Borlak, J. (2004) TBX5 mutations in Non-HoltOram Syndrome (HOS) malformed hearts. Hum. Mutat., 24, 104.[Medline]
-
Reamon-Buettner, S.M. and Borlak, J. (2005) GATA4 zinc finger mutations as a molecular rationale for septation defects of the human heart. J. Med. Genet., 42:e32
[Free Full Text] -
Colgin, L.M., Hackmann, A.F., Emond, M.J. and Monnat, R.J., Jr. (2002) The unexpected landscape of in vivo somatic mutation in a human epithelial cell lineage. Proc. Natl Acad. Sci. USA, 99, 14371442.
[Abstract/Free Full Text] -
Inga, A., Storici, F., Darden, T.A. and Resnick, M.A. (2002) Differential transactivation by the p53 transcription factor is highly dependent on p53 level and promoter target sequence. Mol. Cell. Biol., 22, 86128625.
[Abstract/Free Full Text] - Kennedy, B.K. (2002) Mammalian transcription factors in yeast: strangers in a familiar land. Nat. Rev. Mol. Cell Biol., 3, 4149.[CrossRef][Web of Science][Medline]
- Inga, A., Storici, F. and Resnick, M.A. (2005) In Heitman, J. and Nitiss, J.L. (eds), Kluwer Academic Publishers, in press.
- Resnick, M.A. and Cox, B.S. (2000) Yeast as an honorary mammal. Mutat. Res., 451, 111.[Web of Science][Medline]
-
Storici, F., Durham, C.L., Gordenin, D.A. and Resnick, M.A. (2003) Chromosomal site-specific double-strand breaks are efficiently targeted for repair by oligonucleotides in yeast. Proc. Natl Acad. Sci. USA, 100, 1499414999.
[Abstract/Free Full Text] -
Resnick, M.A., and Inga, A. (2003) Functional mutants of the sequence-specific transcription factor p53 and implications for master genes of diversity. Proc. Natl Acad. Sci. USA, 100, 99349939.
[Abstract/Free Full Text] - Harvey, R.P. (1996) NK-2 homeobox genes and heart development. Dev. Biol., 178, 203216.[CrossRef][Web of Science][Medline]
-
Kasahara, H., Usheva, A., Ueyama, T., Aoki, H., Horikoshi, N. and Izumo, S. (2001) Characterization of homo- and heterodimerization of cardiac Csx/Nkx2.5 homeoprotein. J. Biol. Chem., 276, 45704580.
[Abstract/Free Full Text] -
Chen, C.Y. and Schwartz, R.J. (1995) Identification of novel DNA binding targets and regulatory domains of a murine tinman homeodomain factor, nkx-2.5. J. Biol. Chem., 270, 1562815633.
[Abstract/Free Full Text] - Garg, V., Kathiriya, I.S., Barnes, R., Schluterman, M.K., King, I.N., Butler, C.A., Rothrock, C.R., Eapen, R.S., Hirayama-Yamada, K., Joo, K. et al. (2003) GATA4 mutations cause human congenital heart defects and reveal an interaction with TBX5. Nature, 424, 443447.[CrossRef][Medline]
-
Shiojima, I., Komuro, I., Oka, T., Hiroi, Y., Mizuno, T., Takimoto, E., Monzen, K., Aikawa, R., Akazawa, H., Yamazaki, T., Kudoh, S. and Yazaki, Y. (1999) Context-dependent transcriptional cooperation mediated by cardiac transcription factors Csx/Nkx-2.5 and GATA-4. J. Biol. Chem., 274, 82318239.
[Abstract/Free Full Text] -
Schott, J.J., Benson, D.W., Basson, C.T., Pease, W., Silberbach, G.M., Moak, J.P., Maron, B.J., Seidman, C.E. and Seidman, J.G. (1998) Congenital heart disease caused by mutations in the transcription factor NKX2-5. Science, 281, 108111.
[Abstract/Free Full Text] -
Elliott, D.A., Kirk, E.P., Yeoh, T., Chandar, S., McKenzie, F., Taylor, P., Grossfeld, P., Fatkin, D., Jones, O., Hayes, P., Feneley, M. and Harvey, R.P. (2003) Cardiac homeobox gene NKX2-5 mutations and congenital heart disease: associations with atrial septal defect and hypoplastic left heart syndrome. J. Am. Coll. Cardiol., 41, 20722076.
[Abstract/Free Full Text] -
Goldmuntz, E., Geiger, E. and Benson, D.W. (2001) NKX2.5 mutations in patients with tetralogy of fallot. Circulation, 104, 25652568.
[Abstract/Free Full Text] -
McElhinney, D.B., Geiger, E., Blinder, J., Woodrow, B.D. and Goldmuntz, E. (2003) NKX2.5 mutations in patients with congenital heart disease. J. Am. Coll. Cardiol., 42, 16501655.
[Abstract/Free Full Text] - Hobbs, C.A., Cleves, M.A., Keith, C., Ghaffar, S. and James, S.J. (2005) NKX2.5 and congenital heart defects: A population-based study. Am. J. Med. Genet. A, 134, 223225.[Medline]
-
Gruber, P.J., Epstein, J.A. (2004) Development gone awry: congenital heart disease. Circ. Res., 94, 273283.
[Abstract/Free Full Text] -
Heathcote, K., Braybrook, C., Abushaban, L., Guy, M., Khetyar, M.E., Patton, M.A., Carter, N.D., Scambler, P.J. and Syrris, P. (2005) Common arterial trunk associated with a homeodomain mutation of NKX2.6. Hum. Mol. Genet., 14, 585593.
[Abstract/Free Full Text] - Durocher, D., Charron, F., Warren, R., Schwartz, R.J. and Nemer, M. (1997) The cardiac transcription factors Nkx2-5 and GATA-4 are mutual cofactors. EMBO J., 16, 56875696.[CrossRef][Web of Science][Medline]
-
Lee, Y., Shioi, T., Kasahara, H., Jobe, S.M., Wiese, R.J., Markham, B.E. and Izumo, S. (1998) The cardiac tissue-restricted homeobox protein Csx/Nkx2.5 physically associates with the zinc finger protein GATA4 and cooperatively activates atrial natriuretic factor gene expression. Mol. Cell. Biol., 18, 31203129.
[Abstract/Free Full Text] -
Sepulveda, J.L., Belaguli, N., Nigam, V., Chen, C.Y., Nemer, M. and Schwartz, R.J. (1998) GATA-4 and Nkx-2.5 coactivate Nkx-2 DNA binding targets: role for regulating early cardiac gene expression. Mol. Cell. Biol., 18, 34053415.
[Abstract/Free Full Text] -
Okubo, A., Miyoshi, O., Baba, K., Takagi, M., Tsukamoto, K., Kinoshita, A., Yoshiura, K., Kishino, T., Ohta, T., Niikawa, N. and Matsumoto, N. (2004) A novel GATA4 mutation completely segregated with atrial septal defect in a large Japanese family. J. Med. Genet., 41, e97.
[Free Full Text] -
Flaman, J.M., Frebourg, T., Moreau, V., Charbonnier, F., Martin, C., Chappuis, P., Sappino, A.P., Limacher, I.M., Bron, L., Benhattar, J. Tada, M., van Meir, E.G., Estricher, A., and Iggo, R.D (1995) A simple p53 functional assay for screening cell lines, blood, and tumors. Proc. Natl Acad. Sci. USA, 92, 39633967.
[Abstract/Free Full Text] -
deVere White, R.W., Deitch, A.D., Gumerlock, P.H. and Shi, X.B. (1999) Use of a yeast assay to detect functional alterations in p53 in prostate cancer: review and future directions. Prostate, 41, 134142.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
S. M. Reamon-Buettner, Y. Ciribilli, I. Traverso, B. Kuhls, A. Inga, and J. Borlak A functional genetic study identifies HAND1 mutations in septation defects of the human heart Hum. Mol. Genet., October 1, 2009; 18(19): 3567 - 3578. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Reamon-Buettner, Y. Ciribilli, A. Inga, and J. Borlak A loss-of-function mutation in the binding domain of HAND1 predicts hypoplasia of the human hearts Hum. Mol. Genet., May 15, 2008; 17(10): 1397 - 1405. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. G. Jegga, A. Inga, D. Menendez, B. J. Aronow, and M. A. Resnick Functional evolution of the p53 regulatory network through its target response elements PNAS, January 22, 2008; 105(3): 944 - 949. [Abstract] [Full Text] [PDF] |
||||
![]() |
S M Reamon-Buettner and J Borlak Somatic mutations in cardiac malformations. J. Med. Genet., August 1, 2006; 43(8): e45 - e45. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





