Next Article in Journal
The Gut Microbiota Profile in Children with Prader–Willi Syndrome
Previous Article in Journal
Anaphase Bridges: Not All Natural Fibers Are Healthy
Previous Article in Special Issue
Pitt-Hopkins Syndrome: Clinical and Molecular Findings of a 5-Year-Old Patient
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 8
10.3390/genes11080903
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate

by
Beau Sylvester
1,
Frederick Brindopke
2,
Akiko Suzuki
3,
Melissa Giron
4,
Allyn Auslander
1,5,
Richard L. Maas
6,
Becky Tsai
7,
Hanlin Gao
7,
William Magee III
1,
Timothy C. Cox
3,8,† and
Pedro A. Sanchez-Lara
9,*,†
1
Division of Plastic and Maxillofacial Surgery, Children’s Hospital Los Angeles, Los Angeles, CA 90027, USA
2
Operation Smile Global Headquarters, Virginia Beach, WV 23453, USA
3
Department of Oral & Craniofacial Sciences, University of Missouri-Kansas City School of Dentistry, Kansas City, MO 64108, USA
4
Operación Sonrisa Honduras, Tegucigalpa 11101, Honduras
5
Department of Preventive Medicine, Keck School of Medicine of the University of Southern California, Los Angeles, CA 90033, USA
6
Division of Genetics, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA
7
Fulgent Genetics, Temple City, CA 91780, USA
8
Department of Pediatrics, University of Missouri-Kansas City School of Medicine, Kansas City, MO 64108, USA
9
Department of Pediatrics, Cedars-Sinai Medical Center, David Geffen School of Medicine at UCLA, Los Angeles, CA 90048, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally.
Genes 2020, 11(8), 903; https://doi.org/10.3390/genes11080903
Submission received: 19 June 2020 / Revised: 28 July 2020 / Accepted: 4 August 2020 / Published: 7 August 2020
(This article belongs to the Special Issue Molecular Genetics of Facial Traits and Malformations)
Browse Figures
Review Reports Versions Notes

Abstract

:
Missense, nonsense, splice site and regulatory region variants in interferon regulatory factor 6 (IRF6) have been shown to contribute to both syndromic and non-syndromic forms of cleft lip and/or palate (CL/P). We report the diagnostic evaluation of a complex multigeneration family of Honduran ancestry with a pedigree structure consistent with autosomal-dominant inheritance with both incomplete penetrance and variable expressivity. The proband’s grandmother bore children with two partners and CL/P segregates on both sides of each lineage. Through whole-exome sequencing of five members of the family, we identified a single shared synonymous variant, located in the middle of exon 7 of IRF6 (p.Ser307Ser; g.209963979 G>A; c.921C>T). The variant was shown to segregate in the seven affected individuals and through three unaffected obligate carriers, spanning both sides of this pedigree. This variant is very rare, only being found in three (all of Latino ancestry) of 251,352 alleles in the gnomAD database. While the variant did not create a splice acceptor/donor site, in silico analysis predicted it to impact an exonic splice silencer element and the binding of major splice regulatory factors. In vitro splice assays supported this by revealing multiple abnormal splicing events, estimated to impact >60% of allelic transcripts. Sequencing of the alternate splice products demonstrated the unmasking of a cryptic splice site six nucleotides 5′ of the variant, as well as variable utilization of cryptic splice sites in intron 6. The ectopic expression of different splice regulatory proteins altered the proportion of abnormal splicing events seen in the splice assay, although the alteration was dependent on the splice factor. Importantly, each alternatively spliced mRNA is predicted to result in a frame shift and prematurely truncated IRF6 protein. This is the first study to identify a synonymous variant as a likely cause of NS-CL/P and highlights the care that should be taken by laboratories when considering and interpreting variants.

1. Introduction

Clefts of the lip and/or palate (CL/P) are among the most common and recognizable inborn errors of development, presenting in one in 700 live births [1]. Orofacial clefting often occurs as a component of various genetic syndromes and sequences, but approximately 70% of CL/P cases arise in apparent isolation [1]. These individuals with non-syndromic CL/P (NS-CL/P) display complex epidemiology, varying in prevalence across geographic regions and ethnic groups, and complex etiology, with numerous genetic and environmental risk factors contributing to disease [1].
CL/P and isolated cleft palate have traditionally been treated as etiologically independent [2]. However, families in which some individuals present with cleft lip and others with cleft palate only have been noted with shared mutations in a handful of genes, such as interferon regulatory factor 6 (IRF6) [3], tumor protein p63 (TP63) [4,5], fibroblast growth factor receptor 1 (FGFR1) [6], and msh homeobox 1 (MSX1) [7]. Furthermore, studies of large, multiplex families with cases of NS-CL/P have demonstrated that there are likely sub-clinical phenotypes in otherwise “unaffected” relatives, suggesting that the spectrum of disease may be broader than previously recognized [8,9,10,11].
The identification of gene variants underlying NS-CL/P has mainly been accomplished through studies of individuals, families, and larger population samples with NS-CL/P using genome-wide association studies (GWASs), linkage analysis, and whole genome and exome sequencing. Additional gene candidates have been identified through developmental and gene expression analyses in animal models [1] and inference from individual syndromes and genes disrupted in structural chromosomal anomalies [12]. The best example of this is the IRF6 gene, which was first identified as the cause of Van der Woude syndrome [3], a disorder defined by the co-presentation of CL/P and lower lip pits, and later shown to contribute to NS-CL/P in multiple family-based studies and GWASs [13,14,15,16,17]. Animal model studies have demonstrated that a deficiency of IRF6 is associated with defects in epithelial adhesion and differentiation [18,19,20,21,22], providing potential clues for disease pathophysiology.
In this study, we performed whole exome sequencing (WES) on several family members from a large, multiplex Honduran family with NS-CL/P, identifying a novel segregating synonymous variant in IRF6. We demonstrate that this synonymous variant impacts an exonic splice silencer element that results in the utilization of multiple cryptic exonic and intronic splice sites, each of which is predicted to result in a frame shift and loss of function of IRF6. Although similar types of mutations have been reported in other diseases, this represents the first example underlying the presentation of CL/P.

2. Materials and Methods

2.1. Editorial Policies and Ethical Considerations

Ethical approval for this study was obtained from the Institutional Review Board of Children’s Hospital Los Angeles under CCI-08-00150. Informed consent was obtained from all participating human subjects prior to their inclusion in the study. This study was performed in compliance with the Declaration of Helsinki and the U.S. Federal Policy for the Protection of Human Subjects.

2.2. Patient Recruitment and Clinical Evaluation

The proband (individual III-6; Figure 1) was first identified through Operación Sonrisa (Operation Smile) Honduras during regular pre-surgical programmatic activities for screening children with cleft lip and palate. Operación Sonrisa was established in 1997 and is a country-based foundation of the international NGO Operation Smile, which has been providing free cleft care globally for almost four decades. Upon recognition that this family included multiple affected individuals across several generations, the organization flagged the family as a strong candidate for research. A three-generation pedigree was obtained, individuals were consented for research and phenotypically screened for other comorbidities or subtle signs of an underlying syndrome.

2.3. Samples and DNA Extraction

Saliva samples were collected from available family members using Oragene DNA Genotek 500 adult and 575 saliva collection kits. Genomic DNA was extracted from 400 μL peripheral blood using a QIAamp DNA Mini Kit (QIAGEN, Düsseldorf, Germany) following the manufacturer’s instructions. The quality and quantity of the DNA samples were assessed using a NanoDrop 8000 (Thermo Scientific™, Waltham, MA, USA) and agarose gel electrophoresis.

2.4. Whole Exome Sequencing

Genomic DNA libraries were prepared using TruSight One (Illumina Inc., San Diego, CA, USA) and targets were captured using a custom probe set mixture based on TruSight One. Massively parallel DNA sequencing was initially performed on five individuals (I-2, II-1, II-2, III-3, III-6) using HiSeq 4000 instruments and chemistry (Illumina Inc.). Sequence reads were aligned to the human genome using Novoalign (Novocraft, Selangor, Malaysia). Single nucleotide variants and small indel variants were called using VarScan2 and annotated using a combination of in-house tools and Alamut Batch (Interactive Biosciences division of Sophia Genetics).

2.5. Mutation Screening and Sanger Confirmation

Based on the family history and pedigree structure, autosomal dominant Mendelian inheritance with incomplete penetrance was considered the best model and variants shared among affected individuals were assessed with greater scrutiny. Variants with a sequencing depth of coverage <10x or genotype quality (Q < 20) were excluded from analysis. Only rare variants with a minor allele frequency of <1% in the 1000 Genomes Project (www.1000genomes.org) or Exome Sequencing Project (ESP) (esp.gs.washington.edu/drupal/) reference populations were included for analysis. Each variant was evaluated for the potential to contribute to the orofacial cleft presentation based on a single-gene Mendelian model. The literature was reviewed for individual variants and candidate genes with potential significant relationships to the phenotype and possible deleterious effects on craniofacial structures.
Sanger sequencing with forward (5′-GCCTCCTGCCACCACT-3′) and reverse (5′-GGCTAGATTCCCTAGACCA-3′) primers designed using Primer5 ver. 0.4.0 was applied to the candidate variant. The Sanger sequencing product was separated using an ABI 3730 genetic analyzer, and the data were analyzed using Chromas V1.0.0.1.

2.6. In Silico Analysis

The online software Human Splicing Finder (version 3.1; https://www.genomnis.com/the-system-1) [23] was used to analyze the genomic sequence in and around exon 7 of human IRF6. Specifically, this tool serves two functions: (1) to identify potential exonic and intronic splicing enhancer sequences, silencer sequences and cryptic splice donor and acceptor sites in a given sequence in part based on the similarity to the respective consensus binding sites of major splice regulatory factors; and (2) determine the possible impact on any such sequences by a given variant. Only those sites above the default thresholds of the software were considered as potentially relevant.

2.7. Minigene Splice Assay

A fragment of the human IRF6 gene, comprising exon 7 (of nine) and 309 bp of intron 6 and 428 bp of intron 7, was amplified by polymerase chain reaction (PCR) using HiFi HotStart DNA Polymerase (Kapa Biosystems) and human genomic DNA isolated from an unaffected control individual with the following primer pair containing restriction site linkers: hIRF6-I6 forward (5′-CAATCTCGAGCCGACTCAGTCAG TATAGCGTGGG-3′; XhoI site underlined) and hIRF6-I7 reverse (5′-AGCTCTAGATGTAACTCCCTTC TTTGTTGCC-3′; XbaI site underlined). The purified product was double digested with XhoI and XbaI and cloned into a similarly digested pET01 exon trap vector (MoBiTec, Goettingen, Germany) to generate the wildtype (WT; c.921C) construct.
  • The c.921C>T variant was introduced by overlap PCR with high fidelity polymerase and the following primer sets: hIRF6-I6 forward and hIRF6-synon reverse (5′-CATGACCACTGACCTCCAGGATCAG-3′), and hIRF6-synon forward (5′-GGAGGTCAGTGGTCATGCCATTTATG-3′) and hIRF-I7 reverse. All constructs were verified by sequencing.
  • Wildtype and mutant vectors were transfected into COS7 cells (CRL-1651, ATCC) by Lipofectamine 3000 (Thermo Fisher Scientific) according to the manufacturer’s protocol (n = 3 per group). Total RNA was extracted two days post transfection using QIAshredder and an RNeasy Micro Kit (Qiagen) and then cDNA was synthesized using Superscript IV (Thermo Fisher Scientific). To detect splicing patterns, PCR was carried out using pET01-ETPR04 forward (5′-GGATTCTTCTACACACCC-3′) and pET01-ETPR05 reverse (5′-GTTGACCTCGACCCACCT-3′). As a control for DNA contamination, PCR was also performed using the RNA as a template. An aliquot of each PCR was separated on a 1.2% agarose gel, stained with ethidium bromide, and the gel image capture on an C400 imager (Azure Biosystems, Dublin, CA, USA) for the quantification of individual bands. Quantification was performed in the linear range using the Azure software and manually adjusted for band size. The remainder of the amplification reaction was run on a preparative 1.2% agarose gel and the individual products were extracted using a Qiagen gel extraction kit (Qiagen) and then Sanger sequenced (Genewiz, South Plainfield, NJ, USA) using the amplifying primers.
Plasmids for the expression of Splice Regulatory Factor 1 (pcDNA-FLAG-SF2; Addgene 99021) and Splice Regulatory Factor 2 (pcDNA3.1-SC35-cMyc; Addgene 44721) were separately co-transfected with the empty pET01 vector, or either the WT (c.921C) or mutant (c.921T) vectors into COS7 cells and total RNA was extracted two days post transfection (n = 6 per group), as described above. Following cDNA synthesis, PCR products were amplified as previously described, separated on a 1.2% agarose gel and quantified as above.

3. Results

3.1. Clinical Findings

The pedigree of the family identified by Operación Sonrisa Honduras contained seven affected individuals and two unaffected obligate carriers, as well as the grandmother (individual I-2) who bore offspring with two separate partners, with clefting segregating on both sides of the pedigree (Figure 1). The proband (individual III-6) presented with unilateral cleft lip and cleft palate as did one of his paternal uncles (individual II-3). An additional paternal uncle (individual II-10) and a cousin (individual III-11) presented with bilateral cleft lip and cleft palate, while two additional cousins (individuals III-7 and III-12) and a sibling (III-3) presented with cleft palate only. No family member exhibited lip pits or other clinical signs suggestive of a syndrome. Strikingly, all overtly affected members of the family are males.

3.2. Exome Sequencing and Variant Analysis

Exome sequencing was performed on DNA samples from individuals I-2, II-1, II-2, III-3, and III-6, and the data from each individual were filtered for quality and coverage and then selected for variant allele frequencies of <1%. Under a dominant inheritance model, a total of 411 variants were found to be shared between affected individuals III-6, III-3, and the presumed obligate carriers II-1 and I-2. To further narrow down candidates, we applied the machine learning algorithm, DOMINO, to select only those variants in genes determined to be most likely to result in a dominant condition [24] and then focused on only those variants deemed to be very rare (i.e., minor allele frequencies of <1 × 10−4). This analysis narrowed the number of shared gene variants to 27 (Table 1), which included 13 intronic variants (between 3 and 20 bp from a splice site), three variants located in 3′UTRs, one six-base pair in-frame duplication, four synonymous changes, and six missense variants. The six-base pair duplication (in SPEN) was subsequently determined to be a much more common variant but had been bypassed in the commercial bioinformatic pipeline because it is located in a low complexity region and thus was excluded from further consideration. Of the 26 remaining variants, fourteen of the genes have previously been implicated in diseases. The genes associated with ten of these diseases were excluded as the diseases involved principally non-craniofacial tissues. Three of the disease genes are associated with variable facial dysmorphism as part of a broader spectrum of defining features and thus were deemed less likely to play a role in our family. The remaining disease gene was IRF6, the most prominent CL/P gene and thus an excellent candidate for the family’s presentation. After exclusion of the 13 genes implicated in other diseases, as well as the remaining variants located in introns and untranslated regions, we were left with seven initial candidate genes: SYT11, KIAA1468, ILF3, IRF6, PDE7A, MCL1, and ITGB1. The variant in KIAA1468 was a single nucleotide duplication in a splice site and was deemed unlikely to be pathogenic. The missense variant in ILF3 (encoding interleukin enhancer binding factor 3) alters one of a stretch of five glycine residues that is positioned within the undefined C-terminus of the protein. As the affected glycine was one of 82 glycines in the final 285 amino acids of the protein (i.e., a very glycine-rich region), we deemed this variant unlikely to be causal. The missense variant in ITGB1 is located in a rarely utilized exon which harbors a relatively high proportion of variants compared to other exons and so was excluded. The variants in the remaining four genes, SYT11, PDE7A, MCL1, and IRF6, were synonymous changes. SYT11 encodes synaptotagmin X1, a protein involved in synaptic vesicle regulation. This gene is considered a candidate for schizophrenia. PDE7A encodes a high-affinity cAMP-specific phosphodiesterase that has been shown to be important for heart, muscle, and T-cell function. MCL1 encodes a regulator of apoptosis that is a target of an approved drug treatment for myeloid cell leukemia. The remaining synonymous variant in IRF6, the most prominent CL/P gene, was therefore deemed the top candidate. The IRF6 variant also had the highest scaled combined annotation dependent depletion (CADD) score of the four synonymous changes: 14.48 (https://cadd.gs.washington.edu/) [25].
The synonymous IRF6 variant is located in the middle of exon 7 (p.Ser307Ser; g.209963979G>A; c.921C>T) and is only found in three of 251,352 alleles in gnomAD (minor allele frequency: 1.19 × 10−5), with all three of Latino ancestry. Importantly, Sanger sequencing validated the presence of the IRF6 variant and confirmed segregation in an additional affected individual (III-7), thus confirming the segregation of the variant with the phenotype across both sides of the pedigree. Unfortunately, DNA samples were not available from any of the other affected members of the family.

3.3. Functional Testing of the IRF6 Variant

Initial visual inspection of the synonymous IRF6 variant suggested it did not create a splice acceptor/donor site. However, in silico analyses using the Human Splicing Finder online tool, which uses multiple distinct algorithms, revealed that three separate algorithms predicted the variant would create a new exonic splice silencer (ESS) element (Figure 2A), thus potentially disrupting normal mRNA splicing.
Using an established in vitro minigene splice assay supported the prediction model of disrupted mRNA splicing. Specifically, four distinct splicing events were detected when the variant minigene construct was transfected into COS7 cells, one of which corresponded to the product size seen when the wildtype minigene construct was used (Figure 2(Bi)). The fluorescence intensity of each PCR product was determined using an Azure Biosystems C400 gel imaging system and was then adjusted for band size, which estimated the alternate products to constitute approximately 60% of the allelic transcripts (Figure 2(Bii)).
The direct sequencing of each band following gel purification revealed the nature of each abnormal splicing event and confirmed that the upper band represented a correctly spliced product. Each of the three lower bands shared the use of a cryptic splice donor located six nucleotides 5′ of the variant (AG915/GT), i.e., the c.921T variant located at the +5 position of a perfect consensus splice donor that is normally masked. In the lower band (band 4), this was the only difference in splicing from the expected wildtype band (band 1). For bands 2 and 3, they each shared a new cryptic splice acceptor site 135 nucleotides 5′ of the canonical exon 7 acceptor (i.e., within intron 6). This cryptic site conforms well to the splice acceptor consensus: 5′-(Y)11-13-N-(Y)-AG-3′ (where Y = a pyrimidine and N = any base): 5′-TCTCCCTTTATTGCAG|GA-3′; the vertical line marking the splice site [26]. The distinction between bands 2 and 3 is that band 3 contained only 44 nucleotides of the intronic sequence, whereas band 2 contained all 135 nucleotides up to the original start of exon 7. The inclusion of just part of the intron 6 sequence in band 3 is because of the utilization of a cryptic donor site (5′-AG|GTGGC-3′; vs. consensus: 5′-AG|GTRAG-3′) [26] that is then spliced to the canonical exon 7 acceptor site. The splicing pattern revealed by the sequencing of each product is schematically represented in Figure 2C. Importantly, each alternatively spliced mRNA is predicted to result in a frame shift and prematurely truncated IRF6 protein.
Somewhat surprisingly, alignment of the orthologous IRF6 sequence from 100 vertebrate species found that all 59 mammalian species examined had the C allele, whereas 22 of 40 non-mammalian vertebrates (including marsupials, birds, reptiles and fishes) did not have a C at position c.921 or its equivalent position (see Supplemental Figure S1). Note: the sequence over this region was not available for one vertebrate species. In all but two of these species, the C is substituted for a T (as seen in our family). Yet, notably, in seven of these species, the T substitution is associated with a corresponding substitution of the either G915 or G916 in the cryptic splice donor site (see Supplemental Figure S1). Under the splicing models, this corresponding substitution is predicted to destroy the cryptic exonic splice site, thus countering the potential impact of the C>T change at position c.921.

3.4. The Splicing Pattern Is Affected by Ectopic Expression of Major Splice Regulatory Factors

In addition to the prediction that the synonymous variant created an ESS element, the Human Splice Finder tool also identified consensus binding sites for two major splice regulatory proteins, SRSF1 and SRSF2. Each of these factors can have both positive and negative effects on splicing depending on the context of their binding. We therefore decided to test whether the ectopic expression of either of these factors would impact the pattern or degree of alternatively spliced products seen with the variant minigene. As anticipated, the expression of either factor had no noticeable impact on the splicing of the wildtype minigene. With the variant minigene, however, the co-expression of SRSF1 improved the differential percentage between normal and abnormal splice products (mean difference = 32.93% (SD = 2.41) without SRSF1 vs. 25.53% (SD = 1.47) with SRSF1; two-tailed t-test: p = 0.000159). In contrast, the co-expression of SRSF2 exacerbated the difference between normal and abnormal splice products (mean difference = 43.57% (SD = 2.28) with SRSF2 vs. 32.93% without; p = 0.000030) (see Figure 3).

4. Discussion

Synonymous coding region changes, especially those away from utilized canonical splice sites, are often overlooked as primary causes of disease, as their impact has historically been difficult to predict and/or prove. However, in the recent medical literature, in large part due to improved bioinformatic algorithms and more sensitive analytical techniques, the number of examples where identified synonymous changes have been implicated in various human diseases has been growing. Diseases in which such changes have been shown to contribute include cancer [27,28,29,30,31,32,33,34,35,36], cystic fibrosis [37,38], Crohn’s [39,40], temporomandibular joint disorder [41], spinal muscular atrophy [42], congenital myasthenic syndrome [43], and osteoporosis [44]. Researchers have adopted a range of approaches to investigate the consequences of synonymous mutations, including population genetics models that demonstrate the negative selection of synonymous mutations at exonic splice enhancers (ESEs) [45,46,47,48], RNA sequencing and transcript expression analyses [27,29,30,31,32,34,35,36,37,39,43,44], and functional models that examine resulting protein and mRNA secondary structure and molecular interactions [29,31,35,37,38,39,41,44,49,50,51,52]. As a result, coding-independent mutations have been shown to disrupt gene regulation through a variety of mechanisms [28], including the disruption of splice sites and subsequent exon skipping [27,30,32,34,35,42,50,51,52], the creation of new splice sites [33], the loss of microRNA (miRNA) binding sites [29,39], protein misfolding [31,37,38], and mRNA secondary structure defects and instability [36,41].
Here, we describe the results of whole exome sequencing conducted on a multi-generational family of Honduran ancestry, in which multiple family members present with either isolated cleft palate or cleft lip with cleft palate. In what we believe is the first example for NS-CL/P, we identify a synonymous variant (c.921C>T; S307S) in a well-established cleft-causing gene, IRF6, that segregates with the orofacial clefting phenotype in all members that were available for testing. While the synonymous variant does not destroy a canonical splice site or create a new splice site, bioinformatic analysis suggested it created an exonic splice silencer (ESS) element. ESS elements are thought to bind heterogeneous nuclear ribonucleoproteins (hnRNPs) and generally function to inhibit neighboring splice site usage by interfering with components of the splicing complex. In doing so, they can promote the use of alternate splice sites. However, predicting the actual impact on splicing using bioinformatic algorithms is far from reliable [53]. Consequently, we chose to directly assess the impact of the variant using an established minigene splice assay [53]. This assay demonstrated a clear and significant disruption to the normal splicing of exon 7, detecting three principle aberrant splicing products. Sequencing of the individual products showed that each shared the unmasking of a cryptic splice site six nucleotides 5′ of the variant and suppression of the normal 3′ exon 7 splice site. Two of the products were also involved the concomitant utilization of a cryptic splice acceptor sequence located in intron 6, 135 nucleotides 5′ of the native exon 7 splice acceptor. One of these two products also utilized a cryptic donor site 44 nucleotides 3′ of the cryptic acceptor to join with the established exon 7 acceptor site. Each of the three abnormal splice products is predicted to result in a frame shift and premature truncation, and thus likely to result in nonsense-mediated decay (NMD) [54]. Together, the anomalous splicing impacted over 60% of the allelic splice products in the splicing assay. An important caveat to this result is that this assay is somewhat artificial and may not reflect the actual percentage of abnormal splicing of the endogenous IRF6 locus. This would be most accurately determined using quantitative or real-time PCR using RNA from suitable tissue samples of the affected individuals [53], which unfortunately were not available. Nevertheless, these data do emphasize that the abnormal products represent a significant fraction of the total amount of IRF6 transcript and that the loss of functional IRF6 protein in heterozygous individuals would also likely be significant.
ESS elements can also work in concert with exon splice enhancers (ESEs), as well as intronic silencers and enhancers (ISSs and ISEs) [53,55]. Indeed, some ESS elements can be juxtaposed with ESE elements to finely control splicing events, as in the case of exon 2 splicing of the HIV-1 tat gene [56]. Notably, binding sites for a number of serine/arginine splice factors (SRSF), which are splice regulatory proteins that bind ESEs, were also predicted to be centered over c.921 of IRF6. We therefore tested whether the ectopic expression of these SRSFs would mitigate the impact of the c.921T variant. Consistent with its known role in influencing splice site selection in a concentration-dependent manner [53], the ectopic expression of SRSF1 (ASF/SF2) led to a slightly improved ratio of normal to abnormal splice products. In contrast, SRSF2 (Sc35), which has a role in forming a bridge between the splicing components bound to the donor and acceptor sites, slightly exacerbated the abnormal splicing associated with the c.921T variant. These studies reinforce the complex nature of pre-mRNA splicing and highlight that much more work is needed to dissect out the impact of not only the splice regulatory factors but also the cis elements that govern their binding.
A number of synonymous variants have previously been linked to orofacial clefting. Letra et al. identified a c.471C>T (rs8061351; p.Pro157Pro) variant in the CRISPLD2 gene through family-based association and found that the C allele was significantly overrepresented in CL/P cases [57]. In silico analysis predicted the alteration of an ESE motif, and also created a putative binding site for transcription factor AP-2α, an important regulator of craniofacial development. However, no functional studies were performed to verify either possibility. In a case–control association study of NS-CL/P in an Indian population, Kumari et al. found an association with a synonymous c.348C>T variant (rs34165410; p.Gly116Gly) in the MSX1 gene [58]. In silico modeling of the NS-CL/P-associated variant suggested the variant may create a new splice site and/or affect expression through decreased codon usage. Functional testing of the variant in a minigene splicing assay failed to demonstrate an impact on splicing. More recently, Gaczkowska et al. identified nucleotide variants in the PAX7 gene that were overrepresented in NS-CL/P through a GWAS of a Polish population [59]. In a subsequent evaluation of the PAX7 gene in isolated NS-CL/P cases, they identified a synonymous substitution (c.87G>A, p.Val29Val) in a single patient. In silico modeling predicted the disruption of a potential ESE motif in the NS-CL/P-associated variant but this was not functionally tested. Therefore, we believe our case represents the first such instance of a demonstrable functional impact of a synonymous variant in NS-CL/P. Functionally impactful synonymous variants, including those affecting exonic splice regulatory elements, have been previously described in a number of disorders, including cancer [30,35,60,61], neurofibromatosis type 1 [62], severe combined immunodeficiency [63,64], familial dilated cardiomyopathy [65], and hemophilia B [66]. Given the broad range of phenotypes and disease severity in NS-CL/P, we believe that synonymous variants should not be quickly dismissed but rather more actively pursued, as they could underlie important variation in gene expression throughout the development of the face. Our studies also highlight the necessity to functionally test the impact of any such variants, whether non-synonymous or synonymous, on splicing, as highlighted by studies on BRCA1 [67].
The family described in this study lives in a rural part of Honduras, and we have been unable to make contact with them, despite attempts over the last five years, to inquire whether we might be able to follow up the study. While we would have preferred to collect saliva samples from all family members presenting to the initial Operación Sonrisa screening, several were minors accompanied by non-guardian family members, and thus could not be properly consented. The lack of sequencing data from additional family members is a limitation to this study. Nevertheless, we were able to confirm the variant in all affected individuals and obligate carriers on both sides of the pedigree for whom samples were available. Furthermore, sequencing and confirmation of the variant in unaffected individuals, while interesting, would not necessarily contradict our findings, since the phenotype is incompletely penetrant.
Practical difficulties notwithstanding, if the family members can be re-assembled and consented for further research, direct analysis of their IRF6 transcript sequences and expression levels would provide valuable validation of our in vitro findings. In the future, it may be possible to create an animal model of the variant described in this study. While animal model studies are a crucial step toward discovering underlying mechanisms of disease treatments, few animal models of synonymous change-driven disease have been reported. In a number of studies of spinal muscular atrophy animal models, in which a synonymous mutation was associated with the loss of an ESE and the antagonistic gain of a silencer, the in vivo administration of antisense oligonucleotides that mask the silencer vastly improved symptoms [42,50,51,52]. This study highlights the future potential of similar approaches in utero, in conjunction with prenatal genetic screening, for ensuring proper development and preventing synonymous change-driven NS-CL/P.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/8/903/s1, Figure S1: Alignment of the orthologous IRF6 sequence from vertebrate species

Author Contributions

The contributions of each of the authors is as listed: Conceptualization, P.A.S.-L., F.B., M.G., and W.M.III; Patient recruitment and clinical assessment, F.B., M.G., and P.A.S.-L.; Sequencing, data curation and segregation analysis, B.T.; Validation, A.S., B.S., and R.L.M.; Data Analysis, R.L.M., B.S., H.G., A.S., and T.C.C.; Resources, A.A., W.M.III; Writing—Original Draft Preparation, B.S., F.B., R.L.M., P.A.S.-L., A.S., and T.C.C.; Writing—Review and Editing, B.S., A.A., P.A.S.-L., A.S., R.L.M., and T.C.C.; Supervision, F.B., A.A., and W.M.III; Project Administration, A.A.; Funding Acquisition, W.M.III, R.L.M., P.A.S.-L., and T.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

Support for this research was provided by Operation Smile through its International Family Study with private funding and grants from the Marguerite Foundation, the Sorenson Legacy Foundation, the Moore family, and the Rallis Foundation, as well as by funding from the National Institute of Dental and Craniofacial Research’s FaceBase Consortium (U01 DE024443) to R.L.M. and P.A.S.-L. and from the UMKC School of Dentistry recruitment funds to T.C.C.

Acknowledgments

We thank this family for allowing us to learn from them and for providing consent for presentation and publication. We would also like to thank additional members of Operation Smile, Inc (Kathleen Magee, Melissa DiBona, and Lili Arakaki) and Operation Smile Honduras (Jeanie Barjum, Yanuario Garcia, and Seydi Sanchez) for their valuable assistance.

Conflicts of Interest

M.G. and W.M. are affiliated with Operation Smile and/or its associated foundations who have provided some of the financial support for this research and originally flagged the family for research. The authors declare no other conflicts of interest.

References

  1. Dixon, M.J.; Marazita, M.L.; Beaty, T.H.; Murray, J.C. Cleft lip and palate: Understanding genetic and environmental influences. Nat. Rev. Genet. 2011, 12, 167–178. [Google Scholar] [CrossRef] [Green Version]
  2. Fraser, F.C. Thoughts on the etiology of clefts of the palate and lip. Acta Genet. Stat. Med. 1955, 5, 358–369. [Google Scholar] [CrossRef] [PubMed]
  3. Kondo, S.; Schutte, B.C.; Richardson, R.J.; Bjork, B.C.; Knight, A.S.; Watanabe, Y.; Howard, E.; de Lima, R.L.; Daack-Hirsch, S.; Sander, A.; et al. Mutations in IRF6 cause Van der Woude and popliteal pterygium syndromes. Nat. Genet. 2002, 32, 285–289. [Google Scholar] [CrossRef] [PubMed]
  4. Celli, J.; Duijf, P.; Hamel, B.C.; Bamshad, M.; Kramer, B.; Smits, A.P.; Newbury-Ecob, R.; Hennekam, R.C.; Van Buggenhout, G.; van Haeringen, A.; et al. Heterozygous germline mutations in the p53 homolog p63 are the cause of EEC syndrome. Cell 1999, 99, 143–153. [Google Scholar] [CrossRef] [Green Version]
  5. McGrath, J.A.; Duijf, P.H.; Doetsch, V.; Irvine, A.D.; de Waal, R.; Vanmolkot, K.R.; Wessagowit, V.; Kelly, A.; Atherton, D.J.; Griffiths, W.A.; et al. Hay-Wells syndrome is caused by heterozygous missense mutations in the SAM domain of p63. Hum. Mol. Genet. 2001, 10, 221–229. [Google Scholar] [CrossRef] [Green Version]
  6. Dode, C.; Levilliers, J.; Dupont, J.M.; De Paepe, A.; Le Du, N.; Soussi-Yanicostas, N.; Coimbra, R.S.; Delmaghani, S.; Compain-Nouaille, S.; Baverel, F.; et al. Loss-of-function mutations in FGFR1 cause autosomal dominant Kallmann syndrome. Nat. Genet. 2003, 33, 463–465. [Google Scholar] [CrossRef] [Green Version]
  7. Van den Boogaard, M.J.; Dorland, M.; Beemer, F.A.; van Amstel, H.K. MSX1 mutation is associated with orofacial clefting and tooth agenesis in humans. Nat. Genet. 2000, 24, 342–343. [Google Scholar] [CrossRef]
  8. Neiswanger, K.; Deleyiannis, F.W.; Avila, J.R.; Cooper, M.E.; Brandon, C.A.; Vieira, A.R.; Noorchashm, N.; Weinberg, S.M.; Bardi, K.M.; Murray, J.C.; et al. Candidate genes for oral-facial clefts in Guatemalan families. Ann. Plast. Surg. 2006, 56, 518–521. [Google Scholar] [CrossRef]
  9. Neiswanger, K.; Walker, K.; Klotz, C.M.; Cooper, M.E.; Bardi, K.M.; Brandon, C.A.; Weinberg, S.M.; Vieira, A.R.; Martin, R.A.; Czeizel, A.E.; et al. Whorl patterns on the lower lip are associated with nonsyndromic cleft lip with or without cleft palate. Am. J. Med. Genet. Part A 2009, 149a, 2673–2679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Leslie, E.J.; Marazita, M.L. Genetics of cleft lip and cleft palate. Am. J. Med. Genet. C Semin Med. Genet. 2013, 163, 246–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Boyce, J.O.; Raj, S.; Sanchez, K.; Marazita, M.L.; Morgan, A.T.; Kilpatrick, N. Speech Phenotyping in Unaffected Family Members of Individuals With Nonsyndromic Cleft Lip With or Without Palate. Cleft Palate Craniofac J. 2019, 56, 867–876. [Google Scholar] [CrossRef] [PubMed]
  12. Seto-Salvia, N.; Stanier, P. Genetics of cleft lip and/or cleft palate: Association with other common anomalies. Eur. J. med. Genet. 2014, 57, 381–393. [Google Scholar] [CrossRef] [PubMed]
  13. Ghassibe, M.; Bayet, B.; Revencu, N.; Verellen-Dumoulin, C.; Gillerot, Y.; Vanwijck, R.; Vikkula, M. Interferon regulatory factor-6: A gene predisposing to isolated cleft lip with or without cleft palate in the Belgian population. Eur. J. Hum. Genet. EJHG 2005, 13, 1239–1242. [Google Scholar] [CrossRef] [PubMed]
  14. Park, J.W.; McIntosh, I.; Hetmanski, J.B.; Jabs, E.W.; Vander Kolk, C.A.; Wu-Chou, Y.H.; Chen, P.K.; Chong, S.S.; Yeow, V.; Jee, S.H.; et al. Association between IRF6 and nonsyndromic cleft lip with or without cleft palate in four populations. Genet. Med. Off. J. Am. Coll. Med. Genet. 2007, 9, 219–227. [Google Scholar] [CrossRef] [Green Version]
  15. Rahimov, F.; Marazita, M.L.; Visel, A.; Cooper, M.E.; Hitchler, M.J.; Rubini, M.; Domann, F.E.; Govil, M.; Christensen, K.; Bille, C.; et al. Disruption of an AP-2alpha binding site in an IRF6 enhancer is associated with cleft lip. Nat. Genet. 2008, 40, 1341–1347. [Google Scholar] [CrossRef] [Green Version]
  16. Scapoli, L.; Palmieri, A.; Martinelli, M.; Pezzetti, F.; Carinci, P.; Tognon, M.; Carinci, F. Strong evidence of linkage disequilibrium between polymorphisms at the IRF6 locus and nonsyndromic cleft lip with or without cleft palate, in an Italian population. Am. J. Hum. Genet. 2005, 76, 180–183. [Google Scholar] [CrossRef] [Green Version]
  17. Zucchero, T.M.; Cooper, M.E.; Maher, B.S.; Daack-Hirsch, S.; Nepomuceno, B.; Ribeiro, L.; Caprau, D.; Christensen, K.; Suzuki, Y.; Machida, J.; et al. Interferon regulatory factor 6 (IRF6) gene variants and the risk of isolated cleft lip or palate. N. Engl. J. Med. 2004, 351, 769–780. [Google Scholar] [CrossRef]
  18. Ingraham, C.R.; Kinoshita, A.; Kondo, S.; Yang, B.; Sajan, S.; Trout, K.J.; Malik, M.I.; Dunnwald, M.; Goudy, S.L.; Lovett, M.; et al. Abnormal skin, limb and craniofacial morphogenesis in mice deficient for interferon regulatory factor 6 (Irf6). Nat. Genet. 2006, 38, 1335–1340. [Google Scholar] [CrossRef] [Green Version]
  19. Iwata, J.; Suzuki, A.; Pelikan, R.C.; Ho, T.V.; Sanchez-Lara, P.A.; Urata, M.; Dixon, M.J.; Chai, Y. Smad4-Irf6 genetic interaction and TGFbeta-mediated IRF6 signaling cascade are crucial for palatal fusion in mice. Development (Camb. Engl.) 2013, 140, 1220–1230. [Google Scholar] [CrossRef] [Green Version]
  20. Richardson, R.J.; Dixon, J.; Malhotra, S.; Hardman, M.J.; Knowles, L.; Boot-Handford, R.P.; Shore, P.; Whitmarsh, A.; Dixon, M.J. Irf6 is a key determinant of the keratinocyte proliferation-differentiation switch. Nat. Genet. 2006, 38, 1329–1334. [Google Scholar] [CrossRef]
  21. Richardson, R.J.; Dixon, J.; Jiang, R.; Dixon, M.J. Integration of IRF6 and Jagged2 signalling is essential for controlling palatal adhesion and fusion competence. Hum. Mol. Genet. 2009, 18, 2632–2642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Parada-Sanchez, M.T.; Chu, E.Y.; Cox, L.L.; Undurty, S.S.; Standley, J.M.; Murray, J.C.; Cox, T.C. Disrupted IRF6-NME1/2 Complexes as a Cause of Cleft Lip/Palate. J. Dent. Res. 2017, 96, 1330–1338. [Google Scholar] [CrossRef] [PubMed]
  23. Desmet, F.O.; Hamroun, D.; Lalande, M.; Collod-Béroud, G.; Claustres, M.; Béroud, C. Human Splicing Finder: An online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009, 37, e67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Quinodoz, M.; Royer-Bertrand, B.; Cisarova, K.; Di Gioia, S.A.; Superti-Furga, A.; Rivolta, C. DOMINO: Using Machine Learning to Predict Genes Associated with Dominant Disorders. Am. J. Hum. Genet. 2017, 101, 623–629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Rentzsch, P.; Witten, D.; Cooper, G.M.; Shendure, J.; Kircher, M. CADD: Predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019, 47, D886–D894. [Google Scholar] [CrossRef]
  26. Bretschneider, H.; Gandhi, S.; Deshwar, A.G.; Zuberi, K.; Frey, B.J. COSSMO: Predicting competitive alternative splice site selection using deep learning. Bioinformatics 2018, 34, i429–i437. [Google Scholar] [CrossRef]
  27. Anczukow, O.; Buisson, M.; Salles, M.J.; Triboulet, S.; Longy, M.; Lidereau, R.; Sinilnikova, O.M.; Mazoyer, S. Unclassified variants identified in BRCA1 exon 11: Consequences on splicing. Genes Chromosomes Cancer 2008, 47, 418–426. [Google Scholar] [CrossRef]
  28. Diederichs, S.; Bartsch, L.; Berkmann, J.C.; Frose, K.; Heitmann, J.; Hoppe, C.; Iggena, D.; Jazmati, D.; Karschnia, P.; Linsenmeier, M.; et al. The dark matter of the cancer genome: Aberrations in regulatory elements, untranslated regions, splice sites, non-coding RNA and synonymous mutations. EMBO Mol. Med. 2016, 8, 442–457. [Google Scholar] [CrossRef]
  29. Gartner, J.J.; Parker, S.C.; Prickett, T.D.; Dutton-Regester, K.; Stitzel, M.L.; Lin, J.C.; Davis, S.; Simhadri, V.L.; Jha, S.; Katagiri, N.; et al. Whole-genome sequencing identifies a recurrent functional synonymous mutation in melanoma. Proc. Natl. Acad. Sci. USA 2013, 110, 13481–13486. [Google Scholar] [CrossRef] [Green Version]
  30. Hansen, T.V.; Steffensen, A.Y.; Jonson, L.; Andersen, M.K.; Ejlertsen, B.; Nielsen, F.C. The silent mutation nucleotide 744 G --> A, Lys172Lys, in exon 6 of BRCA2 results in exon skipping. Breast Cancer Res. Treat. 2010, 119, 547–550. [Google Scholar] [CrossRef] [Green Version]
  31. Kimchi-Sarfaty, C.; Oh, J.M.; Kim, I.W.; Sauna, Z.E.; Calcagno, A.M.; Ambudkar, S.V.; Gottesman, M.M. A “silent” polymorphism in the MDR1 gene changes substrate specificity. Science (New York N.Y.) 2007, 315, 525–528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Montera, M.; Piaggio, F.; Marchese, C.; Gismondi, V.; Stella, A.; Resta, N.; Varesco, L.; Guanti, G.; Mareni, C. A silent mutation in exon 14 of the APC gene is associated with exon skipping in a FAP family. J. Med. Genet. 2001, 38, 863–867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Pecina-Slaus, N.; Majic, Z.; Musani, V.; Zeljko, M.; Cupic, H. Report on mutation in exon 15 of the APC gene in a case of brain metastasis. J. Neuro-Oncol. 2010, 97, 143–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Raponi, M.; Kralovicova, J.; Copson, E.; Divina, P.; Eccles, D.; Johnson, P.; Baralle, D.; Vorechovsky, I. Prediction of single-nucleotide substitutions that result in exon skipping: Identification of a splicing silencer in BRCA1 exon 6. Hum. Mutat. 2011, 32, 436–444. [Google Scholar] [CrossRef] [Green Version]
  35. Supek, F.; Minana, B.; Valcarcel, J.; Gabaldon, T.; Lehner, B. Synonymous mutations frequently act as driver mutations in human cancers. Cell 2014, 156, 1324–1335. [Google Scholar] [CrossRef] [Green Version]
  36. Toscano, C.; Raimundo, S.; Klein, K.; Eichelbaum, M.; Schwab, M.; Zanger, U.M. A silent mutation (2939G>A, exon 6; CYP2D6*59) leading to impaired expression and function of CYP2D6. Pharmacogn. Genom. 2006, 16, 767–770. [Google Scholar] [CrossRef]
  37. Bartoszewski, R.A.; Jablonsky, M.; Bartoszewska, S.; Stevenson, L.; Dai, Q.; Kappes, J.; Collawn, J.F.; Bebok, Z. A synonymous single nucleotide polymorphism in DeltaF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein. J. Biolog. Chem. 2010, 285, 28741–28748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Lazrak, A.; Fu, L.; Bali, V.; Bartoszewski, R.; Rab, A.; Havasi, V.; Keiles, S.; Kappes, J.; Kumar, R.; Lefkowitz, E.; et al. The silent codon change I507-ATC->ATT contributes to the severity of the DeltaF508 CFTR channel dysfunction. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2013, 27, 4630–4645. [Google Scholar] [CrossRef] [Green Version]
  39. Brest, P.; Lapaquette, P.; Souidi, M.; Lebrigand, K.; Cesaro, A.; Vouret-Craviari, V.; Mari, B.; Barbry, P.; Mosnier, J.F.; Hebuterne, X.; et al. A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn’s disease. Nat. Genet. 2011, 43, 242–245. [Google Scholar] [CrossRef] [PubMed]
  40. Oh, S.H.; Baek, J.; Liany, H.; Foo, J.N.; Kim, K.M.; Yang, S.C.; Liu, J.; Song, K. A Synonymous Variant in IL10RA Affects RNA Splicing in Paediatric Patients with Refractory Inflammatory Bowel Disease. J. Crohns Colitis 2016, 10, 1366–1371. [Google Scholar] [CrossRef] [Green Version]
  41. Nackley, A.G.; Shabalina, S.A.; Tchivileva, I.E.; Satterfield, K.; Korchynskyi, O.; Makarov, S.S.; Maixner, W.; Diatchenko, L. Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science (New York, N.Y.) 2006, 314, 1930–1933. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Sterne-Weiler, T.; Sanford, J.R. Exon identity crisis: Disease-causing mutations that disrupt the splicing code. Genome Biol. 2014, 15, 201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Richard, P.; Gaudon, K.; Fournier, E.; Jackson, C.; Bauche, S.; Haddad, H.; Koenig, J.; Echenne, B.; Hantai, D.; Eymard, B. A synonymous CHRNE mutation responsible for an aberrant splicing leading to congenital myasthenic syndrome. Neuromuscul. Disord. NMD 2007, 17, 409–414. [Google Scholar] [CrossRef]
  44. Vidal, C.; Cachia, A.; Xuereb-Anastasi, A. Effects of a synonymous variant in exon 9 of the CD44 gene on pre-mRNA splicing in a family with osteoporosis. Bone 2009, 45, 736–742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Parmley, J.L.; Chamary, J.V.; Hurst, L.D. Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers. Mol. Biol. Evolut. 2006, 23, 301–309. [Google Scholar] [CrossRef] [Green Version]
  46. Chamary, J.V.; Parmley, J.L.; Hurst, L.D. Hearing silence: Non-neutral evolution at synonymous sites in mammals. Nat. Rev. Genet. 2006, 7, 98–108. [Google Scholar] [CrossRef]
  47. Savisaar, R.; Hurst, L.D. Exonic splice regulation imposes strong selection at synonymous sites. Genome Res. 2018, 28, 1442–1454. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, Z.; Nielsen, R. Mutation-selection models of codon substitution and their use to estimate selective strengths on codon usage. Mol. Biol. Evolut. 2008, 25, 568–579. [Google Scholar] [CrossRef] [Green Version]
  49. Drummond, D.A.; Wilke, C.O. Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution. Cell 2008, 134, 341–352. [Google Scholar] [CrossRef] [Green Version]
  50. Hua, Y.; Sahashi, K.; Hung, G.; Rigo, F.; Passini, M.A.; Bennett, C.F.; Krainer, A.R. Antisense correction of SMN2 splicing in the CNS rescues necrosis in a type III SMA mouse model. Genes Dev. 2010, 24, 1634–1644. [Google Scholar] [CrossRef] [Green Version]
  51. Hua, Y.; Sahashi, K.; Rigo, F.; Hung, G.; Horev, G.; Bennett, C.F.; Krainer, A.R. Peripheral SMN restoration is essential for long-term rescue of a severe spinal muscular atrophy mouse model. Nature 2011, 478, 123–126. [Google Scholar] [CrossRef] [PubMed]
  52. Passini, M.A.; Bu, J.; Richards, A.M.; Kinnecom, C.; Sardi, S.P.; Stanek, L.M.; Hua, Y.; Rigo, F.; Matson, J.; Hung, G.; et al. Antisense oligonucleotides delivered to the mouse CNS ameliorate symptoms of severe spinal muscular atrophy. Sci. Transl. Med. 2011, 3, 72ra18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Abramowicz, A.; Gos, M. Splicing mutations in human genetic disorders: Examples, detection, and confirmation. J. Appl. Genet. 2018, 59, 253–268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lindeboom, R.G.; Supek, F.; Lehner, B. The rules and impact of nonsense-mediated mRNA decay in human cancers. Nat. Genet. 2016, 48, 1112–1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Glisovic, T.; Bachorik, J.L.; Yong, J.; Dreyfuss, G. RNA-binding proteins and post-transcriptional gene regulation. FEBS Lett. 2008, 582, 1977–1986. [Google Scholar] [CrossRef] [Green Version]
  56. Zahler, A.M.; Damgaard, C.K.; Kjems, J.; Caputi, M. SC35 and heterogeneous nuclear ribonucleoprotein A/B proteins bind to a juxtaposed exonic splicing enhancer/exonic splicing silencer element to regulate HIV-1 tat exon 2 splicing. J. Biolog. Chem. 2004, 279, 10077–10084. [Google Scholar] [CrossRef] [Green Version]
  57. Letra, A.; Menezes, R.; Cooper, M.E.; Fonseca, R.F.; Tropp, S.; Govil, M.; Granjeiro, J.M.; Imoehl, S.R.; Mansilla, M.A.; Murray, J.C.; et al. CRISPLD2 variants including a C471T silent mutation may contribute to nonsyndromic cleft lip with or without cleft palate. Cleft Palate Craniofac J. 2011, 48, 363–370. [Google Scholar] [CrossRef] [Green Version]
  58. Kumari, P.; Singh, S.K.; Raman, R. TGFbeta3, MSX1, and MMP3 as Candidates for NSCL+/-P in an Indian Population. Cleft Palate Craniofac J. 2019, 56, 363–372. [Google Scholar] [CrossRef]
  59. Gaczkowska, A.; Biedziak, B.; Budner, M.; Zadurska, M.; Lasota, A.; Hozyasz, K.K.; Dabrowska, J.; Wojcicki, P.; Szponar-Zurowska, A.; Zukowski, K.; et al. PAX7 nucleotide variants and the risk of non-syndromic orofacial clefts in the Polish population. Oral. Dis. 2019, 25, 1608–1618. [Google Scholar] [CrossRef]
  60. Ke, S.; Shang, S.; Kalachikov, S.M.; Morozova, I.; Yu, L.; Russo, J.J.; Ju, J.; Chasin, L.A. Quantitative evaluation of all hexamers as exonic splicing elements. Genome Res. 2011, 21, 1360–1374. [Google Scholar] [CrossRef] [Green Version]
  61. Sharma, Y.; Miladi, M.; Dukare, S.; Boulay, K.; Caudron-Herger, M.; Groß, M.; Backofen, R.; Diederichs, S. A pan-cancer analysis of synonymous mutations. Nat. Commun. 2019, 10, 2569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Xu, W.; Yang, X.; Hu, X.; Li, S. Fifty-four novel mutations in the NF1 gene and integrated analyses of the mutations that modulate splicing. Int. J. Mol. Med. 2014, 34, 53–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Yamashita, M.; Wakatsuki, R.; Kato, T.; Okano, T.; Yamanishi, S.; Mayumi, N.; Tanaka, M.; Ogura, Y.; Kanegane, H.; Nonoyama, S.; et al. A synonymous splice site mutation in IL2RG gene causes late-onset combined immunodeficiency. Int. J. Hematol. 2019, 109, 603–611. [Google Scholar] [CrossRef] [PubMed]
  64. Platt, C.D.; Massaad, M.J.; Cangemi, B.; Schmidt, B.; Aldhekri, H.; Geha, R.S. Janus kinase 3 deficiency caused by a homozygous synonymous exonic mutation that creates a dominant splice site. J. Allergy Clin. Immunol. 2017, 140, 268–271.e266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Ito, K.; Patel, P.N.; Gorham, J.M.; McDonough, B.; DePalma, S.R.; Adler, E.E.; Lam, L.; MacRae, C.A.; Mohiuddin, S.M.; Fatkin, D.; et al. Identification of pathogenic gene mutations in LMNA and MYBPC3 that alter RNA splicing. Proc. Natl. Acad. Sci. USA 2017, 114, 7689–7694. [Google Scholar] [CrossRef] [Green Version]
  66. Odaira, K.; Tamura, S.; Suzuki, N.; Kakihara, M.; Hattori, Y.; Tokoro, M.; Suzuki, S.; Takagi, A.; Katsumi, A.; Hayakawa, F.; et al. Apparent synonymous mutation F9 c.87A>G causes secretion failure by in-frame mutation with aberrant splicing. Thromb. Res. 2019, 179, 95–103. [Google Scholar] [CrossRef]
  67. Thomassen, M.; Blanco, A.; Montagna, M.; Hansen, T.V.; Pedersen, I.S.; Gutiérrez-Enríquez, S.; Menéndez, M.; Fachal, L.; Santamariña, M.; Steffensen, A.Y.; et al. Characterization of BRCA1 and BRCA2 splicing variants: A collaborative report by ENIGMA consortium members. Breast Cancer Res. Treat. 2012, 132, 1009–1023. [Google Scholar] [CrossRef]
Genes 11 00903 g001 550
Figure 1. Pedigree of Honduran family with non-syndromic cleft lip and/or palate (NS-CL/P). Pedigree shows inheritance from individual I-2 through both sides of the pedigree: II-1, III-3, III-6 & II-3, III-7, II-4, III-11, III-12, and II-10. Arrow marks the proband (individual III-6).
Figure 1. Pedigree of Honduran family with non-syndromic cleft lip and/or palate (NS-CL/P). Pedigree shows inheritance from individual I-2 through both sides of the pedigree: II-1, III-3, III-6 & II-3, III-7, II-4, III-11, III-12, and II-10. Arrow marks the proband (individual III-6).
Genes 11 00903 g001
Genes 11 00903 g002 550
Figure 2. The synonymous IRF6 c.921C>T variant impacts splicing. (A) In silico analysis predicted the creation of an exonic splice silencer (ESS) element (score: 71) and enhanced putative recognition sites for the splicing factors, including SRSF1 and SRSF2. The scores are indicated in brackets. (Bi) In vitro minigene splice assay detected multiple abnormal splice products (bands 2–4) in addition to the normal-sized product (band 1). Note: the slower migrating bands were seen with both the wildtype and mutant and are thought to represent pre-spliced forms, as they were absent from control PCRs using the RNA (without reverse transcriptase) as a template. (Bii) Quantification of splice products revealed that abnormal products constituted ~67% of allelic transcripts. (C) Sequencing of purified splice products resolved the nature of the underlying splice events. The arrowheads depict the location of nucleotide position c.921: C = green; T = purple. The variant unmasks the cryptic exonic splice site at position c.915, as well as a site within intron 6 (135 nucleotide 5′ of the canonical 5′ splice acceptor of exon 7). The canonical acceptor site is used differentially in bands 2 and 3.
Figure 2. The synonymous IRF6 c.921C>T variant impacts splicing. (A) In silico analysis predicted the creation of an exonic splice silencer (ESS) element (score: 71) and enhanced putative recognition sites for the splicing factors, including SRSF1 and SRSF2. The scores are indicated in brackets. (Bi) In vitro minigene splice assay detected multiple abnormal splice products (bands 2–4) in addition to the normal-sized product (band 1). Note: the slower migrating bands were seen with both the wildtype and mutant and are thought to represent pre-spliced forms, as they were absent from control PCRs using the RNA (without reverse transcriptase) as a template. (Bii) Quantification of splice products revealed that abnormal products constituted ~67% of allelic transcripts. (C) Sequencing of purified splice products resolved the nature of the underlying splice events. The arrowheads depict the location of nucleotide position c.921: C = green; T = purple. The variant unmasks the cryptic exonic splice site at position c.915, as well as a site within intron 6 (135 nucleotide 5′ of the canonical 5′ splice acceptor of exon 7). The canonical acceptor site is used differentially in bands 2 and 3.
Genes 11 00903 g002
Genes 11 00903 g003 550
Figure 3. The splice regulatory factors impact the extent of abnormal splicing caused by the synonymous IRF6 c.921C>T variant. The total amount of PCR splice products (WT: band 1, Mut: bands 2–4, as per Figure 2B) were quantified and expressed as a percentage of the total. Six replicates were performed. Statistical analysis was carried out using two-tailed Student’s t-tests; ** p < 0.05. Actual p-values from the tests are presented within the main text.
Figure 3. The splice regulatory factors impact the extent of abnormal splicing caused by the synonymous IRF6 c.921C>T variant. The total amount of PCR splice products (WT: band 1, Mut: bands 2–4, as per Figure 2B) were quantified and expressed as a percentage of the total. Six replicates were performed. Statistical analysis was carried out using two-tailed Student’s t-tests; ** p < 0.05. Actual p-values from the tests are presented within the main text.
Genes 11 00903 g003
Table
Table 1. Shared variants identified by exome sequencing. From the original list of 411 shared variants, further filtering based on DOMINO class consistent with a dominant mode of inheritance and minor allele frequency of less than 1 × 10−4. (VLD = very likely dominant; LD = likely dominant; EDR = either dominant or recessive). Note: Allele count and minor allele frequency (MAF) were determined from the gnomAD database (June, 2020).
Table 1. Shared variants identified by exome sequencing. From the original list of 411 shared variants, further filtering based on DOMINO class consistent with a dominant mode of inheritance and minor allele frequency of less than 1 × 10−4. (VLD = very likely dominant; LD = likely dominant; EDR = either dominant or recessive). Note: Allele count and minor allele frequency (MAF) were determined from the gnomAD database (June, 2020).
Genomic CoordinateGene NameLocationTypeProtein ChangeCADD ScaledDOMINO Class (Probability)gnomAD
Allele Count (MAF)		Associated Disease
chr8:g.13251073G>CDLC1CodingMissensep.Pro435Ala23.4LD (0.706)0Nephrosis
chr1:g.155838558G>ASYT11CodingSynonymousp.Arg279Arg11.89EDR (0.478)0
chr1:g.114940211G>TTRIM333’UTRSNP 0.742VLD (0.997)0
chr3:g.149684103G>APFN23’UTRSNP 16.21VLD (0.948)0
chr1:g.16262494_16262499dupTGTCCCSPENCodingDuplicationp.Val3254_Pro3255dup VLD (1)0
chr5:g.149922536G>ANDST1IntronicSplice 14.37VLD (0.964)0Mental retardation
chr11:g.156907027_156907033delCCTGCTTARHGEF11IntronicDeletion LD (0.794)0
chr7:g.137597556G>CCREB3L2IntronicSplice 5.642LD (0.616)0Fibromyxoid sarcoma
chr8:g.38947723delTADAM9IntronicDeletion EDR (0.584)0Cone-rod dystrophy
chr18:g.59936640dupCKIAA1468IntronicSplice VLD (0.959)0
chr12:g.64882264delTTBK1IntronicSplice VLD (0.882)0Frontotemporal dementia/ALS
chr9:g.73458046dupATRPM3IntronicSplice EDR (0.585)0Intellectual disability/epilepsy
chr19:g.10798342G>AILF3CodingMissensep.Gly798Arg24.2LD (0.677)2 / 247274 (8.09 × 10−6)
chr1:g.209963979G>AIRF6CodingSynonymousp.Ser307Ser14.48LD (0.733)3 / 251352 (1.19 × 10−5)Orofacial clefting
chr8:g.66647047A>GPDE7ACodingSynonymousp.Cys226Cys12.85EDR (0.407)3 / 250766 (1.20 × 10−5)
chr12:g.6094200G>AVWFIntronicSNP 5.659LD (0.617)6 / 251458 (2.39 × 10−5)von Willebrand disease
chr2:g.240085511C>THDAC4CodingMissensep.Arg200His25VLD (0.999)7 / 251428 (2.78 × 10−5)Brachydactyly, mental retardation
chr8:g.26484831G>ADPYSL2IntronicSNP 0.659VLD (0.995)9 / 282698 (3.18 × 10−5)
chr1:g.16268440_16268443delTAAASPEN/ZBTB173’UTRDeletion VLD (1)7 / 218068 (3.21 × 10−5)
chr7:g.98555591G>TTRRAPIntronicSNP 0.096VLD (1)10 / 248972 (4.02 × 10−5)Developmental delay/autism, dysmorphic facies
chr1:g.150551341G>AMCL1CodingSynonymousp.Arg222Arg12.07VLD (0.997)13 / 281774 (4.61 × 10−5)
chr17:g.78113791C>AEIF4A3IntronicSNP 0.049VLD (0.999)12 / 244608 (4.91 × 10−5)Robin sequence, cleft mandible, limb anomalies
chr1:g.154527884G>AUBE2Q1IntronicSNP 4.602LD (0.727)13 / 251300 (5.17 × 10−5)
chr1:g.155948197T>GARHGEF2CodingMissensep.Thr8Pro22.4VLD (0.916)14 / 167764 (8.35 × 10−5)Neurodevelopmental disorder, brain malformations
chr10:g.33196031T>CITGB1codingMissensep.Lys791Arg22.6VLD (0.996)26 / 279598 (9.30 × 10−5)
chr12:g.52911558C>TKRT5intronicSNP 5.636EDR (0.589)28 / 282824 (9.90 × 10−5)epidermolysis bullosa
chr17:g.72759561G>ASLC9A3R1codingMissensep.Arg220Lys14.76VLD (0.865)25 / 250294 (9.99 × 10−5)Nephrolithiasis/osteoporosis, hypophosphatemia

Share and Cite

MDPI and ACS Style

Sylvester, B.; Brindopke, F.; Suzuki, A.; Giron, M.; Auslander, A.; Maas, R.L.; Tsai, B.; Gao, H.; Magee, W., III; Cox, T.C.; et al. A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate. Genes 2020, 11, 903. https://doi.org/10.3390/genes11080903

AMA Style

Sylvester B, Brindopke F, Suzuki A, Giron M, Auslander A, Maas RL, Tsai B, Gao H, Magee W III, Cox TC, et al. A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate. Genes. 2020; 11(8):903. https://doi.org/10.3390/genes11080903

Chicago/Turabian Style

Sylvester, Beau, Frederick Brindopke, Akiko Suzuki, Melissa Giron, Allyn Auslander, Richard L. Maas, Becky Tsai, Hanlin Gao, William Magee, III, Timothy C. Cox, and et al. 2020. "A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate" Genes 11, no. 8: 903. https://doi.org/10.3390/genes11080903

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop