Articles | Volume 65, issue 3
Arch. Anim. Breed., 65, 267–274, 2022
https://doi.org/10.5194/aab-65-267-2022
Arch. Anim. Breed., 65, 267–274, 2022
https://doi.org/10.5194/aab-65-267-2022
Original study
28 Jul 2022
Original study | 28 Jul 2022

Parentage testing and looking for single nucleotide markers associated with antler quality in deer (Cervus elaphus)

Parentage testing and looking for single nucleotide markers associated with antler quality in deer (Cervus elaphus)
Edith Elblinger1, Julianna Bokor2, Árpád Bokor1, Vilmos Altbäcker1, János Nagy2, József Szabó2, Bertalan Sárdi2, Adrian Valentin Bâlteanu3, Zsolt Rónai4, László Rózsa5, József Rátky6, István Anton5, and Attila Zsolnai5,7 Edith Elblinger et al.
  • 1Kaposvár Campus, Hungarian University of Agriculture and Life Sciences, Kaposvár, 7400, Hungary
  • 2Game Management Landscape Center, Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, Bőszénfa, 7475, Hungary
  • 3Institute of Life Sciences, University of Agricultural Sciences and Veterinary Medicine, Cluj-Napoca, Romania
  • 4Department of Medical Chemistry, Molecular Biology and Pathobiochemistry, Eötvös Loránd University, Budapest, 1053, Hungary
  • 5Kaposvár Campus, Hungarian University of Agriculture and Life Sciences, Herceghalom, 2053, Hungary
  • 6Department of Obstetrics and Food Animal Medicine Clinic, University of Veterinary Medicine Budapest, Budapest, 1078, Hungary
  • 7Institute for Farm Animal Gene Conservation, National Centre for Biodiversity and Gene Conservation, Gödöllő, 2100, Hungary

Correspondence: Attila Zsolnai (attila.zsolnai@gmail.com)

Abstract

To provide a cost-efficient parentage testing kit for red deer (Cervus elaphus), a 63 SNP set has been developed from a high-density Illumina BovineHD BeadChip containing 777 962 SNPs after filtering of genotypes of 50 stags. The successful genotyping rate was 38.6 % on the chip. The ratio of polymorphic loci among effectively genotyped loci was 6.5 %. The selected 63 SNPs have been applied to 960 animals to perform parentage control. Thirty SNPs out of the 63 had worked on the OpenArray platform. Their combined value of the probability of identity and exclusion probability was 4.9×10-11 and 0.99803, respectively.

A search for loci linked with antler quality was also performed on the genotypes of the above-mentioned stags. Association studies revealed 14 SNPs associated with antler quality, where low-quality antlers with short and thin main beam antlers had values from 1 to 2, while high-quality antlers with long and strong main beams had values between 4 and 5. The chance for a stag to be correctly identified as having high-value antlers is expected to be over 88 %.

1 Introduction

Molecular genetic information was already used successfully in different species, e.g. for the determination of the origin of modern cattle (Beja-Pereira et al., 2006), for analysis of genetic diversity of horses (Aberle et al., 2004), and classification purposes in pigs (Zsolnai et al., 2006). Availability of SNP chip technology triggered its application in population genetics (Lancioni et al., 2016; Zsolnai et al., 2020a) and in association studies (Zsolnai et al., 2020b).

In deer, comprehensive studies have been conducted in search of quantitative trait loci (Slate et al., 2002) and for the determination of the inbreeding effect on breeding success (Slate et al., 2000). Molecular data have helped to confirm the genetic integrity of Carpathian red deer (Feulner et al., 2004) and were used to seek key genes in antler development (Molnár et al., 2007). The mechanism of antler development is also described by Gyurján et al. (2007) and Stéger et al. (2010).

One aspect of methodologies is the application of a paternity test, which is already routinely implemented for example in cattle and sheep (Glowatzki-Mullis et al., 1995, 2007) in many countries. There is an obvious demand from wildlife management and deer breeders to follow up the lineage, especially when important traits (antler characteristics, body weight, etc.) are logged in the herd book. The molecular genetic approach has already been implemented for completing this task previously. Bonnet et al. (2002) applied 11 microsatellites in three multiplex reactions, and Haanes et al. (2005) used 25 loci in six reactions for genotyping in deer.

In Hungary, two groups have developed multiplex microsatellite tests for parentage testing. Zsolnai et al. (2009) used nine loci in one reaction, and Szabolcsi et al. (2014) applied 10 loci in two reactions. Their obtained probability of identity was 1×10-11 and 1×10-12, respectively. Since deer is an attractive game animal and research methodologies are becoming more advanced, studies – including fitness variation, performance prediction, or comparative population analysis – have been performed using Bovine BeadChip on deer (Kasarda et al., 2014). A total of 136 SNPs selected from the commercially unavailable Illumina 50K CervusSNP50 chip have already been developed for the New Zealand deer industry (Rowe et al., 2015). Such chip development or SNP discovery is on its way in Hungary and can be based on a genome map of red deer, CerEla1.0 (Bana et al., 2018).

Here, we aimed to select and test markers suitable for parentage testing in red deer based on Illumina Bovine HD SNP chip and OpenArray platforms to reduce the costs and labour of microsatellite genotyping. We also aimed to perform a genome-wide association study, to look for loci associated with antler quality, which can contribute to a more effective breeding strategy.

The selection of the Illumina Bovine SNP chip was based on its commercial availability and the high similarity of the bovine and deer genomes (Bana et al., 2018).

2 Materials and methods

We collected 49 red deer blood samples of stags from 10 places, including Baltacím (461446′′ N, 174932′′ E), Cserhát (461400′′ N, 174820′′ E), Égeres (461434′′ N, 174859′′ E), Homokos gödör (461420′′ N, 174906′′ E), Koronafürt felső (461437′′ N, 174947′′ E), Kuszó lucerna jobb alsó (461411′′ N, 174824′′ E), Kuszó lucerna jobb felső (461421′′ N, 174817′′ E), Szent-Tamás (461419′′ N, 174850′′ E), Templom Dél Kelet (461417′′ N, 174841′′ E), and Templom Dél Nyugat (461414′′ N, 174833′′ E) to perform SNP chip genotyping on the Illumina BovineHD BeadChip containing 777 962 SNPs. Sampling on different sites offers elimination of the formation of subgroups due to family structure, and it decreases the chance of population stratification. Selection criteria for stags were to represent high or low antler quality values within each sampled area to search for and elucidate the difference in antler quality. Sampling was performed by trained veterinarians during routine sample collection organized independently from this study. Sampling followed standard procedures and relevant national guidelines to ensure appropriate animal care.

Antler quality values were determined by a scoring method applied in the Bőszénfa deer farm, Hungary. Low-quality antlers (short and thin main beam antlers) reached a value from 1 to 2, while high-quality antlers (long and strong main beam) had values between 4 and 5. The low- and high-quality groups consisted of 24 and 25 animals, respectively. The preferred phenotype for the antler is to have long beams, high mass, and symmetry. Samples were prepared and genotyped by Neogen Europe Ltd. (Scotland, UK).

For parentage testing, the hair roots of 960 animals were collected. DNA from hair root samples was extracted by the Chelex (Bio-Rad, USA) procedure; 8–12 hair root samples were immersed into a 100 µL, 5 % Chelex mixture, incubated overnight at 56 C, and treated 10 min at 96 C. The 63 SNPs, selected from BovineHD SNP chip, were genotyped on the OpenArray (Life Technologies, USA) platform. The surrounding sequences of the SNPs used for the design of primers and probes are available from https://webdata.illumina.com/downloads/productfiles/bovinehd/bovinehd-b1-annotation-file.zip (last access: 22 July 2022) and available in Supplement Table S1.

After filtering (part a) 777 962 loci – presented on the Illumina BovineHD BeadChip – for a call rate >0.95, there were two alleles with MAF >0.05. Additional filtering (part b) was included in the case of selection of SNPs suitable for parentage testing: Hardy–Weinberg equilibrium p value higher than 0.001, and MAF had to be over 0.3. For the association study, only part a of filtering was applied. The sequences of the loci found in GWAS are in Supplement Table S2.

Statistical analyses were performed by SVS SNP & Variation Suite 8.8.1. software (Golden Helix, Bozeman, MT, USA).

The performance of the SNPs selected for parentage testing was characterized by GenAlEx 6.5 (Peakall et al., 2012). For correction of population structure, the genomic kinship matrix was used in a multi-locus mixed model (Segura et al., 2012). The model was

(1) y = X β + Z u + e ,

where y is the antler quality score, X is the matrix of fixed effects composed of SNPs and covariates (date, place of birth, and the father), Z is the matrix of random animal effects, e means the residual effects, and β and u are vectors representing coefficients of fixed and random effects, respectively.

SVS and the PLINK software v.1.9 (Purcell et al., 2007) were used to build a multidimensional scaling (MDS) plot using a genome-wide identity-by-state pairwise distance matrix (mds-plot 2 and cluster options). Population stratification was calculated by SVS as described by Price et al. (2006). Python 3.6 programming language (Van Rossum and Drake, 2009) and the matplotlib 3.2.1 library (Hunter, 2007) were used for the visualization of PCA data.

Classification procedures were performed by Geneclass2 software (Piry et al., 2004).

3 Results

To achieve the goal to select SNPs suitable for pedigree control, we tested 777 962 SNPs on 50 red deer stags. After filtering for call rate and minor allele frequency, the number of retained loci was 103 562 and 25 919, respectively (filtering part a). After additional filtering (filtering part b) for LD pruning and Hardy–Weinberg equilibrium (7146 and 583 loci, respectively), the SNPs proved to be adequate for statistical analyses.

The genotypic principal component analysis of both 25 919 and 583 SNPs did not reveal population stratification according to antler quality (Fig. 1a, b, 583 loci). Eigenvalues of axes 1, 2, and 3 were 5.422, 3.116, and 2.339, respectively.

https://aab.copernicus.org/articles/65/267/2022/aab-65-267-2022-f01

Figure 1Principal component analysis of 49 stags based on their 583 loci. The blue circle shows animals with low antler quality, and the red triangle shows animals with high antler quality. Eigenvalues of axes C1 and C2 are 5.422 and 3.116, respectively.

Download

We selected 63 SNPs for parentage testing. The probability level for the identity of this set was 7×10-24. The surrounding sequence of selected SNPs was used to design primers and probes by OpenArray Product Configurator (Life Technologies, USA) for genotyping 960 samples on the OpenArray platform (Life Technologies, USA). The callable 30 SNPs (Supplement Table S1) out of the selected 63 SNPs reached the 4.9×10-11 value of the probability of identity on 960 animals.

In GWAS analysis date and place of birth and the father of the animal were included in the model. A genomic kinship matrix was used to correct possible relatedness, and 14 SNPs (Supplement Table S2) out of 25 919 displayed significant differentiation between the two groups (high and low antler quality). Eigenvalues of PCA coordinates were 9.857, 5.637, and 5.004 for axes 1, 2, and 3. The selected SNPs were able to position high- and low-antler-value animals into two distinct groups (Fig. 2a, b).

https://aab.copernicus.org/articles/65/267/2022/aab-65-267-2022-f02

Figure 2Principal component analysis of 49 stags based on their 14 loci selected to differentiate by antler quality. The blue circle shows animals with low antler quality, and the red triangle shows animals with high antler quality. Eigenvalues of axes C1 and C2 are 9.857 and 5.637, respectively.

Download

Depending on the algorithm chosen within GeneClass2 software (Piry et al., 2004) the number of misclassified animals (from the low-value group) ranged from 3 to 6, which is 6 %–12 % of the total number of individuals.

4 Discussion

Haynes and Latch (2012) have achieved a 38.7 % successful genotyping rate on Odocoileus hemionus and Odocoileus virginianus using Illumina Bovine BeadChip and they found 5 % polymorphic loci among the successfully genotyped SNPs on the Bovine BeadChip. In our case on Cervus elaphus, 13.3 % of the genotyped loci were effective, while the ratio of polymorphic loci among effectively genotyped loci was 6.5 %. The selected 63 SNPs for parentage testing had an identity value (7×10-24) that was similar to the value of SNP sets used in white pigs (Rohrer et al., 2007) or in Mangalitza pigs (Zsolnai et al., 2013) (1×10-23). Genotyping 960 animals on the OpenArray platform, only 30 SNPs out of 63 were callable. The 4.9×10-11 value of the probability of identity of these 30 SNPs is similar to the previous 1×10-11 (Zsolnai et al., 2009) or 1×10-12 values (Szabolcsi et al., 2014) of deer microsatellite sets. The exclusion probability of the reported 30 SNPs is 0.99803. As for microsatellite sets, this value is slightly better than the previously reported 0.99630 (Zsolnai et al., 2009) and lower than 0.99999 (Szabolcsi et al., 2014). The price of typing 63 SNPs was about a third of the price of microsatellite genotyping and required less working time. By switching from OpenArray (Life Technologies, USA) to another genotyping platform, the number of callable SNPs and the exclusion probability value could be increased to 0.99999.

For the genome-wide association study, we used data acquired from 50 animals where antler quality was known. The genotypic principal component analysis of 583 SNPs did not divide the population by antler quality (Fig. 1a, b). The lambda value calculated by SVS was 1.02, indicating that there was no hidden population stratification among the sampled animals.

After filtering (part a) we found 14 SNPs that are associated with antler values (Supplement Table S2). Previously such a GWA-aided search was successfully applied in the case of Mangalitza pigs to find trait-associated loci for checking quality and quantity parameters in the meat industry (Szántó-Egész et al., 2013, 2016).

Evaluating our 14 SNPs on red deer by PCA (Fig. 2) and classification algorithms implemented in GeneClass2 software (Piry et al., 2004), we have found that animals with a high antler value have never fallen into the low-value group. However, several animals (4 %–10 %) from the low-value group have been assigned to the high-value group.

Several candidate genes have been found in the close vicinity of the antler-quality-associated markers (Supplement Table S2). Among the candidate genes, we are mentioning those where bone- or antler-related functions are known. Sentrin-specific protease 1 (SENP1) markedly enhances androgen-receptor-mediated transcription in males, mediates cell growth and differentiation, and maintains male reproductive functions (Cheng et al., 2004). Collagen alpha-1(II) chain (COL2A1) has key roles in chondrogenesis and osteogenesis (Jia et al., 2021); it is upregulated in the main beams of antlers (Yao et al., 2020). Vitamin D receptor (VDR) alleles are associated with postmenopausal osteoporosis susceptibility and bone mineral density (Zhang et al., 2018). Dynamin 1 like (DNM1L) belongs to a superfamily of GTPases, which are related to osteoclast differentiation and bone loss (Jeong et al., 2021). Tyrosyl-TRNA synthetase 2 (YARS2) is coupled with the impaired ability of the bone marrow to produce normal red blood cells (Riley et al., 2013). Plakophilin 2 (PKP2) overexpression was able to stop the proliferation of osteosarcoma (He et al., 2021). Zinc finger protein 518B (ZNF518B) is a putative ageing modulator (Sleiman et al., 2020). WD repeat-containing protein 1 (WDR1) is part of a deer antler extract (Yao et al., 2019). Solute carrier family 2 member 9 (SLC2A9) plays a role in gout, a form of rheumatic arthritis (Merriman and Dalbeth, 2011), and has a significant role in glucose homeostasis. Among the enriched targets of miRNAs there is zinc finger and BTB domain containing 49 (ZBTB49) after a growth factor-beta 1 stimulation (Ong et al., 2017). The latter compound has a high impact on bone formation and resorption (Bonewald et al., 1990). Molybdenum cofactor sulfurase (MOCOS) is involved in purine metabolism (Kurzawski et al., 2012), and extracellular purine is vital for bone homeostasis (Agrawal and Jørgensen, 2021). Elongation protein 2 (ELP2) has a role in osteogenesis and osteogenic differentiation (Wu et al., 2021). Exocyst complex component 6B (EXOC6B) is associated with joint dislocation syndrome (Girisha et al., 2016). Nucleoporin 58 (NUP58) mediates molecular trafficking to and from the nucleus and can be tied to delayed abscission during mitosis (Hartono et al., 2019). Spermatogenesis associated 13 (SPATA13) can weaken the migration of type I collagen via activation of GTPase (Jean et al., 2013). Zinc finger homeobox 4 (ZFHX4) orchestrates endochondral bone formation (Nakamura et al., 2021). Lysine demethylase 4C (KDM4C) regulates condensin-dependent heterochromatin reorganization and is connected with deterioration and premature bone ageing (Huang et al., 2019). Glycine decarboxylase (GLDC) locus is associated with overall survival in patients with osteosarcoma (Lin et al., 2020). Glutathione peroxidase (GPX) level is higher at the upper sections of the velvet antler and associated with lower Ca content (Cheng et al., 2017). SMAD family member 7 (SMAD7) directly interacts with microtubule actin crosslinking factor 1, initiating downstream osteogenic pathways (Zhao et al., 2020). Dymeclin (DYM) is an important part of the extracellular matrix and has a role in bone development (Denais et al., 2011). Among the differentially abundant metabolites of ageing bone marrow mesenchymal stem cells, acetyl-CoA acyltransferase 2 (ACAA2) is related to senescence and lipid metabolism (Yu et al., 2022).

The approach for finding selection markers for antler quality can be enhanced by (i) incorporating more animals into the study, (ii) development of a deer SNP chip based on red deer populations existing in Hungary and on existing whole-genome sequencing (Bana et al., 2018), and (iii) collecting blood at different stages of antler development and merging genotypic, phenotypic, and blood-metabolite data for deeper analysis.

As for parentage testing, the identified 30 SNPs produce similar discriminating power to that of a microsatellite set (Zsolnai et al., 2009). The advantage of SNP typing lays in its cost-effectiveness compared to microsatellites.

As for antler quality predictions, the listed 14 SNPs coupled with antler quality can be used to identify those animals at birth that are supposed to develop good-quality antlers in their life. The chance that an animal (preselected by SNP typing) has a desirable high-valued antler is expected to be over 88 %.

Data availability

The surrounding sequences of the SNPs described in the MS are available in the Supplement.

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/aab-65-267-2022-supplement.

Author contributions

Conceptualization was formulated by JB, ÁB, IA, and AZ. Resources were ensured by ÁB, VA, JN, JS, BS, LR, JR, and IA. Data were curated by EE, JB, ZR, and AZ. Analyses were performed by EE, AVB, IA, and AZ. The first draft of the manuscript was written by EE, VA, IA, and AZ. All authors commented on and approved the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Ethical statement

Sampling was performed by trained veterinarians during routine sample collection organized independently from this study. Sampling followed standard procedures and relevant national guidelines to ensure appropriate animal care.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

We are grateful to our reviewers for providing their valuable notes, comments, and questions. Many thanks to Copernicus Publications' staff members for their efforts to correct and improve the article.

Financial support

This research has been supported by the Emberi Eroforrások Minisztériuma (grant no. TÁMOP-4.2.2.A11/1/KONV-2012-0039).

Review statement

This paper was edited by Antke-Elsabe Freifrau von Tiele-Winckler and reviewed by Gabriela M. Iglesias and one anonymous referee.

References

Aberle, K. S., Hamann, H., Drögemüller, C., and Distl, O.: Genetic diversity in German draught horse breeds compared with a group of primitive, riding and wild horses by means of microsatellite DNA markers, Anim. Genet., 35, 270–277, https://doi.org/10.1111/j.1365-2052.2004.01166.x, 2004. 

Agrawal, A. and Jørgensen, N. R.: Extracellular purines and bone homeostasis, Biochem. Pharmacol., 187, 114425, https://doi.org/10.1016/j.bcp.2021.114425, 2021. 

Bana, Á. N., Nyíri, A., Nagy, J., Frank, K., Nagy, T., Stéger, V., Schiller, M., Lakatos, P., Sugár, L., Horn, P. Barta, E., and Orosz, L.: The red deer Cervus elaphus genome CerEla1.0: sequencing, annotating, genes, and chromosomes, Mol. Genet. Genomics, 293, 665–684, https://doi.org/10.1007/s00438-017-1412-3, 2018. 

Beja-Pereira, A., Caramelli, D., Lalueza-Fox, C., Vernesif, C., Ferranda, N., Casolih, A., Goyachei, F., Royoi, I.J., Contid, S., Larid, M., Martinij, A., Ouraghk, L., Magidl, A., Atashl, A., Zsolnai, A., Boscaton, P., Triantaphylidiso, C., Ploumip, K., Sineoq, L., Mallegnir, F., Taberlet, P., and Erhardts, G.: The origin of European cattle: Evidence from modern and ancient DNA, P. Natl. Acad. Sci. USA, 103, 8113–8118, https://doi.org/10.1073/pnas.0509210103, 2006. 

Bonewald, L. F. and Mundy, G. R.: Role of transforming growth factor-beta in bone remodeling, Clin. Orthop. Relat. R., 250, 261–276, 1990. 

Bonnet, A., Thévenon, S., Maudet F., and Maillard, J. C.: Efficiency of semi-automated fluorescent multiplex PCRs with 11 microsatellite markers for genetic studies of deer populations, Anim. Genet., 33, 343–350, https://doi.org/10.1111/j.1365-2052.2006.01550.x, 2002. 

Cheng J., Wang, D., Wang, Z., and Yeh, E. T. H.: SENP1 enhances androgen receptor-dependent transcription through desumoylation of histone deacetylase 1, Mol. Cell. Biol., 24, 6021–6028, https://doi.org/10.1128/MCB.24.13.6021-6028.2004, 2004. 

Cheng, S. L., Jian, Y. L., Chen, C. M., and Liu, B. T.: Relationships between antioxidants and quality characteristics from velvet antlers of formosan sambar deer, Korean J. Food Sci. An., 37, 542–551, https://doi.org/10.5851/kosfa.2017.37.4.542, 2017. 

Denais, C., Dent, C. L., Southgate, L., Hoyle, J., Dafou, D., Trembath, R. C., and Machado, R. D.: Dymeclin, the gene underlying Dyggve-Melchior-Clausen syndrome, encodes a protein integral to extracellular matrix and golgi organization and is associated with protein secretion pathways critical in bone development, Human Mutation, 32, 231–239, https://doi.org/10.1002/humu.21413, 2011. 

Feulner, P. G. D., Bielfeldt, W., Zachos, F. E., Bradvarovic, J., Eckert, I., and Hartl, G. B.: Mitochondrial DNA and microsatellite analyses of the genetic status of the presumed subspecies Cervus elaphus montanus (Carpathian red deer), Heredity, 93, 299–306, https://doi.org/10.1038/sj.hdy.6800504, 2004. 

Girisha, K. M., Kortüm, F., Shah, H., Alawi, M., Dalal, A., Bhavani, G. S. L., and Kutsche, K.: A novel multiple joint dislocation syndrome associated with a homozygous nonsense variant in the EXOC6B gene, Eur. J. Hum. Genet., 24, 1206–1210, https://doi.org/10.1038/ejhg.2015.261, 2016. 

Glowatzki-Mullis, M. L., Gaillard, C., Wigger, G., and Fries, R.: Microsatellite-based parentage control in cattle, Anim. Genet., 26, 7–12, https://doi.org/10.1111/j.1365-2052.1995.tb02612.x, 1995. 

Glowatzki-Mullis, M. L., Muntwyler, J., and Gaillard, C.: Cost-effective parentage verification with 17-plex PCR for goats and 19-plex PCR for sheep, Anim. Genet., 38, 86–88, 2007. 

Gyurján, I., Molnár, A., Borsy, A., Stéger, V., Hackler, L., Zomborszky, Z., Papp, P., Duda, E., Deák, F., Lakatos, P., Puskás, L. G., and Orosz, L.: Gene expression dynamics in deer antler: mesenchymal differentiation toward chondrogenesis, Mol. Genet. Genomics, 277, 221–235, https://doi.org/10.1007/s00438-006-0190-0, 2007. 

Haanes, H., Rosef, O., Veiberg, V., and Røed, K. H.: Microsatellites with variation and heredity applicable to genetic studies of Norwegian red deer (Cervus elaphus atlanticus), Anim. Genet., 36, 454–455, https://doi.org/10.1111/j.1365-2052.2005.01343.x, 2005. 

Hartono, Hazawa, M., Lim, K. S., Dewi, F. R. P., Kobayashi, A., and Wong R. W.: Nucleoporin Nup58 localizes to centrosomes and mid-bodies during mitosis, Cell Division, 14, 7, https://doi.org/10.1186/s13008-019-0050-z, 2019. 

Haynes, G. D. and Latch, E. K.: Identification of Novel Single Nucleotide Polymorphisms (SNPs) in Deer (Odocoileus spp.) Using the BovineSNP50 BeadChip, PLoS ONE, 7, e36536, https://doi.org/10.1371/journal.pone.0036536, 2012. 

He, X. L., Xiao, Q., Zhou, Z. P., and Hui, C. Y.: PPM1D accelerates proliferation and metastasis of osteosarcoma by activating PKP2, Eur. Rev. Med. Pharmaco., 25, 78–85, https://doi.org/10.26355/eurrev_202101_24351, 2021. 

Huang, B., Wang, B., Lee, W. Y. W., Pong, K., Leung, K. T., Li, X., Liu, Z., Chen R., Lin, J. C., Tsang, L. L., Liu, B., Ruan, Y. C., Chan, H. C., Li, G., and Jiang, X.: KDM3A and KDM4C Regulate Mesenchymal Stromal Cell Senescence and Bone Aging via Condensin-mediated Heterochromatin Reorganization, iScience, 21, 375–390, https://doi.org/10.1016/j.isci.2019.10.041, 2019. 

Hunter, J. D.: Matplotlib: A 2D Graphics Environment, Comput. Sci. Eng., 9, 90–95, 2007. 

Jean, L., Majumdar, D., Shi, M., Hinkle, L. E., Diggins, N. L., Ao, M., Broussard, J. A., Evans, J. C., Choma, D. P., and Webb, D. J.: Activation of Rac by Asef2 promotes myosin II-dependent contractility to inhibit cell migration on type I collagen, J. Cell. Sci., 126, 5585–5597, https://doi.org/10.1242/jcs.131060, 2013. 

Jeong, S., Seong, J. H., Kang, J. H., Lee, D. S., and Yim, M.: Dynamin-related protein 1 positively regulates osteoclast differentiation and bone loss, FEBS Letters, 595, 58–67, https://doi.org/10.1002/1873-3468.13963, 2021. 

Jia, B., Zhang, L., Zhang, Y., Ge, C., Yang, F., Du, R., and Ba, H.: Integrated analysis of miRNA and mRNA transcriptomic reveals antler growth regulatory network, Mol. Genet. Genomics, 296, 689–703, https://doi.org/10.1007/s00438-021-01776-z, 2021. 

Kasarda, R., Moravčíková, N., and Trakovická, A.: Advances in genomic sequencing using Bovine SNP BeadChip in Deer, Acta Fytotechnica et Zootechnica, 17, 65–71, https://doi.org/10.15414/afz.2014.17.02.65-71, 2014. 

Kurzawski, M., Dziewanowski, K., Safranow, K., and Drozdzik, M.: Polymorphism of genes involved in purine metabolism (XDH, AOX1, MOCOS) in kidney transplant recipients receiving azathioprine, Therapeutic Drug Monitoring, 34, 266–274, https://doi.org/10.1097/FTD.0b013e31824aa681, 2012. 

Lancioni, H., Di Lorenzo, P., Cardinali, I., Ceccobelli, S., Capodiferro, M. R., Fichera, A., Grugni, V., Semino, O., Ferretti, L., Gruppetta, A., Attard, G., Achilli, A., and Lasagna, E.: Survey of uniparental genetic markers in the Maltese cattle breed reveals a significant founder effect but does not indicate local domestication, Anim. Genet., 47, 267–269, https://doi.org/10.1111/age.12408, 2016. 

Lin, O., Han, J., Sun, Q., Wen, L., and Wang, S.: Functional variant of IL33 is associated with survival of osteosarcoma patients, J. Bone Oncol., 2, 100270, https://doi.org/10.1016/j.jbo.2019.100270, 2020. 

Merriman, T. R. and Dalbeth, N.: The genetic basis of hyperuricaemia and gout, Joint Bone Spine, 78, 35–40, https://doi.org/10.1016/j.jbspin.2010.02.027, 2011. 

Molnár, A., Gyurján, I., Korpos, E., Borsy, A., Stéger, V., Buzás, Z., Kiss, I., Zomborszky, Z., Papp, P., Deák, F., and Orosz, L.: Identification of differentially expressed genes in the developing antler of red deer Cervus elaphus, Mol. Genet. Genomics, 277, 237–248, https://doi.org/10.1007/s00438-006-0193-x, 2007. 

Nakamura, E., Hata, K., Takahata, Y., Kurosaka, H., Abe, M., Abe, T., Kihara, M., Komori, T., Kobayashi, S., Murakami, T., Inubushi, T., Yamashiro, T., Yamamoto, S., Akiyama, H., Kawaguchi, M., Sakata, N., and Nishimura, R.: Zfhx4 regulates endochondral ossification as the transcriptional platform of Osterix in mice, Commun. Biol., 4, 1258, https://doi.org/10.1038/s42003-021-02793-9, 2021. 

Ong, J., Timens, W., Rajendran, V., Algra, A., Spira, A., Lenburg, M. E., Campbell, J. D., van den Berge, M., Postma, D. S., van den Berg, A., Kluiver, J., and Brandsma, C. A.: Identification of transforming growth factor-beta-regulated microRNAs and the microRNA-targetomes in primary lung fibroblasts, PLoS One, 12, e0183815, https://doi.org/10.1371/journal.pone.0183815, 2017. 

Peakall, R. and Smouse, P. E.: GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update, Bioinformatics, 28, 2537–2539, https://doi.org/10.1093/bioinformatics/bts460, 2012. 

Piry, S., Alapetite, A., Cornuet, J. M., Paetkau, D., Baudouin, L., and Estoup, A.: GeneClass2: A Software for Genetic Assignment and First-Generation Migrant Detection, J. Hered., 95, 536–539, https://doi.org/10.1093/jhered/esh074, 2004. 

Price, A. L., Patterson, N. J., Plenge, R. M., Weinblatt, M. E., Shadick, N. A., and Reich, D.: Principal Components Analysis corrects for stratification in Genome-Wide Association Studies, Nature Genet., 38, 904–909, https://doi.org/10.1038/ng1847, 2006. 

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., Maller, J., Sklar, P., de Bakker, P. I. W., Daly, M. J., and Sham, P. C.: PLINK: a tool set for whole-genome association and population-based linkage analyses, Am. J. Human Genet., 81, 559–575, https://doi.org/10.1086/519795, 2007. 

Riley, L. G., Menezes, M. J., Rudinger-Thirion, J., Duff, R., de Lonlay, P., Rotig, A., Tchan, M. C., Davis, M., Cooper, S. T., and Christodoulou, J.: Phenotypic variability and identification of novel YARS2 mutations in YARS2 mitochondrial myopathy, lactic acidosis and sideroblastic anaemia, Orphanet J. Rare Dis., 8, 193, https://doi.org/10.1186/1750-1172-8-193, 2013. 

Rohrer, G. A., Freking, B. A., and Nonneman, D.: Single nucleotide polymorphisms for pig identification and parentage exclusion, Anim. Genet., 38, 253–258, https://doi.org/10.1111/j.1365-2052.2007.01593.x, 2007. 

Rowe, S. J., Clarke, S. M., Van Stijn, T. C., Hyndman, D. L., Ward, J. F., McEwan, K. M., Dodds, K. G., McEwan, J. C., Newman, A. N., and Asher, G. W.: Developing genomic tools in the New Zealand, Deer Industry Proceedings of the New Zealand Society of Animal Production, 75, 91–93, 2015. 

Segura, V., Vihjalmsson, B. J., Platt, A., Korte, A., Seren, Ü., Long, Q., and Nordborg, M.: An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations, Nature Genet., 44, 825–830, https://doi.org/10.1038/ng.2314, 2012. 

Slate, J., Kruuk, L. E. B., Marshall, C., Pemberton, J. M., and Clutton-Brock, T. H.: Inbreeding depression influences lifetime breeding success in a wild population of red deer (Cervus elaphus), Proceedings of the Royal Society Biological Sciences, 267, 1657–1662, https://doi.org/10.1098/rspb.2000.1192, 2000. 

Slate, J., Visscher, P. M., Macgregor, S., Stevens, D., Tate, M. L., and Pemberton, J. M.: A Genome Scan for Quantitative Trait Loci in a Wild Population of Red Deer (Cervus elaphus), Genetics, 162, 1863–1873, 2002. 

Sleiman, M. B., Jha, P., Houtkooper, R., Williams, R. W., Wang, X., and Auwerx, J.: The Gene-Regulatory Footprint of Aging Highlights Conserved Central Regulators, Cell Rep., 32, 108203, https://doi.org/10.1016/j.celrep.2020.108203, 2020. 

Stéger, V., Molnár, A., Borsy, A., Gyurján, I., Szabolcsi, Z., Dancs, G., Molnár, J., Papp, P., Nagy, J., Puskás, L., Barta, E., Zomborszky, Z., Horn, P., Podani, J., Semsey, S., Lakatos, P., and Orosz, L.: Antler development and coupled osteoporosis in the skeleton of red deer Cervus elaphus: expression dynamics for regulatory and effector genes, Mol. Genet. Genomics, 284, 273–287, https://doi.org/10.1007/s00438-010-0565-0, 2010. 

Szabolcsi, Z., Egyed, B., Zenke, P., Pádár, Z., Borsy, A., Stéger, V., Pásztor, E., Csányi, S., Buzás, Z., and Orosz, L.: Constructing STR multiplexes for individual identification of Hungarian red deer, J. Forensic Sci., 59, 1090–1099, https://doi.org/10.1111/1556-4029.12403, 2014. 

Szántó-Egész, R., Ujhelyi, G., Koppányné Szabó, E., Jánosi, A., Mohr, A., Sipos, R., Dallmann, K., Micsinai, A., Zsolnai, A., Egerszegi, I., Anton, I., Tóth, G., Molnár, J., Stéger, V., Marincs, F., Tóth, P., and Rátky, J.: Mangalica termékek kimutatására alkalmas real-time PCR módszer fejlesztése (Development of a realtime PCR method for the detection of Mangalitza products), Élelmiszer Tudomány és Technológia, 1, 14–21, 2013. 

Szántó-Egész, R., Jánosi, A., Mohr, A., Szalai, G., Koppányné Szabó, E., Micsinai, A., Sipos, R., Rátky, J., Anton, I., and Zsolnai, A.: Breed-Specific Detection of Mangalica Meat in Food Products, Food Anal. Meth., 9, 889–894, https://doi.org/10.1007/s12161-015-0261-0, 2016. 

Van Rossum, G. and Drake, F. L.: Python 3 Reference Manual, Scotts Valley, CA, Publisher: CreateSpace, ISBN 978-1-4414-1269-0, 2009. 

Wu, W. J., Xia, C. L., Ou, S. J., Yang, Y., Ma, Y. F., Hou, Y. L., Yang, Q. P., Zhang, J., Li, J. W., Qi, Y., and Xu, C. P.: Novel elongator protein 2 inhibitors mitigating tumor necrosis factor-α induced osteogenic differentiation inhibition, Biomed Res. Int., 2021, 3664564, https://doi.org/10.1155/2021/3664564, 2021. 

Yao, B., Gao, H., Liu, J., Zhang, M., Leng, X., and Zhao, D.: Identification of potential therapeutic targets of deer antler extract on bone regulation based on serum proteomic analysis, Mol. Biol. Rep., 46, 4861–4872, https://doi.org/10.1007/s11033-019-04934-0, 2019. 

Yao, B., Wang, C., Zhou, Z., Zhang, M., Zhao, D., Bai, X., and Leng, X.: Comparative transcriptome analysis of the main beam and brow tine of sika deer antler provides insights into the molecular control of rapid antler growth, Cell. Mol. Biol. Lett., 25, 42, https://doi.org/10.1186/s11658-020-00234-9, 2020. 

Yu, X., Sun H., Gao, X., Zhang, C., Sun, Y., Wang, H., Zhang, H., Shi, Y., and He, X.: A comprehensive analysis of age-related metabolomics and transcriptomics reveals metabolic alterations in rat bone marrow mesenchymal stem cells, Aging, 14, 1014–1032, https://doi.org/10.18632/aging.203857, 2022. 

Zhang, L., Yin, X., Wang, J., Xu, D., Wang, Y., Yang, J., Tao, Y., Zhang, S., Feng, X., and Yan, C.: Associations between VDR Gene Polymorphisms and Osteoporosis Risk and Bone Mineral Density in Postmenopausal Women: A systematic review and Meta-Analysis, Sci. Rep., 8, 981, https://doi.org/10.1038/s41598-017-18670-7, 2018. 

Zhao, F., Ma, X., Qiu, W., Wang, P., Zhang, R., Chen, Z., Su, P., Zhang, Y., Li, D., Ma, J., Yang, C., Chen, L., Yin, C., Tian, Y., Hu, L., Li, Y., Zhang, G., Wu, X., and Qian, A.: Mesenchymal MACF1 facilitates SMAD7 nuclear translocation to drive bone formation, Cells, 9, 616, https://doi.org/10.3390/cells9030616, 2020. 

Zsolnai, A., Radnóczy, L., Fésüs, L., and Anton, I.: Do Mangalica Pigs of Different Colours Really Belong to Different Breeds?, Arch. Anim. Breed., 49, 477–483, https://doi.org/10.5194/aab-49-477-2006, 2006. 

Zsolnai, A., Lehoczky, I., Gyurmán, A., Nagy, J., Sugár, L., Anton, I., Horn, P., and Magyary, I.: Development of eight-plex microsatellite PCR for parentage control in deer, Arch. Anim. Breed., 52, 143–149, https://doi.org/10.5194/aab-52-143-2009, 2009. 

Zsolnai, A., Tóth, G., Molnár, J., Stéger, V., Marincs, F., Jánosi, A., Ujhelyi, G., Koppányné Szabó, E., Mohr, A., Anton, I., Szántó-Egész, R., Egerszegi, I., Dallmann, K., Tóth, P., Micsinai, A., Brüssow, K. P., and Rátky, J.: Looking for breed differentiating SNP loci and for a SNP set for parentage testing in Mangalica, Archives für Tierzucht, 56, 200–207, https://doi.org/10.7482/0003-9438-56-019, 2013.  

Zsolnai, A., Maróti-Agóts, Á., Kovács, A., Bâlteanu, A., Kaltenecker, E., and Anton, I.: Genetic position of Hungarian Grey among European cattle and identification of breed-specific markers, Animal, 14, 1–7, https://doi.org/10.1017/s1751731120000634, 2020a. 

Zsolnai, A., Kovács, A., Kaltenecker, E., and Anton, I.: Identification of markers associated with estimated breeding value and horn colour in Hungarian Grey cattle, Asian Austral. J. Anim., 34, 482–488, https://doi.org/10.5713/ajas.19.0881, 2020b. 

Download
Short summary
We report one set of the identified molecular markers linked to antler quality in red deer. These markers can be used in the selection procedure to increase the frequency of stags that develop high-quality antlers with elevated probability. Another set of markers has also been introduced that is capable of effective paternity testing.