Case Report ARTICLE
Identification and Characterization of Known Biallelic Mutations in the IFT27 (BBS19) Gene in a Novel Family With Bardet-Biedl Syndrome
- 1Laboratoire de Génétique Médicale, Institut de Génétique Médicale d’Alsace, INSERM U1112, Fédération de Médecine Translationnelle de Strasbourg, Université de Strasbourg, Strasbourg, France
- 2Service de Génétique Médicale, Hôpitaux Universitaires de Strasbourg, Strasbourg, France
- 3Laboratoire de Diagnostic Génétique, Hôpitaux Universitaires de Strasbourg, Strasbourg, France
- 4Service d’Exploration de la Vision et Neuro-Ophtalmologie, Hôpital Roger Salengro, CHRU Lille, Lille, France
- 5Clinical Genetics Department, Guy’s Hospital, London, United Kingdom
- 6Service de Génétique Clinique, Hôpital Jeanne de Flandre, CHRU Lille, Lille, France
Bardet-Biedl syndrome (BBS; MIM 209900) is a rare ciliopathy characterized by retinitis pigmentosa, postaxial polydactyly, obesity, hypogonadism, cognitive impairment and kidney dysfunction. Mutations in 22 BBS genes have been identified to cause the disease. We report a family with typical BBS features (retinitis pigmentosa, postaxial polydactyly, obesity, cognitive impairment, and atrioventricular septal defect) mutated in IFT27/BBS19. IFT27 is part of the Intraflagellar transport (IFT), a bidirectional mechanism allowing the protein motility within the cilia. Using whole exome sequencing, two compound heterozygous mutations were found in the proband (NM_006860.4:c.[104A > G];[349+1G > T], p.[Tyr35Cys];[?]) consistent with the expected autosomal recessive inheritance mode. These two mutations have already been reported but independently in other families and lacking either familial segregation or functional validation. This is the third report of IFT27 mutations in BBS patients confirming IFT27 as a BBS gene (BBS19). Mutations in IFT genes (IFT27, IFT172 and IFT74) confirm the IFT-pathway as a pathomechanism for BBS.
Bardet-Biedl syndrome (BBS; OMIM 209900) is a recessive and genetically heterogeneous ciliopathy defined by the association of retinitis pigmentosa (RP), postaxial polydactyly, obesity, hypogonadism, kidney dysfunction and cognitive impairment. Mutations in 22 genes have been identified of which WDPCP/BBS15, LZTFL1/BBS17, BBIP1/BBS18 IFT27/BBS19, IFT172/BBS20, C8Orf37/BBS21, and IFT74/BBS22 only account for one or a few families (Khan et al., 2016). Among the last 4 genes identified, 3 (IFT27/BBS19, IFT172/BBS20, and the candidate locus IFT74/BBS22) are members of the intraflagellar transport machinery (IFT) delineating a possible novel mechanisms for the BBS (Aldahmesh et al., 2014; Bujakowska et al., 2015; Lindstrand et al., 2016). The IFT machinery is composed of 2 large protein complexes named IFT-A and IFT-B required, respectively, for retrograde and anterograde transport in the cilia. This machinery is part of an essential process for the assembly and the maintenance of the cilia (Taschner and Lorentzen, 2016). Most of the reported mutations in these genes, and in particular for IFT172, have been identified in patients with skeletal forms of ciliopathies including the Jeune asphyxiating thoracic dystrophy (OMIM 208500), Mainzer–Saldino syndrome (OMIM 266920), the Sensenbrenner syndrome (OMIM 218330).
In this study, we report and validate 2 compound heterozygous mutations in IFT27 confirming the 19th BBS locus in a family with a clear BBS phenotype.
Materials and Methods
Bardet-Biedl syndrome studies have been approved by our Institutional Review Board “Comité de Protection des Personnes” (EST IV, N°DC-20142222) and written informed consent for research and publication was obtained from each participant. We declare that the present research complies with the declaration of Helsinki.
Whole exome sequencing (WES) was performed by the IGBMC Microarray and Sequencing platform. Genomic DNA (2 μg) was sheared to obtain a mean fragment size of 150 bp using the Covaris E210 (KBioscience, Herts, United Kingdom) followed by library preparation using the Agilent SureSelect XT Human all exon V6 (Agilent Technologies, Santa Clara, CA, United States PN G7530-90000, protocol B4). Sequencing was performed on an Illumina HiSeq 4000 (Illumina, San Diego, CA, United States) to generate 100-bp paired-end reads following the manufacturer’s protocols.
Image analysis and base calling were performed using CASAVA v1.8.2 (Illumina) and the 109,658,782 reads were mapped to the reference human genome (GRCh37/hg19) using BWA v0.7.5a (Li and Durbin, 2009) leading to 95.89% of the bases covered at least by 20x. GATK UG v3.2-2 was used to call SNV and indel variations (DePristo et al., 2011). Annotation and ranking of SNV and indel were performed by VaRank (Geoffroy et al., 2015) in combination with the Alamut Batch software (Interactive Biosoftware, Rouen, France). Very stringent filtering criteria were applied to filter out non-pathogenic variants (Supplementary Table S1): (i) variants represented with an allele frequency of more than 1% in public variation databases including the 1000Genomes (Genomes Project et al., 2015), the gnomAD database (Lek et al., 2016), the DGV database (MacDonald et al., 2014) or our internal exome database, (ii) variants in 5′ and 3′ UTR, downstream, upstream or intronic locations and synonymous without pathogenic prediction of local splice effect, (iii) variants not in the ciliary genes (Nevers et al., 2017). Structural variants were predicted using CANOES (Backenroth et al., 2014) and annotated by AnnotSV (Geoffroy et al., 2018). Our analysis was focused on compound heterozygous and homozygous variants consistent with a recessive mode of transmission. The IFT27 nomenclature is based on the RefSeq accession number NM_006860.4 (O’Leary et al., 2016).
RNA Extraction, cDNA Synthesis, and Sanger Sequencing
RNA was extracted from the patient’s blood using the PAXgene Blood RNA Kit (PreAnalytiX GmbH, Hombrechtikon, Switzerland). Reverse transcription was performed on 200 ng of RNA using the BioRad iScript cDNA® Synthesis Kit (BioRad, Hercules, CA, United States). PCR was performed on the cDNA using the Mastercycler ep. Gradient S thermocycler (Eppendorf, Germany). Specific fragment analysis was done by cutting the agarose gel around the fragment of interest (e.g., F3) and eluted in Tris–EDTA buffer 10.1 overnight. The elution product (5 μL) was reamplified by PCR using the same conditions to obtain a unique fragment. Bidirectional sequencing of the purified PCR products was performed by the GATC Sequencing Facilities (Konstanz, Germany). Primers used are summarized in Supplementary Table S2.
The patient is the only child of an unrelated couple without personal or familial medical history. The patient was born at 39 weeks of amenorrhea with the following parameters: weight at 3210 g, height at 51 cm and head circumference at 35 cm. Mesoaxial polydactyly of the right hand with a Y-shaped metacarpian and syndactyly between the 5th and the 6th fingers and postaxial polydactyly of the right foot was noticed at birth. Partial atrioventricular septal defect was also diagnosed and operated at 5 weeks old. Mitral insufficiency persisted after operation and was operated at 2 years old. Renal ultrasound was normal. At birth, cerebral ultrasound showed isolated thin corpus callosum, not confirmed on cerebral MRI. Secondarily, the patient presented with delayed psychomotor development: he walked at 25-month-old and had delayed language (5 words at 2 years; 10 words at 3 years and sentences at 5 years). Audition was normal. The patient received specialized education. Progressively, he developed obesity: at 2 years old, 16.8 kg (+ 3SD) for 93 cm (+ 2SD) and normal head circumference (49 cm); at 3 years old, 22.9 kg (> +3SD) for 101.5 cm (+ 2SD). At 7 years old, his BMI was 25. Initially, ophthalmologic examination revealed myopia at 2 years old and alternate divergent strabismus at 3 years old. Myopia was scalable was noticed at 4 years old. A cone-rod dystrophy was diagnosed at 7 years old.
Several diagnostic hypotheses led to the non-conclusive exploration of GLI3 (Pallister-Hall), OFD1 (Orofaciodigital syndrome I) and the known BBS genes at that time (BBS1-BBS18) using Sanger and targeted exome sequencing (Muller et al., 2010; Redin et al., 2012). This also included array-CGH screening. Finally, WES was performed on the index case and led to the identification of 2 already reported variations in the IFT27/BBS19 gene co-segregating in the family (Figure 1A). The first variation, c.104A > G (p.Tyr35Cys), is affecting a highly conserved position in the protein localized in the switch I region of this Rab-like GTPase (Figures 2A,B). The variation is absent from the variation databases (i.e., 1000 genomes and gnomAD) and predicted as affecting the protein function (Supplementary Table S3). Given its relatively close position toward the end of exon 2 (Figure 2A), a possible splice effect was assessed on the patient’s RNA extracted from blood and failed to reveal any alteration (Supplementary Figure S1).
Figure 1. Identification of two mutations in IFT27. (A) Pedigree of the reported family with one affected individual and segregation analysis of the two IFT27 mutations. Example of Sanger sequencing profiles for the heterozygous individuals. (B) PCR amplification was performed on RNA extracted from blood of individual II.1 and a healthy unrelated control amplified between exon 2 and exon 7. (C) IFT27 cDNA scheme representing the obtained fragments with size and expected composition. PCR primers are positioned. (D) Sanger sequencing of normal F1 in a healthy unrelated control (left side showing each exon boundaries from exon 3 to 7) and cut and eluted F3 band in individual II.1 demonstrating the absence of exons 4 to 6 (right side).
Figure 2. IFT27 gene and protein. (A) Schematic of the IFT27 locus and (B) IFT27 protein and the position and/or nature of known mutations. Nomenclature is given according to the following RefSeq identifiers: NM_006860.4 and NP_006851.1. (C) Multiple sequence alignment of IFT27 proteins from different species and of the human RAB8a and RAB10 proteins. The secondary structures derived from the IFT27 structure are shown in black above the sequences (Bhogaraju et al., 2011). The mutation Tyr35 is highlighted using a black “∗” and the Switch I region is shown using a black line below the sequences. Residues are presented and colored according to the Jalview program (Waterhouse et al., 2009). Hs, Homo sapiens; Mm, Mus musculus; Xl, Xenopus laevis; Dr, Danio rerio; Cr, Chlamydomonas reinhardtii.
Given its localization, the second variation, c.349 + 1G > T, is predicted to abolish the donor splice site of exon 5 (Figure 2A and Supplementary Table S3). It is already described in gnomAD with 21 heterozygous alleles in the European (Non-Finnish) cohort leading to an allele frequency of 0.0085% compatible with the BBS rare condition. Analysis of the patient’s RNA extracted from blood revealed a mix of alternatively spliced isoforms not found in controls (Figures 1B–D). This includes removal of exons 5 + 6 or 4 + 5 + 6 that are predicted to cause an in frame deletion of a significant part of the protein (76 or 96aa out of 185aa).
Bardet-Biedl syndrome is an emblematic ciliopathy characterized by a clinical and genetic heterogeneity. Twenty-two BBS genes have been identified so far and mutations are found in about 80% of the known BBS patients (Redin et al., 2012; Khan et al., 2016). The last genes identified account only for a small fraction of the mutation load, emphasizing the importance to screen more families and increase our knowledge for these genes (e.g., type of mutations and clinical manifestations).
In our study, we identified two compound heterozygous mutations in IFT27 (NM_006860.4: c.[104A > G];[349 + 1G > T]) by WES in a child presenting with the classical BBS signs. IFT27 is part of the RAb like GTPAse that is predicted to have low GTPAse activity. The gene is composed of 7 exons coding for a 186 aa protein.
Our first identified mutation, the c.104A > G is located in exon 2 leading to the drastic change of a tyrosine to a cysteine (p.Tyr35Cys). Localized in the connecting region between α-helix 1 and β-sheet 2 in the so called “switch I” region (Figure 2C), a highly flexible region upon nucleotide binding pocket not crystallized (Bhogaraju et al., 2011) with an important role in GDP/GTP exchange and with interacting partners (Pylypenko et al., 2018), the missense is highly susceptible to affect the function of IFT27. The second mutation, c.349 + 1G > T, located in the donor splice site of exon 5 was predicted to affect the splicing of IFT27 that we have demonstrated with a skipping of several exons of the gene. This mutation has been already reported in a fetus with a severe phenotype but the effect could not be assessed at that time (Quelin et al., 2018).
Initially the first identified mutations in IFT27 were reported in a consanguineous BBS family presenting with RP, obesity, polydactyly of all extremities, mild intellectual disability, renal failure and hypogenitalism. The patients carried a homozygous missense (c.296G > A, p.Cys99Tyr) predicted to affect the stability of the protein (Aldahmesh et al., 2014).
Very recently, another BBS case with RP (cone rod dystrophy), obesity, polydactyly, maturation and learning delay and chronic renal failure has been reported but they could not to assess the segregation and the effect on the protein (c.[104A>G(;)350-2A>G], p.[Tyr35(;)?]) (Sanchez-Navarro et al., 2018).
Lastly, 2 loss of function mutations (c.[115_122del];[349+1G > T], p.[Thr39Glyfs∗11];[?]) have been reported in a fetus presenting with a severe ciliopathy with short ribs polydactyly type II (SRPII or Majewski syndrome) and/or Pallister-Hall syndrome (Quelin et al., 2018). Interestingly, Pallister-Hall syndrome was also evocated for our patient given the association of cardiopathy and polydactyly with Y-Shaped metacarpian. The authors made the hypothesis that the loss of function mutations could explain the severe phenotype.
Comparing the reported patients including the fetus (Supplementary Table S4), all had postaxial polydactyly. All living patients had RP, learning difficulties and obesity. Moreover, one patient of the initial BBS family has deafness and hypogenitalism. The same family and the fetus have renal anomalies (renal hypoplasia to renal agenesis) contrasting with our patient who has no renal involvement. Finally, our patient and one patient of the initial family presented with a cardiopathy, and unexpectedly not found in the more severe fetal form. Interestingly some patients have unusual clinical presentation including mesoaxial polydactyly with Y-shaped metacarpian associated with postaxial polydactyly (our patient) or preaxial (fetus). The Y-shaped metacarpian has been already correlated with mutation in BBS17/LZTFL1 (Schaefer et al., 2014) that is now extended to IFT27. One can notice also that 2 patients had cone-rod dystrophy (our patient and family 2), a rare condition already described in BBS (Scheidecker et al., 2015). These features show the clinical variability associated to IFT27 mutations from severe lethal forms to classical BBS that might be revised when new patients will be described in the future.
Mutations in IFT genes were recently reported in several BBS patients (Bujakowska et al., 2015; Lindstrand et al., 2016; Schaefer et al., 2016). IFT27, IFT74, and IFT172 encode components of the IFT-B complex required for anterograde transport of ciliary proteins (Taschner and Lorentzen, 2016). Mutations in IFT genes are implicated in different ciliopathies, principally in skeletal ciliopathies (Jeune asphyxiating thoracic dystrophy, Sensenbrenner syndrome and Mainzer-Saldino syndrome) but also in isolated RP and Senior-Løken syndrome (associating RP to nephronopthisis). This report confirms the IFT-pathway as a new pathomechanism for BBS.
We identified the third BBS family mutated in IFT27. We confirmed IFT27 as the 19th BBS gene delineating the typical and complete BBS phenotype for those rarely mutated patients and the implication of IFT-pathway in the BBS pathogenesis. This report illustrated the usefulness of WES sequencing to identify mutations in highly heterogeneous genetic disorders. Finally, our observation is a good illustration of the clinical and molecular continuum between the different ciliopathies.
All variants have been submitted to ClinVar using the following range of accessions numbers SCV000839871 and SCV000839872 (https://www.ncbi.nlm.nih.gov/clinvar/).
HD and JM conceived the study. CM-D, JG, and MH-E collected the clinical information. CD, CS, ES, and LM analyzed the data. VG conducted the bioinformatics analysis. ES and JM drafted the manuscript. HD and JM supervised the entire study. All the authors approved the final version of manuscript.
CD was supported by a FRM grant (ECO20170637509). Whole exome sequencing was performed by the IGBMC Microarray and Sequencing platform, a member of the ‘France Génomique’ consortium (ANR-10-INBS-0009) and funded by “La Fondation Maladie Rare” through the programme «High throughput sequencing and rare diseases».
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank the patient’s family for their participation. We would like also to acknowledge the members of the diagnostic laboratories (molecular and cytogenetic units, in particular Manuela Antin, Anne-Sophie Leuvrey and Elsa Nourisson) at the University Hospital in Strasbourg for technical assistance. The computing resources for this work were provided by the BICS and BISTRO bioinformatics platforms in Strasbourg.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00021/full#supplementary-material
Aldahmesh, M. A., Li, Y., Alhashem, A., Anazi, S., Alkuraya, H., Hashem, M., et al. (2014). IFT27, encoding a small GTPase component of IFT particles, is mutated in a consanguineous family with Bardet-Biedl syndrome. Hum. Mol. Genet. 23, 3307–3315. doi: 10.1093/hmg/ddu044
Backenroth, D., Homsy, J., Murillo, L. R., Glessner, J., Lin, E., Brueckner, M., et al. (2014). CANOES: detecting rare copy number variants from whole exome sequencing data. Nucleic Acids Res. 42:e97. doi: 10.1093/nar/gku345
Bujakowska, K. M., Zhang, Q., Siemiatkowska, A. M., Liu, Q., Place, E., Falk, M. J., et al. (2015). Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome. Hum. Mol. Genet. 24, 230–242. doi: 10.1093/hmg/ddu441
DePristo, M. A., Banks, E., Poplin, R., Garimella, K. V., Maguire, J. R., Hartl, C., et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498. doi: 10.1038/ng.806
Geoffroy, V., Herenger, Y., Kress, A., Stoetzel, C., Piton, A., Dollfus, H., et al. (2018). AnnotSV: an integrated tool for structural variations annotation. Bioinformatics 34, 3572–3574. doi: 10.1093/bioinformatics/bty304
Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., et al. (2016). Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285–291. doi: 10.1038/nature19057
Lindstrand, A., Frangakis, S., Carvalho, C. M., Richardson, E. B., McFadden, K. A., Willer, J. R., et al. (2016). Copy-number variation contributes to the mutational load of bardet-biedl syndrome. Am. J. Hum. Genet. 99, 318–336. doi: 10.1016/j.ajhg.2015.04.023
MacDonald, J. R., Ziman, R., Yuen, R. K., Feuk, L., and Scherer, S. W. (2014). The database of genomic variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 42, D986–D992. doi: 10.1093/nar/gkt958
Muller, J., Stoetzel, C., Vincent, M. C., Leitch, C. C., Laurier, V., Danse, J. M., et al. (2010). Identification of 28 novel mutations in the Bardet-Biedl syndrome genes: the burden of private mutations in an extensively heterogeneous disease. Hum. Genet. 127, 583–593. doi: 10.1007/s00439-010-0804-9
Nevers, Y., Prasad, M. K., Poidevin, L., Chennen, K., Allot, A., Kress, A., et al. (2017). Insights into ciliary genes and evolution from multi-level phylogenetic profiling. Mol. Biol. Evol. 34, 2016–2034. doi: 10.1093/molbev/msx146
O’Leary, N. A., Wright, M. W., Brister, J. R., Ciufo, S., Haddad, D., McVeigh, R., et al. (2016). Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 44, D733–D745. doi: 10.1093/nar/gkv1189
Pylypenko, O., Hammich, H., Yu, I. M., and Houdusse, A. (2018). Rab GTPases and their interacting protein partners: structural insights into Rab functional diversity. Small GTPases 9, 22–48. doi: 10.1080/21541248.2017.1336191
Quelin, C., Loget, P., Boutaud, L., Elkhartoufi, N., Milon, J., Odent, S., et al. (2018). Loss of function IFT27 variants associated with an unclassified lethal fetal ciliopathy with renal agenesis. Am. J. Med. Genet. A 176, 1610–1613. doi: 10.1002/ajmg.a.38685
Redin, C., Le Gras, S., Mhamdi, O., Geoffroy, V., Stoetzel, C., Vincent, M. C., et al. (2012). Targeted high-throughput sequencing for diagnosis of genetically heterogeneous diseases: efficient mutation detection in Bardet-Biedl and Alstrom syndromes. J. Med. Genet. 49, 502–512. doi: 10.1136/jmedgenet-2012-100875
Sanchez-Navarro, R. J. I, da Silva, L., Blanco-Kelly, F., Zurita, O., Sanchez-Bolivar, N., Villaverde, C., et al. (2018). Combining targeted panel-based resequencing and copy-number variation analysis for the diagnosis of inherited syndromic retinopathies and associated ciliopathies. Sci. Rep. 8:5285. doi: 10.1038/s41598-018-23520-1
Schaefer, E., Lauer, J., Durand, M., Pelletier, V., Obringer, C., Claussmann, A., et al. (2014). Mesoaxial polydactyly is a major feature in Bardet-Biedl syndrome patients with LZTFL1 (BBS17) mutations. Clin. Genet. 85, 476–481. doi: 10.1111/cge.12198
Schaefer, E., Stoetzel, C., Scheidecker, S., Geoffroy, V., Prasad, M. K., Redin, C., et al. (2016). Identification of a novel mutation confirms the implication of IFT172 (BBS20) in Bardet-Biedl syndrome. J. Hum. Genet. 61, 447–450. doi: 10.1038/jhg.2015.162
Scheidecker, S., Hull, S., Perdomo, Y., Studer, F., Pelletier, V., Muller, J., et al. (2015). Predominantly cone-system dysfunction as rare form of retinal degeneration in patients with molecularly confirmed Bardet-Biedl syndrome. Am. J. Ophthalmol. 160, 364.e1–372.e1. doi: 10.1016/j.ajo.2015.05.007
Keywords: Bardet-Biedl syndrome, IFT27 gene, BBS19, ciliopathy, whole exome sequencing
Citation: Schaefer E, Delvallée C, Mary L, Stoetzel C, Geoffroy V, Marks-Delesalle C, Holder-Espinasse M, Ghoumid J, Dollfus H and Muller J (2019) Identification and Characterization of Known Biallelic Mutations in the IFT27 (BBS19) Gene in a Novel Family With Bardet-Biedl Syndrome. Front. Genet. 10:21. doi: 10.3389/fgene.2019.00021
Received: 29 October 2018; Accepted: 15 January 2019;
Published: 30 January 2019.
Edited by:Musharraf Jelani, King Abdulaziz University, Saudi Arabia
Reviewed by:Muhammad Tariq, University of Tabuk, Saudi Arabia
Sulman Basit, Quaid-i-Azam University, Pakistan
Copyright © 2019 Schaefer, Delvallée, Mary, Stoetzel, Geoffroy, Marks-Delesalle, Holder-Espinasse, Ghoumid, Dollfus and Muller. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jean Muller, email@example.com