DNA-based methods for pedigree reconstruction and kinship analysis in natural populations

Structured Incorporation of Prior Information in Relationship Identification Problems

The objective of this paper is to show how various sources of information can be modelled and integrated to address relationship identification problems. Applications come from areas as diverse as evolution and conservation research, genealogical research in human, plant and animal populations, and forensic problems including paternity cases, identification following disasters, family reunions and immigration issues. We propose assigning a prior probability distribution to the sample space of pedigrees, calculating the likelihood based on DNA data using available software and posterior probabilities using Bayes' Theorem. Our emphasis here is on the modelling of this prior information in a formal and consistent manner. We introduce the distinction between local and global prior information, whereby local information usually applies to particular components of the pedigree and global prior information refers to more general features. When it is difficult to decide on a prior distribution, robustness to various choices should be studied. When suitable prior information is not available, a flat prior can be used which will then correspond to a strict likelihood approach. In practice, prior information is often considered for these problems, but in a generally ad hoc manner. This paper offers a consistent alternative. We emphasise that many practical problems can be addressed using freely available software.

A graphical approach to relatedness inference

The estimation of relatedness structure in natural populations using molecular marker data has become an important tool in population biology, resulting in a variety of estimation procedures for specific sampling scenarios. In this article a general approach is proposed, in which the detailed relationship structure, typically a pedigree graph or partition, is considered to be the object of inference. This makes available tools used in complex model selection theory which have demonstrated effectiveness. An important advantage of this approach is that it permits a fully Bayesian approach to the problem, providing a principled and accessible way to measure statistical error. The approach is demonstrated by applying the minimum description length principle. This technique is used in model selection to provide a rational way of comparing models of varying complexity. We show how the resulting score may be interpreted and applied as a Bayesian posterior density.

Life history correlates of inbreeding depression in mandrills (Mandrillus sphinx)

Inbreeding depression reflects the negative consequences of increased homozygosity at genes that affect fitness. We investigate inbreeding depression in a semi-free-ranging colony of mandrills (Mandrillus sphinx), using high-quality pedigree data, comprising five maternal generations and 20 years of morphological and demographic data. We examine the relationship between inbreeding coefficients and four fitness correlates: two growth parameters (mass and height for age) and longevity in both sexes, and age at first conception in females. Inbreeding was correlated with both growth parameters, but only in females, with inbred females being smaller than noninbred females. Inbreeding was also correlated significantly with age at first conception, with inbred females giving birth earlier in life than noninbred females. We suggest that sex-biased maternal investment may explain this sex-differential response to inbreeding, although the lack of a significant association between inbreeding and growth in males may also be due to the provisioned nature of the colony. The surprising relationship between age at first conception and inbreeding may be related to smaller adult size in inbred females, or to their being less able to escape from male sexual coercion.

Dispersal, philopatry and intergroup relatedness: fine-scale genetic structure in the white-breasted thrasher, Ramphocinclus brachyurus

Dispersal is a fundamental process influencing evolution, social behaviour, and the long-term persistence of populations. We use both observational and genetic data to investigate dispersal, kin-clustering and intergroup relatedness in the white-breasted thrasher, Ramphocinclus brachyurus, a cooperatively breeding bird that is globally endangered. Mark-resighting data suggested sex-biased dispersal, with females dispersing over greater distances while males remained philopatric. Accordingly, spatial autocorrelation analysis showed highly significant fine-scale genetic structure among males, but not among females. This fine-scale genetic structuring of the male population resulted in very high levels of relatedness between dominant males at neighbouring nests, similar to that seen within cooperative groups in many species where kin selection is cited as a cause of cooperation. By implication, between-group as well as within-group cooperation may be important, potentially creating a feedback loop in which short-distance dispersal by males leads to the formation of male kin clusters that in turn facilitate nepotistic interactions and favour further local recruitment. The strength of spatial autocorrelation, as measured by the autocorrelation coefficient, r, was approximately two to three times greater than that reported in previous studies of animals. Relatively short dispersal distances by both males and females may have a negative impact on the white-breasted thrasher's ability to colonize new areas, and may influence the long-term persistence of isolated populations. This should be taken into account when designating protected areas or selecting sites for habitat restoration.

Unravelling first-generation pedigrees in wild endangered salmon populations using molecular genetic markers

There is increasing interest in the use of molecular genetic data to infer genealogical relationships among individuals in the absence of parental information. Such analyses can provide insight into mating systems and estimations of heritability in the wild. In addition, accurate pedigree reconstruction among the founders of endangered populations being reared in captivity would be invaluable. Many breeding programs for endangered species attempt to minimize loss of genetic variation and inbreeding through strategies designed to minimize global co-ancestry, but they assume a lack of relatedness among the founders. Yet populations that are the target of such programs are generally in serious demographic decline, and many of the available founders may be closely related. Here we demonstrate determination of full and half-sib relationships among the wild founders of a captive breeding program involving two endangered Atlantic salmon populations using two different approaches and associated software, pedigree and colony. A large portion of the juveniles collected in these two rivers appear to be derived from surprisingly few females mating with a large number of males, probably small precocious parr. Another group of potential founders, obtained from a local hatchery, clearly originated from a small number of full-sib crosses. These results allowed us to prioritize individuals on the basis of conservation value, and are expected to help minimize loss of genetic variation through time. In addition, insight is provided into the number of contributing parents and the mating systems that produced this last generation of endangered wild Atlantic salmon.

Towards unbiased parentage assignment: combining genetic, behavioural and spatial data in a Bayesian framework

Inferring the parentage of a sample of individuals is often a prerequisite for many types of analysis in molecular ecology, evolutionary biology and quantitative genetics. In all but a few cases, the method of parentage assignment is divorced from the methods used to estimate the parameters of primary interest, such as mate choice or heritability. Here we present a Bayesian approach that simultaneously estimates the parentage of a sample of individuals and a wide range of population-level parameters in which we are interested. We show that joint estimation of parentage and population-level parameters increases the power of parentage assignment, reduces bias in parameter estimation, and accurately evaluates uncertainty in both. We illustrate the method by analysing a number of simulated test data sets, and through a re-analysis of parentage in the Seychelles warbler, Acrocephalus sechellensis. A combination of behavioural, spatial and genetic data are used in the analyses and, importantly, the method does not require strong prior information about the relationship between nongenetic data and parentage.

Estimating relatedness and relationships using microsatellite loci with null alleles

Relatedness is often estimated from microsatellite genotypes that include null alleles. When null alleles are present, observed genotypes represent one of several possible true genotypes. If null alleles are detected, but analyses do not adjust for their presence (ie, observed genotypes are treated as true genotypes), then estimates of relatedness and relationship can be incorrect. The number of loci available in many wildlife studies is limited, and loci with null alleles are commonly a large proportion of data that cannot be discarded without substantial loss of power. To resolve this problem, we present a new approach for estimating relatedness and relationships from data sets that include null alleles. Once it is recognized that the probability of the observed genotypes is dependent on the probabilities of a limited number of possible true genotypes, the required adjustments are straightforward. The concept can be applied to any existing estimators of relatedness and relationships. We review established maximum likelihood estimators and apply the correction in that setting. In an application of the corrected method to data from striped hyenas, we demonstrate that correcting for the presence of null alleles affect results substantially. Finally, we use simulated data to confirm that this method works better than two common approaches, namely ignoring the presence of null alleles or discarding affected loci.

Performance of marker-based relatedness estimators in natural populations of outbred vertebrates

Knowledge of relatedness between pairs of individuals plays an important role in many research areas including evolutionary biology, quantitative genetics, and conservation. Pairwise relatedness estimation methods based on genetic data from highly variable molecular markers are now used extensively as a substitute for pedigrees. Although the sampling variance of the estimators has been intensively studied for the most common simple genetic relationships, such as unrelated, half- and full-sib, or parent-offspring, little attention has been paid to the average performance of the estimators, by which we mean the performance across all pairs of individuals in a sample. Here we apply two measures to quantify the average performance: first, misclassification rates between pairs of genetic relationships and, second, the proportion of variance explained in the pairwise relatedness estimates by the true population relatedness composition (i.e., the frequencies of different relationships in the population). Using simulated data derived from exceptionally good quality marker and pedigree data from five long-term projects of natural populations, we demonstrate that the average performance depends mainly on the population relatedness composition and may be improved by the marker data quality only within the limits of the population relatedness composition. Our five examples of vertebrate breeding systems suggest that due to the remarkably low variance in relatedness across the population, marker-based estimates may often have low power to address research questions of interest.

Genealogy of wine grape cultivars: "Pinot" is related to "Syrah"

Since the domestication of wild grapes ca 6000 years ago, numerous cultivars have been generated by spontaneous or deliberate crosses, and up to 10 000 are still in existence today. Just as in human paternity analysis, DNA typing can reveal unexpected parentage of grape cultivars. In this study, we have analysed 89 grape cultivars with 60 microsatellite markers in order to accurately calculate the identity-by-descent (IBD) and relatedness (r) coefficients among six putatively related cultivars from France ("Pinot", "Syrah" and "Dureza") and northern Italy ("Teroldego", "Lagrein" and "Marzemino"). Using a recently developed likelihood-based approach to analyse kinship in grapes, we provide the first evidence of a genetic link between grapes across the Alps: "Dureza" and "Teroldego" turn out to be full-siblings (FS). For the first time in grapevine genetics we were able to detect FS without knowing one of the parents and identify unexpected second-degree relatives. We reconstructed the most likely pedigree that revealed a third-degree relationship between the worldwide-cultivated "Pinot" from Burgundy and "Syrah" from the Rhone Valley. Our finding was totally unsuspected by classical ampelography and it challenges the commonly assumed independent origins of these grape cultivars. Our results and this new approach in grape genetics will (a) help grape breeders to avoid choosing closely related varieties for new crosses, (b) provide pedigrees of cultivars in order to detect inheritance of disease-resistance genes and (c) open the way for future discoveries of first- and second-degree relationships between grape cultivars in order to better understand viticultural migrations.

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Seasonal changes in pollinator activity influence pollen dispersal and seed production of the alpine shrub Rhododendron aureum (Ericaceae)

In alpine ecosystems, microscale variation in snowmelt timing often causes different flowering phenology of the same plant species and seasonal changes in pollinator activity. We compared the variations in insect visitation, pollen dispersal, mating patterns, and sexual reproduction of Rhododendron aureum early and late in the flowering season using five microsatellites. Insects visiting the flowers were rare early in the flowering season (mid-June), when major pollinators were bumblebee queens and flies. In contrast, frequent visitations by bumblebee workers were observed late in the season (late July). Two-generation analysis of pollen pool structure demonstrated that quality of pollen-mediated gene flow was more diverse late in the season in parallel with the high pollinator activity. The effective number of pollen donors per fruit (N(ep)) increased late in the season (N(ep) = 2.2-2.7 early, 3.4-4.4 late). However, both the outcrossing rate (t(m)) and seed-set ratio per fruit were smaller late in the season (t(m) = 0.89 and 0.71, seed-set ratio = 0.52 and 0.18, early and late in the season, respectively). In addition, biparental inbreeding occurred only late in the season. We conclude that R. aureum shows contrasting patterns of pollen movement and seed production between early and late season: in early season, seed production can be high but genetically less diverse and, during late season, be reduced, possibly due to higher inbreeding and inbreeding depression, but have greater genetic diversity. Thus, more pollinator activity does not always mean more pollen movement.

Kin distribution of amphibian larvae in the wild

According to kin selection theory, the location of an individual with respect to its relatives can have important ramifications for its fitness. Perhaps more than any other vertebrate group, anuran amphibian larvae have been the subject of many experiments on this topic. Some anuran species have been shown in the laboratory to recognize and associate with their siblings and half-siblings. However, due to the difficulty of identifying sibships, no kinship studies with anuran larvae have been conducted in the wild. Here, we use microsatellite analysis to show that wood frog (Rana sylvatica) tadpoles were nonrandomly distributed in two ponds with respect to their relatives. In one pond, the tadpoles were significantly clumped with their siblings or half-siblings as expected from other published laboratory studies on this species. However, in another pond, the tadpoles were significantly nonrandomly dispersed from their siblings or half-siblings. This is the first example of kin repulsion of nonreproductive animals in the wild and the first time a species has been shown to display both aggregation and repulsion under different circumstances. These results suggest that kin distribution is context dependent and demonstrate the importance of testing kin selection hypotheses under natural conditions.

Estimating relatedness between individuals in general populations with a focus on their use in conservation programs

Relatedness estimators are widely used in genetic studies, but effects of population structure on performance of estimators, criteria to evaluate estimators, and benefits of using such estimators in conservation programs have to date received little attention. In this article we present new estimators, based on the relationship between coancestry and molecular similarity between individuals, and compare them with existing estimators using Monte Carlo simulation of populations, either panmictic or structured. Estimators were evaluated using statistical criteria and a diversity criterion that minimized relatedness. Results show that ranking of estimators depends on the population structure. An existing estimator based on two-gene and four-gene coefficients of identity performs best in panmictic populations, whereas a new estimator based on coancestry performs best in structured populations. The number of marker alleles and loci did not affect ranking of estimators. Statistical criteria were insufficient to evaluate estimators for their use in conservation programs. The regression coefficient of pedigree relatedness on estimated relatedness (beta2) was substantially lower than unity for all estimators, causing overestimation of the diversity conserved. A simple correction to achieve beta2 = 1 improves both existing and new estimators. Using relatedness estimates with correction considerably increased diversity in structured populations, but did not do so or even decreased diversity in panmictic populations.

Interspecies mating in sympatric species of Syngnathus pipefish

Hybridization is thought to be an important source of novel genetic variation, and interspecific hybridization may increase the adaptive potential of wild populations. While hybridization has not been previously reported in syngnathid fishes (seahorses and pipefish), the sympatric occurrence of closely related species at high densities increases the probability of interspecies mating in this group. Southern California is home to five species of Syngnathus pipefish, and these species frequently co-occur in near-shore eelgrass beds along the California coast. Recent work has identified exceptionally high levels of genetic diversity in southern populations of Syngnathus leptorhynchus, a widespread species which ranges from Mexico to Alaska. Microsatellite genotyping and mitochondrial sequence data are used here to study the population genetics of S. leptorhynchus and S. auliscus at a site in San Diego Bay where they are found to co-occur at high densities. While no adult hybrids were detected in the study population, analysis of male broods indicates that interspecies mating is occurring between the two species. The lack of premating isolating mechanisms between these two relatives suggests that hybridization may be common in sympatric species of Syngnathus.

The dilemma of female mate selection in the brown bear, a species with sexually selected infanticide

Because of differential investment in gametes between sexes, females tend to be the more selective sex. Based on this concept, we investigate mate selection in a large carnivore: the brown bear (Ursus arctos). We hypothesize that, in this species with sexually selected infanticide (SSI), females may be faced with a dilemma: either select a high-quality partner based on phenotypic criteria, as suggested by theories of mate choice, or rather mate with future potentially infanticidal males as a counter-strategy to SSI. We evaluated which male characteristics were important in paternity assignment. Among males available in the vicinity of the females, the largest, most heterozygous and less inbred and also the geographically closest males were more often the fathers of the female's next litter. We suggest that female brown bears may select the closest males as a counter-strategy to infanticide and exercise a post-copulatory cryptic choice, based on physical attributes, such as a large body size, reflecting male genetic quality. However, male-male competition either in the form of fighting before copulation or during the post-copulatory phase, in the form of sperm competition, cannot entirely be ruled out.

Informativeness of genetic markers for pairwise relationship and relatedness inference

Measuring the information content of markers in relationship/relatedness inferences is important in selecting highly informative markers to attain a given statistical power with the minimal genotyping effort. Using information-theoretic principles, I introduce the informativeness for relationship (I(R)) and the informativeness for relatedness (I(r)) to measure the amount of information provided by markers in inferring pairwise relationships (R) and relatedness (r), respectively. I also propose a fast and accurate algorithm to calculate the power (PW(R)) of a set of markers in differentiating two candidate relationships, and the reciprocal of the mean squared deviations of relatedness estimates (RMSD) to measure the amount of information of markers actually used by an estimator in estimating relatedness. All of the four measurements (I(R), I(r), PW(R), RMSD) apply to dominant and codominant markers, haploid and diploid individuals, and take into account of mutations and typing errors in data. The statistical properties of the four measurements and their relationships are investigated analytically and are examined by applying these methods to simulated and empirical data.

Phenotypic approaches for understanding patterns of intracemetery biological variation

This paper reviews studies of phenotypic inheritance and microevolutionary processes in archaeological populations using data on cranial and dental phenotypic variation, often referred to as paleogenetics or biodistance analysis. The estimation of biological distances between populations, or among individuals within populations, is one component of bioarchaeological research on past populations. In this overview, five approaches that focus on morphological variation within cemeteries are summarized: kinship and cemetery structure analysis, postmarital residence analysis, sample aggregate phenotypic variability, temporal microchronology, and age-structured phenotypic variation. Previous research, theoretical justifications, and methods are outlined for each topic. Case studies are presented that illustrate these theoretical and methodological bases, as well as demonstrate the kinds of inferences possible using these approaches. Kinship and cemetery structure analysis seeks to identify the members of family groups within larger cemeteries or determine whether cemeteries were kin-structured. Analysis of sex-specific phenotypic variation allows estimation of postmarital residence practices, which is important for understanding other aspects of prehistoric social organization. Analysis of aggregate phenotypic variability can be used to infer site formation processes or cemetery catchment area. The study of temporal microchronologies can be used to evaluate provisional archaeological chronologies or study microevolutionary processes such as adaptive selection or changing patterns of gene flow. Finally, age-structured phenotypic variation can be reflective of selection processes within populations or it can be used as a measure of morbidity, growth arrest, and early mortality within past populations. Use of phenotypic data as a genotypic proxy is theoretically sound, even at small scales of analysis.

Genetic Structure Inside a Declining Red Oak Community in Old-Growth Forest

Problems with oak regeneration have been documented in the last 50 years at numerous sites in the Midwestern United States. We applied nuclear microsatellites to examine the demographic and fine-scale spatial genetic structure of red oaks in two old-growth stands in Indiana. Oaks in one stand have declined in numbers over the past several decades whereas oaks in the other, smaller stand have increased. Large amounts of genetic variation were maintained within stands, and there was slight but significant differentiation among stands. There was significant but weak isolation by distance genetic structure within the large stand, likely reflecting family structure. No significant differences exist in allele frequencies or in levels of genetic diversity between cohorts that remain well represented within each stand, even between medium-sized adults and those antedating European settlement of the area. However, a virtual absence of smaller size classes in the forest interior of the large stand represents the early stages of a genetic bottleneck in what had been the core habitat of this stand. Whether future generations of this old-growth stand will retain the present genetic character depends on the oaks regenerating at the forest margins, absent any major changes in disturbance regimes. Similar demographic and genetic dynamics are likely occurring in a large number of remnant oak forests across the Midwest.

Paternity and kinship patterns in polyandrous moustached tamarins (Saguinus mystax)

We studied patterns of genetic relatedness and paternity in moustached tamarins, small Neotropical primates living in groups of 1-4 adult males and 1-4 adult females. Generally only one female per group breeds, mating with more than one male. Twin birth are the norm. In order to examine the genetic consequences of this mating pattern, DNA was extracted from fecal samples collected from two principal and six neighboring groups. DNA was characterized at twelve microsatellite loci (average: seven alleles/locus). We addressed the following questions: Do all adult males have mating access to the reproductive female of the group? How is paternity distributed across males in a group? Can polyandrous mating lead to multiple paternity? Are nonparental animals more closely related to the breeders than to the population mean? And, are mating partners unrelated? Breeding females mated with all nonrelated males. In at least one group the father of the older offspring did not sire the youngest infant although he was still resident in the group. We also found evidence for multiple paternity in a supposed twin pair. Yet, within each group the majority (67-100%) of infants had the same father, suggesting reproductive skew. Relatedness within groups was generally high (average R = 0.31), although both nonrelated males and females occurred, i.e., immigrations of both sexes are possible. Mating partners were never found to be related, hence inbreeding seems to be uncommon. The results suggest that while the social mating system is polyandry, paternity is often monopolized by a single male per group.

Migratory charr schools exhibit population and kin associations beyond juvenile stages

Few studies have critically investigated the genetic composition of wild fish schools. Yet, such investigations may have profound implications for the understanding of social organization and population differentiation in both fundamental and applied research. Using 20 microsatellite loci, we investigated the composition of 53 schools (total n = 211) of adult and subadult migratory brook charr (Salvelinus fontinalis) sampled from the known feeding areas of two populations inhabiting Mistassini Lake (Québec, Canada). We specifically tested whether (i) school members originated from the same population, (ii) individuals from the same population within schools were kin (half- or full-siblings), and (iii) kin schooling relationships differed between sexes. Randomization tests revealed a tendency for most schools to be population specific, although some schools were population mixtures. Significantly more kin were found within schools than expected at random for both populations (approximately 21-34% of the total number of school members). This result, combined with the observed size range of individuals, indicated that stable associations between kin may occur beyond juvenile stages for up to 4 years. Nevertheless, a high proportion of school members were non-kin (approximately 66-79%). No differences were detected between sexes in the propensity to school with kin. We discuss the hypothesis that the stable kin groups, rather than arising from kin selection, may instead be a by-product of familiarity based on individual selection for the maintenance of local adaptations related to migration (natal and feeding area philopatry). Our results are noteworthy because they suggest that there is some degree of permanence in the composition of wild fish schools. Additionally, they support the hypothesis that schools can be hierarchically structured (from population members down to family groups) and are thus nonrandom genetic entities.

Multilocus estimation of pairwise relatedness with dominant markers

Estimators for pairwise relatedness designed for dominant markers are derived, based on a genetic model that accounts for the full structure of pairwise relatedness between two individuals at a diploid locus with dominance. They jointly estimate 'relatedness' and 'fraternity', in which case the estimators are inherently multilocus, as at least two loci of differing gene frequency are required. Extensions to cases of zero fraternity and isolation by distance (inbreeding) are also examined. Properties of estimators are examined by simulation and compared to the estimator of Hardy. The most statistical power for pairwise relatedness occurs when roughly half of individuals are the recessive phenotype. Estimation procedures are implemented in the computer program mark.

Modified SIMPSON O(n3) algorithm for the full sibship reconstruction problem

MOTIVATION

The problem of reconstructing full sibling groups from DNA marker data remains a significant challenge for computational biology. A recently published heuristic algorithm based on Mendelian exclusion rules and the Simpson index was successfully applied to the full sibship reconstruction (FSR) problem. However, the so-called SIMPSON algorithm has an unknown complexity measure, questioning its applicability range.

RESULTS

We present a modified version of the SIMPSON (MS) algorithm that behaves as O(n(3)) and achieves the same or better accuracy when compared with the original algorithm. Performance of the MS algorithm was tested on a variety of simulated diploid population samples to verify its complexity measure and the significant improvement in efficiency (e.g. 100 times faster than SIMPSON in some cases). It has been shown that, in theory, the SIMPSON algorithm runs in non-polynomial time, significantly limiting its usefulness. It has been also verified via simulation experiments that SIMPSON could run in O(n(a)), where a > 3.

AVAILABILITY

Computer code written in Java is available upon request from the first author.

CONTACT

Dmitry.Konovalov@jcu.edu.au.

Measuring polyandry in wild populations: a case study using promiscuous crickets

Mating rate has important implications for patterns of sexual selection and sexual conflict and hence for issues such as speciation and the maintenance of genetic diversity. Knowledge of natural mating rates can provide insights into the factors driving female mating behaviour. We investigated the level of polyandry in a Spanish population of the field cricket Gryllus bimaculatus using microsatellite markers. Two approaches were employed: (i) genotyping the offspring of wild-caught gravid females to determine the number of males siring the brood and (ii) genotyping sperm stored in the spermathecae of females mated in the wild to estimate the number of mating partners. We compared existing methods for inferring the minimum and probable number of fathers and described a novel probabilistic technique estimating the number of mates by genotyping stored sperm. Using the most conservative allele-counting method, 71% of females produced offspring sired by at least two males (a minimum mean of 2.4 fathers per clutch), and all females had mated to at least two males with minimum mean estimates of 2.7-5.1 mates per female. Our study reveals high levels of polyandry in the wild and suggests that females mate with more males than sire their offspring.

The estimation of genetic relationships using molecular markers and their efficiency in estimating heritability in natural populations

Molecular marker data collected from natural populations allows information on genetic relationships to be established without referencing an exact pedigree. Numerous methods have been developed to exploit the marker data. These fall into two main categories: method of moment estimators and likelihood estimators. Method of moment estimators are essentially unbiased, but utilise weighting schemes that are only optimal if the analysed pair is unrelated. Thus, they differ in their efficiency at estimating parameters for different relationship categories. Likelihood estimators show smaller mean squared errors but are much more biased. Both types of estimator have been used in variance component analysis to estimate heritability. All marker-based heritability estimators require that adequate levels of the true relationship be present in the population of interest and that adequate amounts of informative marker data are available. I review the different approaches to relationship estimation, with particular attention to optimizing the use of this relationship information in subsequent variance component estimation.

Molecular ecology of parasites: elucidating ecological and microevolutionary processes

We review studies that have used molecular markers to address ecological and microevolutionary processes in parasites. Our goal is to highlight areas of research that may be of particular interest in relation to the parasitic lifestyle, and to draw attention to areas that require additional study. Topics include species identification, phylogeography, host specificity and speciation, population genetic structure, modes of reproduction and transmission patterns, and searching for loci under selection.

Mountain gorilla tug-of-war: silverbacks have limited control over reproduction in multimale groups

To determine who fathers the offspring in wild mountain gorilla groups containing more than one adult male silverback, we genotyped nearly one-fourth (n = 92) of the mountain gorillas (Gorilla beringei beringei) living in the Virunga Volcanoes region of Africa. Paternity analysis of 48 offspring born into four groups between 1985 and 1999 revealed that, although all infants were sired by within-group males, the socially dominant silverback did not always monopolize reproduction within his group. Instead, the second-ranking male sired an average of 15% of group offspring. This result, in combination with previous findings that second-ranking males fare best by not leaving the group but by staying and waiting to assume dominance even if no reproduction is possible while waiting, is not consistent with expectations from a reproductive skew model in which the silverback concedes controllable reproduction to the second-ranking male. Instead, the data suggest a "tug-of-war" scenario in which neither the dominant nor the second-ranking male has full control over his relative reproductive share. The two top-ranked males were typically unrelated and this, in combination with the mixed paternity of group offspring, means that multimale gorilla groups do not approximate family groups. Instead, as long-term assemblages of related and unrelated individuals, gorilla groups are similar to chimpanzee groups and so offer interesting possibilities for kin-biased interactions among individuals.

Backward simulation of ancestors of sampled individuals

If the population is large and the sampling mechanism is random, the coalescent is commonly used to model the haplotypes in the sample. Ordered genotypes can then be formed by random matching of the derived haplotypes. However, this approach is not realistic when (1) there is departure from random mating (e.g., dominant individuals in breeding populations or monogamy in humans), or (2) the population is small and/or the individuals in the sample are ascertained by applying some particular non-random sampling scheme, as is usually the case when considering the statistical modeling and analysis of pedigree data. For such situations, we present here a data generation method where an ancestral graph with non-overlapping generations is first generated backwards in time, using ideas from coalescent theory. Alleles are randomly assigned to the founders, and subsequently the gene flow over the entire genome is simulated forwards in time by dropping alleles down the graph according to recombination model without interference. The parameters controlling the mating behavior of generated individuals in the graph (degree of monogamy) can be tuned in order to match a particular demographic situation, without restriction to simple random mating. The performance of the approach is illustrated with a simulation example. The software (written in C-language) is freely available for research purposes at http://www.rni.helsinki.fi/~dag/.

Identifying nineteenth century genealogical links from genotypes

We have developed a likelihood method to identify moderately distant genealogical relationships from genomewide scan data. The aim is to compare the genotypes of many pairs of people and identify those pairs most likely to be related to one another. We have tested the algorithm using the genotypes of 170 Tasmanians with multiple sclerosis recruited into a haplotype association study. It is estimated from genealogical records that approximately 65% of Tasmania's current population of 470,000 are direct descendants of the 13,000 female founders living in this island state of Australia in the mid-nineteenth century. All cases and four to five relatives of each case have been genotyped with microsatellite markers at a genomewide average density of 4 cM. Previous genealogical research has identified 51 pairwise relationships linking 56 of the 170 cases. Testing the likelihood calculation on these known relative pairs, we have good power to identify relationships up to degree eight (e.g. third cousins once removed). Applying the algorithm to all other pairs of cases, we have identified a further 61 putative relative pairs, with an estimated false discovery rate of 10%. The power to identify genealogical links should increase when the new, denser sets of SNP markers are used. Except in populations where there is a searchable electronic database containing virtually all genealogical links in the past six generations, the algorithm should be a useful aid for genealogists working on gene-mapping projects, both linkage studies and association studies.

How to use molecular marker data to measure evolutionary parameters in wild populations

Estimating the genetic basis of phenotypic traits and the selection pressures acting on them are central to our understanding of the evolution and conservation of wild populations. However, obtaining such evolutionary-related parameters is not an easy task as it requires accurate information on both relatedness among individuals and their breeding success. Polymorphic molecular markers are very useful in estimating relatedness between individuals and parentage analyses are now extensively used in most taxa. The next step in the application of molecular data to wild populations is to use them to derive estimates of evolutionary-related parameters for quantitative traits, such as quantitative genetic parameters (e.g. heritability, genetic correlations) and measures of selection (e.g. selection gradients). Despite their great appeal and potential, the optimal use of molecular tools is still debated and it remains unclear how they should best be used to obtain reliable estimates of evolutionary parameters in the wild. Here, we review the methods available for estimating quantitative genetic and selection parameters and discuss their merits and shortcomings, to provide a tool that summarizes the potential uses of molecular data to obtain such parameters in wild populations.

Multiplexed microsatellite markers for the genetic analysis of Eucalyptus leucoxylon (Myrtaceae) and their utility for ecological and breeding studies in other eucalyptus species

Eucalyptus leucoxylon is a widespread woodland tree species found in southeastern Australia that has suffered from, and continues to be, threatened by the impacts of habitat clearance and degradation. Populations now consist predominantly of scattered individuals, and their conservation status is of increasing concern. We report the development and characterization of a set of eight highly polymorphic microsatellite loci for E. leucoxylon. The loci can be amplified in three PCR multiplexes and electrophoresed in a single lane, allowing rapid throughput of large numbers of samples. A total of 111 alleles were detected in 68 individuals with an average of 12.3 alleles per locus, a mean expected heterozygosity of 0.83, and a mean observed heterozygosity of 0.72. The combined probabilities of identity and probabilities of paternity exclusion allow an extremely precise level of individual identification, indicating that these microsatellite markers will be ideal for population genetic and parentage-type studies in E. leucoxylon. The markers also exhibited an average of 76% conservation within the subgenus Symphyomyrtus, to which E. leucoxylon belongs, and 53% conservation across other subgenera of Eucalyptus, demonstrating the potential of these markers in ecological and breeding studies in a wide range of Eucalyptus species.

SNPs in ecological and conservation studies: a test in the Scandinavian wolf population

Single nucleotide polymorphisms (SNPs) have the potential to become the genetic marker of choice in studies of the ecology and conservation of natural populations because of their capacity to access variability across the genome. In this study, we provide one of the first demonstrations of SNP discovery in a wild population in order to address typical issues of importance in ecology and conservation in the recolonized Scandinavian and neighbouring Finnish wolf Canis lupus populations. Using end sequence from BAC (bacterial artificial chromosome) clones specific for dogs, we designed assays for 24 SNP loci, 20 sites of which had previously been shown to be polymorphic in domestic dogs and four sites were newly identified as polymorphic in wolves. Of the 24 assayed loci, 22 SNPs were found to be variable within the Scandinavian population and, importantly, these were able to distinguish individual wolves from one another (unbiased probability of identity of 4.33 x 10(-8)), providing equivalent results to that derived from 12 variable microsatellites genotyped in the same population. An assignment test shows differentiation between the Scandinavian and neighbouring Finnish wolf populations, although not all known immigrants are accurately identified. An exploration of the misclassification rates in the identification of relationships shows that neither 22 SNP nor 20 microsatellite loci are able to discriminate across single order relationships. Despite the remaining obstacle of SNP discovery in nonmodel organisms, the use of SNPs in ecological and conservation studies is encouraged by the advent of large scale screening methods. Furthermore, the ability to amplify extremely small fragments makes SNPs of particular use for population monitoring, where faecal and other noninvasive samples are routinely used.

Kin selection and cooperative courtship in wild turkeys

In the few species of birds in which males form display partnerships to attract females, one male secures most or all of the copulations. This leads to the question of why subordinate males help in the absence of observable reproductive benefits. Hamilton's concept of kin selection, whereby individuals can benefit indirectly by helping a relative, was a crucial breakthrough for understanding apparently altruistic systems. However in the only direct test of kin selection in coordinated display partnerships, partners were unrelated, discounting kin selection as an explanation for the evolution of cooperation. Here I show, using genetic measures of relatedness and reproductive success, that kin selection can explain the evolution of cooperative courtship in wild turkeys. Subordinate (helper) males do not themselves reproduce, but their indirect fitness as calculated by Hamilton's rule more than offsets the cost of helping. This result confirms a textbook example of kin selection that until now has been controversial and also extends recent findings of male relatedness on avian leks by quantifying the kin-selected benefits gained by non-reproducing males.

A simple method for the calculation of microsatellite genotype distances irrespective of ploidy level

Abstract Microsatellites are powerful molecular markers, used commonly to estimate intraspecific genetic distances. With the exception of band sharing similarity index, available distance measures were developed specifically for diploid organisms and are unsuited for comparisons of polyploids. Here, we present a simple method for calculation of microsatellite genotype distances, which takes into account mutation processes and permits comparison of individuals with different ploidy levels. This method should provide a valuable tool for intraspecific analyses of polyploid organisms, which are widespread among plants and some animal taxa. An illustration is given using data from the planarian flatworm Schmidtea polychroa (Platyhelminthes).

Use of microsatellites for parentage and kinship analyses in animals

Microsatellite markers are quickly becoming the molecular marker of choice for studies of parentage and kinship in animals. In this chapter, we review methods and give protocols for screening potential microsatellite markers, as well as protocols for genotyping individuals with useful markers once they have been identified. In addition, we explain how microsatellites can be used to assess parentage and kinship, give basic analytical methods, and briefly review more sophisticated approaches that can be used to circumvent many of the problems that arise in any real empirical study.

The expansion of conservation genetics

The 'crisis discipline' of conservation biology has voraciously incorporated many technologies to speed up and increase the accuracy of conservation decision-making. Genetic approaches to characterizing endangered species or areas that contain endangered species are prime examples of this. Technical advances in areas such as high-throughput sequencing, microsatellite analysis and non-invasive DNA sampling have led to a much-expanded role for genetics in conservation. Such expansion will allow for more precise conservation decisions to be made and, more importantly, will allow conservation genetics to contribute to area- and landscape-based decision-making processes.

Ex situ population management in the absence of pedigree information

For captive breeding to play a significant role in conservation, ex situ populations must be scientifically managed to meet objective goals for retaining representative genetic variation. Imperfect genealogical information requires fundamental assumptions to be made that may bias downstream measures of genetic importance, upon which management decisions are based. The impacts of such assumptions are most pronounced within breeding programmes characterized by a high proportion of individuals of unknown ancestry, as exemplified by the large captive population of the St Vincent parrot (Amazona guildingii). The degree to which microsatellite-based estimates of relatedness may improve upon the assumptions of conventional pedigree-based management was investigated using genotypic data collected at eight microsatellite loci and two marker-based relatedness estimators. The measure, rxyLR, was found to explain the highest amount of variation in true relatedness. Integration of pairwise estimates of founder relatedness with studbook data transformed current understanding of the relatedness structure of the A. guildingii population from two subgroups characterized by a high and low degree of relatedness, respectively, to a situation where all 72 individuals are prioritized for breeding according to their estimated mean kinships. Furthermore, the discovery of opposing, directional bias exhibited by rxyLR and rxyQG in assigning dyads to a given relationship category suggests that an approach that utilizes a combination of pairwise relatedness estimators may provide the most genetic information for balancing the dual considerations of maximizing gene diversity and minimizing inbreeding in developing breeding recommendations.