Abstract
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.
Introduction
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/BBS18IFT27/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
Subjects
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
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.
Bioinformatics Analysis
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.
Results
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 (Figure 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
FIGURE 2
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 (Figure 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).
Discussion
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.
Conclusion
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.
Statements
Data availability statement
All variants have been submitted to ClinVar using the following range of accessions numbers SCV000839871 and SCV000839872 (https://www.ncbi.nlm.nih.gov/clinvar/).
Author contributions
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.
Funding
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».
Acknowledgments
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.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00021/full#supplementary-material
References
1
AldahmeshM. A.LiY.AlhashemA.AnaziS.AlkurayaH.HashemM.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.233307–3315. 10.1093/hmg/ddu044
2
BackenrothD.HomsyJ.MurilloL. R.GlessnerJ.LinE.BruecknerM.et al (2014). CANOES: detecting rare copy number variants from whole exome sequencing data.Nucleic Acids Res.42:e97. 10.1093/nar/gku345
3
BhogarajuS.TaschnerM.MorawetzM.BasquinC.LorentzenE. (2011). Crystal structure of the intraflagellar transport complex 25/27.EMBO J.301907–1918. 10.1038/emboj.2011.110
4
BujakowskaK. M.ZhangQ.SiemiatkowskaA. M.LiuQ.PlaceE.FalkM. J.et al (2015). Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome.Hum. Mol. Genet.24230–242. 10.1093/hmg/ddu441
5
DePristoM. A.BanksE.PoplinR.GarimellaK. V.MaguireJ. R.HartlC.et al (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data.Nat. Genet.43491–498. 10.1038/ng.806
6
Genomes ProjectC.AutonA.BrooksL. D.DurbinR. M.GarrisonE. P.KangH. M.et al (2015). A global reference for human genetic variation.Nature52668–74. 10.1038/nature15393
7
GeoffroyV.HerengerY.KressA.StoetzelC.PitonA.DollfusH.et al (2018). AnnotSV: an integrated tool for structural variations annotation.Bioinformatics343572–3574. 10.1093/bioinformatics/bty304
8
GeoffroyV.PizotC.RedinC.PitonA.VasliN.StoetzelC.et al (2015). VaRank: a simple and powerful tool for ranking genetic variants.PeerJ3e796. 10.7717/peerj.796
9
KhanS. A.MuhammadN.KhanM. A.KamalA.RehmanZ. U.KhanS. (2016). Genetics of human Bardet-Biedl syndrome, an updates.Clin. Genet.903–15. 10.1111/cge.12737
10
LekM.KarczewskiK. J.MinikelE. V.SamochaK. E.BanksE.FennellT.et al (2016). Analysis of protein-coding genetic variation in 60,706 humans.Nature536285–291. 10.1038/nature19057
11
LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform.Bioinformatics251754–1760. 10.1093/bioinformatics/btp324
12
LindstrandA.FrangakisS.CarvalhoC. M.RichardsonE. B.McFaddenK. A.WillerJ. R.et al (2016). Copy-number variation contributes to the mutational load of bardet-biedl syndrome.Am. J. Hum. Genet.99318–336. 10.1016/j.ajhg.2015.04.023
13
MacDonaldJ. R.ZimanR.YuenR. K.FeukL.SchererS. W. (2014). The database of genomic variants: a curated collection of structural variation in the human genome.Nucleic Acids Res.42D986–D992. 10.1093/nar/gkt958
14
MullerJ.StoetzelC.VincentM. C.LeitchC. C.LaurierV.DanseJ. 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.127583–593. 10.1007/s00439-010-0804-9
15
NeversY.PrasadM. K.PoidevinL.ChennenK.AllotA.KressA.et al (2017). Insights into ciliary genes and evolution from multi-level phylogenetic profiling.Mol. Biol. Evol.342016–2034. 10.1093/molbev/msx146
16
O’LearyN. A.WrightM. W.BristerJ. R.CiufoS.HaddadD.McVeighR.et al (2016). Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation.Nucleic Acids Res.44D733–D745. 10.1093/nar/gkv1189
17
PylypenkoO.HammichH.YuI. M.HoudusseA. (2018). Rab GTPases and their interacting protein partners: structural insights into Rab functional diversity.Small GTPases922–48. 10.1080/21541248.2017.1336191
18
QuelinC.LogetP.BoutaudL.ElkhartoufiN.MilonJ.OdentS.et al (2018). Loss of function IFT27 variants associated with an unclassified lethal fetal ciliopathy with renal agenesis.Am. J. Med. Genet. A1761610–1613. 10.1002/ajmg.a.38685
19
RedinC.Le GrasS.MhamdiO.GeoffroyV.StoetzelC.VincentM. 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.49502–512. 10.1136/jmedgenet-2012-100875
20
Sanchez-NavarroR. J. Ida SilvaL.Blanco-KellyF.ZuritaO.Sanchez-BolivarN.VillaverdeC.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. 10.1038/s41598-018-23520-1
21
SchaeferE.LauerJ.DurandM.PelletierV.ObringerC.ClaussmannA.et al (2014). Mesoaxial polydactyly is a major feature in Bardet-Biedl syndrome patients with LZTFL1 (BBS17) mutations.Clin. Genet.85476–481. 10.1111/cge.12198
22
SchaeferE.StoetzelC.ScheideckerS.GeoffroyV.PrasadM. K.RedinC.et al (2016). Identification of a novel mutation confirms the implication of IFT172 (BBS20) in Bardet-Biedl syndrome.J. Hum. Genet.61447–450. 10.1038/jhg.2015.162
23
ScheideckerS.HullS.PerdomoY.StuderF.PelletierV.MullerJ.et al (2015). Predominantly cone-system dysfunction as rare form of retinal degeneration in patients with molecularly confirmed Bardet-Biedl syndrome.Am. J. Ophthalmol.160364.e1–372.e1. 10.1016/j.ajo.2015.05.007
24
TaschnerM.LorentzenE. (2016). The Intraflagellar Transport Machinery.Cold Spring Harb. Perspect. Biol.8:a028092. 10.1101/cshperspect.a028092
25
WaterhouseA. M.ProcterJ. B.MartinD. M.ClampM.BartonG. J. (2009). Jalview version 2–a multiple sequence alignment editor and analysis workbench.Bioinformatics251189–1191. 10.1093/bioinformatics/btp033
Summary
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
Volume
10 - 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
Updates
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, jeanmuller@unistra.fr
This article was submitted to Genetic Disorders, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.