Ministry of Education Key Laboratory for Biodiversity and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
2.
College of Life Sciences, Leshan Normal University, Leshan, Sichuan 614000, China
The Emei Shan Liocichla (Liocichla omeiensis) is an endemic bird species to southwestern China with a small geographic range. However, little was known about the genetic status of this threatened species.
Methods
We applied restriction-site-associated DNA sequencing (RAD-Seq) for rapid mass identification of microsatellite markers of the Emei Shan Liocichla.
Results
A total of 11, 564 microsatellite sequences were obtained, 600 random loci were designed for screening and 24 polymorphic microsatellite loci were selected for further validation. The average allele number, average observed heterozygosity and average expected heterozygosity were relatively low in our samples, which were 6.08, 0.6618 and 0.7048, respectively, indicating that the Emei Shan Liocichla might have lost some genetic diversity. Further analyses suggested that the populations distributed on two mountains (Daxiangling and Xiaoliangshan) showed a modest degree of genetic differentiation.
Conclusions
These novel microsatellite markers provided valuable preliminary knowledge regarding the genetic status of the Emei Shan Liocichla and can be useful in further studies, as well as in the management and conservation of this species.
Over the past 40 years there have been intensive theoretical and experimental investigations of how natural selection moulds the sex ratio. Fisher's (1930) classic frequency-dependent model of sex allocation predicts that the offspring sex ratio should not deviate from parity as long as the costs of producing males and females are equal. Modern empirical sex ratio research began when Hamilton (1967) observed that many insects and mites have highly female-biased sex ratios and that this trait is associated with high levels of brother-sister mating. Hamilton explained this result by his theory of local mate competition, which has been extended to other cases of interactions between siblings. Briefly, if one sex suffers a greater reduction in fitness through competition with siblings of the same sex, the mother is programmed in her selection of the sex ratio to favor the sex that suffers less competition. Fisher's frequency-dependent argument is thus modified in populations where relatives compete. A second major development was by Trivers and Willard (1973) who originally proposed that maternal quality should affect the direction of this investment for species, exhibiting high variation in male reproductive success. Modern classic sex allocation theory predicts that selection should act on parents to sex-specifically vary the level of investment in offspring when the fitness returns differ for the two sexes (Charnov, 1982; Frank, 1990; Hardy, 2002).
Birds have long been considered incapable of facultatively shifting their primary sex ratio because their sex determination is chromosomal (Clutton-Brock, 1986). However, recent studies have suggested that sex ratio adjustment in birds is not as constrained as thought earlier (West and Sheldon, 2002). In birds, the sex-ratio adjustment can occur at laying period (primary sex-ratio adjustment) and during the period of provisioning for the young (secondary sex-ratio adjustment) (Ellegren and Sheldon, 1997; Kilner, 1998).
Mechanisms for both primary and secondary sex-ratio adjustment in birds remain largely unknown (Emlen, 1997; Sheldon, 1998; Krackow, 1999; Martins, 2004), but preovulatory mechanisms have been suggested as the most efficient means of primary sex-ratio adjustment for some species (Komdeur et al., 2002). Pike and Petrie (2003) presented an excellent review of the potential mechanisms of avian sex manipulation, such as steroid hormones, paternal attractiveness, maternal condition, photoperiod and phenotypic sex and others. In addition, much new experimental and relevant evidence have appeared since then, providing new and interesting perspectives (Pike and Petrie, 2006; Rutkowska and Badyaev, 2008).
Sex allocation patterns among raptor populations have attracted considerable research attention because most raptorial birds exhibit reversed sexual size dimorphism (females are the larger sex) (Newton and Marquiss, 1979; Olsen and Cockburn, 1991; Appleby et al., 1997; Brommer et al., 2003; McDonald et al., 2005). One consequence of this size difference is a disparity in food requirements: daughters require more food than sons do during the nestling period (Anderson et al., 1993a; Krijgsveld et al., 1998; Riedstra et al., 1998; Laaksonen et al., 2004). Rearing larger daughters is more costly since they require more parental investment in the form of food requirements as opposed to sons. Theoretically, the offspring sex ratio in raptors should be male-biased. However, there seems no consistent pattern in the offspring sex ratios of raptors; the proportion of males is either increasing (Newton and Marquiss, 1979; Olsent and Cockburn, 1991; Korpimäki et al., 2000) or decreasing (Tella et al., 1996) and no clear adaptive explanations for these trends have been found so far.
The Eurasian Kestrel (Falco tinnunculus), a sexually size-dimorphic raptor (females ca. 20% heavier than males; Village, 1990), widespread in open country throughout the Palearctic, Afrotropical and Oriental regions (Cramp and Simmons, 1980). A total of 12 subspecies have been identified (Dickinson, 2003). The subspecies in our study is F. tinnunculus interstinctus. The female Eurasian Kestrel might be able to adjust her reproductive effort by promoting sex differences in the mass of first-laid eggs (Blanco et al., 2003), where sexual differences in the size of first-laid eggs defines the pattern of within-brood mass hierarchies of fledglings. In addition, some researchers have documented that females have the potential ability to adjust resources allocated to different-sex eggs laid in varying order (Anderson et al., 1997; Blanco et al., 2003; Martínez-Padilla and Fargallo, 2007).
We investigated the primary and secondary offspring sex ratio of the Eurasian Kestrel at the Zuojia Nature Reserve in northeastern China. Our objectives were: 1) to test whether its offspring sex ratio demonstrates a general deviation from a balanced sex ratio and 2) whether a pattern of biased sex ratio exists in relation to laying order, egg mass, hatching order, laying date and clutch size.
Study area and methods
Study area
We conducted the study in the Zuojia Nature Reserve in Jilin Province, northeastern China during the breeding seasons (from late March to late July) in 2007 and 2008. The Zuojia Nature Reserve stretches from the eastern Changbai mountains to the western plains (126°1′–127°2′N, 44°6′–45°5′E), where the elevation ranges from 200 m to 530 m. The region is subject to an eastern monsoon climate, characterized by hot, dry summers and cold, snowy winters; the area is covered by approximately 35% open habitat and 65% forests. Vegetation in the area is diverse and the existing forests are secondary.
Data collection
The Eurasian Kestrels do not build their own nests, but rely mainly on the stick nests of the Magpie (Pica pica) in our study area (Xiang et al., 1991; Deng et al., 2006; Zhou et al., 2009). Most magpie nests used by the kestrels were inaccessible. We therefore installed nest boxes to develop the study and most kestrels then shifted to breed in these boxes.
The nest boxes we used in the study were constructed from rough-cut boards and put together with exterior nails or deck screws. The internal dimensions of the boxes were all the same: 50 cm deep, a 35 cm × 35 cm bottom area and with a 12.5 cm or 15 cm diameter entrance hole near the top. These nest boxes were placed randomly about 8–13 m above ground in various tree species with distances ranging from 50 m to 150 m between the nearest two boxes. The number of boxes monitored was 68 in 2007 and 77 in 2008.
We monitored nest boxes used by the Eurasian Kestrels every day to record the laying date and marked the laying sequence of each egg accurately during the egg laying period. We used an electronic balance to weigh eggs to the nearest 0.01 g on the day of being laid. In order to assign all marked eggs to their corresponding nestlings, and record the hatching order, we checked nests hourly (between 08:00 and 19:00 hours) on the days of expected hatching. We distinguished nestlings by drawing some simple symbols on their head with indelible ink. We removed unhatched eggs from the nest after several days when the rest of the clutch had hatched. We discarded eggs containing no visible embryos, and stored dead embryos in 100% ethanol. After hatching, we visited nests at an interval of two days to collect body tissues of the dead nestlings as a source of DNA and to monitor fledging success. We collected a blood sample of each nestling (20–50 μL) from brachial veins with a sterilized hypodermic needle, 9 to 11 days after hatching. We stored blood samples in EDTA for determination of sex in the laboratory.
Sex determination
We sexed the nestlings using a polymerase chain reaction (PCR) amplification based on the technique used by Fridolfsson and Ellegren (1999). We isolated DNA from the red blood cells using proteinase K digestion and phenol-chloroform extraction. We ran PCR amplification using a particular set of primers (2550F and 2718R), as suggested by Fridolfsson and Ellegren (1999). The temperature profile in PCR was as follows: hotstart at 94℃, initial denaturation at 94℃ for 5 min, followed by a "touch-down" scheme (Don et al., 1991) where the annealing temperature was lowered 1℃ per cycle, starting from 58℃, until a temperature of 48℃ was reached. Then 25–35 additional cycles were run at a constant annealing temperature of 48℃. Denaturation was at 94℃ for 30 s, annealing for 45 s and extension at 72℃ for 1 min. A final extension step of 5 min was added after the last cycle. PCR products, separated by electrophoresis on 2% agarose gels containing ethidium bromide, were visualized under UV light. We assigned the sex of individuals from five adult pairs correctly using this technique, so we could be confident the method worked successfully in this species.
Statistics analysis
We used two data sets to analyze the offspring sex ratio. The data set for analyzing the primary sex ratio included broods that all laid eggs had hatched successfully, including dead embryos. We sampled all nestlings (both blood and tissue samples) for a total of 81 nestlings in 15 broods. The data set for analyzing the secondary sex ratio included 105 nestlings in 23 successfully fledged broods. We examined the brood sex ratio (hereafter referred to as the proportion of males) using the binomial test (Zar, 1984). We used the Wilcoxon signed rank test in S-PLUS 6.0 to test the primary and secondary sex ratios at population level. In order to study the relationship between primary sex ratio, laying order, clutch size and egg mass, we performed linear mixed model (LMM) in SPSS 14.0 (SPSS Inc., Chicago, Illinois, USA). In this model, we entered the sex ratio as the response variable, clutch size and laying order as fixed factors. The kestrel hatches asynchronously with a hatching spread of 3–5 days and therefore we divided the laying order into three categories: initial, middle and last. The initial category corresponded to the first eggs and the last category to the last eggs. The other eggs were defined as the middle category (Hardy, 2002). Thus, the sex ratio used in the model stood for the sex ratio of each category. Moreover, egg mass in the model also represented the egg mass of each category. We encoded the laying dates into five stages: March, early April (April Ⅰ), middle April (April Ⅱ), late April (April Ⅲ) and May and introduced it as a random factor, with egg mass as covariate. By introducing the laying date as a higher level random effect in addition to the residual error term, our LMM permitted an estimation of the fixed effects and interaction parameters while accounting for random variation in the sex ratio between broods. We analyzed the relationship between secondary sex ratio, laying order, clutch size, hatching order and egg mass using the same model. For all statistical tests, we used a significance level of 0.05. We presented values as means ± SD.
Results
Sex ratio at brood and population level
The brood sex ratio of the Eurasian Kestrel in eggs was 47.0±4.3%, not statistically biased from 0.5 (Wilcoxon signed rank test, Z15 = –0.402, p = 0.342). However, the brood sex ratio in fledglings was 40.8±4.5%, significantly different from 0.5 (Wilcoxon signed rank test, Z23 = –1.882, p = 0.029).
The primary sex ratio at population level was 48.2%, which did not depart from parity (binomial test, n = 81, p = 0.824). The sex ratio in eggs for the two study years showed a similar pattern (binomial test, 48.9%, n = 45, p = 1.000 in 2007; 47.2%, n = 36, p = 0.868 in 2008). The secondary sex ratio at population level was 42.9%, which did not depart significantly from parity (binomial test, n = 105, p = 0.172). Moreover, we detected the same tendency in nestlings (binomial test, 41.9%, n = 62, p = 0.253 in 2007; 44.2%, n = 43, p = 0.542 in 2008).
Effect of laying order, egg mass, clutch size on primary brood sex ratio
The sex ratio of eggs for different categories of laying order did not differ from parity (binomial tests: initial, 40.0%, n = 15, p = 0.607; middle, 43.1%, n = 51, p = 0.401; last, 73.3%, n = 15, p = 0.118). Results of LMM showed that the sex ratio in eggs was not associated with laying order, clutch size and their interaction (Table 1). Our results indicate that the laying date had a significant impact on primary offspring sex (Table 1), but there was no significant correlation (r = –0.025, p = 0.929) (Fig. 1a). Egg mass did not have a significant effect on primary sex ratio, but the egg mass of males (23.15±1.36, n = 42) was slightly heavier than that of females (22.89±1.28, n = 38), although the difference was not statistically significant (Mann-Whitney test U = 712, p = 0.407). In addition, egg mass was not associated linearly with laying order; the eggs in the middle position in the laying order were the heaviest (Fig. 2).
Table
1.
LMM analyses of the effects of clutch size and laying order on the sex of the eggs in Eurasian Kestrels
Source
NDF
DDF
FNDF, DDF
p
Clutch size
3
8
1.570
0.271
Laying order
1
8
0.004
0.951
Egg mass
1
8
1.796
0.217
Clutch size × Laying order
1
8
0.000
0.983
Covariance parameter estimates:
Estimate
SE
Z
p
Laying date
0.262
0.130
2.000
0.046
The model accounts for random variation between laying date. NDF represents the degrees of freedom of the numerator and DDF that of the denominator. FNDF, DDF stands for Fisher's test statistics, p represents the significance level, "estimate" represents the estimated covariance parameters, SE stands for standard error, Z represents Z statistic value of the Wald test.
Effects of laying order, egg mass, clutch size and hatching order on secondary brood sex ratio
The secondary brood sex ratio was significantly related to clutch size (F = 5.297, p = 0.017) and egg mass (F = 8.099, p = 0.016). With increasing clutch size and egg mass, the proportion of males became gradually higher, but there were no significantly positive correlations (Figs. 3 and 4). Laying order and hatching order did not have a significant effect on the proportion of males in fledglings; however, when they interacted with clutch size, both factors had a significant effect on the secondary brood sex ratio (Table 2). In addition, the LMM showed that laying date had a significant effect on secondary offspring sex (Table 1), but there was no significant correlation (r = –0.267, p = 0.218) (Fig. 1b).
Figure
3.
Relationship between secondary brood sex ratio and clutch size in Eurasian Kestrel in northeastern China
Table
2.
LMM analyses of the effects of clutch size, laying order and hatching order on the sex of the fledgings in Eurasian Kestrel
Source
NDF
DDF
FNDF, DDF
p
Clutch size
3
11
5.297
0.017
Laying order
1
11
0.851
0.376
Hatching order
1
11
0.128
0.728
Egg mass
1
11
8.099
0.016
Clutch size × Laying order
1
11
6.447
0.028
Clutch size × Hatching order
1
11
5.508
0.039
Covariance parameter estimates:
Estimate
SE
Z
p
Laying date
0.194
0.008
2.345
0.019
The model accounts for random variation between laying date. NDF is the degree of freedom of the numerator and DDF of the denominator, FNDF, DDF is Fisher's test statistic, p represents significance level, "estimate" represents the estimated covariance parameters, SE stands for standard error, and Z represents the Z statistic value of the Wald test.
In this study, the sex ratio at population level and the primary brood sex ratio did not depart from parity, which was similar to the result obtained by Ankney (1982). We found that the secondary brood sex ratio of Eurasian Kestrels was consistently female-biased across the two-year study. A number of studies have shown that females can adjust sex ratio of their offspring by adjusting the laying and hatching orders (Ankney, 1982; Kilner, 1998; Cichoń et al., 2003; Ležalová et al., 2005). However, our results, as another study (Martínez-Padilla and Fargallo, 2007) on the species, indicated that the laying order did not affect the sex ratio of Eurasian Kestrels. Interestingly, our results showed that clutch size and egg mass only affected secondary brood sex ratio, while the laying date affected both primary and secondary sex ratios.
For many bird species, there is no a consistent relationship between egg mass and sex ratio, because females might adjust egg mass either positively or negatively to correlate with the laying order in assigning resources to different sexual offspring (Fiala, 1981; Ankney, 1982; Bancroft, 1984). Anderson et al. (1997) suggested that the sex ratio of offspring was significantly associated with the laying order and that male eggs were significant larger than female eggs in the American Kestrel (F. sparverius). In our study, egg mass had a significant effect on the secondary brood sex ratio and showed the same pattern as the American Kestrel in spite of a lack of significant difference between female and male eggs. Therefore, we speculate that the sexually size-dimorphic Eurasian Kestrel was able to adjust the sex ratio of its offspring by adopting different strategies of egg mass provisioning according to egg sex and laying order.
Some studies showed that the fledgling sex ratios of several raptor species are skewed toward males early in the season while late broods had an excess of females (Dijkstra et al., 1990; Daan et al., 1996; Tella et al., 1996; Griggio et al., 2002). In our study, the laying date had a significant effect on the primary and secondary brood sex ratios. Moreover, both sex ratios showed a slightly decreasing trend as the season progressed in spite of a lack of significant shift. Eurasian Kestrel nestlings are sexually dimorphic, with daughters larger than sons. The larger daughters may need more foods and have an advantage during sibling competition for food (Anderson et al., 1993a, 1993b). Sex-specific characteristics of nestlings may affect the reproductive biology of the species, while nestling mortality should fall more heavily on sons than on daughters. We suggest that the seasonal variation in the sex ratio of the Eurasian Kestrel might be associated with favouritism for the male function on the part of the parents, both to counter the risk of excess mortality of sons and because sons require less parental care per capita during earlier breeding seasons.
A large number of studies show that clutch size does not affect the sex ratio of birds, except for the study by Wegge (1980), who indicated that in a declining population of Capercaillie (Tetrao urogallus), the sex ratio of fledglings declined with clutch size, but that this was probably a consequence of differential mortality between hatchings and fledgings. In our study, clutch size had a significant effect on the secondary sex ratio in spite of a lack of significant correlation. The highest male proportion occurred with a clutch size of seven, with the highest probability of mortality of hatchings from the seventh egg. Two potential reasons may explain this conclusion. First, we found that the Eurasian Kestrel, bred in an artificial nest box, had a larger clutch size than kestrels bred in natural nests (5.77±0.15 vs 4.86± 0.26), where clutch size of the species was not more than 6 in natural nests (W. Gao unpublished data). As a result, it was difficult for the parents to rear all the nestlings in this condition. We therefore considered that the kestrel may adopt a "brood reduction strategy" in order to adjust the secondary brood sex ratio (Slagsvold et al., 1984). Secondly, food availability might also constrain the survival rate of nestlings. Kestrels feed mainly on small mammals, especially voles, but birds, frogs and insects also form part of their diet during the nestling period, which might reflect a reduced availability of rodents due to their being preyed upon continuously by the kestrels and several other raptor species in the study area (Geng et al., 2009).
Our results revealed relatively complex patterns of sex allocation in the Eurasian Kestrel that appears to be associated with both maternal ability to modify the sex of the progeny by adapting to its environment. Future work on the kestrel should include experiments that evaluate the ultimate mechanism of sex ratio allocation and parental provisioning patterns to different sex eggs with egg mass, adjusting for reproductive efforts.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant No. 30400047), the Training Fund of Northeast Normal University Scientific Innovation Project (Grant No. 07013) and the Jilin Provincial Science and Technology Department (Grant No. 20070558 and 20090575). The authors thank Shi Jiao, Bo Qin and Yufang Qin for assistance in the field and anonymous reviewers for providing comments to improve the manuscript.
Andrews KR, Good JM, Miller MR, Luikart G, Hohenlohe PA. Harnessing the power of RADseq for ecological and evolutionary genomics. Nat Rev Genet. 2016;17:81-92.
Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, Selker EU, Cresko WA, Johnson EA. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE. 2008;3:e3376.
Barchi L, Lanteri S, Portis E, Acquadro A, Valè G, Toppino L, Rotino GL. Identification of SNP and SSR markers in eggplant using RAD tag sequencing. BMC Genom. 2011;12:1-9.
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 Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier (France); 2004.
BirdLife International. Liocichla omeiensis. The IUCN Red List of Threatened Species 2017. Cambridge: BirdLife International; 2017.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114-20.
Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Human Genet. 1980;32:314-31.
Bruford MW, Wayne RK. Microsatellites and their application to population genetic studies. Curr Opin Genet Dev. 1993;3:939-43.
Carvajal-Rodriguez A, de Uña-Alvarez J. Assessing significance in high-throughput experiments by sequential goodness of fit and q-value estimation. PLoS ONE. 2011;6:e24700.
Castoe TA, Poole AW, de Koning APJ, Jones KL, Tomback DF, Oyler-McCance SJ, Fike JA, Lance SL, Streicher JW, Smith EN, Pollock DD. Rapid microsatellite identification from illumina paired-end genomic sequencing in two birds and a snake. PLoS ONE. 2012;7:e30953.
Chapuis M-P, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24:621-31.
Chapuis MP, Plantamp C, Streiff R, Blondin L, Piou C. Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations. Mol Ecol. 2015;24:6107-19.
Collar N, Robson C, de Juana E. Grey-cheeked Liocichla (Liocichla omeiensis). In: del Hoyo J, Elliott A, Sargatal J, Christie DA, de Juana E, editors. Handbook of the birds of the world alive. Barcelona: Lynx Edicions; 2016.
Dawson DA, Horsburgh GJ, Küpper C, Stewart IRK, Ball AD, Durrant KL, Hansson B, Bacon IDA, Bird S, Klein Á, Krupa AP, Lee JW, Martín-Gálvez D, Simeoni M, Smith G, Spurgin LG, Burke T. 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-94.
Duncan CJ, Worth JRP, Jordan GJ, Jones RC, Vaillancourt RE. Genetic differentiation in spite of high gene flow in the dominant rainforest tree of southeastern Australia, Nothofagus cunninghamii. Heredity. 2016;116:99-106.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564-7.
Faria J, Pita A, Rivas M, Martins GM, Hawkins SJ, Ribeiro P, Neto AI, Presa P. A multiplex microsatellite tool for conservation genetics of the endemic limpet Patella candei in the Macaronesian archipelagos. Aquat Conserv. 2016;26:775-81.
Frankham R, Briscoe DA, Ballou JD. Introduction to conservation genetics. Cambridge: Cambridge University Press; 2002.
Fu Y, Zhang Z. The influence of abnormal low-temperature of spring and summer on the breeding of Omei Shan Liocichla (Liocichla omeiensis). Beijing Norm Univ (Nat Sci). 2011;47:292-5 (in Chinese).
Fu Y, Dowell SD, Zhang Z. Breeding ecology of the Emei Shan Liocichla (Liocichla omeiensis). Wilson J Ornithol. 2011;123:748-54.
Gill FB. Ornithology. New York: W.H. Freeman & Company; 2007.
Gu LY, Liu Y, Wang N, Zhang ZW. A panel of polymorphic microsatellites in the Blue Eared Pheasant (Crossoptilon auritum) developed by cross-species amplification. Chin Birds. 2012;3:103-7.
Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, Lepoittevin C, Malausa T, Revardel E, Salin F, Petit RJ. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11:591-611.
Hua F, Wang X, Zheng X, Fisher B, Wang L, Zhu J, Tang Y, Douglas WY, Wilcove DS. Opportunities for biodiversity gains under the world's largest reforestation programme. Nat Commun. 2016;7:12717.
Huang J, Zhu D, Song X, Chen B, Zeng C, Moermond T, Zhang X, Yue B. High-throughput microsatellite markers discovery for the Sichuan Hill Partridge (Arborophila rufipectus) and assessment of genetic diversity in the Laojunshan population. Biochem Syst Ecol. 2015;60:266-72.
Jamieson IG. Founder effects, inbreeding, and loss of genetic diversity in four avian reintroduction programs. Conserv Biol. 2011;25:115-23.
Jarne P, Lagoda PJL. Microsatellites, from molecules to populations and back. Trends Ecol Evol. 1996;11:424-9.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870-4.
Lalitha S. Primer premier 5. Biotech Softw Internet Rep. 2000;1:270-2.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754-60.
Li Y, Korol AB, Fahima T, Beiles A, Nevo E. Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review. Mol Ecol. 2002;11:2453-65.
Marshall TC, Slate J, Kruuk LEB, Pemberton JM. Statistical confidence for likelihood-based paternity inference in natural populations. Mol Ecol. 1998;7:639-55.
McKay BD, Mays HL, Peng Y-W, Kozak KH, Yao C-T, Yuan H-W. Recent range-wide demographic expansion in a Taiwan endemic montane bird, Steere's Liocichla (Liocichla steerii). BMC Evol Biol. 2010;10:71.
Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 2007;17:240-8.
Nagy S, Poczai P, Cernák I, Gorji AM, Hegedűs G, Taller J. PICcalc: an online program to calculate polymorphic information content for molecular genetic studies. Biochem Genet. 2012;50:670-2.
Pfender WF, Saha MC, Johnson EA, Slabaugh MB. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theor Appl Genet. 2011;122:1467-80.
Ribon R, Simon JE, Theodoro De Mattos G. Bird extinctions in Atlantic forest fragments of the Viçosa region, southeastern Brazil. Conserv Biol. 2003;17:1827-39.
Ryman N, Leimar O. Effect of mutation on genetic differentiation among nonequilibrium populations. Evolution. 2008;62:2250-9.
Sudheer PDVN, Mastan SG, Rahman H, Prakash CR, Singh S, Reddy MP. Cross species amplification ability of novel microsatellites isolated from Jatropha curcas and genetic relationship with sister taxa. Mol Biol Rep. 2010;38:1383-8.
Tregenza T, Wedell N. Polyandrous females avoid costs of inbreeding. Nature. 2002;415:71-3.
Twedt DJ, Wilson RR, Henne-Kerr JL, Grosshuesch DA. Avian response to bottomland hardwood reforestation: the first 10 years. Restor Ecol. 2002;10:645-55.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535-8.
Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163:1177-91.
Zane L, Bargelloni L, Patarnello T. Strategies for microsatellite isolation: a review. Mol Ecol. 2002;11:1-16.
Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821-9.
Zhang L, Zhang ZH, Shen FJ, Hou R, Zhang WP, Liu YL, Tu KY, Yang AL. Identification and characterization of polymorphic microsatellite loci in the red-crowned crane. Genet Mol Res. 2015;14:15169-76.
John A. Nash, Richard C. Harrington, Kristof Zyskowski, et al. Species status and phylogenetic relationships of the enigmatic Negros Fruit Dove (Ptilinopus arcanus). Ibis, 2024.
DOI:10.1111/ibi.13305
2.
Andrew Ross T. Reintar, Lisa J. Paguntalan, Philip Godfrey C. Jakosalem, et al. Habitat preference and population density of threatened Visayan hornbills Penelopides panini and Rhabdotorrhinus waldeni in the Philippines. Journal of Threatened Taxa, 2022, 14(3): 20713.
DOI:10.11609/jott.7261.14.3.20713-20720
Table
1.
LMM analyses of the effects of clutch size and laying order on the sex of the eggs in Eurasian Kestrels
Source
NDF
DDF
FNDF, DDF
p
Clutch size
3
8
1.570
0.271
Laying order
1
8
0.004
0.951
Egg mass
1
8
1.796
0.217
Clutch size × Laying order
1
8
0.000
0.983
Covariance parameter estimates:
Estimate
SE
Z
p
Laying date
0.262
0.130
2.000
0.046
The model accounts for random variation between laying date. NDF represents the degrees of freedom of the numerator and DDF that of the denominator. FNDF, DDF stands for Fisher's test statistics, p represents the significance level, "estimate" represents the estimated covariance parameters, SE stands for standard error, Z represents Z statistic value of the Wald test.
Table
2.
LMM analyses of the effects of clutch size, laying order and hatching order on the sex of the fledgings in Eurasian Kestrel
Source
NDF
DDF
FNDF, DDF
p
Clutch size
3
11
5.297
0.017
Laying order
1
11
0.851
0.376
Hatching order
1
11
0.128
0.728
Egg mass
1
11
8.099
0.016
Clutch size × Laying order
1
11
6.447
0.028
Clutch size × Hatching order
1
11
5.508
0.039
Covariance parameter estimates:
Estimate
SE
Z
p
Laying date
0.194
0.008
2.345
0.019
The model accounts for random variation between laying date. NDF is the degree of freedom of the numerator and DDF of the denominator, FNDF, DDF is Fisher's test statistic, p represents significance level, "estimate" represents the estimated covariance parameters, SE stands for standard error, and Z represents the Z statistic value of the Wald test.