- Research article
- Open Access
Genome resequencing and bioinformatic analysis of SNP containing candidate genes in the autoimmune vitiligo Smyth line chicken model
BMC Genomics volume 15, Article number: 707 (2014)
The Smyth line (SL) chicken is the only animal model for autoimmune vitiligo that spontaneously displays all clinical and biological manifestations of the human disorder. To understand the genetic components underlying the susceptibility to develop SL vitiligo (SLV), whole genome resequencing analysis was performed in SLV chickens compared with non-vitiliginous parental Brown line (BL) chickens, which maintain a very low incidence rate of vitiligo.
Illumina sequencing technology and reference based assembly on Red Jungle Fowl genome sequences were used. Results of genome resequencing of pooled DNA of each 10 BL and SL chickens reached 5.1x and 7.0x coverage, respectively. The total number of SNPs was 4.8 and 5.5 million in BL and SL genome, respectively. Through a series of filtering processes, a total of ~1 million unique SNPs were found in the SL alone. Eventually of the 156 reliable marker SNPs, which can induce non-synonymous-, frameshift-, nonsense-, and no-start mutations in amino acid sequences in proteins, 139 genes were chosen for further analysis. Of these, 14 randomly chosen SNPs were examined for SNP verification by PCR and Sanger sequencing to detect SNP positions in 20 BL and 70 SL chickens. The results of the analysis of the 14 SNPs clearly showed differential frequencies of nucleotide bases in the SNP positions between BL and SL chickens. Bioinformatic analysis showed that the 156 most reliable marker SNPs included genes involved in dermatological diseases/conditions such as ADAMTS13, ASPM, ATP6V0A2, BRCA2, COL12A1, GRM5, LRP2, OBSCN, PLAU, RNF168, STAB2, and XIRP1. Intermolecular gene network analysis revealed that candidate genes identified in SLV play a role in networks centered on protein kinases (MAPK, ERK1/2, PKC, PRKDC), phosphatase (PPP1CA), ubiquitinylation (UBC) and amyloid production (APP).
Various potential genetic markers showing amino acid changes and potential roles in vitiligo development were identified in the SLV chicken through genome resequencing. The genetic markers and bioinformatic interpretations of amino acid mutations found in SLV chickens may provide insight into the genetic component responsible for the onset and the progression of autoimmune vitiligo and serve as valuable markers to develop diagnostic tools to detect vitiligo susceptibility.
Vitiligo is an acquired hypomelanotic disorder characterized by circumscribed depigmented macules in the skin resulting from the loss of melanocytes. Autoimmunity has been identified as the major etiological factor in vitiligo, although many other factors including infections, stress, neural abnormalities, aberrant melanocyte function, and genetic susceptibility have been implicated . The Smyth line (SL) chicken is the only animal model for autoimmune vitiligo that spontaneously displays all clinical and biological manifestations of the human disorder [2, 3]. Like other autoimmune diseases, SL vitiligo (SLV) is multi-factorial in nature and involves the interplay of genetic, immune system, and environmental-factors. SLV susceptibility is manifested in part in an inherent melanocyte defect and loss of melanocytes is due to melanocyte-specific cell-mediated and humoral immune activities [2, 3].
Recent genome wide association studies (GWAS) in humans to understand the role of genetic components in a variety of autoimmune diseases including vitiligo have identified hundreds of loci harboring risk alleles . Several GWAS results identified vitiligo susceptible loci in human populations [5–10]. However, most susceptible loci identified by GWAS results were found in regulatory regions of gene expression, therefore the identified associations were not sufficient to identify the causal gene or deduce alterations caused by risk variants, which generally do not induce profound changes to genes (e.g. coding sequence changes, deletions, or duplications). Recently, the encyclopedia of DNA elements (ENCODE) of mammalian species suggested that ~90% of disease associated genetic variations in human lie in noncoding regions, while only ~10% of variations in coding regions were causative mutations associated with human disease [11, 12]. Nevertheless, the identification of potential coding mutations that alter protein functionalities is a prerequisite process to understand disease etiology. Moreover, the functional study of candidate genetic risk factors is almost impossible without appropriate model systems.
The SL chicken is an excellent model to conduct a functional verification study of candidate genes that underlie genetic susceptibility for vitiligo due to the tractable, definite phenotype, the high vitiligo incidence in the population (80-90%), the feasibility of in vivo characterization and the relatively short generation time. Recently, microarray analysis examined global gene expression during SLV development and provided comprehensive information at the transcriptome level that supported the multifactorial etiology of vitiligo . In this study, whole genome resequencing analysis using an Illumina platform was performed to more deeply investigate the genetic aspects of SLV expression in comparison with the parental Brown line (BL) of chickens from which the SL was originally derived. BL chickens retain vitiligo susceptibility but with a very low (0 – 2%) incidence rate of vitiligo development [2, 3], although none of the BL chickens used in this study had vitiligo. Millions of single nucleotide polymorphisms (SNP) were identified by genome resequencing and only potentially causal genes containing non-synonymous mutations that can induce amino acid changes in proteins were focused on in this study.
Results and discussion
Genome resequencing for BL and SL chickens
Genome sequencing of pooled DNA from 10 non-vitiliginous BL and SL chickens each with confirmed SLV produced ~63 and 89 million sequence reads of 200 bp, respectively (Table 1). Of those, ~80% of the reads were used for sequence alignment, while 20% of sequence reads were not aligned. Therefore, genome coverage for BL and SL reached 5.1x and 7.0x, respectively, of the Red Jungle Fowl chicken genome. The total number of SNPs was 4.8 and 5.5 million (~0.5% of template genome) for the BL and SL genome, respectively. The large number of SNPs per examined chicken line was based on data of at least 2 read coverage depths (number of read counts per nucleotide location). Most SNPs were found on the larger chromosomes (Chr), including Chr 1 through 5 and Z (sex chromosome) (data not shown). To identify genetic biomarkers that are responsible for the incidence of SLV, unique SNPs that are found in SL only were selected by removing SNPs that overlapped with those found in the parental BL. Then, mutations with ≥75% SNP rates were chosen as reliable marker SNPs. Since the objective of this study was to identify mutational SNPs uniquely found in SL compared to the parental BL, the filtering process used did not involve a typical SNP calling and filtering method based on quality score (Q call column in Additional file 1: Table S1) . Instead, SNP filtering was conducted by removing overlapping SNPs found in both BL and SL, and applying fixed %SNP rates (≥75%) as described in Methods. As a result, a total of ~1 million unique SNPs were identified throughout the SL chicken genome (Figure 1). Over 100,000 SNPs were found in Chr 1, 2, 3, and Z while Chr 32 did not contain any unique SNPs for SL. When SNPs were grouped by the feature types of chromosome regions, ~50% of SNPs were in the intergenic (heterochromatic) regions and 13,710 SNPs were found in CDS sequences (protein coding regions) (Figure 2). Most genes containing SNPs in regulatory regions, not in CDS regions, identified by human GWAS studies were also observed to contain unique SNPs in the current SL study (data not shown). Around 60% of SL SNPs in protein coding CDS regions were synonymous mutations that did not induce amino acid changes. To identify potentially causal mutations that induce protein coding alterations, SNP analysis focused on SNPs leading to changes in amino acid sequences. Using this approach, a total of 3518 SNPs were identified that could induce non-synonymous-, frameshift-, nonsense-, and no-start-changes in the CDS region (Figure 2 and Additional file 1: Table S1), suggesting that the 3518 SNPs are part of the genetic components in functional protein coding regions that may drive the high incidence rate of SLV. Of the 3518 candidate SNPs that are associated with amino acid changes, SNPs showing ≥10 read depths (considered to be more reliable candidate genetic markers) were chosen for the further analysis. Using this approach, 195 SNPs remained (data not shown). To reduce false positives due to possible errors in the assembly process, re-scanning of each SNP position for the 195 potentially more reliable protein coding SNPs was conducted using the Seqman-Pro viewer program. This process yielded 156 more reliable SNPs that were chosen as candidate marker SNPs for further analysis (Table 2).
SNP validation using PCR and Sanger sequencing
Since pooled DNA samples of 10 chickens for each line were used for genome sequencing, individual SNPs were subjected to the verification process with larger bird populations. For this, 14 SNPs were randomly chosen from the 156 candidate marker SNPs showing ≥10 read depths and were subjected to SNP verification analysis using PCR and Sanger sequencing to detect SNP positions with larger numbers of birds; specifically 20 non-vitiliginous BL chickens and 70 SL chickens exhibiting vitiligo. The results clearly showed differential frequencies of nucleotide bases in the 14 SNP positions between BL and SL chickens (Table 3). Thus, the 156 SNPs known to induce amino acid changes can become potential genetic biomarkers for vitiligo in SL chickens.
Bioinformatic analyses of genes containing amino acid change SNPs
Amino acid changes may have impacts on the functional interpretations for vitiligo induction in SL chickens. The Ingenuity Pathway Analysis (IPA) program generated bioinformatics data sets including functional groups (gene ontology; GO) and gene networks for genes containing amino acid changes in SL chicken. The 156 SNPs were found in 139 genes encompassing known- and unknown functions, chromosomal open reading frames, and hypothetical proteins (Additional file 2: Table S2).
Genes were categorized in 76 functional groups (Additional file 3: Table S3). Of these, six functional groups are of particular interest to autoimmune vitiligo development, including dermatological diseases/conditions, inflammatory response, inflammatory disease, immunological disease, immune cell trafficking, and infectious disease (Table 4). The functional group of genes for dermatological diseases/conditions contained the following genes: ADAMTS13 (ADAM metallopeptidase with thrombospondin type 1 motif 13); ASPM [asp (abnormal spindle) homolog, microcephaly associated; Drosophila)]; ATP6V0A2 (ATPase, H + transporting, lysosomal V0 subunit A2); BRCA2 (breast cancer 2, early onset); COL12A1 (collagen, type XII, alpha 1); GRM5 (glutamate receptor, metabotropic 5); LRP2 (low density lipoprotein receptor-related protein 2); MKI67 (marker of proliferation Ki-67); OBSCN (obscurin, cytoskeletal calmodulin and titin-interacting RhoGEF); PLAU (plasminogen activator, urokinase); RNF168 (ring finger protein 168, E3 ubiquitin protein ligase); STAB2 [stabilin 2 or FEEL2 (fasciclin EGF-like, laminin-type EGF-like, and link domain-containing scavenger receptor 2)]; and XIRP1 (xin actin-binding repeat containing 1). General and dermatological disease related functions for these genes are summarized in Table 5.
Interestingly, a recent report by Nikolaev et al. (2012)  indicated that amino acid changes found in ASPM, LRP2, STAB2, and XIRP1 proteins were associated with human melanoma by exome sequencing. Melanocytes in vitiligo also exhibit morphological and biological melanocyte defects/alterations compared to melanocytes from individuals with normal pigmentation . While these alterations may be different from those observed in melanoma, e.g. slower growth, and higher sensitivity to oxidative stress of cultured melanocytes [30, 31], alterations in amino acid sequences found in homolog proteins but different residues may result in opposite phenotypes of dermatological diseases/conditions . In addition to these molecules, BRCA2, GRM5, MKI67, and OBSCN associated with SLV are also known to be associated with human melanoma. The relationship between candidate genes and other dermatological diseases including melanoma is summarized in Table 5.
Gene network analysis, which represents the intermolecular connections among interacting genes based on functional knowledge inputs, was performed on genes with amino acid changes in SLV chickens using the IPA program. The gene network analysis was carried out using the simplest setting of 35 focus molecules to facilitate and summarize the intermolecular connections (Table 6 and Figures 3, 4, 5, 6, 7, 8, and 9). A discussion of the 7 gene networks is provided below and gene information for focus molecules in each network is listed in Additional file 4: Table S4.
Candidate genes in Network #1 are associated with signaling pathways of the mitogen activated protein kinase (MAPK; also ERK1/2) and protein kinase C (Pkc) connected to VEGF (vascular endothelial growth factor) and PDGF (platelet derived growth factor) with PLAU in the center (Figure 3). The top functions related to network #1 are cardiovascular disease, hematological disease, and cardiac infarction. Interestingly, molecules including LRP2, PLAU, ADAMTS13, and GRM5 that are part of Network #1 were also identified as functional factors for melanoma as described above . Additionally, mutations in the amino acid sequence of LRP2, altered function of MAP2K1 and MAP2K2 induced by genetic mutations in melanoma patients , and mutations in GRM5 in mouse melanoma models  were also reported. The connections in Network #1 therefore suggest genetic mutations that generated amino acid changes in LRP2, PLAU, ADAMTS13, and GRM5 may influence dermatological diseases, including vitiligo, through MAPK and ERK1/2 signaling pathways.
The top functions of Network #2 include Cell Cycle, DNA Replication, Recombination and Repair, and Developmental Disorder (Figure 4) and Network #3 is involved in Developmental Disorder, Endocrine System Disorders, Hereditary Disorder (Figure 5). Most molecules in Networks #2 and #3 directly bind to UBC (ubiquitin C). Ubiquitinylation, the covalent attachment of ubiquitin to proteins, regulates numerous cellular processes such as protein degradation and signal transduction. Recently, many ubiquitinylated proteins and their lysine ubiquitinylation site were identified using proteomic technologies in mammalian species [33–37]. Indeed, in SLV, the amino acid changes in lysine residues for CEP192 (centrosomal protein 192 kDa; pK169R) and API5 (apoptosis inhibitor 5; pK1219Q) were identified (Table 2), suggesting that the various cellular functions including protein degradation by altered ubiquitinylation properties of proteins may play a significant role in the induction of vitiligo.
Molecules in Network #4 are involved in Cell Morphology, Cellular Function and Maintenance, Hair and Skin Development and Function (Figure 6). Molecules in network #4 mainly interact with IFNG (interferon gamma), IL4 (interleukin 4), MAPK8, and calcium signaling pathways in addition to protein phosphatase (PPP1CA; protein phosphatase 1 catalytic subunit alpha isozyme) functions. MYO16 (myocin 16), KIF18A (kinesin family member 18A), and CASC1 (cancer susceptibility candidate 1) are known to directly bind to PPP1CA [38–40], suggesting that amino acid changes in those proteins may induce alterations (hyper vs hypo) in protein phophorylation states in SL chickens, resulting in vitiligo development. Molecules interacting with IFNG, IL4, MAPK8 and calcium signaling pathways showed an indirect relationship, not a direct relationship with each other making it difficult to explain how amino acid changes in these molecules affect vitiligo induction in SL chicken.Molecules in Network #5 also mainly bind to UBC as discussed in Networks #2 and #3 and the functions include Lipid Metabolism, Small Molecule Biochemistry, Digestive System Development and Function (Figure 7).
Network #6 contains molecules involved in Hereditary Disorder, Ophthalmic Disease, Neurological Disease (Figure 8). In this network, ZC2HC1A (zinc finger, C2HC-type containing 1A) and RUFY3 (RUN and FYVE domain containing 3) directly bind to APP [amyloide beta (A4) precursor protein]. APP is a precursor protein for beta-amyloide, which is the main constituent of amyloid plaques in the brains of Alzheimer disease patients . RUFY3 (also known as single axon-related; singar1), which is a brain specific protein, regulates neuronal polarity by suppressing formation of surplus axons . Though the binding of ZC2HC1A and RUFY3 to APP was found during the progression of Alzheimer disease , the functional roles for this binding in the progression of this disease has not been characterized. Similarly, amino acid changes found in ZC2HC1A and RUFY3 are implicated in SLV development possibly as a result of altered APP binding properties. USH2A [Usher syndrome 2A (autosomal recessive, mild)] was included in network #6. Various mutations in USH2A have been identified in patients of Usher syndrome type II, which is characterized by moderate to severe sensorineural hearing loss and progressive retinitis pigmentosa . Vitiliginous SL chickens may also develop severe visual impairment and blindness due to autoimmune activity directed against choroidal melanocytes and subsequent damage to the retinal pigment epithelium. [44, 45]. Taken together, the amino acid change in USH2A also may affect vitiligo progression and retinal depigmentation.
Molecules in Network #7 mainly function in Cellular Response to Therapeutics, Cellular Assembly and Organization, DNA Replication, Recombination, and Repair. PRKDC (protein kinase, DNA-activated, catalytic polypeptide) and its interacting protein kinases are mainly involved in this network (Figure 9). In addition to knock-out and inactive mutations, alteration of autophosphorylation capability by single amino acid change of PRKDC has been known to influence rejoining of DNA double stranded breaks in mammalian cells , suggesting that vitiligo development may be affected by aberrant PRKDC kinase activity due to an observed amino acid change.
Vitiligo susceptible loci in human populations identified by several GWAS [5–10] also showed several loci that induced amino acid changes in various proteins such as STRN3 (Striatin, calmodulin binding protein 3), DNAH5 (dynein, axonemal, heavy chain 5), KIAA1005 (immunoglobulin heavy variable 3), TYR (tyrosinase), OCA2 (oculocutaneous albinism II), PTPN22 [protein tyrosine phosphatase, non-receptor type 22 (lymphoid)], IFIH1 (interferon induced with helicase C domain 1), SLA (SRC-like-adaptor), CD44, MC1R [melanocortin 1 receptor (alpha melanocyte stimulating hormone receptor)], UBASH3A (ubiquitin associated and SH3 domain containing A), C1QTNF6 (C1Q and tumor necrosis factor related protein 6), CASP7 (caspase 7), and GZMB (granzyme B). One similar mutation in UBASH3A protein coding region from human  was also found in the SLV chicken model (Table 2). Additionally, when the long list of 3518 candidate amino acid altering SNPs (read depth <10) was considered, several genes, including IFIH1 (interferon-induced helicase C domain-containing protein 1), CD44 antigen and DNAH5 (dynein, axonemal, heavy chain 5), matched those identified by human GWAS studies [Additional file 1: Table S1 and [5, 6]], although the amino acid position and alterations did not match. UBASH3A is one of the two family members belonging to the T cell ubiquitin ligand (TULA) family and can negatively regulate T cell signaling . Together with the UBC molecule discussed elsewhere in this paper, functions for UBASH3A related to ubiquitinylation and T cell signaling pathway may be important for SLV development. IFIH1 encodes an interferon-induced RNA helicase involved in antiviral innate immune responses, associated with type 1 diabetes, Graves’ disease, multiple sclerosis, psoriasis, and perhaps lupus [48–53]. CD44 encodes a cell surface glycoprotein with various functions, including a role in T cell development , and is associated with lupus . DNAH5 gene mutation is found in patients with primary ciliary dyskinesia (PCD) , a rare disease transmitted as an autosomal recessive trait and characterized by recurrent airway infections due to abnormal ciliary structure and function. Primary defects in the structure and function of sensory and motile cilia result in multiple ciliopathies .
In this study, various potential genetic markers showing amino acid changes were identified in the SLV model through genome re-sequencing. When considering functionality based on the interpretation of factors involved, development of vitiligo appeared to be associated with the interactions among cytoskeletal factors (OBSCN, ASPM, XIRP1, ADAMTS13), protein kinases (MAPK, ERK1/2, PKC, PRKDC), phosphatase (PPP1CA), ubiquitinylation (UBC) and amyloid (APP) production. Further functional validation study, such as allele specific expression of the candidate genes with candidate SNPs at the target tissues involved in SLV development will be carried out using the SL chicken model for spontaneous autoimmune vitiligo.
Animals and Illumina sequencing
Adult SL chickens with vitiligo and parental non-vitiliginous BL chickens, maintained by G. Erf at the University of Arkansas (Fayetteville, AR), were selected from the breeder populations. Blood (3 ml) was collected from 12 birds each following an animal use protocol approved by the University of Arkansas Institutional Animal Care and Use Committee (IACUC; approval number: 11019). Genomic DNA was isolated from each whole blood sample using QiaAmp DNA mini kit (Qiagen, Hilden, Germany) following manufacturer’s instructions. DNA quality was determined by agarose gel electrophoresis and 10 samples having the highest quality in each line were pooled to represent each line. Library preparation and Illumina genome sequencing for the pooled DNA samples were performed by the National Center for Genome Resources (NCGR; Santa Fe, NM). Illumina HiSeq system 2x100 bp paired end read technology was used for genome sequencing.
Genome sequence assembly and data analysis
Illumina sequencing data received from NCGR was aligned to the chicken reference genome sequence for Red Jungle Fowl (GBK 4.0) that was retrieved from NCBI. For the reference based genome alignment, the NGen genome sequence assembly program of the Lasergene software package (DNAStar, Madison, WI) was used. Assembly parameters were as follows: file format, BAM; mer Size, 21; mer skip query, 2; minimum match percentage, 93; maximum gap size, 6; minimum aligned length, 35; match score, 10; mismatch penalty, 20; gap penalty, 30; SNP calculation method, diploid bayesian; minimum SNP percentage, 5; SNP confidence threshold, 10; minimum SNP count, 2; minimum base quality score, 5. After assembly, the SeqMan Pro program of the Lasergene package was used for further analyses including SNP data.
SNP detection and analysis
JMP genomics (SAS Institute, Inc., Cary, NC) program was used for filtering unique SNPs for vitiligo SL chickens. SNPs occurring in both SL and BL lines were filtered out, leaving behind unique SNPs for each line. To identify highly fixed and homozygous SNPs, the SNPs were filtered based on SNP percentages (SNP%). SNPs with a SNP% of ≥75 (%) (for example, number of SNP = 3 of read depth = 4) were chosen. The 75% cutoff for SNP selection was set by considering potential sequencing errors that can be generated by the massively parallel sequencing method. Potential causal SL SNPs that induce non-synonymous changes in CDS regions were chosen for further analysis. Since the read depth of many SL SNPs was low, unique SNPs showing ≥10 read depths were considered as reliable SNPs. Reliable and causal SNPs, which were chosen by criteria described above were confirmed by double-checking the raw assembly data with alignment view to reduce false positives.
SNP validation using PCR and Sanger sequencing
Fourteen randomly chosen SNPs, which induce amino acid changes in the CDS region, were subjected to validation using PCR and Sanger sequencing with larger numbers of SL and BL chickens. Twenty BL and 70 SL chickens that were verified phenotypically to be non-vitiliginous and vitiliginous, respectively, were used for blood sampling. Approximately 100 μL of blood was collected from each bird by wing vein puncture into tubes containing citrate (anticoagulant). Genomic DNA was isolated from whole blood using the Wizard SV 96 Genomic DNA Purification System (Promega; Madison, WI) following manufacturer’s instructions with modifications. Isolated DNA was quantified using a Nanodrop 1000 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA) and a dilution of 1 ng/μL was prepared in 96 well PCR format for all samples. For PCR reaction, forward and reverse primers were designed based on the RJF genome sequence (GenBank assembly ID: GCA_000002315.2) using Primer 3 online software (Table 7). The sequencing primers were designed to anneal at least 50 bp upstream of the SNP position and forward/reverse primers were chosen at the flanking regions of the seq primer and the SNP position. All primers were commercially synthesized by Integrated DNA Technology (Ames, IA). PCR was carried out as 25 μL reaction volumes in 96 well plates with cyle conditions as follows: denaturation 95°C for 1 min, 40 cycles of amplification (95°C for 30 s, 55-63°C for 1 min, 72°C for 1 min), final extension 72°C for 10 min. Verification of PCR was performed by 1% agarose gel electrophoresis. PCR products were purified using the Wizard SV 96 PCR Clean-Up System (Promega; Madison, WI) following manufacturer’s instructions. Briefly, four plates (four different PCR products) were pooled into one plate and were subjected to PCR clean-up. Cross-specificity of seq-primers to the four pooled PCR products was examined by BLAST function (NCBI) and only products that were not cross-specific with other seq primers were pooled. Purified PCR products were subjected to Sanger sequencing performed by the University of Arkansas DNA Resources Center (Fayetteville, AR). Results were analyzed using ABI Sequence scanner software (Life Technologies, Carlsbad, CA). Ratios of bases occurring at SNP locations were recorded.
Functional interpretation of 139 genes showing ≥75 SNP%, ≥10 read depths and non-synonymous changes was analyzed in the context of gene ontology and molecular networks using Ingenuity Pathways Analysis (IPA; Ingenuity Systems®; http://www.ingenuity.com). Since IPA is based on human and mouse bioinformatics, functionalities for differentially expressed genes in the chicken were interpreted based primarily on mammalian biological mechanisms. The limit of number of molecules in the network was set to 35, leaving only the most important molecules based on the number of connections for each focus gene (a subset of uploaded significant genes having direct interactions with other genes in the database) to other significant genes .
Availability of supporting data
All sequence reads described in the manuscript are available under BioProject accession PRJNA256208. Illumina sequence reads have been deposited at NCBI’s SRA archive under following numbers (SL: Sample: SRS670088, Experiment: SRX665272, Read: SRR1531502; BL: Sample: SRS670098, Experiment: SRX665286, Reads: SRR1531503).
Spritz RA: Six decades of vitiligo genetics: genome-wide studies provide insights into autoimmune pathogenesis. J Invest Dermatol. 2012, 132 (2): 268-273.
Erf GF: Ainmal Models. Vitiligo. Edited by: Picardo M, Taieb A. 2010, Berlin Heidelberg, Germany: Springer-Verlag GmbH, 205-218.
Erf GF: Autoimmune Diseases of Poultry. Avian Immunology. Edited by: Davison F, Kaspers B, Schat K. 2008, London: Elsevier
Cotsapas C, Hafler DA: Immune-mediated disease genetics: the shared basis of pathogenesis. Trends Immunol. 2013, 34 (1): 22-26.
Cheong KA, Kim NH, Noh M, Lee AY: Three new single nucleotide polymorphisms identified by a genome-wide association study in Korean patients with vitiligo. J Korean Med Sci. 2013, 28 (5): 775-779.
Jin Y, Birlea SA, Fain PR, Ferrara TM, Ben S, Riccardi SL, Cole JB, Gowan K, Holland PJ, Bennett DC, Luiten RM, Wolkerstorfer A, van der Veen JP, Hartmann A, Eichner S, Schuler G, van Geel N, Lambert J, Kemp EH, Gawkrodger DJ, Weetman AP, Taieb A, Jouary T, Ezzedine K, Wallace MR, McCormack WT, Picardo M, Leone G, Overbeck A, Silverberg NB, et al: Genome-wide association analyses identify 13 new susceptibility loci for generalized vitiligo. Nat Genet. 2012, 44 (6): 676-680.
Jin Y, Birlea SA, Fain PR, Gowan K, Riccardi SL, Holland PJ, Bennett DC, Herbstman DM, Wallace MR, McCormack WT, Kemp EH, Gawkrodger DJ, Weetman AP, Picardo M, Leone G, Taieb A, Jouary T, Ezzedine K, van Geel N, Lambert J, Overbeck A, Spritz RA: Genome-wide analysis identifies a quantitative trait locus in the MHC class II region associated with generalized vitiligo age of onset. J Invest Dermatol. 2011, 131 (6): 1308-1312.
Jin Y, Ferrara T, Gowan K, Holcomb C, Rastrou M, Erlich HA, Fain PR, Spritz RA: Next-generation DNA re-sequencing identifies common variants of TYR and HLA-A that modulate the risk of generalized vitiligo via antigen presentation. J Invest Dermatol. 2012, 132 (6): 1730-1733.
Tang J, Liu JL, Zhang C, da Hu Y, He SM, Zuo XB, Wang PG, Sun LD, Zhang XJ, Yang S: The association between a single nucleotide polymorphism rs11966200 in MHC region and clinical features of generalized vitiligo in Chinese Han population. Mol Biol Rep. 2013, 40 (6): 4097-4100.
Tang XF, Zhang Z, Hu DY, Xu AE, Zhou HS, Sun LD, Gao M, Gao TW, Gao XH, Chen HD, Xie HF, Tu CX, Hao F, Wu RN, Zhang FR, Liang L, Pu XM, Zhang JZ, Han JW, Pan GP, Wu JQ, Li K, Su MW, Du WD, Zhang WJ, Liu JJ, Xiang LH, Yang S, Zhou YW, Zhang XJ: Association analyses identify three susceptibility loci for vitiligo in the Chinese Han population. J Invest Dermatol. 2013, 133 (2): 403-410.
Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009, 106 (23): 9362-9367.
Maurano MT, Humbert R, Rynes E, Thurman RE, Haugen E, Wang H, Reynolds AP, Sandstrom R, Qu H, Brody J, Shafer A, Neri F, Lee K, Kutyavin T, Stehling-Sun S, Johnson AK, Canfield TK, Giste E, Diegel M, Bates D, Hansen RS, Neph S, Sabo PJ, Heimfeld S, Raubitschek A, Ziegler S, Cotsapas C, Sotoodehnia N, Glass I, Sunyaev SR, et al: Systematic localization of common disease-associated variation in regulatory DNA. Science. 2012, 337 (6099): 1190-1195.
Shi F, Kong BW, Song JJ, Lee JY, Dienglewicz RL, Erf GF: Understanding mechanisms of vitiligo development in Smyth line of chickens by transcriptomic microarray analysis of evolving autoimmune lesions. BMC Immunol. 2012, 13: 18-2172-13-18-
Nielsen R, Paul JS, Albrechtsen A, Song YS: Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011, 12 (6): 443-451.
Balakrishnan A, Bleeker FE, Lamba S, Rodolfo M, Daniotti M, Scarpa A, van Tilborg AA, Leenstra S, Zanon C, Bardelli A: Novel somatic and germline mutations in cancer candidate genes in glioblastoma, melanoma, and pancreatic carcinoma. Cancer Res. 2007, 67 (8): 3545-3550.
Tamura Y, Adachi H, Osuga J, Ohashi K, Yahagi N, Sekiya M, Okazaki H, Tomita S, Iizuka Y, Shimano H, Nagai R, Kimura S, Tsujimoto M, Ishibashi S: FEEL-1 and FEEL-2 are endocytic receptors for advanced glycation end products. J Biol Chem. 2003, 278 (15): 12613-12617.
Nikolaev SI, Rimoldi D, Iseli C, Valsesia A, Robyr D, Gehrig C, Harshman K, Guipponi M, Bukach O, Zoete V, Michielin O, Muehlethaler K, Speiser D, Beckmann JS, Xenarios I, Halazonetis TD, Jongeneel CV, Stevenson BJ, Antonarakis SE: Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma. Nat Genet. 2011, 44 (2): 133-139.
Spoelgen R, Hammes A, Anzenberger U, Zechner D, Andersen OM, Jerchow B, Willnow TE: LRP2/megalin is required for patterning of the ventral telencephalon. Development. 2005, 132 (2): 405-414.
Kouprina N, Pavlicek A, Collins NK, Nakano M, Noskov VN, Ohzeki J, Mochida GH, Risinger JI, Goldsmith P, Gunsior M, Solomon G, Gersch W, Kim JH, Barrett JC, Walsh CA, Jurka J, Masumoto H, Larionov V: The microcephaly ASPM gene is expressed in proliferating tissues and encodes for a mitotic spindle protein. Hum Mol Genet. 2005, 14 (15): 2155-2165.
Choi KY, Chang K, Pickel JM, Badger JD 2nd, Roche KW: Expression of the metabotropic glutamate receptor 5 (mGluR5) induces melanoma in transgenic mice. Proc Natl Acad Sci U S A. 2011, 108 (37): 15219-15224.
Jensen RB, Carreira A, Kowalczykowski SC: Purified human BRCA2 stimulates RAD51-mediated recombination. Nature. 2010, 467 (7316): 678-683.
The Breast Cancer Linkage Consortium: Cancer risks in BRCA2 mutation carriers. The Breast Cancer Linkage Consortium. J Natl Cancer Inst. 1999, 91 (15): 1310-1316.
Otten C, van der Ven PF, Lewrenz I, Paul S, Steinhagen A, Busch-Nentwich E, Eichhorst J, Wiesner B, Stemple D, Strahle U, Furst DO, Abdelilah-Seyfried S: Xirp proteins mark injured skeletal muscle in zebrafish. PLoS One. 2012, 7 (2): e31041-
Stewart GS, Panier S, Townsend K, Al-Hakim AK, Kolas NK, Miller ES, Nakada S, Ylanko J, Olivarius S, Mendez M, Oldreive C, Wildenhain J, Tagliaferro A, Pelletier L, Taubenheim N, Durandy A, Byrd PJ, Stankovic T, Taylor AM, Durocher D: The RIDDLE syndrome protein mediates a ubiquitin-dependent signaling cascade at sites of DNA damage. Cell. 2009, 136 (3): 420-434.
Pimanda JE, Maekawa A, Wind T, Paxton J, Chesterman CN, Hogg PJ: Congenital thrombotic thrombocytopenic purpura in association with a mutation in the second CUB domain of ADAMTS13. Blood. 2004, 103 (2): 627-629.
Kornak U, Reynders E, Dimopoulou A, van Reeuwijk J, Fischer B, Rajab A, Budde B, Nurnberg P, Foulquier F, Lefeber D, Urban Z, Gruenewald S, Annaert W, Brunner HG, van Bokhoven H, Wevers R, Morava E, Matthijs G, Van Maldergem L, Mundlos S, ARCL Debre-type Study Group: Impaired glycosylation and cutis laxa caused by mutations in the vesicular H+-ATPase subunit ATP6V0A2. Nat Genet. 2008, 40 (1): 32-34.
Kiian I, Tkachuk N, Haller H, Dumler I: Urokinase-induced migration of human vascular smooth muscle cells requires coupling of the small GTPases RhoA and Rac1 to the Tyk2/PI3-K signalling pathway. Thromb Haemost. 2003, 89 (5): 904-914.
Bolon I, Zhou HM, Charron Y, Wohlwend A, Vassalli JD: Plasminogen mediates the pathological effects of urokinase-type plasminogen activator overexpression. Am J Pathol. 2004, 164 (6): 2299-2304.
Le Poole IC, Luiten RM: Autoimmune etiology of generalized vitiligo. Current directions in autoimmunity: dermatologic immunity. Volume 22. Edited by: Nickloff BJ, Nestle FO. 2008, Basel: Karger, 566-569.
Medrano EE, Nordlund JJ: Successful culture of adult human melanocytes obtained from normal and vitiligo donors. J Invest Dermatol. 1990, 95 (4): 441-445.
Bowers RR, Nguyen B, Buckner S, Gonzalez Y, Ruiz F: Role of antioxidants in the survival of normal and vitiliginous avian melanocytes. Cell Mol Biol (Noisy-le-grand). 1999, 45 (7): 1065-1074.
Spritz RA: The genetics of vitiligo. J Invest Dermatol. 2011, 131 (E1): E18-E20.
Wagner SA, Beli P, Weinert BT, Nielsen ML, Cox J, Mann M, Choudhary C: A proteome-wide, quantitative survey of in vivo ubiquitylation sites reveals widespread regulatory roles. Mol Cell Proteomics. 2011, 10 (10): M111.013284-
Wagner SA, Beli P, Weinert BT, Scholz C, Kelstrup CD, Young C, Nielsen ML, Olsen JV, Brakebusch C, Choudhary C: Proteomic analyses reveal divergent ubiquitylation site patterns in murine tissues. Mol Cell Proteomics. 2012, 11 (12): 1578-1585.
Danielsen JM, Sylvestersen KB, Bekker-Jensen S, Szklarczyk D, Poulsen JW, Horn H, Jensen LJ, Mailand N, Nielsen ML: Mass spectrometric analysis of lysine ubiquitylation reveals promiscuity at site level. Mol Cell Proteomics. 2011, 10 (3): M110.003590-
Udeshi ND, Mani DR, Eisenhaure T, Mertins P, Jaffe JD, Clauser KR, Hacohen N, Carr SA: Methods for quantification of in vivo changes in protein ubiquitination following proteasome and deubiquitinase inhibition. Mol Cell Proteomics. 2012, 11 (5): 148-159.
Lee KA, Hammerle LP, Andrews PS, Stokes MP, Mustelin T, Silva JC, Black RA, Doedens JR: Ubiquitin ligase substrate identification through quantitative proteomics at both the protein and peptide levels. J Biol Chem. 2011, 286 (48): 41530-41538.
Colland F, Jacq X, Trouplin V, Mougin C, Groizeleau C, Hamburger A, Meil A, Wojcik J, Legrain P, Gauthier JM: Functional proteomics mapping of a human signaling pathway. Genome Res. 2004, 14 (7): 1324-1332.
Patel KG, Liu C, Cameron PL, Cameron RS: Myr 8, a novel unconventional myosin expressed during brain development associates with the protein phosphatase catalytic subunits 1alpha and 1gamma1. J Neurosci. 2001, 21 (20): 7954-7968.
Hendrickx A, Beullens M, Ceulemans H, Den Abt T, Van Eynde A, Nicolaescu E, Lesage B, Bollen M: Docking motif-guided mapping of the interactome of protein phosphatase-1. Chem Biol. 2009, 16 (4): 365-371.
Olah J, Vincze O, Virok D, Simon D, Bozso Z, Tokesi N, Horvath I, Hlavanda E, Kovacs J, Magyar A, Szucs M, Orosz F, Penke B, Ovadi J: Interactions of pathological hallmark proteins: tubulin polymerization promoting protein/p25, beta-amyloid, and alpha-synuclein. J Biol Chem. 2011, 286 (39): 34088-34100.
Mori T, Wada T, Suzuki T, Kubota Y, Inagaki N: Singar1, a novel RUN domain-containing protein, suppresses formation of surplus axons for neuronal polarity. J Biol Chem. 2007, 282 (27): 19884-19893.
Eudy JD, Weston MD, Yao S, Hoover DM, Rehm HL, Ma-Edmonds M, Yan D, Ahmad I, Cheng JJ, Ayuso C, Cremers C, Davenport S, Moller C, Talmadge CB, Beisel KW, Tamayo M, Morton CC, Swaroop A, Kimberling WJ, Sumegi J: Mutation of a gene encoding a protein with extracellular matrix motifs in Usher syndrome type IIa. Science. 1998, 280 (5370): 1753-1757.
Boissy RE: The melanocyte. Its structure, function, and subpopulations in skin, eyes, and hair. Dermatol Clin. 1988, 6 (2): 161-173.
Boissy RE, Moellmann G, Trainer AT, Smyth JR, Lerner AB: Delayed-amelanotic (DAM or Smyth) chicken: melanocyte dysfunction in vivo and in vitro. J Invest Dermatol. 1986, 86 (2): 149-156.
Chan DW, Chen BP, Prithivirajsingh S, Kurimasa A, Story MD, Qin J, Chen DJ: Autophosphorylation of the DNA-dependent protein kinase catalytic subunit is required for rejoining of DNA double-strand breaks. Genes Dev. 2002, 16 (18): 2333-2338.
Collingwood TS, Smirnova EV, Bogush M, Carpino N, Annan RS, Tsygankov AY: T-cell ubiquitin ligand affects cell death through a functional interaction with apoptosis-inducing factor, a key factor of caspase-independent apoptosis. J Biol Chem. 2007, 282 (42): 30920-30928.
Kato H, Takeuchi O, Sato S, Yoneyama M, Yamamoto M, Matsui K, Uematsu S, Jung A, Kawai T, Ishii KJ, Yamaguchi O, Otsu K, Tsujimura T, Koh CS, Reis e S C, Matsuura Y, Fujita T, Akira S: Differential roles of MDA5 and RIG-I helicases in the recognition of RNA viruses. Nature. 2006, 441 (7089): 101-105.
Smyth DJ, Cooper JD, Bailey R, Field S, Burren O, Smink LJ, Guja C, Ionescu-Tirgoviste C, Widmer B, Dunger DB, Savage DA, Walker NM, Clayton DG, Todd JA: A genome-wide association study of nonsynonymous SNPs identifies a type 1 diabetes locus in the interferon-induced helicase (IFIH1) region. Nat Genet. 2006, 38 (6): 617-619.
Sutherland A, Davies J, Owen CJ, Vaikkakara S, Walker C, Cheetham TD, James RA, Perros P, Donaldson PT, Cordell HJ, Quinton R, Pearce SH: Genomic polymorphism at the interferon-induced helicase (IFIH1) locus contributes to Graves’ disease susceptibility. J Clin Endocrinol Metab. 2007, 92 (8): 3338-3341.
Martinez A, Santiago JL, Cenit MC, de Las HV, de la Calle H, Fernandez-Arquero M, Arroyo R, de la Concha EG, Urcelay E: IFIH1-GCA-KCNH7 locus: influence on multiple sclerosis risk. Eur J Hum Genet. 2008, 16 (7): 861-864.
Li Y, Liao W, Cargill M, Chang M, Matsunami N, Feng BJ, Poon A, Callis-Duffin KP, Catanese JJ, Bowcock AM, Leppert MF, Kwok PY, Krueger GG, Begovich AB: Carriers of rare missense variants in IFIH1 are protected from psoriasis. J Invest Dermatol. 2010, 130 (12): 2768-2772.
Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, Ortmann W, Kosoy R, Ferreira RC, Nordmark G, Gunnarsson I, Svenungsson E, Padyukov L, Sturfelt G, Jonsen A, Bengtsson AA, Rantapaa-Dahlqvist S, Baechler EC, Brown EE, Alarcon GS, Edberg JC, Ramsey-Goldman R, McGwin G, Reveille JD, Vila LM, Kimberly RP, Manzi S, Petri MA, Lee A, Gregersen PK, et al: A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet. 2009, 41 (11): 1228-1233.
Baaten BJ, Li CR, Bradley LM: Multifaceted regulation of T cells by CD44. Commun Integr Biol. 2010, 3 (6): 508-512.
Ramos PS, Williams AH, Ziegler JT, Comeau ME, Guy RT, Lessard CJ, Li H, Edberg JC, Zidovetzki R, Criswell LA, Gaffney PM, Graham DC, Graham RR, Kelly JA, Kaufman KM, Brown EE, Alarcon GS, Petri MA, Reveille JD, McGwin G, Vila LM, Ramsey-Goldman R, Jacob CO, Vyse TJ, Tsao BP, Harley JB, Kimberly RP, Alarcon-Riquelme ME, Langefeld CD, Moser KL: Genetic analyses of interferon pathway-related genes reveal multiple new loci associated with systemic lupus erythematosus. Arthritis Rheum. 2011, 63 (7): 2049-2057.
Omran H, Haffner K, Volkel A, Kuehr J, Ketelsen UP, Ross UH, Konietzko N, Wienker T, Brandis M, Hildebrandt F: Homozygosity mapping of a gene locus for primary ciliary dyskinesia on chromosome 5p and identification of the heavy dynein chain DNAH5 as a candidate gene. Am J Respir Cell Mol Biol. 2000, 23 (5): 696-702.
Zariwala MA, Knowles MR, Omran H: Genetic defects in ciliary structure and function. Annu Rev Physiol. 2007, 69: 423-450.
Kong BW, Lee JY, Bottje WG, Lassiter K, Lee J, Foster DN: Genome-wide differential gene expression in immortalized DF-1 chicken embryo fibroblast cell line. BMC Genomics. 2011, 12: 571-2164-12-571-
This work was supported by Arkansas Bioscience Institute and Arkansas Agricultural Experimental Station.
The authors declare that they have no competing interests.
HMJ performed experiments, analyzed the data, and wrote the manuscript. GFE maintains the animal model, prepared the BL and SL blood samples, and edited the manuscript. KCR prepared DNA samples from blood, designed primers and performed PCR and Sanger sequencing reactions. BWK designed the experiments, performed experiments and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Gene name and functions of 139 gene containing amino acid changes showing ≥10 depth counts in SL chicken. Gene symbol, Entrez gene name, cellular location and type were provided. (XLSX 18 KB)
Additional file 4: List of focus molecules in gene networks. Gene symbols and Entrez gene names were displayed for the illustrations of network analysis. (XLSX 17 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Jang, HM., Erf, G.F., Rowland, K.C. et al. Genome resequencing and bioinformatic analysis of SNP containing candidate genes in the autoimmune vitiligo Smyth line chicken model. BMC Genomics 15, 707 (2014). https://doi.org/10.1186/1471-2164-15-707
- Amino Acid Change
- Single Nucleotide Polymorphism
- Ingenuity Pathway Analysis
- Read Depth