This study aimed to find the molecular basis of Bardet-Biedl syndrome (BBS) in Pakistani consanguineous families. A total of 12 affected families were enrolled. Clinical investigations were performed to access the BBS-associated phenotypes. Whole exome sequencing was conducted on one affected individual from each family. The computational functional analysis predicted the variants' pathogenic effects and modeled the mutated proteins. Whole-exome sequencing revealed 9 pathogenic variants in six genes associated with BBS in 12 families. The BBS6/MKS was the most common BBS causative gene identified in five families (5/12, 41.6%), with one novel (c.1226G>A, p.Gly409Glu) and two reported variants. c.774G>A, Thr259LeuTer21 was the most frequent BBS6/MMKS allele in three families 3/5 (60%). Two variants, c.223C>T, p.Arg75Ter and a novel, c. 252delA, p.Lys85STer39 were detected in the BBS9 gene. A novel 8bp deletion c.387_394delAAATAAAA, p. Asn130GlyfsTer3 was found in BBS3 gene. Three known variants were detected in the BBS1, BBS2, and BBS7 genes. Identification of novel likely pathogenic variants in three genes reaffirms the allelic and genetic heterogeneity of BBS in Pakistani patients. The clinical differences among patients carrying the same pathogenic variant may be due to other factors influencing the phenotype, including variants in other modifier genes.