Informativeness of genetic markers for pairwise relationship and relatedness inference

Parentage Analysis in Giant Grouper (Epinephelus lanceolatus) Using Microsatellite and SNP Markers from Genotyping-by-Sequencing Data

Pedigree information is necessary for the maintenance of diversity for wild and captive populations. Accurate pedigree is determined by molecular marker-based parentage analysis, which may be influenced by the polymorphism and number of markers, integrity of samples, relatedness of parents, or different analysis programs. Here, we described the first development of 208 single nucleotide polymorphisms (SNPs) and 11 microsatellites for giant grouper (Epinephelus lanceolatus) taking advantage of Genotyping-by-sequencing (GBS), and compared the power of SNPs and microsatellites for parentage and relatedness analysis, based on a mixed family composed of 4 candidate females, 4 candidate males and 289 offspring. CERVUS, PAPA and COLONY were used for mutually verification. We found that SNPs had a better potential for relatedness estimation, exclusion of non-parentage and individual identification than microsatellites, and > 98% accuracy of parentage assignment could be achieved by 100 polymorphic SNPs (MAF cut-off < 0.4) or 10 polymorphic microsatellites (mean H_o = 0.821, mean PIC = 0.651). This study provides a reference for the development of molecular markers for parentage analysis taking advantage of next-generation sequencing, and contributes to the molecular breeding, fishery management and population conservation.

Hare's affairs: Lessons learnt from a noninvasive genetic monitoring for tracking mountain hare individuals

Systematic monitoring of individuals and their abundance over time has become an important tool to provide information for conservation. For genetic monitoring studies, noninvasive sampling has emerged as a valuable approach, particularly so for elusive or rare animals. Here, we present the 5-year results of an ongoing noninvasive genetic monitoring of mountain hares (Lepus timidus) in a protected area in the Swiss Alps. We used nuclear microsatellites and a sex marker to identify individuals and assign species to noninvasively collected feces samples. Through including a marker for sex identification, we were able to assess sex ratio changes and sex-specific demographic parameters over time. Male abundance in the area showed high fluctuations and apparent survival for males was lower than for females. Generally, males and females showed only little temporary migration into and out of the study area. Additionally, using genotyped tissue samples from mountain hares, European hares (Lepus europaeus) and their hybrids, we were able to provide evidence for the first occurrence of a European hare in the study area at an elevation of 2,300 m a.s.l. in spring 2016. For future monitoring studies, we suggest to include complementary analysis methods to reliably infer species identities of the individuals analyzed and, thus, not only monitor mountain hare individual abundance, but also assess the potential threats given through competitive exclusion by and hybridization with the European hare.

Panmixia in a sea ice-associated marine mammal: evaluating genetic structure of the Pacific walrus (Odobenus rosmarus divergens) at multiple spatial scales

The kin structure of a species at relatively fine spatial scales impacts broad-scale patterns in genetic structure at the population level. However, kin structure rarely has been elucidated for migratory marine mammals. The Pacific walrus (Odobenus rosmarus divergens) exhibits migratory behavior linked to seasonal patterns in sea ice dynamics. Consequently, information on the spatial genetic structure of the subspecies, including kin structure, could aid wildlife managers in designing future studies to evaluate the impacts of sea ice loss on the subspecies. We sampled 8,303 individual walruses over a 5-year period and used 114 single-nucleotide polymorphisms to examine both broad-scale patterns in genetic structure and fine-scale patterns in relatedness. We did not detect any evidence of genetic structure at broad spatial scales, with low FST values (≤ 0.001) across all pairs of putative aggregations. To evaluate kin structure at fine spatial scales, we defined a walrus group as a cluster of resting individuals that were less than one walrus body length apart. We found weak evidence of kin structure at fine spatial scales, with 3.72% of groups exhibiting mean relatedness values greater than expected by chance, and a significantly higher overall observed mean value of relatedness within groups than expected by chance. Thus, the high spatiotemporal variation in the distribution of resources in the Pacific Arctic environment likely has favored a gregarious social system in Pacific walruses, with unrelated animals forming temporary associations.

Show me you care: female mate choice based on egg attendance rather than male or territorial traits

Female mate choice is often based on male traits, including signals or behaviors, and/or the quality of a male’s territory. In species with obligate paternal care, where care directly affects offspring survival, females may also base their mate choices on the quality of a sire’s care. Here, we quantified male reproductive success in a natural population of the glass frog Hyalinobatrachium cappellei, a species with male parental care, to determine the influence of territory quality, male traits, and paternal care behaviors on female mate choice. We found that attending males have a higher chance of gaining new clutches than nonattending males. Our results indicate that females do not select males based only on body condition, calling persistence, or territory traits. Instead, our findings support the hypothesis that females choose males based on care status. Indeed, males already attending a clutch were 70% more likely to obtain another clutch, and the time to acquire an additional clutch was significantly shorter. We also found that males adjust their parental care effort in response to genetic relatedness by caring only for their own offspring; however, remaining close to unrelated clutches serves as a strategy to attract females and increase chances of successful mating. Thus, males that establish territories that already contain clutches benefit from the signal eggs provide to females.

Genetic monogamy despite frequent extrapair copulations in "strictly monogamous" wild jackdaws

"Monogamy" refers to different components of pair exclusiveness: the social pair, sexual partners, and the genetic outcome of sexual encounters. Avian monogamy is usually defined socially or genetically, whereas quantifications of sexual behavior remain scarce. Jackdaws (Corvus monedula) are considered a rare example of strict monogamy in songbirds, with lifelong pair bonds and little genetic evidence for extrapair (EP) offspring. Yet jackdaw copulations, although accompanied by loud copulation calls, are rarely observed because they occur visually concealed inside nest cavities. Using full-day nest-box video surveillance and on-bird acoustic bio-logging, we directly observed jackdaw sexual behavior and compared it to the corresponding genetic outcome obtained via molecular parentage analysis. In the video-observed nests, we found genetic monogamy but frequently detected forced EP sexual behavior, accompanied by characteristic male copulation calls. We, thus, challenge the long-held notion of strict jackdaw monogamy at the sexual level. Our data suggest that male mate guarding and frequent intrapair copulations during the female fertile phase, as well as the forced nature of the copulations, could explain the absence of EP offspring. Because EP copulation behavior appeared to be costly for both sexes, we suggest that immediate fitness benefits are an unlikely explanation for its prevalence. Instead, sexual conflict and dominance effects could interact to shape the spatiotemporal pattern of EP sexual behavior in this species. Our results call for larger-scale investigations of jackdaw sexual behavior and parentage and highlight the importance of combining social, sexual, and genetic data sets for a more complete understanding of mating systems.

Genetic evidence supports polygamous mating system in a wild population of Prochilodus lineatus (Characiformes: Prochilodontidae), a Neotropical shoal spawner fish

ABSTRACT Behavioral observations made on fish have revealed remarkably diverse reproductive strategies, including polygamy by both sexes. Still, to date, most Neotropical species remain unstudied as to whether the observed reproductive behavior in natural populations correlates with their genetic mating systems. Here, we investigated the genetic mating system of a wild population of Prochilodus lineatus settled in the Middle Uruguay River basin. By using sibship reconstruction and parental inference methods based on microsatellites’ genotypes, we inferred 45 females and 47 males as potential parents of the 87 larvae analyzed. We found evidence supporting polygamous mating in both sexes: while a high percentage of males (44.7%) fertilized the eggs of one female, 55.3% of the inferred males fertilized eggs of up to four females. Likewise, while 44.5% of the inferred females had their eggs fertilized by one only male, 55.5% of females were fertilized by multiple males. The estimated proxy of the effective population size (Nb) was 126, exhibiting moderate to high levels of genetic diversity. The genetic evidence contributed in this study complements earlier behavioral observations of formation of spawning nuclei of aggregating breeders, which may be promoting a polygamous mating strategy in this long-distance migratory fish.

Pedigree reconstruction from poor quality genotype data

Marker genotype data could suffer from a high rate of errors such as false alleles and allelic dropouts (null alleles) in situations such as SNPs from low-coverage next-generation sequencing and microsatellites from noninvasive samples. Use of such data without accounting for mistyping properly could lead to inaccurate or incorrect inferences of family relationships such as parentage and sibship. This study shows that markers with a high error rate are still informative. Simply discarding them could cause a substantial loss of precious information, and is impractical in situations where virtually all markers (e.g. SNPs from low-coverage next-generation sequencing, microsatellites from noninvasive samples) suffer from a similarly high error rate. This study also shows that some previous error models are valid for markers of low error rates, but fail for markers of high error rates. It proposes an improved error model and demonstrates, using simulated and empirical data of a high error rate (say, >0.5), that it leads to more accurate sibship and parentage inferences than previous models. It suggests that, in reality, markers of high error rates should be used rather than discarded in pedigree reconstruction, so long as the error rates can be estimated and used properly in the analyses.

Genetic constraints of population expansion of the Carpathian lynx at the western edge of its native distribution range in Central Europe

Even though populations of many large carnivores are expanding throughout Europe, the Eurasian lynx population in the Western Carpathians seems unable to spread beyond the western boundaries of its current distributional range. Many factors, both extrinsic and intrinsic, can influence the potential for range expansion: landscape fragmentation, natal philopatry, low natural fecundity and high mortality, and low and sex-biased dispersal rates. In this study we used non-invasive genetic sampling to determine population size fluctuation, sub-structuring and social organisation of the peripheral lynx population at the Czech-Slovak border. Even though the population size has been relatively stable over the period studied (2010-2016), the individual inbreeding coefficients of residents at the end of the study were much higher than those of founders at the beginning of the study. While non-resident individuals (predominantly males) occurred regularly in the study population, only resident individuals with well-established home ranges participated in breeding and produced offspring. Almost half the offspring detected in the study (predominantly females) settled in or near the natal area. Subsequent incestuous mating resulted in production of inbred individuals, reduction of effective population size of the population, and sub-structuring of the population through formation of two distinct family lineages. Our study illustrates how social constraints, such as territoriality, breeding of residents and natal philopatry of females, lead to incestuous mating in small-sized populations, especially at the periphery of their distribution. This threat should be taken into account in planning of conservation and population recovery of species with similar social structure.

New perspectives on multilocus ancestry informativeness

We present an axiomatic approach for multilocus informativeness measures for determining the amount of information that a set of polymorphic genetic markers provides about individual ancestry. We then reveal several surprising properties of a decision-theoretic based measure that is consistent with the set of proposed criteria for multilocus informativeness. In particular, these properties highlight the interplay between information originating from population priors and the information extractable from the population genetic variants. This analysis then reveals a certain deficiency of mutual information based multilocus informativeness measures when such population priors are incorporated. Finally, we analyse and quantify the inevitable inherent decrease in informativeness due to learning from finite population samples.

Effects of Sample Size and Full Sibs on Genetic Diversity Characterization: A Case Study of Three Syntopic Iberian Pond-Breeding Amphibians

Accurate characterization of genetic diversity is essential for understanding population demography, predicting future trends and implementing efficient conservation policies. For that purpose, molecular markers are routinely developed for nonmodel species, but key questions regarding sampling design, such as calculation of minimum sample sizes or the effect of relatives in the sample, are often neglected. We used accumulation curves and sibship analyses to explore how these 2 factors affect marker performance in the characterization of genetic diversity. We illustrate this approach with the analysis of an empirical dataset including newly optimized microsatellite sets for 3 Iberian amphibian species: Hyla molleri, Epidalea calamita, and Pelophylax perezi. We studied 17-21 populations per species (total n = 547, 652, and 516 individuals, respectively), including a reference locality in which the effect of sample size was explored using larger samples (77-96 individuals). As expected, FIS and tests for Hardy-Weinberg equilibrium and linkage disequilibrium were affected by the presence of full sibs, and most initially inferred disequilibria were no longer statistically significant when full siblings were removed from the sample. We estimated that to obtain reliable estimates, the minimum sample size (potentially including full sibs) was close to 20 for expected heterozygosity, and between 50 and 80 for allelic richness. Our pilot study based on a reference population provided a rigorous assessment of marker properties and the effects of sample size and presence of full sibs in the sample. These examples illustrate the advantages of this approach to produce robust and reliable results for downstream analyses.

Female vervet monkeys fine-tune decisions on tolerance versus conflict in a communication network

Group living promotes opportunities for both cooperation and competition. Selection on the ability to cope with such opposing social opportunities has been proposed as a driving force in the evolution of large brains in primates and other social species. However, we still know little about the degree of complexity involved in such social strategies. Here, we report advanced social strategies in wild vervet monkeys. Building on recent experimental evidence that subordinate females trade grooming for tolerance from higher-ranking individuals during foraging activities, we show that the audience composition strongly affects this trade. First, tolerance was lower if the audience contained individuals that outranked the subordinate partner, independently of audience size and kinship relationships. Second, we found a significant interaction between previous grooming and relative rank of bystanders: dominant subjects valued recent grooming by subordinates while intermediate ranked subjects valued the option to aggress subordinate partners in the presence of a dominant audience. Aggressors were also more likely to emit coalition recruitment calls if the audience contained individuals that outranked the subordinate partner. In conclusion, vervet monkeys include both recent grooming and knowledge about third-party relationships to make complex decisions when trading grooming for tolerance, leading to a finely balanced trade-off between reciprocation and opportunities to reinforce rank relationships.

Genome-wide SNPs reveal the drivers of gene flow in an urban population of the Asian Tiger Mosquito, Aedes albopictus

Aedes albopictus is a highly invasive disease vector with an expanding worldwide distribution. Genetic assays using low to medium resolution markers have found little evidence of spatial genetic structure even at broad geographic scales, suggesting frequent passive movement along human transportation networks. Here we analysed genetic structure of Aedes albopictus collected from 12 sample sites in Guangzhou, China, using thousands of genome-wide single nucleotide polymorphisms (SNPs). We found evidence for passive gene flow, with distance from shipping terminals being the strongest predictor of genetic distance among mosquitoes. As further evidence of passive dispersal, we found multiple pairs of full-siblings distributed between two sample sites 3.7 km apart. After accounting for geographical variability, we also found evidence for isolation by distance, previously undetectable in Ae. albopictus. These findings demonstrate how large SNP datasets and spatially-explicit hypothesis testing can be used to decipher processes at finer geographic scales than formerly possible. Our approach can be used to help predict new invasion pathways of Ae. albopictus and to refine strategies for vector control that involve the transformation or suppression of mosquito populations.

Aedes albopictus, the Asian Tiger Mosquito, is a highly invasive disease vector with a growing global distribution. Designing strategies to prevent invasion and to control Ae. albopictus populations in invaded regions requires knowledge of how Ae. albopictus disperses. Studies comparing Ae. albopictus populations have found little evidence of genetic structure even between distant populations, suggesting that dispersal along human transportation networks is common. However, a more specific understanding of dispersal processes has been unavailable due to an absence of studies using high-resolution genetic markers. Here we present a study using high-resolution markers, which investigates genetic structure among 152 Ae. albopictus from Guangzhou, China. We found that human transportation networks, particularly shipping terminals, had an influence on genetic structure. We also found genetic distance was correlated with geographical distance, the first such observation in this species. This study demonstrates how high-resolution markers can be used to investigate ecological processes that may otherwise escape detection. We conclude that strategies for controlling Ae. albopictus will have to consider both passive reinvasion along human transportation networks and active reinvasion from neighbouring regions.

Microsatellite Mutation Rate in Atlantic Sturgeon (Acipenser oxyrinchus)

Understanding mutation rates can greatly extend the utility of population and conservation genetic analyses. Herein, we present an estimate of genome-wide microsatellite mutation rate in Atlantic sturgeon (Acipenser oxyrinchus) based on parent-offspring transmission patterns. We screened 307 individuals for parentage and mutation-rate analysis applying 43 variable markers. Out of 13228 allele transfers, 11 mutations were detected, producing a mutation rate of 8.3 × 10-4 per locus per generation (95% confidence interval: 1.48 × 10-3, 4.15 × 10-4). Single-step mutations predominated and there were trends toward mutations in loci with greater polymorphism and allele length. Two of the detected mutations were most probably cluster mutations, being identified in 12 and 28 sibs, respectively. Finally, we observed evidences of polyploidy based on the sporadic presence of 3 or 4 alleles per locus in the genotyped individuals, supporting previous reports of incomplete diploidization in Atlantic sturgeon.

Use of noninvasive 'bug-eggs' to enable comparative inferences on genetic mating system with and without parental information: A study in a cattle egret colony

Colonial waterbirds such as herons, egrets and spoonbills exhibit ecological characteristics that could have promoted the evolution of conspecific brood parasitism and extra-pair copulation. However, an adequate characterization of the genetic mating systems of this avian group has been hindered by the lack of samples of elusive candidate parents which precluded conducting conventional parentage allocation tests. Here, we investigate the genetic mating system of the invasive cattle egret using hematophagous insects contained in fake eggs to collect blood from incubating adults in a wild breeding colony. We tested a protocol with a previously unused Neotropical Triatominae, Panstrongylus megistus, obtained blood samples from males and females in 31 nests built on trees, drew blood from 89 nestlings at those nests, and genotyped all samples at 14 microsatellite loci, including six new species-specific loci. We comparatively addressed the performance of parentage allocation versus kinship classification of nestlings to infer the genetic mating system of cattle egrets. In line with previous behavioral observations, we found evidence in support of a non-monogamous genetic mating system, including extra-pair paternity (EPP) and conspecific brood parasitism (CBP). Parentage allocation tests detected a higher percentage of nests with alternative reproductive tactics (EPP: 61.7%; CBP: 64.5%) than the kinship classification method (EPP: 50.0%; CBP: 43.3%). Overall, these results indicate that rates of alternative reproductive tactics inferred in the absence of parental genetic information could be underestimated and should be interpreted with caution. This study highlights the importance of incorporating samples from candidate parents to adequately determine the genetic mating system of a species. We expand knowledge on the reproductive tactics of colonial waterbirds, contributing novel data on the genetic mating system of the cattle egret, valuable for the design of management strategies for this invasive bird.

Novel efficient genome-wide SNP panels for the conservation of the highly endangered Iberian lynx

BACKGROUND

The Iberian lynx (Lynx pardinus) has been acknowledged as the most endangered felid species in the world. An intense contraction and fragmentation during the twentieth century left less than 100 individuals split in two isolated and genetically eroded populations by 2002. Genetic monitoring and management so far have been based on 36 STRs, but their limited variability and the more complex situation of current populations demand more efficient molecular markers. The recent characterization of the Iberian lynx genome identified more than 1.6 million SNPs, of which 1536 were selected and genotyped in an extended Iberian lynx sample.

METHODS

We validated 1492 SNPs and analysed their heterozygosity, Hardy-Weinberg equilibrium, and linkage disequilibrium. We then selected a panel of 343 minimally linked autosomal SNPs from which we extracted subsets optimized for four different typical tasks in conservation applications: individual identification, parentage assignment, relatedness estimation, and admixture classification, and compared their power to currently used STR panels.

RESULTS

We ascribed 21 SNPs to chromosome X based on their segregation patterns, and identified one additional marker that showed significant differentiation between sexes. For all applications considered, panels of autosomal SNPs showed higher power than the currently used STR set with only a very modest increase in the number of markers.

CONCLUSIONS

These novel panels of highly informative genome-wide SNPs provide more powerful, efficient, and flexible tools for the genetic management and non-invasive monitoring of Iberian lynx populations. This example highlights an important outcome of whole-genome studies in genetically threatened species.

Accurate recapture identification for genetic mark-recapture studies with error-tolerant likelihood-based match calling and sample clustering

Error-tolerant likelihood-based match calling presents a promising technique to accurately identify recapture events in genetic mark-recapture studies by combining probabilities of latent genotypes and probabilities of observed genotypes, which may contain genotyping errors. Combined with clustering algorithms to group samples into sets of recaptures based upon pairwise match calls, these tools can be used to reconstruct accurate capture histories for mark-recapture modelling. Here, we assess the performance of a recently introduced error-tolerant likelihood-based match-calling model and sample clustering algorithm for genetic mark-recapture studies. We assessed both biallelic (i.e. single nucleotide polymorphisms; SNP) and multiallelic (i.e. microsatellite; MSAT) markers using a combination of simulation analyses and case study data on Pacific walrus (Odobenus rosmarus divergens) and fishers (Pekania pennanti). A novel two-stage clustering approach is demonstrated for genetic mark-recapture applications. First, repeat captures within a sampling occasion are identified. Subsequently, recaptures across sampling occasions are identified. The likelihood-based matching protocol performed well in simulation trials, demonstrating utility for use in a wide range of genetic mark-recapture studies. Moderately sized SNP (64+) and MSAT (10-15) panels produced accurate match calls for recaptures and accurate non-match calls for samples from closely related individuals in the face of low to moderate genotyping error. Furthermore, matching performance remained stable or increased as the number of genetic markers increased, genotyping error notwithstanding.

Postdispersal nepotism in male long-tailed macaques (Macaca fascicularis)

Cooperative behaviors are promoted by kin selection if the costs to the actor are smaller than the fitness benefits to the recipient, weighted by the coefficient of relatedness. In primates, cooperation occurs primarily among female dyads. Due to male dispersal before sexual maturity in many primate species, however, it is unknown whether there are sufficient opportunities for selective tolerance and occasional coalitionary support for kin selection to favor male nepotistic support. We studied the effect of the presence of male kin on correlates of male reproductive success (residence time, duration of high dominance rank) in non‐natal male long‐tailed macaques (Macaca fascicularis). We found that “related” (i.e., related at the half‐sibling level or higher) males in a group have a significantly higher probability to remain in the non‐natal group compared to males without relatives. Moreover, males stayed longer in a group when a relative was present at group entry or joined the same group within 3 months upon arrival. Males with co‐residing relatives also maintained a high rank for longer than those without. To our knowledge, this is the first demonstration of a potential nepotistic effect on residence and rank maintenance among non‐natal males in a social system without long‐term alliances.

The use and abuse of genetic marker-based estimates of relatedness and inbreeding

Genetic marker-based estimators remain a popular tool for measuring relatedness (r xy ) and inbreeding (F) coefficients at both the population and individual level. The performance of these estimators fluctuates with the number and variability of markers available, and the relatedness composition and demographic history of a population. Several methods are available to evaluate the reliability of the estimates of r xy and F, some of which are implemented in the program COANCESTRY. I used the simulation module in COANCESTRY since assess the performance of marker-based estimators of r xy and F in a species with very low genetic diversity, New Zealand's little spotted kiwi (Apteryx owenii). I also conducted a review of published papers that have used COANCESTRY as its release to assess whether and how the reliability of the estimates of r xy and F produced by genetic markers are being measured and reported in published studies. My simulation results show that even when the correlation between true (simulated) and estimated r xy or F is relatively high (Pearson's r = 0.66-0.72 and 0.81-0.85, respectively) the imprecision of the estimates renders them highly unreliable on an individual basis. The literature review demonstrates that the majority of studies do not report the reliability of marker-based estimates of r xy and F. There is currently no standard practice for selecting the best estimator for a given data set or reporting an estimator's performance. This could lead to experimental results being interpreted out of context and render the robustness of conclusions based on measures of r xy and F debatable.

How Well Do Molecular and Pedigree Relatedness Correspond, in Populations with Diverse Mating Systems, and Various Types and Quantities of Molecular and Demographic Data?

Kinship analyses are important pillars of ecological and conservation genetic studies with potentially far-reaching implications. There is a need for power analyses that address a range of possible relationships. Nevertheless, such analyses are rarely applied, and studies that use genetic-data-based-kinship inference often ignore the influence of intrinsic population characteristics. We investigated 11 questions regarding the correct classification rate of dyads to relatedness categories (relatedness category assignments; RCA) using an individual-based model with realistic life history parameters. We investigated the effects of the number of genetic markers; marker type (microsatellite, single nucleotide polymorphism SNP, or both); minor allele frequency; typing error; mating system; and the number of overlapping generations under different demographic conditions. We found that (i) an increasing number of genetic markers increased the correct classification rate of the RCA so that up to >80% first cousins can be correctly assigned; (ii) the minimum number of genetic markers required for assignments with 80 and 95% correct classifications differed between relatedness categories, mating systems, and the number of overlapping generations; (iii) the correct classification rate was improved by adding additional relatedness categories and age and mitochondrial DNA data; and (iv) a combination of microsatellite and single-nucleotide polymorphism data increased the correct classification rate if <800 SNP loci were available. This study shows how intrinsic population characteristics, such as mating system and the number of overlapping generations, life history traits, and genetic marker characteristics, can influence the correct classification rate of an RCA study. Therefore, species-specific power analyses are essential for empirical studies.

Nest suitability, fine-scale population structure and male-mediated dispersal of a solitary ground nesting bee in an urban landscape

Bees are the primary pollinators of flowering plants in almost all ecosystems. Worldwide declines in bee populations have raised awareness about the importance of their ecological role in maintaining ecosystem functioning. The naturally strong philopatric behavior that some bee species show can be detrimental to population viability through increased probability of inbreeding. Furthermore, bee populations found in human-altered landscapes, such as urban areas, can experience lower levels of gene flow and effective population sizes, increasing potential for inbreeding depression in wild bee populations. In this study, we investigated the fine-scale population structure of the solitary bee Colletes inaequalis in an urbanized landscape. First, we developed a predictive spatial model to detect suitable nesting habitat for this ground nesting bee and to inform our field search for nests. We genotyped 18 microsatellites in 548 female individuals collected from nest aggregations throughout the study area. Genetic relatedness estimates revealed that genetic similarity among individuals was slightly greater within nest aggregations than among randomly chosen individuals. However, genetic structure among nest aggregations was low (Nei's GST = 0.011). Reconstruction of parental genotypes revealed greater genetic relatedness among females than among males within nest aggregations, suggesting male-mediated dispersal as a potentially important mechanism of population connectivity and inbreeding avoidance. Size of nesting patch was positively correlated with effective population size, but not with other estimators of genetic diversity. We detected a positive trend between geographic distance and genetic differentiation between nest aggregations. Our landscape genetic models suggest that increased urbanization is likely associated with higher levels of inbreeding. Overall, these findings emphasize the importance of density and distribution of suitable nesting patches for enhancing bee population abundance and connectivity in human dominated habitats and highlights the critical contribution of landscape genetic studies for enhanced conservation and management of native pollinators.

Multiple populations of pantropical spotted dolphins in Hawaiian waters

Understanding gene flow and dispersal patterns is important for predicting effects of natural events and anthropogenic activities on animal populations. In Hawaii, most species of odontocetes are managed as single populations. Recent exceptions include false killer whales, spinner dolphins, and common bottlenose dolphins, for which studies have shown fidelity to individual islands or groups of islands. Our study focused on pantropical spotted dolphins. We analyzed mitochondrial control region and 11 microsatellite loci from 101 individuals from 4 areas: Hawaii, Maui/Lanai, Oahu, and Kauai/Niihau. We examined F ST, F' ST, R ST, Jost's D, and ΦST and used TESS to estimate number of populations and assignment probabilities. Our results support genetic differentiation among Hawaii, Maui/Lanai, and Oahu and suggest that pantropical spotted dolphins near Kauai/Niihau are likely transient and in low numbers. Between island regions, F ST for microsatellites ranged from 0.016 to 0.045 and for mtDNA, from 0.011 to 0.282. F ' ST, ranged from 0.098 to 0.262 for microsatellites and 0.019 to 0.415 for mtDNA. R ST and ΦST showed similar results to F ST for microsatellites and mtDNA respectively, and Jost's D fell between F ST and F ' ST. TESS supported 3 populations, and greatest mean assignment probability by island region ranged from 0.50 to 0.72. The private alleles method indicated migration rates among regions from 1.49 to 3.45, and effective population size of the island of Hawaii was estimated to be 220. There was no strong evidence to support sex-biased dispersal or group fidelity. Considering this study in the larger context of other odontocete population studies and studies of connectivity, we suggest genetic differentiation may be mediated by behavior adapted to differing habitat types and niches.

Fine-scale genetic assessment of sex-specific dispersal patterns in a multilevel primate society

Like humans, hamadryas baboons (Papio hamadryas) are unusual among primates in having a multilevel social system and stable pair bonds, and are thus a useful model for the evolution of human sociality. While the kinship structure and sex-biased dispersal patterns that underlie human social organization have been extensively elucidated, the impact of these factors on the social system of hamadryas baboons is currently unclear. Here we use genetic analysis of individuals to elucidate the patterns of male and female dispersal across multiple levels of society in a wild population of hamadryas baboons. We characterized 244 members of five hamadryas bands at Filoha, Ethiopia by genotyping one Y-linked and 23 autosomal microsatellite loci and sequencing part of the mitochondrial hypervariable control region I. We found both male and female dispersal to be limited at the level of the band, with more movement of females than males among bands. By integrating long-term behavioral data for Band 1, we also found evidence for male and female philopatry at the clan level. We speculate that male philopatry at the clan level and female dispersal across one-male units and clans may enable both kin-based cooperation among males and the maintenance of kin bonds among females after dispersal. This would mean that, as in humans, kin bonds within both sexes are a core feature of the hamadryas social system, thus contributing to our understanding of the evolution of social organization in humans.

A simulation module in the computer program COLONY for sibship and parentage analysis

A simulation module is built into the software package COLONY to simulate marker genotype data of individuals with a predefined parentage and sibship structure. The simulated data can then be used to compare the accuracy, robustness and computational efficiency of different methods for sibship and parentage reconstruction, to examine the impact of different parameter options in a software on its accuracy and computational efficiency and to assess the information sufficiency of a given set of markers for a sibship and parentage analysis. This computer note describes the method used for simulating genotype data with a pedigree and its possible applications. The method can quickly generate genotype data for a one- or two-generation pedigree of virtually any complexity with up to 30k offspring, at up to 30k codominant or dominant loci with an arbitrary degree of linkage and a user-defined mistyping rate. The data can be fed directly into the COLONY program for analysis by three sibship and parentage reconstruction methods and can also be imported into other programs such as Excel and R. With slight modification, the data can be analysed by other relationship analysis software.

Conspecific brood parasitism in the white-faced ibis Plegadis chihi (Aves: Pelecaniformes) revealed by microsatellites' based kinship-reconstruction

The white-faced ibis Plegadis chihi Vieillot, 1817 (Pelecaniformes: Threskiornithidae) is a socially monogamous colonially breeding bird in which behavioral and ecological observations suggest the occurrence of conspecific brood parasitism (CBP). We inferred aspects of the genetic mating system of P. chihi in nature, using a genetic approach in the absence of parental information. We used five heterologous microsatellite loci and a multiple-step methodological approach to infer kinship patterns among 104 pairs of nestlings sampled inside 80 nests in a breeding colony from southern Brazil. The estimated effective population size was 69 white-faced ibises (95% CI: 50-98), enough to ensure long-term population survival. Kinship patterns were identified for 38% of the analyzed pairs: 60% of the diagnosed pairs were identified as full-siblings, 2.5% as half-siblings and 37.5% as unrelated individuals. CBP could explain the presence of unrelated nestlings within broods, in agreement with available non-genetic evidence. The presence of half-siblings within broods could indicate extra-pair paternity. Results suggest that a non-strictly monogamous genetic mating system may be present in the white-faced ibis. This study is the first molecular approach to better characterize the reproductive behavior of P. chihi in the wild. Our findings set the stage for further research to investigate the possible causes and consequences of alternative reproductive strategies in this species.

A comparison of four methods for detecting weak genetic structure from marker data

Genetic structure is ubiquitous in wild populations and is the result of the processes of natural selection, genetic drift, mutation, and gene flow. Genetic drift and divergent selection promotes the generation of genetic structure, while gene flow homogenizes the subpopulations. The ability to detect genetic structure from marker data diminishes rapidly with a decreasing level of differentiation among subpopulations. Weak genetic structure may be unimportant over evolutionary time scales but could have important implications in ecology and conservation biology. In this paper we examine methods for detecting and quantifying weak genetic structures using simulated data. We simulated populations consisting of two putative subpopulations evolving for up to 50 generations with varying degrees of gene flow (migration), and varying amounts of information (allelic diversity). There are a number of techniques available to detect and quantify genetic structure but here we concentrate on four methods: F_ST, population assignment, relatedness, and sibship assignment. Under the simple mating system simulated here, the four methods produce qualitatively similar results. However, the assignment method performed relatively poorly when genetic structure was weak and we therefore caution against using this method when the analytical aim is to detect fine-scale patterns. Further work should examine situations with different mating systems, for example where a few individuals dominate reproductive output of the population. This study will help workers to design their experiments (e.g., sample sizes of markers and individuals), and to decide which methods are likely to be most appropriate for their particular data.

Monitoring the effective population size of a brown bear (Ursus arctos) population using new single-sample approaches

The effective population size (N(e) ) could be the ideal parameter for monitoring populations of conservation concern as it conveniently summarizes both the evolutionary potential of the population and its sensitivity to genetic stochasticity. However, tracing its change through time is difficult in natural populations. We applied four new methods for estimating N(e) from a single sample of genotypes to trace temporal change in N(e) for bears in the Northern Dinaric Mountains. We genotyped 510 bears using 20 microsatellite loci and determined their age. The samples were organized into cohorts with regard to the year when the animals were born and yearly samples with age categories for every year when they were alive. We used the Estimator by Parentage Assignment (EPA) to directly estimate both N(e) and generation interval for each yearly sample. For cohorts, we estimated the effective number of breeders (N(b) ) using linkage disequilibrium, sibship assignment and approximate Bayesian computation methods and extrapolated these estimates to N(e) using the generation interval. The N(e) estimate by EPA is 276 (183-350 95% CI), meeting the inbreeding-avoidance criterion of N(e) > 50 but short of the long-term minimum viable population goal of N(e) > 500. The results obtained by the other methods are highly consistent with this result, and all indicate a rapid increase in N(e) probably in the late 1990s and early 2000s. The new single-sample approaches to the estimation of N(e) provide efficient means for including N(e) in monitoring frameworks and will be of great importance for future management and conservation.

Population connectivity: dam migration mitigations and contemporary site fidelity in arctic char

BACKGROUND

Animal feeding and spawning migrations may be limited by physical barriers and behavioral interactions. Dam constructions (e.g. hydropower) commonly include gateways for fish migrations to sustain ecological connectivity. Relative genetic impacts of fish passage devices versus natural processes (e.g. hybrid inferiority) are, however, rarely studied. We examined genetic (i.e. microsatellite) population connectivity of highly migrating lake-dwelling Arctic char (Salvelinus alpinus), introduced 20 generations ago, across and within two subalpine lakes separated by a dam with a subterranean tunnel and spill gates after 7 generations. Due to water flow regime, the time window for fish migration is highly restricted.

RESULTS

Char populations, with similar genetic structuring and diversity observed across and within lakes, were admixed across the dam with fishways during feeding. For spawning, however, statistically significant, but very low population differentiation (θ; 0.002 - 0.013) was found in nine out of ten reproductive site comparisons, reflecting interactions between extensive migration (mean first generation (F0) = 10.8%) and initial site fidelity. Simulations indicated that genetic drift among relatively small effective populations (mean N(e) = 62) may have caused the observed contemporary differentiation. Novel Bayesian analyses indicated mean contributions of 71% F0 population hybrids in spawning populations, of which 76% had maternal or paternal native origin.

CONCLUSIONS

Ecological connectivity between lakes separated by a dam has been retained through construction of fishways for feeding migration. Considerable survival and homing to ancestral spawning sites in hybrid progeny was documented. Population differentiation despite preceding admixture is likely caused by contemporary reduced reproductive fitness of population hybrids. The study documents the beginning stages of population divergence among spatial aggregations with recent common ancestry.

Effects of genotyping errors on parentage exclusion analysis

Genetic markers are widely used to determine the parentage of individuals in studies of mating systems, reproductive success, dispersals, quantitative genetic parameters and in the management of conservation populations. These markers are, however, imperfect for parentage analyses because of the presence of genotyping errors and undetectable alleles, which may cause incompatible genotypes (mismatches) between parents and offspring and thus result in false exclusions of true parentage. Highly polymorphic markers widely used in parentage analyses, such as microsatellites, are especially prone to genotyping errors. In this investigation, I derived the probabilities of excluding a random (related) individual from parentage and the probabilities of Mendelian-inconsistent errors (mismatches) and Mendelian-consistent errors (which do not cause mismatches) in parent-offspring dyads, when a marker having null alleles, allelic dropouts and false alleles is used in a parentage analysis. These probabilities are useful in evaluating the impact of various types of genotyping errors on the information content of a set of markers in and thus the power of a parentage analysis, in determining the threshold number of genetic mismatches that is appropriate for a parentage exclusion analysis and in estimating the rates of genotyping errors and frequencies of null alleles from observed mismatches between known parent-offspring dyads. These applications are demonstrated by numerical examples using both hypothetical and empirical data sets and discussed in the context of practical parentage exclusion analyses.

Detecting populations in the 'ambiguous' zone: kinship-based estimation of population structure at low genetic divergence

Identifying population structure is one of the most common and important objectives of spatial analyses using population genetic data. Population structure is detected either by rejecting the null hypothesis of a homogenous distribution of genetic variation, or by estimating low migration rates. Issues arise with most current population genetic inference methods when the genetic divergence is low among putative populations. Low levels of genetic divergence may be as a result of either high ongoing migration or historic high migration but no current, ongoing migration. We direct attention to recent developments in the use of the tempo-spatial distribution of closely related individuals to detect population structure or estimate current migration rates. These 'kinship-based' approaches complement more traditional population-based genetic inference methods by providing a means to detect population structure and estimate current migration rates when genetic divergence is low. However, for kinship-based methods to become widely adopted, formal estimation procedures applicable to a range of species life histories are needed.

On the use of large marker panels to estimate inbreeding and relatedness: empirical and simulation studies of a pedigreed zebra finch population typed at 771 SNPs

In recent years there has been a dramatic increase in the availability of high density genetic marker data for both model and non-model organisms. A potential application of these data is to infer relatedness in the absence of a complete pedigree. Using a marker panel of 771 SNPs genotyped in three generations of an extensive zebra finch pedigree, correlations between pedigree relatedness and seven marker-based estimates of relatedness were examined, as was the relationship between heterozygosity and inbreeding. Although marker-based and pedigree relatedness were highly correlated, the variance in estimated relatedness was high. Further, the correlation between heterozygosity and inbreeding was weak, even though mean inbreeding coefficient is typical of that seen in wild vertebrate pedigrees; the weak relationship was in part due to the small variance in inbreeding in the pedigree. Our data suggest that using marker information to reconstruct the pedigree, and then calculating relatedness from the pedigree, is likely to give more accurate relatedness estimates than using marker-based estimators directly.

Estimation of bumblebee queen dispersal distances using sibship reconstruction method

Dispersal ability is a key determinant of the propensity of an organism to cope with habitat fragmentation and climate change. Here we quantify queen dispersal in two common bumblebee species in an arable landscape. Dispersal was measured by taking DNA samples from workers in the spring and summer, and from queens in the following spring, at 14 sites across a landscape. The queens captured in the spring must be full sisters of workers that were foraging in the previous year. A range of sibship reconstruction methods were compared using simulated data sets including or no genotyping errors. The program Colony gave the most accurate reconstruction and was used for our analysis of queen dispersal. Comparison of queen dispersion with worker foraging distances was used to take into account an expected low level of false identification of sister pairs which might otherwise lead to overestimates of dispersal. Our data show that Bombus pascuorum and B. lapidarius queens can disperse by at least 3 and 5 km, respectively. These estimates are consistent with inferences drawn from studies of population structuring in common and rare bumblebee species, and suggest that regular gene flow over several kilometres due to queen dispersal are likely to be sufficient to maintain genetic cohesion of ubiquitous species over large spatial scales whereas rare bumblebee species appear unable to regularly disperse over distances greater than 10 km. Our results have clear implications for conservation strategies for this important pollinator group, particularly when attempting to conserve fragmented populations.