{'fraction_non_alphanumeric': 0.034269680633059046, 'fraction_numerical': 0.01803335540212491, 'mean_word_length': 5.006438174588146, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes

Strigiformes are affected by a substantial decline mainly caused by habitat loss and destruction, poaching, and trapping. Moreover, the increasing trend in bird trade and the growing interest in wild-caught rather than captive-bred birds are expected to encourage illegal trade. The biomolecular investigation represents a valuable tool to track illegal trade and to explore the genetic variability to preserving biodiversity. Microsatellite loci (STRs) are the most used markers to study genetic variability. Despite the availability of species-specific microsatellite loci in Strigiformes, a unique panel permitting the description of the genetic variability across species has not been identified yet. We tested 32 highly polymorphic microsatellite markers to evaluate the reliability of a unique microsatellite panel in different species of Strigiformes and its use for conservation and forensic purposes. We included in the study 84 individuals belonging to 28 parental groups and 11 species of Strigiformes. After screening polymorphic microsatellite loci, the description of genetic variability, and the kinship assessment, we characterized a final panel of 12 microsatellite loci able to identify individuals in 9 Strigiformes species. This STR panel might support the authorities in the forensic investigation for suspected smugglers and false parental claims; moreover, it can be useful to evaluate relatedness among individuals in captive-bred populations and to implement research projects finalized to the description of the genetic variability in wild populations.


Introduction
Habitat loss and fragmentation represent the main threats to the survival of wildlife species. Several studies have well documented their effect on biodiversity loss [1]. Moreover, also international wildlife trade contributes to the depletion of natural resources; hence, the over-exploitation of natural populations, such as the withdrawal of individuals from the wild, is considered one of the leading causes driving the species into extinction [2][3][4]. It is known that seafood, furniture, and fashion are the main categories requested from the international trade; in addition to this, the commerce of pets affects many individuals. At least 5% of the import/export requests regards parrots and bird of prey (raptors and owls), a quote that equals the sum of all other commercialized birds [1].
Baker and colleagues [5] calculated that animals' demand for pets and entertainment purposes contributed to at least one-fifth of the wildlife trade [6][7][8], and they showed that the removal of wild individuals from native populations was accounted in migrating species for the second threat that pushes species into extinction, lower only to the habitat loss [9].
Species over-exploitation increases the negative effects produced by habitat unsuitability, causing the reduction of individuals in wild populations that, in turn, can produce a fast loss of genetic variability in a short timescale. Domínguez et al. [10] showed that decreasing populations and pet trade could have generated a bottleneck, reducing genetic variability in the yellow cardinal, Gubernatrix cristata, an endangered passerine from the southern region of South America. Harris et al. [11] found that 14 species routinely traded in Indonesian wildlife markets had undergone population decline differently from the other 22 untraded species, which did not suffer any number reduction. Evans and Sheldon [12] found a correlation between the decline of mean heterozygosity with an increasing extinction risk and identified the smaller population sizes as the main cause of this correlation.
Valuable tools to rule or forbid the trade of threatened species to save them from extinction have been provided since the application in 1976 of the Washington Convention. All the countries belonging to the European Community joined the convention that has been applied through several Council (EU) Regulations. Species subjected to maximum protection are listed in Annex A that includes most of the species listed in Appendix I and several species of Appendix II recorded as critically endangered in the European Community.
One of the orders suffering a strong decline in Europe is Strigiformes, which includes two main families: Tytonidae and Strigidae. While the first family is split into two subfamilies, each of which including one genus (Tyto and Phodilus), the taxonomy of Strigidae is quite controversial [13]. Differently from traditional systematics, Wink et al. [14] recognized the existence of three subfamilies (Ninoxinae, Striginae, Surninae) with several tribes in the last two through a molecular phylogenetic approach. Concerning the substantial decline of Strigiformes in Europe, it is caused mainly by pesticides, changing agricultural techniques with the loss of the structures for nesting and lower rodent availability than in traditional farming, and caught of individuals from the wild. Several main factors threatening owl survival in the Gaza Strip were identified by Abdel Rabou [15], including habitat loss and destruction, poaching and trapping, myths and superstitions, secondary poisoning, road kills, and fence agricultural lands. Wan et al. [16] found a negative influence of fragmentation on the population size and genetic diversity. Despite their inclusion in CITES Appendix II, they enjoy the greatest protection and are listed in Annex A. As for other species included in the CITES Appendix II, only individuals imported into the EU before the Convention of Washington entered into force can be traded. This permission is also applicable to these individuals' descendants, but only if their birth in captivity is proved.
The demand for pets among birds gradually increased in the past decades [17]. Panter, Atkinson and White [7] recorded a wild-caught export of 18,948 individuals belonging to 86 owl species from 1975 to 2015. In 2019, Ribeiro and colleagues [18] denounced a future trend in legal bird trade driven by sociocultural motivations with an increased demand that will interest wild-caught rather than captive-bred birds. This preference is expected to drive towards an increase in the illegal trade of the rarest or more popular species. Also, the reduced number of wild individuals increases their commercial value, encouraging, even more, their illegal trade. Remarkably, the illicit traffic involves the smuggling of eggs and the laundering through captive breeding facilities of wild-caught animals [7].
Owls trade increased in the last decades and has also been influenced by a change in habits; for example, Nijman and Nekaris [19] recorded an increased owl trade, from <0.06% before 2002 to >0.43% post 2008, in wildlife markets in Java and Bali and suggested a delayed "Harry Potter effect". Siriwat and colleagues [20] identified the same increasing trade in Thailand without finding a correlation with the "Harry Potter effect"; nevertheless, they found an association between the increasing owl request, the novel online market and more popularity of owls as a pet. Moreover, they related the higher market price of some species to lower individual availability.
Even though the main actions aim to prevent and hamper over-exploitation and illegal trade represent the first steps in preserving biodiversity, a biomolecular investigation can give valuable hints for either verifying the depletion of genetic variability or tracking illegal trade. This can be carried out through individual identification using molecular techniques and kinship analyses. One of the most used genetic markers to study genetic variability is microsatellite loci (i.e., short tandem repeats or STRs). They are multiallelic and PCR multiplexable, allowing the description of genetic variability with reduced costs. Although the comparability between different laboratories has often proved to represent a limit for their use, this hampering can be easily resolved by using an allelic ladder for each marker-that identifies the position of all the alleles-and by sharing reference samples.
Despite the availability of species-specific panels of microsatellite loci in Strigiformes [21][22][23][24][25][26] and the cross-amplification of some markers in species lacking their markers [27], a unique panel permitting the description of the genetic variability across species has not been found yet. Such a protocol would permit the recording of comparable variability indices among different species and speed the process, reducing molecular analyses costs; it could also be informative in detecting illegal trade of individuals collected from the wild.
Thus, we tested 32 highly polymorphic microsatellite markers available from literature in 11 species of Strigiformes with the following aims:

•
To evaluate the reliability of a unique panel of microsatellite loci for several species of Strigiformes; • To test its use for conservation and forensic purposes; • To assess the use of the panel in confirming phylogenetic relationships among species.

Materials and Methods
The overall experimental strategy used in this work was based on the following steps: (i) Thirty-two microsatellite markers selected from the literature were tested on two unrelated individuals of a restricted panel of three species (Athene noctua, Strix aluco, Bubo bubo) (Table 1). Table 1. List of species and individuals used in each step: family groups from 1 to 9 were tested in the first and second screening; additional species and family groups (from 10 to 28) were tested in the third step. (ii) Microsatellite markers giving polymorphic profiles in all the species tested at step (i) were used for genotyping three family groups per species (Table 1).
(iii) Microsatellite primers giving polymorphic amplification patterns at step (ii) were then used to evaluate the STR panel for forensic purposes on a larger panel of family groups, comprising different subfamilies, tribes and species, aiming at identifying a unique shared panel of polymorphic loci ( Table 1). The selected species are the most commonly traded in the Italian market.
According to the manufacturer's instructions, DNA was isolated from feathers using DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany). After digestion in 180 µL ATL buffer, 20 µL proteinase K and 20 µL DTT and incubation overnight at 56 • C, the lysate was loaded in a QIAcube HT robotic station (Qiagen, Hilden, Germany) for further purification steps.
DNA was amplified in an 8 µL final volume reaction. The PCR reactions were carried out as follows: in 8 µL of the final volume, we added 1x reaction buffer, 0.02% BSA, 1.5 mM MgCl 2 , 0.125 mM of dNTPs, Hot Start Taq Polymerase 0.025 U (Qiagen, Hilden, Germany), 0.125 mM primer (forward and reverse), 1 µL DNA template and nuclease-free water to reach the final volume.
Twenty-five out of 32PCR primers used in this work were originally isolated from five species of Strigidae (Bubo bubo, Otus elegans, Strix occidentalis, Strix nebulosa, Glaucidium brasilianum), one from the goshawk (Accipiter gentilis). Six were previously designed to and/or shown to display high cross utility throughout many genetically distant passerine and shorebirds and were found polymorphic in Tyto alba. (Table 2); each forward primer was labelled with fluorescent ABI dyes. Amplification reactions were performed in simplex to test for the right allelic range and no specific amplification. According to melting temperature and reference bibliography, we chose the following thermal profile: 94 • C for 15 , followed by 35 cycles at 94 • C for 40 , 60 • C for 40 , ending with a final extension at 72 • C for 10 .  Oe129 GTCACCTCTTGACATCCGAGTAGC GCTAAGAGTCCATTTGCCCATCTG Otus elegans [29] Amplicons were separated through capillary electrophoresis in an ABI 3130xl genetic analyzer (Thermo Fisher Scientific); alleles were scored in GeneMapper 4.0 using GeneScan 500 ROX size standard (Thermo Fisher Scientific, Waltham, MA, USA). An allelic ladder was constructed in Genemapper, merging the sample electropherograms at each locus. Scoring has been reported for each peak, allowing the data comparison among laboratories (Supplementary Material Figure S1: allelic ladders).
Genotypes resulting from the screening were analyzed for the selected loci. Allele number (NA), effective allele number (NE), number of private alleles, observed and expected heterozygosities (Ho, He) were computed using GenAlEx 6.41 [32]. Gimlet 1.3.3 [33] was used to evaluate the probability of identity (PID) among individuals of each cluster species, to calculate the minimum number of markers necessary to achieve a reliable individual genotyping for unrelated (PID) and related samples (PID _sib , sibling).
A factorial correspondence analysis (FCA) was carried out in Genetix 4.05 [34] to check the genetic distances between species, tribes and subfamilies.
Parentage relationships among individuals within the family cluster were reconstructed in Colony 2.0 [35]. This analysis allows to check if an individual that is claimed to be bred in captivity truly descends from the declared parents. A single analysis was performed including all genotypes of all species together. We used non-inbreeding data and monogamy models, and set the genotyping error rate value 0.0001. The same analysis was conducted in Cervus [36,37] and the LOD score was computed using a proportion of loci typed = 0.60 and mistyped = 0.01. An additional computation with the same parameters was performed using the Delta score defined as the difference in LOD scores between the most and the second most likely candidate parents. The confidence levels were set to 80% (relaxed) and 95% (strict) in both analyses.

Evaluation of the STR Potential in Family Groups
Three family groups (father, mother and one offspring) for each species (n = 9, Table 1) were chosen from the CITES database and analyzed at the twenty-one loci selected in the previous section.

Evaluation of the STR Potential in Other Subfamilies of Strigidae
In order to evaluate the potential application of the STR panel for forensic purposes, the 12 microsatellites were used on a larger panel of samples (Table 1), which included: 1. Additional species belonging to the already analyzed tribes of Bubonini and Strigini (Bubo scandiacus, Strix uralensis, Strix nebulosa); 2. Two species belonging to the tribes of Asionini (Asio otus) and Otini (Otus scops); 3. Two species belonging to the Surninae subfamily, Tribes Surnini (Surnia ulula, Glaucidium passerinum).
4. Tyto alba, belonging to the Tytonidae family, subfamily Tytoninae. Three family groups were chosen for each species except for A. otus, S. uralensis and G. passerinum because only two confirmed parental nuclei were available in the database. The final dataset consisted of 81 individuals belonging to 27 families.
The DNA of each of the 81 individuals was amplified using the twelve primer sets. Data obtained revealed that the DNA from the additional species included in the analysis (O. scops, A. otus, B. scandiacus, S. ulula, S. uralensis, S. nebulosa, G. passerinum, T. alba) was amplified at the examined loci with the following exceptions: 1. S. ulula showed fixed genotypes at four loci (Oe054, Oe129, Oe053, Oe321); 2. Ten loci were amplified in G. passerinum. However, only four loci were polymorphic. Because of lack of variable loci, G. passerinum was removed from the further analysis; 3. S. uralensis and S. nebulosa showed no polymorphisms at Oe321 and Oe054 loci; 4. In A. otus, a unique fixed allele was recorded at locus Oe142; 5. Six loci were amplified in T. alba (FePo42, Oe54, Oe128, Oe129, Oe321, Tgu06), but only 2 were variable (FePo42, Tgu06). Thus, this species was discarded from further analyses because two loci are too few to be able to distinguish individuals reliably.
All the percentage values of polymorphic loci per species varied between 100% (n = 12 in A. noctua and O. scops) and 66.67% (n = 8 in S. ulula), if we do not take into account the value in G. passerinum and T. alba (Table 3).

Genetic Variability between and within Species
Statistical analysis was carried out on 9 species at 12 loci ( Table 3). The mean allele number (NA) was 3.7 (±0.2), with the highest values recorded respectively in O. scops (6.3 ± 0.6) and the lowest one in S. ulula (2.2 ± 0.3). Mean expected and observed heterozygosity (He and Ho) ranged from 0.750 ± 0.027 and 0.733 ± 0.049 in O. scops to 0.318 ± 0.074 and 0.356 ± 0.088 in S. ulula with an average value of 0.525 ± 0.025 and 0.566 ± 0.030, respectively. The probability of identity resulted in different thresholds depending on the species and if it was estimated for unrelated or related individuals (Figure 1). A PID value lower than 0.001 was reached in all the species using an average of 3.9 markers, while the same value was achieved with PID _sib using at least an average of 9.2 loci. S. ulula, S. nebulosa and S. uralensis were not included in this last computation because their threshold resulted higher using all the loci (minimum values respectively of 0.0155, 0.00162 and 0.00246) ( Table 4). A principal component analysis was then carried out, revealing that, as expected, the three species of the genus Strix (S. uralensis, S. nebulosa, S. aluco) overlapped in the PCA graphic (Figure 2). A. otus and O. scops showed a low distance from each other and resulted not very distant from the Strix species. S. ulula exhibited the greatest genetic distance from the other species while A. noctua was in a mid-range position (Figure 2). B. scandiacus and B. bubo did not overlap, but the genetic distance of the two species is low.  A principal component analysis was then carried out, revealing that, as expected, the three species of the genus Strix (S. uralensis, S. nebulosa, S. aluco) overlapped in the PCA graphic ( Figure 2). A. otus and O. scops showed a low distance from each other and resulted not very distant from the Strix species. S. ulula exhibited the greatest genetic distance from the other species while A. noctua was in a mid-range position (Figure 2). B. scandiacus and B. bubo did not overlap, but the genetic distance of the two species is low. A principal component analysis was then carried out, revealing that, as expected, the three species of the genus Strix (S. uralensis, S. nebulosa, S. aluco) overlapped in the PCA graphic ( Figure 2). A. otus and O. scops showed a low distance from each other and resulted not very distant from the Strix species. S. ulula exhibited the greatest genetic distance from the other species while A. noctua was in a mid-range position (Figure 2). B. scandiacus and B. bubo did not overlap, but the genetic distance of the two species is low.

Evaluation of the STR Potential Panel for Parentage Analysis in Family Groups
The paternity test yielded inconsistent results and different reliability values depending on the species (Tables 5 and 6). In Colony, the correct parent pair was assigned with the maximum probability value (1.00) but decreased in S. nebulosa (0.992 and 0.907), S. uralensis (0.997), and S. ulula (0.983 and 0.753) when the computation has been limited to association with putative mother or father. In two individuals, respectively, of S. uralensis (S_ur6) and S. ulula (S_ul3), the probability of assignment to the right parent is reduced to 0.500.
In Cervus, using the LOD computation, all the individuals have been correctly associated with the right parents with trio confidence values higher than 95%, except for S_ul3 that has been associated with the wrong father (S_ul5 instead of S_ul2). In S. uralensis (S_ur6), the assignment to the right mother did not reach a significant value. Using the Delta calculation, the probability values decreased in S. nebulosa, S. uralensis and S. ulula (see Table 6 for more details). Again, S_ul3 was associated with the wrong father.

Discussion
International wildlife trafficking today is recognized as one of the largest organized transnational crimes [38], which equals the trafficking of drugs, arms and humans (World Wildlife Report, United Nations: Office on Drugs and Crime).
Genotyping assays through microsatellites are a quick, informative and low-cost approach for linking items of evidence to crimes in forensic investigations [39]. Since microsatellites have high mutation rates, they are used within the context of monitoring illegal wildlife trade primarily to identify individuals, assign individuals to specific populations or for relatedness testing. Oklander et al. [39] generated a multilocus microsatellite genotype reference database of the black and gold howler monkey (Alouatta caraya), a neotropical primate threatened by habitat loss and capture for illegal trade in Argentina, to assign confiscated individuals to localities of origin, illustrating the applicability of genotype databases for inferring hotspots of illegal capture. Potoczniak et al. [39] developed a STR genotyping assay able to associate a biological sample to Asian elephant (Elephas maximus) or African elephant (Loxodonta africana). Fitzsimmons et al. [40] designed a panel of 26 microsatellite loci for Crocodylus spp. to answer questions regarding population assignment, mating system and geneflow. Jan and Fumagalli [41] isolated DNA microsatellite markers in seven parrot species threatened with extinction and subjected to illegal trafficking, characterized a total of 106 polymorphic microsatellite markers and tested them for individual identification and parental analyses. Mucci and colleagues [42] developed a panel of 16 de novo sequenced microsatellites (STRs) for Testudo graeca and tested its effectiveness for parentage analysis in two other species of endangered tortoises, T. hermanni and T. marginata.
Given the utility of STR-based approaches to answer questions related to wildlife crime investigations as forensic genetics, efforts should be devoted to the characterization of microsatellite primers for species threatened by illegal trafficking. However, while in human forensics, the selection of around 20 core STR loci allowed the standardization around the globe for human identity testing [43,44], accomplishing this same achievement is more challenging for the hyper-diverse animal assemblage encountered in wildlife forensics [38].
The lack of species-specific molecular markers or their inadequate representation in genetic databases is major limitation in wildlife forensics [45].
Cross-amplification is a widely used approach permitting to avoid investing time and money in the development of new markers, and many studies have proven the efficiency of microsatellite loci developed in closely related species [42,[46][47][48][49].
The first aim of this study was to test for the presence of a minimum number of microsatellite loci reliable for individual identification and parentage analysis for forensic purposes in 11 species of Strigiformes listed in the CITES Appendix II and regularly traded in the Italian national market.
Literature data show that cross-species transferability is unevenly distributed across taxa. Barbará et al. [50] reviewed 64 primer notes and found more than 40% transfer success in mammals, more than 25% in fishes and more than 10% in birds. Our results permitted to define a unique panel of 12 out of 32 highly polymorphic microsatellites (37%), able to identify individuals in nine species of two subfamilies of Strigiformes (Striginae and Surninae) belonging to Family Strigidae (A. otus, A. noctua, B. bubo, B. scandiacus, O. scops, S. aluco, S. nebulosa, S. uralensis, S. ulula).
The test was non-efficient in individuals of T. alba in which only six markers yield a positive result but four of them resulted monomorphic. Such results could rely on the fact that T. alba is the only species of this study belonging to a different family [51]. The two monophyletic families of Strigiformes, Tytonidae and Strigidae, diverged in the middle of the Eocene [52]. The phylogenetic divergence of Tytonidae from Strigidae could be justified by the retrieved inefficiency of markers panel tested on species. Since only two subfamilies represent Tytonidae family with a single genus each [53]-Tytoninae with Tyto and Phodilinae with Phodilus-it is not possible to verify the discriminant power of tested markers on other species of the same taxon.
The common barn owl T. alba is one of the most cosmopolitan species and represents a taxon-rich species complex with several subspecies [54]. The interest in its conservation status clears the need for integrating this set with more polymorphic loci for this species. Twenty-one microsatellite loci already isolated and characterized in T. alba [26] could be tested for cross-amplification in CITES species to implement this panel.
The principal component analysis results reflect the most recent knowledge about taxonomy and systematics of Strigiformes [53].
A slightly better but similar result was obtained for G. passerinum, for which only 4 markers were polymorphic. The high number of monomorphic loci in this species could be due to several cumulative factors: the low number of related individuals analyzed, the high inbreeding of captive-bred individuals, or the inefficiency of the selected markers for this species. Unlike T. alba, microsatellites were developed only for a species of the same genera: Glaucidium brasilianum [23].
Evaluating the STR panel potential in identifying family groups, we found that in two out of the nine species analyzed (S. uralensis and S. ulula), the probability value associated with parent pairs was reduced. The PID _sib almost reached the 0.001 threshold value in S. ulula, S. nebulosa and S. uralensis (see Table 4: Probability of identity values, PID and PID_ sibs ). However, PID and PID _sib values are subjected to bias due to the low number of tested individuals per species and bottleneck in captive breeding facilities. For these three species, it could be helpful to increase the number of samples and markers to obtain a more confident individual identification and association to parent pairs.
Besides reducing time and costs when adopting a cross-amplification approach, another advantage of using a shared panel among several species is the possibility of comparing genetic variability values among species. Nevertheless, even if microsatellites loci are very useful genetic markers in studying the mating system, population genetics, and conservation of owls, many studies focus on species belonging to the same genus. Dial et al. [55] screened many markers developed in strigids but found only four polymorphic pairs in the great horned owl (Bubo virginianus), short-eared owl (Asio flammeus) and the snowy owl (Bubo scandiacus). In addition, another eight reliably amplified polymorphic fragments only in the great horned owl, eleven in the short-eared owl, and ten in the snowy owl.
Hsu et al. [22] developed six new microsatellite markers containing tetranucleotide repeat motifs (GATA/CTAT) for Lanyu scops owl (Otus elegans botelensis). They tested them, and additional further microsatellite primer pairs previously developed from O. elegans on four other species of owls (O. lettia, O. spilocephalus, O. scops and Ninox scutulata). Data obtained showed a reduced degree of polymorphism with most of the loci resulting polymorphic in the three Otus owls but only five in N. scutulata.
Our study gives a valuable tool to implement research involving most Strigidae threatened by illegal trafficking, habitat loss and fragmentation in Italy and other countries.
Delport et al. [56] tested 19 loci originally developed for Vidua and Geospiza for cross-amplification in Nesospiza buntings. They detected a degree of polymorphism and heterozygosity lower in loci developed for Vidua than those explicitly developed for Nesospiza. These data demonstrate that microsatellite markers isolated in the reference species are frequently less variable in related species. Moreover, cross-species amplification is usually limited to the loci that were found polymorphic in the referent species. In this study, we selected from the literature the most variable loci in each reference species: most polymorphic loci were discarded only when they were found monomorphic or did not give amplification products in all the target species, with only three exceptions as mentioned in Results (Paragraph ii), because of their utility in other species. We are aware that this a priori selection can cause an ascertainment bias, as suggested by Delport et al. [56]; however, we found high levels of polymorphisms in nine out 11 analyzed species.
Variability indices found in different species using this panel were not discordant from the ones found in the natural populations with different markers: Pellegrino et al. [57] found in the A. noctua European populations an average and an effective number of alleles = 5.6 and 3.5, respectively, and observed and expected heterozygosities equal to 0.59 and 0.61, respectively. Pertoldi et al. [58] found lower alleles in the Danish population (effective number of alleles = 2.8; Ho = 0.51 and He = 0.60), probably caused by a population bottleneck in the last decades.
Microsatellite loci represent reliable molecular markers to describe genetic variability or its drastic reduction, as demonstrated by Macías-Duarte et al. [59] that found different values in three different populations of Athene cunicularia, with the average number of alleles varying from 2.7 in Clarion Islands to 5.1 in Florida and 22.5 in Western North America.
Though their high polymorphism makes them adequate for conservation and forensic genetics purposes, the main difficulties are represented by comparing samples between laboratories.
This problem that has been resolved through the exchange of reference samples has been recently fixed by the set-up of an allelic ladder [60]. According to these authors, we constructed an allelic ladder for each locus to standardize a protocol between laboratories for conservation and forensic purposes.
Comparability between laboratories is now also possible thanks to high-throughput sequencing (HTS) technologies [61]. This application that has been developed and used in human forensics [62][63][64] has been already applied also in conservation genetics [65][66][67][68]. This method permits the sequencing of STRs, allowing the identification of the correct number of repeats. The possibility of multiplexing several dozen of markers from a single individual will allow cost and time reduction. De Barba et al. [69] used this protocol in the study of a brown bear population yielded reliable results of parentage analysis also from low quality DNA, confirming a broader application in conservation genetics and forensics.

Conclusions
The newly optimized 12 loci can provide the authorities with the ability to investigate suspected smugglers and false parental claims or establish a link between evidence and an individual (e.g., identifying when a bird is illegally transferred between different breeding facilities). In addition to this, the selected markers can be used to assess relatedness among captive-bred individuals, which can be crucial in designing optimal breeding protocols to avoid genetic variability loss and minimize inbreeding. Moreover, these markers are useful for applications in research projects on wild populations where individual genotype identification allows to describe the mating system, genetic variability, population structure and geneflow among nine different species of Strigiformes with conservation concern.