Vitiligo is an autoimmune complex pigmentation disease characterized by non-pigmented patches on the surface of the skin that affect approximately 0.5–2% population worldwide. The exact etiology is still unknown; however, vitiligo is hypothesized to be a multifactorial and genetically heterogeneous condition. Therefore, the current study is designed to investigate the anthropometric presentation and genetic spectrum of vitiligo in fifteen consanguineous Pakistani families. The clinical evaluation of participating individuals revealed varying degrees of disease severity, with 23 years as the average age of disease onset. The majority of the affected individuals had non-segmental vitiligo (NSV). Whole exome sequencing analysis revealed clustering of rare variants of known vitiligo-associated genes. For instance, in the affected individuals of family VF-12, we identified three novel rare variants ofPTPN22(c.1108C>A),NRROS(c.197C>T) andHERC2(c.10969G>A) genes. All three variants replaced evolutionarily conserved amino acid residues in encoded proteins, which are predicted to impact the ionic interactions in the secondary structure. Although various in silico algorithms predicted low effect sizes for these variants individually, the clustering of them in affected individuals increases the polygenic burden of risk alleles. To our knowledge, this is the first study that highlights the complex etiology of vitiligo and genetic heterogeneity in multiplex consanguineous Pakistani families.