TLE6 mutation causes the earliest known human embryonic lethality
- Anas M. Alazami†1,
- Salma M. Awad†2,
- Serdar Coskun†3,
- Saad Al-Hassan4,
- Hadia Hijazi2,
- Firdous M. Abdulwahab1,
- Coralie Poizat2Email author and
- Fowzan S. Alkuraya1, 5Email authorView ORCID ID profile
© Alazami et al. 2015
Received: 18 August 2015
Accepted: 28 September 2015
Published: 5 November 2015
Embryonic lethality is a recognized phenotypic expression of individual gene mutations in model organisms. However, identifying embryonic lethal genes in humans is challenging, especially when the phenotype is manifested at the preimplantation stage.
In an ongoing effort to exploit the highly consanguineous nature of the Saudi population to catalog recessively acting embryonic lethal genes in humans, we have identified two families with a female-limited infertility phenotype. Using autozygosity mapping and whole exome sequencing, we map this phenotype to a single mutation in TLE6, a maternal effect gene that encodes a member of the subcortical maternal complex in mammalian oocytes. Consistent with the published phenotype of mouse Tle6 mutants, embryos from female patients who are homozygous for the TLE6 mutation fail to undergo early cleavage, with resulting sterility. The human mutation abrogates TLE6 phosphorylation, a step that is reported to be critical for the PKA-mediated progression of oocyte meiosis II. Furthermore, the TLE6 mutation impairs its binding to components of the subcortical maternal complex.
In this first report of a human defect in a member of the subcortical maternal subcritical maternal complex, we show that the TLE6 mutation is gender-specific and leads to the earliest known human embryonic lethality phenotype.
Thousands of genes are known to exert specific phenotypic effects when mutated individually and are referred to as ‘Mendelian’. The phenotypic consequences of mutations in Mendelian genes involve virtually every known aspect of human physiology and anatomy . Embryonic lethality sometimes represents the severe end of the phenotypic expression of Mendelian genes that are otherwise better known for their postnatal phenotypic consequences . On the other hand, embryonic lethality may be the only phenotypic expression of genes that play a fundamental role in early development . Identifying these genes, therefore, can provide novel insights into basic biological processes during development [1, 4]. The fraction of Mendelian genes in humans that are ‘embryonic lethal’ is unknown but extrapolation from model organisms, for example, 19 % of yeast genes and 30 % of mouse genes are embryonic lethal when knocked out, suggests that they are at least in the hundreds [5, 6].
Consanguineous families facilitate the occurrence of autosomal recessive Mendelian diseases by virtue of rendering founder mutations homozygous through autozygosity, and have greatly aided the discovery of novel Mendelian genes through autozygosity mapping [7, 8]. More recently, combining autozygosity mapping with genomic sequencing has markedly accelerated novel disease gene discovery . Due to practical reasons, the entry phenotypes in these studies is almost always a postnatal recognizable pattern rather than embryonic lethality. We have recently shown that consanguineous families that segregate embryonic lethal phenotypes provide an exceptionally high yield for novel gene discovery where seven novel genes were identified from a study of less than 20 families . However, embryonic lethality in that study was limited to late first trimester and beyond, such that very early embryonic lethal events could not be assessed. In this study, we report the identification of the genetic cause of the earliest known embryonic lethal phenotype through the study of consanguineous families that segregate female-limited infertility and failure of embryonic development beyond the zygote formation.
Materials and methods
Families were recruited through the Reproductive Medicine Unit at KFSHRC based on the observation of abnormal embryo development during regular in vitro fertilization treatment of consanguineous couples. Eligible families and controls were enrolled after signing a KFSHRC IRB-approved written informed consent (RAC #2121053). Venous blood was collected in EDTA and, when possible, in Na-heparin tubes for DNA extraction and lymphoblastoid cell line establishment, respectively. All methods comply with the Helsinki Declaration.
Autozygome and linakge analysis
Determination of the entire set of autozygous intervals per genome (autozygome) was through genomewide SNP genotyping (Axiom SNP chip, Affymetrix) followed by mapping of runs of homozygosity as surrogates of autozygosity using AutoSNPa v4, as described before [10, 11]. Overlap in the autozygome of affected individuals was employed as a strategy to determine the critical disease locus. Statistical confirmation of the critical locus was achieved by linkage analysis using easyLINKAGE .
Whole exome sequencing
Exome capture was performed using TruSeq Exome Enrichment kit (Illumina) following the manufacturer’s protocol. Samples were prepared as an Illumina sequencing library, and in the second step, the sequencing libraries were enriched for the desired target using the Illumina Exome Enrichment protocol. The captured libraries were sequenced using Illumina HiSeq 2000 Sequencer. The reads are mapped against UCSC hg19  by BWA ver.0.5.9rc1 , without unordered sequences and alternate haplotypes. The Picard-tools suite (ver.1.59) was then utilized to sort by mapping coordinates, and BEDtools (ver. 2.15.0) filtered out any reads not present in the targeted exonic regions. SNPs and Indels were detected by SAMTOOLS ver.0.1.18  and annotated using ANNOVAR ver.Nov 2011 . The candidacy of the resulting variants was based on their physical location within the autozygome of the affected individual, their population frequency and predicted effect on the protein as described before . Data used in this paper come from a small and well-defined family. To protect the identity of individuals, these confidential data are not publicly available.
Western blot and phosphorylation analysis
Epstein Barr Virus (EBV) transformed cell lines were produced from three healthy donors (controls) and from three individuals who are homozygous for the TLE6 S510Y mutation (patients, see below). Western blot analysis was performed as described . Briefly, cells were harvested by centrifugation and resuspended in lysis buffer (20 mM Tris pH 7.5, 350 mM NaCl, 0.05 % β-mercaptoethanol) supplemented with a protease inhibitor cocktail. After sonication and centrifugation, 30 μg of total cell lysates were analyzed by SDS-PAGE on 10 % acrylamide or on 8 % Phospho-tag acrylamide gels (Wako, TX, USA), followed by transfer of the proteins onto nitrocellulose membrane. After blocking in 5 % milk in TBS-Tween, the membranes were incubated with anti-TLE6, anti-KDHC3L/Ecat1 from (Abcam, Cambridge, MA, USA) or anti-OOEP, anti-Flag and anti-GAPDH from (Santa Cruz, CA, USA). After washing, secondary reactions were carried out with biotin conjugated secondary antibodies followed by anti-avidin-HRP conjugated antibody. Signals were visualized using an LAS 4000 mini (GE Healthcare, UK) and quantified using ImageQuant software (GE Healthcare, UK).
Phosphatase inhibitor treatment
A total of 30 μg of whole cell lysates from the two control individuals were incubated with 40 μM of calf intestine alkaline phosphatase (CIP) (Promega, Madison, WI, USA) and equimolar amount of PKI (5–24), PKA Inhibitor (Santa Cruz, CA, USA) in a reaction buffer (50 mM Tris–HCl, pH 9.3; 1 mM MgCl2; 0.1 mM ZnCl2; 1 mM spermidine). Extracts were then analyzed by western blot analysis as described before .
Production of TLE6 S510Y protein
Mutant TLE6 S510Y was produced using the Stratagene site-directed mutagenesis system (Stratagene, La Jolla, CA, USA) according to manufacturer’s instructions; using Flag tagged TLE6 plasmid (Origene, Rockville, MD, USA). Primers used to generate the mutation are sense: 5′-TCCTGAGCGTCAAGTTCT(A)CCCCTTTGGCCAGTGGTG-3′ and anti-sense: 5′-AGGACTCGCAGTTCAAGATGGGGAAACCGGTCACCAC-3. The base change was from C to A, shown between brackets. After verification of the mutation by sequencing, control transformed lymphocytes were transfected with wild-type and mutant plasmids using lipofectamine (Invitrogen, Waltham, MA, USA). Twenty-four hours post transfection, total cell lysates were prepared and analyzed by western blotting or immunoprecipitation.
Cell extracts were re-suspended in a buffer containing 50 mM Tris pH7.6, 150 mM or 500 mM NaCl, 1 mM EDTA, 1%Triton X-100, 1 mM PMSF supplemented with phosphatase and protease inhibitor cocktail (Sigma). A total of 60 μg of total protein was incubated with Flag-IgG sepharose beads (Sigma-Aldrich, St. Louis, MO, USA) or Protein A dynabeads (Invitrogen, Waltham, MA, USA) coupled to anti-OOEP antibody (Santa Cruz, CA, USA) or control IgG in a pull-down buffer (50 mM HEPES pH 7.5, 1 mM EDTA, 150 mM NaCl, 10 % glycerol, 0.1 % Tween 20, 0.5 mM DTT, 1 mM PMSF, 2 μg/mL leupeptin, and 2 μg/mL pepstatin A). Extracts were incubated with the beads for at least 2 h at 4 °C while mixing on a rotating wheel. After collection of the supernatants, the beads were washed with pull-down buffer and left as a 50 % slurry after a final wash. After elution, proteins were loaded on 10 % SDS-PAGE gel and immunoblotting was performed using the indicated antibodies.
In vitro phosphorylation
Wild-type and TLE6-S510Y proteins were expressed in HEK293 and immunoprecipitated using anti-Flag beads (Sigma-Aldrich, St. Louis, MO, USA). A total of 20 μg of purified proteins was added to a standard PKA mixture containing 20 mM Hepes (pH 7.5), 5 mM MgCl2, 1 mM unlabeled ATP, 1 mM 1,4-dithiothreitol, 100 mM NaCl, and 1 mM [γ32P]ATP (0.5 Ci/mmol- Perkin Elmer, Waltham, MA, USA), with different amounts of PKA (Promega, Madison, WI, USA) as indicated and incubated for 1 h at 37 °C. Reactions were stopped by adding 20 μL of 2X Laemmli sample buffer and proteins were separated by SDS-PAGE electrophoresis followed by autoradiography and immunoblotting.
Identification of a preimplantation embryonic lethal phenotype
Identification of a novel disease locus for early embryonic lethality defined by TLE6 mutation
Although the phenotype in the three women is primarily failure of zygote formation, the occasional formation of zygotes that did not divide suggested that the phenotype in these patients can be considered preimplantation embryonic lethality. The observation of this phenotype only in certain female members of each family suggests that this phenotype is female-limited, perhaps due to a recessive mutation of a maternal effect gene. To identify the likely causal mutation we performed whole-exome sequencing in two of the three affected women and only considered homozygous, coding/splicing variants within the autozygome of each individual. Only one novel variant remained after applying these filters in both exomes, a homozygous substitution in TLE6 (transducin-like enhancer of split 6, OMIM ID: 612399; NM_001143986.1:c.1529C > A:p.S510Y) (Fig. 1b, d). Reassuringly, this variant was present in an autozygous interval that is shared by the three females (chr19: 2712016–3918047, hg19) (Fig. 1c). The identical haplotype indicated the presence of a common ancestor, consistent with the shared geographic location of the two families. Furthermore, linkage analysis showed that this haplotype is the only haplotype to achieve genomewide significance with LOD of 4.17 (Additional file 1: Figure S1).
This variant was confirmed by Sanger sequencing to be homozygous in the three female patients. Consistent with this mutation exerting a gender-limited phenotype, we also had access to one fertile brother who was also homozygous for this variant (Fig. 1a). This mutation was only observed once in the heterozygous state in 615 in-house Saudi exomes (allele frequency <0.001) and was absent in the 1000 Genomes and Exome Variant Server databases. The affected residue appears to be universally conserved among mammalian TLE6 orthologs (Fig. 1e), and is present within a stretch of seven WD40 domain repeats, which dominate the protein’s C-terminus (Fig. 1f).
TLE6 mutation causes impaired phosphorylation
To confirm that PKA significantly phosphorylates the S510 site in TLE6, we performed in vitro kinase assays with increasing concentrations of purified PKA, and wild-type or mutant TLE6-S510Y in the presence of [γ32P] ATP. PKA successfully phosphorylated wild-type TLE6 (Fig. 2e, lanes 2–5) while mutant TLE6-S510Y had a significantly lower phosphorylation (Fig. 2e, lanes 7–10). The reduced phosphorylation of TLE6-S510Y by PKA was also observed by immunoblotting (Fig. 2e, lower panel). Thus, these results indicate that the TLE6 mutation at S510 impairs its phosphorylation by PKA.
TLE6 mutation affects binding to SCMC proteins
Maternal genes play a critical role in the very early stages of embryonic development because of the lag in transcribing genes derived from the male pronucleus. Although maternal effect genes are very well studied in other model organisms, little is known about them in mammalian oocytes due to practical limitations . Recently, TLE6 was identified as an essential member of the four-protein, murine SCMC along with MATER, Filia, and FLOPED . Mice that are knocked out for Tle6 exhibit complete infertility that is female-specific, which places TLE6 as one of the less than 30 maternal effect genes identified to date in mammals . In sheep, expression of the four homologous proteins was recently linked to oocyte development potential .
The phenotypic overlap between our patients with TLE6 mutation and the corresponding mouse model is noteworthy. Ovulation proceeds normally in Tle6 mutants and the retrieved oocytes appear normal just as we observed in our three patients. The infertility phenotype in Tle6 mutants specifically involves the zygote stage and beyond. While a comparable number of one-cell zygotes was retrieved from Tle6 mutant females following hormonal stimulation and natural mating, there was a protracted delay in undergoing the first cleavage and the few that did, succumbed to disintegration and never formed morula, a phenotype similar to that observed in the few zygotes from our patients with homozygous TLE6 mutation. On the other hand, we note that zygote formation, which appears to have occurred normally in the mouse model, was severely impaired in the human patients. The human homologues of the four SCMC proteins (NLRP5, OOEP, KHDC3L and TLE6) share 39–46 % identity with their mouse counterparts , so the mechanism for this difference may lie in interspecies divergence. Nonetheless, these observations suggest that the phenotype is not caused by an impaired receptor mediated sperm-egg physical interaction as shown for Juno mutants but rather a maternal gene that exerts its effect after the introduction of the sperm .
The mechanism through which TLE6 exerts its maternal effect is not completely understood. Tle6 −/− oocytes lack the normal subcortical distribution of the SCMC despite normal abundance of RNA from their component genes, which clearly shows an essential role for TLE6 in stabilizing this complex . Indeed, the Tle6 mutant phenotype is very similar to that observed in knockouts of the other SCMC components indicating that stabilized SCMC is necessary for the very early stages of embryonic development [22, 26, 27]. Interestingly, this early developmental role of TLE6 was independently discovered through a study that examined the dynamic function of PKA in meiotic maturation of mammalian oocytes . Suppressing PKA resulted in oocytes that were blocked in meiosis II and TLE6 was found to be a major substrate of PKA in this process. TLE6 phosphorylation occurs in a very narrow window within meiosis II and when this was blocked, oocytes failed to mature.
Our data, based on patient cells as well as the mutant Flag-TLE6 plasmid, demonstrate that the S510Y mutation results in markedly reduced phosphorylation of TLE6 and that this phosphorylation is mostly catalyzed by PKA. Mutant TLE6 also exhibits impaired binding ability to at least two SCMC component proteins. Taken together, these data reveal two, potentially related, phenotypic mechanisms: reduced potential for PKA-catalyzed phosphorylation and reduced potential to bind SCMC. Although it is unclear whether the patient phenotype is due to the absence of phosphorylation of S510, or its substitution with tyrosine, we note that both PKA-catalyzed phosphorylation and binding to SCMC have been proposed in the literature to explain the critical role of maternal TLE6 in supporting peri-implantation embryogenesis [18, 20, 23]. Since serine is a canonical phosphorylation site, it is difficult to separate the two roles and, in fact, our data suggest they are related.
To the best of our knowledge, this is the first reported mutation in an SCMC gene in humans. Our data suggest that TLE6 mutations are a rare cause of human female-limited infertility and represent the earliest known human embryonic lethality that is explicable by a single gene mutation. Our ongoing analysis of embryonic lethal phenotypes is likely to reveal additional novel genes and contribute to the understanding of the fundamental biological processes that control early human development.
We thank the families for their enthusiastic participation. We are grateful to Salma Wakil, Dorota Monies, Faisal Alotaibi, the KFSHRC Sequencing and Genotyping Core Facilities, and to the Logistics Office at KFSHRC for their invaluable help. This work was supported in part by KACST 13-BIO1113-20 (FSA) and KACST 13-MED456-20 (CP).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Alkuraya FS. Human knockout research: new horizons and opportunities. Trends Genet. 2015;31:108–15.View ArticlePubMedGoogle Scholar
- Shaheen R, Faqeih E, Shamseldin HE, Noche RR, Sunker A, Alshammari MJ, et al. POC1A truncation mutation causes a ciliopathy in humans characterized by primordial dwarfism. Am J Hum Genet. 2012;91:330–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Shamseldin HE, Tulbah M, Kurdi W, Nemer M, Alsahan N, Al Mardawi E, et al. Identification of embryonic lethal genes in humans by autozygosity mapping and exome sequencing in consanguineous families. Genome Biol. 2015;16:116.PubMed CentralView ArticlePubMedGoogle Scholar
- Khalak HG, Wakil SM, Imtiaz F, Ramzan K, Baz B, Almostafa A, et al. Autozygome maps dispensable DNA and reveals potential selective bias against nullizygosity. Genet Med. 2012;14:515–9.View ArticlePubMedGoogle Scholar
- Ayadi A, Birling M-C, Bottomley J, Bussell J, Fuchs H, Fray M, et al. Mouse large-scale phenotyping initiatives: Overview of the European mouse disease clinic (EUMODIC) and of the Wellcome Trust Sanger Institute Mouse Genetics Project. Mamm Genome. 2012;23:600–10.PubMed CentralView ArticlePubMedGoogle Scholar
- Giaever G, Chu AM, Ni L, Connelly C, Riles L, Véronneau S, et al. Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002;418:387–91.View ArticlePubMedGoogle Scholar
- Alkuraya FS. Discovery of rare homozygous mutations from studies of consanguineous pedigrees. Curr Protoc Hum Genet. 2012;6:6.12.Google Scholar
- Alkuraya FS. Autozygome decoded. Genet Med. 2010;12:765–71.View ArticlePubMedGoogle Scholar
- Alkuraya FS. The application of next-generation sequencing in the autozygosity mapping of human recessive diseases. Hum Genet. 2013;132:1197–211.View ArticlePubMedGoogle Scholar
- Carr IM, Flintoff KJ, Taylor GR, Markham AF, Bonthron DT. Interactive visual analysis of SNP data for rapid autozygosity mapping in consanguineous families. Hum Mutat. 2006;27:1041.View ArticlePubMedGoogle Scholar
- Alazami AM, Patel N, Shamseldin HE, Anazi S, Al-Dosari MS, Alzahrani F, et al. Accelerating novel candidate gene discovery in neurogenetic disorders via whole-exome sequencing of prescreened multiplex consanguineous families. Cell Rep. 2015;10:148–61.View ArticlePubMedGoogle Scholar
- Lindner TH, Hoffmann K. easyLINKAGE: a PERL script for easy and automated two-/multi-point linkage analyses. Bioinformatics. 2005;21:405–7.View ArticlePubMedGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.PubMed CentralView ArticlePubMedGoogle Scholar
- Awad S, Al-Dosari MS, Al-Yacoub N, Colak D, Salih MA, Alkuraya FS, et al. Mutation in PHC1 implicates chromatin remodeling in primary microcephaly pathogenesis. Hum Mol Genet. 2013;22:2200–13.View ArticlePubMedGoogle Scholar
- Duncan FE, Padilla-Banks E, Bernhardt ML, Ord TS, Jefferson WN, Moss SB, et al. Transducin-like enhancer of split-6 (TLE6) is a substrate of protein kinase A activity during mouse oocyte maturation 1. Biol Reprod. 2014;90:63.PubMed CentralView ArticlePubMedGoogle Scholar
- Kinoshita E, Kinoshita-Kikuta E, Takiyama K, Koike T. Phosphate-binding tag, a new tool to visualize phosphorylated proteins. Mol Cell Proteomics. 2006;5:749–57.View ArticlePubMedGoogle Scholar
- Zhu K, Yan L, Zhang X, Lu X, Wang T, Yan J, et al. Identification of a human subcortical maternal complex. Mol Hum Reprod. 2015;21:320–9.View ArticlePubMedGoogle Scholar
- Li L, Zheng P, Dean J. Maternal control of early mouse development. Development. 2010;137:859–70.PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Baibakov B, Dean J. A subcortical maternal complex essential for preimplantation mouse embryogenesis. Dev Cell. 2008;15:416–25.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu XJ, Yi Z, Gao Z, Qin D, Zhai Y, Chen X, et al. The subcortical maternal complex controls symmetric division of mouse zygotes by regulating F-actin dynamics. Nat Commun. 2014;5:4887.View ArticlePubMedGoogle Scholar
- Bebbere D, Ariu F, Bogliolo L, Masala L, Murrone O, Fattorini M, et al. Expression of maternally derived KHDC3, NLRP5, OOEP and TLE6 is associated with oocyte developmental competence in the ovine species. BMC Dev Biol. 2014;14:40.PubMed CentralView ArticlePubMedGoogle Scholar
- Bianchi E, Doe B, Goulding D, Wright GJ. Juno is the egg Izumo receptor and is essential for mammalian fertilization. Nature. 2014;508:483–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Tong Z-B, Bondy CA, Zhou J, Nelson LM. A human homologue of mouse Mater, a maternal effect gene essential for early embryonic development. Hum Reprod. 2002;17:903–11.View ArticlePubMedGoogle Scholar
- Zheng P, Dean J. Role of Filia, a maternal effect gene, in maintaining euploidy during cleavage-stage mouse embryogenesis. Proc Natl Acad Sci. 2009;106:7473–8.PubMed CentralView ArticlePubMedGoogle Scholar