- Research article
- Open Access
Identification and analysis of KLF13 variants in patients with congenital heart disease
BMC Medical Genetics volume 21, Article number: 78 (2020)
The protein Kruppel-like factor 13 (KLF13) is a member of the KLF family and has been identified as a cardiac transcription factor that is involved in heart development. However, the relationship between KLF13 variants and CHDs in humans remains largely unknown. The present study aimed to screen the KLF13 variants in CHD patients and genetically analyze the functions of these variants.
KLF13 variants were sequenced in a cohort of 309 CHD patients and population-matched healthy controls (n = 200) using targeted sequencing. To investigate the effect of variants on the functional properties of the KLF13 protein, the expression and subcellular localization of the protein, as well as the transcriptional activities of downstream genes and physical interactions with other transcription factors, were assessed.
Two heterozygous variants, c.487C > T (P163S) and c.467G > A (S156N), were identified in two out of 309 CHD patients with tricuspid valve atresia and transposition of the great arteries, respectively. No variants were found among healthy controls. The variant c.467G > A (S156N) had increased protein expression and enhanced functionality compared with the wild type, without affecting the subcellular localization. The other variant, c.487C > T (P163S), did not show any abnormalities in protein expression or subcellular localization; however, it inhibited the transcriptional activities of downstream target genes and physically interacted with TBX5, another cardiac transcription factor.
Our results show that the S156N and P163S variants may affect the transcriptional function of KLF13 and physical interaction with TBX5. These results identified KLF13 as a potential genetic risk factor for congenital heart disease.
Congenital heart disease (CHD) is a term that refers to any anatomical defect in the heart or major blood vessels that is present at birth. CHD affects 1% of newborn babies, and approximately one to two per 1000 newborn babies have critical CHDs that can cause death during the neonatal period, such as pulmonary atresia and transposition of the great arteries [1,2,3]. Infants and children with CHDs often suffer from delays in growth and serious progressive cardiac complications . Despite significant advances in our understanding of cardiac morphogenesis, emerging evidence suggests that genetic risk factors play a crucial role in the etiology of CHD. Further elucidation of the molecular and genetic intricacies of CHDs will likely be an important step in improving treatments for CHD patients.
The heart is the first functional organ to develop during vertebrate embryogenesis, which is a complex process in which multiple cell lineages contribute to organ development. Epithelial–mesenchymal transformation (EMT) is an important step during heart development and is regulated by myocardial signals to generate mesenchymal cells that migrate into the cardiac jelly and proliferate to cellularize the endocardial cushions [5, 6]. Malformation of the atrioventricular canal cushions can lead to atrioventricular septal defects (AVSDs) or other endocardial cushion defects. Kruppel-like factor 13 (KLF13) protein has been identified as a novel collaborator of the cardiac GATA4 and TBX5 factors, which are highly expressed in the endocardium and myocardium during embryonic development. Knocking down KLF13 alleles in Xenopus embryos leads to valve immaturity, septal defects, and hypotrabeculation [7,8,9]. However, the mechanisms of KLF13 dysfunction causing CHDs are still poorly understood.
A spatiotemporal expression analysis of KLF13 in murine hearts showed that the earliest sign of expression was at E9.5, and subsequent expression was detected in the developing atrial myocardium, ventricular trabeculae, AV cushions, and truncus arteriosus, which then decreased in the postnatal heart . Additionally, in Xenopus embryos, the highest levels of KLF13 were found in the interventricular septum and the atrioventricular valves, and knocking down KLF13 induced ventricular hypotrabeculation, atrial septal defects (ASDs) and delayed AV cushion formation [ ].
KLF13 is a member of the KLF family of transcription factors that have been categorized along with the specificity protein (Sp) family of transcription factors . To date, twenty-one KLF members have been identified in humans. These transcription factors are characterized by C-terminal Cys2His2 (C2H2) zinc-finger motifs that confer preferential binding to GC/GT rich sequences in gene promoters and enhancer regions to activate or inhibit downstream target expression [11, 12]. ANF and BNP are natriuretic peptides that are expressed in the heart and are developmentally regulated. Their levels rise continuously during embryonic cardiac development in both atria and ventricles when cells differentiate into cardiomyocytes. Their levels then drop after postnatal development. At the genetic level, analysis of the ANP and BNP promoters has led to the characterization of key cardiac transcription factors that govern cardiac growth and differentiation . Studies by Lavallee et al. revealed that KLF13 could activate the BNP promoter [ ].
Additionally, analysis of the binding interactions of transcription factor variants helps us better understand how the disruption of combinatorial interactions can lead to specific congenital disabilities. TBX5 is a member of the T-box transcription factor family which has crucial roles in regulating early cellular commitment, differentiation, and heart development . TBX5 cooperates with other transcription factors, such as GATA4 and NKX2.5, to synergistically regulate downstream targets during cardiac development [ [15, 16]]. Previously, we showed that KLF13 is a TBX5 cofactor and that KLF13 may be a genetic modifier of Holt-Oram Syndrome and possibly other congenital heart diseases .
However, since current studies of KLF13 have focused on animal models [ [7, 17]], it remains unclear whether genetic variants are involved in the mechanisms of CHDs in humans. In the present study, we identified two KLF13 heterozygous variants in a cohort of patients with complex CHDs and compared them with those of healthy controls to evaluate the prevalence of KLF13 variants in sporadic CHDs. Our results demonstrated that these variants altered protein expression, changed the transcriptional activation of BNP and impaired the genetic interaction of KLF13 with TBX5.
In this study, we recruited a total of 309 patients with complex CHDs. These patients were diagnosed by echocardiography or cardiac catheterization or underwent cardiac surgery at Shanghai Xinhua Hospital. The patients included 191 males and 118 females (Table 1). Patients with known syndromic CHDs or chromosomal abnormalities were excluded from our study. The controls were 200 population-matched healthy children without heart disease. The study protocol was reviewed and approved by the Xinhua Hospital Ethics Committee. Both parents and legal guardians of the patients and healthy controls provided signed informed consent. Subsequently, peripheral blood was collected for DNA extraction. Genomic DNA was extracted using the QIAamp DNA Blood Mini Kit (Qiagen, Germany) and stored at − 80 °C.
Target sequencing and analysis
Genomic DNA was sequenced by target sequencing technology using the Illumina HiSeq 2000 platform for variants of KLF13 (GenBank accession number NC_000015.10, NM_015995.3) and several other cardiac transcriptional factors involved in cardiovascular development (GATA4, TBX5, TBX1, GATA6, GATA5 and so on). Then, Sanger sequencing was performed to validate all the candidate variants. To evaluate the protein characteristics of nonsynonymous variants, we used SIFT (http://sift.jcvi.org/), PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/), and Mutation Taster (www.mutationtaster.org/). Amino acid substitutions were predicted as damaging when the score was ≤0.05 in SIFT or ≥ 0.85 in PolyPhen-2. KLF13 protein sequences from Homo sapiens (human), Mus musculus (house mouse), Bos taurus (cattle), Capra hircus (goat), Pan troglodytes (chimpanzee), Xenopus laevis (frog) and Sus scrofa (swine) were downloaded from the Universal Protein (UniProt) database (http://uniprot.org/) and aligned with Clustal X software.
Plasmid construction and site-directed mutagenesis
The KLF13 and TBX5 cDNA plasmids were purchased from Genomeditech. Site-directed mutagenesis for the KLF13 point mutations c.467G > A (S156N) and c.487C > T (P163S) was performed according to the protocol provided for the Site-Directed Mutagenesis Kit (Stratagene, USA). Then, the mutated sites were confirmed with Sanger sequencing. The luciferase reporter with human B-type natriuretic peptide (BNP) promoter was constructed as previously described [ ].
Cell cultures and transfection
For cell culture experiments, 293 T cells were used for protein extraction and immunofluorescent staining. NIH 3 T3 cells were used for luciferase assays and were maintained in growth medium (Dulbecco’s modified Eagle’s medium) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin. Plasmids were transfected with FuGene HD (Promega, USA) according to the manufacturer’s protocol 24 h after the cells were seeded.
NIH 3 T3 cells were seeded onto a 24-well plate, and 600 ng of BNP luciferase reporter vector, different dosages (25, 50, 100, 200, and 300 ng) of the wild-type/variant KLF13 plasmid, and pCMV-Tag2B vector with or without TBX5 plasmids were cotransfected with FuGene HD. The luciferase activity was measured by the Dual-Glo Luciferase Assay System (Promega) following the manufacturer’s protocol after 48 h of transfection. Firefly luciferase activities are reported as fold activation. All experiments were repeated at least three times in duplicate.
For myc staining, 293 T cells were harvested 48 h after transfection with wild-type KLF13 or its variants. After fixation with 4% paraformaldehyde at room temperature for 10 min, permeabilization with 0.3% Triton X-100 for 15 min, and blocking with 5% BSA for 1 h, cells were incubated with the primary antibody (1:500, rabbit polyclonal, anti-myc, 16,286–1-AP, Proteintech, USA) overnight at a temperature of 4 °C, followed by incubation with the conjugated secondary antibody, anti-rabbit Alexa Fluor® 488 (1:500, R37116, Life Technologies, USA) at room temperature in the dark for 1 h. DAPI (Vector Laboratories, USA) was used for nuclear staining. The images were acquired using fluorescence microscopy at 200X magnification with Image-Pro Plus software.
Western blotting analysis
293 T cells were harvested at 48 h after transfection with wild-type KLF13 or its variants. The protein was extracted as previously described . A total of 50 μg of the extract was subjected to 10% SDS-PAGE, transferred onto nitrocellulose membranes and blocked with skim milk (5%) in TBST (0.1%) at room temperature for 2 h, using gentle agitation during the process. Then, the membranes were incubated with the primary antibody (myc: 1:1000, rabbit polyclonal anti-myc, 16,286–1-AP, Proteintech, USA; GAPDH: 1:1000, rabbit monoclonal anti-GAPDH, ab181602, Abcam, USA) overnight at a temperature of 4 °C. Following incubation with the primary antibody, horseradish peroxidase-conjugated secondary antibody (1:1000, 111–005-003, Jackson ImmunoResearch, USA) and Immobilon Western Chemiluminescent HRP Substrate (Millipore, USA) were used to visualize chemiluminescent immunodetection.
293 T cells were transfected with wild-type KLF13 or variants and TBX5 plasmids using FuGene HD according to the manufacturer’s guidelines. Protein extracts were incubated overnight with anti-myc antibody (1:100, rabbit polyclonal anti-myc, 16,286–1-AP, Proteintech, USA) coupled with magnetic beads. Bound proteins were revealed with anti-TBX5 (1:1000, rabbit polyclonal anti-TBX5, 42–6500, Invitrogen, USA) antibody by Western blotting.
Data are reported as the means ± SEMs. Two-way analysis of variance (ANOVA) was used to compare differences among groups by SPSS. A P value < 0.05 was selected to indicate statistical significance.
KLF13 sequencing in CHD patients
In this study, a cohort of 309 patients with CHDs was established. Additionally, a total of 200 ethnically matched healthy subjects were enrolled and used as controls. The baseline clinical characteristics of the 309 CHD patients are summarized in Table 1. The entire coding region of the KLF13 gene was sequenced in all patients with CHDs by targeted sequencing. We identified two heterozygous variants of KLF13 in two out of 309 patients. A transversion of proline into serine at amino acid position 163 (exon 1 c.487C > T P163S) was identified in one patient with TA VSD ASD. Additionally, a transversion of serine into asparagine at the first nucleotide of codon 156 of the KLF13 gene (exon 1: c.467G > A S156N) was found in the D-TGA patient (Fig. 1a-d). These two patients did not have any variants in other sequenced genes, and no variants of any kind were found among the control individuals. SIFT, PolyPhen-2, and Mutation Taster were used to perform bioinformatic predictions for these two variants. Details of the specific sequence variants are summarized in Table 2.
In silico analyses of KLF13 mutant protein
Human KLF13 has 6825 bp on chromosome 15, is composed of two exons and two introns and encodes 288 amino acids. Song et al. analyzed the KLF13 protein structure and revealed that it is composed of three zinc-finger domains (Zns), distinct transcriptional activation and repression domains, and two potent independent nuclear localization signals (NLS) . Both S156N and P163S were located in the NLS1 region that was close to the DNA-binding domain (Fig. 1e). Multiple sequence alignments revealed that the S156N and P163S variants were highly conserved in vertebrates (Fig. 1f), indicating that these variants are important and might result in alteration of gene function.
Subcellular localization and expression of KLF13 variants
We analyzed the protein expression of these two variants in the nucleus and cytoplasm. The subcellular localization was determined via immunofluorescence (Fig. 2a). Although both mutations were located in the NLS1 region, the variants had an equal protein distribution in the two cellular compartments compared to the wild type. However, measurements of the expression of the myc-tagged-KLF13 variant using total protein extracts prepared from transfected 293 T cells (Fig. 2b) showed that the expression of the S156N variant was distinctly higher than that of wild-type KLF13 (P < 0.05), while the P163S variant had a similar expression to that of the wild type, which suggested that the S156N variant might be a gain-of-function variant.
Transcriptional activity of KLF13 variant protein
To evaluate the functional ability of the KLF13 variants to regulate downstream genes, the transcriptional activities of the variants were analyzed using a BNP-luciferase reporter, and cells were cotransfected with wild type or KLF13 variants. As shown in Fig. 2c, a dose response study (25, 50, 100, 200, and 300 ng) showed that, in comparison to wild-type KLF13, the P163S variant had decreased transcriptional activity. In contrast, the S156N variant had increased transcriptional activity, which indicated that the KLF13 variants affected BNP promoter transactivation (P < 0.001). Therefore, the variants S156N and P163S may affect KLF13 downstream target gene expression during cardiogenesis, for which variant S156N was a gain-of-function mutant, while variant P163S disrupted the function of KLF13.
KLF13 variants affect the physical interactions with TBX5
Previously, we reported that KLF13 is a TBX5 cofactor; these two proteins were found to colocalize in several cardiac cells, and they physically and functionally interact to synergistically activate cardiac promoters . Therefore, to verify the impact of KLF13 variants on physical interactions with TBX5, wild-type or variant KLF13 constructs were cotransfected with TBX5 using the BNP-luciferase reporter. As shown in Fig. 3a, when NIH 3 T3 cells were cotransfected with KLF13 and TBX5, the BNP-luciferase reporter was activated up to 75-fold. Decreased synergistic transcriptional activation was detected with variant P163S, and increased synergy was observed with variant S156N. To test whether the altered synergy was due to altered physical interactions, the binding of wild-type KLF13 or KLF13 variants with TBX5 was tested by co-IP. Figure 3b shows the analysis of the co-IP assay data. P163S showed reduced synergy with TBX5 and had a decreased physical interaction with the TBX5 protein. In contrast, S156N had a significantly increased ability to physically and functionally interact with TBX5.
KLF13, a member of the KLF family of transcription factors, has been identified as a cardiac factor that is involved in heart development [7, 17, 19, 20]. Previous research by Lavallee et al. showed that KLF13 is expressed during cardiac morphogenesis in mice and that knocking it down induced ventricular hypotrabeculation, atrial septal defects (ASDs), and delayed AV cushion formation in Xenopus embryos . However, variants of KLF13 in CHD patients have remained unidentified and uncharacterized. In the present study, two heterozygous KLF13 variants were discovered using targeted sequencing in 309 patients with complex CHDs. Variant c.467G > A (S156N), identified in patients with TA, VSD, and ASD, had altered protein expression and functionality compared with the wild-type protein, without affecting the subcellular localization. The other variant, c.487C > T (P163S), found in a D-TGA patient, did not show any abnormalities in protein expression or subcellular localization. However, this variant had altered transcriptional activities of downstream genes and altered physical interactions with another cardiac factor, TBX5.
Both variants, c.467G > A and c.487C > T, were highly preserved in vertebrates based on multiple sequence alignment analyses and were located in the NLS1 region. Song et al. identified a bipartite NLS in the KLF13 protein based on sequence searching and found that AA 147–168 (NLS1) was not the only NLS region for protein subcellular localization because deletion of this signal did not affect nuclear transport . In our present results, although the changed sequences of these two variants were located in the NLS1 region, the variants were also expressed in the nucleus, as was the wild-type KLF13 protein. However, the locations of the mutations in these two variants were close to the Cys2/His2 zinc-finger motifs that bind to gene promoters and enhancer regions to activate or inhibit transcription , suggesting that these two variants might have crucial functions in mediating the expression of downstream target genes. Lavallee et al. used the BNP promoter and showed that KLF13 could activate the promoter by binding to an evolutionarily conserved CACCC box . Although the subcellular localization of the variant protein did not change, which is in contrast to most previously reported CHD-associated variants [22, 23], we found that the S156N variant was a gain-of-function mutation. The mutation at this site may promote DNA binding activity by increasing the protein levels, which may explain the upregulated BNP in this study. Although a gain-of-function mutation was detected in our experiments, the protein that might be less stable in vivo could also result in a loss-of-function phenotype [24, 25]. Consistent with this observation, human heterozygous microdeletions and duplication of the chromosomal band harboring KLF13 (15q13.3) have been demonstrated to be associated with a wide range of cardiac defects . Unlike the S156N variant, the c.487C > T variant decreased the transcriptional activity of BNP in the present study, but we did not find any disruption of expression or subcellular localization. However, a proline-to-serine substitution may result in impaired interactions with DNA, and the location of the P163S variant in the tertiary structure that is near the bound DNA might also lead to a decrease in the DNA-binding activity.
The reduced binding ability of our identified variants to their targets may be primarily related to decreased DNA binding. However, structural instability of the variants when bound to cofactors is another possibility . The mechanism of the function of cardiac transcription factors in the developing heart involves cooperative interactions with other conserved factors, and any structural defect is often linked to more than one candidate gene. Additionally, variable expressivity of the phenotype is often observed in cases where CHDs are linked to a specific gene mutation . Thus, variants of these factors may eliminate any cooperative interactions between these cofactors [ [28, 29]]. Previously, we showed that KLF13 and TBX5 exhibit physical and functional interactions and that a combined heterozygous loss of TBX5 and its associated KLF13 leads to decreased postnatal viability and a higher incidence of septal defects compared to that of Tbx5 heterozygous mice. Additionally, several TBX5 mutations have been associated with CHDs, as indicated by impaired functional and physical interactions with KLF13, which confirms that TBX5 is a genetic cofactor of KLF13 . In the present study, the increased transcriptional activity of BNP was more pronounced when wild-type KLF13 was cotransfected with TBX5, which indicated that KLF13 has functional interactions with TBX5 to activate BNP promoters. In our previous study of the structure-function of TBX5-KLF13 synergy , the N-terminal domain of KLF13 was mainly found to support synergy with TBX5. The two variants of KLF13 that were detected in the 309 CHD patients were in the N-terminal domain of the protein. However, the two variants had different effects on the interactions with TBX5. Significantly enhanced synergistic activity was detected in the S156N variant, while the P163S variant eliminated this synergistic activity, which suggested that the loss of function observed in KLF13 was associated with disruption of the physical and functional interactions with cofactors (TBX5), while the gain of function could enhance the interactions with cofactors.
This study has some limitations. Corrections between genotype and phenotype remain challenging in exploring CHDs. In this study, although the variant P163S eliminated synergistic activity with TBX5, the patient did not have typical Holt-Oram Syndrome or other TBX5-related features. It is uncertain what other consequences have not been anticipated, such as other transcription factors not investigated in this study that interacted with KLF13. Thus, larger cohort retrospective and prospective studies of CHDs are still needed to pave the way for future elucidation of its etiology.
This study suggested that CHD patients harbored KLF13 gene variants, and KLF13 variants may contribute to the occurrence of congenital heart disease. Additional studies that test a larger population will help further our understanding of the role of KLF13 in CHD molecular pathogenesis.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Congenital heart disease
Kruppel-like factor 13
Atrioventricular septal defects
Atrial septal defects
Ventricular septal defects
Tricuspid valve atresia transposition of the great arteries,
Transposition of the great arteries
Nuclear localization signals
Zinc fingers domain
vander Linde D, Konings EE, Slager MA, Witsenburg M, Helbing WA, Takkenberg JJ, Roos-Hesselink JW. Birth prevalence of congenital heart disease worldwide: a systematic review and meta-analysis. J Am Coll Cardiol. 2011;58:2241–7.
Zhao QM, Ma XJ, Jia B, Huang GY. Prevalence of congenital heart disease at live birth: anaccurate assessment by echocardiographic screening. Acta Paediatr. 2013;102:397–402.
Chang RK, Gurvitz M, Rodriguez S. Missed diagnosis of critical congenital heart disease. Arch Pediatr Adolesc Med. 2008;162:969–74.
Cameron JW, Rosenthal A, Olson AD. Malnutrittionin hospitalized children with congenital heart disease. Arch Pediatr Adolesc Med. 1995;149(10):1098–102.
Eisenberg LM, Markwald RR. Molecular regulation of atrioventricular valvuloseptal morphogenesis. Circ Res. 1995;77:1–6.
Anderson RH, Webb S, Brown NA, Lamers W, Moorman A. Development of the heart: (3) formation of the ventricular outflow tracts, arterial valves, and intrapericardial arterial trunks. Heart. 2003b;89:1110–8.
Lavallée G, Andelfinger G, Nadeau M, Lefebvre C, Nemer G, Horb ME, et al. The Kruppel-like transcription factor KLF13 is a novel regulator of heart development. EMBO J. 2006;25(21):5201–13.
Nemer M, Horbi ME. The KLF of family of transcriptional regulators in cardiomyocyte proliferation and differentiation. Cell Cycle. 2007;6(2):117–21.
Darwich R, Li W, Yamak A, Komati H, Andelfinger G, Sun K, et al. KLF13 is a genetic modifier of the Holt-Oram syndrome gene TBX5. Hum Mol Genet. 2017;26(5):942–54.
Lomberk G, Urrutia R. The family feud: turning off Sp1 by Sp1-like KLF proteins. Biochem J. 2005;392(1):1–11.
Song A, Patel A, Thamatrakoln K, Liu C, Feng D, Clayberger C, et al. Functional domains and DNA-binding sequences of RFLAT-1/KLF13, a Krüppel-like transcription factor of activated T lymphocytes. J Biol Chem. 2002;277(33):30055–65.
Pearson R, Fleetwood J, Eaton S, Crossley M, Bao S. Krüppel-like transcription factors: a functional family. Int J Biochem Cell Biol. 2008;40(10):1996–2001.
Hayek S, Nemer M. Cardiac natriuretic peptides: from basic discovery to clinical practice. Cardiovasc Ther. 2011;29(6):362–76.
Hariri F, Nemer M, Nemer G. T-box factors: insights into the evolutionary emergence of the complex heart. Ann Med. 2012;44:680–93.
Garg V, Kathiriya IS, Barnes R, Schluterman MK, King IN, Butler CA, et al. GATA4 mutations cause human congenital heart defects and reveal an interaction With TBX5. Nature. 2003;424:443–7.
Hiroi Y, Kudoh S, Monzen K, Ikeda Y, Yazaki Y, Nagai R, et al. Tbx5 associates with Nkx2–5 and synergistically promotes cardiomyocyte differentiation. Nat Genet. 2001;28:276–80.
Martin KM, Metcalfe JC, Kemp PR. Expression of Klf9 and Klf13 in mouse development. Mech Dev. 2001;103(1–2):149–51.
Zhang E, Hong N, Chen S, Fu Q, Li F, Yu Y, Sun K. Target sequencing identifies novel GATA6 variants in a large cohort of patients with conotruncal heart disease. Gene. 2018;641:341–8.
Chen CP, Chen CY, Chern SR, Wu PS, Chen SW, Wu FT, et al. Detection of de novo del (18)(q22.2) and a familial of 15q13.2-q13.3 microduplication in a fetus with congenital heart defects. Taiwan J Obstet Gynecol. 2019;58(5):704–8.
Derwińska K, Bartnik M, Wiśniowiecka-Kowalnik B, Jagła M, Rudziński A, Pietrzyk JJ, et al. Assessment of the role of copy-number variants in 150 patients with congenital heart defects. Med Wieku Rozwoj. 2012;16(3):175–82.
Toshio Nagashima FH, Umehara T, Yokoyama AS. In: Nagai RF, Scott L, Kasuga M, editors. Molecular Structures of Krüppel-like Factors, in The Biology of Krüppel-like Factors, vol. 268. Tokyo: Springer; 2009.
Huang RT, Wang J, Xue S, Qiu XB, Shi HY, Li RG, et al. TBX20 loss-of-function mutation responsible for familial tetralogy of Fallot or sporadic persistent truncus arteriosus. Int J Med Sci. 2017;14(4):323–32.
Yang YQ, Gharibeh L, Li RG, Xin YF, Wang J, Liu ZM, et al. GATA4 loss-of-function mutations underlie familial tetralogy of fallot. Hum Mutat. 2013;34(12):1662–71.
Posch MG, Gramlich M, Sunde M, Schmitt KR, Lee SH, Richter S, et al. A gain-of function TBX20 mutation causes congenital atrial septal defects, patent foramen ovaleand cardiac valve defects. J Med Genet. 2010;47(4):230–5.
Bouveret R, Waardenberg AJ, Schonrock N, Ramialison M, Doan T, de Jong D, et al. NKX2-5 mutations causative for congenital heart disease retain functionality and are directed to hundreds of target. Elife. 2015;6:4.
Van Bon BW, Mefford HC, Menten B, Koolen DA, Sharp AJ, Nillesen WM, et al. Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical sepectrumvarying from non-pathogenic to a severe outcome. J Med Genet. 2009;46(8):511–23.
Kasinath V, Faini M, Poepsel S, Reif D, Feng XA, Stjepanovic G, et al. Structures of human PRC2 with its cofactors AEBP2 and JARID2. Science. 2018;359(6378):940–4.
Nemer M. Genetic insights into normal and abnormal heart development. Cardiovasc Pathol. 2008;17(1):48–54.
Temsah R, Nemer M. GATA factors and transcriptional regulation of cardiac natriuretic peptide genes. Regul Pept. 2005;3:177–85.
The authors think Pro. Rang Xu for advice and LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.
This works was supported by two grants from the National Basic Research Program of China, the “Study on the molecular mechanisms of genetics interaction between GATA5 and KLF13 in aortic valvulogenesis and bicuspid aortic valve pathogenesis” (grant number 81720108003), and the “Study on the mechanisms of PEX1 inhibition of DOX-induced cardiomyocytes apoptosis” (grant number 81700226).
Ethics approval and consent to participate
Research was approved by the Xinhua Hospital Ethics Committee affiliated to Shanghai Jiaotong University School of Medicine (XHEC-C-2017-105). All parents were informed of the purpose and significance of the experiment and signed an informed consent form.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, W., Li, B., Li, T. et al. Identification and analysis of KLF13 variants in patients with congenital heart disease. BMC Med Genet 21, 78 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-020-01009-x