A new method for estimating effective population sizes from a single sample of multilocus genotypes

Equations for the effective size ( N_e ) of a population were derived in terms of the frequencies of a pair of offspring taken at random from the population being sibs sharing the same one or two parents. Based on these equations, a novel method (called sibship assignment method) was proposed to infer N_e from the sibship frequencies estimated from a sibship assignment analysis, using the multilocus genotypes of a sample of offspring taken at random from a single cohort in a population. Comparative analyses of extensive simulated data and some empirical data clearly demonstrated that the sibship assignment method is much more accurate [measured by the root mean squared error, RMSE, of 1/(2 N_e )] than other methods such as the heterozygote excess method, the linkage disequilibrium method, and the temporal method. The RMSE of 1/(2 N_e ) from the sibship assignment method is typically a small fraction of that from other methods. The new method is also more general and flexible than other methods. It can be applied to populations with nonoverlapping generations of both diploid and haplodiploid species under random or nonrandom mating, using either codominant or dominant markers. It can also be applied to the estimation of N_e for a subpopulation with immigration. With some modification, it could be applied to monoecious diploid populations with self-fertilization, and to populations with overlapping generations.     

Three new blood groups of rhesus monkeys

Three new blood group systems, called “T,” “U,” and “V,” have been identified in the rhesus monkey (Macaca mulatta). Each system consists of a single antigenic factor (blood group) detected by a monospecific alloimmune reagent that agglutinates erythrocytes. The antisera that detect these blood groups were obtained following a series of alloimmunizations and absorption fractionizations of the resulting antisera to produce operationally monospecific typing reagents. Analyses of family data indicated that each blood group was controlled by an autosomal dominant gene and that each system was independent of previously defined systems. With the addition of these new blood groups, we can identify 16 different blood group systems and well over one hundred million possible phenotypes in this species.     

COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients

The software package COANCESTRY implements seven relatedness estimators and three inbreeding estimators to estimate relatedness and inbreeding coefficients from multilocus genotype data. Two likelihood estimators that allow for inbred individuals and account for genotyping errors are for the first time included in this user-friendly program for PCs running Windows operating system. A simulation module is built in the program to simulate multilocus genotype data of individuals with a predefined relationship, and to compare the estimators and the simulated relatedness values to facilitate the selection of the best estimator in a particular situation. Bootstrapping and permutations are used to obtain the 95% confidence intervals of each relatedness or inbreeding estimate, and to test the difference in averages between groups.     

Genetic variability in the Skyros pony and its relationship with other Greek and foreign horse breeds

In Greece, seven native horse breeds have been identified so far. Among these, the Skyros pony is outstanding through having a distinct phenotype. In the present study, the aim was to assess genetic diversity in this breed, by using different types of genetic loci and available genealogical information. Its relationships with the other Greek, as well as foreign, domestic breeds were also investigated. Through microsatellite and pedigree analysis it appeared that the Skyros presented a similar level of genetic diversity to the other European breeds. Nevertheless, comparisons between DNA-based and pedigree-based results revealed that a loss of genetic diversity had probably already occurred before the beginning of breed registration. Tests indicated the possible existence of a recent bottleneck in two of the three main herds of Skyros pony. Nonetheless, relatively high levels of heterozygosity and Polymorphism Information Content indicated sufficient residual genetic variability, probably useful in planning future strategies for breed conservation. Three other Greek breeds were also analyzed. A comparison of these with domestic breeds elsewhere, revealed the closest relationships to be with the Middle Eastern types, whereas the Skyros itself remained isolated, without any close relationship, whatsoever.     

Do marker‐based paternity assignments favour heterozygous and unrelated males?

Genetic marker‐based parentage analyses are widely applied to studies of natural populations in the fields of evolutionary biology, conservation biology and ecology. When the same markers used in a parentage analysis are used together with the inferred parentage in a downstream analysis, such as the analysis of mate choice in terms of heterozygosity or relatedness, a bias may be incurred because a subset of the genotypes are favoured in parentage assignments or non‐exclusions. A previous simulation study shows that exclusion‐based paternity analyses are biased in favour of heterozygous males, and males less related to the mothers than expected under random mating. In this study, I investigated the biases of genetic paternity analyses achieved by both exclusion‐ and likelihood‐based methods, using both analytical and simulation approaches. It is concluded that while both exclusion‐ and likelihood‐based methods can lead to biased paternity assignments or non‐exclusions in favour of a subset of genotypes, the bias is not consistently towards heterozygous males or males apparently less related to mothers. Both the direction and extent of the bias depend heavily on the allele frequency distribution and the number of markers, the methods used for paternity assignments, and the estimators of relatedness. There exist important differences in the patterns of the biases between exclusion‐ and likelihood‐based paternity analysis methods. It is concluded that the markers, except when they are highly informative to yield accurate paternity assignments or exclusions, should be split into two subsets which are used separately in the paternity and downstream analyses.     

Implementation of a parentage control system in Portuguese beef-cattle with a panel of microsatellite markers

A study was conducted to assess the feasibility of applying a panel of 10 microsatellite markers in parentage control of beef cattle in Portugal. In the first stage, DNA samples were collected from 475 randomly selected animals of the Charolais, Limousin and Preta breeds. Across breeds and genetic markers, means for average number of alleles, effective number of alleles, expected heterozygosity and polymorphic information content, were 8.20, 4.43, 0.733 and 0.70, respectively. Enlightenment from the various markers differed among breeds, but the set of 10 markers resulted in a combined probability above 0.9995 in the ability to exclude a random putative parent. The marker-set thus developed was later used for parentage control in a group of 140 calves from several breeds, where there was the suspicion of possible faulty parentage recording. Overall, 76.4% of the calves in this group were compatible with the recorded parents, with most incompatibilities due to misidentification of the dam. Efforts must be made to improve the quality of pedigree information, with particular emphasis on information recorded at the calf's birth.     

Power of exclusion for parentage verification and probability of match for identity in American Kennel Club breeds using 17 canine microsatellite markers

DNA analysis of microsatellite markers has become a common tool for verifying parentage in breed registries and identifying individual animals that are linked to a database or owner. Panels of markers have been developed in canines, but their utility across and within a wide range of breeds has not been reported. The American Kennel Club (AKC) authorized a study to determine the power to exclude non-parents and identify individuals using DNA genotypes of 17 microsatellite markers in two panels. Cheek swab samples were voluntarily collected at Parent Breed Club National Specialty dog shows and 9561 samples representing 108 breeds were collected, averaging 88.5 dogs per breed. The primary panel of 10 markers exceeded 99% power of exclusion for canine parentage verification of 61% of the breeds. In combination with the secondary panel of seven markers, 100% of the tested breeds exceeded 99% power of exclusion. The minimum probability match rate of the first panel was 3.6 x 10(-5) averaged across breeds, and with the addition of the second panel, the probability match rate was 3.2 x 10(-8); thus the probability of another random, unrelated dog with the same genotype is very low. The results of this analysis indicated that, on average, the primary panel meets the AKC's needs for routine parentage testing, but that a combination of 10-15 genetic markers from the two panels could yield a universal canine panel with enhanced processing efficiency, reliability and informativeness.     

Correcting for bias in estimation of quantitative trait loci effects

Estimates of quantitative trait loci (QTL) effects derived from complete genome scans are biased, if no assumptions are made about the distribution of QTL effects. Bias should be reduced if estimates are derived by maximum likelihood, with the QTL effects sampled from a known distribution. The parameters of the distributions of QTL effects for nine economic traits in dairy cattle were estimated from a daughter design analysis of the Israeli Holstein population including 490 marker-by-sire contrasts. A separate gamma distribution was derived for each trait. Estimates for both the α and β parameters and their SE decreased as a function of heritability. The maximum likelihood estimates derived for the individual QTL effects using the gamma distributions for each trait were regressed relative to the least squares estimates, but the regression factor decreased as a function of the least squares estimate. On simulated data, the mean of least squares estimates for effects with nominal 1% significance was more than twice the simulated values, while the mean of the maximum likelihood estimates was slightly lower than the mean of the simulated values. The coefficient of determination for the maximum likelihood estimates was five-fold the corresponding value for the least squares estimates.     

一种基于高密度遗传标记的亲子鉴定方法及其应用

系谱是人类遗传及动植物育种研究与实践的重要信息来源之一。系谱记录错误是育种生产中普遍存在的一种记录错误,影响基因定位、遗传值及表型值预测等相关研究结果的可靠性。现有的方法软件可以利用遗传标记信息对疑似亲子进行亲子鉴定,但这些软件方法操作复杂,限制标记数量,如Cervus。针对当前高密度SNP标记在人类及动植物研究中广泛应用的现状,文章提出了一种基于全基因组高密度SNP数据的亲子鉴定新方法,命名为EasyPC。对EasyPC及Cervus的运行效率进行了对比,并用中国荷斯坦牛(n=2180)和杜洛克猪(n=191)的全基因组SNP芯片数据对EasyPC进行了验证。结果表明：EasyPC运行效率高于Cervus,牛和猪群体系谱错误率分别为20%和6%,与相关研究报道相符。通过使用全基因组SNP标记对群体孟德尔错误率的经验分布进行分析,该方法不仅可以简单、快速、准确地判别系谱的正确性,而且还可以对错误系谱进行校正。EasyPC为解决全基因组研究中基因型及系谱数据前处理过程中的系谱校正问题提供了一种新的途径。     

Microsatellite analysis of maternity and the mating system in the Gulf pipefish Syngnathus scovelli, a species with male pregnancy and sex-role reversal

Highly variable microsatellite loci were employed to study the mating system of the sexually dimorphic Gulf pipefish Syngnathus scovelli . In this species, like others in the family Syngnathidae, 'pregnant' males provide all parental care. Gulf pipefish were collected from one locale in the northern Gulf of Mexico, and internally carried broods of 40 pregnant males were analysed genetically. By comparing multilocus microsatellite fingerprints for the inferred mothers against expected genotypic distributions from the population sample, it was determined that: (i) only one male had received eggs from more than a single female; and (ii) on two separate occasions, two different males had received eggs from the same female. Given the high power to detect multiple matings by males, the first finding indicates that only rarely are individual males impregnated by multiple females during the course of a pregnancy. Conversely, given the lower power to detect multiple matings by females due to sampling constraints, the second finding suggests a high frequency of multiple successful matings by females. Thus, this population of Gulf pipefish displays a polyandrous genetic mating system. The relevance of these genetic findings is discussed with regard to the evolution of secondary sex traits in this species, and in other syngnathids.     

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

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

Precision of methods for calculating identity-by-descent matrices using multiple markers

A rapid, deterministic method (DET) based on a recursive algorithm and a stochastic method based on Markov Chain Monte Carlo (MCMC) for calculating identity-by-descent (IBD) matrices conditional on multiple markers were compared using stochastic simulation. Precision was measured by the mean squared error (MSE) of the relationship coefficients in predicting the true IBD relationships, relative to MSE obtained from using pedigree only. Comparisons were made when varying marker density, allele numbers, allele frequencies, and the size of full-sib families. The precision of DET was 75–99% relative to MCMC, but was not simply related to the informativeness of individual loci. For situations mimicking microsatellite markers or dense SNP, the precision of DET was ≥ 95% relative to MCMC. Relative precision declined for the SNP, but not microsatellites as marker density decreased. Full-sib family size did not affect the precision. The methods were tested in interval mapping and marker assisted selection, and the performance was very largely determined by the MSE. A multi-locus information index considering the type, number, and position of markers was developed to assess precision. It showed a marked empirical relationship with the observed precision for DET and MCMC and explained the complex relationship between relative precision and the informativeness of individual loci.     

A bayesian framework for parentage analysis: the value of genetic and other biological data

We develop fractional allocation models and confidence statistics for parentage analysis in mating systems. The models can be used, for example, to estimate the paternities of candidate males when the genetic mother is known or to calculate the parentage of candidate parent pairs when neither is known. The models do not require two implicit assumptions made by previous models, assumptions that are potentially erroneous. First, we provide formulas to calculate the expected parentage, as opposed to using a maximum likelihood algorithm to calculate the most likely parentage. The expected parentage is superior as it does not assume a symmetrical probability distribution of parentage and therefore, unlike the most likely parentage, will be unbiased. Second, we provide a mathematical framework for incorporating additional biological data to estimate the prior probability distribution of parentage. This additional biological data might include behavioral observations during mating or morphological measurements known to correlate with parentage. The value of multiple sources of information is increased accuracy of the estimates. We show that when the prior probability of parentage is known, and the expected parentage is calculated, fractional allocation provides unbiased estimates of the variance in reproductive success, thereby correcting a problem that has previously plagued parentage analyses. We also develop formulas to calculate the confidence interval in the parentage estimates, thus enabling the assessment of precision. These confidence statistics have not previously been available for fractional models. We demonstrate our models with several biological examples based on data from two fish species that we study, coho salmon (Oncorhychus kisutch) and bluegill sunfish (Lepomis macrochirus). In coho, multiple males compete to fertilize a single female's eggs. We show how behavioral observations taken during spawning can be combined with genetic data to provide an accurate calculation of each male's paternity. In bluegill, multiple males and multiple females may mate in a single nest. For a nest, we calculate the fertilization success and the 95% confidence interval of each candidate parent pair.