Next Article in Journal
CMAX3: A Robust Statistical Test for Genetic Association Accounting for Covariates
Next Article in Special Issue
On the Identification of Body Fluids and Tissues: A Crucial Link in the Investigation and Solution of Crime
Previous Article in Journal
Effect of High Hydrostatic Pressure on Stress-Related dnaK, hrcA, and ctsR Expression Patterns in Selected Lactobacilli Strains
Previous Article in Special Issue
Precision DNA Mixture Interpretation with Single-Cell Profiling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes

1
Conservation Genetics Area, Institute for Environmental Protection and Research, Via Ca’ Fornacetta 9, Ozzano dell’Emilia, 40064 Bologna, Italy
2
Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, Piazzale delle Cascine 18, 50144 Florence, Italy
3
Department of Biology, University of Florence, Via Madonna del Piano 6, Sesto Fiorentino, 50019 Florence, Italy
*
Author to whom correspondence should be addressed.
Genes 2021, 12(11), 1721; https://doi.org/10.3390/genes12111721
Submission received: 5 October 2021 / Revised: 20 October 2021 / Accepted: 25 October 2021 / Published: 28 October 2021
(This article belongs to the Special Issue Advances in Forensic Genetics)

Abstract

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

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

2. 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).
(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 MgCl2, 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′.
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.

3. Results

3.1. Preliminary Screening of STR Polymorphic Loci

Thirty-two microsatellite markers were chosen to evaluate their cross-species amplification potential on three species (Athene noctua, Strix aluco, Bubo bubo) (Table 1), belonging to two different sub-families, Striginae and Surninae, and three tribes (Strigini, Bubonini, Athenini) respectively [14].
Two unrelated individuals per species were selected from the CITES sample database (managed since 1995 by the Italian Institute of Environmental Protection and Research –ISPRA-on behalf of the Italian Environmental Ministry). Data obtained revealed that twenty-one out of 32 markers used (15a6, Age5, Bb101, Bb126, BbuS027, Bbus116, Calex05, Fepo42, Oe054, Oe128, Oe129, Oe142, Oe149, Oe321, Oe53, Oe81, Oe84, SneD113, SneD202, SneD218, Tgu06) gave amplicons in A. noctua, S. aluco and B. bubo DNA samples. The remaining 11 microsatellites were discarded because they either did not give any amplification product or led to unreliable PCR amplicons.

3.2. 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.
Data obtained revealed that 12 of them resulted polymorphic in at least 2 species (15a6, Fepo42, Oe053, Oe054, Oe128, Oe129, Oe142, Oe149, Oe321, SneD113, SneD218, Tgu06) and were retained for further analyses. Fepo42 and Oe054 resulted monomorphic in S. aluco, and Oe129 did not give any amplification product in B. bubo but were retained.

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

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

3.5. 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 (Table 5 and Table 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.

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

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

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12111721/s1, Figure S1: allelic ladders; Figure S2: individual genotypes.

Author Contributions

Conceptualization, N.M., G.A., R.F., M.Z.; methodology: N.M., C.M., C.G., P.G.; Validation: P.G., N.M., C.G.; formal analysis, P.G., L.N., A.P.; writing—original draft preparation, N.M., P.G.; writing—review and editing, P.G., C.M., A.P., C.G., M.Z., R.F., G.A.; supervision, N.M., R.F.; project administration, N.M., C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Ethical review and approval were waived for this study since they are not officially required by national regulations. In any case, the study abided by local policies and regulations. The methods applied in this study were carried out in strict accordance with the Nature and Sea Protection Directorate (DPNM, which belongs to the Ministry for the Ecological Transition) and the Italian Institute for Environmental Protection and Research (ISPRA). The breeders authorized the feather’s sampling of their own animals. The handling of animals and feather’s samplings were carried out by qualified veterinary personnel in accordance with animal welfare procedures and in the presence and under the jurisdiction of the Italian forest police units (Carabinieri Command). The Carabinieri Command exercises management, coordination and control functions of the CITES operational units, in charge for the application of CITES directives. The CITES units checked that the breeders fulfilled requirements of the art. 54 of Reg. (CE) n. 865/2006 for investigations of birth in captivity.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the Italian Ministry for Ecological Transition (MITE), particularly the Italian CITES Management Authority, which since 1997 granted support to the Unit for Conservation Genetics (BIO-CGE) at ISPRA, and allowed us to develop molecular procedures for improving wildlife forensics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anderson, S.C.; Elsen, P.R.; Hughes, B.B.; Tonietto, R.K.; Bletz, M.C.; Gill, D.A.; Holgerson, M.A.; Kuebbing, S.E.; McDonough MacKenzie, C.; Meek, M.H.; et al. Trends in ecology and conservation over eight decades. Front. Ecol. Environ. 2021, 19, 274–282. [Google Scholar] [CrossRef]
  2. Bishop, J.M.; Leslie, A.J.; Bourquin, S.L.; O’Ryan, C. Reduced effective population size in an overexploited population of the Nile crocodile (Crocodylus niloticus). Biol. Conserv. 2009, 142, 2335–2341. [Google Scholar] [CrossRef]
  3. Betancourth-Cundar, M.; Palacios-Rodríguez, P.; Mejía-Vargas, D.; Paz, A.; Amézquita, A. Genetic differentiation and overexploitation history of the critically endangered Lehmann’s Poison Frog: Oophaga lehmanni. Conserv. Genet. 2020, 21, 453–465. [Google Scholar] [CrossRef]
  4. Righi, T.; Splendiani, A.; Fioravanti, T.; Casoni, E.; Gioacchini, G.; Carnevali, O.; Barucchi, V.C. Loss of mitochondrial genetic diversity in overexploited mediterranean swordfish (Xiphias gladius, 1759) population. Diversity 2020, 12, 170. [Google Scholar] [CrossRef]
  5. Baker, S.E.; Cain, R.; Kesteren, F.V.; Zommers, Z.A.; D’Cruze, N.; MacDonald, D.W. Rough trade: Animal welfare in the global wildlife trade. BioScience 2013, 63, 928–938. [Google Scholar] [CrossRef] [Green Version]
  6. Tingley, M.W.; Harris, J.B.C.; Hua, F.; Wilcove, D.S. The pet trade’s role in defaunation. Science 2017, 356, 1–916. [Google Scholar] [CrossRef] [PubMed]
  7. Panter, C.T.; Atkinson, E.D.; White, R.L. Retraction: Quantifying the global legal trade in live CITES-listed raptors and owls for commercial purposes over a 40-year period. Avocetta 2019, 43, 23–36. [Google Scholar] [CrossRef]
  8. Vall-Llosera, M.; Su, S. Trends and characteristics of imports of live CITES-listed bird species into Japan. Ibis 2019, 161, 590–604. [Google Scholar] [CrossRef]
  9. Brochet, A.L.; Van Den Bossche, W.; Jbour, S.; Ndang’Ang’A, P.K.; Jones, V.R.; Abdou, W.A.L.I.; Al-Hmoud, A.R.; Asswad, N.G.; Atienza, J.C.; Atrash, I.; et al. Preliminary assessment of the scope and scale of illegal killing and taking of birds in the Mediterranean. Bird Conserv. Int. 2016, 26, 1–28. [Google Scholar] [CrossRef] [Green Version]
  10. Domínguez, M.; Tiedemann, R.; Reboreda, J.C.; Segura, L.; Tittarelli, F.; Mahler, B. Genetic structure reveals management units for the yellow cardinal (Gubernatrix cristata), endangered by habitat loss and illegal trapping. Conserv. Genet. 2017, 18, 1131–1140. [Google Scholar] [CrossRef]
  11. Harris, J.B.C.; Green, J.M.H.; Prawiradilaga, D.M.; Giam, X.; Giyanto; Hikmatullah, D.; Putra, C.A.; Wilcove, D.S. Using market data and expert opinion to identify overexploited species in the wild bird trade. Biol. Conserv. 2015, 187, 51–60. [Google Scholar] [CrossRef]
  12. Evans, S.R.; Sheldon, B.C. Interspecific patterns of genetic diversity in birds: Correlations with extinction risk. Conserv. Biol. 2008, 22, 1016–1025. [Google Scholar] [CrossRef] [PubMed]
  13. Enríquez, P.L.; Eisermann, K.; Mikkola, H.; Motta-Junior, J.C. A review of the systematics of Neotropical owls (Strigiformes). Neotrop. Owls Divers. Conserv. 2017, 1–670. [Google Scholar] [CrossRef]
  14. Wink, M.; El-Sayed, A.A.; Sauer-Gürth, H.; Gonzalez, J. Molecular phylogeny of owls (Strigiformes) inferred from DNA sequences of the mitochondrial cytochrome b and the nuclear RAG-1 gene. Ardea 2009, 97, 581–591. [Google Scholar] [CrossRef] [Green Version]
  15. Abdel Rabou, A.F.N. On the Owls (Order Strigiformes) Inhabiting the Gaza Strip—Palestine. JOJ Wildl. Biodivers. 2020, 3, 1–11. [Google Scholar] [CrossRef]
  16. Wan, H.Y.; Ganey, J.L.; Vojta, C.D.; Cushman, S.A. Managing emerging threats to spotted owls. J. Wildl. Manag. 2018, 82, 682–697. [Google Scholar] [CrossRef]
  17. Li, L.; Jiang, Z. International trade of CITES listed bird species in China. PLoS ONE 2014, 9, e85012. [Google Scholar] [CrossRef]
  18. Ribeiro, J.; Reino, L.; Schindler, S.; Strubbe, D.; Vall-llosera, M.; Araújo, M.B.; Capinha, C.; Carrete, M.; Mazzoni, S.; Monteiro, M.; et al. Trends in legal and illegal trade of wild birds: A global assessment based on expert knowledge. Biodivers. Conserv. 2019, 28, 3343–3369. [Google Scholar] [CrossRef]
  19. Nijman, V.; Nekaris, K.A.I. The Harry Potter effect: The rise in trade of owls as pets in Java and Bali, Indonesia. Glob. Ecol. Conserv. 2017, 11, 84–94. [Google Scholar] [CrossRef]
  20. Siriwat, P.; Nekaris, K.A.I.; Nijman, V. Digital media and the modern-day pet trade: A test of the ‘Harry Potter effect‘ and the owl trade in Thailand. Endanger. Species Res. 2020, 41, 7–16. [Google Scholar] [CrossRef] [Green Version]
  21. Isaksson, M.; Tegelström, H. Characterization of polymorphic microsatellite markers in a captive population of the eagle owl (Bubo bubo) used for supportive breeding. Mol. Ecol. Notes 2002, 2, 91–93. [Google Scholar] [CrossRef]
  22. Hsu, Y.C.; Li, S.H.; Lin, Y.S.; Severinghaus, L.L. Microsatellite loci from Lanyu scops owl (Otus elegans botelensis) and their cross-species application in four species of strigidae. Conserv. Genet. 2006, 7, 161–165. [Google Scholar] [CrossRef]
  23. Proudfoot, G.; Honeycutt, R.; Douglas Slack, R. Development and characterization of microsatellite DNA primers for ferruginous pygmy-owls (Glaucidium brasilianum). Mol. Ecol. Notes 2005, 5, 90–92. [Google Scholar] [CrossRef]
  24. Hull, J.M.; Keane, J.J.; Tell, L.A.; Ernest, H.B. Development of 37 microsatellite loci for the great gray owl (Strix nebulosa) and other Strix spp. owls. Conserv. Genet. 2008, 9, 1357–1361. [Google Scholar] [CrossRef]
  25. Klein, A.; Horsburgh, G.J.; Kupper, C.; Major, A.; Lee, P.L.M.; Hoffmann, G.; Matics, R.; Dawson, D.A. Microsatellite markers characterized in the barn owl (Tyto alba) and high utility in the other owls (Strigiformes: AVES). Mol. Ecol. Resour. 2009, 9, 1513–1519. [Google Scholar] [CrossRef]
  26. Burri, R.; Antoniazza, S.; Siverio, F.; Klein, Á.; Roulin, A.; Fumagalli, L. Isolation and characterization of 21 microsatellite markers in the barn owl (Tyto alba). Mol. Ecol. Resour. 2008, 8, 977–979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Hogan, F.E.; Cooke, R.; Norman, J.A. Reverse ascertainment bias in microsatellite allelic diversity in owls (Aves, Strigiformes). Conserv. Genet. 2009, 10, 635–638. [Google Scholar] [CrossRef]
  28. Kleven, O.; Dawson, D.A.; Gjershaug, J.O.; Horsburgh, G.J.; Jacobsen, K.O.; Wabakken, P. Isolation, characterization and predicted genome locations of Eurasian eagle-owl (Bubo bubo) microsatellite loci. Conserv. Genet. Resour. 2013, 5, 723–727. [Google Scholar] [CrossRef] [Green Version]
  29. Hsu, Y.C.; Severinghaus, L.L.; Lin, Y.S.; Li, S.H. Isolation and characterization of microsatellite DNA markers from the Lanyu scops owl (Otus elegans botelensis). Mol. Ecol. Notes 2003, 3, 595–597. [Google Scholar] [CrossRef] [Green Version]
  30. Thode, A.B.; Maltbie, M.; Hansen, L.A.; Green, L.D.; Longmire, J.L. Microsatellite markers for the Mexican spotted owl (Strix occidentalis lucida). Mol. Ecol. Notes 2002, 2, 446–448. [Google Scholar] [CrossRef]
  31. Topinka, J.R.; May, B. Development of polymorphic microsatellite loci in the northern goshawk (Accipiter gentilis) and cross-amplification in other raptor species. Conserv.Genet. 2004, 5, 861–864. [Google Scholar] [CrossRef]
  32. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research--an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Valière, N. GIMLET: A computer program for analysing genetic individual identification data. Mol. Ecol. Notes 2002, 2, 377–379. [Google Scholar] [CrossRef]
  34. Jones, O.R.; Wang, J. COLONY: A program for parentage and sibship inference from multilocus genotype data. Mol. Ecol. Resour. 2010, 10, 551–555. [Google Scholar] [CrossRef]
  35. Marshall, T.C.; Slate, J.; Kruuk, L.E.B.; Pemberton, J.M. Statistical confidence for likelihood-based paternity inference in natural populations. Mol. Ecol. 1998, 7, 639–655. [Google Scholar] [CrossRef] [Green Version]
  36. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 2007, 16, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  37. Belkhir, K.; Borsa, P.; Chikhi, L.; Raufaste, N.; Bonhomme, F. Genetix 4.05, Logiciel Sous Windows Tm Pour La Génétique Des Populations; Laboratoire Genome, Populations, Interactions, CNRS UMR 5000, Universite de Montpellier II: Montpellier, France, 1996–2004.
  38. Smart, U.; Cihlar, J.C.; Budowle, B. International Wildlife Trafficking: A perspective on the challenges and potential forensic genetics solutions. Forensic Sci. Int. Genet. 2021, 54, 102551. [Google Scholar] [CrossRef] [PubMed]
  39. Potoczniak, M.J.; Chermak, M.; Quarino, L.; Tobe, S.S.; Conte, J. Development of a multiplex, PCR-based genotyping assay for African and Asian elephants for forensic purposes. Int. J. Leg. Med. 2020, 134, 55–62. [Google Scholar] [CrossRef]
  40. Fitzsimmons, N.N.; Tanksley, S.; Forstner, M.R.J.; Louis, E.E.; Daglish, R.; Gratten, J.; Davis, S. Microsatellite markers for Crocodylus: New genetic tools for population genetics, mating system studies and forensics. Crocodilian Biol. Evol. 2001, 51–57. [Google Scholar]
  41. Jan, C.; Fumagalli, L. Polymorphic DNA microsatellite markers for forensic individual identification and parentage analyses of seven threatened species of parrots (family Psittacidae). PeerJ 2016, 2016, 1–7. [Google Scholar] [CrossRef] [Green Version]
  42. Mucci, N.; Giangregorio, P.; Cirasella, L.; Isani, G.; Mengoni, C. A new STR panel for parentage analysis in endangered tortoises. Conserv. Genet. Resour. 2020, 12, 67–75. [Google Scholar] [CrossRef]
  43. Zhou, Z.; Shao, C.; Xie, J.; Xu, H.; Liu, Y.; Zhou, Y.; Liu, Z.; Zhao, Z.; Tang, Q.; Sun, K. Genetic polymorphism and phylogenetic analyses of 21 non-CODIS STR loci in a Chinese Han population from Shanghai. Mol. Genet. Genom. Med. 2020, 8, 1–8. [Google Scholar] [CrossRef] [PubMed]
  44. Kumar, A.; Kumar, R.; Kumawat, R.K.; Shrivastava, P.; Chaubey, G. Genomic diversity at 22 STR loci (extended CODIS STR) in the population of Rajasthan, India. Gene Rep. 2021, 23, 101150. [Google Scholar] [CrossRef]
  45. Rocco, F.D.; Anello, M. The Use of Forensic DNA on the Conservation of Neotropical Mammals. In Molecular Ecology and Conservation Genetics of Neotropical Mammals; Springer Nature Switzerland AG: Cham, Switzerland, 2021; pp. 85–98. [Google Scholar]
  46. Gebhardt, K.J.; Waits, L.P. Cross-species amplification and optimization of microsatellite markers for use in six Neotropical parrots. Mol. Ecol. Resour. 2008, 8, 835–839. [Google Scholar] [CrossRef] [PubMed]
  47. Dawson, D.A.; Horsburgh, G.J.; Küpper, C.; Stewart, I.R.K.; Ball, A.D.; Durrant, K.L.; Hansson, B.; Bacon, I.; Bird, S.; Klein, Á.; et al. New methods to identify conserved microsatellite loci and develop primer sets of high cross-species utility—As demonstrated for birds. Mol. Ecol. Resour. 2010, 10, 475–494. [Google Scholar] [CrossRef] [PubMed]
  48. Maulidi, A.; Fatchiyah, F.; Hamidy, A.; Kurniawan, N. Microsatellite Marker for Cross-Species Amplification: Study Case for Indonesian Sundaland Python (Serpentes: Pythonidae). J. Exp. Life Sci. 2018, 8, 61–65. [Google Scholar] [CrossRef] [Green Version]
  49. Du Toit, Z.; Dalton, D.L.; du Plessis, M.; Jansen, R.; Grobler, J.P.; Kotzé, A. Isolation and characterization of 30 STRs in Temminck’s ground pangolin (Smutsia temminckii) and potential for cross amplification in other African species. J. Genet. 2020, 99. [Google Scholar] [CrossRef]
  50. Barbará, T.; Palma-Silva, C.; Paggi, G.M.; Bered, F.; Fay, M.F.; Lexer, C. Cross-species transfer of nuclear microsatellite markers: Potential and limitations. Mol. Ecol. 2007, 16, 3759–3767. [Google Scholar] [CrossRef]
  51. König, C.; Weick, F. Owls of the World; AandC Black Publishers Ltd.: London, UK, 2008. [Google Scholar]
  52. Prum, R.O.; Berv, J.S.; Dornburg, A.; Field, D.J.; Townsend, J.P.; Lemmon, E.M.; Lemmon, A.R. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature 2015, 526, 569–573. [Google Scholar] [CrossRef] [PubMed]
  53. Wink, M.; Sauer-Gurth, H. Molecular taxonomy and systematics of owls (Strigiformes)—An update. Airo 2021, 29, 487–500. [Google Scholar]
  54. Aliabadian, M.; Alaei-Kakhki, N.; Mirshamsi, O.; Nijman, V.; Roulin, A. Phylogeny, biogeography, and diversification of barn owls (Aves: Strigiformes). Biol. J. Linn. Soc. 2016, 119, 904–918. [Google Scholar] [CrossRef]
  55. Dial, C.R.; Talbot, S.L.; Sage, G.K.; Seidensticker, M.T.; Holt, D.W. Cross-species Amplification of Microsatellite Markers in the Great Horned Owl Bubo virginianus, Short-eared Owl Asio flammeus and Snowy Owl B. scandiacus for Use in Population Genetics, Individual Identification and Parentage Studies. J. Yamashina Inst. Ornithol. 2012, 44, 1–12. [Google Scholar] [CrossRef]
  56. Delport, W.; Grant, T.J.; Ryan, P.G.; Bloomer, P. Ten microsatellite loci for evolutionary research on Nesospiza buntings. Mol. Ecol. Notes 2006, 6, 1180–1183. [Google Scholar] [CrossRef]
  57. Pellegrino, I.; Negri, A.; Boano, G.; Cucco, M.; Kristensen, T.N.; Pertoldi, C.; Randi, E.; Šálek, M.; Mucci, N. Evidence for strong genetic structure in European populations of the little owl Athene noctua. J. Avian Biol. 2015, 46, 462–475. [Google Scholar] [CrossRef]
  58. Pertoldi, C.; Pellegrino, I.; Cucco, M.; Mucci, N.; Randi, E.; Laursen, J.T.; Sunde, P.; Loeschcke, V.; Kristensen, T.N.; du Toit, Z.; et al. Genetic consequences of population decline in the Danish population of the little owl (Athene noctua). Mol. Ecol. 2020, 14, 358–367. [Google Scholar] [CrossRef] [Green Version]
  59. Macías-Duarte, A.; Conway, C.J.; Holroyd, G.L.; Valdez-Gómez, H.E.; Culver, M. Genetic Variation among Island and Continental Populations of Burrowing Owl (Athene cunicularia) Subspecies in North America. J. Raptor Res. 2019, 53, 127–133. [Google Scholar] [CrossRef] [Green Version]
  60. Biello, R.; Zampiglia, M.; Corti, C.; Deli, G.; Biaggini, M.; Crestanello, B.; Delaugerre, M.; Di Tizio, L.; Leonetti Francesco, L.; Stefano, C.; et al. Mapping the geographic origin of captive and confiscated Hermann’s tortoises: A genetic toolkit for conservation and forensic analyses. Forensic Sci. Int. Genet. 2021, 51, 102447. [Google Scholar] [CrossRef] [PubMed]
  61. Glenn, T.C. Field guide to next-generation DNA sequencers. Mol. Ecol. Resour. 2011, 11, 759–769. [Google Scholar] [CrossRef]
  62. Fordyce, S.L.; Ávila-Arcos, M.C.; Rockenbauer, E.; Børsting, C.; Frank-Hansen, R.; Petersen, F.T.; Willerslev, E.; Hansen, A.J.; Morling, N.; Gilbert, T.P. High-throughput sequencing of core STR loci for forensic genetic investigations using the Roche Genome Sequencer FLX platform. BioTechniques 2011, 51, 127–133. [Google Scholar] [CrossRef] [PubMed]
  63. Bornman, D.M.; Hester, M.E.; Schuetter, J.M.; Kasoji, M.D.; Minard-Smith, A.; Barden, C.A.; Faith, S.A. Short-read, high-throughput sequencing technology for STR genotyping. BioTechniques. Rapid Dispatches 2012, 23, 1–7. [Google Scholar] [CrossRef]
  64. Van Neste, C.; Van Nieuwerburgh, F.; Van Hoofstat, D.; Deforce, D. Forensic STR analysis using massive parallel sequencing. Forensic Sci. Int. Genet. 2012, 6, 810–818. [Google Scholar] [CrossRef]
  65. Stoeckle, B.C.; Theuerkauf, J.; Rouys, S.; Gula, R.; Lorenzo, A.; Lambert, C.; Kaeser, T.; Kuehn, R. Identification of polymorphic microsatellite loci for the endangered Kagu (Rhynochetos jubatus) by high-throughput sequencing. J. Ornithol. 2012, 153, 249–253. [Google Scholar] [CrossRef] [Green Version]
  66. Pimentel, J.S.M.; Carmo, A.O.; Rosse, I.C.; Martins, A.P.V.; Ludwig, S.; Facchin, S.; Pereira, A.H.; Brandão-Dias, P.F.P.; Abreu, N.L.; Kalapothakis, E. High-throughput sequencing strategy for microsatellite genotyping using neotropical fish as a model. Front. Genet. 2018, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  67. Qi, W.H.; Lu, T.; Zheng, C.L.; Jiang, X.M.; Jie, H.; Zhang, X.Y.; Yue, B.S.; Zhao, G.J. Distribution patterns of microsatellites and development of its marker in different genomic regions of forest musk deer genome based on high throughput sequencing. Aging 2020, 12, 4445–4462. [Google Scholar] [CrossRef] [PubMed]
  68. Song, C.; Feng, Z.; Li, C.; Sun, Z.; Gao, T.; Song, N.; Liu, L. Profile and development of microsatellite primers for Acanthogobius ommaturus based on high-throughput sequencing technology. J. Oceanol. Limnol. 2020, 38, 1880–1890. [Google Scholar] [CrossRef]
  69. De Barba, M.; Miquel, C.; Lobréaux, S.; Quenette, P.Y.; Swenson, J.E.; Taberlet, P. High-throughput microsatellite genotyping in ecology: Improved accuracy, efficiency, standardization and success with low-quantity and degraded DNA. Mol. Ecol. Resour. 2017, 17, 492–507. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graph of the PID_sib trend in the species.
Figure 1. Graph of the PID_sib trend in the species.
Genes 12 01721 g001
Figure 2. Principal Component Analysis plot obtained using Genetix and visualized in Excel. Strix sp. overlapped while Bubo bubo and Bubo scandiacus plotted very closed each other. Surnia ulula individuals diverged consistently from other species. Ovals in different colors represent the distribution of individuals belonging to six tribes of Striginae and Surninae, in the blue and green squares, respectively.
Figure 2. Principal Component Analysis plot obtained using Genetix and visualized in Excel. Strix sp. overlapped while Bubo bubo and Bubo scandiacus plotted very closed each other. Surnia ulula individuals diverged consistently from other species. Ovals in different colors represent the distribution of individuals belonging to six tribes of Striginae and Surninae, in the blue and green squares, respectively.
Genes 12 01721 g002
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.
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.
FamilySubfamilyTribeSpeciesIDGenderRelationshipsFamily Group
StrigidaeStriginaeStriginiStrix alucoS_al 1MaleFather1
S_al 2FemaleMother
S_al 3unknownOffspring
S_al 4MaleFather2
S_al 5FemaleMother
S_al 6unknownOffspring
S_al 7MaleFather3
S_al 8FemaleMother
S_al 9unknownOffspring
BuboniniBubo buboB_bu 1MaleFather4
B_bu 2FemaleMother
B_bu 3unknownOffspring
B_bu 4MaleFather5
B_bu 5FemaleMother
B_bu 6unknownOffspring
B_bu 7MaleFather6
B_bu 8FemaleMother
B_bu 9unknownOffspring
SurninaeAtheniniAthene noctuaA_no 1MaleFather7
A_no 2FemaleMother
A_no 3unknownOffspring
A_no 4MaleFather8
A_no 5FemaleMother
A_no 6unknownOffspring
A_no 7MaleFather9
A_no 8FemaleMother
A_no 9unknownOffspring
StriginaeStriginiStrix nebulosaS_ne 1MaleFather10
S_ne 2FemaleMother
S_ne 3unknownOffspring
S_ne 4MaleFather11
S_ne 5FemaleMother
S_ne 6unknownOffspring
S_ne 7MaleFather12
S_ne 8FemaleMother
S_ne 9unknownOffspring
Strix uralensisS_ur 1MaleFather13
S_ur 2FemaleMother
S_ur 3unknownOffspring
S_ur 4MaleFather14
S_ur 5FemaleMother
S_ur 6unknownOffspring
BuboniniBubo scandiacusB_sc 1MaleFather15
B_sc 2FemaleMother
B_sc 3unknownOffspring
B_sc 4MaleFather16
B_sc 5FemaleMother
B_sc 6unknownOffspring
B_sc 7MaleFather17
B_sc 8FemaleMother
B_sc 9unknownOffspring
OtiniOtus scopsO_sc 1MaleFather18
O_sc 2FemaleMother
O_sc 3unknownOffspring
O_sc 4MaleFather19
O_sc 5FemaleMother
O_sc 6unknownOffspring
O_sc 7MaleFather20
O_sc 8FemaleMother
O_sc 9unknownOffspring
AsioniniAsio otusA_ot 1MaleFather21
A_ot 2FemaleMother
A_ot 3unknownOffspring
A_ot 4MaleFather22
A_ot 5FemaleMother
A_ot 6unknownOffspring
SurninaeSurniniSurnia ululaS_ul 1MaleFather23
S_ul 2FemaleMother
S_ul 3unknownOffspring
S_ul 4MaleFather24
S_ul 5FemaleMother
S_ul 6unknownOffspring
S_ul 7MaleFather25
S_ul 8FemaleMother
S_ul 9unknownOffspring
Glaucidium passerinumG_pa 1MaleFather26
G_pa 2FemaleMother
G_pa 3unknownOffspring
G_pa 4MaleFather27
G_pa 5FemaleMother
G_pa 6unknownOffspring
TytonidaeTytoninae Tyto albaT_al 1MaleFather28
T_al 2FemaleMother
T_al 3unknownOffspring
Table 2. Loci chosen from literature and tested for the cross-species amplification study.
Table 2. Loci chosen from literature and tested for the cross-species amplification study.
Primer IDForward Primer SequenceReverse Primer SequenceOriginal SpeciesReferences
54f2GAGAGATGTTGGGCTCTTGTGTGTTAATTGCATTGAATGTCAGCPasserine and shorebird species[25]
BOOW19GGAAACTTACTAGAAAATAAATGACTGGCTTTCTAACTTTCCCATGCAACPasserine and shorebird species[25]
Calex05TCCAGCTGAAGTCTTCCGTGAATTCCACACCTGTTCGACAGTTCAATAPasserine and shorebird species[25]
HvoB1AAGCAAGGACTTTCCTTCCAGTCTCAAATTGGAACAGAGAAAGGPasserine and shorebird species[25]
Tgu06CGAGTAGCGTATTTGTAGCGAAGGAGCGGTGATTGTTCAGTPasserine and shorebird species[25]
TG04-061GACAATGGCTATGAAATAAATTAGGCAGAAGGGCATTGAAGCACACPasserine and shorebird species[25]
BbuS027TCATGAGGAACTTTCAGTGCTCGAAGAAAGGCAGCTCTCACCBubo bubo[28]
BbuS102AACTGATTTGGAAACCACCATCCTGGAACACCCAGTGTTTGTCBubo bubo[28]
Bbus116GTTTCTGCAGCTGGGTCAGAAACAGTTTCCATGCCTTACGBubo bubo[28]
BbuS132TCATTGTAGGTCCCATCCAACCCATATCTATCAAGCAACCTTGGBubo bubo[28]
Oe050AGAGTTGTCCTTGGTTGGTTCTAGTAACCTCATCTGCOtus elegans[29]
Oe054TCAGAAAGAAAACTTCAGCAACCCATATATGTATACACAGGCACATGCOtus elegans[29]
15A6ACCTCAGAAGCAGACAGAACCCCTTTGCGATTGCTGTAACStrix occidentalis[30]
Age5ACGTTACAGACACCGATTACTTCCAGCCACGCGTCTGATACTTTAccipiter gentilis[31]
Bb101AATAACCCCAATAGAAGCACCAGAAGGAGATGAGACCBubo bubo[21]
Bb126TCTCCAGAAGGGTTGTCATCTGCTAAAACCTTACAGAATAACAGBubo bubo[21]
Oe142TGAATCAGCAAACCTGTGCCTGAGCTAACCTAGAGTCAGCCAGCOtus elegans[29]
Oe171TTTTACAAACTACTAGTGCATGTCTCCAGATGTTGTATTTCAGTGTCAGOtus elegans[29]
Oe3-21GATTAGAGACCCGATTCCACATTATCTGAGTGGAAGGGTAGTGCOtus elegans[29]
Oe3-7GTGGGTTTATTGCCCCCTCGCAGATGAATGAATGGATAGATGGOtus elegans[22]
Oe149CACACATCCATTTTGGGGTGCGGATGCTGGAACTGACCTGCOtus elegans[22]
Oe081GTAAGGGAAGTAGACGTCTGTTGGCAACTTTGTGTCATCTGAAAGGOtus elegans[22]
Oe084GGGCATAGTTAGACCTTTGCAGCACATCTGTTGTTCTCCGTGTTACCOtus elegans[22]
SneD105AGCCTTGGGAGGTTAAGTCCTACACGCAATCACTGAAGCAGTStrix nebulosa[24]
SneD113AGCTCTTTCCCGACAGTGTCTAGCCAAAAAACCCCATATCCTAGStrix nebulosa[24]
SneD202GCTGGGCTTAGAAATTACATGGCTCTGCCCTAATCAGGAACACTStrix nebulosa[24]
SneD211TAGCCCTGTGGATTTGCATTAGTCACCAGAAGTTCAAAGCAGGTAGStrix nebulosa[24]
SneD218GGGTGTGAGAATGACCTACCTGAGCAGACAGGAGATGGCTTTTAStrix nebulosa[24]
Fepo42CGTATACATCGAAATAAATACCCGAATAAAACATCCCTAACCGlaucidium brasilianum[23]
Oe053CTCTGCATCTTAACGCACAGGACCCTCCAAGTGGACAGGAAAAGCOtus elegans[29]
Oe128CGTTGTAAATGATGAATCGCCTAGTGCATGCATGTATACATACAAACCTGGOtus elegans[29]
Oe129GTCACCTCTTGACATCCGAGTAGCGCTAAGAGTCCATTTGCCCATCTGOtus elegans[29]
Table 3. Variability indices: Allele number (NA), effective allele number (NE); observed heterozygosity (HO), expected heterozygosity (HE), unbalanced expected heterozygosity (uHE), percentage of polymorphic loci (P%). Standard errors are shown in brackets.
Table 3. Variability indices: Allele number (NA), effective allele number (NE); observed heterozygosity (HO), expected heterozygosity (HE), unbalanced expected heterozygosity (uHE), percentage of polymorphic loci (P%). Standard errors are shown in brackets.
SpeciesNANEHOHEuHEP%
Strix aluco3.6 (±0.5)2.8 (±0.4)0.528 (±0.091)0.525 (±0.083)0.556 (±0.087)83.33%
Strix nebulosa3.0 (±0.3)2.2 (±0.2)0.565 (±0.096)0.471 (±0.074)0.499 (±0.079)83.33%
Strix uralensis3.0 (±0.5)2.3 (±0.4)0.528 (±0.119)0.429 (±0.0920.468 (±0.100)75.00%
Bubo bubo4.2 (±0.6)2.6 (±0.4)0.539 (±0.090)0.545 (±0.068)0.578 (±0.072)91.67%
Bubo scandiacus3.3 (±0.4)2.6 (±0.3)0.515 (±0.092)0.550 (±0.061)0.584 (±0.065)91.67%
Otus scops6.3 (±0.6)4.6 (±0.5)0.733 (±0.049)0.750 (±0.027)0.798 (±0.030)100.00%
Athene noctua4.1 (±0.5)3.3 (±0.5)0.620 (±0.061)0.604 (±0.060)0.641 (±0.064)100.00%
Asio otus3.2 (±0.4)2.5 (±0.3)0.681 (±0.097)0.529 (±0.062)0.590 (±0.070)91.67%
Surnia ulula2.2 (±0.3)1.70 (±0.21)0.356 (±0.088)0.318 (±0.074)0.338 (±0.079)66.67%
Mean3.7 (±0.2)2.8 (±0.14)0.566 (±0.030)0.525 (±0.025)0.561 (±0.026)87.04 (±3.70)%
Table 4. PID and PID_sib values: In the table, values lower than 0.001 were bold-evidenced. In the PID_sib columns, this threshold wasindicated with (*). The PID_sib never reached this value in Surnia ulula (0.0155), Strix nebulosa (0.00162) and Strix uralensis (0.0246). When PID or PID_sib values stop decreasing it means that the markers are not variable and informative (underlined).
Table 4. PID and PID_sib values: In the table, values lower than 0.001 were bold-evidenced. In the PID_sib columns, this threshold wasindicated with (*). The PID_sib never reached this value in Surnia ulula (0.0155), Strix nebulosa (0.00162) and Strix uralensis (0.0246). When PID or PID_sib values stop decreasing it means that the markers are not variable and informative (underlined).
Strix alucoStrix nebulosaStrix uralensis
locusPIDPID_siblocusPID locusPIDPID_sib
Oe1498.16 × 10−023.78 × 10−0115a61.40 × 10−014.33 × 10−01Oe1288.16 × 10−023.82 × 10−01
15a66.98 × 10−031.46 × 10−01Oe532.38 × 10−021.99 × 10−01Oe539.28 × 10−031.57 × 10−01
SneD2187.10 × 10−045.87 × 10−02SneD2184.82 × 10−039.54 × 10−02Oe1421.25 × 10−036.72 × 10−02
Oe537.22 × 10−052.36 × 10−02Oe1491.01 × 10−034.68 × 10−02Oe1491.86 × 10−042.96 × 10−02
Oe1427.77 × 10−069.59 × 10−03Fepo422.28 × 10−042.35 × 10−02SneD2183.23 × 10−051.38 × 10−02
Oe1281.45 × 10−064.56 × 10−03Oe1425.18 × 10−051.18 × 10−0215a67.33 × 10−066.93 × 10−03
Oe1293.52 × 10−072.42 × 10−03SneD1131.24 × 10−056.07 × 10−03Oe1292.14 × 10−063.82 × 10−03
SneD1139.97 × 10−081.33 × 10−03Oe1293.17 × 10−063.20 × 10−03Fepo421.20 × 10−062.87 × 10−03
Oe3213.77 × 10−087.97 × 10−04 *Oe1281.14 × 10−061.99 × 10−03Tgu068.74 × 10−072.46 × 10−03
Tgu062.11 × 10−085.99 × 10−04Tgu067.53 × 10−071.62 × 10−03Oe0548.74 × 10−072.46 × 10−03
Fepo422.11 × 10−085.99 × 10−04Oe0547.53 × 10−071.62 × 10−03Oe3218.74 × 10−072.46 × 10−03
Oe0542.11 × 10−085.99 × 10−04Oe3217.53 × 10−071.62 × 10−03SneD1138.74 × 10−072.46 × 10−03
Bubo buboBubo scandiacusOtus scops
locusPIDPID_siblocusPIDPID_siblocusPIDPID_sib
SneD2187.65 × 10−023.74 × 10−01Oe538.79 × 10−023.86 × 10−01Oe0542.99 × 10−023.22 × 10−01
Oe537.26 × 10−031.46 × 10−01SneD2188.94 × 10−031.55 × 10−01Oe1491.18 × 10−031.08 × 10−01
15a68.09 × 10−046.01 × 10−02SneD1131.01 × 10−036.48 × 10−02SneD1134.70 × 10−053.60 × 10−02
Oe1429.11 × 10−052.56 × 10−02Oe1421.36 × 10−042.76 × 10−0215a62.73 × 10−061.28 × 10−02
Fepo421.48 × 10−051.16 × 10−02Oe1282.52 × 10−051.31 × 10−02Fepo421.85 × 10−074.70 × 10−03
SneD1132.70 × 10−065.77 × 10−03Oe0545.84 × 10−066.70 × 10−03Oe1421.48 × 10−081.79 × 10−03
Oe1288.03 × 10−073.16 × 10−03Tgu061.50 × 10−063.49 × 10−03SneD2181.47 × 10−097.19 × 10−04*
Oe1492.23 × 10−071.78 × 10−03Oe1494.68 × 10−072.03 × 10−03Oe1281.48 × 10−102.93 × 10−04
Oe0546.17 × 10−081.00 × 10−03Fepo421.68 × 10−071.26 × 10−03Oe532.49 × 10−111.34 × 10−04
Tgu062.33 × 10−085.98 × 10−04 *15a66.83 × 10−087.92 × 10−04 *Oe3214.46 × 10−126.45 × 10−05
Oe3211.55 × 10−084.89 × 10−04Oe1292.78 × 10−084.99 × 10−04Tgu068.29 × 10−133.13 × 10−05
Oe129\\Oe3212.78 × 10−084.99 × 10−04Oe1291.76 × 10−131.54 × 10−05
Surnia ululaAthene noctuaAsio otus
locusPIDPID_siblocusPIDPID_siblocusPIDPID_sib
SneD2181.34 × 10−014.26 × 10−01SneD1133.94 × 10−023.34 × 10−01Tgu061.14 × 10−014.10 × 10−01
Oe1423.56 × 10−022.27 × 10−01Oe1492.42 × 10−031.20 × 10−01SneD2181.36 × 10−021.72 × 10−01
15a61.34 × 10−021.35 × 10−01Oe1291.51 × 10−044.30 × 10−02Fepo421.62 × 10−037.20 × 10−02
Fepo424.63 × 10−038.06 × 10−02SneD2181.46 × 10−051.71 × 10−02Oe1292.41 × 10−043.17 × 10−02
Oe1281.77 × 10−034.86 × 10−02Oe1421.96 × 10−067.27 × 10−03Oe1284.70 × 10−051.50 × 10−02
SneD1135.87 × 10−042.95 × 10−0215a63.66 × 10−073.38 × 10−0315a69.91 × 10−067.33 × 10−03
Tgu062.86 × 10−042.06 × 10−02Fepo426.59 × 10−081.59 × 10−03Oe1493.24 × 10−064.31 × 10−03
Oe1491.60 × 10−041.55 × 10−02Oe3211.34 × 10−087.88 × 10−04 *SneD1131.24 × 10−062.60 × 10−03
Oe0541.60 × 10−041.55 × 10−02Oe0545.22 × 10−094.80 × 10−04Oe3215.04 × 10−071.64 × 10−03
Oe1291.60 × 10−041.55 × 10−02Oe1282.59 × 10−093.39 × 10−04Oe0542.00 × 10−071.06 × 10−03
Oe3211.60 × 10−041.55 × 10−02Tgu061.36 × 10−092.50 × 10−04Oe539.20 × 10−087.18 × 10−04 *
Oe531.60 × 10−041.55 × 10−02Oe537.63 × 10−101.87 × 10−04Oe1429.20 × 10−087.18 × 10−04
Table 5. Parental assignment values from Colony: the lowest values were recovered when the one-parent assignment was tested. In the table, misallocation or low assignment probabilities are highlighted in bold.
Table 5. Parental assignment values from Colony: the lowest values were recovered when the one-parent assignment was tested. In the table, misallocation or low assignment probabilities are highlighted in bold.
Parent Pair ProbabilitySingle Parent Probability
Offspring Father MotherProbability ValueMain ProbabilitySecond Probability
S_al 3S_al 1S_al 21.00
S_al 6S_al 4S_al 51.00
S_al 9S_al 7S_al 81.00
S_ne 3S_ne 1S_ne 21.00S_ne 20.907S_ne 8
S_ne 6S_ne 4S_ne 51.00S_ne 40.992S_ne 1
S_ne 9S_ne 7S_ne 81.00
S_ur 3S_ur 1S_ur 21.00S_ur 20.997S_ur 5
S_ur 6S_ur 4S_ur 51.00S_ur 50.500S_ur 2
B_bu 3B_bu 1B_bu 21.00
B_bu 6B_bu 4B_bu 51.00
B_bu 9B_bu 7B_bu 81.00
B_sc 3B_sc 1B_sc 21,00
B_sc 6B_sc 4B_sc 51.00B_sc 60.999
B_sc 9B_sc 7B_sc 81,00
O_sc 3O_sc 1O_sc 21.00
O_sc 6O_sc 4O_sc 51.00
O_sc 9O_sc 7O_sc 81.00
A_no 3A_no 1A_no 21.00
A_no 6A_no 4A_no 51.00
A_no 9A_no 7A_no 81.00
S_ul 3S_ul 1S_ul 21.00S_ul 10.500S_ul 4
S_ul 6S_ul 4S_ul 51.00S_ul 40.983S_ul 7
S_ul 50.753S_ul 2
S_ul 9S_ul 7S_ul 81.00
A_ot 3A_ot 1A_ot 21.00
A_ot 6A_ot 7A_ot 81.00
Table 6. Parentage analysis in Cervus using the LOD and Delta computation. The different confidence values are indicated by the following codes: * = confidence level higher than 95%. In the table, misallocation or low assignment probabilities are highlighted in bold.
Table 6. Parentage analysis in Cervus using the LOD and Delta computation. The different confidence values are indicated by the following codes: * = confidence level higher than 95%. In the table, misallocation or low assignment probabilities are highlighted in bold.
LODDelta
OffspringCandidate Mother Pair ConfidenceCandidate FatherPair ConfidenceTrio ConfidenceOffspringCandidate MotherPair ConfidenceCandidate FatherPair Confidence
S_al 3S_al 1*S_al 2**S_al 3S_al 1*S_al 2*
S_al 6S_al 3*S_al 4**S_al 6S_al 3*S_al 4*
S_al 9S_al 7*S_al 8**S_al 9S_al 7*S_al 8*
S_ne 3S_ne 1*S_ne 2**S_ne 3S_ne 1*S_ne 2*
S_ne 6S_ne 4*S_ne 5**S_ne 6S_ne 4*S_ne 5*
S_ne 9S_ne 7*S_ne 8**S_ne 9S_ne 7*S_ne 8*
S_ur 3S_ur 1*S_ur 2**S_ur 3S_ur 1*S_ur 2*
S_ur 6S_ur 4 S_ur 5**S_ur 6S_ur 4 S_ur 5*
B_bu 3B_bu 1*B_bu 2**B_bu 3B_bu 1*B_bu 2*
B_bu 6B_bu 4*B_bu 5**B_bu 6B_bu 4*B_bu 5*
B_bu 9B_bu 7*B_bu 8**B_bu 9B_bu 7*B_bu 8*
B_sc 3B_sc 1*B_sc 2**B_sc 3B_sc 1*B_sc 2*
B_sc 6B_sc 4*B_sc 5**B_sc 6B_sc 4*B_sc 5*
B_sc 9B_sc 7*B_sc 8**B_sc 9B_sc 7*B_sc 8*
O_sc 3O_sc 1*O_sc 2**O_sc 3O_sc 1*O_sc 2*
O_sc 6O_sc 4*O_sc 5**O_sc 6O_sc 4*O_sc 5*
O_sc 9O_sc 7*O_sc 8**O_sc 9O_sc 7*O_sc 8*
A_no 3A_no 1*A_no 2**A_no 3A_no 1*A_no 2*
A_no 6A_no 4*A_no 5**A_no 6A_no 4*A_no 5*
A_no 9A_no 7*A_no 8**A_no 9A_no 7*A_no 8*
S_ul 3S_ul 1*S_ul 5**S_ul 3S_ul 1*S_ul 5*
S_ul 6S_ul 4*S_ul 5**S_ul 6S_ul 4*S_ul 5*
S_ul 9S_ul 7*S_ul 8**S_ul 9S_ul 7*S_ul 8*
A_ot 3A_ot 1*A_ot 2**A_ot 3A_ot 1*A_ot 2*
A_ot 6A_ot 7*A_ot 8**A_ot 6A_ot 7*A_ot 8*
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giangregorio, P.; Naldi, L.; Mengoni, C.; Greco, C.; Padula, A.; Zaccaroni, M.; Fani, R.; Argenti, G.; Mucci, N. Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes. Genes 2021, 12, 1721. https://doi.org/10.3390/genes12111721

AMA Style

Giangregorio P, Naldi L, Mengoni C, Greco C, Padula A, Zaccaroni M, Fani R, Argenti G, Mucci N. Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes. Genes. 2021; 12(11):1721. https://doi.org/10.3390/genes12111721

Chicago/Turabian Style

Giangregorio, Patrizia, Lorenzo Naldi, Chiara Mengoni, Claudia Greco, Anna Padula, Marco Zaccaroni, Renato Fani, Giovanni Argenti, and Nadia Mucci. 2021. "Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes" Genes 12, no. 11: 1721. https://doi.org/10.3390/genes12111721

APA Style

Giangregorio, P., Naldi, L., Mengoni, C., Greco, C., Padula, A., Zaccaroni, M., Fani, R., Argenti, G., & Mucci, N. (2021). Cross-Amplification in Strigiformes: A New STR Panel for Forensic Purposes. Genes, 12(11), 1721. https://doi.org/10.3390/genes12111721

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop