Using blocks of linked single nucleotide polymorphisms as highly polymorphic genetic markers for parentage analysis

Single nucleotide polymorphisms (SNPs) are plentiful in most genomes and amenable to high throughput genotyping, but they are not yet popular for parentage or paternity analysis. The markers are bi-allelic, so individually they contain little information about parentage, and in nonmodel organisms the process of identifying large numbers of unlinked SNPs can be daunting. We explore the possibility of using blocks of between three and 26 linked SNPs as highly polymorphic molecular markers for reconstructing male genotypes in polyandrous organisms with moderate (five offspring) to large (25 offspring) clutches of offspring. Haplotypes are inferred for each block of linked SNPs using the programs Haplore and Phase 2.1. Each multi-SNP haplotype is then treated as a separate allele, producing a highly polymorphic, 'microsatellite-like' marker. A simulation study is performed using haplotype frequencies derived from empirical data sets from Drosophila melanogaster and Mus musculus populations. We find that the markers produced are competitive with microsatellite loci in terms of single parent exclusion probabilities, particularly when using six or more linked SNPs to form a haplotype. These markers contain only modest rates of missing data and genotyping or phasing errors and thus should be seriously considered as molecular markers for parentage analysis, particularly when the study is interested in the functional significance of polymorphisms across the genome.     

Double‐digest RAD sequencing outperforms microsatellite loci at assigning paternity and estimating relatedness: A proof of concept in a highly promiscuous bird

Information on genetic relationships among individuals is essential to many studies of the behaviour and ecology of wild organisms. Parentage and relatedness assays based on large numbers of single nucleotide polymorphism (SNP) loci hold substantial advantages over the microsatellite markers traditionally used for these purposes. We present a double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) analysis pipeline that, as such, simultaneously achieves the SNP discovery and genotyping steps and which is optimized to return a statistically powerful set of SNP markers (typically 150–600 after stringent filtering) from large numbers of individuals (up to 240 per run). We explore the trade‐offs inherent in this approach through a set of experiments in a species with a complex social system, the variegated fairy‐wren (Malurus lamberti) and further validate it in a phylogenetically broad set of other bird species. Through direct comparisons with a parallel data set from a robust panel of highly variable microsatellite markers, we show that this ddRAD‐seq approach results in substantially improved power to discriminate among potential relatives and considerably more precise estimates of relatedness coefficients. The pipeline is designed to be universally applicable to all bird species (and with minor modifications to many other taxa), to be cost‐ and time‐efficient, and to be replicable across independent runs such that genotype data from different study periods can be combined and analysed as field samples are accumulated.     

The fitness consequences of kin-biased dispersal in a cooperatively breeding bird

Cooperative alliances among kin may not only lead to indirect fitness benefits for group-living species, but can also provide direct benefits through access to mates or higher social rank. However, the immigrant sex in most species loses any potential benefits of living with kin unless immigrants disperse together or recruit relatives into the group in subsequent years. To look for evidence of small subgroups of related immigrants within social groups (kin substructure), we used microsatellites to assess relatedness between immigrant females of the cooperatively breeding superb starling, Lamprotornis superbus. We determined how timing of immigration led to kin subgroup formation and if being part of one influenced female fitness. Although mean relatedness in groups was higher for males than females, 26% of immigrant females were part of a kin subgroup with a sister. These immigrant sibships formed through kin recruitment across years more often than through coalitions immigrating together in the same year. Furthermore, females were more likely to breed when part of a kin subgroup than when alone, suggesting that female siblings form alliances that may positively influence their fitness. Ultimately, kin substructure should be considered when determining the role of relatedness in the evolution of animal societies.     

Single nucleotide polymorphisms for pig identification and parentage exclusion

Single nucleotide polymorphisms (SNPs) have become an important type of marker for commercial diagnostic and parentage genotyping applications as automated genotyping systems have been developed that yield accurate genotypes. Unfortunately, allele frequencies for public SNP markers in commercial pig populations have not been available. To fulfil this need, SNP markers previously mapped in the USMARC swine reference population were tested in a panel of 155 boars that were representative of US purebred Duroc, Hampshire, Landrace and Yorkshire populations. Multiplex assay groups of 5-7 SNP assays/group were designed and genotypes were determined using Sequenom's massarray system. Of 80 SNPs that were evaluated, 60 SNPs with minor allele frequencies >0.15 were selected for the final panel of markers. Overall identity power across breeds was 4.6 x 10(-23), but within-breed values ranged from 4.3 x 10(-14) (Hampshire) to 2.6 x 10(-22) (Yorkshire). Parentage exclusion probability with only one sampled parent was 0.9974 (all data) and ranged from 0.9594 (Hampshire) to 0.9963 (Yorkshire) within breeds. Sire exclusion probability when the dam's genotype was known was 0.99998 (all data) and ranged from 0.99868 (Hampshire) to 0.99997 (Yorkshire) within breeds. Power of exclusion was compared between the 60 SNP and 10 microsatellite markers. The parental exclusion probabilities for SNP and microsatellite marker panels were similar, but the SNP panel was much more sensitive for individual identification. This panel of SNP markers is theoretically sufficient for individual identification of any pig in the world and is publicly available.     

Probability of random sire exclusion using microsatellite markers for parentage verification

Many microsatellite sequences have been described in the bovine genome. Being highly polymorphic these have been suggested as markers for parentage verification and individual identification in cattle. We have evaluated the use of five highly polymorphic microsatellite markers for parentage verification in 14 breeds of cattle in the UK. Three of the microsatellite loci occur within introns in genes: BoLA DRB3 , steroid 21-hydroxylase, and the beta subunit of the follicle-stimulating hormone. The other two are anonymous sites ETH131 and HEL6. Results were analysed by a statistical approach that takes in to account deviations from Hardy-Wienberg equilibrium and linkage disequilibrium for multiple loci. The method of determining the probability of random sire exclusion uses observed genotype frequencies instead of allele frequencies. Independently, the markers used have a probability of between 0.72 and 0.62 of identifying a parentage error, while used together the five markers give, on average across breeds, a probability of 0.99 of excluding an incorrect sire.     

Validation of a single‐nucleotide polymorphism panel for parentage testing of farmed red deer

Red deer behaviour prevents the accurate physical matching of dams to fawns and, therefore, breeders have to rely on DNA‐based parentage testing. A panel consisting of 100 single nucleotide polymorphism markers, with an average minor allele frequency of 0.25, was able to resolve 92% of fawns to both parents. In comparison, an existing 12‐marker microsatellite panel was able to resolve 68% of fawns to both parents. When excluding dam DNA information, the single nucleotide polymorphism panel matched 81% of the fawns to their sires and the microsatellite panel 71%.     

Estimation of the number of SNP genetic markers required for parentage verification

Formulae were developed to compute exclusion probabilities for parentage confirmation for any number of diallelic markers under the assumption that the minor allele frequency (MAF) varied among markers, but has a uniform distribution. Three scenarios were analysed: a progeny with (1) a single putative parent; (2) two putative parents; and (3) one actual parent and one putative parent. Exclusion probabilities were computed for minimum values for the MAFs of 0.1, 0.2 and 0.3, and required either one or at least two conflicts for exclusion. The numbers of markers required to obtain 99% exclusion probabilities based on a single conflict for the three minimum MAFs were 54, 45 and 39 for scenario 1; 17, 16 and 15 for scenario 2; and 28, 25 and 24 for scenario 3. The requirement of at least two conflicts for exclusion increased the number of markers required by approximately 45% for all three scenarios and all three minimum MAFs. The results obtained by the analytical formulae were very close to results obtained by simulation and to values in the literature for specific marker sets.     

Mating system, philopatry and patterns of kinship in the cooperatively breeding subdesert mesite Monias benschi

In the first molecular study of a member of the threatened avian family, Mesitornithidae, we used nine polymorphic microsatellite loci to elucidate parentage, patterns of within-group kinship and occurrence of extra-group paternity in the subdesert mesite Monias benschi, of southwest Madagascar. We found this cooperatively breeding species to have a very fluid mating system. There was evidence of genetic monogamy and polygynandry: of the nine groups with multiple offspring, six contained one breeding pair with unrelated helpers and three contained multiple male and female breeders with related helpers. Although patterns of within-group kinship varied, there was a strong positive relationship between group size and relatedness, suggesting that groups form by natal philopatry. There was also a strong positive correlation between within-sex and between-sex relatedness, indicating that unlike most cooperatively breeding birds, philopatry involved both sexes. In contrast to predictions of kin selection and reproductive skew models, all monogamous groups contained unrelated individuals, while two of the three polygynandrous groups were families. Moreover, although between-group variation in seasonal reproductive success was related to within-group female relatedness, relatedness among males and between the sexes had no bearing on a group's reproductive output. While kin selection may underlie helping behaviour in females, factors such as direct long-term fitness benefits of group living probably determine helping in males. Of the 14 offspring produced by fully sampled groups, at least two were sired by males from neighbouring groups: one by a breeding male and one by a nonbreeding male, suggesting that males may augment their reproductive success through extra-group paternity.     

Characterization of eight single nucleotide polymorphism markers in Aedes aegypti

Sequencing of part of seven genes from Aedes aegypti collected in 16 Brazilian cities revealed the existence of 53 single nucleotide polymorphisms (SNPs), representing one SNP every 52 base pairs. From these 53 SNPs, we selected eight that are independent and highly polymorphic. We describe the use of these markers for differentiation of Brazilian populations of A. aegypti. These are the first SNPs developed for delineating population structure in A. aegypti, and will be a useful complement to epidemiological studies.     

Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High‐throughput high‐density genotyping arrays continue to be a fast, accurate, and cost‐effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high‐density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole‐genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives—honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.     

Characterization of 15 single nucleotide polymorphism markers for chimpanzees (Pan troglodytes)

We report the characterization of 15 new single nucleotide polymorphism markers for a threatened species, the chimpanzee (Pan troglodytes), developed using a targeted gene approach. These markers are derived from the Y chromosome and autosomal regions of the genome and show frequency differences between chimpanzee subspecies from central and western Africa. These single nucleotide polymorphism markers are the first to be designed for the genotyping of wild chimpanzee populations and will provide a useful addition to the genetic tools employed for the conservation management of this threatened species.     

Exclusion probabilities of 22 bovine microsatellite markers in fluorescent multiplexes for semiautomated parentage testing

Six multiplexes developed for semiautomated fluorescence genotyping were evaluated for parentage testing. These multiplexes contained primer pairs for the amplification of 22 microsatellites on 17 bovine autosomes. Exclusion probabilities were determined from genotypes of 1022 Holstein cattle and 311 beef cattle belonging to five breeds. Two cases were considered: case 1, genotypes are known for an alleged parent and an offspring but genotypes of a confirmed parent are unknown; and case 2, genotypes are known for an alleged parent, a confirmed parent and an offspring. If the alleged parent is not the true parent, then the 22 markers will exclude the alleged parent with a probability of >0·9986 for case 1 and with a probability of >0·99999 for case 2. On the basis of these exclusion probabilities, the probability that an alleged parent will be falsely included as a parent is in the range of 1/716 to 1/2845 for case 1 and 1/1·2 million to 1/159753 for case 2. In addition to these results, a rapid and efficient non-organic method for extraction of DNA from semen is described.     

Multi-locus inference of population structure: a comparison between single nucleotide polymorphisms and microsatellites

Although growing numbers of single nucleotide polymorphisms (SNPs) and microsatellites (short tandem repeat polymorphisms or STRPs) are used to infer population structure, their relative properties in this context remain poorly understood. SNPs and STRPs mutate differently, suggesting multi-locus genotypes at these loci might differ in ability to detect population structure. Here, we use coalescent simulations to measure the power of sets of SNPs and STRPs to identify population structure. To maximize the applicability of our results to empirical studies, we focus on the popular STRUCTURE analysis and evaluate the role of several biological and practical factors in the detection of population structure. We find that: (1) fewer unlinked STRPs than SNPs are needed to detect structure at recent divergence times <0.3 N_e generations; (2) accurate estimation of the number of populations requires many fewer STRPs than SNPs; (3) for both marker types, declines in power due to modest gene flow (N_em=1.0) are largely negated by increasing marker number; (4) variation in the STRP mutational model affects power modestly; (5) SNP haplotypes (θ=1, no recombination) provide power comparable with STRP loci (θ=10); (6) ascertainment schemes that select highly variable STRP or SNP loci increase power to detect structure, though ascertained data may not be suitable to other inference; and (7) when samples are drawn from an admixed population and one of its parent populations, the reduction in power to detect two populations is greater for STRPs than SNPs. These results should assist the design of multi-locus studies to detect population structure in nature.     

Six diagnostic single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trouts

Ten primer pairs were screened to develop single nucleotide polymorphism (SNP) TaqMan assays that will distinguish California golden trout and some rainbow trouts (Oncorhynchus mykiss sspp., O. m. aguabonita) from the Paiute and Lahontan cutthroat trouts (Oncorhynchus clarkii seleniris, O. c. henshawi). From these 10 primer pairs, one mitochondrial and five nuclear fixed SNP differences were discovered and developed into TaqMan assays. These six assays will be useful for characterizing and monitoring hybridization between these groups. Additional Oncorhynchus clarkii sspp. and Oncorhynchus mykiss sspp. were assayed to determine if these assays are useful in closely related species.     

Isolation and characterization of polymorphic microsatellite loci in the plural cooperatively breeding superb starling,Lamprotornis superbus

Superb starlings (Lamprotornis superbus) are separate nesting plural cooperative breeders endemic to East Africa that show high levels of reproductive conflict and low reproductive skew. I isolated 31 polymorphic microsatellite loci from the superb starling genome and characterized them in approximately 31 individuals. The numbers of alleles per locus ranged from two to 16 and observed heterozygosities ranged from 0.03 to 0.88. Many of these loci amplified in other passerine species including the European starling, Sturnus vulgaris, and a variety of other starlings and mockingbirds.     

Bayesian parentage analysis reliably controls the number of false assignments in natural populations

Parentage analysis in natural populations is a powerful tool for addressing a wide range of ecological and evolutionary questions. However, identifying parent–offspring pairs in samples collected from natural populations is often more challenging than simply resolving the Mendelian pattern of shared alleles. For example, large numbers of pairwise comparisons and limited numbers of genetic markers can contribute to incorrect assignments, whereby unrelated individuals are falsely identified as parent–offspring pairs. Determining which parentage methods are the least susceptible to making false assignments is an important challenge facing molecular ecologists. In a recent paper, Harrison et al. (2013a) address this challenge by comparing three commonly used parentage methods, including a Bayesian approach, in order to explore the effects of varied proportions of sampled parents on the accuracy of parentage assignments. Unfortunately, Harrison et al. made a simple error in using the Bayesian approach, which led them to incorrectly conclude that this method could not control the rate of false assignment. Here, I briefly outline the basic principles behind the Bayesian approach, identify the error made by Harrison et al., and provide detailed guidelines as to how the method should be correctly applied. Furthermore, using the exact data from Harrison et al., I show that the Bayesian approach actually provides greater control over the number of false assignments than either of the other tested methods. Lastly, I conclude with a brief introduction to solomon , a recently updated version of the Bayesian approach that can account for genotyping error, missing data and false matching.     

Bateman's principle is reversed in a cooperatively breeding bird

Bateman''s principle is not only used to explain sex differences in mating behaviour, but also to determine which sex has the greater opportunity for sexual selection. It predicts that the relationship between the number of mates and the number of offspring produced should be stronger for males than for females. Yet, it is unclear whether Bateman''s principle holds in cooperatively breeding systems where the strength of selection on traits used in intrasexual competition is high in both sexes. We tested Bateman''s principle in the cooperatively breeding superb starling (Lamprotornis superbus), finding that only females showed a significant, positive Bateman gradient. We also found that the opportunity for selection was on average higher in females, but that its strength and direction oscillated through time. These data are consistent with the hypothesis that sexual selection underlies the female trait elaboration observed in superb starlings and other cooperative breeders. Even though the Bateman gradient was steeper for females than for males, the year-to-year oscillation in the strength and direction of the opportunity for selection likely explains why cooperative breeders do not exhibit sexual role reversal. Thus, Bateman''s principle may not hold in cooperative breeders where both sexes appear to be under mutually strong sexual selection.     

Characterization of single nucleotide polymorphism markers for eelgrass (Zostera marina)

We characterized 37 single nucleotide polymorphism (SNP) makers for eelgrass Zostera marina. SNP markers were developed using existing EST (expressed sequence tag)-libraries to locate polymorphic loci and develop primers from the functional expressed genes that are deposited in The ZOSTERA database (V1.2.1). SNP loci were genotyped using a single-base-extension approach which facilitated high-throughput genotyping with minimal optimization time. These markers show a wide range of variability among 25 eelgrass populations and will be useful for population genetic studies including evaluation of population structure, historical demography, and phylogeography. Potential applications include haplotype inference of physically linked SNPs and identification of genes under selection for temperature and desiccation stress.     

Evaluation of 16 loci to examine the cross‐species utility of single nucleotide polymorphism arrays

Large collections of single nucleotide polymorphisms (SNPs) have recently been identified from a number of livestock genomes. This raises the possibility that SNP arrays might be useful for analysis in related species for which few genetic markers are currently available. To address the likely success of such an approach, the aim of this study was to examine the threshold number and position of flanking mutations which act to prevent genotype calls being produced. Sequence diversity was measured across 16 loci containing SNPs known either to work successfully between species or fail between species. In pairwise comparisons between domestic and wild sheep, sequence divergence surrounding working SNP assays was significantly lower than that surrounding non‐functional assays. In addition, the location of flanking mismatches tended to be closer to the target SNP in loci that failed to generate genotype calls across species. The magnitude of sequence divergence observed for both working and non‐functional assays was compared with the divergence separating domestic sheep from European Mouflon, African Barbary, goat and cattle. The results suggest that the utility of SNP arrays for analysis of shared polymorphism will be restricted to closely related pairs of species. Analysis across more divergent species will, however, be successful for other objectives, such as the identification of the ancestral state of SNPs.