Performance of different SNP panels for parentage testing in two East Asian cattle breeds

Pedigree reconstruction and population structure using SNP markers in Gir cattle

Our objective was to establish a SNPs panel for pedigree reconstruction using microarrays of different densities and evaluate the genomic relationship coefficient of the inferred pedigree, in addition to analyzing the population structure based on genomic analyses in Gir cattle. For parentage analysis and genomic relationship, 16,205 genotyped Gir animals (14,458 females and 1747 males) and 1810 common markers to the four SNP microarrays were used. For population structure analyses, including linkage disequilibrium, effective population size, and runs of homozygosity (ROH), genotypes from 21,656 animals were imputed. Likelihood ratio (LR) approach was used to reconstruct the pedigree, deepening the pedigree and showing it is well established in terms of recent information. Coefficients for each relationship category of the inferred pedigree were adequate. Linkage disequilibrium showed rapid decay. We detected a decrease in the effective population size over the last 50 generations, with the average generation interval around 9.08 years. Higher ROH-based inbreeding coefficient in a class of short ROH segments, with moderate to high values, was also detected, suggesting bottlenecks in the Gir genome. Breeding strategies to minimize inbreeding and avoid massive use of few proven sires with high genetic value are suggested to maintain genetic variability in future generations. In addition, we recommend reducing the generation interval to maximize genetic progress and increase effective population size.

Performance of using opposing homozygotes for paternity testing in Japanese Black cattle

Genome-wide single nucleotide polymorphism (SNP) markers in Japanese Black cattle enable genomic prediction and verifying parent-offspring relationships. We assessed the performance of opposing homozygotes (OH) for paternity testing in Japanese Black cattle, using SNP genotype information of 50 sires and 3,420 fattened animals, 1,945 of which were fathered by the 50 genotyped sires. The number of OH was counted for each sire-progeny pair in 28,764 SNPs with minor allele frequencies of ≥0.05 in this population. Across all pairs of animals, the number of OH tended to increase as the pedigree-based coefficient of relationship decreased. With a threshold of 288 (1% of SNPs) for paternity testing, most sire-progeny pairs were detected as true relationships. The frequency of Mendelian inconsistencies was 2.4%, reflecting the high accuracy of pedigree information in Japanese Black cattle population. The results indicate the utility of OH for paternity testing in Japanese Black cattle.

SNP panels for the estimation of dairy breed proportion and parentage assignment in African crossbred dairy cattle

BACKGROUND

Understanding the relationship between genetic admixture and phenotypic performance is crucial for the optimization of crossbreeding programs. The use of small sets of informative ancestry markers can be a cost-effective option for the estimation of breed composition and for parentage assignment in situations where pedigree recording is difficult. The objectives of this study were to develop small single nucleotide polymorphism (SNP) panels that can accurately estimate the total dairy proportion and assign parentage in both West and East African crossbred dairy cows.

METHODS

Medium- and high-density SNP genotype data (Illumina BovineSNP50 and BovineHD Beadchip) for 4231 animals sampled from African crossbreds, African Bos taurus, European Bos taurus, Bos indicus, and African indigenous populations were used. For estimating breed composition, the absolute differences in allele frequency were calculated between pure ancestral breeds to identify SNPs with the highest discriminating power, and different combinations of SNPs weighted by ancestral origin were tested against estimates based on all available SNPs. For parentage assignment, informative SNPs were selected based on the highest minor allele frequency (MAF) in African crossbred populations assuming two Scenarios: (1) parents were selected among all the animals with known genotypes, and (2) parents were selected only among the animals known to be a parent of at least one progeny.

RESULTS

For the medium-density genotype data, SNPs selected for the largest differences in allele frequency between West African indigenous and European Bos taurus breeds performed best for most African crossbred populations and achieved a prediction accuracy (r2) for breed composition of 0.926 to 0.961 with 200 SNPs. For the high-density dataset, a panel with 70% of the SNPs selected on their largest difference in allele frequency between African and European Bos taurus performed best or very near best across all crossbred populations with r2 ranging from 0.978 to 0.984 with 200 SNPs. In all African crossbred populations, unambiguous parentage assignment was possible with ≥ 300 SNPs for the majority of the panels for Scenario 1 and ≥ 200 SNPs for Scenario 2.

CONCLUSIONS

The identified low-cost SNP assays could overcome incomplete or inaccurate pedigree records in African smallholder systems and allow effective breeding decisions to produce progeny of desired breed composition.

Selection and implementation of single nucleotide polymorphism markers for parentage analysis in crossbred cattle population

Crossbreeding is an essential way of improving herd performance. However, frequent parentage record errors appear, which results in the lower accuracy of genetic parameter estimation and genetic evaluation. This study aims to build a single nucleotide polymorphism (SNP) panel with sufficient power for parentage testing in the crossbred population of Simmental and Holstein cattle. The direct sequencing technique in PCR products of pooling DNA along with matrix-assisted laser desorption/ionization time-of-flight MS method for genotyping the individuals was applied. A panel comprising 50 highly informative SNPs for parentage analysis was developed in the crossbred population. The average minor allele frequency for SNPs was 0.43, and the cumulative probability of exclusion for single-parent and both-parent inference met 0.99797 and 0.999999, respectively. The maker-set for parentage verification was then used in a group of 81 trios with aid of the likelihood-based parentage-assignment program of Cervus software. Reconfirmation with on-farm records showed that this 50-SNP system could provide sufficient and reliable information for parentage testing with the parental errors for mother-offspring and sire-offspring being 8.6 and 18.5%, respectively. In conclusion, a set of low-cost and efficient SNPs for the paternity testing in the Simmental and Holstein crossbred population are provided.

Evaluation of the International Society for Animal Genetics bovine single nucleotide polymorphism parentage panel in South African Bonsmara and Drakensberger cattle

A panel of 200 single nucleotide polymorphisms (SNPs) have been recommended by the International Society for Animal Genetics (ISAG) for use in parentage verification of cattle. While the SNPs included on the ISAG panel are segregating in European Bos taurus and Bos indicus breeds, their applicability in South African (SA) Sanga cattle has never been evaluated. This study, therefore, assessed the usefulness of the ISAG panel in SA Bonsmara (BON) and Drakensberger (DRB) cattle. Genotypes of 185 ISAG SNPs from 64 BON and 97 DRB sire-offspring pairs were available, all of which were validated with 119,375 SNPs. Of the 185 ISAG SNPs, 14 and 18 in the BON and DRB, respectively (9 in common to both breeds), were either monomorphic, exhibited at least one discordance between validated sire-offspring pairs, or had poor call rate or clustering issue. The mean minor allele frequency of the 185 ISAG SNPs was 0.331 in the BON and 0.359 in the DRB. The combined probability of parentage exclusion (P_E) was the same (99.46%) for both breeds, while the probability of identity varied from 1.61 × 10^-48 (BON) to 1.11 × 10^-54 (DRB). Fifteen (23.4%) and 32 (33%) of the already validated sire-offspring pairs for the BON and DRB, respectively, were determined by the ISAG panel to be false-negatives based on a threshold of having at least two discordant SNPs. In comparison to sire discovery using the 119,375 SNPs, sire discovery using only the ISAG panel identified correctly 44 (out of 64 identified using the 119,375 SNPs) unique sire-offspring BON pairs and 62 (out of 97 identified using the 119,375 SNPs) unique sire-offspring DRB when all sires were masked. Five BON and three DRB offspring had > 1 sire nominated. This study demonstrated that the use of the ISAG panel may result in incorrect exclusions and multiple candidate sires for a given animal. Selection of more informative SNPs is, therefore, necessary in the pursuit of a low-cost and effective SNP panel for indigenous cattle breeds in SA.

Genotyping-free parentage assignment using RAD-seq reads

Parentage assignment is defined as the identification of the true parents of one focal offspring among a list of candidates and has been commonly used in zoological, ecological, and agricultural studies. Although likelihood-based parentage assignment is the preferred method in most cases, it requires genotyping a predefined set of DNA markers and providing their population allele frequencies. In the present study, we proposed an alternative method of parentage assignment that does not depend on genotype data and prior information of allele frequencies. Our method employs the restriction site-associated DNA sequencing (RAD-seq) reads for clustering into the overlapped RAD loci among the compared individuals, following which the likelihood ratio of parentage assignment could be directly calculated using two parameters-the genome heterozygosity and error rate of sequencing reads. This method was validated on one simulated and two real data sets with the accurate assignment of true parents to focal offspring. However, our method could not provide a statistical confidence to conclude that the first ranked candidate is a true parent.

The future of parentage analysis: From microsatellites to SNPs and beyond

Parentage analysis is a cornerstone of molecular ecology that has delivered fundamental insights into behaviour, ecology and evolution. Microsatellite markers have long been the king of parentage, their hypervariable nature conferring sufficient power to correctly assign offspring to parents. However, microsatellite markers have seen a sharp decline in use with the rise of next-generation sequencing technologies, especially in the study of population genetics and local adaptation. The time is ripe to review the current state of parentage analysis and see how it stands to be affected by the emergence of next-generation sequencing approaches. We find that single nucleotide polymorphisms (SNPs), the typical next-generation sequencing marker, remain underutilized in parentage analysis but are gaining momentum, with 58 SNP-based parentage analyses published thus far. Many of these papers, particularly the earlier ones, compare the power of SNPs and microsatellites in a parentage context. In virtually every case, SNPs are at least as powerful as microsatellite markers. As few as 100-500 SNPs are sufficient to resolve parentage completely in most situations. We also provide an overview of the analytical programs that are commonly used and compatible with SNP data. As the next-generation parentage enterprise grows, a reliance on likelihood and Bayesian approaches, as opposed to strict exclusion, will become increasingly important. We discuss some of the caveats surrounding the use of next-generation sequencing data for parentage analysis and conclude that the future is bright for this important realm of molecular ecology.

SNP Data Quality Control in a National Beef and Dairy Cattle System and Highly Accurate SNP Based Parentage Verification and Identification

A major use of genetic data is parentage verification and identification as inaccurate pedigrees negatively affect genetic gain. Since 2012 the international standard for single nucleotide polymorphism (SNP) verification in Bos taurus cattle has been the ISAG SNP panels. While these ISAG panels provide an increased level of parentage accuracy over microsatellite markers (MS), they can validate the wrong parent at ≤1% misconcordance rate levels, indicating that more SNP are needed if a more accurate pedigree is required. With rapidly increasing numbers of cattle being genotyped in Ireland that represent 61 B. taurus breeds from a wide range of farm types: beef/dairy, AI/pedigree/commercial, purebred/crossbred, and large to small herd size the Irish Cattle Breeding Federation (ICBF) analyzed different SNP densities to determine that at a minimum ≥500 SNP are needed to consistently predict only one set of parents at a ≤1% misconcordance rate. For parentage validation and prediction ICBF uses 800 SNP (ICBF800) selected based on SNP clustering quality, ISAG200 inclusion, call rate (CR), and minor allele frequency (MAF) in the Irish cattle population. Large datasets require sample and SNP quality control (QC). Most publications only deal with SNP QC via CR, MAF, parent-progeny conflicts, and Hardy-Weinberg deviation, but not sample QC. We report here parentage, SNP QC, and a genomic sample QC pipelines to deal with the unique challenges of >1 million genotypes from a national herd such as SNP genotype errors from mis-tagging of animals, lab errors, farm errors, and multiple other issues that can arise. We divide the pipeline into two parts: a Genotype QC and an Animal QC pipeline. The Genotype QC identifies samples with low call rate, missing or mixed genotype classes (no BB genotype or ABTG alleles present), and low genotype frequencies. The Animal QC handles situations where the genotype might not belong to the listed individual by identifying: >1 non-matching genotypes per animal, SNP duplicates, sex and breed prediction mismatches, parentage and progeny validation results, and other situations. The Animal QC pipeline make use of ICBF800 SNP set where appropriate to identify errors in a computationally efficient yet still highly accurate method.

Single nucleotide polymorphisms for DNA typing in the domestic horse

Genetic markers are important resources for individual identification and parentage assessment. Although short tandem repeats (STRs) have been the traditional DNA marker, technological advances have led to single nucleotide polymorphisms (SNPs) becoming an attractive alternative. SNPs can be highly multiplexed and automatically scored, which allows for easier standardization and sharing among laboratories. Equine parentage is currently assessed using STRs. We obtained a publicly available SNP dataset of 729 horses representing 32 diverse breeds. A proposed set of 101 SNPs was analyzed for DNA typing suitability. The overall minor allele frequency of the panel was 0.376 (range 0.304-0.419), with per breed probability of identities ranging from 5.6 × 10^-35 to 1.86 × 10^-42 . When one parent was available, exclusion probabilities ranged from 0.9998 to 0.999996, although when both parents were available, all breeds had exclusion probabilities greater than 0.9999999. A set of 388 horses from 35 breeds was genotyped to evaluate marker performance on known families. The set included 107 parent-offspring pairs and 101 full trios. No horses shared identical genotypes across all markers, indicating that the selected set was sufficient for individual identification. All pairwise comparisons were classified using ISAG rules, with one or two excluding markers considered an accepted parent-offspring pair, two or three excluding markers considered doubtful and four or more excluding markers rejecting parentage. The panel had an overall accuracy of 99.9% for identifying true parent-offspring pairs. Our developed marker set is both present on current generation SNP chips and can be highly multiplexed in standalone panels and thus is a promising resource for SNP-based DNA typing.

Genetic tests for estimating dairy breed proportion and parentage assignment in East African crossbred cattle

BACKGROUND

Smallholder dairy farming in much of the developing world is based on the use of crossbred cows that combine local adaptation traits of indigenous breeds with high milk yield potential of exotic dairy breeds. Pedigree recording is rare in such systems which means that it is impossible to make informed breeding decisions. High-density single nucleotide polymorphism (SNP) assays allow accurate estimation of breed composition and parentage assignment but are too expensive for routine application. Our aim was to determine the level of accuracy achieved with low-density SNP assays.

METHODS

We constructed subsets of 100 to 1500 SNPs from the 735k-SNP Illumina panel by selecting: (a) on high minor allele frequencies (MAF) in a crossbred population; (b) on large differences in allele frequency between ancestral breeds; (c) at random; or (d) with a differential evolution algorithm. These panels were tested on a dataset of 1933 crossbred dairy cattle from Kenya/Uganda and on crossbred populations from Ethiopia (N = 545) and Tanzania (N = 462). Dairy breed proportions were estimated by using the ADMIXTURE program, a regression approach, and SNP-best linear unbiased prediction, and tested against estimates obtained by ADMIXTURE based on the 735k-SNP panel. Performance for parentage assignment was based on opposing homozygotes which were used to calculate the separation value (sv) between true and false assignments.

RESULTS

Panels of SNPs based on the largest differences in allele frequency between European dairy breeds and a combined Nelore/N'Dama population gave the best predictions of dairy breed proportion (r² = 0.962 to 0.994 for 100 to 1500 SNPs) with an average absolute bias of 0.026. Panels of SNPs based on the highest MAF in the crossbred population (Kenya/Uganda) gave the most accurate parentage assignments (sv = -1 to 15 for 100 to 1500 SNPs).

CONCLUSIONS

Due to the different required properties of SNPs, panels that did well for breed composition did poorly for parentage assignment and vice versa. A combined panel of 400 SNPs was not able to assign parentages correctly, thus we recommend the use of 200 SNPs either for breed proportion prediction or parentage assignment, independently.

Development of a SNP panel dedicated to parentage assignment in French sheep populations

BACKGROUND

The efficiency of breeding programs partly relies on the accuracy of the estimated breeding values which decreases when pedigrees are incomplete. Two reproduction techniques are mainly used by sheep breeders to identify the sires of lambs: animal insemination and natural matings with a single ram per group of ewes. Both methods have major drawbacks, notably time-consuming tasks for breeders, and are thus used at varying levels in breeding programs. As a consequence, the percentage of known sires can be very low in some breeds and results in less accurate estimated breeding values.

RESULTS

In order to address this issue and offer an alternative strategy for obtaining parentage information, we designed a set of 249 SNPs for parentage assignment in French sheep breeds and tested its efficiency in one breed. The set was derived from the 54 K SNP chip that was used to genotype the thirty main French sheep populations. Only SNPs in Hardy-Weinberg equilibrium, displaying the highest Minor Allele Frequency across all the thirty populations and not associated with Mendelian errors in verified family trios were selected. The panel of 249 SNPs was successfully used in an on-farm test in the BMC breed and resulted in more than 95% of lambs being assigned to a unique sire.

CONCLUSION

In this study we developed a SNP panel for assignment that achieved good results in the on-farm testing. We also raised some conditions for optimal use of this panel: at least 180 SNPs should be used and a minute preparation of the list of candidate sires. Our panel also displays high levels of MAF in the SheepHapMap breeds, particularly in the South West European breeds.

A method for single nucleotide polymorphism selection for parentage assessment in goats

Accurate pedigrees are essential to optimize genetic improvement and conservation of animal genetic resources. In goats, the use of mating groups and kidding management procedures hamper the identification of parentage. Small panels of single nucleotide polymorphisms (SNP) have been proposed in other species to substitute microsatellites for parentage assessment. Using data from the current GoatSNP50 chip, we developed a new 3-step procedure to identify a low-density SNP panel for highly accurate parentage assessment. Methodologies for SNP selection used in other species are less suitable in the goat because of uncertainties in the genome assembly. The procedure developed in this study is based on parent-offspring identification and on estimation of Mendelian errors, followed by canonical discriminant analysis identification and stepwise regression reduction. Starting from a reference sample of 109 Alpine goats with known pedigree relationships, we first identified a panel of 200 SNP that was further reduced to 2 final panels of 130 and 114 SNP with random coincidental match inclusion of 1.51×10(-57) and 2.94×10(-34), respectively. In our reference data set, all panels correctly identified all parent-offspring combinations, revealing a 40% pedigree error rate in the information provided by breeders. All reference trios were confirmed by official tests based on microsatellites. Panels were also tested on Saanen and Teramana breeds. Although the testing on a larger set of breeds in the reference population is still needed to validate these results, our findings suggest that our procedure could identify SNP panels for accurate parentage assessment in goats or in other species with unreliable marker positioning.

Recent development of allele frequencies and exclusion probabilities of microsatellites used for parentage control in the German Holstein Friesian cattle population

Background

Methods for parentage control in cattle have changed since their initial implementation in the late 1950’s from blood group typing to more current single nucleotide polymorphism determination. In the early 1990’s, 12 microsatellites were selected by the International Society for Animal Genetics based on their informativeness and robustness in a variety of different cattle breeds. Since then this panel is used as standard in cattle herd book breeding and its application is accompanied by recurrent international comparison tests ensuring permanent validity for the most common commercial dairy and beef cattle breeds for example Holstein Friesian, Simmental, Angus, and Hereford. Although, nearly every parentage can be resolved using these microsatellites, cases with very close relatives became an emerging resolution problem during recent years. This is mainly due to an increase of monomorphism and a trend to the fixation of alleles, although no direct selection against their variability was applied. Thus other effects must be presumed resulting in a loss of polymorphism information content, heterozygosity, and exclusion probabilities.

Results

To determine changes of allele frequencies and exclusion probabilities, we analyzed the development of these parameters for the 12 microsatellites from 2004 to 2014. One hundred sixty eight thousand recorded Holstein Friesian cattle genotypes were evaluated. During this period certain alleles of nine microsatellites increased significantly (t-values >5). When calculating the exclusion probabilities for 11 microsatellites, reduction was determined for the three situations, i.e. one parent is wrongly identified (p = 0.01), both parents are wrongly identified (p = 0.005), and the genotype of one parent is missing (p = 0.048). With the addition of BM1818 to the marker set in 2009, this development was corrected leading to significant increases in exclusion probabilities. Although, the exclusion probabilities for the three family situations using the 12 microsatellites are >99 %, the clarification of 142 relationships in 40,000 situations where one parent is missing will still be impossible.

Twenty-five sires were identified that are responsible for the most significant microsatellite allele increases in the population. The corresponding alleles are mainly associated with milk protein and fat yield, body weight at birth and weaning, as well as somatic cell score, milk fat percentage, and longissimus muscle area.

Conclusions

Our data show that most of the microsatellites used for parentage control in cattle show directional changes in allele frequencies consistent with the history of artificial selection in the German Holstein population.

How many markers are enough? Factors influencing parentage testing in different livestock populations

Reliability of parentage test panels is usually based on its power to exclude wrong parentage assignments based on allele frequencies. We evaluated the rates of false exclusions and inclusions in parentage assignments, and how these results are affected by allele frequencies, panel sizes and the number of allowed mismatches. We also evaluated the reliability of parentage testing by comparing populations with distinct genetic backgrounds using pure and composite families of cattle and sheep. Allowing for 1% genotype mismatches in true parent-offspring relations provided the best compromise between false-positive and false-negative assignments. Pure breeds needed at least 200-210 single-nucleotide polymorphism (SNP) markers to correctly assign relations, but between 700 and 890 markers to avoid assigning incorrect relationships. Composite breeds needed between 220 (sheep) and 500 (cattle) markers for correct assignment; 680 (cattle) to 4400 (sheep) SNPs were needed to eliminate false-positive assignments. Allowing 0% genotype mismatches decreased false-positive but increased false-negative assignments, whilst a higher threshold of 2% showed the opposite effects. Panels with high minor allele frequencies (0.35-0.45) provided the best chance for correct parentage resolutions requiring fewer markers. Further, we propose that a dynamic threshold would allow adapting to population specific error rates. A comparison to the performance of the official International Society for Animal Genetics SNP panel for cattle and a recently published SNP panel for sheep showed that randomly selected markers performed only slightly worse for the applied parentage test based on opposing homozygotes. This suggests that even with carefully selected panels, only marginal assignment improvements are obtainable for a particular number of SNPs. The main point for improvement is the number of markers used. We recommend using at least 200 SNP markers for parentage testing if the aim is to reduce false-negative results. To fully exclude false positives at least 700 markers are required.

Analytical and statistical consideration on the use of the ISAG-ICAR-SNP bovine panel for parentage control, using the Illumina BeadChip technology: example on the German Holstein population

Background

Parentage control is moving from short tandem repeats- to single nucleotide polymorphism (SNP) systems. For SNP-based parentage control in cattle, the ISAG-ICAR Committee proposes a set of 100/200 SNPs but quality criteria are lacking. Regarding German Holstein-Friesian cattle with only a limited number of evaluated individuals, the exclusion probability is not well-defined. We propose a statistical procedure for excluding single SNPs from parentage control, based on case-by-case evaluation of the GenCall score, to minimize parentage exclusion, based on miscalled genotypes. Exclusion power of the ISAG-ICAR SNPs used for the German Holstein-Friesian population was adjusted based on the results of more than 25 000 individuals.

Results

Experimental data were derived from routine genomic selection analyses of the German Holstein-Friesian population using the Illumina BovineSNP50 v2 BeadChip (20 000 individuals) or the EuroG10K variant (7000 individuals). Averages and standard deviations of GenCall scores for the 200 SNPs of the ISAG-ICAR recommended panel were calculated and used to calculate the downward Z-value. Based on minor allelic frequencies in the Holstein-Friesian population, one minus exclusion probability was equal to 1.4×10⁻¹⁰ and 7.2×10⁻²⁶, with one and two parents, respectively. Two monomorphic SNPs from the 100-SNP ISAG-ICAR core-panel did not contribute. Simulation of 10 000 parentage control combinations, using the GenCall score data from both BeadChips, showed that with a Z-value greater than 3.66 only about 2.5% parentages were excluded, based on the ISAG-ICAR recommendations (core-panel: ≥ 90 SNPs for one, ≥ 85 SNPs for two parents). When applied to real data from 1750 single parentage assessments, the optimal threshold was determined to be Z = 5.0, with only 34 censored cases and reduction to four (0.2%) doubtful parentages. About 70 parentage exclusions due to weak genotype calls were avoided, whereas true exclusions (n = 34) were unaffected.

Conclusions

Using SNPs for parentage evaluation provides a high exclusion power also for parent identification. SNPs with a low GenCall score show a high tendency towards intra-molecular secondary structures and substantially contribute to false exclusion of parentages. We propose a method that controls this error without excluding too many parent combinations from the evaluation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0085-1) contains supplementary material, which is available to authorized users.