Molecular Vision 2017; 23:482-494
Received 31 December 2016 | Accepted 19 July 2017 | Published 21 July 2017
Asmat Ullah,1 Muhammad Umair,1,4,5 Maryam Yousaf,2 Sher Alam Khan,6 Muhammad Nazim-ud-din,3 Khadim Shah,1 Farooq Ahmad,1 Zahid Azeem,2 Ghazanfar Ali,3 Bader Alhaddad,4,5 Afzal Rafique,1 Abid Jan,1,6 Tobias B. Haack,4,5 Tim M. Strom,4,5 Thomas Meitinger,4,5 Tahseen Ghous,2 Wasim Ahmad1
The first two authors contributed equally to this work.
1Department of Biochemistry, Faculty of Biological Sciences, Quaid-i-Azam University, Islamabad, Pakistan; 2Department of Chemistry, University of Azad Jammu and Kashmir, Muzaffarabad, Pakistan; 3Department of Biotechnology, University of Azad Jammu and Kashmir, Muzaffarabad, Pakistan; 4Institute of Human Genetics, Technische Universitat Munchen, Munchen, Germany; 5Institute of Human Genetics, Helmholtz Zentrum Munchen, Neuherberg, Germany; 6Kohat University of Science and Technology (KUST), Khyber Pakhtunkhwa Province, Pakistan
Correspondence to: Wasim Ahmad, Department of Biochemistry, Faculty of Biological Sciences, Quaid-i-Azam University, Islamabad, Pakistan; Phone: +92-51-90643003; FAX: +92-51-90643003; email: email@example.com
Purpose: To investigate the molecular basis of Bardet-Biedl syndrome (BBS) in five consanguineous families of Pakistani origin.
Methods: Linkage in two families (A and B) was established to BBS7 on chromosome 4q27, in family C to BBS8 on chromosome 14q32.1, and in family D to BBS10 on chromosome 12q21.2. Family E was investigated directly with exome sequence analysis.
Results: Sanger sequencing revealed two novel mutations and three previously reported mutations in the BBS genes. These mutations include two deletions (c.580_582delGCA, c.1592_1597delTTCCAG) in the BBS7 gene, a missense mutation (p.Gln449His) in the BBS8 gene, a frameshift mutation (c.271_272insT) in the BBS10 gene, and a nonsense mutation (p.Ser40*) in the MKKS (BBS6) gene.
Conclusions: Two novel mutations and three previously reported variants, identified in the present study, further extend the body of evidence implicating BBS6, BBS7, BBS8, and BBS10 in causing BBS.
Bardet-Biedl syndrome (BBS; OMIM 209900) is a rare genetically heterogeneous developmental disorder with primary features of retinitis pigmentosa, postaxial polydactyly, obesity, renal dysfunctions, hypogonadism, and intellectual disability. Additional features reported in cases of BBS include strabismus, nystagmus, brachydactyly, syndactyly, truncal obesity, hydronephrosis, pyelonephritis, cryptorchidism, small penis/hypospadias [1,2], delay in reaching puberty, hypoplastic fallopian tubes, poor articulation, anxiety, depression, obsessive-compulsive disorder, autism spectrum disorder or psychosis, hypertension, anosmia, dental anomalies (micrognathia, malocclusion, and microdontia), ataxia, and Hirschsprung disease [1-6]. Variation in the phenotypes within and among different families has been reported [4,5,7]. Facial similarities reported among patients with BBS include deep-set eyes, hypertelorism with downward slanting palpebral fissures, a flat nasal bridge with anteverted nares and prominent nasolabial folds, a long philtrum, a thin upper lip, and apathetic facial appearance [6,8].
BBS belongs to the group of ciliopathies that share partial-overlapping phenotypes and common genes . Genetic studies revealed different variations within the genotype and phenotype of the disease . To date, 21 BBS genes (BBS1–21) have been identified, which includes BBS1 located on chromosome 11q13 , BBS2 on chromosome 16q21 , BBS3 on chromosome 3p13-p12 , BBS4 on chromosome 15q22.3-q23 [11,14], BBS5 on chromosome 2q31 , BBS6 on chromosome 20p12 , BBS7 on chromosome 4q27 , BBS8 on chromosome 14q32.11 , BBS9 on chromosome 7p14 , BBS10 on chromosome 12q21.2 , BBS11 (TRIM32) on chromosome 9q33.1 , BBS12 on chromosome 4q27 , BBS13 (MKS1) on chromosome 17q23 , BBS14 (CEP290) on chromosome 12q21.3 , BBS15 (C2ORF86) on chromosome 2p15 , BBS16 on chromosome 1q43 , BBS17 (LZTFL1) on chromosome 3p21.31 , BBS18 (BBIP1) on chromosome 10q25.2 , BBS19 (IFT27) on chromosome 22q12 , BBS20 (IFT172) on chromosome 2p23.3 , and BBS21 (C8ORF37) on chromosome 8q22.1 . We report two novel mutations and three previously reported variants in the BBS genes in five consanguineous families of Pakistani origin segregating BBS in an autosomal recessive pattern.
Permission to conduct the present research work was obtained from the Institutional Review Board (IRB) of Quaid-i-Azam University, Islamabad, Pakistan and Technical University Munich, Germany. This study adhered to the ARVO statement on human subjects, all those who participated in the study signed informed written consent forms and the research followed the tenets of the Declaration of Helsinki.
In total, five consanguineous families (A–E) segregating Bardet-Biedl syndrome in an autosomal recessive manner were investigated in the present study. Families A and C live in a remote area of the Khyber Pakhtunkhwa (KPK) province of Pakistan. Two other families (B and D) originate from a remote village in Azad Jammu and Kashmir, bordering India and Pakistan. Family E was recruited from a remote area of Nawab Shah City, Sindh province of Pakistan. Information provided by the family elders was used for the construction of the pedigrees (Figure 1, Figure 2A). All five pedigrees convincingly support the autosomal recessive inheritance pattern of the disease.
Affected members in the families underwent clinical examinations at local government hospitals. Venous blood samples from 38 members were collected in vacutainers (BD Biosciences, Franklin Lakes, NJ) containing EDTA.
Genomic DNA was extracted from the collected blood samples using the Nucleospin® Blood kit (Macherey-Nagel, Germany). A NanoDrop-1000 Spectrophotometer (Thermal Scientific, Wilmington, NJ) was used for DNA quantification, measuring optical density at 260 nm and diluted to 40–50 ng/μl for amplification with PCR. PCR was performed in 25 μl reaction volume containing 40 ng genomic DNA, 20 pmol of each primer, 200 mM of each deoxynucleoside triphosphate (dNTP), 2.5 μl reaction buffer (KCl 50 mM, Tris–HCl pH 8.3, and MgCl2 1.5 mM), and 1 U Taq DNA polymerase (MBI Fermentas, Life Sciences, York, UK). PCR was performed using the GeneAmp® PCR System 9700 obtained from Applied Biosystems (Applera Corp, Foster City, CA). PCR conditions (1st cycle: 95 °C for 5 min; 39 cycles: 95 °C for 30 s, 53-58 °C for 30 s, 72 °C for 30 s; 72 °C for 5 min). Amplified PCR products were resolved on 8% non-denaturing polyacrylamide gel followed by ethidium bromide staining to score the alleles.
To identify the causative genes, homozygosity mapping followed by Sanger sequencing was performed in four families (A–D). In family E, exome sequencing was used to search for the causative gene. Linkage in the four families (A–D) was tested using several microsatellite markers linked to 21 candidate genes (BBS1–BBS21) known for causing Bardet-Biedl Syndrome. After linkage was established in the families, the candidate regions were further saturated by typing additional markers linked to three genes (BBS7, BBS8, and BBS10).
In family E, DNA from a proband (IV-3) was subjected to whole exome sampling (WES; Figure 2A). Exome enrichment was accomplished using the SureSelect XT Human All Exon 50 Mb kit, version 5 (Agilent Technologies, San Clara, CA), and sequencing was performed on HiSeq 2500 systems (Illumina, San Diego, CA). All the reads were aligned against the human assembly hg19 (GRCh37), using Burrows-Wheeler Aligner (BWA v 0.7.5). Subsequently, the Exome Depth (v1.0.0), PINDEL (v 0.2.4t), and SAM tools (v 0.1.18) were used for variant calling. Filtering of the variants was performed with the help of the SAM tools varFilter script and custom scripts. All the variants obtained after filtering were inserted in an in-house database for further analysis.
After linkage was established in the four families (A–D), three genes (BBS7, BBS8, and BBS10) were Sanger sequenced. In family E, WES identified a pathogenic homozygous variant in the MKKS gene. All exons and flanking intronic sequences of the genes were PCR amplified using gene-specific primers. The primer sequences used to screen the BBS7 and BBS8 genes are listed in Table 1, while those used for BBS10 were the same as reported previously . The PCR products were purified using the Rapid PCR Purification System 9700 (Axygen, Union City, CA) and sequenced following the dideoxy chain termination method using the BigDye Terminator v3.1 Cycle Sequencing Kit and the Applied Biosystems 310 DNA Analyzer (Applied Biosystems Inc.) according to the manufacturer’s instructions. To identify sequence variants, the BioEdit sequence alignment editor version 6.0.7 was used. PCR primers were designed using the Primer3 program ()  and checked for specificity using BLAST.
All three affected members presented clinical features of obesity, retinitis pigmentosa, postaxial polydactyly (hexadactyly), mental delay, hypertension, and cutaneous syndactyly. Radiographs of the affected individual (IV-3) revealed normal carpals, metacarpals, and phalanges while the feet radiographs revealed postaxial polydactyly (PAP) type A, with an extra toe originating from a two-headed metatarsal. Tarsals and metatarsals were hypoplastic and distorted (Figure 3P,Q). Carpals were hypoplastic in the second affected member (IV-4) of the family (Figure 3R,S).
Fundus examination in two affected individuals (IV-3 and IV-4) showed typical features of sine pigmentosa, a variant of retinitis pigmentosa that have characteristic features, including atrophy of the RPE and vessel (arteriolar) attenuation. The vessels were thinner than normal, and the optic nerve had pallor. Features such as microaneurysms, edema, exudates, neovascularization, and hypo- and hyperpigmentation of macula were not observed in either affected individual (Figure 3W,X).
Affected individuals showed mental delay, hypertension, retinitis pigmentosa, reproductive tract/organ anomalies, obesity, postaxial polydactyly, and cutaneous syndactyly.
Affected members displayed weak analytical ability, vision impairment, obesity, hypogonadism, and postaxial polydactyly. Radiographs of both hands of affected individual (IV-6) displayed complete PAP type A. An extra digit, originated from the fifth metacarpal, was present in the left hand. In the right hand, the extra digit had no connection with the fifth metacarpal. Both extra digits had fix flexion deformity. Feet radiographs of the affected individuals (IV-1 and IV-6) showed PAP type A, and an extra toe originated from the fifth metatarsal (Figure 3T,U).
Both affected individuals (V-1 and V-2) showed clinical features of cognitive impairment, obesity, vision impairment, polysyndactylism, and hypogonadism. Camptodactyly was also observed in the right hand of affected individual V-2. Clinical reports of an affected individual (V-2) suggested the right kidney is either missing or severely hypoplastic, and the affected individual had a small penis with cryptorchidism. Presence of cognitive impairment was assessed through interviews with parents of affected members. In affected members of all four families (A–D), retinitis pigmentosa debuted with night blindness started at the age of 5 years with progressive decreased visual acuity at the age of 7–10 years.
Clinical evaluation of all three affected individuals in family E demonstrated common BBS phenotypes, such as obesity, learning disability, speech difficulties, mild hearing problems, and slight mental retardation. PAP was observed in the left hand and feet of affected individual IV-3 (Figure 2C–E) and the right hand of affected individual IV-4 (Figure 2G,H). An extra digit was surgically removed from affected individual III-7.
Selected microsatellite markers, mapped on chromosomal regions harboring genes for BBS1–BBS21, were used to search genetic linkage based upon the homozygosity mapping technique in four families (A–D). Haplotype analysis showed linkage in two families (A and B) to chromosome 4q27 harboring BBS7, in family C to chromosome 14q32 harboring BBS8, and in family D to chromosome 12q21.2 harboring BBS10. All affected members were heterozygous with microsatellite markers linked to other BBS genes.
In family E, direct WES was performed using DNA from an affected individual (IV-3) at the Institute of Human Genetics, Helmholtz Zentrum Munchen, Germany . All the variants were filtered and validated according to minor allele frequency (MAF) >0.001 in the Single Nucleotide Polymorphism database (dbSNP) and in the 1000 Genome Project, Exome Variant Server (EVS), and Exome Aggregation Consortium (ExAC), 7,000 in-house exome database (Appendix 1).
Sanger sequencing of the BBS7 gene identified a novel 3 bp deletion (c.580_582delGCA) in family A (Figure 1B) and a previously reported 6 bp in-frame deletion (c.1592_1597delTTCCAG) in family B (Figure 1D). In family C, sequence analysis of the BBS8 gene revealed a homozygous G to C transition at nucleotide position 1347(c.1347G>C) resulting in a missense variant (p.Gln449His) in all affected individuals (Figure 1F). In family D, sequence analysis of the BBS10 gene revealed a homozygous frameshift mutation (c.271_272insT; Figure 1H). In family E, a novel nonsense mutation (c.119C>G, p.Ser40*) was detected in the MKKS gene (Figure 2L). Unaffected members in the respective families were either heterozygous for a mutant allele or had a wild-type allele. Sanger sequencing validated cosegregation of the variants with the disease phenotype in all five families (A–E). The frequency of the identified variants in the in-house 7,000 exomes (IHG; Germany) and ExAc is presented in Appendix 2. The nonpathogenic nature of the four variants (c.580_582delGCA and c.1592_1597delTTCCAG in the BBS7 gene, c.1395G>C in the BBS8 gene, c.271_272insT in the BBS10 gene, and c.119C>G in the MKKS gene) were excluded in 175 ethnically matched control individuals.
Bardet-Biedl syndrome is a clinically pleiotropic disorder segregating in an autosomal recessive pattern . However, in rare cases, the triallelic nature of BBS involving three mutated alleles in two genes have been reported as well . In the study presented here, we identified disease-causing alleles in five consanguineous families (A–E) of Pakistani origin. The disease was inherited in an autosomal recessive pattern. Clinical features, observed in the five families, were similar to those reported previously . Linkage and WES analysis followed by Sanger sequencing revealed five variants, including two novel mutations in BBS6 (MKKS) and BBS7 and three previously reported mutations in the BBS7, BBS8, and BBS10 genes.
The BBS7 gene spans a 60.06 kb genomic region on chromosome 4q27. It is composed of 19 exons encoding 715 amino acids for the BBS7 protein. The BBS7 protein shares structural features with the BBS1 and BBS2 proteins . A group of seven BBS proteins (BBS1, BBS2, BBS4, BBS5, BBS7, BBS8, and BBS9) form a complex called the BBSome that promotes ciliogenesis. The cilia are microtubule-based structures that play important roles in the development of several tissues. Two of these proteins, BBS7 and BBS8, although structurally different, function as intraflagellar transport (IFT) cohesion factors . Three other proteins (BBS6, BBS10, and BBS12) function as chaperonins and arbitrate assembly of BBSome. Any sequence variant causing disruption in the seven genes results in defective ciliogenesis.
To date, 31 mutations, including 16 missense/nonsense, four splice sites, eight deletions, one small indel (insertion/deletion),
and one complex rearrangement, have been reported in the BBS7 gene
The BBS8 gene is one of the 21 genes identified for Bardet-Biedl syndrome. This gene encodes the 515 amino acid BBS8 protein that
interacts with PCM1 in ciliary biogenesis . Thus far, 15 mutations (including four missense/nonsense, five splice sites, five small deletions, and one small insertion)
have been reported in the BBS8 gene
BBS10 with two exons, encoding a 723 amino acid protein, mapped to chromosome 12q21.2 . Stoetzel et al. reported the most common mutation (271dupT, C91fsX95) in the BBS10 gene in several families with Bardet-Biedl syndrome. The same mutation was detected in family D in the present study . In a study involving five fetuses and a child, one fetus was found to be carrying homozygous 271dupT, 3 compound heterozygous with another mutation in the BBS10, and a fifth with homozygous 271dupT in addition to a truncating variant in the BBS6 gene . To date, 88 mutations, including 54 missense/nonsense, one splice site, 25 deletions, seven insertions, and one indel mutation, have been found in the BBS10 gene (HGMD) , which account for 20% of BBS cases.
MKKS/BBS6, mapped on chromosome 20p12, encodes a 570 amino acid protein. MKKS has been reported to cause the phenotypically overlapping McKusick-Kaufmann syndrome. Mutations in critical ciliary-regulating proteins result in an increase in rhodopsin in the inner segment and cause eventual photoreceptor cell death . This results in phenotypes that include tunnel vision, peripheral vision loss, and blindness. The BBS chaperonin complex is formed at the base of the primary cilia in the photoreceptor cell and comprises the three BBS proteins (MKKS/BBS6, BBS10, and BBS12) that form the BBSome complex [42,43]. Mice that lack Bbs6 and humans who have mutations in the MKKS/BBS6, BBS10, or BBS12 phenotypically resemble defects caused by BBSome genes . However, mutations in the MKKS/BBS6 gene give rise to variable phenotypes, yet similar to other BBS syndromes . The mutation (p.Ser40*) identified in family E resides within the predicted equatorial domain (Figure 2K) and most likely results in loss of function of the MKKS protein either through nonsense-mediated mRNA decay (NMD) or resulting in the production of a truncated MKSS protein. To date, 52 mutations have been identified in the MKKS gene, including 43 nonsense/missense, one splice site, six deletions, one insertion, and one indel (HGMD) .
BBS genes show expression in ciliated cells. Proteins, encoded by these genes, are divided into two groups, including BBSome, a complex formed by the assembly of BBS1, BBS2, BBS4, BBS5, BBS7, BBS8, and BBS9, involved in promoting ciliogenesis by recruiting Rab8GTP to enter the cilium and type II chaperonin-like proteins, a complex formed by BBS6, BBS10, and BBS12, which plays a role in the regulation of BBSome assembly . In conclusion, we have identified mutations in four BBS genes that cause BBS phenotypes. This study will support genetic testing of patients with BBS in Pakistan.
Appendix 1. Filtering steps followed to search for the candidate variant in affected individual of family E (IV-3).
Appendix 2. Frequency of the identified variants identifies in five families (A-E) in both in-house exoms and ExAc.
We highly appreciate cooperation and participation of the family members in this study. The work presented in the study was supported by Pakistan Academy of Science (PAS), Islamabad Pakistan. This work was funded by the Higher Education Commission (HEC), Islamabad, Pakistan. Muhammad Umair was supported by an Indigenous Ph.D. Fellowship and International Research Support Initiative Program (IRSIP) from HEC, Islamabad, Pakistan. Funding: The work presented in the study was supported by Pakistan Academy of Science (PAS), Islamabad Pakistan and funded by the Higher Education Commission (HEC), Islamabad, Pakistan.