Next Article in Journal
Association of Surfactant Protein D Single Nucleotide Polymorphisms rs721917, rs2243639, rs3088308 with Recurrent Aphthous Stomatitis in Pakistani Population
Previous Article in Journal
Evolutionary Adaptation of Genes Involved in Galactose Derivatives Metabolism in Oil-Tea Specialized Andrena Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 5
10.3390/genes14051118
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole Exome Sequencing Reveals Clustering of Variants of Known Vitiligo Genes in Multiplex Consanguineous Pakistani Families

by
Rafaqat Ishaq
1,2,
Muhammad Ilyas
1,3,
Umme Habiba
1,
Muhammad Noor ul Amin
1,
Sadia Saeed
1,4,
Ghazala Kaukab Raja
1,
Pakeeza Arzoo Shaiq
1,* and
Zubair M. Ahmed
2,5,6,*
1
University Institute of Biochemistry and Biotechnology, Pir Mehr Ali Shah Arid Agriculture University, Rawalpindi 46300, Pakistan
2
Department of Otorhinolaryngology-Head and Neck Surgery, School of Medicine, University of Maryland, Baltimore, MD 20742, USA
3
Department of Medical Laboratory Technology, Riphah International University, Malakand Campus, Malakand 23010, Pakistan
4
Department of Clinical Molecular Biology, EpiGen, Institute of Clinical Medicine, University of Oslo, 0313 Oslo, Norway
5
Department of Biochemistry and Molecular Biology, School of Medicine, University of Maryland, Baltimore, MD 20742, USA
6
Department of Ophthalmology and Visual Sciences, School of Medicine, University of Maryland, Baltimore, MD 20742, USA
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(5), 1118; https://doi.org/10.3390/genes14051118
Submission received: 13 April 2023 / Revised: 12 May 2023 / Accepted: 18 May 2023 / Published: 22 May 2023
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
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 of PTPN22 (c.1108C>A), NRROS (c.197C>T) and HERC2 (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.

1. Introduction

Vitiligo is a common skin pigmentation disorder with remarkable genetic and clinical heterogeneity. Vitiligo is characterized by the appearance of colorless patches on the epidermal surface of the skin. Discrepancies in the production and distribution of melanin in affected areas result in the onset of this socially stigmatizing and psychologically devastating disease. Vitiligo may affect one or both halves of the body. The vitiligo European task force assessment (VETFa) has broadly classified vitiligo into two main types, i.e., non-segmental vitiligo (NSV) or generalized vitiligo (GV) and Segmental Vitiligo (SV). Vitiligo is known to be progressive in most patients, which eventually results in complete loss of skin coloration, referred to as vitiligo Universalis (VU) in elderly patients. Various studies have reported that the worldwide incidence of vitiligo ranges between 0.5–2% [1,2]. However, there is no prior comprehensive study documenting the prevalence of vitiligo in the Pakistani population.
Vitiligo is a multifactorial, polygenic disorder that follows a non-Mendelian inheritance pattern with considerable incomplete penetrance due to complex genetic and environmental etiology [3]. A number of theories have been proposed to explain vitiligo pathogenesis that include genetics, autoimmune, auto-cytotoxic, and gene-environment interactions, which indicates that gene-environment interactions are predominantly involved in disease pathogenesis [4]. Multiplex vitiligo families partially support genetic hypotheses [5]. A report on the Korean population suggests a polygenic etiology of the disease [6]. The possible involvement of different sets of genes in melanin production, epistatic interactions, multistep regulation of the autoimmune system, and control of auto-toxicity has been implicated in disease pathogenesis [4]. The presence of certain environmental triggers, such as oxidative stress, ultraviolet radiation, and exposure to phenolic derivative compounds, may initiate vitiligo macule formation [7,8,9].
Advancement in genetic screening technologies has tremendously helped in the identification of rare variants associated with vitiligo. In the current study, we identified fifteen clinically diagnosed multiplex vitiligo families from the Punjab province of Pakistan. Whole exome sequencing was performed on seven of these families, which identified rare variants in vitiligo-associated genes. Our study also points towards the clinical and genetic complexity and locus heterogeneity of vitiligo in the inbred Pakistani population.

2. Materials and Methods

2.1. Family Enrollment and Clinical Classification

The proceedings of the current study were approved by committees of the Institutional Review Board (IRB) at the University of Maryland, School of Medicine, Baltimore, MD, USA (HP-00061036) and Pir Mher Ali Shah Arid Agriculture University Rawalpindi, Pakistan. The tenets of the Declaration of Helsinki were followed for the use of human subjects in the current study. Informed consent in writing was obtained from all participants. Fifteen vitiligo families were recruited from different regions of the Punjab province of Pakistan in this study.
Detailed family histories were obtained to document disease onset, affected area of skin, disease progression status, smoking/non-smoking status, chemical exposure, general lifestyle, eating habits, treatment etc. Family pedigrees were drawn based on the information received and verified by multiple participating family members.

2.2. Sample Collection and Isolation of DNA

From all participating individuals, 3 mL of venous blood samples were collected in an EDTA Orsin vacutainer (formerly imuMed®) and stored at 4 °C until further use. The genomic DNA was isolated from the blood samples of participants using the Thermo Scientific™ GeneJET genomic DNA purification kit (Catalog #: K0722). The concentration of DNA was quantified by NanoDrop® ND-1000 UV-Vis Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and DNA quality was assessed by 260/280 and 260/230 ratios.

2.3. Whole Exome Sequencing

Whole exome sequencing was performed on the probands of families VF-01, VF-02-, VF-03, VF-09, VF-10, VF-12 and VF-15. Genomic libraries were recovered for exome enrichment by using the Agilent SureSelect Human Expanded All Exon V5 kit and sequenced on an Illumina HiSeq4000 (Illumina, San Diego, CA, USA) with an average of 100× coverage [10]. Broad Institute’s Toolkit for Genome Analysis was used to analyze the sequencing data, as described previously [11,12], with the following modifications: (a) allele frequencies cut off was set to ≤0.003; (b) potential pathogenic prediction was achieved by at least 2 in silico algorithms; and (c) CADD score ≥ 6.0. Primers were designed using the online resource Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/, accessed on 23 October 2018) for the variants that passed our filtration criteria; targeted DNA amplification by conventional PCR followed by Sanger sequencing was carried out for segregation studies [13].

2.4. Molecular Modeling

The Phyre2 (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index, accessed on 15 March 2019) and protein variation predictor resource HOPE (https://www3.cmbi.umcn.nl/hope/, accessed on 15 March 2019) were applied for 3D modeling of protein structures. Visualization of the mutation impact on the 3D structure and change of ionic interaction between the residues of interest was visualized by the UCSF CHIMERA tool (https://www.cgl.ucsf.edu/chimera/, accessed on 19 March 2019). MUpro (http://mupro.proteomics.ics.uci.edu/, accessed on 27 March 2023) and I-Mutant 2.0 (http://folding.biofold.org/i-mutant/i-mutant2.0.html, accessed on 27 March 2023) were used to predict the effects of these mutations on protein structure stability on the bases of changes in Gibbs free energy (G). The evolutionary conservation status of multiple amino acids across different species was performed by Clustal Omega multiple sequence alignment (https://www.ebi.ac.uk/Tools/msa/clustalo/, accessed on 10 March 2019).

3. Results

3.1. Clinical Evaluation of Vitiligo in Multiplex Pakistani Families

After IRB approval, a total of 15 multiplex consanguineous families segregating vitiligo (Figure 1) were enrolled in this study from five different cities (Rawalpindi, Chakwal, Mianwali, Abbottabad, Islamabad) in Pakistan (Table 1). A total of 130 individuals from these families participated in our study, including 72 males (55%) and 58 females (45%). Of these enrolled subjects, 53 (33 males and 20 females) exhibited the vitiligo phenotype, while 77 individuals had no skin problem and, thus, were used as controls in our downstream genetic studies. The age range for onset vitiligo macules was between 2 to 68 years, with an average of 23 years. Among the enrolled females, non-segmental vitiligo (NSV) was the most noticeable clinical phenotype noted in forty-eight females (91%), three females (6%) had vitiligo universalis (VU), while two females (3%) reported complete recovery from their vitiligo phenotype. Almost all affected individuals reported having no other autoimmune diseases or any other comorbidity besides vitiligo, except two affected individuals that had type 1 diabetes mellitus (Table 1).
Inter- and intra-familial variability in the areas of skin discoloration among the affected individuals was identified (Table 1). For instance, family VF-02 is an extended vitiligo family having nine affected individuals (Figure 1). Three of the affected individuals (VF02_15, _16, _17) had 90–100% of their body with vitiligo patches and thus present VU (Table 1), while the other affected individuals had much fewer patches (Table 1). Intriguingly, according to family history interviews, one affected female (VH02_07) recovered from vitiligo without any medication. In contrast, the three affected individuals of family VF-04 had relatively less progressive disease severity (Table 1).

3.2. Genetic Analysis of Families with Vitiligo

To decipher the genetic risk factors contributing to the vitiligo phenotype in the enrolled families, we performed whole exome sequencing on index patients from seven families along with ethnically matched normal control samples. After quality control, we used in silico panel to filter variants in the known vitiligo-associated genes [14]. Next, we selected variants of the genes that had (a) an allele frequency of ≤0.001; (b) a CADD score of ≥7.0; and (c) predicted damage by at least one in silico algorithm. Table 2 lists all variants of known vitiligo genes that passed these criteria. Although various in silico algorithms predicted low effect size for these variants individually, the clustering of these variants in affected individuals increases the polygenic burden of risk alleles (Table 2); we found potentially pathogenic variants of three known vitiligo-associated genes, PTPN22, HERC2 and NRROS, in a varying combination of zygosity among affected individuals of family VF-12 (Figure 2). Proband III-2, having the most severe phenotype, was found heterozygous for the rare missense variants (Figure 2) of PTPN22 (c.1108C>A; p.(His370Asn)), HERC2 (c.10969G>A; p. (Val3657Ile)), and NRROS (c.197C>T; p.(Ala66Val)). Different family members inherited different combinations of these potentially pathogenic variants, and hence disease severity and allelic burden vary among these individuals (Figure 2 and Table 1). Variants found in the PTPN22, HERC2 and NRROS genes have been previously known for their association with NSV in different ethnicities around the world [15,16,17,18].
Amino acid conservation studies have shown that the variants identified in PTPN22, HERC2 and NRROS genes are present in the evolutionarily conserved regions of the encoded proteins (Figure 2C). Next, we performed molecular modeling analysis using the Phyre2 online resource for retrieving structural information based on protein homology present in the protein database. Using this approach, we were able to model the two missense variants found in the HERC2 and NROSS proteins (Figure 3). The lack of available crystalline structural information for PTPN22 restrained further analysis from predicting the effect of the mutation on protein structure and function. A different pattern of interresidual bonding interaction between wild type and mutant protein was observed for the p.Val 3657Ile variant of HERC2 (Figure 3A). The modeling data predicted that p.Val3657Ile substitution might impact protein folding due to the formation of an additional H-bond between p.Ile3657 and p.Leu3661 (Figure 3A). Similarly, the p.Ala66Val variant of NRROS is predicted to greatly hamper the microscopic interaction in the surrounding area of 10 Å by differential H-bonding patterns (Figure 3B). Surrounding the p.Ala66 residue in WT NRROS, there are four prominent hydrogen bonds: between p.Gly45 and p.Asp65 (2.837 Å); p.Cys91 and Leu64 (0.3129 Å); p.Cys91 and p.His89 (0.374 Å); and p.Ala66 and p.His89 (2.921 Å; Figure 3B). In contrast, the p.Ala66Val substitution exhibits a different set of H-bonding capabilities, resulting in the loss of H-bonds between p.Gly45 and p.Asp65 and the formation of two new H-bond bonds between p.Val66 and p.His89 (2.962 Å), and between p.Gly113 and p.Ser90 (2.293 Å; Figure 3B). This new pattern of interaction may induce torsion, leading to improper protein folding and properties, which may have an effect on function. A noticeable decrease in Gibbs free energy was recorded in PTPN22 (−0.73516404), HERC2 (−0.52853749) and NRROS (−0.22229897) by MUpro, indicating a decrease in protein 3D structure stability, as corroborated by the I-Mutant prediction program. Overall, the three identified variants have been anticipated to have deleterious effects on protein structures, rendering the proteins unstable (Figure 3C).
Using this in silico approach, we also modeled six missense variants found in the affected individuals of our vitiligo families (Table 2) and predicted their potential pathogenic impact on the encoded proteins (Figure 4). These include p.Asn201Thr (GSTP1), p.Arg63Gln (CDH1), p.Arg161Trp (CASP7), p.Pro695Leu (NOS2), p.Glu518Lys (UBASH3A), and p.Asp442Gly (PTPRC) variants (Figure 4). The p.Asn201Thr variant of GSTP1 is located within the GST C-terminal domain, and due to different sizes and ionic properties, it might impact the folding of the domain (Figure 4). A similar impact was predicted for the p.Arg63Gln variant of CDH1, p.Arg161Trp of CASP7, and p.Asp442Gly of PTPRC (Figure 4). The p.Pro695Leu of NOS2 is predicted to form a new H-bond with p.Ile692 (Figure 4). Moreover, the WT Proline residue is known to be very rigid and therefore induce a special backbone conformation; therefore, substitution with Leucine would affect this conformation. Finally, the p.Glu518Lys variant of UBASH3A is predicted to alter the salt bridge with Arginine at position 406 and Arginine at position 485 (Figure 4), and thus it ultimately affects protein folding and function. Overall, the above missense variants are predicted to have an impact on protein folding and secondary structure, which might lead to impaired function.

4. Discussion

Vitiligo is a disease of complex etiology, whereby the involvement of genetic contributors in disease pathogenesis has been emphasized greatly since the advent of modern genetic technologies. Previously, several genetic association studies have identified risk alleles responsible for various types of vitiligo in different ethnicities [17,19,20,21,22,23,24]. In the present study, we reported clinical and genetic findings in fifteen highly inbred consanguineous Pakistani families. Our study highlights the genetic heterogeneity of the disease in terms of risk allele clustering as well as associated risk allele burden in affected individuals. We report the identification of several rare variants of known vitiligo-associated genes, including three novel missense variants in PTPN22 (c.1108C>A; p.(His370Asn)), NRROS (c.197C>T; p.(Ala66Val)) and HERC2 (c.10969G>A; p. (Val3657Ile)) in family VF-12. In silico analysis of the risk variants for pathogenicity determination predicted a low to moderate pathogenic impact of all the identified variants; however, clustering of mutated alleles in an individual collectively contributed to the onset of vitiligo.
Autoimmune-mediated melanocyte destruction is a common phenomenon in vitiligo pathogenesis that involves intrinsic defects within melanocytes and an autoimmune system that targets these cells [25]. PTPN22, which encrypts genetic information for lymphoid-specific intracellular phosphatase, is the main inhibitor of T-cell activation [19]. Many autoimmune diseases, including type 1 diabetes, arthritis, systemic lupus and vitiligo, have also been previously associated with PTPN22 variants [4]. NSV was found to be strongly associated with PTPN22 (c.1858C>T; p. (Arg620Trp)) in the cohorts from the United Kingdom, Gujarat, India and Romania in three different case-control studies [15,26]. Similarly, HERC2, which encodes for the E3 ubiquitin-protein, provides protection against ionizing radiation and DNA damage related to radiation. Previous studies have shown that an elevated risk of vitiligo is associated with HERC2 variations, which have a significant impact on determining skin and iris color besides the OCA2. Previous studies on Caucasian and Chinese Han cohorts also reported a significant association between the T allele of rs6583331 SNP of NRROS and vitiligo [17]. However, in our cohort, a novel vitiligo-associated variant of NRROS was identified. Our outcomes are in agreement with previous findings, thus indicating an increased risk for intra-familial NSV development linked to NRROS.
Besides the above-mentioned genes, as shown in Table 2, the majority of the variants identified in our cohort are present in proteins that function in innate immune responses, inflammation reaction, control of cell apoptosis, anti-oxidative damage, and have a role in detoxification reactions. For instance, IFIH1 is an interferon-induced helicase C domain-containing protein, an innate immune receptor that acts as a cytoplasmic sensor of viral nucleic acids and plays a major role in sensing viral infection and in the activation of a cascade of antiviral responses, including the induction of type I interferons and pro-inflammatory cytokines [25]. The IL1RAPL1-Interleukin-1 receptor accessory protein-like 1 regulates secretion and presynaptic differentiation through inhibition of the activity of the N-type voltage-gated calcium channel and is also involved in the activation of the MAP kinase JNK. TLR4 is a transmembrane protein, a member of the Toll-like receptor family, which belongs to the pattern recognition receptor (PRR) family. TLR4 participates in intracellular signaling NF-κB and inflammatory cytokine production, which is responsible for activating the innate immune system. CASP7 encodes a member of the Cysteine-Aspartic acid protease (caspase) family and plays a central role in the execution phase of cell apoptosis and survival of melanocytes [4]. Similarly, the BCL2L12 protein family and their interactors include inhibitors and inducers of cell death, regulate and mediate the process by which mitochondria contribute to cell death known as the intrinsic apoptosis pathway, and thus regulate the survival of pigmented cells.
In conclusion, the convergence theory of vitiligo development bridges all previously known etiologies and describes vitiligo as a multivariate, complex, and multifactorial disorder [9]. Vitiligo is not caused merely by the presence of an underlying genetic variation or only by exposure to some specific external factor such as stimuli or defective melanocytes and an internal autoimmune condition but is a result of an intricate combination of all these etiological contributors. Therefore, it is important to acknowledge that the underlying molecular pathways affected in one patient do not always lead to the same level of pathogenicity in other individuals. Such attributes of vitiligo define that identification of individual-specific pathways is significant for the targeted development of more effective therapeutics [26]. Polygenic Risk Score (PRS) estimation is a trending approach for the genetic description of complex disorders, including vitiligo. However, there always remains quite a fair chance for the involvement of additional genes or novel variants in the known vitiligo-associated genes. In this regard, it has been recently described that the phenomenon of clustering of risk alleles in an inbred family increases polygenic risk burden [5], and our study further supports this notion.

Author Contributions

Conceptualization, P.A.S. and Z.M.A.; methodology, P.A.S.; software, R.I. and U.H.; validation, U.H., M.I. and R.I.; formal analysis, R.I.; investigation, R.I. and M.N.u.A.; resources, Z.M.A., S.S. and G.K.R.; data curation, R.I. and U.H.; writing—original draft preparation, R.I. and Z.M.A.; writing—review and editing, R.I. and Z.M.A.; visualization, R.I.; supervision, P.A.S. and Z.M.A.; project administration, P.A.S.; funding acquisition, P.A.S. and Z.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially supported by International Research Support Initiative Program, Higher Education Commission of Pakistan and by National Institute of Arthritis and Musculoskeletal and Skin Disease, National Institutes of Health, Bethesda, MD, USA (R01AR077563 to ZMA).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board/Ethics committee of Pir Mehr Ali Shah Arid Agriculture University, Rawalpindi (No. PMAS-AAUR/IEC/15 dated 11 October 2017).

Informed Consent Statement

Informed consent was obtained from all subjects included in the study.

Data Availability Statement

The data presented in this study are available on request from Corresponding authors. The data is not publically available due to confidentiality of patient’s information.

Acknowledgments

We are thankful to all the participants of the study and health care professionals for their cooperation in clinical investigations. We are very thankful to Muhammad Naeem, Fahmida (Quaid-e-Azam University, Islamabad, Pakistan) and Ehsanullah (National Institute of eye, NIH, USA) for their technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berqvist, C.; Ezzedine, K. Vitiligo: A Review. Dermatology 2020, 236, 571–592. [Google Scholar] [CrossRef]
  2. Garcia-Melendez, M.E.; Salinas-Santander, M.; Sanchez-Dominguez, C.; Gonzalez-Cardenas, H.; Cerda-Flores, R.M.; Ocampo-Candiani, J.; Ortiz-López, R. Protein tyrosine phosphatase PTPN22 +1858C/T polymorphism is associated with active vitiligo. Exp. Med. 2014, 8, 1433–1437. [Google Scholar] [CrossRef] [PubMed]
  3. Spritz, R.A. Recent progress in the genetics of generalized vitiligo. J. Genet. Genom. 2011, 38, 271–278. [Google Scholar] [CrossRef] [PubMed]
  4. Czajkowski, R.; Mecinska-Jundzill, K. Current aspects of vitiligo genetics. Postep. Derm. Alergol. 2014, 31, 247–255. [Google Scholar] [CrossRef] [PubMed]
  5. Roberts, G.H.L.; Paul, S.; Yorgov, D.; Santorico, S.A.; Spritz, R.A. Family Clustering of Autoimmune Vitiligo Results Principally from Polygenic Inheritance of Common Risk Alleles. Am. J. Hum. Genet. 2019, 105, 364–372. [Google Scholar] [CrossRef]
  6. Kim, S.M.; Chung, H.S.; Hann, S.-K. The genetics of vitiligo in Korean patients. Int. J. Dermatol. 1998, 37, 908–910. [Google Scholar] [CrossRef]
  7. Laddha, N.C.; Dwivedi, M.; Mansuri, M.S.; Gani, A.R.; Ansarullah, M.; Ramachandran, A.V.; Dalai, S.; Begum, R. Vitiligo: Interplay between oxidative stress and immune system. Exp. Derm. 2013, 22, 245–250. [Google Scholar] [CrossRef]
  8. Yu, S.; Doycheva, D.M.; Gamdzyk, M.; Yang, Y.; Lenahan, C.; Li, G.; Li, D.; Lian, L.; Tang, J.; Lu, J.; et al. Activation of MC1R with BMS-470539 attenuates neuroinflammation via cAMP/PKA/Nurr1 pathway after neonatal hypoxic-ischemic brain injury in rats. J. Neuroinflamm. 2021, 18, 26. [Google Scholar] [CrossRef]
  9. Alikhan, A.; Felsten, L.M.; Daly, M.; Petronic-Rosic, V. Vitiligo: A comprehensive overview Part I. Introduction, epidemiology, quality of life, diagnosis, differential diagnosis, associations, histopathology, etiology, and work-up. J. Am. Acad. Derm. 2011, 65, 473–491. [Google Scholar] [CrossRef]
  10. Noman, M.; Ishaq, R.; Bukhari, S.A.; Ahmed, Z.M.; Riazuddin, S. Delineation of Homozygous Variants Associated with Prelingual Sensorineural Hearing Loss in Pakistani Families. Genes 2019, 10, 1031. [Google Scholar] [CrossRef]
  11. Richard, E.M.; Santos-Cortez, R.L.P.; Faridi, R.; Rehman, A.U.; Lee, K.; Shahzad, M.; Acharya, A.; Khan, A.A.; Imtiaz, A.; Chakchouk, I.; et al. Global genetic insight contributed by consanguineous Pakistani families segregating hearing loss. Hum. Mutat. 2019, 40, 53–72. [Google Scholar] [CrossRef] [PubMed]
  12. Riazuddin, S.; Hussain, M.; Razzaq, A.; Iqbal, Z.; Shahzad, M.; Polla, D.L.; Song, Y.; van Beusekom, E.; Khan, A.A.; Tomas-Roca, L.; et al. Exome sequencing of Pakistani consanguineous families identifies 30 novel candidate genes for recessive intellectual disability. Mol. Psychiatry 2017, 22, 1604–1614. [Google Scholar] [CrossRef] [PubMed]
  13. Jaworek, T.J.; Kausar, T.; Bell, S.M.; Tariq, N.; Maqsood, M.I.; Sohail, A.; Ali, M.; Iqbal, F.; Rasool, S.; Riazuddin, S.; et al. Molecular genetic studies and delineation of the oculocutaneous albinism phenotype in the Pakistani population. Orphanet. J. Rare Dis. 2012, 7, 44. [Google Scholar]
  14. Debasis, d.; Vivek; Manchanda; Singh, P.; Ankita; Gupta, I. A Compendium of Genes and Variations Associated with Vitiligo. 2023. Available online: https://vitivar.igib.res.in/genes.php (accessed on 12 March 2023).
  15. Canton, I.; Akhtar, S.; Gavalas, N.G.; Gawkrodger, D.J.; Blomhoff, A.; Watson, P.F.; Weetman, A.P.; Kemp, E.H. A single-nucleotide polymorphism in the gene encoding lymphoid protein tyrosine phosphatase (PTPN22) confers susceptibility to generalised vitiligo. Genes Immun. 2005, 6, 584–587. [Google Scholar] [CrossRef] [PubMed]
  16. Passeron, T.; Ortonne, J.P. Activation of the unfolded protein response in vitiligo: The missing link? J. Investig. Derm. 2012, 132, 2502–2504. [Google Scholar] [CrossRef]
  17. Tang, X.; Fang, F.; Yang, J.; Zheng, X.; Fan, M.; Wang, L.; Zhang, A. Association Study Reveals One Susceptibility Locus with Vitiligo in the Chinese Han Population. Genet. Test. Mol. Biomark. 2019, 23, 791–796. [Google Scholar] [CrossRef]
  18. Liu, J.; Zhang, Z.; Chai, L.; Che, Y.; Min, S.; Yang, R. Identification and characterization of a unique leucine-rich repeat protein (LRRC33) that inhibits Toll-like receptor-mediated NF-κB activation. Biochem. Biophys. Res. Commun. 2013, 434, 28–34. [Google Scholar] [CrossRef]
  19. Rajendiran, K.S.; Rajappa, M.; Chandrashekar, L.; Thappa, D.M. Association of PTPN22 gene polymorphism with non-segmental vitiligo in South Indian Tamils. Postep. Derm. Alergol. 2018, 35, 280–285. [Google Scholar] [CrossRef]
  20. Cheong, K.A.; Kim, N.H.; Noh, M.; Lee, A.Y. Three new single nucleotide polymorphisms identified by a genome-wide association study in Korean patients with vitiligo. J. Korean Med. Sci. 2013, 28, 775–779. [Google Scholar] [CrossRef]
  21. Rodrigues, M.; Ezzedine, K.; Hamzavi, I.; Pandya, A.G.; Harris, J.E.; Vitiligo Working Group. New discoveries in the pathogenesis and classification of vitiligo. J. Am. Acad. Dermatol. 2017, 77, 1–13. [Google Scholar] [CrossRef]
  22. Laddha, N.C.; Dwivedi, M.; Shajil, E.M.; Prajapati, H.; Marfatia, Y.S.; Begum, R. Association of PTPN22 1858C/T polymorphism with vitiligo susceptibility in Gujarat population. J. Derm. Sci. 2008, 49, 260–262. [Google Scholar] [CrossRef] [PubMed]
  23. Bekker-Jensen, S.; Rendtlew Danielsen, J.; Fugger, K.; Gromova, I.; Nerstedt, A.; Lukas, C.; Bartek, J.; Lukas, J.; Mailand, N. HERC2 coordinates ubiquitin-dependent assembly of DNA repair factors on damaged chromosomes. Nat. Cell. Biol. 2010, 12, 80–86. [Google Scholar] [CrossRef] [PubMed]
  24. Jin, Y.; Birlea, S.A.; Fain, P.R.; Ferrara, T.M.; Ben, S.; Riccardi, S.L.; Cole, J.B.; Gowan, K.; Holland, P.J.; Bennett, D.C.; et al. Genome-wide association analyses identify 13 new susceptibility loci for generalized vitiligo. Nat. Genet. 2012, 44, 676–680. [Google Scholar] [CrossRef] [PubMed]
  25. Jin, Y.; Andersen, G.H.L.; Santorico, S.A.; Spritz, R.A. Multiple Functional Variants of IFIH1, a Gene Involved in Triggering Innate Immune Responses, Protect against Vitiligo. J. Invest. Derm. 2017, 137, 522–524. [Google Scholar] [CrossRef]
  26. Kundu, R.V.; Mhlaba, J.M.; Rangel, S.M.; Le Poole, I.C. The convergence theory for vitiligo: A reappraisal. Exp. Dermatol. 2019, 28, 647–655. [Google Scholar] [CrossRef] [PubMed]
Genes 14 01118 g001 550
Figure 1. Pedigrees of vitiligo families. Black-filled and empty symbols represent affected and unaffected individuals, respectively. The double line shows only consanguineous marriages. The arrow sign indicate probands that were used for whole exome sequencing. Clinical manifestations of the vitiligo patches for representative affected members are also shown below each pedigree.
Figure 1. Pedigrees of vitiligo families. Black-filled and empty symbols represent affected and unaffected individuals, respectively. The double line shows only consanguineous marriages. The arrow sign indicate probands that were used for whole exome sequencing. Clinical manifestations of the vitiligo patches for representative affected members are also shown below each pedigree.
Genes 14 01118 g001
Genes 14 01118 g002 550
Figure 2. Pedigree of vitiligo family and identified causative mutation. (A) Polygenic pattern of disease-causing alleles segregating in multiplex Pakistani vitiligo family. Black-filled and empty symbols represent affected and unaffected individuals, respectively. Symbols filled in gray represent individuals that recovered (<5%) from vitiligo. Double line shows only consanguineous marriages. The genotypes (wild type, heterozygous) of the identified mutant alleles are also shown for each of the participating family members. Arrow sign represents the proband used for whole exome sequencing. (B) Vitiligo patches were observed in the affected individuals of VF-12 family. (C) Clustal-W multiple amino acid sequence alignments of orthologous proteins show evolutionarily conserved mutated residues across different species.
Figure 2. Pedigree of vitiligo family and identified causative mutation. (A) Polygenic pattern of disease-causing alleles segregating in multiplex Pakistani vitiligo family. Black-filled and empty symbols represent affected and unaffected individuals, respectively. Symbols filled in gray represent individuals that recovered (<5%) from vitiligo. Double line shows only consanguineous marriages. The genotypes (wild type, heterozygous) of the identified mutant alleles are also shown for each of the participating family members. Arrow sign represents the proband used for whole exome sequencing. (B) Vitiligo patches were observed in the affected individuals of VF-12 family. (C) Clustal-W multiple amino acid sequence alignments of orthologous proteins show evolutionarily conserved mutated residues across different species.
Genes 14 01118 g002
Genes 14 01118 g003 550
Figure 3. Three-dimensional protein structural modeling encoded from variants identified in the current study. Structures of both wild type (shown in green) and mutated residues (shown in red) are highlighted. Solid black lines represent the interactions with neighboring amino acids. Stick representation in green shows wild type residue, while red color shows mutated amino acid. Helices are shown in magenta, sheets in light green and loops in cyan color. (A) Missense variant in HERC2 p.(Val3657Ile), wild type residue is relatively smaller in size than the mutant. Mutant residue p.(Ile 3657), shown in red, forms H-bonds with neighboring p.(Leu 3661), shown in blue, and bond length is 0.2561 Å, which may have an impact on protein folding and functionality. (B) Similarly, for missense variant of NRROS p.(Ala66Val), wild type residue p.(Ala66), shown in green, forms H-bonds with p.(His89), shown in blue. Neighboring residue, shown in yellow with solid black line, represents H-bonds in the wild type. Mutant model represents different sets of interactions in the adjacent region of mutation. p.(Val66), shown in red, forms H-bonds with p.(His89), shown in blue, with different bond lengths causing loss of H-bonding between p.(Gly45) and p.(Asp65), shown in yellow, in wild type. Formation of new intra residual H-bonding between p.(Ser90) and p.(Gly113) is shown in violet color. (C) The underlying mutations negatively impacted protein structures, and the encoded proteins were predicted to become unstable by MUpro and I-mutant.
Figure 3. Three-dimensional protein structural modeling encoded from variants identified in the current study. Structures of both wild type (shown in green) and mutated residues (shown in red) are highlighted. Solid black lines represent the interactions with neighboring amino acids. Stick representation in green shows wild type residue, while red color shows mutated amino acid. Helices are shown in magenta, sheets in light green and loops in cyan color. (A) Missense variant in HERC2 p.(Val3657Ile), wild type residue is relatively smaller in size than the mutant. Mutant residue p.(Ile 3657), shown in red, forms H-bonds with neighboring p.(Leu 3661), shown in blue, and bond length is 0.2561 Å, which may have an impact on protein folding and functionality. (B) Similarly, for missense variant of NRROS p.(Ala66Val), wild type residue p.(Ala66), shown in green, forms H-bonds with p.(His89), shown in blue. Neighboring residue, shown in yellow with solid black line, represents H-bonds in the wild type. Mutant model represents different sets of interactions in the adjacent region of mutation. p.(Val66), shown in red, forms H-bonds with p.(His89), shown in blue, with different bond lengths causing loss of H-bonding between p.(Gly45) and p.(Asp65), shown in yellow, in wild type. Formation of new intra residual H-bonding between p.(Ser90) and p.(Gly113) is shown in violet color. (C) The underlying mutations negatively impacted protein structures, and the encoded proteins were predicted to become unstable by MUpro and I-mutant.
Genes 14 01118 g003
Genes 14 01118 g004 550
Figure 4. WT residue is shown in dodger blue color, mutant in firebrick (red) and hydrogen-bond (orange color line) associated amino acids in dark magenta. Amino acid atoms are shown in ball and stick forms. Secondary structure: Helix: Green strands: pink-purple and coils: golden-yellow in color.
Figure 4. WT residue is shown in dodger blue color, mutant in firebrick (red) and hydrogen-bond (orange color line) associated amino acids in dark magenta. Amino acid atoms are shown in ball and stick forms. Secondary structure: Helix: Green strands: pink-purple and coils: golden-yellow in color.
Genes 14 01118 g004
Table
Table 1. Clinical phenotype of Vitiligo families.
Table 1. Clinical phenotype of Vitiligo families.
Family IDFamily EthnicityIndividual IDAge (Year)Age of OnsetGenderPhenotypeDoctoral Description of the DiseaseAffected Area of Skin (%)Comorbidity
VF-01RajputVF01_012818MaleAffectedNSV, Generalized vitiligo20–30%Nil
VF01_022621FemaleAffectedNSV, Generalized Vitiligo20–30%Nil
VF01_041714FemaleAffectedNSV, Generalized vitiligo10–20%Nil
VF-02MughalVF02_032820FemaleAffectedNSV, Generalized vitiligo10–20%Nil
VF02_061818FemaleAffectedNSV, Generalized VitiligoLess than 5%Nil
VF02_074239MaleRecoveredGeneralized Vitiligo (NSV*)N/ANil
VF02_113217MaleAffectedNSV, Generalized vitiligo 30–40%Nil
VF02_158014FemaleAffectedVitiligo Universalis 90–100%Diabetes
VF02_167919MaleAffectedNSV, Generalized vitiligo 90–100%Nil
VF02_175618FemaleAffectedNSV, Generalized vitiligo 90–100%Nil
VF02_182311FemaleAffectedNSV, Generalized vitiligo 10–20%Nil
VF02_192010FemaleRecoveredVitiligo recovered (NSV*)N/ANil
VF-03AwanVF03_014624MaleAffectedNSV, Generalized vitiligo 30–40%Nil
VF03_07704MaleAffectedNSV, Acrofacial vitiligo 5–10%Nil
VF03_083027FemaleAffectedNSV, Acrofacial vitiligo Less than 5%Nil
VF-04SardarVF04_013024MaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF04_023831MaleAffectedNSV, Generalized vitiligo Less than 5%Nil
VF04_03137FemaleAffectedNSV, Generalized vitiligo Less than 5%Nil
VF-05MughalVF05_014733MaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF05_042221MaleAffectedVitiligo started in neck region Less than 5%Nil
VF-06RajputVF06_013111MaleAffectedNSV, Generalized vitiligo 70–80%Nil
VF06_022513MaleAffectedNSV, Generalized vitiligo 40–50%Nil
VF-07NiaziVF07_015523MaleAffectedNSV, Generalized vitiligo 40–50%Nil
VF07_032119MaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF-08ShiekhVF08_016919FemaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF08_042117MaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF-09AwanVF09_016459FemaleAffectedNSV, Generalized vitiligo20–30%Nil
VF09_051715MaleAffectedNSV, Generalized vitiligo Less than 5%Nil
VF09_071915FemaleAffectedNSV, Generalized vitiligo Less than 5%Nil
VF09_1092FemaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF-10MalikVF10_017933MaleAffectedNSV, Generalized vitiligo 60–70%Nil
VF10_048342MaleAffectedNSV, Generalized vitiligo 70–80%Nil
VF10_078040FemaleAffectedNSV, Generalized vitiligo 40–50%Nil
VF10_103527MaleAffectedNSV, Generalized vitiligo 40–50%Nil
VF10_115133MaleAffectedNSV, Generalized vitiligo50–60%Nil
VF10_14127MaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF10_162817MaleAffectedNSV, Generalized vitiligo 40–50%Nil
VF-11QaziVF11_015611FemaleAffectedVitiligo Universalis 90–100%Nil
VF11_025853MaleAffectedVitiligo just started at the skull Less than 1%Nil
VF11_031911MaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF-12MalikVF12_015743MaleAffectedNSV, Generalized vitiligo 30–40%Diabetes
VF12_052615MaleAffectedNSV, Generalized vitiligo 30–40%Nil
VF12_062322MaleAffectedNSV, Generalized vitiligo Less than 5%Nil
VF12_073524MaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF12_097968FemaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF-13RajputVF13_017466MaleAffectedNSV, Generalized vitiligo 30–40%Nil
VF13_04139MaleAffectedNSV, Generalized vitiligo 5–10%Nil
VF-14ChaudaryVF14_016731FemaleAffectedNSV, Generalized vitiligo 15–20%Nil
VF14_034738MaleAffectedNSV, Generalized vitiligo 30–40%Nil
VF14_04159FemaleAffectedNSV, Generalized vitiligo 10–20%Nil
VF-15RajputVF15_017719MaleAffectedVitiligo Universalis90–100%Nil
VF15_042114MaleAffectedNSV, Generalized vitiligo 20–30%Nil
VF15_072113MaleAffectedNSV, Generalized vitiligo 30–40%Nil
NSV: Non segmental vitiligo; NSV*: Non-segmental vitiligo (Recovered); N/A: Not applicable.
Table
Table 2. Genetic analysis of index patients for known vitiligo gene variants.
Table 2. Genetic analysis of index patients for known vitiligo gene variants.
FamilyChrGeneRefAltc.DNA ChangeProtein ChangeSNP IDrs NumberExACCADDSIFTPolyphenMutation
Taster		Mutation
Assessor		Genotyping of Proband
VF-012IFIH1TCc.A1046Gp.K349RNM_022168rs726506640.003514.59TBDNHet
11GSTP1ACc.A601Cp.N201HNM_000852rs145957405-12.59DPDMHet
11GSTP1ACc.A602Cp.N201TNM_000852rs3760742800.000016777.554TBDMHet
12ATXN2TGC-c.534_536delp.178_179delNM_002973rs1939229270-----Het
15HERC2CT-c.5780_5781dep.K1927fsNM_004667-------Het
16RPGRIP1LTAc.A278Tp.E93VNM_015272rs7658768390.0000082518.56TPDLHet
XIL1RAPL1TG-c.63_64delp.V21fsNM_014271rs147274241......Het
VF-0219BCL2L12GAc.G542Ap.R181QNM_138639rs7622862810.000319.56TDDLHet
VF-038TGCTc.C3902Tp.P1301LNM_003235rs5491842030.000221.1DDDMHet
11GSTP1AGc.A1Gp.M1VNM_000852rs1448870282-10.05DB--Het
VF-093GPX1CAc.G442Tp.A148SNM_001329455rs6446261-6.871TBDLHet
3GPX1GCC-c.33_35delp.11_12delNM_001329455rs17838762------Het
9TLR4AGc.A1084Gp.K362ENM_003266rs5391537080.000212.84TPNNHet
10CASP7CTc.C481Tp.R161WNM_033340rs18564998200TBN-Het
21UBASH3AGAc.G1666Ap.E556KNM_018961rs1167523270.000115.6DDDMHet
VF-101PTPRCTAc.T566Ap.I189NNM_002838rs2017151570.00137.444DBNLHet
2FARP2AGc.A134Gp.H45RNM_014808rs1503124580.00348.045TBNNHet
8TGGAc.G2963Ap.R988HNM_003235rs168933320.006211.02TBPNHet
14STRN3GTc.C53Ap.P18HNM_001083893rs1036923688.15.33DBDNHet
16CDH1GAc.G188Ap.R63QNM_004360rs5877801170.0000576512.03TBNMHet
17NOS2GAc.C2084Tp.P695LNM_000625rs7678073660.0000495419.85TPDMHet
VF-121PTPN22GTc.C1108Ap.H370NNM_015967rs726506710.00236.752TPNMHet
10ACTA2CT-Splicing -rs112687898-26--D-Het
15HERC2CTc.G10969Ap.V3657INM_004667rs1399533760.00253.685TBNNHet
20RALYGAc.G694Ap.G232SNM_016732rs7797450090.00857.476TBNNHet
VF-151PTPRCAGc.A1808Gp.D603GNM_002838rs7546992790.0000496712.77TPDLHet
3TGFBR2GAc.G1283Ap.R428HNM_001024847rs1430957460.000417.14DBDLHet
3GPX1GCC-c.33_35delp.11_12delNM_001329455rs17838762------Het
CADD: Combined Annotation Dependent Depletion, https://cadd.gs.washington.edu/, accessed on 4 February 2019; ExAC: Exome Aggregation Consortium, http://exac.broadinstitute.org/, accessed on 4 February 2019; SIFT: Sorting intolerant from tolerant, https://sift.bii.a-star.edu.sg/, accessed on 4 February 2019; Mutation Taster: Mutation impact testing, http://www.mutationtaster.org/, accessed on 4 February 2019.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ishaq, R.; Ilyas, M.; Habiba, U.; Amin, M.N.u.; Saeed, S.; Raja, G.K.; Shaiq, P.A.; Ahmed, Z.M. Whole Exome Sequencing Reveals Clustering of Variants of Known Vitiligo Genes in Multiplex Consanguineous Pakistani Families. Genes 2023, 14, 1118. https://doi.org/10.3390/genes14051118

AMA Style

Ishaq R, Ilyas M, Habiba U, Amin MNu, Saeed S, Raja GK, Shaiq PA, Ahmed ZM. Whole Exome Sequencing Reveals Clustering of Variants of Known Vitiligo Genes in Multiplex Consanguineous Pakistani Families. Genes. 2023; 14(5):1118. https://doi.org/10.3390/genes14051118

Chicago/Turabian Style

Ishaq, Rafaqat, Muhammad Ilyas, Umme Habiba, Muhammad Noor ul Amin, Sadia Saeed, Ghazala Kaukab Raja, Pakeeza Arzoo Shaiq, and Zubair M. Ahmed. 2023. "Whole Exome Sequencing Reveals Clustering of Variants of Known Vitiligo Genes in Multiplex Consanguineous Pakistani Families" Genes 14, no. 5: 1118. https://doi.org/10.3390/genes14051118

APA Style

Ishaq, R., Ilyas, M., Habiba, U., Amin, M. N. u., Saeed, S., Raja, G. K., Shaiq, P. A., & Ahmed, Z. M. (2023). Whole Exome Sequencing Reveals Clustering of Variants of Known Vitiligo Genes in Multiplex Consanguineous Pakistani Families. Genes, 14(5), 1118. https://doi.org/10.3390/genes14051118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop