Molecular Sexing and Species Detection of Antlered European Hunting Game for Forensic Purposes

Simple Summary The reasons behind illegal hunting can be widely different. There are also tricky methods that hunters use in the attempt to legalize their illegally-acquired trophies, specifically, introducing them in later seasons, and registering the eviscerated corpses as hinds. During certain periods of the year, hunters are only able to acquire a license for the shooting of female deer exclusively, with the male-hunting season beginning later. The eviscerated and decapitated carcass of the animal delivered to the wild game processing house cold store can be falsely registered as a female based on phenotype. If the hunting association suspects that the animal’s sex has been falsely reported, a forensic genetic investigation can be proposed. In other criminal events, there are no carcasses at all. When several biological remains left on the crime scene need to be tested, a fast and cost-effective detection of a given species or species-group might be substantial before subsequent analysis. Therefore, sex and species detection can provide an early-stage credibility to the resolution of illegal activities related to trophy animals, and additionally may disclose potential poaching disputes. Abstract Molecular sexing techniques are widely applied in conservation biology, although the range of forensically validated methods is fairly limited. The primary aim of this work was to develop forensically validated assays, using two PCR panels for sex and species assignment for the abundant antlered European game species: red deer (Cervus elaphus), roe deer (Capreolus capreolus) and fallow deer (Dama dama). Segments of the SRY and Amelogenin X/Y genes for sex determination, additionally species-specific cytochrome b regions for species detection were targeted and separately amplified in two multiplex reactions. These assays can reliably analyze trace amounts of DNA. The results of both can easily be visualized and interpreted practically, either on agarose gel or by capillary electrophoresis. These simple, fast molecular assays are able to affect the early-stage resolution of disputed or unsolved poaching cases, without the need of individualization or sequencing of forensic samples.


Introduction
Species identification and sexing methods have significant relevance in conservation biology, population ecology and also in wildlife forensics. It is reasonable that the genetic methods represent a major part of the current analytical toolbox for these tasks, particularly when partial, less recognizable, eviscerated remains of the animal are the single evidence for a given offense. This is not a rare scenario at wildlife crime scenes, especially in cases Several recent methods of sex determination have been developed for the detection of single Y-chromosome specific sequences, completed with nuclear or mitochondrial internal controls to confirm the adequacy of reaction [12,[25][26][27][28][29][30]. However, the absence of a signal male specific marker does not necessarily mean that the sample is of female origin. The allelic drop-out caused by sample degradation or technical imbalances can also lead to false negative results [35,42]. To avoid this type of misinterpretation, the simultaneous detection of at least two different Y-chromosome specific sequences is the most reliable method. This is the case not only for humans or primates, but for other species as well. A complementary sexing protocol developed for detection of the X and Y chromosome linked Amelogenin gene or STR markers can disambiguate the confirmation of analytical results on the SRY gene [43,44].
For species identification, two genetic markers on the mitochondrial genome are mostly analyzed by sequencing: the cytochrome oxidase I (COI) gene, which is the basis of DNA barcoding initiative and BOLD database (Bold Systems v4; http://boldsystems.org/, accessed on 22 November 2021.) and the cytochrome b (Cytb) gene. Both meet the criteria for species assignment, such as high inter-species variability and low intra-species variation [45], and therefore commonly applied in forensic cases [12]. Multiplexing several primers together in one reaction mixture enables simultaneous identification of target species and improves the success rate of species detection when using only a minute amount of DNA. In contrast to the sequencing-based methods, combining species-specific primers with fragment length analysis, real-time PCR detection or multiplex hybridization assay allow the parallel identification of more species and can be applied to a DNA mixture [3,[46][47][48][49]. As many crime scenes are outdoor "living crime scenes", evidential samples are often exposed to environmental breakdown, and DNA shows considerable degradation over time, which can inhibit the subsequent DNA analysis. Success of evaluating trace amounts of DNA depends on the length of the amplified fragment. Methods, using shortened amplicons are able to significantly increase the success rate of analysis from degraded remnants [50].
Although a combination of rapid molecular assays for species detection and sex determination in the case of swamp deer (Rucervus duvaucelii duvaucelii) has previously been described for conservation purpose [31], the criteria of forensic validation as well as extending of testing to European deer species are not mentioned or referred to. Red deer (Cervus elaphus), fallow deer (Dama dama) and roe deer (Capreolus capreolus) are not only among the most important game species in European and Near Eastern countries, having a wide distribution in many regions, but they can also be found on other continents. The objective of this work was to create a relatively fast and validated assays for molecular sexing and species detection for these three above mentioned species, which is suited for simultaneous visualization by capillary electrophoresis (CE) and is easy to use in poaching or illegal trading cases.

Sample Collection and DNA Extraction
All samples, such as muscle, hide, hair and feces were collected from registered shootings of the given species. Red deer (n = 30), fallow deer (n = 47) and roe deer (n = 38) with known sex were provided by the Pilis Park Forest Hunting Company and by hunters with a license (Table S1, Figure S2).
Genomic DNA was isolated from the muscle, hide and hair samples using a DNeasy ® Blood and Tissue Kit (Qiagen, Hilden, Germany), and a QIAamp ® Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) was used for the feces samples, following the procedural guidelines as instructed. Extracted DNA quality was measured on 1% agarose gel using a GelRed TM Nucleic Acid Gel Stain (Biotium, Fremont, CA, USA), and concentration was measured by a Qubit 2.0 Fluorometer (Life Technologies Corporation, Carlsbad, CA, USA). Extracted DNA from the tissue samples were stored at 4 • C until subsequent analysis.

Marker Selection and Primer Design
DNA sequences of Sex-determining Region Y (SRY), Amelogenin X and Y (AmelX/Y) and cytochrome b (Cytb) genes of the three species were obtained from the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/genbank/, accessed on 28 November 2021).
The selected homologue sequences of SRY, Amel X, Amel Y and Cytb genes were aligned in MEGA (Molecular Evolutionary Genetics Analysis across computing platforms) [35] to identify conservative regions of similarity in the case of AmelX, AmelY and SRY genes. In case of the Cytb gene, the species-specific SNP (single nucleotide polymorphism) variations were identified and compared for developing species-specific forward primers, the reverse primer was designed universally for all the three species. For gender analysis, the universal primers were created for the three given species using the conserved parts of the AmelX, AmelY and SRY markers. All primers were designed by Primer Designer 4 software (http://www.scied.com, accessed on 28 November 2021). Selected primers have the length of 18-36 bases and no self-complementarity (hairpins) or complementarity to other primers (more than four base pairs at the 3 end) were observed. Moreover, the primers were designed considering the reduced amplicon size of <220 bp to assure higher amplification success and the varying amplicon sizes among markers enable the parallel detection with capillary electrophoresis (Table 1). In order to examine possible cross-amplification with the goal of discovering limitations or widened utilization fields, DNA sequences of each marker from a great variety of related cervid and other species were obtained from the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/genbank/, accessed on 28 November 2021) for further in silico analysis. Intra-population diversity at the primer binding sites for each marker and investigated species were also tested by MEGA platform.

Singleplex Reactions and Sequencing Analysis
Amelogenin, SRY and cytochrome b gene fragments were amplified with male and female red deer, roe deer and fallow deer samples in singleplex reactions to determine functionality. PCR (25 µL in volume) containing 12,5 µL DreamTaq™ Green PCR Master Mix (ThermoFisher Scientific, Waltham, MA, USA), 1 µL of BSA (10 mg/mL, Sigma-Aldrich, St. Louis, MO, USA), 0.5 µM unlabeled forward and 0.5 µM unlabeled reverse primers, 5 ng DNA template and PCR grade H 2 O to volume. PCR was carried out in an AppliedBiosystems 2720 Thermal Cycler with the following conditions: an initial 95 • C for 30 s followed by 36 cycles of 25 s at 94 • C, annealing of 25 s at the given marker specific temperature (as described in Section 2.4) and 25 s at 72 • C, and final extension for 20 min at 72 • C, for subsequent sequencing reactions. Amelogenin X and Y alleles were separated on 2.5% agarose gel, the selected allele bands were cut from the gel and purified with Gel Advanced TM Gel Extraction System (Viogene, Sunnyvale, CA, USA). The reamplification parameters were the same as singleplex PCR (see above) using 26 cycles. All PCR amplicons were purified using GenElute TM PCR Clean-Up Kit (Sigma-Aldrich, St. Louis, MO, USA). Both DNA strands were sequenced using the BigDye ® Terminator v.1.1 Cycle Sequencing Kit (Thermo Fisher Scientific, Waltham, MA, USA) in the manner recommended by the manufacturer. For sequence detection, an ABI Prism 3130XL Genetic Analyzer (Applied Biosystems, Waltham, MA, USA) was applied, according to manufacturer's guidelines. Sequence analyses were performed using Sequencing Analysis Software 5.1 (Applied Biosystems, Waltham, MA, USA) and aligned by Sequencher TM 4.1.2 software (Gene Codes Corp, Ann Arbor, MI USA). A homology search of resulted sequences was accomplished using the GenBank (BLAST ® , Basic Local Alignment Search Tool).

Multiplex Amplification and Detection
Annealing temperatures were optimized for a diplex use of the sex determining loci, hereafter referred to as "DeerSex-plex", and in a separate triplex reaction for the species detection markers, hereafter referred to as "DeerSpec-plex". PCR conditions were optimized as follows, using female and male genomic DNA obtained from each of the three above mentioned species. The primer concentrations were combined equimolar first, then based on the peak heights of the capillary electropherograms, variable concentrations were set to obtain inter-locus balance (Table 1). Negative (no-template) reaction controls were applied in each step of analysis to avoid/detect the risk of contamination.
Amplification of AmelX/Y and SRY were performed in DeerSex-plex in a 10 µL reaction volume, containing 5 µL DreamTaq™ Green PCR Master Mix (Thermo Fisher Scientific, Waltham, MA, USA), 0.5 µL of BSA (10 mg/mL), 1.5 µL of primer-mix, 2 ng of DNA template and PCR grade H 2 O to volume. Touch-down PCR was carried out on Applied Biosystems 2720 Thermal Cycler using the following conditions: an initial step at 95 • C for For species detection, three specific forward primers and one universal reverse primer were designed to amplify the Cytb gene segment in a 10 µL reaction volume in DeerSpec-plex, containing 5 µL DreamTaq™ Green PCR Master Mix (ThermoFisher Scientific, Waltham, MA, USA), 0.5 µL of BSA (10 mg/mL), 1.5 µL of primer-mix (three species-specific forward and the universal reverse primer), 2 ng of DNA template and PCR grade H 2 O to volume. PCR was carried out on Applied Biosystems 2720 Thermal Cycler using the following conditions: an initial step at 95 • C for 30 s followed by 34 cycles of 25 s at 94 • C, annealing of 25 s at 64 • C and 25 s at 72 • C, and final extension for 20 min at 72 • C.
From both (sex-and species detection) assays, the amplified products were separated and visualized on a 2% agarose gel using GelRed ® nucleic acid stain (Biotium, Fremont, CA, USA). The two reactions were then analyzed by capillary electrophoresis on an ABI Prism 3130XL Genetic Analyzer using GeneScanTM-500 LIZ TM Size Standard (Thermo Fisher Scientific, Waltham, MA, USA). The minimum detection threshold was set at 150 relative fluorescence unit (RFU) during fragment analysis using GeneMapper ® ID-X software version 1.4. Bin-width used was 1 bp, which was set for allele designation for each investigated marker based on allele sizing of sequencing verified alleles.

Basic Validation
To estimate the reliability and limitations of the two assays, the steps of basic developmental validation were implemented according to recommendations [32,51].
Cross-species amplification of both plexes (DeerSex-plex and DeerSpec-plex) were tested on male and female samples (input amount of 2 ng) from related species such as cattle (Bos taurus), takin (Budorcas taxicolor), bison (Bison bison), mouflon (Ovis musimon) and sheep (Ovis aries); and on samples from forensically relevant species such as human, swine (Sus scrofa), dog (Canis lupus familiaris) and cat (Felis catus). The success rate of intra-species amplification was tested by both plexes on 115 samples-red deer (n = 30), fallow deer (n = 47) and roe deer (n = 38)-collected from different regions of Hungary.
Sensitivity was evaluated using 1.0, 0.5, 0.25, 0.125 and 0.062 ng DNA templates in each reaction. Each DNA dilution was amplified with both assays two times on different days using the same dilution series.
Specific alleles of the sensitivity and reproducibility studies detected by CE were used for calculating peak height ratios. Heterozygote balance in the case of AmelX and AmelY marker in male animals was calculated by dividing the peak height of the higher peak by the peak height of the lower peak.
Both plexes were tested on forensically relevant sample types, such as muscle, hide, hair and feces. Five previous casework DNA samples, including blood stain, blood swab and hairs from different poaching cases involving red deer and fallow deer collected in 2019 and 2020, were reanalyzed. Conditions of DNA preparation and amplification are described in Section 2.4.
For the study mixture, equal DNA amounts (1-1 ng) were used to prepare the red/roe/fallow deer DNA mixture. Additionally, 5 ng human DNA was also mixed with 0.5, 0.25, and 0.125 ng red-, roe-, and fallow deer DNA in separate tubes. These mixtures were examined with the DeerSpec-plex assay described in Section 2.4. To estimate the sensitivity of male specific markers' detection 5 ng female red-, roe-, and fallow deer DNA were mixed with 1, 0.5, and 0.33 ng male red-, roe-, and fallow deer DNA in separate tubes and examined with the DeerSex-plex assay described in Section 2.4.
The accuracy was established through repeatability and reproducibility studies. The multiplexes were tested by one operator on two runs, and by different operators on separate days on two runs.

Field Samples-Mock Cases
In two cases, five field samples-blood, vessel, muscle, bone and hair samples-from two red deer corpses were provided by two official veterinarian examiners of Hungarian hunting companies ( Figure S3). The preliminary hypothesis in both cases was that the samples had to come from red deer hinds, which were delivered into processing houses.

Results
DNA was successfully extracted, and measurable amounts (0.8-96.4 ng/µL) of DNA were obtained from blood stain, muscle, hide, hair and feces.

Marker Selection and Primer Design
The SRY and AmelX, AmelY loci were selected for sex determination and four primers were designed to amplify both gene regions for the three species of interest. One reverse primer and three species-specific forward primers based on the SNP pattern of the Cytb gene of red deer, roe deer, and fallow deer has been created for species detection (Table 1). Intra-and inter-species specificity were estimated using as many sequences as available from GenBank for in silico analysis (Table S2, Figures S4-S6). Regarding the four sexing primers (SRY forward and reverse, AmelX/Y forward and reverse), only one SNP position within the AmelX/Y reverse primer was observed in the downloaded sequences (in roe deer species). Within the Cytb reverse primer, only one SNP affected the primer binding site in fallow deer and two SNP in red and roe deer samples. The forward primer sequence of Capreolus capreolus Cytb contains three possible SNPs from 135 downloaded roe deer sequences. However, based on 435 downloaded red deer sequences, a relatively high number of polymorphic base positions (nine SNPs) can be existing within the Cytb Cervus elaphus forward primer binding site, five of them have only <1% frequency. Within a single sequence, the most common SNP numbers were two in primer binding sites of the Cytb marker in red-and roe deer samples. No SNP positions were observed at the 3 end position within the six developed primer sequences from neither species.
Having sequences available from a great variety of species found in GenBank, similarities among species were found based on the in silico analysis using MEGA software ( Figures S4-S6). The results show that the designed primers for sexing could be used in several other cervid species in the case of SRY and Amelogenin X/Y genes (Figures S4 and S5).

Singleplex Reactions, Sequencing and GenBank Submission
The PCR products from each marker were confirmed by sequencing. The sequences were submitted into GenBank (accession no. MW876187-MW876198), which indicated homology with the Amelogenin X and Y, SRY and Cytb genes of red deer, roe deer and fallow deer. Possible SNP positions within the primer sequences can cause a partial homology in these species, but the specificity of the result was not influenced. Primers do not require full matching in the middle section (Table 1).

PCR Based Multiplex Design and Detection
Our original goal was to amplify all selected markers for sex-and species detection in one multiplex reaction, but this attempt failed due to the significant differences among the optimal annealing temperatures, which induced allele/locus drop-out and artefact peaks. Accordingly, primers of sex-and species detection markers were developed and optimized successfully in two separate multiplexes. The specific fragments in question were robustly amplified and detected on agarose gel from control samples with known sexand species controls ( Figure 1). As expected, both Y-chromosomal (SRY and AmelY) and X-chromosomal (AmelX, as a reaction control) segments were detected in male samples, while only the AmelX was detected in female samples using the DeerSex-plex assay. Although some aspecific fragments were also produced, the artefacts did not influence the judgment of the sex, as their sizes are longer than the specific Amelogenin X allele. All of the three species-specific primers in the DeerSpec-plex designed for the cytochrome b gene amplified only the target cervid species using the optimized multiplex PCR profile (Figure 1). To avoid cross species or nonspecific amplifications during testing, annealing temperature should reach 64 • C.
Due to the different amplicon lengths, followed by PCR amplification and agarose gel electrophoresis, a parallel detection of each PCR product from both assays was feasible by capillary electrophoresis using fluorescent labeled primers with 6-FAM dye (Table 1, Figure 2). The CE detection can be implemented in optionally mixed form, in which based on the strength of the amplicons' bands on the semiquantitative agarose gel an approximate 1:1 mixing ratio of the two plexes can be appropriate. The males show two Y-chromosomal peaks (SRY, AmelY) and the AmelX as a reaction control, plus a species specific peak (Cytb); females show a single sex chromosomal peak (AmelX) as a reaction control, plus a species specific peak (Cytb). Based on the size of the amplified cytochrome b gene, the species origin can be clearly determined. Some artefact peaks were also detected out of the bin sets, these did not influence the reliability of the evaluation of results. As expected, no-template reaction control samples were negative in each reaction by gel electrophoresis and by capillary electrophoresis.

Species Specificity
Species specificity was separately tested on male and female samples in the case of the DeerSex-plex to know whether detected specific or aspecific PCR products have relatedness with the sex, and on female samples regarding the DeerSpec-plex assay (Figures S7 and S8).
DeerSex-plex assay: in the case of mouflon and sheep, three amplicons were detected at SRY, AmelY and AmelX bin positions on male samples, and one amplicon was detected at AmelX bin on female samples. No specific AmelX allele was observed using samples from other species, neither from males nor females. Using male cattle sample, PCR product at SRY and AmelY bin could be found specifically to the male sex. Using bison, dog and cat samples, PCR product at SRY bin could be found specifically to the male sex, as other amplicons at AmelY bin position or out of range could be observed on female samples as well. In cases of human, swine and horse only, one amplicon fell into (or very close to) bin position of AmelY marker, but from both sexes. No amplicons within the specified bin sets on the capillary electropherogram were observed using takin DNA, neither from the male nor female animal ( Figure S7).
DeerSpec-plex assay: The designed primers for cytochrome b gene segments showed limited species specificity since two PCR products manifested very close to the marker specific bin settings in the case of the mouflon and sheep, and one amplicon using bison DNA ( Figure S8). However, no PCR product has been detected using human, cat, horse, swine, takin, cattle or dog DNA.
Intra-species specificity study based on the analysis of collected population samples shows that all investigated samples amplified with DeerSpec-plex provide a clear PCR product on agarose gel ( Figure S9). Using the DeerSex-plex assay, all but one sample (Dd_16) could be effectively determined for the sex by agarose gel detection ( Figure S10).

Sensitivity, Mixtures and Case-Type Samples
Sensitivity of multiplexes were evaluated using template concentration in dilution interval of 1 to 0.06 ng. Clearly visible specific bands were observed in each three species until 0.125 ng dilution using the DeerSpec-plex assay and until 0.5 ng dilution using the DeerSex-plex assay. Below this DNA amount, AmelX/Y or SRY locus drop-out could be observed ( Figures S11 and S12).
Using DNA mixtures from red deer, roe deer and fallow deer, all three species could be detected until 0.06 ng dilution using the species identification assay. The detectability of male specific markers proved to be sensitive enough to see them on agarose gel using 5 ng female and 0.33 ng male red and fallow deer DNA. In the case of roe deer, 5 ng female and 0.5 ng male DNA gave a visible result of each allele involved in DeerSex-plex assay ( Figure S13).
Extracted DNA from blood swab, blood stain, hairs, muscle and hide could be successfully amplified, resulting specific amplicons clearly visible on agarose gel and on capillary electropherogram as well (Figures S14-S16).
The electropherograms obtained in the course of the sensitivity study provided the basis for the calculation of intra-locus peak height ratios. The peak height ratios for AmelX/Y ranged from 0.66, 1.11 and 1.51 to 9.26, 2.29 and 11.26, with an average value of 2.35, 1.67 and 3.85 in roe deer, red deer and fallow deer (respectively).

Field Samples-Mock Cases
The field sample (vessel) in case 1 resulted in red deer specific Cytb, and male-specific fragments, so the initial hypothesis, namely the sample originated from a female was excluded ( Figure S15). The field sample (corium) in case 2 resulted in red deer specific Cytb, and female-specific fragments, so the initial hypothesis was supported by analysis ( Figure S16).

Discussion
Poaching is a criminal activity, which is strictly defined by laws regulating the taking of wildlife. The literature is filled with studies of poaching cases worldwide, since the illegal killing of wildlife has been a social problem for centuries [10]. Parallel with the changes in law, new, highly reliable, molecular techniques and markers for sex-and species differentiation are also needed. In compliance with the SWGDAM [51], SWFS [52] and ISFG [32] recommendations, standards and guidelines for animal and wildlife forensic analysis, these methods can be incorporated into forensic practice only when systematic validation steps are performed.
For this purpose, we set up and validated two assays for sexing-and for species assignment, which generated reproducible and reliable data from a range of sample types. The validation process regarding the sensitivity, specificity, mixture study and repeatability was accomplished to the forensics standards; therefore the assays are applicable to the enforcement for the detection of both species and sex for red-, roe-and fallow deer species.
The sexing panel simultaneously amplifies two conserved regions of the Y-chromosome, the SRY and the Amelogenin Y gene, and one X-chromosomal locus (AmelX) providing an internal control. Double testing of Y-chromosome loci provides a result with higher accuracy, thus decreasing the chance of a false-positive or false-negative result. The assumption that the Amelogenin gene exists on both X and Y chromosomes with a deletion region on the Y chromosome in these species was supported by our results, as DNA fragments of different length between the male and female animals can be detected. As only red deer from the three tested species has a nine base pair deletion within the amplified AmelY region (Table 1, Figure 2), the sexing assay by itself can also indicate the species origin. However, based on GenBank sequences, this deletion exists also in Sika deer (Cervus nippon). In addition to the investigated cervid species, specific amplicons were detected from the nearest taxonomically related species (mouflon and sheep) for each marker, but with quite low signal intensity at the AmelY marker. However, since some other minor PCR products outside of the bin settings were also manifested based on sequences aligned in MEGA, the possible utilization may extend to these species and other cervid species as well ( Figures S4-S6). Major artefacts affecting the interpretation of the sex were not observed.
Amplifying short mitochondrial DNA regions with species-specific primers, using different point mutations within the forward primer 3 binding sites, is an effective method for simultaneous detection of mammalian species [3]. Species-specific primers can be designed based on inter species SNPs, provided those SNPs have little or no intraspecies variation [53]. In this study, all the four sexing primers (AmelX/Y forward and reverse, SRY forward and reverse) and the Cytb Dama dama forward primer proved to be designed to a highly conserved region within species, however limited number of GenBank sequences were accessed (Table 1 and Table S2). Although within the Cytb primers binding sites one to four SNPs can occur, drop out of specific fragment was not observed in our study ( Figure S9). In cases where there is a well-defined investigational question i.e., a nominated set of wellcharacterized species, an assay using primers to target SNPs can be a cost-effective and rapid alternative to sequencing, followed by agarose or/and capillary electrophoresis [1].
Based on this assumption, there have been multiplex tests developed to detect a small number of deer species [54,55]. As it cannot be excluded that unforeseen species may produce a signal that is similar in length but not homologous with the target taxon, these tests are often considered presumptive methods of detection and not identification [56]. These markers, however, can be used as a screening tool to categorize samples for sequencing and are particularly suited for situations where sequencing commonly fails, such as when species mixtures are expected (e.g., common food meats or animal based traditional medicine) [50,57,58]. Such assays are best suited to situations where there is a limited number of expected species; they also require more extensive validation by sequencing to ensure the specificity of results [56].
Among ungulates, the evidence of hybridization varies in different countries [59][60][61][62][63]. The volume and effects of anthropogenic hybridization and wildlife translocations can differ region to region, as well as the efficiency of efforts of wildlife conservation or national prevention [21]. However, where there is a possibility of the presence of hybrids, the single utilization of maternal inherited mt DNA marker is not appropriate. The limitation using mitochondrial genes for species assignment can be solved with the expanded examination tool of paternal and/or biparental markers, for instance microsatellites or SNPs [63][64][65][66][67][68]. However, it is noted that the determination of species origin of hybrids requires more complex genetic methodology involving statistical analysis [69,70].
In our developed DeerSex-plex and DeerSpec-plex assays, the different fragment sizes were well separated in agarose gel electrophoresis; thus, different sex and species origin were easy to identify. However, fluorescent gel electrophoresis adaptation was implemented for better resolution by labeling the PCR products with fluorescent dyes. Combining the two panels in one capillary electrophoresis detection gives a rapid and prompt result (attention should be paid to use a proper mixing ratio of the two reactions) and makes sample identification more reliable and cost-effective. These approaches facilitate fast screening of evidential samples of various qualities collected from crime scenes, reduce time and cost in comparison with traditional sequencing or restriction fragment length polymorphism-based methods. Such molecular approaches are definitely advantageous when there is a need to test several samples for whether they come from the relevant species or sex in connection with the given case. In this study, we obtained no misidentification regarding the sex or species of the three tested cervid species; however, some imbalance among the specific peaks as well as some non-specific peaks could be observed.

Conclusions
Simple and precise methods for sex and species assignment in animals are a prerequisite for several applications in animal production, wildlife conservation and forensics. This study targeted species and sex detection for forensic purpose of three coexisting European cervid species as hunting games. Our methods, which match to forensic standards, could contribute credibility to the resolution of potential poaching disputes and confine the illegal activities related to red deer, roe deer and fallow deer.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/ani12030246/s1, Figure S1: Eviscerated and decapitated carcasses in processing house cold store. Table S1: Sample database from three antlered species used for validation studies. Figure S2: Distribution of collected samples within Hungary from three antlered species. Figure S3: Field sample types. Table S2: In silico intra-species variation study of developed primers. Figure S4: Conservative regions of the SRY gene identified among species by MEGA. Figure S5: Primer binding regions marked out with yellow color within the highly conservative regions of the AmelogeninX/Y gene among species identified by MEGA. Figure S6: The species-specific forward and the universal reverse primer binding regions within the cytochrome b gene among species identified by MEGA. Figure S7: Species specificity tested on male and female samples using the DeerSex-plex assay, detected on agarose gel and by capillary electrophoresis. Figure S8: Species specificity tested on female samples using the DeerSpec-plex assay, detected by agarose gel and capillary electrophoresis. Figure S9: Intra-species specificity study, population samples from three species amplified with DeerSex-plex and detected on agarose gel. Figure S10: Intra-species specificity study, population samples from three species amplified with DeerSex-plex and detected on agarose gel. Figure S11: Sensitivity tested on male and female samples using the DeerSex-plex and DeerSpec-plex assays on agarose gel. Figure S12: The SRY and AMELX/Y loci in red (a), roe (b) and fallow (c) deer samples amplified with varying amounts of input DNA with the DeerSex-plex assay. Figure S13: Mixture study: red deer, roe deer and fallow deer samples amplified with varying amounts of input male DNA with the DeerSex-plex assay. Figure S14: Case type samples amplified with both assays. Figure S15: CE result of field sample (vessel) in mock case 1. Figure S16: CE result of field sample (corium) in mock case 2.

Funding:
The project was supported by the Hungarian Hunters' National Chamber and by the European Union and co-financed by the European Social Fund: EFOP-3.6.3-VEKOP-16-2017-00005-"Strengthening the scientific replacement by supporting the academic workshops and programs of students, developing a mentoring process".
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: The data will be available with the corresponding author upon request.

Conflicts of Interest:
The authors declare that they have no interests that could have appeared to influence the work reported in this paper.