- Open Access
Genome-wide association study for posthitis in the free-living population of European bison (Bison bonasus)
Biology Direct volume 10, Article number: 2 (2015)
About 5–6% of the European bison (Bison bonasus) males are affected by posthitis (necrotic inflammation of the prepuce) and die in the wild forest. Despite many years of study, pathogenesis of this disease has not yet been determined. The main aim of the study was to find SNP markers significantly associated with the incidence of posthitis and mine the genome for candidate genes potentially involved in the development of the disease.
It was shown that relatively small number of SNPs effects reached genome-wide significance after false discovery rate (FDR) correction. Among 25 significant markers, the highest effects were found for two SNPs (rs110456748 and rs136792896) located at the distance of 23846 bp and 37742 bp, respectively, from OR10A3 gene (olfactory receptor genes), known to be involved in atopic dermatitis in humans. It was also observed that five other significant SNP markers were located in the proximity of candidate genes involved in severe diseases of skin tissue and cancer/tumour development of epithelial or testicular germ cells, which suggest their potential participation in the posthitis. The 25 investigated SNPs showed marked differences in allelic and genotypic frequencies between the healthy and affected bison groups.
The 2 Mb region of the BTA15 chromosome is involved in genetic background of posthitis and should be closer examined to find causal mutations helpful in better understanding of the disease ethology and to control its incidence in the future.
This article was reviewed by Prof. Lev Klebanov and Dr. Fyodor Kondrashov.
Genome Wide Association Studies (GWAS) have become one of the major tools in explaining the ethology of human diseases of genetic origin . The technique has also been applied in studies on model organisms  and domestic animals [3,4] as well as in crops . Performing GWAS in wild living species is usually difficult, as all the commercially available microarray sets are designed for model or domestic animals. It has been shown that cattle SNP microarrays work well in bison species [6-8]. The European Bison (Bison bonasus), (EB), after extinction in the wild before World War II, was restored based on progeny of just 12 individuals preserved in private collections and zoological gardens . The restored population was divided into two isolated genetic lines: the Lowland or the Bialowieza line (in which the founders were seven pure Lowland bison) and the Lowland-Caucasian line (in which the founders were 11 Lowland bison individuals and one Caucasian bison bull). Both populations of the EB are extremely inbred .
In the last four decades, a severe disease called posthitis affected EB males. This disease was discovered in the 1980s in the Bialowieza Forest [11,12]. At the end of the 1990s, similar symptoms were also observed in five young EB of the Lowland line from the Bayerisher Wald National Park, Germany . Posthitis is defined as necrotic inflammation of the prepuce. About 5–6% of the males are affected each year and die in the wild forest or are eliminated by officially approved hunting . Despite many years of study, pathogenesis of posthitis in bison has not yet been determined. The affected bulls do not exhibit changes in the general physiological mechanism as indicated in the long-term studies of 30 physiological indices . Posthitis is thought to be a disease of an endemic character since it is only observed in males living in Bialowieza National Park or transported from it. Apparently, no environmental or genetic factors which would explain the susceptibility of some individuals to develop the disease have been identified . New tools for whole genome analysis (SNP microarrays) available for the Bovidae species since 2009 (www.illumina.com) have opened new possibilities to improve the genetic analysis of posthitis.
Among the commercially available microarrays, the Illumina Bovine HD 777 K microarray is the most suitable tool for bison genome studies. Its usefulness was evaluated recently by showing that the Lowland line is a genetically homogenous population, with only a small amount of polymorphic loci, which can be helpful in pedigree analysis [7,16]. It was also shown that only 2.8% out of 54,000 bovine point mutations are polymorphic in ten EB representing the Caucasian and Lowland lines .
In this paper, we used high density bovine SNP microarray to investigate the entire bison genome in order to identify differences in SNP allele and genotype frequencies between the group of individuals affected by posthitis and the non-affected group. The differences in allelic frequency found between these two groups were used to perform a GWAS to identify genes potentially involved in the posthitis disease.
A global view of all SNPs effects for posthitis incidence is presented in the Manhattan plot in Figure 1. It is shown that a relatively small number of SNPs effects reached genome-wide significance after FDR correction. Among the 25 top significant markers, most of them (18) were located on chromosome 15. Single significant markers were also found on chromosome 3, 9, 13, 26 and three markers on chromosome 12.
Detailed characterization of the top 25 markers are given in Table 1, including the SNP position on UMD 3.1 genome map, the name of the closest candidate gene and its distance from the significant marker and the reference confirming the potential involvement of a candidate gene in posthitis. Eleven significant markers were located on chromosome 15 in the close vicinity of olfactory receptor genes.
To better show how the number of SNP genotypes change between affected and unaffected bison, the allelic and genotypic frequencies of the 25 markers were estimated for the two bison groups and are listed in Table 2.
For two randomly chosen significant markers, SNP genotypes designated by Illumina Bead Studio were checked by sequencing (rs136792896and rs137427505). No discrepancies between the alternative genotypes (AA vs BB) for these SNPs were found (Figure 2).
European bison is an example of an extinct species, which was successfully reintroduced to wild nature. Due to the founder effect, the current free-living bison population of the Lowland line (about 1000 individuals in 2011 ) is extremely inbred and, therefore special attention should be paid to any threats decreasing the reproductive potential of the population and, consequently, lowering the genetic diversity of the species .
The inbreeding level of EB was also considered as a potential factor increasing the risk of disease, but the difference between the inbreeding coefficient of healthy (32 bison) and affected (12 bison) was not found to be significant . Another analysis of 69 affected males did not find significant associations between MHC DRB3 alleles and the susceptibility to posthitis .
In this study, the utilization of an HD (High Density) microarray, allowed a detailed screening of 18,079 polymorphic markers which were used in a GWAS analysis to identify those significantly associated with the incidence of posthitis. In Figure 1 and Table 1 it is shown that a group of markers located on chromosome 15 have the highest effect. They are located within OLFML1 and PPFIBP2 and in the close vicinity of the ZNF214 gene. Most significant effects were found for two SNPs (rs110456748 and rs136792896) located at a distance of 23846 bp and 37742 bp from OR10A3 gene, known as being involved in atopic dermatitis, a common inflammatory skin disease .
Several significant markers were mapped within two genes: PPFIBP2 and OLFML1. PPFIBP2 encodes protein-tyrosine phosphatase-interacting protein, which is differentially expressed in endometrial cancer  and is also involved in tumour cell migration and invasiveness of extracellular signal-regulated kinase depleted cells . OLFML1 (olfactomedin-like protein 1) encodes a secreted protein, is expressed in many tissues and may play a significant role in the regulation of cell proliferation in vitro . Another gene located closely to significant marker rs134495844 is coding zing finger protein 214 (ZNF214) - a potential tumour suppressor gene whose expression is increased in tumours derived from testicular germ cells . The next candidate genes, COG6 and LHFP, were found to be associated with susceptibility to psoriasis [25,26]. Among the genes listed in Table 1, five candidate genes are associated with severe skin diseases (atopic dermatitis, psoriasis) or cancer/tumour development of epithelial or testicular germ cells, which suggest their potential participation in the posthitis.
The fact that the remaining five candidate genes (Table 1) do not show a direct link with an inflammation/necrosis specific function does not exclude the possibility that these genes are contributing to pothistis susceptibility. It is possible that their contribution is indirect, through immunological pathways, or has not yet been reported. It can be assumed that, as for any complex trait, many genes whose names are not “trait-specific” participate in the overall genetic variance of the trait .
Being aware that in typical GWAS, most of the significant markers are located outside the functional genes and are in linkage disequilibrium with other genes, we additionally used the MapViewer service to identify genes located in the broader neighborhood of the significant marker (at a distance of less than 1 Mb). The size of the inspected region is difficult to be justified objectively since the distribution of the functional gene along the chromosomes is very different and specific to each autosome . One Mb distance from a significant marker seems to be a manageable region to manually screen all genes included. All such genes were checked for functional association with inflammation/tumour/cancer/necrosis. In this way, we were able to find further candidate genes related to the trait of interest. For example, COL17A1 encoding collagen type XVII, alpha 1 (located 529,5 Kb from the significant marker rs110088444, on BTA 26) is involved in epidermolysis bullosa, causing separation of the basal keratinocyte from the underlying basement membrane and consequently skin atrophy . TRIM6 coding tripartite motif 6 located 3799 bp from the significant marker rs1986306 has been reported as playing a role in an IFN-induced antiviral state against retrovirus infection . ILK gene encoding integrin-linked kinase (mapped 477,9 K from significant marker rs134156378, BTA15) interacts with pathogenic bacteria reinforcing host cell adherence . Another gene associated with skin pathology, FERMT1 encoding fermitin family member 1 is located 34,7Kb from significant marker rs137458105 (BTA13). Its mutation causes Kindler syndrome, whose symptoms include skin blistering, skin atrophy, photosensitivity, colonic inflammation and mucosal stenosis . Proliferating cell nuclear antigen (PCNA) is located about 787Kb from marker rs137458105. It was shown that the percentage of PCNA-positive cells in malignant and some non-malignant skin diseases (atopic dermatitis, psoriasis) is from 2- to 5- times higher than in normal skin cells .
In the close vicinity of significant marker rs134255411 (BTA 3) CSF3R gene encoding colony stimulating factor 3 receptor is located. Its mutation causes chronic neutrophilia . It was reported that periostin, osteoblast specific factor, a protein encoded by the POSTN gene exacerbates the pathogenesis of atopic dermatitis in mice [35,36]. POSTN is located 618Kb from the significant marker rs137427505 on BTA12. Another gene located near the significant marker rs137024050 is forkhead box O1 transcription factor, whose loss of expression in lesional psoriatic tissues has a potential contribution to the development of psoriasis .
All of the above-mentioned reports support the main outcome of our work – that the most significant markers for posthitis indicate regions containing genes involved in the immune system and disorders of skin and epithelial cells. Unlike in Johnston , we used marker set of extremely high coverage of the cattle genome – over 770,000 loci. Attempts to use standard microchips (of approx. 54,000 SNPs) failed as we could not find any marker associated with the disease symptoms (data not shown). In future research, GWAS should be validated for a larger population of posthitis-affected bison and candidate genes should be screened to find causal mutations enabling a better understanding of the disease and, prospectively, to control its incidence. The fact that the 25 SNPs investigated showed marked differences in allelic and genotypic frequencies between the healthy and affected bison groups is a clear indication of an excess of the genotype BB in the affected population together with a higher frequency of the B allele in the affected compared to the non-affected bison population. In addition, the fact that many alleles showed the same allelic frequency indicate a strong linkage between several SNPs on the same chromosome. This suggests that an appropriate breeding strategy based on identity by descent (IBD) could be developed for the European bison. Information from the SNP chip could be used as an accurate tool for guiding which individuals should mate, and which should be excluded from the mating as they are bearing the alleles associated with incidence of posthitis. Given the fact that many of the SNPs which showed associations are strongly linked, it would be enough to simply screen one of the SNPs which are linked to each other. At the same time it would also be desirable to optimize the best contribution of animals to the next generation, in order to reduce the pace at which the inbreeding coefficient is increasing every generation. This can be in the form of a specific list with suggested matings or guidelines on the number of matings that a given bison should be engaged in during a given number of generations. Developing a SNP chip with a subset of SNPs found to be associated with the posthitis incidence and with a subset of SNPs informative for paternity analysis and to estimate the kinship between individuals could be useful for this purpose and make these goals more realistic from an economic point of view.
Allele frequencies of 25 SNP markers are significantly different in health and posthitis-affected animals. Genes localized in the closest vicinity of these markers show connections with immunological processes. 18 out of 25 markers are localized on chromosome BTA15 within the space of 2 Mb, which suggests this region may play an important role in the ethology of the posthitis. This subregion should be closer examined to find causal mutations helpful in better understanding of the genetic background of the disease and to control its incidence in the future.
The study included 90 adult males of the Lowland line of the EB living in the Bialowieza National Park. Bison DNA samples were provided by the Mammal Research Institute in Bialowieza (Polish Academy of Sciences). The animals were divided into two groups: healthy (36) and affected by posthitis (54). The presence of the disease and the degree of its advancement were ascertained post mortem by an expert veterinarian, on the basis of anatomical and pathological changes  (Figure 3).
Genotyping and statistical analysis
Each individual was genotyped using the Illumina BovineHD 777 K microarray, which consists of 777,962 SNP markers. Associations were tested using GoldenHelix SVS7 analysis software (Bozeman, MT). Initial data clean-up was performed to remove poorly performing SNPs. SNPs in the control population with significant deviations from the Hardy–Weinberg equilibrium (HWE) at the significance level of P < 0.0001 were removed prior to association analysis. Additionally, the SNP selection criteria were applied based on a minor allele frequency (MAF) of at least MAF > 0.005, due to highly monomorphic data and the technical quality of the chip – a minimum call rate of 98%. After filtering, 18,079 SNPs were used in the final analysis. The density of SNPs used in the analysis was 1 SNP per 138 kb. The average call rate obtained for our data set amounted to 98.86% (SD = 0.2%) for selected SNPs.
An additive model of the Linear Regression Analysis with a False Discovery Rate, kinship matrix (computed on the autosomal SNPs in SVS7.6) and heterozygosity rate as additional covariate was performed to estimate the effect of an SNP marker for posthitis incidence. SNP association analysis was performed on autosomal SNPs. Following this, the allelic and genotypic frequencies of the significant SNP markers were then estimated for the two bison groups.
Since bison DNA was used on a bovine SNP array, two fragments of bison genomic DNA (containing two SNPs) were amplified and sequenced to check whether the SNP genotype was correct. The following SNPs were PCR-amplified and sequenced: BovineHD1500012672 (rs136792896) and BovineHD1200007116 (rs137427505). Primers were designed and blasted using the NCBI primer-blast website (http://www.ncbi.nlm.nih.gov/tools/primer-blast/):
In the PCR reaction, 20 pM of each primer, 1 U of Taq DNA polymerase (Bioron, Germany), 100 ng of genomic DNA, 2.5 μl of 10× PCR buffer, 2.5 mM of MgCl2 and 800 μM of dNTPs were used. The amplification program consisted of: an initial DNA denaturation (94°C/3 min), 35 cycles of denaturation (94°C/30 s), annealing (61°C/30 s) and elongation (72°C/30 s), final synthesis (72°C/5 min).
PCR yielded two specific amplicons of expected 185 (BovineHD1500012672) and 537 bp (BovineHD1200007116). The PCR products were purified and sequenced by the GENOMED Ltd (Poland). The sequences were analyzed using BioEdit v. 7.2.0 software.
Mining genome for candidate genes
Two approaches were used to locate candidate genes: statistical and functional. In the statistical approach, genes located in the closest vicinity of markers above the cut-off line were chosen (P value after FDR correction was 3.40). In the functional approach, 1 Mb region around the position of the significant SNP marker was inspected in the MapViewer service (www.ncbi.nlm.nih.gov) in order to check whether there are genes functionally associated with inflammation or diseases showing symptoms related to posthitis.
The genomic positions of SNPs included in the Illumina BovineHD 777 K were taken from the Illumina publication (www.illumina.com). Genomic positions of candidate genes were assigned based on the UMD 3.1 bovine genome assembly  through the current Ensembl database (www.ensembl.org).
Statement of ethical approval
Tissue samples were collected during annual culling of the European bison approved by General Directorate for Environmental Protection (Warsaw, Poland) and Regional Directorate for Environmental Protection (Bialystok, Poland) and released for Mammal Research Institute in Bialowieza. Reference numbers of the approvals: 2010: DOP-Pozgiz-4200/IV-57/4139/10/kp; 2011: DOP-OZGTZ.6401.05.5.2012.kp; 2012: WPN.6401.222.2012.EJ; 2013: WPN.6401.266.2013.EJ; 2014: WPN.6401.221.2014.EJ.
We appreciate very much all remarks, suggestions and comments provided by reviewers. We considered them carefully trying to improve the manuscript.
Reviewer 1: Dr. Fyodor Kondrashov
I do not see any issues in the design and implementation of the present study and the results are straightforward. The authors capitalized on the application of tools that have been developed for domestic species to one of conservation value. The study identifies several candidate genes for a specific disease, posthitis, which can now drive a more enlightened breeding effort to reduce its impact. While this study is a good demonstration of how tools developed for specific industries can be taken advantage of for environmental issues, it is unlikely that species without a closely related relative of economic impact can be studied in a similar manner.
Quality of written English: Acceptable
Author respond: We appreciate very much sharing our view on how new tools developed for domestic animals may be efficiently used for related species living in wild conditions.
Reviewer 2: Prof. Lev Klebanov
The manuscript contains an interesting study of finding SNP markers significantly associated with the incidence of posthitis and mine the genome for candidate genes potentially involved in the development of the disease. The methods developed here are based on the use of high-density bovine SNP microarray to investigate the entire bison genome. The study shows that among the 25 top significant markers, 18 were located on chromosome 15. Eleven significant markers were located on chromosome 15 in the close vicinity of olfactory receptor genes. In my opinion, the manuscript on the review is very interesting and worth to be published in Biology Direct. However, I have some questions, and think that the answers them may make the presentation better. Namely, on page 10 there is written: “SNPs with significant deviations from the Hardy–Weinberg equilibrium (HWE) at the significance level of P < 0.0001 were removed prior to association analysis. Additionally, the SNP selection criteria were applied based on a minor allele frequency (MAF) of at least MAF > 0.005, due to highly monomorphic data and the technical quality of the chip – a minimum call rate of 98%.” I did not understood, why such levels of P and MAF were used here.
Author respond: Illumina BovineHD BeadChip is designed especially for the domestic cattle and some out-group Bovidae species (bison, water buffalo, yak). The out-group species are only related with Bos species, therefore we have to be aware that bovine SNPs sometimes may generate spurious results when we use bison DNA. The only way to minimize that is to exclude all samples (bisons) giving low call rate. Moreover careful analysis of SNP cluster quality is also of key value. It was the reason why we used higher cal rate (>98%) than standard value (95%).
Another problem with bison SNPs is their limited number. The Lowland (Bialowieza) line of the European bison was restored from 7 pure Lowland animals – one male and six females. The relationship and the inbred level of all living animals are extremely high. Obviously, such low heterozygosity reduces the number of SNPs available for association study. Our intention was to keep almost each polymorphic SNP to cover the whole genome with markers located with relatively short distance from each other. Therefore we decreased the typical level of MAF from 0.01 to 0.005.
Summing up, both criteria were used to maximize the number of good-quality heterozygotic SNPs.
Analyzing criteria for SNP selection we found that one sentence should be corrected “SNPs in the control population with significant deviations from the Hardy–Weinberg equilibrium (HWE) at the significance level of P < 0.0001 were removed prior to association analysis”.
Genome-Wide Association Study
Identity By Descent
Bos Taurus Autosome
False Discovery Rate
Single Nucleotide Polymorphism
Minor Allele Frequency
Megabases (106 basepairs)
kilobases (103 basepairs)
University of Maryland Database
Bush WS, Moore JH. Chapter 11: genome-wide association studies. PLoS Comput Biol. 2012;8(12):e1002822. doi:10.1371.
Cheng R, Lim JE, Samocha KE, Sokoloff G, Abney M, Skol AD, et al. Genome-wide association studies and the problem of relatedness among advanced intercross lines and other highly recombinant populations. Genetics. 2010;185(3):1033–44.
Johnston SE, McEwan JC, Pickering NK, Kijas JW, Beraldi D, Pilkington JG, et al. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population. Mol Ecol. 2011;20:2555–66.
Zhang H, Wang Z, Wang S, Li H. Progress of genome wide association study in domestic animals. J Anim Sci Biotechnol. 2012;3:26.
Huang X, Han B. Sequencing-based GWAS for cloning QTLs in wild and cultivated rice. San Diego, USA: Plant and Animal Genome XXI; 2013.
Pertoldi C, Tokarska M, Wojcik JM, Demontis D, Loeschcke V, Gregersen VR, et al. Depauperate genetic variability detected in the American and European bison using genomic techniques. Biol Direct. 2009;4:48.
Tokarska M, Marshall T, Kowalczyk R, Wojcik JM, Pertoldi C, Kristensen TN, et al. Effectiveness of microsatellite and SNP markers for parentage and identity analysis in species with low genetic diversity: the case of European bison. Heredity. 2009;103:326–32.
Kamiński S, Olech W, Oleński K, Nowak Z, Ruść A. Single nucleotide polymorphisms between two lines of European bison (Bison bonasus) detected by the use of Illumina Bovine 50 K Bead Chip. Conserv Genet Resour. 2012;4:311–4.
Pucek Z, Belousova IP, Krasinski ZA, Krasinska M, Olech W. European bison status survey and conservation action plan IUCN/SSC bison specialist group. Gland: IUCN; 2004.
Olech W. The analysis of European bison genetic diversity using the pedigree data. In: Kita J, Anusz K, editors. Health threats for the European bison particularly in free-roaming populations in Poland. Warszawa: SGGW Publishers; 2006. p. 205–10.
Korochkina LN, Kochko FP. K voprosu o smertnosti Zubrom v estestvennych uslovijach Belovezhskojj Pushchi. Zapov Bel Issl. 1982;6:96–103 (in Russian).
Shabailo VE, Pererva VI. Issledovanie zabolevanij samcov zubrov Belovezhoskoj Pushchi i Nadvirnyanskogo lespromkhoza. In: Materialy nauchnoj konferencii posvashchennoj 50-letiyu issledovanij v Belovezhoskoj Pushchi; 1989 Dec 19–21; Kamenjuki. p. 219–20.
Wolf O, Jakob W, Lange A, Rudolph M, Borchers K, Glazel PS, et al. Balanoposthitis of the European bison (Bison bonasus). Proceedings of International Symposium “European bison – yesterday, today and tomorrow”, 2000 Dec 2–10; Siauliai, Lithuania; p. 68–73.
Krasinska M, Krasinski ZA. Course and dispersion of the posthitis/balanoposthitis in males of the European bison in the Polish part of Bialowieza Forest. Natl Park Nat Reserv. 2010;29(2):107–28 (in Polish, abstract in English).
Gill J. Zarys fizjologii żubra (Outlines of bison physiology). Warszawa: Severus Publishers; 1999 (in Polish).
Pertoldi C, Wójcik JM, Tokarska M, Kawałko A, Kristensen TN, Loeschcke V, et al. Genome variability in European and American bison detected using the BovineSNP50 BeadChip. Conserv Genet. 2010;11:627–34.
Raczyński J. European bison pedigree book. European Bison Conservation Newsletter. 2008;1:87–104.
Matuszewska M, Olech W, Bielecki W, Osinska B. The influence of inbreeding on the pathological changes occurrence in European bison male reproductive tract. Natl Park Nat Reserv. 2004;23:679–85 (in Polish, abstract in English).
Radwan J, Kawalko A, Wojcik JM, Babik W. MHC-DRB3 variation in a free-living population of the European bison, Bison bonasus. Mol Ecol. 2007;16(3):531–40.
Hirota T, Takahashi A, Kubo M, Tsunoda T, Tomita K, Sakashita M, et al. Genome-wide association study identifies eight new susceptibility loci for atopic dermatitis in the Japanese population. Nat Genet. 2012;44(11):1222–6.
Colas E, Pedrola N, Devis L, Ertekin T, Campoy I, Martínez E, et al. Molecular markers of endometrial carcinoma detected in uterine aspirates. Int J Cancer. 2011;129(10):2435–44.
von Thun A, Birtwistle M, Kalna G, Grindlay J, Strachan D, Kolch W, et al. ERK2 drives tumour cell migration in three-dimensional microenvironments by suppressing expression of Rab17 and liprin-β2. J Cell Sci. 2012;125(6):1465–77.
Wan B, Zhou YB, Zhang X, Zhu H, Huo K, Han ZG. hOLFML1, a novel secreted glycoprotein, enhances the proliferation of human cancer cell lines in vitro. FEBS Lett. 2008;582(21–22):3185–92.
Fedele M, Franco R, Salvatore G, Paronetto MP, Barbagallo F, Pero R, et al. PATZ1 gene has a critical role in the spermatogenesis and testicular tumours. J Pathol. 2008;215(1):39–47.
Liu Y, Helms C, Liao W, Zaba LC, Duan S, Gardner J, et al. A genome-wide association study of psoriasis and psoriatic arthritis identifies new disease loci. PLoS Genet. 2008;28;4(3):e1000041.
Schrodi SJ. Genome-wide association scan in psoriasis: new insights into chronic inflammatory disease. Expert Rev Clin Immunol. 2008;4(5):565–71.
Cole JB, Wiggans GR, Ma L, Sonstegard TS, Lawlor Jr TJ, Crooker BA, et al. Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U.S. Holstein cows. BMC Genomics. 2011;12:408.
Qanbari S, Pimentel ECG, Tetens J, Thaller G, Lichtner P, Sharifi AR, et al. The pattern of linkage disequilibrium in German Holstein cattle. Anim Genet. 2010;41:346–56.
Väisänen L, Has C, Franzke C, Hurskainen T, Tuomi ML, Bruckner-Tuderman L, et al. Molecular mechanisms of junctional epidermolysis bullosa: Col 15 domain mutations decrease the thermal stability of collagen XVII. J Invest Dermatol. 2005;125(6):1112–8.
Asaoka K, Ikeda K, Hishinuma T, Horie-Inoue K, Takeda S, Inoue S. A retrovirus restriction factor TRIM5alpha is transcriptionally regulated by interferons. Biochem Biophys Res Commun. 2005;338(4):1950–6.
Kim M, Ogawa M, Fujita Y, Yoshikawa Y, Nagai T, Koyama T, et al. Bacteria hijack integrin-linked kinase to stabilize focal adhesions and block cell detachment. Nature. 2009;459(7246):578–82.
Techanukul T, Sethuraman G, Zlotogorski A, Horev L, Macarov M, Trainer A, et al. Novel and recurrent FERMT1 gene mutations in Kindler syndrome. Acta Derm Venereol. 2011;91(3):267–70.
Kawahira K. Immunohistochemical staining of proliferating cell nuclear antigen (PCNA) in malignant and nonmalignant skin diseases. Arch Dermatol Res. 1999;291:413–8.
Plo I, Zhang Y, Le Couédic JP, Nakatake M, Boulet JM, Itaya M, et al. An activating mutation in the CSF3R gene induces a hereditary chronic neutrophilia. J Exp Med. 2009;206(8):1701–7.
Shiraishi H, Masuoka M, Ohta S, Suzuki S, Arima K, Taniguchi K, et al. Periostin contributes to the pathogenesis of atopic dermatitis by inducing TSLP production from keratinocytes. Allergol Int. 2012;61:563–72.
Masuoka M, Shiraishi H, Ohta S, Suzuki S, Arima K, Aoki S, et al. Periostin promotes chronic allergic inflammation in response to Th2 cytokines. J Clin Invest. 2012;122(7):2590–600.
Liu Y, Luo W, Chen S. Comparison of gene expression profiles reveals aberrant expression of FOXO1, Aurora A/B and EZH2 in lesional psoriatic skins. Mol Biol Rep. 2011;38:4219–24.
Krasinska M, Krasinski ZA. European bison. The nature monograph. Bialowieza: Mammal Research Institute Polish Academy of Science; 2007.
Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, et al. A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009;10(4):R42.
This work was financially supported by the University of Warmia and Mazury, grant no 105-0804. CP was supported by a grant from Danish Natural Science Research Council (grant numbers: 11-103926, 09-065999, 95095995), the Carlsberg Foundation (grant number 2011-01-0059) and the Aalborg Zoo Conservation Foundation (AZCF).
The authors declare that they have no competing interests.
KO performed SNP genotyping and done statistical analysis, MT provided the biological samples and pathological data, DMH cooperated with KO in genotyping and statistical analysis, PP participated in sequencing of single SNPs, AR prepared DNA samples, CP discussed the design of the experiment and co-edited the manuscript, SK designed the experiment and edit the manuscript. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
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 in a credit line to the data.
About this article
Cite this article
Oleński, K., Tokarska, M., Hering, D.M. et al. Genome-wide association study for posthitis in the free-living population of European bison (Bison bonasus). Biol Direct 10, 2 (2015). https://doi.org/10.1186/s13062-014-0033-6
- Illumina BovineHD 777K