A preliminary study of cross-amplified microsatellite loci using molted feathers from a near-threatened Painted Stork (Mycteria leucocephala) population of north India as a DNA source

In continuation of an earlier study in which we reported the cross-amplification of Wood stork microsatellites on the DNA obtained from molted feathers of Painted stork (Mycteria leucocephala), here we investigated the nature of cross-amplified microsatellites and the effect of non-invasive samples on cross-amplification success. In a limited manner, we also addressed the genetic diversity and differentiation in a north Indian population of the Painted Stork examined over three nesting seasons. Among the nine cross-amplified loci, only 5 were polymorphic. Three and 6 loci exhibited low (< 50%) and high amplification success rates (> 80), respectively. For 36 of 145 samples most of the loci failed to amplify. For genetic diversity, only 3 loci could be used since others exhibited low amplification and linkage disequilibrium. Probability of identity (0.034) was not low enough to develop a confidence that the similar genotypes originate from the same individual. Forty-two unique genotypes were identified. In 3 loci, a low to moderate level of genetic diversity (mean He = 0.435) was reported. Non-significant Fst (0.003, P = 0.230), G’stH (0.005, P = 0.247) and Dest (0.003, P = 0.250) values indicate a lack of structuring in temporally distributed populations of Delhi Zoo. The limitations and uniqueness of this study are discussed.


Introduction
The near threatened Painted Stork (PS) (Mycteria leucocephala) is widespread across south and south-east Asia with a stronghold in India. While in the non-breeding season these waterbirds remain widely dispersed across the countryside, after the cessation of the summer monsoon showers they start congregating at their traditional nesting sites which are trees planted on islands in village tanks or in urban premises. One important nesting site in north India is located within the confines of the Delhi Zoological Park in India's capital city. A wild PS population established itself in 1959-1960 and has been regularly visiting the park for nesting since its introduction [1]. No previous genetic studies of the PS have been conducted in India, yet such studies are key to understating geographic patterns of genetic diversity of these storks in the Indian subcontinent and devising conservation strategies [2].
In continuation of our earlier work [3] in which we reported the successful cross-amplification of Wood Stork (WS) (Mycteria americana) microsatellites on PS DNA obtained exclusively from molted feathers, here we investigate the nature of cross-amplified loci and the effect of non-invasive samples on cross-amplification Open Access BMC Research Notes *Correspondence: ajurfi@rediffmail.com 1 Department of Environmental Studies, University of Delhi, Delhi 110007, India Full list of author information is available at the end of the article success of different loci. In addition, we report results of estimates of genetic diversity obtained across three nesting seasons.

Materials and methods
The present study was conducted by using DNA extracted from molted feathers obtained from the PS nesting colony in the premises of the National Zoological Park (Delhi Zoo), Delhi (28° 36′ 21.11″ N, 77° 14′ 48.47″ E). One hundred and forty-five feather samples (142 molted and 2 plucked from dead individuals) were collected over a 3 year period: (1)  For DNA extraction, barbs of the vane were removed and the calamus region was cut and washed in 70% ethanol for 30 min followed by a 30 min wash in ddH 2 O. The basal tip of the calamus and the blood clot from the superior umbilicus [4] were put in a microcentrifuge tube and finely chopped with sterile scissors. The details of the standardized isopropanol method of DNA extraction [5] have already been described in Sharma et al. [3].

Data analysis
PCR products were sized with GeneMapper Software v4.0 (Applied Biosystems). GenAlEx v 6.5 [7,8] was used for the estimation of allele frequencies, performing Chi square test for deviations from Hardy-Weinberg equilibrium (HWE), and estimating probability of identity (PI), diversity measures (observed and expected heterozygosities, and observed and effective number of alleles), and statistically significant AMOVA Fst, G'stH (Hedrick's standardized Gst) and Dest (Jost's estimate of differentiation) [9] for genetic structuring between temporally separated populations. Linkage disequilibrium (LD) between pair of loci was tested by GenePop web version 4.2 [10,11] by 10,000 dememorisation 500 batches and 10,000 iterations per batch. Polymorphic information content (PIC) was calculated by Cervus version 3.0.7 [12]. An assignment test was performed manually to assign the unique genotypes to different populations.
Only those polymorphic loci that showed high amplification success i.e. WSµ09, WSµ17, WSµ18 and WSµ23 were subjected to population genetic analysis. Since a highly significant LD was detected between WSµ17 and WSµ18 (P = 0.0000, across all populations; according to Fisher's method) and WSµ17 was not in HWE in any of the nesting seasons and showed lower amplification success and PIC-value than WSµ18 (Table 2), WSµ17 was dropped from further analysis. Diversity measures from different loci are given in Table 2 The PI of multilocus genotypes for different breeding seasons is given in Table 2. The multilocus PI value of data (0.034) is much higher than the logical cut-off point value (0.0092) (the reciprocal of sample size, 109) i.e. the level of PI below which we have some confidence that two similar genotypes are likely to originate from the same individual. However, 42 unique genotypes were identified, 27 of which were represented by single samples, while the remaining 15 were represented by multiple samples.

Discussion
The nature and size of cross-amplified loci can vary among the source and the target species. Many of the WS polymorphic microsatellites are reported to be monomorphic in Jabiru Stork (JS) (Jabiru mycteria) [13]. While WSµ03 and WSµ20 are monomorphic in PS, in JS and European White Stork (EWS) (Ciconia ciconia) they are reported to be polymorphic [13,14]. Loci WSµ08 and WSµ24, which are polymorphic in WS, are monomorphic in both PS and JS [13]. DNA replication slippage is considered as a major source of microsatellite variability [15]. Except for WSµ17, variations in allele sizes of polymorphic loci of PS are in agreement with the size of repeat motifs described in WS [6]. WSµ17 was described as a tetranucleotide in WS, but in PS we also observed an allele which varies in length from others by 2 bp. Although WSµ08 shows two alleles of size difference of around 1 bp, we counted it as monomorphic because in WS it was described as dinucleotide marker [6]. Such size differences could be because of mutations outside the stretch of repeat strand. WSµ09, WSµ13, WSµ18 and WSµ23 were as expected tri-, di-, tetra-and tetranucleotide microsatellites, respectively. The number of alleles over four polymorphic loci in PS is more than the number reported over same loci in WS.
Although all of the reported WS microsatellites are known for 100% cross-amplification success in congeneric PS and Milky Stork (Mycteria cinerea) [3,16], the reasons of high amplification failures of some loci could be explained in terms of general problems of insufficient and degraded templates from non-invasive samples [3,17] or in terms of possibility of null alleles due to mutations in primer binding sites in certain lines [18]. Due to the high probability of amplification failure because of the poor DNA quantity and quality from non-invasive samples, the possibility of null alleles could not be tested. Even for the samples which failed to amplify the loci on first attempt, the PCR experiment was performed repeatedly (3-4 times) with minute changes in the concentrations of the reagents and cycling conditions. In some cases the amplification success was seen on repeated Table 2  PCR experiments. Five samples were repeatedly genotyped for 3 times and each time the similar fragment size was reported, which indicated a 0% scoring error. Maximum missing data for WSµ17 could be the reason for its significant deviations from HWE. LD between WSµ17 and WSµ18 was also reported in WS [6] and EWS [14], which indicates that in stork species these two loci are closely linked.
In WS, most of the studies report low-moderate genetic diversity in colonies of Americas [6,19,20], which are comparable with the 3 microsatellite limited genetic diversity encountered in PS of Delhi Zoo. As per the general trend of genetic diversity for waterbirds [21], lowmoderate genetic diversity in storks could be because of being inhabitants of wetland habitat and their larger body mass. Some other factors have also been attributed to low genetic diversity in other storks. For instance, in the case of EWS it has been attributed to historical bottlenecks [14].
Field and molecular studies in WS colonies have indicated a lack of differentiation and low natal site philopatry [19,20,[22][23][24][25][26][27]. Although the assignment test and differentiation results of our study indicate a lack of genetic structuring, due to low number of loci used, the presence of similar genotypes between the nesting seasons is not a strong reason to conclude that it is indeed the same individuals which are visiting the zoo year after year. Field observations indicate that in any given season PS can fluctuate between the short distance separated colonies [28]. So to be more certain if PS display nesting site fidelity, there is a need to test the genetic connectivity of multiple, spatially and temporally distributed PS colonies. Any high degree of fidelity in PS of Delhi Zoo, if indeed such is the case, may be on account of the fact that the zoo premises continue to offer a relatively safe environment for nesting [1].
Being an exclusively piscivorous bird [1], PS is a flagship of wetland ecosystems which are gravely threatened all across the Indian subcontinent. Its wide geographical spread and capacity to disperse far away from natal sites warrants further studies along biogeographical lines. Molecular genetic studies are likely to be extremely relevant in this regard, besides providing important baseline data for conservation of this species [2].

Limitations of study
The limitations of the present study are the exclusive use of non-invasive samples and low number of useful polymorphic microsatellites. The fact that around 25% samples showed high PCR failures at most of the loci, severely limited our study. Except for some initial studies [19,20], the number of microsatellites used in most of the genetic diversity studies on WS are considerably larger [29] than what we have managed. Higher number of polymorphic loci results an increase in confidence that similar genotypes are from same individual. For further improvements, we recognize the need to develop species specific markers of PS.