Ancestral Stories of Ghanaian Bimoba Reflect Millennia-Old Genetic Lineages

Oral history and oral genealogies are mechanisms of collective memory and a main cultural heritage of many populations without a writing system. In the effort to analytically address the correspondence between genetic data and historical genealogies, anthropologists hypothesised that genealogies evolve through time, ultimately containing three parts: literal – where the most recent ancestry is truthfully represented; intended – where ancestry is inferred and reflects political relations among groups; and mythical – that does not represent current social reality. While numerous studies discuss oral genealogies, to our knowledge no genetic studies have been able to investigate to what extent genetic relatedness corresponds to the literal and intended parts of oral genealogies. We report on the correspondence between genetic data and oral genealogies among Bimoba males in a single village in North-Eastern Ghana. We compared the pairwise mismatch distribution of Y chromosome short tandem repeat (Y-STR) haplotypes among all lineages present in this village to the self-reported (oral) relatedness. We found that Bimoba are able to correctly identify unrelated individuals in 92% of the cases. In contrast, they are able to correctly identify related individuals only in 38% of the cases, which can be explained by three processes: (1) the compression of genealogies, leading to increasing inaccuracy with increasing genealogical distance, (2) inclusions into the lineage from intended relations such as clan co-option or adoptions, and (3) false paternities, which in this study were found to have a minor effect on the correspondence between genetic data and oral genealogies. In addition, we observed that 70% of unrelated pairs have from six to eight Y-STR differences, a diversification peak which we attribute to an ancient West African expansion dating around 9454 years ago. We conclude that, despite all caveats, oral genealogies are reflecting ancient lineages more accurately than previously thought.     

Use of microsatellite loci to classify individuals by relatedness

This study investigates the use of microsatellite loci for estimating relatedness between individuals in wild, outbred, vertebrate populations. We measured allele frequencies at 20 unlinked, dinucleotide-repeat microsatellite loci in a population of wild mice ( Mus musculus ), and used these observed frequencies to generate the expected distributions of pairwise relatedness among full sib, half sib, and unrelated pairs of individuals, as would be estimated from the microsatellite data. In this population one should be able to discriminate between unrelated and full-sib dyads with at least 97% accuracy, and to discriminate half-sib pairs from unrelated pairs or from full-sib pairs with better than 80% accuracy. If one uses the criterion that parent-offspring pairs must share at least one allele per locus, then only 15% of full-sib pairs, 2% of half-sib pairs, and 0% of unrelated pairs in this population would qualify as potential parent-offspring pairs. We verified that the simulation results (which assume a random mating population in Hardy-Weinberg and linkage equilibrium) accurately predict results one would obtain from this population in real life by scoring laboratory-bred full- and half-sib families whose parents were wild-caught mice from the study population. We also investigated the effects of using different numbers of loci, or loci of different average heterozygosities ( He ), on misclassification frequencies. Both variables have strong effects on misclassification rate. For example, it requires almost twice as many loci of He = 0.62 to achieve the same accuracy as a given number of loci of He = 0.75. Finally, we tested the ability of UPGMA clustering to identify family groups in our population. Clustering of allele matching scores among the offspring of four sets of independent maternal half sibships (four females, each mated to two different males) perfectly recovered the true family relationships.     

Inference of relationships in population data using identity-by-descent and identity-by-state

It is an assumption of large, population-based datasets that samples are annotated accurately whether they correspond to known relationships or unrelated individuals. These annotations are key for a broad range of genetics applications. While many methods are available to assess relatedness that involve estimates of identity-by-descent (IBD) and/or identity-by-state (IBS) allele-sharing proportions, we developed a novel approach that estimates IBD0, 1, and 2 based on observed IBS within windows. When combined with genome-wide IBS information, it provides an intuitive and practical graphical approach with the capacity to analyze datasets with thousands of samples without prior information about relatedness between individuals or haplotypes. We applied the method to a commonly used Human Variation Panel consisting of 400 nominally unrelated individuals. Surprisingly, we identified identical, parent-child, and full-sibling relationships and reconstructed pedigrees. In two instances non-sibling pairs of individuals in these pedigrees had unexpected IBD2 levels, as well as multiple regions of homozygosity, implying inbreeding. This combined method allowed us to distinguish related individuals from those having atypical heterozygosity rates and determine which individuals were outliers with respect to their designated population. Additionally, it becomes increasingly difficult to identify distant relatedness using genome-wide IBS methods alone. However, our IBD method further identified distant relatedness between individuals within populations, supported by the presence of megabase-scale regions lacking IBS0 across individual chromosomes. We benchmarked our approach against the hidden Markov model of a leading software package (PLINK), showing improved calling of distantly related individuals, and we validated it using a known pedigree from a clinical study. The application of this approach could improve genome-wide association, linkage, heterozygosity, and other population genomics studies that rely on SNP genotype data.     

Estimating heritabilities and genetic correlations with marker-based methods: an experimental test in Mimulus guttatus

The calculation of heritabilities and genetic correlations, which are necessary for predicting evolutionary responses, requires knowledge about the relatedness between individuals. This information is often not directly available, especially not for natural populations, but can be inferred by using molecular markers such as allozymes. Several methods based on inferred relatedness from marker data have been developed to estimate heritabilities and genetic correlations in natural populations. Most methods use maximum-likelihood procedures to assign pairs or groups of individuals to predefined discrete relatedness classes (e.g., half sibs and unrelated individuals). The Ritland method, on the other hand, uses method of moments estimators to estimate pairwise relatedness among individuals as continuous values. We tested both the Ritland method and a maximum-likelihood method by applying them to a greenhouse population consisting of seed families of the herb Mimulus guttatus and comparing the results to the ones from a frequently used standard method based on half-sib families. Estimates of genetic correlations were far from accurate, especially when we used the Ritland method. However, this study shows that even with a few variable allozyme loci, it is possible to get qualitatively good indications about the presence of heritable genetic variation from marker-based methods, even though both methods underestimated it.     

A new method to estimate relatedness from molecular markers

Four major problems can affect the efficiency of methods developed to estimate relatedness between individuals from information of molecular markers: (i) some of them are dependent on the knowledge of the true allelic frequencies in the base population; (ii) they assume that all loci are unlinked and in Hardy-Weinberg and linkage equilibrium; (iii) pairwise methods can lead to incongruous assignations because they take into account only two individuals at a time; (iv) most are usually constructed for particular structured populations (only consider a few relationship classes, e.g. full-sibs vs. unrelated). We have developed a new approach to estimate relatedness that is free from the above limitations. The method uses a 'blind search algorithm' (actually simulated annealing) to find the genealogy that yield a co-ancestry matrix with the highest correlation with the molecular co-ancestry matrix calculated using the markers. Thus (i and ii) it makes no direct assumptions about allelic frequencies or Hardy-Weinberg and linkage equilibrium; (iii) it always provide congruent relationships, as it considers all individuals at a time; (iv) degrees of relatedness can be as complex as desired just increasing the 'depth' (i.e. number of generations) of the proposed genealogies. Computer simulations have shown that the accuracy and robustness against genotyping errors of this new approach is comparable to that of other proposed methods in those particular situations they were developed for, but it is more flexible and can cope with more complex situations.     

基于亲权鉴定的千岛湖社鼠家群遗传结构与亲缘关系特征研究

于2009年7月至2010年11月,对浙江千岛湖两个岛屿上的社鼠(Niviventer confucianus)种群进行标志重捕,并采用8个高多态性的微卫星位点,对两个岛屿的社鼠种群进行家群分析和亲权鉴定,探讨了社鼠家群的亲缘关系特征。结果显示,8个微卫星位点能可靠地对两个岛屿社鼠种群进行亲权鉴定,A岛已确定亲缘关系的71只社鼠分为12个家群,家群中的个体数最多达到19个,B岛已确定的49只社鼠个体共分为11个家群,家群中的个体数最多达到14个。家群内部成员之间的亲缘关系表现为配对繁殖的个体对间亲缘系数最小,揭示了社鼠倾向于选择亲缘关系较远的异性作为配偶。家群中雄性后代个体之间与雌性个体之间的亲缘关系相比,两岛表现情况相反,该结果暗示两岛屿上社鼠扩散行为可能有所不同。通过计算与同一雄性(或同一雌性)交配的个体间的亲缘系数,发现两个岛屿上的社鼠在与不同异性交配时也存在选择性,即避免选择亲缘关系较近的异性作为混交的对象。     

Inference of Unexpected Genetic Relatedness among Individuals in HapMap Phase III

The International Haplotype Map Project (HapMap) has provided an essential database for studies of human population genetics and genome-wide association. Phases I and II of the HapMap project generated genotype data across ∼3 million SNP loci in 270 individuals representing four populations. Phase III provides dense genotype data on ∼1.5 million SNPs, generated by Illumina and Affymetrix platforms in a larger set of individuals. Release 3 of phase III of the HapMap contains 1397 individuals from 11 populations, including 250 of the original 270 phase I and phase II individuals and 1147 additional individuals. Although some known relationships among the phase III individuals have been described in the data release, the genotype data that are currently available provide an opportunity to empirically ascertain previously unknown relationships. We performed a systematic analysis of genetic relatedness and were able not only to confirm the reported relationships, but also to detect numerous additional, previously unidentified pairs of close relatives in the HapMap sample. The inferred relative pairs make it possible to propose standardized subsets of unrelated individuals for use in future studies in which relatedness needs to be clearly defined.     

Inbreeding avoidance in cunningham's skinks (Egernia cunninghami) in natural and fragmented habitat

Habitat fragmentation/alteration has been proposed as a distinct process threatening the viability of populations of many organisms. One expression of its impact may be the disruption of core population processes such as inbreeding avoidance. Using the experimental design outlined in our companion paper, we report on the impact of habitat alteration (deforestation) on inbreeding in the rock-dwelling Australian lizard Egernia cunninghami. Ten microsatellite loci were used to calculate relatedness coefficients of potential and actual breeding pairs, and to examine mate-choice and heterozygosity. Despite significantly less dispersal and higher within-group relatedness between potential mates in deforested than in natural habitats, this did not result in significantly more inbred matings. Average relatedness amongst breeding pairs was low, with no significant difference between natural and fragmented populations in relatedness between breeding pairs, or individual heterozygosity. Active avoidance of close kin as mates was indicated by the substantially and significantly lower relatedness in actual breeding pairs than potential ones. These facts, and heterozygote excesses in all groups of immature lizards from both habitats, show that E. cunninghami maintained outbreeding in the face of increased accumulation of relatives.     

Empirical comparison of Simple Sequence Repeats and single nucleotide polymorphisms in assessment of maize diversity and relatedness

While Simple Sequence Repeats (SSRs) are extremely useful genetic markers, recent advances in technology have produced a shift toward use of single nucleotide polymorphisms (SNPs). The different mutational properties of these two classes of markers result in differences in heterozygosities and allele frequencies that may have implications for their use in assessing relatedness and evaluation of genetic diversity. We compared analyses based on 89 SSRs (primarily dinucleotide repeats) to analyses based on 847 SNPs in individuals from the same 259 inbred maize lines, which had been chosen to represent the diversity available among current and historic lines used in breeding. The SSRs performed better at clustering germplasm into populations than did a set of 847 SNPs or 554 SNP haplotypes, and SSRs provided more resolution in measuring genetic distance based on allele-sharing. Except for closely related pairs of individuals, measures of distance based on SSRs were only weakly correlated with measures of distance based on SNPs. Our results suggest that 1) large numbers of SNP loci will be required to replace highly polymorphic SSRs in studies of diversity and relatedness and 2) relatedness among highly-diverged maize lines is difficult to measure accurately regardless of the marker system.     

Population genetics and relatedness in a critically endangered island raptor, Ridgway’s Hawk Buteo ridgwayi

Many island avian populations are of conservation interest because they have a higher risk of extinction than mainland populations. Susceptibility of island birds to extinction is primarily related to human induced change through habitat loss, persecution, and introduction of exotic species, in combination with genetic factors. We used microsatellite profiles from 11 loci to assess genetic diversity and relatedness in the critically endangered hawk Buteo ridgwayi endemic to the island of Hispaniola in the Caribbean. Using samples collected between 2005 and 2009, our results revealed a relatively high level of heterozygosity, evidence of a recent genetic bottleneck, and the occurrence of inbreeding within the population. Pair relatedness analysis found 4 of 7 sampled breeding pairs to be related similar to that of first cousin or greater. Pedigree estimates indicated that up to 18 % of potential pairings would be between individuals with relatedness values similar to that of half-sibling. We discuss our findings in the context of conservation genetic management suggesting both carefully managed translocations and the initiation of a captive population as a safeguard of the remaining genetic diversity.     

Estimating kinship in admixed populations

Genome-wide association studies (GWASs) are commonly used for the mapping of genetic loci that influence complex traits. A problem that is often encountered in both population-based and family-based GWASs is that of identifying cryptic relatedness and population stratification because it is well known that failure to appropriately account for both pedigree and population structure can lead to spurious association. A number of methods have been proposed for identifying relatives in samples from homogeneous populations. A strong assumption of population homogeneity, however, is often untenable, and many GWASs include samples from structured populations. Here, we consider the problem of estimating relatedness in structured populations with admixed ancestry. We propose a method, REAP (relatedness estimation in admixed populations), for robust estimation of identity by descent (IBD)-sharing probabilities and kinship coefficients in admixed populations. REAP appropriately accounts for population structure and ancestry-related assortative mating by using individual-specific allele frequencies at SNPs that are calculated on the basis of ancestry derived from whole-genome analysis. In simulation studies with related individuals and admixture from highly divergent populations, we demonstrate that REAP gives accurate IBD-sharing probabilities and kinship coefficients. We apply REAP to the Mexican Americans in Los Angeles, California (MXL) population sample of release 3 of phase III of the International Haplotype Map Project; in this sample, we identify third- and fourth-degree relatives who have not previously been reported. We also apply REAP to the African American and Hispanic samples from the Women's Health Initiative SNP Health Association Resource (WHI-SHARe) study, in which hundreds of pairs of cryptically related individuals have been identified.     

The cryptic genetic structure of the North American captive gorilla population

Western lowland gorillas (Gorilla gorilla gorilla) were imported from across their geographical range to North American zoos from the late 1800s through 1974. The majority of these gorillas were imported with little or no information regarding their original provenance and no information on their genetic relatedness. Here, we analyze 32 microsatellite loci in 144 individuals using a Bayesian clustering method to delineate clusters of individuals among a sample of founders of the captive North American zoo gorilla collection. We infer that the majority of North American zoo founders sampled are distributed into two distinct clusters, and that some individuals are of admixed ancestry. This new information regarding the existence of ancestral genetic population structure in the North American zoo population lays the groundwork for enhanced efforts to conserve the evolutionary units of the western lowland gorilla gene pool. Our data also show that the genetic diversity estimates in the founder population were comparable to those in wild gorilla populations (Mondika and Cross River), and that pairwise relatedness among the founders is no different from that expected for a random mating population. However, the relatively high level of relatedness (R = 0.54) we discovered in a pair of known breeding pairs reveals the need for incorporating genetic relatedness estimates in the captive management of western lowland gorillas.     

A comparison of approaches to estimate the inbreeding coefficient and pairwise relatedness using genomic and pedigree data in a sheep population

Genome-wide SNP data provide a powerful tool to estimate pairwise relatedness among individuals and individual inbreeding coefficient. The aim of this study was to compare methods for estimating the two parameters in a Finnsheep population based on genome-wide SNPs and genealogies, separately. This study included ninety-nine Finnsheep in Finland that differed in coat colours (white, black, brown, grey, and black/white spotted) and were from a large pedigree comprising 319 119 animals. All the individuals were genotyped with the Illumina Ovine SNP50K BeadChip by the International Sheep Genomics Consortium. We identified three genetic subpopulations that corresponded approximately with the coat colours (grey, white, and black and brown) of the sheep. We detected a significant subdivision among the colour types (F _ST = 5.4%, P<0.05). We applied robust algorithms for the genomic estimation of individual inbreeding (F _SNP) and pairwise relatedness (Φ _SNP) as implemented in the programs KING and PLINK, respectively. Estimates of the two parameters from pedigrees (F _PED and Φ _PED) were computed using the RelaX2 program. Values of the two parameters estimated from genomic and genealogical data were mostly consistent, in particular for the highly inbred animals (e.g. inbreeding coefficient F>0.0625) and pairs of closely related animals (e.g. the full- or half-sibs). Nevertheless, we also detected differences in the two parameters between the approaches, particularly with respect to the grey Finnsheep. This could be due to the smaller sample size and relative incompleteness of the pedigree for them.We conclude that the genome-wide genomic data will provide useful information on a per sample or pairwise-samples basis in cases of complex genealogies or in the absence of genealogical data.     

HLA genes in populations of the Aleutian islands

We typed a subset of the Aleut population for HLA loci (HLA-A, HLA-B, HLA-DRB1, HLA-DQB1) to obtain an HLA profile, which was compared to other Eurasian and Amerindian populations for studying Aleut origin and its significance on the peopling of the Americas. Allele frequencies at the four loci were identified in an Aleut sample using standard indirect DNA sequencing methods. Genetic distances with Amerindians and Eurasians were obtained by comparing Aleut allele frequencies with a worldwide population database (13,164 chromosomes). The most frequently extended HLA haplotypes were also calculated. We also generated Aleut relatedness dendrograms and calculated correspondence relatedness in a multidimensional scale. Both neighbor-joining dendrograms and correspondence analysis separated Aleuts from Eskimos and Amerindians. Aleuts are closer genetically to Europeans, including Scandinavians and English. Our results are concordant with those obtained by Y-chromosome analysis, suggesting that most male Aleut ancestors of our sample came mainly from Europe.     

Phylogeographic breaks in low-dispersal species: the emergence of concordance across gene trees

Computer simulations were used to investigate population conditions under which phylogeographic breaks in gene genealogies can be interpreted with confidence to infer the existence and location of historical barriers to gene flow in continuously distributed, low-dispersal species. We generated collections of haplotypic gene trees under a variety of demographic scenarios and analyzed them with regard to salient genealogical breaks in their spatial patterns. In the first part of the analysis, we estimated the frequency in which the spatial location of the deepest phylogeographic break between successive pairs of populations along a linear habitat coincided with a spatial physical barrier to dispersal. Results confirm previous reports that individual gene trees can show ‘haphazard’ phylogeographic discontinuities even in the absence of historical barriers to gene flow. In the second part of the analysis, we assessed the probability that pairs of gene genealogies from a set of population samples agree upon the location of a geographical barrier. Our findings extend earlier reports by demonstrating that spatially concordant phylogeographic breaks across independent neutral loci normally emerge only in the presence of longstanding historical barriers to gene flow. Genealogical concordance across multiple loci thus becomes a deciding criterion by which to distinguish between stochastic and deterministic causation in accounting for phylogeographic discontinuities in continuously distributed species.     

A comparison of microsatellite-based pairwise relatedness estimators

Studies of inbreeding depression or kin selection require knowledge of relatedness between individuals. If pedigree information is lacking, one has to rely on genotypic information to infer relatedness. In this study we investigated the performance (absolute and relative) of 10 marker-based relatedness estimators using allele frequencies at microsatellite loci obtained from natural populations of two bird species and one mammal species. Using Monte Carlo simulations we show that many factors affect the performance of estimators and that different sets of loci promote the use of different estimators: in general, there is no single best-performing estimator. The use of locus-specific weights turns out to greatly improve the performance of estimators when marker loci are used that differ strongly in allele frequency distribution. Microsatellite-based estimates are expected to explain between 25 and 79% of variation in true relatedness depending on the microsatellite dataset and on the population composition (i.e. the frequency distribution of relationship in the population). We recommend performing Monte Carlo simulations to decide which estimator to use in studies of pairwise relatedness.     

Estimating genealogies from unlinked marker data: a Bayesian approach

An issue often encountered in statistical genetics is whether, or to what extent, it is possible to estimate the degree to which individuals sampled from a background population are related to each other, on the basis of the available genotype data and some information on the demography of the population. In this article, we consider this question using explicit modelling of the pedigrees and gene flows at unlinked marker loci, but then restricting ourselves to a relatively recent history of the population, that is, considering the genealogy at most some tens of generations backwards in time. As a computational tool we use a Markov chain Monte Carlo numerical integration on the state space of genealogies of the sampled individuals. As illustrations of the method, we consider the question of relatedness at the level of genes/genomes (IBD estimation), using both simulated and real data.     

Relatedness and genetic structure in a socially polymorphic population of the spider Anelosimus studiosus

The evolution of sociality remains a challenge in evolutionary biology and a central question is whether association between kin is a critical factor favouring the evolution of cooperation. This study examines genetic structure of Anelosimus studiosus, a spider exhibiting polymorphic social behaviour. Two phenotypes have been identified: an ‘asocial’ phenotype with solitary female nests and a ‘social’ phenotype with multi‐female/communal nests. To address the questions of whether these phenotypes are differentiated populations and whether cooperative individuals are closely related, we used microsatellites to analyse individuals from both communal and solitary nests. We found no evidence of differentiation between social and solitary samples, implying high rates of interbreeding. This is consistent with the hypothesis that these phenotypes coexist as a behavioural polymorphism within populations. Pairwise relatedness coefficients were used to test whether cooperating individuals are more closely related than expected by chance. Pairwise relatedness of females sharing communal webs averaged 0.25, the level expected for half‐siblings and significantly more closely related than random pairs from the population. Solitary females collected at similar distances to the communal spider pairs were also more closely related than expected by chance (mean relatedness = 0.18), but less related than social pairs. These results imply that low dispersal contributes to increase likelihood of interaction between kin, but relatedness between social pairs is not explained by spatial structure alone. We propose that these phenotypes represent stages in the evolution of sociality, where viscous population structure creates opportunities for kin selection and cooperation is favoured under certain environmental conditions.     

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.     

Sibship within samples of brown trout (<Emphasis Type="Italic">Salmo trutta</Emphasis>) and implications for supportive breeding

We analysed family relationships among brown trout from two small tributary populations that have been suggested as a source of individuals for supportive breeding, using variation at eight microsatellite loci. As a control, we analysed a sample of supposedly unrelated individuals representing a large anadromous population, and we simulated unrelated individuals based on the allelic distributions in all three samples. Two different approaches were used: (1) pairwise estimates of relatedness between individuals and (2) a method for partitioning individuals into half-sib and full-sib families. The anadromous population did not show evidence of a significant number of closely related individuals. In both tributary populations, however, the distributions of pairwise relatedness estimates suggested the presence of several related individuals, and sibship reconstruction suggested fewer families consisting of more individuals than were observed for the simulated individuals. The expected increase of inbreeding coefficient in the two samples due to family structure was 0.026 and 0.030 respectively. Moreover, tests for recent bottlenecks yielded significant outcomes in both populations suggesting a history of low effective population sizes. Depending on the effective population size of captive spawners and past effective population sizes in the populations it could be beneficial to conduct sib-avoidance matings, though this cannot eliminate inbreeding but only delay it. Alternatively, individuals from different populations could be crossed. Sibship reconstruction provided the clearest evidence for family structure, but pairwise relatedness is the best measure for designing mating schemes, as it allows for mating as unrelated individuals as possible rather than just avoiding mating between sibs.