DNA Extraction and Amplification from Contemporary Polynesian Bark-Cloth

Background

Paper mulberry has been used for thousands of years in Asia and Oceania for making paper and bark-cloth, respectively. Museums around the world hold valuable collections of Polynesian bark-cloth. Genetic analysis of the plant fibers from which the textiles were made may answer a number of questions of interest related to provenance, authenticity or species used in the manufacture of these textiles. Recovery of nucleic acids from paper mulberry bark-cloth has not been reported before.

Methodology

We describe a simple method for the extraction of PCR-amplifiable DNA from small samples of contemporary Polynesian bark-cloth (tapa) using two types of nuclear markers. We report the amplification of about 300 bp sequences of the ITS1 region and of a microsatellite marker.

Conclusions

Sufficient DNA was retrieved from all bark-cloth samples to permit successful PCR amplification. This method shows a means of obtaining useful genetic information from modern bark-cloth samples and opens perspectives for the analyses of small fragments derived from ethnographic materials.     

The local origin of the Tibetan pig and additional insights into the origin of Asian pigs

Background

The domestic pig currently indigenous to the Tibetan highlands is supposed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. However, there is no direct genetic evidence of either the local origin or exotic migration of the Tibetan pig.

Methods and Findings

We analyzed mtDNA hypervariable segment I (HVI) variation of 218 individuals from seven Tibetan pig populations and 1,737 reported mtDNA sequences from domestic pigs and wild boars across Asia. The Bayesian consensus tree revealed a main haplogroup M and twelve minor haplogroups, which suggested a large number of small scale in situ domestication episodes. In particular, haplogroups D1 and D6 represented two highly divergent lineages in the Tibetan highlands and Island Southeastern Asia, respectively. Network analysis of haplogroup M further revealed one main subhaplogroup M1 and two minor subhaplogroups M2 and M3. Intriguingly, M2 was mainly distributed in Southeastern Asia, suggesting for a local origin. Similar with haplogroup D6, M3 was mainly restricted in Island Southeastern Asia. This pattern suggested that Island Southeastern Asia, but not Southeastern Asia, might be the center of domestication of the so-called Pacific clade (M3 and D6 here) described in previous studies. Diversity gradient analysis of major subhaplogroup M1 suggested three local origins in Southeastern Asia, the middle and downstream regions of the Yangtze River, and the Tibetan highlands, respectively.

Conclusions

We identified two new origin centers for domestic pigs in the Tibetan highlands and in the Island Southeastern Asian region.     

Control region variability of haplogroup C1d and the tempo of the peopling of the Americas

Background

Among the founding mitochondrial haplogroups involved in the peopling of the Americas, haplogroup C1d has been viewed as problematic because of its phylogeny and because of the estimates of its antiquity, apparently being much younger than other founding haplogroups. Several recent analyses, based on data from the entire mitochondrial genome, have contributed to an advance in the resolution of these problems. The aim of our analysis is to compare the conclusions drawn from the available HVR-I and HVR-II data for haplogroup C1d with the ones based on whole mitochondrial genomes.

Methodology/Principal Findings

HVR-I and HVR-II sequences defined as belonging to haplogroup C1d by standard criteria were gathered from the literature as well as from population studies carried out in Uruguay. Sequence phylogeny was reconstructed using median-joining networks, geographic distribution of lineages was analyzed and the age of the most recent common ancestor estimated using the ρ-statistic and two different mutation rates. The putative ancestral forms of the haplogroup were found to be more widespread than the derived lineages, and the lineages defined by np 194 were found to be widely distributed and of equivalent age.

Conclusions/Significance

The analysis of control region sequences is found to still harbor great potential in tracing microevolutionary phenomena, especially those found to have occurred in more recent times. Based on the geographic distributions of the alleles of np 7697 and np 194, both discussed as possible basal mutations of the C1d phylogeny, we suggest that both alleles were part of the variability of the haplogroup at the time of its entrance. Moreover, based on the mutation rates of the different sites stated to be diagnostic, it is possible that the anomalies found when analyzing the haplogroup are due to paraphyly.     

Evolutionary Pressure on Mitochondrial Cytochrome b Is Consistent with a Role of CytbI7T Affecting Longevity during Caloric Restriction

Background

Metabolism of energy nutrients by the mitochondrial electron transport chain (ETC) is implicated in the aging process. Polymorphisms in core ETC proteins may have an effect on longevity. Here we investigate the cytochrome b (cytb) polymorphism at amino acid 7 (cytbI7T) that distinguishes human mitochondrial haplogroup H from haplogroup U.

Principal Findings

We compared longevity of individuals in these two haplogroups during historical extremes of caloric intake. Haplogroup H exhibits significantly increased longevity during historical caloric restriction compared to haplogroup U (p = 0.02) while during caloric abundance they are not different. The historical effects of natural selection on the cytb protein were estimated with the software TreeSAAP using a phylogenetic reconstruction for 107 mammal taxa from all major mammalian lineages using 13 complete protein-coding mitochondrial gene sequences. With this framework, we compared the biochemical shifts produced by cytbI7T with historical evolutionary pressure on and near this polymorphic site throughout mammalian evolution to characterize the role cytbI7T had on the ETC during times of restricted caloric intake.

Significance

Our results suggest the relationship between caloric restriction and increased longevity in human mitochondrial haplogroup H is determined by cytbI7T which likely enhances the ability of water to replenish the Q_i binding site and decreases the time ubisemiquinone is at the Q_o site, resulting in a decrease in the average production rate of radical oxygen species (ROS).     

Deep Rooting In-Situ Expansion of mtDNA Haplogroup R8 in South Asia

Background

The phylogeny of the indigenous Indian-specific mitochondrial DNA (mtDNA) haplogroups have been determined and refined in previous reports. Similar to mtDNA superhaplogroups M and N, a profusion of reports are also available for superhaplogroup R. However, there is a dearth of information on South Asian subhaplogroups in particular, including R8. Therefore, we ought to access the genealogy and pre-historic expansion of haplogroup R8 which is considered one of the autochthonous lineages of South Asia.

Methodology/Principal Findings

Upon screening the mtDNA of 5,836 individuals belonging to 104 distinct ethnic populations of the Indian subcontinent, we found 54 individuals with the HVS-I motif that defines the R8 haplogroup. Complete mtDNA sequencing of these 54 individuals revealed two deep-rooted subclades: R8a and R8b. Furthermore, these subclades split into several fine subclades. An isofrequency contour map detected the highest frequency of R8 in the state of Orissa. Spearman''s rank correlation analysis suggests significant correlation of R8 occurrence with geography.

Conclusions/Significance

The coalescent age of newly-characterized subclades of R8, R8a (15.4±7.2 Kya) and R8b (25.7±10.2 Kya) indicates that the initial maternal colonization of this haplogroup occurred during the middle and upper Paleolithic period, roughly around 40 to 45 Kya. These results signify that the southern part of Orissa currently inhabited by Munda speakers is likely the origin of these autochthonous maternal deep-rooted haplogroups. Our high-resolution study on the genesis of R8 haplogroup provides ample evidence of its deep-rooted ancestry among the Orissa (Austro-Asiatic) tribes.     

Host Factors Determine Anti-GM1 Response Following Oral Challenge of Chickens with Guillain-Barré Syndrome Derived Campylobacter jejuni Strain GB11

Background

Anti-ganglioside antibodies with a pathogenic potential are present in C. jejuni-associated Guillain-Barré syndrome (GBS) patients and are probably induced by molecular mimicry. Immunization studies in rabbits and mice have demonstrated that these anti-ganglioside antibodies can be induced using purified lipo-oligosaccharides (LOS) from C. jejuni in a strong adjuvant.

Methodology/Principal Findings

To investigate whether natural colonization of chickens with a ganglioside-mimicking C. jejuni strain induces an anti-ganglioside response, and to investigate the diversity in anti-ganglioside response between and within genetically different chicken lines, we orally challenged chickens with different C. jejuni strains. Oral challenge of chickens with a C. jejuni strain from a GBS patient, containing a LOS that mimics ganglioside GM1, induced specific IgM and IgG anti-LOS and anti-GM1 antibodies. Inoculation of chickens with the Penner HS:3 serostrain, without a GM1-like structure, induced anti-LOS but no anti-ganglioside antibodies. We observed different patterns of anti-LOS/ganglioside response between and within the five strains of chickens.

Conclusions

Natural infection of chickens with C. jejuni induces anti-ganglioside antibodies. The production of antibodies is governed by both microbial and host factors.     

Functional Differences between Mitochondrial Haplogroup T and Haplogroup H in HEK293 Cybrid Cells

Background

Epidemiological case-control studies have revealed associations between mitochondrial haplogroups and the onset and/or progression of various multifactorial diseases. For instance, mitochondrial haplogroup T was previously shown to be associated with vascular diseases, including coronary artery disease and diabetic retinopathy. In contrast, haplogroup H, the most frequent haplogroup in Europe, is often found to be more prevalent in healthy control subjects than in patient study groups. However, justifications for the assumption that haplogroups are functionally distinct are rare. Therefore, we attempted to compare differences in mitochondrial function between haplogroup H and T cybrids.

Methodology/Principal Findings

Mitochondrial haplogroup H and T cybrids were generated by fusion of HEK293 cells devoid of mitochondrial DNA with isolated thrombocytes of individuals with the respective haplogroups. These cybrid cells were analyzed for oxidative phosphorylation (OXPHOS) enzyme activities, mitochondrial DNA (mtDNA) copy number, growth rate and susceptibility to reactive oxygen species (ROS). We observed that haplogroup T cybrids have higher survival rate when challenged with hydrogen peroxide, indicating a higher capability to cope with oxidative stress.

Conclusions/Significance

The results of this study show that functional differences exist between HEK293 cybrid cells which differ in mitochondrial genomic background.     

The Druze: a population genetic refugium of the Near East

Background

Phylogenetic mitochondrial DNA haplogroups are highly partitioned across global geographic regions. A unique exception is the X haplogroup, which has a widespread global distribution without major regions of distinct localization.

Principal Findings

We have examined mitochondrial DNA sequence variation together with Y-chromosome-based haplogroup structure among the Druze, a religious minority with a unique socio-demographic history residing in the Near East. We observed a striking overall pattern of heterogeneous parental origins, consistent with Druze oral tradition, together with both a high frequency and a high diversity of the mitochondrial DNA (mtDNA) X haplogroup within a confined regional subpopulation. Furthermore demographic modeling indicated low migration rates with nearby populations.

Conclusions

These findings were enabled through the use of a paternal kindred based sampling approach, and suggest that the Galilee Druze represent a population isolate, and that the combination of a high frequency and diversity of the mtDNA X haplogroup signifies a phylogenetic refugium, providing a sample snapshot of the genetic landscape of the Near East prior to the modern age.     

Geographic Variation in Venom Allelic Composition and Diets of the Widespread Predatory Marine Gastropod Conus ebraeus

Background

Members of the predatory gastropod genus Conus use a venom comprised of a cocktail of peptide neurotoxins, termed conotoxins or conopeptides, to paralyze prey and conotoxin gene family members diversify via strong positive selection. Because Conus venoms are used primarily to subdue prey, the evolution of venoms is likely affected by predator-prey interactions.

Methodology/Principal Findings

To identify the selective forces that drive the differentiation of venoms within species of Conus, we examined the distribution of alleles of a polymorphic O-superfamily conotoxin locus of Conus ebraeus at Okinawa, Guam and Hawaii. Previous analyses of mitochondrial cytochrome oxidase I gene sequences suggest that populations of C. ebraeus, a worm-eating Conus, are not structured genetically in the western and central Pacific. Nonetheless, because the sample size from Guam was relatively low, we obtained additional data from this location and reexamined patterns of genetic variation at the mitochondrial gene at Okinawa, Guam and Hawaii. We also utilized a DNA-based approach to identify prey items of individuals of C. ebraeus from Guam and compared this information to published data on diets at Okinawa and Hawaii. Our results show that conotoxin allelic frequencies differ significantly among all three locations, with strongest differentiation at Hawaii. We also confirm previous inferences that C. ebraeus exhibits no genetic differentiation between Okinawa, Guam and Hawaii at the mitochondrial locus. Finally, DNA-based analyses show that eunicid polychaetes comprise the majority of the prey items of C. ebraeus at Guam; while this results compares well with observed diet of this species at Okinawa, C. ebraeus preys predominantly on nereid polychaetes at Hawaii.

Conclusions/Significance

These results imply that strong selection pressures affect conotoxin allelic frequencies. Based on the dietary information, the selection may derive from geographic variation in dietary specialization and local coevolutionary arms races between Conus and their prey.     

Domestic Animal Hosts Strongly Influence Human-Feeding Rates of the Chagas Disease Vector Triatoma infestans in Argentina

Background

The host species composition in a household and their relative availability affect the host-feeding choices of blood-sucking insects and parasite transmission risks. We investigated four hypotheses regarding factors that affect blood-feeding rates, proportion of human-fed bugs (human blood index), and daily human-feeding rates of Triatoma infestans, the main vector of Chagas disease.

Methods

A cross-sectional survey collected triatomines in human sleeping quarters (domiciles) of 49 of 270 rural houses in northwestern Argentina. We developed an improved way of estimating the human-feeding rate of domestic T. infestans populations. We fitted generalized linear mixed-effects models to a global model with six explanatory variables (chicken blood index, dog blood index, bug stage, numbers of human residents, bug abundance, and maximum temperature during the night preceding bug catch) and three response variables (daily blood-feeding rate, human blood index, and daily human-feeding rate). Coefficients were estimated via multimodel inference with model averaging.

Findings

Median blood-feeding intervals per late-stage bug were 4.1 days, with large variations among households. The main bloodmeal sources were humans (68%), chickens (22%), and dogs (9%). Blood-feeding rates decreased with increases in the chicken blood index. Both the human blood index and daily human-feeding rate decreased substantially with increasing proportions of chicken- or dog-fed bugs, or the presence of chickens indoors. Improved calculations estimated the mean daily human-feeding rate per late-stage bug at 0.231 (95% confidence interval, 0.157–0.305).

Conclusions and Significance

Based on the changing availability of chickens in domiciles during spring-summer and the much larger infectivity of dogs compared with humans, we infer that the net effects of chickens in the presence of transmission-competent hosts may be more adequately described by zoopotentiation than by zooprophylaxis. Domestic animals in domiciles profoundly affect the host-feeding choices, human-vector contact rates and parasite transmission predicted by a model based on these estimates.     

Ashkenazi Jewish centenarians do not demonstrate enrichment in mitochondrial haplogroup J

Background

Association of mitochondrial haplogroup J with longevity has been reported in several population subgroups. While studies from northern Italy and Finland, have described a higher frequency of haplogroup J among centenarians in comparison to non-centenarian, several other studies could not replicate these results and suggested various explanations for the discrepancy.

Methodology/Principal Findings

We have evaluated haplogroup frequencies among Ashkenazi Jewish centenarians using two different sets of matched controls. No difference was observed in the haplogroup J frequencies between the centenarians or either matched control group, despite adequate statistical power to detect such a difference. Furthermore, the lack of association was robust to population substructure in the Ashkenazi Jewish population. Given this discrepancy with the previous reported associations in the northern Italian and the Finnish populations, we conducted re-analysis of these previously published data, which supported one of several possible explanations: i) inadequate matching of cases and controls; ii) inadequate adjustment for multiple comparison testing; iii) cryptic population stratification.

Conclusions/Significance

There does not exist a universal association of mitochondrial haplogroup J with longevity across all population groups. Reported associations in specialized populations may reflect genetic or other interactions specific to those populations or else cryptic confounding influences, such as inadequate matching attributable to population substructure, which are of general relevance to all studies of the possible association of mitochondrial DNA haplogroups with common complex phenotypes.     

Genetic affinities of the central Indian tribal populations

Background

The central Indian state Madhya Pradesh is often called as ‘heart of India’ and has always been an important region functioning as a trinexus belt for three major language families (Indo-European, Dravidian and Austroasiatic). There are less detailed genetic studies on the populations inhabited in this region. Therefore, this study is an attempt for extensive characterization of genetic ancestries of three tribal populations, namely; Bharia, Bhil and Sahariya, inhabiting this region using haploid and diploid DNA markers.

Methodology/Principal Findings

Mitochondrial DNA analysis showed high diversity, including some of the older sublineages of M haplogroup and prominent R lineages in all the three tribes. Y-chromosomal biallelic markers revealed high frequency of Austroasiatic-specific M95-O2a haplogroup in Bharia and Sahariya, M82-H1a in Bhil and M17-R1a in Bhil and Sahariya. The results obtained by haploid as well as diploid genetic markers revealed strong genetic affinity of Bharia (a Dravidian speaking tribe) with the Austroasiatic (Munda) group. The gene flow from Austroasiatic group is further confirmed by their Y-STRs haplotype sharing analysis, where we determined their founder haplotype from the North Munda speaking tribe, while, autosomal analysis was largely in concordant with the haploid DNA results.

Conclusions/Significance

Bhil exhibited largely Indo-European specific ancestry, while Sahariya and Bharia showed admixed genetic package of Indo-European and Austroasiatic populations. Hence, in a landscape like India, linguistic label doesn''t unequivocally follow the genetic footprints.     

Mitochondrial haplogroups and control region polymorphisms in age-related macular degeneration: a case-control study

Background

Onset and development of the multifactorial disease age-related macular degeneration (AMD) are highly interrelated with mitochondrial functions such as energy production and free radical turnover. Mitochondrial dysfunction and overproduction of reactive oxygen species may contribute to destruction of the retinal pigment epithelium, retinal atrophy and choroidal neovascularization, leading to AMD. Consequently, polymorphisms of the mitochondrial genome (mtDNA) are postulated to be susceptibility factors for this disease. Previous studies from Australia and the United States detected associations of mitochondrial haplogroups with AMD. The aim of the present study was to test these associations in Middle European Caucasians.

Methodology/Principal Findings

Mitochondrial haplogroups (combinations of mtDNA polymorphisms) and mitochondrial CR polymorphisms were analyzed in 200 patients with wet AMD (choroidal neovascularization, CNV), in 66 patients with dry AMD, and in 385 controls from Austria by means of multiplex primer extension analysis and sequencing, respectively. In patients with CNV, haplogroup H was found to be significantly less frequent compared to controls, and haplogroup J showed a trend toward a higher frequency compared to controls. Five CR polymorphisms were found to differ significantly in the two study populations compared to controls, and all, except one (T152C), are linked to those haplogroups.

Conclusions/Significance

It can be concluded that haplogroup J is a risk factor for AMD, whereas haplogroup H seems to be protective for AMD.     

Accuracy of estimated breeding values with genomic information on males,females, or both: an example on broiler chicken

Background

As more and more genotypes become available, accuracy of genomic evaluations can potentially increase. However, the impact of genotype data on accuracy depends on the structure of the genotyped cohort. For populations such as dairy cattle, the greatest benefit has come from genotyping sires with high accuracy, whereas the benefit due to adding genotypes from cows was smaller. In broiler chicken breeding programs, males have less progeny than dairy bulls, females have more progeny than dairy cows, and most production traits are recorded for both sexes. Consequently, genotyping both sexes in broiler chickens may be more advantageous than in dairy cattle.

Methods

We studied the contribution of genotypes from males and females using a real dataset with genotypes on 15 723 broiler chickens. Genomic evaluations used three training sets that included only males (4648), only females (8100), and both sexes (12 748). Realized accuracies of genomic estimated breeding values (GEBV) were used to evaluate the benefit of including genotypes for different training populations on genomic predictions of young genotyped chickens.

Results

Using genotypes on males, the average increase in accuracy of GEBV over pedigree-based EBV for males and females was 12 and 1 percentage points, respectively. Using female genotypes, this increase was 1 and 18 percentage points, respectively. Using genotypes of both sexes increased accuracies by 19 points for males and 20 points for females. For two traits with similar heritabilities and amounts of information, realized accuracies from cross-validation were lower for the trait that was under strong selection.

Conclusions

Overall, genotyping males and females improves predictions of all young genotyped chickens, regardless of sex. Therefore, when males and females both contribute to genetic progress of the population, genotyping both sexes may be the best option.     

The Peopling of Korea Revealed by Analyses of Mitochondrial DNA and Y-Chromosomal Markers

Background

The Koreans are generally considered a northeast Asian group because of their geographical location. However, recent findings from Y chromosome studies showed that the Korean population contains lineages from both southern and northern parts of East Asia. To understand the genetic history and relationships of Korea more fully, additional data and analyses are necessary.

Methodology and Results

We analyzed mitochondrial DNA (mtDNA) sequence variation in the hypervariable segments I and II (HVS-I and HVS-II) and haplogroup-specific mutations in coding regions in 445 individuals from seven east Asian populations (Korean, Korean-Chinese, Mongolian, Manchurian, Han (Beijing), Vietnamese and Thais). In addition, published mtDNA haplogroup data (N = 3307), mtDNA HVS-I sequences (N = 2313), Y chromosome haplogroup data (N = 1697) and Y chromosome STR data (N = 2713) were analyzed to elucidate the genetic structure of East Asian populations. All the mtDNA profiles studied here were classified into subsets of haplogroups common in East Asia, with just two exceptions. In general, the Korean mtDNA profiles revealed similarities to other northeastern Asian populations through analysis of individual haplogroup distributions, genetic distances between populations or an analysis of molecular variance, although a minor southern contribution was also suggested. Reanalysis of Y-chromosomal data confirmed both the overall similarity to other northeastern populations, and also a larger paternal contribution from southeastern populations.

Conclusion

The present work provides evidence that peopling of Korea can be seen as a complex process, interpreted as an early northern Asian settlement with at least one subsequent male-biased southern-to-northern migration, possibly associated with the spread of rice agriculture.     

Relationship between seasonal cold acclimatization and mtDNA haplogroup in Japanese

Background

The purpose of this study was to elucidate the interaction between mtDNA haplogroup and seasonal variation that contributes to cold adaptation.

Methods

There were 15 subjects (seven haplotype D subjects and eight haplotype non-D subjects). In summer and winter, the subjects were placed in an environment where the ambient temperature dropped from 27 °C to 10 °C in 30 minutes. After that, they were exposed to cold for 60 minutes.

Results

In summer, the decrease in rectal temperature and increase in oxygen consumption was smaller and cold tolerance was higher in the haplotype non-D group than in the haplotype D group. In winter, no significant differences were seen in rectal temperature or oxygen consumption, but the respiratory exchange ratio decreased in the haplotype D group.

Conclusions

The results of the present study suggest that haplogroup D subjects are a group that changes energy metabolism more, and there appears to be a relationship between differences in cold adaptability and mtDNA polymorphism within the population. Moreover, group differences in cold adaptability seen in summer may decrease in winter due to supplementation by seasonal cold acclimatization.     

Trypanosoma cruzi in the chicken model: Chagas-like heart disease in the absence of parasitism

Background

The administration of anti-trypanosome nitroderivatives curtails Trypanosoma cruzi infection in Chagas disease patients, but does not prevent destructive lesions in the heart. This observation suggests that an effective treatment for the disease requires understanding its pathogenesis.

Methodology/Principal Findings

To understand the origin of clinical manifestations of the heart disease we used a chicken model system in which infection can be initiated in the egg, but parasite persistence is precluded. T. cruzi inoculation into the air chamber of embryonated chicken eggs generated chicks that retained only the parasite mitochondrial kinetoplast DNA minicircle in their genome after eight days of gestation. Crossbreeding showed that minicircles were transferred vertically via the germ line to chicken progeny. Minicircle integration in coding regions was shown by targeted-primer thermal asymmetric interlaced PCR, and detected by direct genomic analysis. The kDNA-mutated chickens died with arrhythmias, shortness of breath, cyanosis and heart failure. These chickens with cardiomyopathy had rupture of the dystrophin and other genes that regulate cell growth and differentiation. Tissue pathology revealed inflammatory dilated cardiomegaly whereby immune system mononuclear cells lyse parasite-free target heart fibers. The heart cell destruction implicated a thymus-dependent, autoimmune; self-tissue rejection carried out by CD45⁺, CD8γδ⁺, and CD8α lymphocytes.

Conclusions/Significance

These results suggest that genetic alterations resulting from kDNA integration in the host genome lead to autoimmune-mediated destruction of heart tissue in the absence of T. cruzi parasites.     

Genome resequencing and bioinformatic analysis of SNP containing candidate genes in the autoimmune vitiligo Smyth line chicken model

Background

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.

Results

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).

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-707) contains supplementary material, which is available to authorized users.