- Open Access
A loss-of-function mutation p.T52S in RIPPLY3 is a potential predisposing genetic risk factor for Chinese Han conotruncal heart defect patients without the 22q11.2 deletion/duplication
Journal of Translational Medicinevolume 16, Article number: 260 (2018)
Conotruncal heart defect (CTD) is a complex congenital heart disease with a complex and poorly understood etiology. The transcriptional corepressor RIPPLY3 plays a pivotal role in heart development as a negative regulator of the key cardiac transcription factor TBX1. A previous study showed that RIPPLY3 contribute to cardiac outflow tract development in mice, however, the relationship between RIPPLY3 and human cardiac malformation has not been reported.
615 unrelated CTD Chinese Han patients were enrolled, we excluded the 22q11.2 deletion/duplication using a modified multiplex ligation-dependent probe amplification method—CNVplex®, and investigated the variants of RIPPLY3 in 577 patients without the 22q11.2 deletion/duplication by target sequencing. Functional assays were performed to testify the potential pathogenicity of nonsynonymous variants found in these CTD patients.
Four rare heterozygous nonsynonymous variants (p.P30L, p.T52S, p.D113N and p.V179D) were identified in four CTD patients, the variant NM_018962.2:c.155C>G (p.T52S) is referred as rs745539198, and the variant NM_018962.2:c.337G>A (p.D113N) is referred as rs747419773. However, variants p.P30L and p.V179D were not found in multiple online human gene variation databases. Western blot analysis and immunofluorescence showed that there were no significant difference between wild type RIPPLY3 and these four variants. Luciferase assays revealed that the p.T52S variant altered the inhibition of TBX1 transcriptional activity in vitro, and co-immunoprecipitation assays showed that the p.T52S variant reduced the physical interaction of RIPPLY3 with TBX1. In addition to the results from pathogenicity prediction tools and evolutionary protein conservation, the p.T52S variant was thought to be a potentially deleterious variant.
Our results provide evidence that deleterious variants in RIPPLY3 are potential molecular mechanisms involved in the pathogenesis of human CTD.
Congenital heart disease (CHD), including heterogeneous anatomy with distinct phenotypic subtypes, is the most common form of birth defect in humans, and it affects approximately 1% of all live births [1, 2]. Among these, conotruncal heart defect (CTD) is a severe malformation. It is characterized by a disordered orchestration of the ventricles, the aorta and the pulmonary artery, and consists of the tetralogy of Fallot (TOF), transposition of the great arteries (TGA), double outlet of right ventricle (DORV), persistent truncus arteriosus (PTA), pulmonary atresia with ventricular septal defect (PA/VSD), and interrupted aortic arch (IAA). CTD is estimated to occur in approximately 1 out of 1000 live births . It is the most common cyanosis CHD, and it usually requires surgical treatment. Patients born with CTD often need lifelong specialized cardiac care.
Although our understanding of the molecular pathways underlying cardiac development has grown remarkably in the past few years, most genetic factors that predispose patients to CTD remain unidentified . Most of the known causes of CTD are sporadic genetic changes such as nucleotide variants or focal copy-number DNA variations . 22q11.2 deletion syndrome is the most common microdeletion syndrome, affecting nearly 1 in 2000–4000 live births . It has been reported to be associated with various types of congenital heart defects, especially CTD . TBX1 (MIM 602054), a member of the phylogenetically conserved T-box gene family of DNA-binding transcription factors, was mapped to the critical 1.5 Mb deletion region of 22q11.2 . TBX1 variants were reported to be responsible for the cardiac phenotype of 22q11.2 deletion syndrome, and it plays a vital role in the development of the cardiac outflow tract [7,8,9,10,11]. However, copy-number or single-nucleotide variants of TBX1 account for only approximately 7% of CTD patients in China [11, 12].
RIPPLY proteins (including RIPPLY1/2/3) have been shown to modulate the transcriptional properties of T-box proteins, which play a pivotal role in embryonic development [13,14,15,16]. They interact with the transcriptional corepressor Groucho/TLE and the T-box proteins through two distinct amino acid sequences: the WRPW motif, which is a highly conserved four-amino-acid stretch in the N-terminal half, and the Ripply homology (RH) domain, which is a conserved ~ 50-amino-acid stretch that interacts with the T-domain [15, 16]. In mouse embryo, RIPPLY1 mRNA was first detected at E8.5, and was continuously expressed until E12.5. During this period, RIPPLY1 was mainly expressed in the anterior presomitic mesoderm (PSM) . Similarly, RIPPLY2 expression was localized to the PSM in wild-type mouse embryos from E9.0 to E11.5, when somites were being formed . Moreover, RIPPLY2-null mice had a partial pedicle of neural arches and separated lamina, RIPPLY1/2 double-null mice showed a complete loss of pedicle and fusion of the lamina in the vertebrae . In ripply1-deficient zebrafish embryos, somite boundaries do not form, and the characteristic gene expression in the presomitic mesoderm (PSM) is not properly terminated . These results revealed that RIPPLY1 and RIPPLY2 play important roles somite segmentation during development [17,18,19,20].
RIPPLY3, also known as Down syndrome critical region 6 (DSCR6; MIM 609892), plays a critical role in the development of cardiac outflow tract. Previous studies suggested that Ripply3 knockout appeared to lead to the abnormal development of the cardiac outflow tract, including hypotrophy of the aorta and incomplete formation of the ventricular septum in mice . Further evidence showed that RIPPLY3 can repress the intrinsic transcriptional property of TBX1 by recruiting the Groucho/TLE co-repressor . TBX1 loss- and/or gain-of-function in mice led to malformations of the cardiac outflow tract [21,22,23], partly resembling the Ripply3-deficient embryos. However, the relationship between RIPPLY3 and human cardiac malformation has not been investigated.
Since copy number or single-nucleotide variants of TBX1 account for only a small proportion of CTD cases, and its corresponding protein partner RIPPLY3 is also required for normal development of cardiac outflow tract, we hypothesized that deleterious variants in RIPPLY3 gene may be potential molecular mechanisms involved in the pathogenesis of human CTD. Therefore, we screened for RIPPLY3 variants in a CTD cohort and assessed potentially deleterious variants. Here, we report four rare nonsynonymous RIPPLY3 variants in four CTD patients. We also show that the RIPPLY3 variant p.T52S confers consistent functional changes. This is the first report of a human CHD phenotype related to genetic variation in RIPPLY3.
Gene nomenclature information
The RIPPLY3 sequence used in this work is accessible as NCBI RefSeq NM_018962.2. Nucleotide numbering was based on the cDNA sequence (#NM_018962.2) with + 1 corresponding to the A of the ATG initiation codon of translation (exon 1) in the reference sequence.
All assessments were done with the approval of the Medical Ethics Committee of Xinhua hospital and Shanghai Children’s Medical Center. And all experiments were carried out in accordance with the approved guidelines. Fully written informed consent was obtained from all participants (or their parents if children were too young to consent by themselves).
Study subjects and DNA isolation
From November 2011 to January 2014, a cohort of 615 unrelated patients of Han ethnicity diagnosed with CTD were recruited from Shanghai Children’s Medical Center affiliated to Shanghai Jiao Tong University School of Medicine, a cardiologist confirmed the CHD diagnosis for all patients by reviewing and evaluating patient history, physical examinations, and medical records (Table 1) . All subjects were unrelated, and the median age at the time of diagnosis was 10 months with a range of 3 days to 18 years. A total of 391 unrelated healthy individuals (median age was 3.5 years), free of CHD, of Han ethnicity were also enrolled as controls. Peripheral blood samples of all cases were obtained, and genomic DNA was extracted using the QIAmp DNA Blood Mini Kit (Qiagen, Hilden, Germany) with standard protocols.
22q11.2 loci in all CTD patients were explored by CNVplex® (a technique for high-throughput detection of sub-chromosomal copynumber aberrations) as described before . Finally, 577 patients without 22q11.2 deletion/duplication were recruited.
For the patients without 22q11.2 deletion/duplication, fragments covering the promoter region, 5′UTR, coding sequence and splicing site of RIPPLY3 were amplified using the EasyTarget® amplification kit (Genesky Biotechnologies Inc, Shanghai, China) which was developed according to the cycled primer extension and ligation-dependent amplification (CPELA) method. The CPELA was a new method developed by Genesky Biotechnologies for fast and simple enrichment of multiple gene regions for massively parallel sequencing .
All candidate variants of RIPPLY3 were confirmed by Sanger sequencing. The primer pairs used to amplify the coding regions contain candidate variants were designed using Primer Premier 5 in Additional file 1: Table S1. Samples were amplified in a volume of 10 μl containing 1× MyTaq™Mix (Bioline USA Inc.), 200 ng of DNA and 0.4 μM of each primer following the suggested PCR protocol. PCR products were sequencing on an ABI 3730 sequencer (Applied Biosystems, USA). The PCR program was as follows: 98 °C for 1 min; 30 cycles of 98 °C for 20 s, 62 °C for 30 s and 72 °C for 30 s; then 72 °C for 10 min. Amplified products were analyzed on 1% agarose gels stained with GelRed (Biotium, USA). The sequence traces were aligned with the reference sequence for RIPPLY3 (RefSeq NM_018962.2) using the GenBank BLAST program (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Furthermore, because mutations in cardiac kernel members, such as GATA4, NKX2-5, and TBX5, underlie a range of CHDs [25, 26], and TBX1 plays a vital role in the development of the cardiac outflow tract [7,8,9,10,11]. To validate the potential function of identified RIPPLY3 variant, we also screened coding sequences and splicing sites of these known CHD pathogenic genes (TBX1, GATA4, NKX2.5 and TBX5) in CTD patients and controls using the EasyTarget® amplification kit (Genesky Biotechnologies Inc, Shanghai, China). All candidate variants were also validated by Sanger sequencing. The primer pairs used to amplify the coding regions contain candidate variants were designed using Primer Premier 5 in Additional file 1: Table S2. The verified sequence variants were queried in the SNP database at National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov), 1000 Genomes database (http://www.1000genomes.org), Exome Aggregation Consortium database (ExAC, http://exac.broadinstitute.org) and Ensembl database (http://asia.ensembl.org).
In silico analysis
To predict the potential pathogenic impact of missense variants, three different pathogenicity prediction tools (PPTs), namely PolyPhen2, SIFT and Mutation Taster were used.
Alignment of multiple RIPPLY3 protein sequences
Multiple RIPPLY3 protein sequences across various species were aligned using the clustalx 2.1 software.
RIPPLY3 expression vector pCMV6-Entry-RIPPLY3 (Myc-DDK-tagged) containing the cDNA of human RIPPLY3 (RefSeq NM_018962.2) was purchased from Origene (Rockville, MD, USA). The identified variants were introduced into wild type pCMV6-Entry-RIPPLY3 using a QuickChange II Site-Directed Mutagenesis kit (Stratagene, La Jolla, CA, USA). Both wild type and variant inserts were sequenced by Sanger sequencing to confirm the variant sequence and exclude any other sequence variations. TBX1C expression vector pCMV6-XL6-TBX1 containing the cDNA of human TBX1C (RefSeq NM_080647.1) was purchased from Origene (Rockville, MD, USA). TBX1C cDNA was digested and inserting into the plasmid pcDNA3.1(+) at the KpnI and XhoI sites. The luciferase reporter Wnt5a-luc, was constructed by inserting the conserved T-box binding sites (TBEs) in the 3′-UTR of human WNT5A, a target of TBX1 , into the pGL3-promoter plasmid (Promega, USA). An FGF10 luc reporter was also constructed according to the described by Agarwal et al. . FGF10 genomic sequence was obtained from the GenBank (http://www.ncbi.nlm.nih.gov). A fragment comprising 4.5 kb upstream of the coding region (− 4.3 kb/+ 176 bp) of FGF10 was amplified from DNA of a healthy individual and inserted into the MluI/SmaI sites of the pGL3-Basic vector (Promega, Madison, Wisconsin, USA) to generate the FGF10-luc reporter construct. The schemes of these two luciferase reporter plasmids were showed in Fig. 3b, c.
Cell culture and transfection
HEK293T and C2C12 cell lines were purchased from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China) and cultured in the DMEM medium (Invitrogen, California, USA) supplemented with 10% foetal calf serum (Invitrogen, California, USA), penicillin (100 units/ml) and streptomycin (100 μg/ml). Cells were incubated at 37 °C with 5% CO2. All the transient transfections were performed with Fugene HD transfection reagent (Promega, Madison, Wisconsin, USA) according to the manufacturer’s protocol for adherent cells.
Dual luciferase reporter assay
HEK293T or C2C12 cells were seeded into 48-well plates 24 h before transfection. Plasmids including the wild-type or variant RIPPLY3 vector, the TBX1 vector, the luciferase reporters and the corresponding renilla luciferase reporter were co-transfected into HEK293T or C2C12 cellsusing Fugene HD following the manufacturer’s protocol. Transfected cells were incubated for 36 h, then washed using DPBS and lysed using 1× passive lysis buffer for 15 min at room temperature (RT) provided by the Dual Luciferase Reporter Assay Kit (Promega). Firefly luciferase and renilla luciferase activities were measured with the Dual-Glo luciferase assay system (Promega) and the Centro XS3 LB 960 Microplate Luminometer (Berthold, Bad Wildbad, Germany) according to the manufacturer’s recommended protocol. The transfection efficiency was normalized to paired renilla luciferase activity. The results are the means of three independent experiments, each done in duplicate.
HEK293T cells were seeded into 12-well plates 24 h before transfection, and cells were transiently transfected with a total 800 ng pCMV6-Entry-RIPPLY3 or variant RIPPLY3 constructs, and pcDNA3.1-TBX1 was also transfected as positive control. Cells were harvested 24 h after transfection, washed three times with PBS, fixed using 4% paraformaldehyde/PBS for 10 min and washed again in PBS. Cells were permeabilized using 0.3% Triton X-100/PBS for 10 min. After washing with PBS, the cells were blocked with 5% BSA/PBS solution for 1 h at RT. Then, the cells were incubated with a 1:100 solution of the rabbit primary anti-TBX1 antibody (ab18530, abcam) or 1:200 solution of the goat primary anti-Myc tag antibody (ab9132, abcam) in 1% BSA/PBS overnight at 4 °C. Next,the cells were washed with PBS and incubated with a 1:500 solution of the secondary anti-Rabbit IgG H&L (Cy3®) (ab6939, abcam) or 1:200 solution of the secondary Anti-Goat IgG H&L (Alexa Fluor® 488) preadsorbed (ab150129, abcam) in 1% BSA/PBS for 2 h at RT. After washing with PBS, the cells were incubated with DAPI for 6–8 min at RT, then washed with PBS and mounted on a rectangular slide containing an anti-fading agent. The slides were examined using the Olympus BX43 microscope (Olympus, Shinjuku-ku, Tokyo, Japan).
To confirm whether these variants affect the expression of RIPPLY3 protein, western blot analysis was carried out. HEK293T cells were transfected with a total of 500 ng of pCMV6-Entry-RIPPLY3 or variant RIPPLY3 constructs in 24-well plates. Cells were harvested 40–48 h after transfection and lysed in 90 μl cell lysis buffer. After boiling for 10 min, cell lysates were analyzed by sodium dodecyl sulfatepolyacrylamide gel electrophoresis and electrotransferred to Nitrocellulose Blotting Membrane (GE Healthcare). Then, we block the membranes with 5% skim milk for 2 h at RT, and incubated the membranes with anti-Myc tag antibody (ab9132, abcam) and anti-actin antibody (ab3280, abcam) overnight at 4 °C. Next, the membranes were incubated with horseradish peroxidase-conjugated secondary antibodies (Proteintech, Chicago, IL, USA) for 2 h at RT. The membranes were washed, and detected using enhanced chemiluminescence (ECL) reagents (Millipore, Billerica, MA, USA). Bands were visualized with a chemiluminescence detection system (BioRad, Philadelphia, PA, USA) and analyzed using the Imagelab program (BioRad, Philadelphia, PA, USA).
For co-immunoprecipitation experiment, HEK293T cells were cotransfected with a total 7.2 μg pcDNA3.1-TBX1 and pCMV6-Entry-RIPPLY3 or variant RIPPLY3 constructs. The washed Protein A/G Sepharose (GE Healthcare) was mixed with 6 μg of rabbit anti-TBX1 antibody (ab18530, abcam) and placed on a shaker at 4 °C for 3 h. Then, 1 mg protein from HEK 293T cell lysates were added and left at 4 °C for 8 h or overnight. After washing 4 times with PBS, the pellets were resuspended in 30 μl 1× Protein Loading buffer (Sangon Biotech), boiled for 10 min at 100 °C, and the products were analyzed by sodium dodecyl sulfatepolyacrylamide gel electrophoresis, next, a regular Western blot was run as previously described. The grayscale value of protein bands was calculated using Image J software.
Results are expressed as mean ± standard deviation. The independent-samples t-test was used to determine statistical significance of unpaired samples. The Chi square test was used to determine statistical significance of allele frequency between CTD patients and controls. The two-sided statistical tests were considered significant with a level of p < 0.05. All statistical analysis in our study was carried out with SPSS 22.0 (SPSS Inc. Chicago, IL, USA).
Four missense variants of RIPPLY3 were identified in 22q11.2 CNV-negative patients
Four heterozygous missense variants in RIPPLY3 were identified in 4 unrelated CTD patients out of the 577 CTD patients without 22q11.2 deletion/duplication, all the variants were absent from the 391 controls in our cohort: p.P30L in TOF, p.T52S in TOF, p.D113N in TGA/PA/VSD, and p.V179D in PA/VSD and patent ductus arteriosus (PDA) (Table 2).
Within the group of 4 missense RIPPLY3 variants, p.T52S (NM_018962.2:c.155C>G, rs745539198) and p.D113N (NM_018962.2:c.337G>A, rs747419773) were found in ExAC database, the allelic frequencies of these two variants were 3.32e−05 and 1.654e−05 respectively, and the p.T52S variant was found only in Asian population. Notably, variants p.P30L and p.V179D were not identified in the NCBI’s SNP database, 1000 Genomes database, ExAC database or Ensembl database. The chromatograms showing the identified heterozygous RIPPLY3 variants in contrast to corresponding control sequences are shown in Fig. 1a–d, and the positions of these variants in the RIPPLY3 protein are illustrated in Fig. 1e.
As shown in Fig. 1f, a cross-species alignment of multiple RIPPLY3 protein sequences showed that both the p.P30L and p.T52S occur in the highly evolutionarily conserved residues, indicating that the amino acids are functionally important. However, the valine 179 of the protein was just conserved in humans and gorillas, and the aspartate 113 of the protein was not conserved.
Several known CHD pathogenic genes (TBX1, GATA4, NKX2.5 and TBX5) were also screened in these four patients. There were no nonsynonymous variants identified in GATA4, NKX2.5 and TBX5, but two known nonsynonymous single nucleotide polymorphisms (SNPs) were identified in TBX1C (Table 3). The nonsynonymous SNP rs72646967 (TBX1C NM_080647.1:c.1189A>C) appeared commonly in population (MAF > 0.05), and the allele frequency was similar between CHD patients and controls . The another nonsynonymous SNP rs41298838 (TBX1C NM_080647.1:c.928G>A) was previously predicted to be a tolerant polymorphism , and was found in both the non-del22q11 CTD patients and healthy controls at similar frequencies . The allele frequency of rs72646967 and rs41298838 was similar between CTD patients and controls in this study (p > 0.05) (Table 3). Therefore, there were no deleterious variants identified in these four CHD pathogenic genes (TBX1, GATA4, NKX2.5 and TBX5) in patients harboring RIPPLY3 variants.
The variants do not affect the intracellular localization of RIPPLY3
The generated plasmids were transfected into HEK293T cells to investigate the effect of the variants on the trafficking of the RIPPLY3 protein. Immunofluorescence staining showed that both the wild-type RIPPLY3 and variant RIPPLY3 proteins (p.P30L, p.T52S, p.D113N and p.V179D) were located in both the cytoplasm and nucleus (Fig. 2). The TBX1 protein was located only in the nucleus (Additional file 1: Figure S1).
The variants do not affect the expression of RIPPLY3
Wild-type and variant RIPPLY3 expression plasmids were separately transfected into HEK293T cells. Western blot results revealed that there was no significant change in dosage between wild-type and variant RIPPLY3 proteins (Fig. 3a).
The p.P30L, p.T52S and p.V179D variants show impaired inhibition of TBX1 transcriptional activity in vitro
RIPPLY3 is able to repress the transcriptional activity of TBX1 in vitro. To further validate the effects of these variants, we performed dual luciferase reporter assays using the Wnt5a-luc and FGF10-luc in HEK293T cells and C2C12 cells. All of the RIPPLY3 variants are heterozygous in patients, so we also detected whether the variants would affect function of the wild-type RIPPLY3 protein by co-transfected the variant and wild-type RIPPLY3 protein together. Compared to the wild-type RIPPLY3 protein, the p.D113N variant showed no significant change in inhibition of TBX1 transcriptional activity in our assay. The others, namely, p.P30L, p.T52S and p.V179D, showed impaired inhibition of TBX1 transcriptional activity in both 293T cells for Wnt5a-luc (p.P30L: p = 0.0144 for homozygote and p = 0.0258 for heterozygote; p.T52S: p = 0.0031 for homozygote and p = 0.0167 for heterozygote; p.V179D: p = 0.0109 for homozygote and p = 0.0105 for heterozygote, respectively) and C2C12 cells for FGF10-luc. (p.P30L: p = 0.0032 for homozygote and p = 0.0102 for heterozygote; p.T52S: p = 0.0036 for homozygote and p = 0.0098 for heterozygote; p.V179D: p = 0.016 for homozygote and p = 0.0102 for heterozygote, respectively) (Fig. 3b, c). Therefore, the three variants (p.P30L, p.T52S and p.V179D) showed impaired inhibition of TBX1 transcriptional activity in both heterozygote and homozygote, although it seems that RIPPLY3 variants do not have significant dominant negative effect on the wild-type protein.
The p.T52S variant affects the interaction of RIPPLY3 with TBX1
Evidence suggested that RIPPLY3 inhibited TBX1 transcriptional activity by physically interacting with TBX1 in a T-domain-dependent manner . To characterize the effect of the variants on the protein–protein interaction, co-immunoprecipitation assays were performed using equal amounts of protein extracts from HEK293T cells overexpressing the wild-type and variant RIPPLY3 proteins. The results showed that the p.T52S variant significantly disrupted the physical interaction between RIPPLY3 and TBX1 (p = 0.0059). However, the p.P30L, p.D113N and p.V179D variants had no effect on the interaction between RIPPLY3 and TBX1 (Fig. 4).
Dysfunctional RIPPLY3 predisposing to congenital cardiovascular defects has been substantiated in animal models . In mice, Ripply3−/− embryos were born but died with cyanosis within 24 h after birth. In addition, an outstanding characteristic of Ripply3−/− embryos is their almost complete lack of the third and fourth pharyngeal arches . Therefore, Ripply3−/− mouse embryos showed abnormal development of the vascular system, including misshapen major blood vessels and deletion of the aortic arch. In addition, an abnormality in cardiac outflow tract development was also seen in the Ripply3−/− embryos, including hypotrophy of the aorta and incomplete formation of the ventricular septum . The phenotypes partly resembled that of mouse embryos overexpressing Tbx1 [21, 32]. Meanwhile, previous studies have shown that RIPPLY3 interacts with TBX1, and inhibits the transcriptional activity of TBX1 [14, 15]. Taken together, these experimental findings from animals have provided strong evidence that RIPPLY3 plays a pivotal role in cardiovascular morphogenesis, especially in cardiac outflow tract development. Functionally compromised RIPPLY3 may underlie a wide variety of congenital cardiac malformations, including CTD in humans.
In the present study, four rare, heterozygous and nonsynonymous variants of RIPPLY3—p.P30L, p.T52S, p.D113N and p.V179D—were identified in 577 unrelated CTD patients without the 22q11.2 deletion/duplication. The variant NM_018962.2:c.155C>G (p.T52S) is referred as rs745539198 with a 3.32e−05 allelic frequency, and the variant NM_018962.2:c.337G>A (p.D113N) is referred as rs747419773 with a 1.654e−05 allelic frequency. Interestingly, in spite of the very low allelic frequency, the p.T52S variant was only found in the Asian population, suggesting the presence of race specificity. However, variants p.P30L and p.V179D were not found in multiple online human gene variation databases. A cross-species alignment of multiple RIPPLY3 protein sequences showed that both the p.P30L and p.T52S occur in the highly evolutionarily conserved residues, while the p.D113N and p.V179D were not conserved. The p.T52S variant was predicted to be “possibly damaging” by the Polyphen2 software. Further functional analysis revealed that three of these variants, p.P30L, p.T52S, and p.V179D, showed impaired inhibition of TBX1 transcriptional activity in both heterozygote and homozygote, and the p.T52S variant affected the interaction of RIPPLY3 with TBX1 in vitro. Moreover, there were no deleterious variants identified in CHD pathogenic genes (TBX1, GATA4, NKX2.5 and TBX5) in patients harboring RIPPLY3 variants.
The human RIPPLY3 gene is composed of four exons that encode a protein of 190 amino acids and maps to chromosome 21q22.13.RIPPLY3 protein has two highly conserved regions found in all members of the Ripply family: a WRPW motif and a C-terminal Ripply homology domain (also called the bowline-DSCR-LedgerLine conserved ‘BDLC’ region) . The WRPW motif (amino acids 39–42) facilitates interaction with GROUCHO [34, 35]. The ‘BDLC’ domain (amino acids 77–112) is predicted to contact with T-box proteins . RIPPLY3 can recruit the Groucho/TLE co-repressor to TBX1 and control its intrinsic transcriptional activity. It has been established that TBX1 is an upstream regulator of select genes expressed during embryogenesis and cardiac morphogenesis, including the genes that encode Wnt family member 5A  and fibroblast growth factor 10 . Therefore, the functional characteristics of the RIPPLY3 variants may be delineated by analyzing the inhibition of TBX1 transcriptional activity on the luciferase reporter Wnt5a-luc and FGF10-luc.
Dual luciferase reporter assays revealed that variants p.P30L, p.T52S and p.V179D, showed impaired inhibition of TBX1 transcriptional activity in both heterozygote and homozygote in vitro, and indicating that heterozygous of these variants could also be deleterious. It is intriguing that protein–protein interaction was hampered by the p.T52S variant, although this variant is not located in the ‘BDLC’ domain, which is predicted to contact with T-box proteins . This may be ascribed to the distortion of the three-dimensional structure of the variant protein, hindering in such away its interaction with the TBX1 protein. The p.P30L and p.V179D variants also showed impaired inhibition of TBX1 transcriptional activity, however, these two variants seemed to have no effect on the physical interaction of TBX1 and RIPPLY3. This goes in parallel with previously reported variants in GATA4 and MESP1 [37,38,39]. These variants reduced the transcriptional activity of GATA4 or MESP1 without affecting their affinity of binding to DNA or another protein partner, as is the case with the p.P30L and p.V179D variants. It is difficult to say whether the impaired inhibition of TBX1 transcriptional activity with these two variant RIPPLY3 proteins (p.P30L and p.V179D) is due to protein instability or the inability of recruiting other TBX1 partners, such as the Groucho/TLE, to modulate its activity, or a combination of both. Further studies are needed to clarify the specific mechanism.
Additionally, immunofluorescence staining revealed that the subcellular localization of these four variants was not affected. Although much is known about the interaction between RIPPLY3 and TBX1 in the pharyngeal apparatus and heart development , relatively less is known of RIPPLY3 playing a pivotal role in the development of other tissues, such as pre-placodal ectoderm and pancreas. Previously studies showed that RIPPLY3, acting downstream of retinoic acid receptor signaling, is a key player in establishing boundaries in the PPE . RIPPLY3 negatively regulates the proliferation of early endocrine cells, acting as an Insm1-regulated gene enriched in the Pdx1-high cell population . RIPPLY3 is also expressed in pancreatic β-cells, partially in the cytoplasm. Therefore, the subcellular localization of RIPPLY3 may not be limited to the nucleus.
This study has some limitations. Firstly, all functional assays in this study were performed in vitro, transgenic animal models would be helpful to confirm the potential function of these variants. Moreover, the specific mechanism of how the p.P30L and p.V179D variants impaired inhibition of TBX1 transcriptional activity need to be elucidated in the future.
In conclusion, our study suggests that deleterious variants in RIPPLY3 are potential molecular mechanisms involved in the pathogenesis of human CTD. This could have critical implications for clinical practice and genetic counseling. Furthermore, these findings reaffirm the importance of the cardiac transcription factor network and demonstrate that protein binding partners could be related to congenital heart disease. Additional studies that test a larger population and more extensive CHD phenotypes will help further our understanding of the role of RIPPLY3 in CHD molecular pathogenesis.
conotruncal heart defect
congenital heart disease
tetralogy of Fallot
transposition of the great arteries
double outlet of right ventricle
persistent truncus arteriosus
pulmonary atresia with ventricular septal defect
interrupted aortic arch
modified from the multiplex ligation-dependent probe amplification
cycled primer extension and ligation-dependent amplification
T-box binding sites
patent ductus arteriosus
pathogenicity prediction tool
copy number variant
single nucleotide polymorphism
Hoffman JI, Kaplan S, Liberthson RR. Prevalence of congenital heart disease. Am Heart J. 2004;147:425–39.
Johnson TR. Conotruncal cardiac defects: a clinical imaging perspective. Pediatr Cardiol. 2010;31:430–7.
O’Malley CD, Shaw GM, Wasserman CR, Lammer EJ. Epidemiologic characteristics of conotruncal heart defects in California, 1987–1988. Teratology. 1996;53:374–7.
Xu Y, Fang S, Zhang E, Pu T, Cao R, Fu Q, et al. A 3 base pair deletion in TBX1 leads to reduced protein expression and transcriptional activity. Sci Rep. 2017;7:44165.
Lindsay EA. Chromosomal microdeletions: dissecting del22q11 syndrome. Nat Rev Genet. 2001;2(11):858–68.
Robin NH, Shprintzen RJ. Defining the clinical spectrum of deletion 22q11.2. J Pediatr. 2005;147:90–6.
Yagi H, Furutani Y, Hamada H, Sasaki T, Asakawa S, Minoshima S, et al. Role of TBX1 in human del22q11.2 syndrome. Lancet. 2003;362:1366–73.
Buckingham M, Meilhac S, Zaffran S. Building the mammalian heart from two sources of myocardial cells. Nat Rev Genet. 2005;6:826–35.
Liao J, Aggarwal VS, Nowotschin S, Bondarev A, Lipner S, Morrow BE. Identification of downstream genetic pathways of Tbx1 in the secondary heart field. Dev Biol. 2008;316:524–37.
Chen L, Fulcoli FG, Tang S, Baldini A. Tbx1 regulates proliferation and differentiation of multipotent heart progenitors. Circ Res. 2009;105:842–51.
Xu YJ, Chen S, Zhang J, Fang SH, Guo QQ, Wang J, et al. Novel TBX1 loss-of-function mutation causes isolated conotruncal heart defects in Chinese patients without 22q11.2 deletion. BMC Med Genet. 2014;15:78.
Zhang X, Xu Y, Liu D, Geng J, Chen S, Jiang Z, et al. A modified multiplex ligation-dependent probe amplification method for the detection of 22q11.2 copy number variations in patients with congenital heart disease. BMC Genomics. 2015;16:364.
Kondow A, Hitachi K, Okabayashi K, Hayashi N, Asashima M. Bowline mediates association of the transcriptional corepressor XGrg-4 with Tbx6 during somitogenesis in Xenopus. Biochem Biophys Res Commun. 2007;359:959–64.
Kawamura A, Koshida S, Takada S. Activator-to-repressor conversion of T-box transcription factors by the Ripply family of Groucho/TLE associated mediators. Mol Cell Biol. 2008;28:3236–44.
Okubo T, Kawamura A, Takahashi J, Yagi H, Morishima M, Matsuoka R, et al. RIPPLY3, a Tbx1 repressor, is required for development of the pharyngeal apparatus and its derivatives in mice. Development. 2011;138:339–48.
Janesick A, Shiotsugu J, Taketani M, Blumberg B. RIPPLY3 is a retinoic acid-inducible repressor required for setting the borders of the pre-placodal ectoderm. Development. 2012;139:1213–24.
Takahashi J, Ohbayashi A, Oginuma M, Saito D, Mochizuki A, Saga Y, et al. Analysis of Ripply1/2-deficient mouse embryos reveals a mechanism underlying the rostro-caudal patterning within a somite. Dev Biol. 2010;342:134–45.
Chan T, Kondow A, Hosoya A, Hitachi K, Yukita A, Okabayashi K, et al. Ripply2 is essential for precise somite formation during mouse early development. FEBS Lett. 2007;581:2691–6.
Kawamura A, Koshida S, Hijikata H, Ohbayashi A, Kondoh H, Takada S. Groucho-associated transcriptional repressor ripply1 is required for proper transition from the presomitic mesoderm to somites. Dev Cell. 2005;9(6):735–44.
Morimoto M, Sasaki N, Oginuma M, Kiso M, Igarashi K, Aizaki K, et al. The negative regulation of Mesp2 by mouse Ripply2 is required to establish the rostro-caudal patterning within a somite. Development. 2007;134:1561–9.
Liao J, Kochilas L, Nowotschin S, Arnold JS, Aggarwal VS, Epstein JA, et al. Full spectrum of malformations in velo-cardio-facial syndrome/DiGeorge syndrome mouse models by altering Tbx1 dosage. Hum Mol Genet. 2004;13:1577–85.
Zhang Z, Baldini A. In vivo response to high-resolution variation of Tbx1 mRNA dosage. Hum Mol Genet. 2008;17:150–7.
Vitelli F, Huynh T, Baldini A. Gain of function of Tbx1 affects pharyngeal and development in the mouse. Genesis. 2009;47:188–95.
Zhang Z, Li C, Wu F, Ma R, Luan J, Yang F, et al. Genomic variations of the mevalonate pathway in porokeratosis. eLife. 2015;4:e06322.
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.
Schott JJ, Benson DW, Basson CT, Pease W, Silberbach GM, Moak JP, et al. Congenital heart disease caused by mutations in the transcription factor NKX2-5. Science. 1998;281:108–11.
Chen L, Fulcoli FG, Ferrentino R, Martucciello S, Illingworth EA, Baldini A. Transcriptional control in cardiac progenitors. Tbx1 interacts with the BAF chromatin remodeling complex and regulates Wnt5a. PLoS Genet. 2012;8:e1002571.
Agarwal P, Wylie JN, Galceran J, Arkhitko O, Li C, Deng C, et al. Tbx5 is essential for forelimb bud initiation following patterning of the limb field in the mouse embryo. Development. 2003;130:623–33.
Guo T, McDonald-McGinn D, Blonska A, Shanske A, Chow E, Bassett AS, et al. Genotype and cardiovascular phenotype correlations with TBX1in 1,022 velo-cardio-facial/DiGeorge/22q11.2 deletion syndromepatients. Hum Mutat. 2011;32(11):1278–89.
Heike CL, Starr JR, Rieder MJ, Cunningham ML, Edwards KL, Stanaway IB, et al. Single nucleotide polymorphism discovery in TBX1 in individuals with and without 22q11.2 deletion syndrome. Birth Defects Res A. 2010;88(1):54–63.
Yue-Juan Xu, Wang Jian, Rang Xu, Zhao Peng-Jun, Wang Xi-Ke, Sun Heng-Juan, et al. Detecting 22q11.2 deletion in Chinese children with conotruncal heart defects and single nucleotide polymorphisms in the haploid TBX1 locus. BMC Med Genet. 2011;12:169.
Merscher S, Funke B, Epstein JA, Heyer J, Puech A, Lu MM, et al. TBX1 is responsible for cardiovascular defects in velo-cardio-facial/DiGeorge syndrome. Cell. 2001;104:619–29.
Kondow A, Hitachi K, Ikegame T, Asashima M. Bowline, a novel protein localized to the presomitic mesoderm, interacts with Groucho/TLE in Xenopus. Int J Dev Biol. 2006;50:473–9.
Paroush Z, Finley RL Jr, Kidd T, Wainwright SM, Ingham PW, Brent R, et al. Groucho is required for Drosophila neurogenesis, segmentation, and sex determination and interacts directly with hairy-related bHLH proteins. Cell. 1994;79:805–15.
Fisher AL, Ohsako S, Caudy M. The WRPW motif of the hairy-related basic helix-loop-helix repressor proteins acts as a 4-amino-acid transcription repression and protein-protein interaction domain. Mol Cell Biol. 1996;16:2670–7.
Kawamura A, Koshida S, Hijikata H, Ohbayashi A, Kondoh H, Takada S. Groucho-associated transcriptional repressor ripply1 is required for proper transition from the presomitic mesoderm to somites. Dev Cell. 2005;9:735–44.
Nemer G, Fadlalah F, Usta J, Nemer M, Dbaibo G, Obeid M, et al. A novel mutation in the GATA4 gene in patients with Tetralogy of Fallot. Hum Mutat. 2006;27:293–4.
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:1662–71.
Werner P, Latney B, Deardorff MA, Goldmuntz E. MESP1 mutations in patients with congenital heart defects. Hum Mutat. 2016;37:308–14.
Osipovich AB, Long Q, Manduchi E, Gangula R, Hipkens SB, Schneider J, et al. Insm1 promotes endocrine cell differentiation by modulating the expression of a network of genes that includes Neurog3 and RIPPLY3. Development. 2014;141:2939–49.
KS and YX designed the study. NH drafted the manuscript and all authors contributed to the revision of the manuscript. All authors read and approved the final manuscript.
We thank all the participants for their contribution to this work. We would like to thank our study team in the Department of Pediatric Cardiology, Xinhua hospital, affiliated to Shanghai Jiao Tong University School of Medicine, for their efforts and excellent work. We are also indebted to the senior management and Medical Ethics Committee of Xinhua hospital and Shanghai Children’s Medical Center for their unwavering support and valuable recommendations. Moreover, we thank postdoctor Alice Lau from university of Ottawa (Canada) for helping us to re-edit our manuscript. Finally, we are grateful to grants from the National Natural Science Foundation of China, Shanghai Municipal Commission of Health and Family Planning and Shanghai Hospital Development Center for financial support of this Research Project.
The authors declare that they have no competing interests.
Availability of data and materials
The data sets analyzed within the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
All assessments were done with the approval of the Medical Ethics Committee of Xinhua hospital and Shanghai Children’s Medical Center. And all experiments were carried out in accordance with the approved guidelines. Fully written informed consent was obtained from all participants (or their parents if children were too young to consent by themselves).
This work was supported by grants from the National Natural Science Foundation of China (81670285, 81300068, 81771623, 81741066), a Grant (GWIV-23) from Shanghai Municipal Commission of Health and Family Planning, and a Grant (SHDC12015102) from Shanghai Hospital Development Center.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.