Optimizing the use of shed feathers for genetic analysis

Large-scale genotypic identification reveals density-dependent natal dispersal patterns in an elusive bird of prey

BACKGROUND

Natal dispersal, the distance between site of birth and site of first breeding, has a fundamental role in population dynamics and species' responses to environmental changes. Population density is considered a key driver of natal dispersal. However, few studies have been able to examine densities at both the natal and the settlement site, which is critical for understanding the role of density in dispersal. Additionally, the role of density on natal dispersal remains poorly understood in long-lived and slowly reproducing species, due to their prolonged dispersal periods and often elusive nature. We studied the natal dispersal of the white-tailed eagle (Haliaeetus albicilla) in response to local breeder densities. We investigated the effects of the number of active territories around the natal site on (a) natal dispersal distance and (b) the difference between natal and settlement site breeder density. We were interested in whether eagles showed tendencies of conspecific attraction (positive density-dependence) or intraspecific competition (negative density-dependence) and how this related to settlement site breeder density.

METHODS

We used a combination of long-term visual and genotypic identification to match individuals from their breeding site to their natal nest. We identified natal dispersal events for 355 individuals hatched between 1984 and 2015 in the Baltic Sea coast and Arctic areas of Finland. Of those, 251 were identified by their genotype.

RESULTS

Individuals born in high-density areas dispersed shorter distances than those born in low-density areas, but settled at lower density breeding sites in comparison to their natal site. Eagles born in low natal area densities dispersed farther but settled in higher density breeding sites compared to their natal site.

CONCLUSIONS

We show that eagles might be attracted by conspecifics (positive density-dependence) to identify high-quality habitats or find mates, but do not settle in the most densely populated areas. This indicates that natal dispersal is affected by an interplay of conspecific attraction and intraspecific competition, which has implications for population dynamics of white-tailed eagles, but also other top predators. Furthermore, our study demonstrates the value of long-term collection of both nestling and (non-invasive) adult DNA samples, and thereafter using genotype matching to identify individuals in long-lived and elusive species.

Effect of storage temperature and duration on direct PCR amplification of various feather types and DBS matrices

The use of direct PCR has been pioneered over the last decade for DNA analysis of biological specimens of distinct origins. The information on how longer these specimens can be stored and amplified by direct PCR is however scanty. Such a piece of information could expedite research and diagnostic studies without compromising the reliability of results. The current study was therefore designed to analyze the effect of storage temperature and duration on direct PCR amplification of biological specimens having either low quantity or high quantity of DNA. Whole blood, dried blood spots (DBS), and feathers from chicken were stored for five years at three different temperatures, viz. room temperature (∼25 °C), 4 °C, and -20 °C. These samples were subjected to crude DNA extraction by diluting them in PBS buffer and heating at 98 °C after 1 day, 7 days, 15 days, 1 month, 3 months, 6 months, 1 year, 3 years and 5 years of storage. The crude DNA was PCR-amplified with the use of DNA sexing primers as well as DNA barcoding primers. Incubation at 98 °C for 10 min of any type of sample in PBS buffer was sufficient for crude DNA extraction. There was irrelevant impact of feather type, DBS matrix nature and storage temperature on amplification success over the period of analysis. It was possible to successfully accomplish the amplification of 96 samples with the use of routine PCR reagents within 3.5-6.0 hrs. In short, economical and fast genetic analysis of commonly used avian samples is feasible after their storage for longer time at room temperature.

Leveraging genomics to understand threats to migratory birds

Understanding how risk factors affect populations across their annual cycle is a major challenge for conserving migratory birds. For example, disease outbreaks may happen on the breeding grounds, the wintering grounds, or during migration and are expected to accelerate under climate change. The ability to identify the geographic origins of impacted individuals, especially outside of breeding areas, might make it possible to predict demographic trends and inform conservation decision-making. However, such an effort is made more challenging by the degraded state of carcasses and resulting low quality of DNA available. Here, we describe a rapid and low-cost approach for identifying the origins of birds sampled across their annual cycle that is robust even when DNA quality is poor. We illustrate the approach in the common loon (Gavia immer), an iconic migratory aquatic bird that is under increasing threat on both its breeding and wintering areas. Using 300 samples collected from across the breeding range, we develop a panel of 158 single-nucleotide polymorphisms (SNP) loci with divergent allele frequencies across six genetic subpopulations. We use this SNP panel to identify the breeding grounds for 142 live nonbreeding individuals and carcasses. For example, genetic assignment of loons sampled during botulism outbreaks in parts of the Great Lakes provides evidence for the significant role the lakes play as migratory stopover areas for loons that breed across wide swaths of Canada, and highlights the vulnerability of a large segment of the breeding population to botulism outbreaks that are occurring in the Great Lakes with increasing frequency. Our results illustrate that the use of SNP panels to identify breeding origins of carcasses collected during the nonbreeding season can improve our understanding of the population-specific impacts of mortality from disease and anthropogenic stressors, ultimately allowing more effective management.

A DNA toolbox for non-invasive genetic studies of sambar deer (Rusa unicolor)

Invasive sambar deer (Rusa unicolor) are having significant detrimental impacts on natural environments in south-eastern Australia. Little, however, is known about their ecology, limiting evidence-based management strategies directed at reducing deer impacts. Genetic data, generated from DNA isolated from deer scats, can be used to fill ecological knowledge gaps. This study outlines a non-invasive genetic sampling strategy by which good-quality DNA from a single deer scat can be used to determine (1) species of origin, (2) sex and (3) a unique DNA profile. DNA from deer tissue and sambar deer scat samples were used to develop and optimise molecular methods to collect reliable genetic information. A DNA toolbox is presented that describes how to find, collect and store scat samples, isolate DNA and use molecular markers to generate informative genetic data. Generating genetic data using this approach will support studies aimed at acquiring ecological knowledge about sambar deer. Such knowledge will be critical for developing evidence-based recommendations to improve on-ground management decisions for sambar deer.

A novel method to optimise the utility of underused moulted plumulaceous feather samples for genetic analysis in bird conservation

Non-invasive sampling methods are increasingly being used in conservation research as they reduce or eliminate the stress and disturbance resulting from invasive sampling of blood or tissue. Here we present a protocol optimised for obtaining usable genetic material from moulted plumulaceous feather samples. The combination of simple alterations to a ‘user-developed’ method, comprised of increased incubation time and modification of temperature and volume of DNA elution buffer, are outlined to increase DNA yield and significantly increase DNA concentration (W = 81, p < 0.01, Cohens’s d = 0.89). We also demonstrate that the use of a primerless polymerase chain reaction (PCR) technique increases DNA quality and amplification success when used prior to PCR reactions targeting avian mitochondrial DNA (mtDNA). A small amplicon strategy proved effective for mtDNA amplification using PCR, targeting three overlapping 314–359 bp regions of the cytochrome oxidase I barcoding region which, when combined, aligned with target-species reference sequences. We provide evidence that samples collected non-invasively in the field and kept in non-optimal conditions for DNA extraction can be used effectively to sequence a 650 bp region of mtDNA for genetic analysis.

DNA based method for determining source country of the short beaked echidna (Tachyglossus aculeatus) in the illegal wildlife trade

The illegal trade in wild animals being sold as 'captive bred' is an emerging issue in the pet and zoo industry and has both animal welfare and conservation implications. DNA based methods can be a quick, inexpensive, and definitive way to determine the source of these animals, thereby assisting efforts to combat this trade. The short beaked echidna (Tachyglossus aculeatus) is currently one of the species suspected to be targeted in this trade. As this species is distributed throughout Australia and in New Guinea (currently comprising of five recognised sub-species), this project aimed to develop a DNA based method to definitively determine the source country of an echidna and explore the use of non-invasive sampling techniques. Here we use non-invasively sampled echidna quills and demonstrate the extraction of mitochondrial DNA and amplification of a region of the mitochondrial genome. Phylogenetically informative markers for analysis of a 322bp segment of the D-loop region were developed, and subsequently validated, using animals with known source locations allowing us to reliably distinguish between echidnas from New Guinea, and Australia. This research presents the first validated forensic protocols for short beaked echidnas and will be an integral tool in understanding the movement of animals in this emerging trade.

Integrated tool for microsatellite isolation and validation from the reference genome and their application in the study of breeding turnover in an endangered avian population

Accurate individual identification is required to estimate survival rates in avian populations. For endangered species, non-invasive methods of obtaining individual identification, such as using molted feathers as a source of DNA for microsatellite markers, are preferred because of less disturbance, easy sample preparation and high efficiency. With the availability of many avian genomes, a few pipelines isolating genome-wide microsatellites have been published, but it is still a challenge to isolate microsatellites from the reference genome efficiently. Here, we have developed an integrated tool comprising a bioinformatic pipeline and experimental procedures for microsatellite isolation and validation based on the reference genome. We have identified over 95 000 microsatellite loci and established a system comprising 10 highly polymorphic markers (PIC value: 0.49-0.93, mean: 0.79) for an endangered species, saker falcon (Falco cherrug). These markers (except 1) were successfully amplified in 126 molted feathers, exhibiting high amplification success rates (83.9-99.7%), high quality index (0.90-0.97) and low allelic dropout rates (1-9.5%). To further assess the efficiency of this marker system in a population study, we identified individual sakers using these molted feathers (adult) and 146 plucked feathers (offspring). The use of parent and offspring samples enabled us to infer the genotype of missing samples (N = 28), and all adult genotypes were used to ascertain that breeding turnover is a useful proxy for survival estimation in sakers. Our study presents a cost-effective tool for microsatellite isolation based on publicly available reference genomes and demonstrates the power of this tool in estimating key parameters of avian population dynamics.

Estimating demographic contributions to effective population size in an age-structured wild population experiencing environmental and demographic stochasticity

A population's effective size (N_e ) is a key parameter that shapes rates of inbreeding and loss of genetic diversity, thereby influencing evolutionary processes and population viability. However, estimating N_e , and identifying key demographic mechanisms that underlie the N_e to census population size (N) ratio, remains challenging, especially for small populations with overlapping generations and substantial environmental and demographic stochasticity and hence dynamic age-structure. A sophisticated demographic method of estimating N_e /N, which uses Fisher's reproductive value to account for dynamic age-structure, has been formulated. However, this method requires detailed individual- and population-level data on sex- and age-specific reproduction and survival, and has rarely been implemented. Here, we use the reproductive value method and detailed demographic data to estimate N_e /N for a small and apparently isolated red-billed chough (Pyrrhocorax pyrrhocorax) population of high conservation concern. We additionally calculated two single-sample molecular genetic estimates of N_e to corroborate the demographic estimate and examine evidence for unobserved immigration and gene flow. The demographic estimate of N_e /N was 0.21, reflecting a high total demographic variance (σ2dg) of 0.71. Females and males made similar overall contributions to σ2dg. However, contributions varied among sex-age classes, with greater contributions from 3 year-old females than males, but greater contributions from ≥5 year-old males than females. The demographic estimate of N_e was ~30, suggesting that rates of increase of inbreeding and loss of genetic variation per generation will be relatively high. Molecular genetic estimates of N_e computed from linkage disequilibrium and approximate Bayesian computation were approximately 50 and 30, respectively, providing no evidence of substantial unobserved immigration which could bias demographic estimates of N_e . Our analyses identify key sex-age classes contributing to demographic variance and thus decreasing N_e /N in a small age-structured population inhabiting a variable environment. They thereby demonstrate how assessments of N_e can incorporate stochastic sex- and age-specific demography and elucidate key demographic processes affecting a population's evolutionary trajectory and viability. Furthermore, our analyses show that N_e for the focal chough population is critically small, implying that management to re-establish genetic connectivity may be required to ensure population viability.

Cooled Propylene Glycol as a Pragmatic Choice for Preservation of DNA From Remote Field-Collected Diptera for Next-Generation Sequence Analysis

Next-generation sequencing (NGS)-based methods can now be applied to large population-scale studies, but this demands very high-quality DNA. For specimens collected from remote field locations, DNA degradation can be a problem, requiring logistically challenging preservation techniques. Simpler preservation techniques are therefore required. Prior to collection of exotic fruit fly (Tephritidae) species, a number of readily available preservatives with storage at either 4°C or room temperature were trialed here to determine the DNA quality for three locally available Diptera species, Fannia canicularis (L.), Musca domestica L., and Lucilia sericata Meigen. Considerable variation was observed between the different preservatives, species, and temperatures, but several preservatives at 4°C were favored. Chilled propylene glycol was subsequently used for the storage and carriage of Australian field-collected Bactrocera fruit fly specimens to New Zealand. When processed up to 20 d later, DNA fragments of ∼10-20 kb were obtained for successful genotyping by sequencing analysis. This protocol is therefore recommended as a logistically simple and safe approach for distant collection of dipteran samples for NGS population genomic studies.

Evidence of the phenotypic expression of a lethal recessive allele under inbreeding in a wild population of conservation concern

Deleterious recessive alleles that are masked in outbred populations are predicted to be expressed in small, inbred populations, reducing both individual fitness and population viability. However, there are few definitive examples of phenotypic expression of lethal recessive alleles under inbreeding conditions in wild populations. Studies that demonstrate the action of such alleles, and infer their distribution and dynamics, are required to understand their potential impact on population viability and inform management responses. The Scottish population of red-billed choughs (Pyrrhocorax pyrrhocorax), which currently totals <60 breeding pairs and is of major conservation concern, has recently been affected by lethal blindness in nestlings. We used family data to show that the pattern of occurrence of blindness within and across affected families that produced blind nestlings was exactly 0·25, matching that expected given a single-locus autosomal lethal recessive allele. Furthermore, the observed distribution of blind nestlings within affected families did not differ from that expected given Mendelian inheritance of such an allele. Relatedness estimates showed that individuals from affected families were not more closely related to each other than they were to individuals from unaffected families that did not produce blind nestlings. Blind individuals tended to be less heterozygous than non-blind individuals, as expected if blindness was caused by the expression of a recessive allele under inbreeding. However, there was no difference in the variance in heterozygosity estimates, suggesting that some blind individuals were relatively outbred. These results suggest carriers of the blindness allele may be widely distributed across contemporary families rather than restricted to a single family lineage, implying that the allele has persisted across multiple generations. Blindness occurred at low frequency (affecting 1·6% of observed nestlings since 1981). However, affected families had larger initial brood sizes than unaffected families. Such high fecundity of carriers of a lethal recessive allele might reflect overdominance, potentially reducing purging and increasing allele persistence probability. We thereby demonstrate the phenotypic expression of a lethal recessive allele in a wild population of conservation concern, and provide a general framework for inferring allele distribution and persistence and informing management responses.

SkydancerPlex: A novel STR multiplex validated for forensic use in the hen harrier (Circus cyaneus)

The hen harrier (Circus cyaneus) is a bird of prey which is heavily persecuted in the UK because it preys on the game bird red grouse (Lagopus lagopus scoticus). To help investigations into illegal killings of hen harrier, a STR multiplex kit containing eight short tandem repeat (STR) markers and a chromohelicase DNA binding protein 1 (CHD 1) sexing marker was developed. The multiplex kit was tested for species specificity, sensitivity, robustness, precision, accuracy and stability. Full profiles were obtained with as little as 0.25 ng of template DNA. Concurrent development of an allelic ladder to ensure reliable and accurate allele designation across laboratories makes the SkydancerPlex the first forensic DNA profiling system in a species of wildlife to be fully validated according to SWGDAM and ISFG recommendations. An average profile frequency of 3.67 × 10(-8), a PID estimate of 5.3 × 10(-9) and a PID-SIB estimate of 9.7 × 10(-4) make the SkydancerPlex an extremely powerful kit for individualisation.

The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica)

Landscape genomics promises to provide novel insights into how neutral and adaptive processes shape genome-wide variation within and among populations. However, there has been little emphasis on examining whether individual-based phenotype-genotype relationships derived from approaches such as genome-wide association (GWAS) manifest themselves as a population-level signature of selection in a landscape context. The two may prove irreconcilable as individual-level patterns become diluted by high levels of gene flow and complex phenotypic or environmental heterogeneity. We illustrate this issue with a case study that examines the role of the highly prevalent gastrointestinal nematode Trichostrongylus tenuis in shaping genomic signatures of selection in red grouse (Lagopus lagopus scotica). Individual-level GWAS involving 384 SNPs has previously identified five SNPs that explain variation in T. tenuis burden. Here, we examine whether these same SNPs display population-level relationships between T. tenuis burden and genetic structure across a small-scale landscape of 21 sites with heterogeneous parasite pressure. Moreover, we identify adaptive SNPs showing signatures of directional selection using F(ST) outlier analysis and relate population- and individual-level patterns of multilocus neutral and adaptive genetic structure to T. tenuis burden. The five candidate SNPs for parasite-driven selection were neither associated with T. tenuis burden on a population level, nor under directional selection. Similarly, there was no evidence of parasite-driven selection in SNPs identified as candidates for directional selection. We discuss these results in the context of red grouse ecology and highlight the broader consequences for the utility of landscape genomics approaches for identifying signatures of selection.

Genome-wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco-evolutionary processes in natural populations. Host-parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade-offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome-wide SNPs. We first use genome-wide association to identify individual SNPs associated with nematode burden. We then partition genome-wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome-wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large-effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).

Self-care behaviors of adults with type 2 diabetes mellitus in Greece

The purpose of this study was to examine self-care behaviors of adults with type 2 diabetes mellitus living in the Metropolitan Area of Thessaloniki in Northern Greece. The Summary of Diabetes Self-Care behaviors measurement was administered to 215 patients, out of which 177 were eligible to participate (87 males). Patients, aged 30 years or more, were recruited through a university hospital day-clinic. Older patients (>65 years), as well as those with "higher educational level" did not distribute their daily carbohydrate intake equally. Nevertheless, they were more likely to adapt to their physician's recommendations regarding medication and to regularly perform suggested blood glucose checking. Exercise patterns were more often found for higher educated, earlier diagnosed males. Younger patients were less likely to follow their healthcare professional's recommendations, regarding diet, medication intake, blood glucose checking, foot care and exercise compared to older patients. These results pose a higher risk for complications and morbidity in younger patients with type 2 diabetes mellitus, who most possibly will require intensive treatment in the future.

Fine-scale population epigenetic structure in relation to gastrointestinal parasite load in red grouse (Lagopus lagopus scotica)

Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north-east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome-wide and locus-specific methylation states. Following methylation-sensitive AFLP (MSAP), 129 bands, representing 73 methylation-susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome-wide methylation levels and were also significantly epigenetically (F_SC = 0.0227; P < 0.001) and genetically (F_SC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome-wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus-specific methylation states (F_ST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam) and Bayesian environmental analysis (bayenv2). Following Sanger sequencing, genome mapping and geneontology (go) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite-associated cytosine methylation.

Techniques for optimising long-term, large-scale capture–mark–resighting raptor studies: climb hard, band fast

Context Efficiency of large-scale capture–mark–recapture (CMR) studies can be improved by developing accurate methods for predicting the window of opportunity in which banding can occur. Aims This research aimed to investigate potential efficiency improvements in a long-term CMR raptor study. The research focussed on: (1) developing selection processes for adopting CMR protocols; (2) testing methods for increasing the number of nestlings successfully banded; and (3) assessing the efficacy of visual identification (VID) bands for collecting re-sight data. Methods Ten selection criteria were developed into a robust CMR-technique selection process and used to assess marking techniques commonly applied to birds. Optimising banding effort by predicting banding dates using two different techniques a priori and a posteriori were tested against a traditional approach to the timing of banding. The cost (in time) to collect resight data at an active nest site was also measured. Key results The CMR selection criteria and parameters provided a transparent selection process and scored metal VID bands the highest for the study design. This provided individual recognition of marked birds up to the expected life-span of 14 years. Both techniques for predicting banding dates improved the proportion of whole clutches banded by 40%. The average time to identify both peregrine falcon adults of a breeding pair wearing VID bands was 30 min. Conclusions The two methods described here for predicting preferred banding dates are of particular value as efficient approaches to banding large numbers of nestlings are key to the success of CMR studies. All of the methods developed in this research can be applied to CMR studies of almost any bird species with a predictable seasonal breeding system. Implications Optimisation and cost effectiveness of CMR studies for seasonal breeding birds can be significantly improved by accurately predicting the window of opportunity in which banding of nestlings can be carried out, and also utilising VID colour bands for rapid collection of recapture data.

Identification of protected avian species using a single feather barb

We report on the unambiguous identification of protected avian species from as little as one barb of a feather. Many avian species are protected by international agreements and national legislation, yet they are traded illegally because of their high value. Two sections of the avian mitochondrial genome were chosen to identify bird species, being a 561-bp section of ND2 gene and a 921-bp section of the ND5 gene. Two different DNA extraction methods were compared for their ability to reliably isolate sufficient DNA to be detected in a subsequent PCR. Using a commercial kit supplied by QIAGEN, a complete sequence was obtained from one barb for the ND2 gene, whereas two barbs were required to reliably sequence the 921-bp section of the ND5 gene. The process worked on all species tested using feathers from archival museum specimens, resulted in minimal damage to the specimen and can readily be adopted by a forensic science laboratory.

Feather barbs as a good source of mtDNA for bird species identification in forensic wildlife investigations

Background

The ability to accurately identify bird species is crucial for wildlife law enforcement and bird-strike investigations. However, such identifications may be challenging when only partial or damaged feathers are available for analysis.

Results

By applying vigorous contamination controls and sensitive PCR amplification protocols, we found that it was feasible to obtain accurate mitochondrial (mt)DNA-based species identification with as few as two feather barbs. This minimally destructive DNA approach was successfully used and tested on a variety of bird species, including North American wild turkey (Meleagris gallopavo), Canada goose (Branta canadensis), blue heron (Ardea herodias) and pygmy owl (Glaucidium californicum). The mtDNA was successfully obtained from 'fresh' feathers, historic museum specimens and archaeological samples, demonstrating the sensitivity and versatility of this technique.

Conclusions

By applying appropriate contamination controls, sufficient quantities of mtDNA can be reliably recovered and analyzed from feather barbs. This previously overlooked substrate provides new opportunities for accurate DNA species identification when minimal feather samples are available for forensic analysis.

Description of microsatellite markers and genotyping performances using feathers and buccal swabs for the Ivory gull (Pagophila eburnea)

We report 22 new polymorphic microsatellites for the Ivory gull (Pagophila eburnea), and we describe how they can be efficiently co-amplified using multiplexed polymerase chain reactions. In addition, we report DNA concentration, amplification success, rates of genotyping errors and the number of genotyping repetitions required to obtain reliable data with three types of noninvasive or nondestructive samples: shed feathers collected in colonies, feathers plucked from living individuals and buccal swabs. In two populations from Greenland (n=21) and Russia (Severnaya Zemlya Archipelago, n=21), the number of alleles per locus varied between 2 and 17, and expected heterozygosity per population ranged from 0.18 to 0.92. Twenty of the markers conformed to Hardy-Weinberg and linkage equilibrium expectations. Most markers were easily amplified and highly reliable when analysed from buccal swabs and plucked feathers, showing that buccal swabbing is a very efficient approach allowing good quality DNA retrieval. Although DNA amplification success using single shed feathers was generally high, the genotypes obtained from this type of samples were prone to error and thus need to be amplified several times. The set of microsatellite markers described here together with multiplex amplification conditions and genotyping error rates will be useful for population genetic studies of the Ivory gull.