Ultra-low-density genotype panels for breed assignment of Angus and Hereford cattle

Angus and Hereford beef is marketed internationally for apparent superior meat quality attributes; DNA-based breed authenticity could be a useful instrument to ensure consumer confidence on premium meat products. The objective of this study was to develop an ultra-low-density genotype panel to accurately quantify the Angus and Hereford breed proportion in biological samples. Medium-density genotypes (13 306 single nucleotide polymorphisms (SNPs)) were available on 54 703 commercial and 4042 purebred animals. The breed proportion of the commercial animals was generated from the medium-density genotypes and this estimate was regarded as the gold-standard breed composition. Ten genotype panels (100 to 1000 SNPs) were developed from the medium-density genotypes; five methods were used to identify the most informative SNPs and these included the Delta statistic, the fixation (F_st) statistic and an index of both. Breed assignment analyses were undertaken for each breed, panel density and SNP selection method separately with a programme to infer population structure using the entire 13 306 SNP panel (representing the gold-standard measure). Breed assignment was undertaken for all commercial animals (n=54 703), animals deemed to contain some proportion of Angus based on pedigree (n=5740) and animals deemed to contain some proportion of Hereford based on pedigree (n=5187). The predicted breed proportion of all animals from the lower density panels was then compared with the gold-standard breed prediction. Panel density, SNP selection method and breed all had a significant effect on the correlation of predicted and actual breed proportion. Regardless of breed, the Index method of SNP selection numerically (but not significantly) outperformed all other selection methods in accuracy (i.e. correlation and root mean square of prediction) when panel density was ⩾300 SNPs. The correlation between actual and predicted breed proportion increased as panel density increased. Using 300 SNPs (selected using the global index method), the correlation between predicted and actual breed proportion was 0.993 and 0.995 in the Angus and Hereford validation populations, respectively. When SNP panels optimised for breed prediction in one population were used to predict the breed proportion of a separate population, the correlation between predicted and actual breed proportion was 0.034 and 0.044 weaker in the Hereford and Angus populations, respectively (using the 300 SNP panel). It is necessary to include at least 300 to 400 SNPs (per breed) on genotype panels to accurately predict breed proportion from biological samples.     

Accuracy‐driven loci selection and assignment of individuals

We present two programs: gafs for optimal selection of loci for use in individual assignment tests, and mlc , a program for individual classification using maximum likelihood and k‐nearest neighbour decision rules. gafs software employs a genetic algorithm to heuristically search multilocus subsets with several objective functions to maximize predictive accuracy of the assignments.     

Individual assignment using microsatellite DNA reveals unambiguous breed identification in the domestic dog

Modern individual clustering methods utilising hypervariable nuclear microsatellite DNA polymorphisms are being increasingly applied in the field of population genetics. This study explores the efficiency of the clustering methods in identifying the breeds of origin of 250 domestic dog (Canis familiaris) individuals based on 10 microsatellite loci. An allele sharing distance (DAS) matrix and the corresponding neighbour-joining tree of individuals revealed monophyletic assemblages that corresponded perfectly with the breeds of origin of the dogs. Individual assignment tests using a Bayesian statistical approach, an allele frequency based method, and a DCE genetic distance based method were all extremely powerful. Most strikingly, the Bayesian method provided 100% assignment success of individuals into their correct breeds of origin and 100% exclusion success of individuals from all alternate reference populations with a high level of statistical confidence (P < 0.0001). A Bayesian Markov Chain Monte Carlo clustering approach revealed clear distinction of individuals into groups according to their breeds of origin, with a near-zero level of 'genetic admixture' among breeds. The results demonstrate that an FST of 0.18, mean expected gene diversity of 0.6 across 10 loci, and approximately 50 individuals per reference population suffice to provide maximum individual assignment success in C. familiaris. This refutes the traditional view that DNA based dog breed identification is not feasible at the individual level of resolution.     

Genetic diversity of dog breeds: between-breed diversity, breed assignation and conservation approaches

Genetic relationships between 61 dog breeds were investigated, using a sampling of 1514 animals and a panel of 21 microsatellite markers. Based on the results from distance-based and Bayesian methods, breed constituted the main genetic structure, while groups including genetically close breeds showed a very weak structure. Depending on the method used, between 85.7% and 98.3% of dogs could be assigned to their breed, with large variations according to the breed. However, breed heterozygosity influenced assignment results differently according to the method used. Within-breed and between-breed diversity variations when breeds were removed were highly negatively correlated ( r  = −0.963, P  <   0.0001), because of the genetic structure of the breed set.     

Breed demarcation and potential for breed allocation of horses assessed by microsatellite markers

Population demarcation of eight horse breeds was investigated using genotype information of 306 horses from 26 microsatellite loci. The breeds include the indigenous Norwegian breeds Fjord Horse, Nordland/Lyngen Horse, Døle Horse and Coldblooded Trotter together with Icelandic Horse, Shetland Pony, Standardbred and Thoroughbred. Both phylogenetic analysis and a maximum likelihood method were applied to examine the potential for breed allocation of individual animals. The phylogenetic analysis utilizing simple allele sharing statistics revealed clear demarcation among the breeds; 95% of the individuals clustered together with animals of the same breed in the phylogenetic tree. Even breeds with a short history of divergence like Døle Horse and Coldblooded Trotter formed distinct clusters. Implementing the maximum likelihood method allocated 96% of the individuals to their source population, applying an assignment stringency of a log of the odds ratio larger than 2. Lower allocation stringency assigned nearly all the horses. Only three individuals were wrongly allocated a breed by both methods. In conclusion, the study demonstrates clear distinction among horse breeds, and by combining the two assignment methods breed allocation could be determined for more than 99% of the individuals.     

Cattle breed as a source of variation in embryo transfer

The records of 1596 embryo collections were retrospectively analysed to investigate the influence of breed on embryo production and pregnancy rates. The breed of the donor cow was a significant source of variation in the results of embryo transfer. Total embryo production per collection ranged between breeds from 6.0 to 16.2, number of transferable embryos from 2.8 to 6.6 and percent transferable from 37% to 68%. The percent of pregnancies per collection ranged from 2.0 to 4.0 and pregnancy rates from 37% to 75%. The percent of transferable embryos and the pregnancy rate was independent of the total embryos collected. Care should be taken when interpreting the breed differences, since the donor cows were not randomly representative of the breeds and were selected on different bases. Brangus (16.2), Simbrah (15.8) and Beefmaster (13.0) produced the most embryos. The most transferable embryos were collected from Simmentals (6.6), Brangus (6.6), Chianina (6.2), Beefmaster (6.1), Simbrah (6.1) and Saler (6.0) cows. Angus (68%) had the highest percent transferable, followed by Saler (64%), Chianina (60%), Limousin (52%) and Simmental (51%). Pregnancy rates were highest in the Herefords (75%), Saler (57%), Zebu (56%), Charolais (53%), Longhorns (52%) and Simbrah (50%). Estrus response to superovulation varied between breeds, but this did not account for all the breed differences in embryo production.     

Use of biosynthetic fractional 13C-labeling for backbone NMR assignment of proteins

We describe a simple approach to classify amino acid residue types in NMR spectra of proteins for supporting the backbone resonance assignments. It makes use of the differences in biosynthetic pathways of the 20 amino acids in Escherichia coli. Therefore, it is distinct from the parameters routinely exploited in the backbone resonance assignment such as chemical shifts and spin topology information. The combination of biosynthetically directed fractional 13C-labeling and uniform 15N-labeling enables us to obtain both residue-type specific information and sequential connectivities from a single protein sample. The residue-type classification exploiting biosynthetic pathways can be used for accelerating the conventional backbone assignment procedure.     

p-loci: a computer program for choosing the most efficient set of loci for parentage assignment

Determining how many and which codominant marker loci are required for accurate parentage assignment is not straightforward because levels of marker polymorphism, linkage, allelic distributions among potential parents and other factors produce differences in the discriminatory power of individual markers and sets of markers. p-loci software identifies the most efficient set of codominant markers for assigning parentage at a user-defined level of success, using either simulated or actual offspring genotypes of known parentage. Simulations can incorporate linkage among markers, mating design and frequencies of null alleles and/or genotyping errors. p-loci is available for windows systems at http://marineresearch.oregonstate.edu/genetics/ploci.htm.     

Assessing the power of informative subsets of loci for population assignment: standard methods are upwardly biased

It is well known that statistical classification procedures should be assessed using data that are separate from those used to train the classifier. This principle is commonly overlooked when the classification procedure in question is population assignment using a set of genetic markers that were chosen specifically on the basis of their allele frequencies from amongst a larger number of candidate markers. This oversight leads to a systematic upward bias in the predicted accuracy of the chosen set of markers for population assignment. Three widely used software programs for selecting markers informative for population assignment suffer from this bias. The extent of this bias is documented through a small set of simulations. The relative effect of the bias is largest when screening many candidate loci from poorly differentiated populations. Simple unbiased methods are presented and their use encouraged.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.     

Frequency and origin of haplotypes associated with the beta-globin gene cluster in individuals with trait and sickle cell anemia in the Atlantic and Pacific coastal regions of Colombia

Sickle cell anemia is a genetic disease with high prevalence in people of African descent. There are five typical haplotypes associated with this disease and the haplotypes associated with the beta-globin gene cluster have been used to establish the origin of African-descendant people in America. In this work, we determined the frequency and the origin of haplotypes associated with hemoglobin S in a sample of individuals with sickle cell anemia (HbSS) and sickle cell hemoglobin trait (HbAS) in coastal regions of Colombia. Blood samples from 71 HbAS and 79 HbSS individuals were obtained. Haplotypes were determined based on the presence of variable restriction sites within the β-globin gene cluster. On the Pacific coast of Colombia the most frequent haplotype was Benin, while on the Atlantic coast Bantu was marginally higher than Benin. Eight atypical haplotypes were observed on both coasts, being more diverse in the Atlantic than in the Pacific region. These results suggest a differential settlement of the coasts, dependent on where slaves were brought from, either from the Gulf of Guinea or from Angola, where the haplotype distributions are similar. Atypical haplotypes probably originated from point mutations that lost or gained a restriction site and/or by recombination events.     

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.     

Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution

Through thousands of years of breeding and strong human selection, the dog (Canis lupus familiaris) exists today within hundreds of closed populations throughout the world, each with defined phenotypes. A singular geographic region with broad diversity in dog breeds presents an interesting opportunity to observe potential mechanisms of breed formation. Italy claims 14 internationally recognized dog breeds, with numerous additional local varieties. To determine the relationship among Italian dog populations, we integrated genetic data from 263 dogs representing 23 closed dog populations from Italy, seven Apennine gray wolves, and an established dataset of 161 globally recognized dog breeds, applying multiple genetic methods to characterize the modes by which breeds are formed within a single geographic region. Our consideration of each of five genetic analyses reveals a series of development events that mirror historical modes of breed formation, but with variations unique to the codevelopment of early dog and human populations. Using 142,840 genome‐wide SNPs and a dataset of 1,609 canines, representing 182 breeds and 16 wild canids, we identified breed development routes for the Italian breeds that included divergence from common populations for a specific purpose, admixture of regional stock with that from other regions, and isolated selection of local stock with specific attributes.     

bels: backward elimination locus selection for studies of mixture composition or individual assignment

Methods of evaluating loci in studies of mixture composition or individual assignment are largely based on performance characteristics of individual loci. Synergisms between loci are not exploited. Loci are often evaluated based on their ability to resolve individual populations, even though multipopulation aggregations are more commonly of interest. In addition, measures of locus performance may indirectly relate to investigative objectives. A new computer program, bels, offers an alternative that addresses these limitations and may be preferable to existing methods in some applications. The algorithm is illustrated using Yukon River chum salmon (Oncorhynchus keta) data.     

Twenty-three microsatellite DNA loci for population genetic studies and parentage assignment in orangethroat darter, Etheostoma spectabile

The genus Etheostoma is a species-rich and ecologically important group of fishes in North America. The orangethroat darter (Etheostoma spectabile) is widely distributed and abundant in headwater streams throughout the central Midwest, and is an excellent model for ecological and mating system studies. We developed 23 novel, polymorphic, and independent microsatellite loci for E. spectabile. We found from two to 14 alleles per locus, and observed heterozygosities ranged from 0.39 to 1.0. These markers, in combination with others isolated from Etheostoma taxa, will be useful for ecological and evolutionary studies in the genus.     

Evaluation of factors affecting individual assignment precision using microsatellite data from horse breeds and simulated breed crosses

Assignment tests have been utilized to investigate population classification, measure genetic diversity and to solve forensic questions. Using microsatellite data from 26 loci genotyped in eight horse breeds we examined how population differentiation, number of scored loci, number of scored animals per breed and loci variability affected individual assignment precision applying log likelihood methods. We found that both genetic differentiation and number of scored loci were highly important for recognizing the breed of origin. When comparing two and two breeds, a proportion of 95% of the most differentiated breeds (0.200 < or = FST < or = 0.259) could be identified scoring only three loci, while the corresponding number was six for the least differentiated breeds (0.080 < or = FST < or = 0.139). An identical proportion of simulated breed crosses, differentiated from their parental breeds by FST estimates in the range 0.050-0.069, was identified when scoring 12 loci. This level of source identification was not obtained for the less differentiated breed crosses. The current data further suggested that population sample size and locus variability were not critical for the assignment precision as long as moderately large sample sizes (> or = 20 animals per population) and fairly variable loci were used.     

Use of autosomal loci for clustering individuals and populations of East Asian origin

We studied the genetic relationships among East Asian populations based on allele frequency differences to clarify the relative similarities of East Asian populations with a specific focus on the relationships among the Koreans, the Japanese, and the Chinese populations known to be genetically similar. The goal is to find markers appropriate for differentiating among the specific populations. In this study, no prior data existed for Koreans and the markers were selected to differentiate Chinese and Japanese. We typed, using AB TaqMan assays, single-nucleotide polymorphisms (SNPs) at 43 highly selected mostly independent diallelic sites, on 386 individuals from eight East Asian populations (Han Chinese from San Francisco, Han Chinese from Taiwan, Hakka, Koreans, Japanese, Ami, Atayal, and Cambodians) and one Siberian population (Yakut). We inferred group membership of individuals using a model-based clustering method implemented by the STRUCTURE program and population clustering by using computer programs DISTANCE, NEIGHBOR, LSSEARCH, and DRAWTREE, respectively, calculating genetic distances among populations, calculating neighbor-joining and least-squares trees, and drawing the calculated trees. On average 52% of individuals in the three Chinese groups were assigned into one cluster, and, respectively, 78 and 69% of Koreans and Japanese into a different cluster. Koreans differentiated from the Chinese groups and clustered with the Japanese in the principal component analysis (PCA) and in the best least-squares tree. The majority of Koreans were difficult to distinguish from the Japanese. This study shows that a relatively few highly selected markers can, within limits, differentiate between closely related populations.Electronic Supplementary Material Supplementary material is available for this article at     

Evidence for a geographical cline of casein haplotypes in Portuguese cattle breeds

Genetic variants of bovine milk proteins have been intensively used to characterize breeds and as markers for population/QTL studies throughout the world. However, a large number of cattle breeds including those found in Portugal, remain unstudied. In this work, we have analysed the genetic variation of six milk protein loci in 10 Portuguese cattle breeds by isoelectric focusing. High genetic diversities were generally found across breeds, with the exception of Mirandesa that showed a trend to fixation of the most common alleles in five loci, as well as of the rarer CSN3B allele. The casein haplotype BA2A was often the most frequent, followed by haplotypes BA2B and BA1A. Remarkably, CA2A was found to be the second most frequent haplotype in Southern breeds, supporting a geographical cline between Central-Northern European breeds and Bos indicus populations. Our data suggest that high genetic similarity among neighbouring Portuguese breeds is mainly caused by gene flow, and that the geographical distribution of particular casein haplotypes may indicate an influence of African cattle.     

Use of the canonical discriminant analysis to select SNP markers for bovine breed assignment and traceability purposes

Several market research studies have shown that consumers are primarily concerned with the provenance of the food they eat. Among the available identification methods, only DNA‐based techniques appear able to completely prevent frauds. In this study, a new method to discriminate among different bovine breeds and assign new individuals to groups was developed. Bulls of three cattle breeds farmed in Italy – Holstein, Brown, and Simmental – were genotyped using the 50K SNP Illumina BeadChip. Multivariate canonical discriminant analysis was used to discriminate among breeds, and discriminant analysis (DA) was used to assign new observations. This method was able to completely identify the three groups at chromosome level. Moreover, a genome‐wide analysis developed using 340 linearly independent SNPs yielded a significant separation among groups. Using the reduced set of markers, the DA was able to assign 30 independent individuals to the proper breed. Finally, a set of 48 high discriminant SNPs was selected and used to develop a new run of the analysis. Again, the procedure was able to significantly identify the three breeds and to correctly assign new observations. These results suggest that an assay with the selected 48 SNP could be used to routinely track monobreed products.     

An evaluation of allowing for mismatches as a way to manage genotyping errors in parentage assignment by exclusion

In parentage assignment by exclusion, using multiple and very polymorphic loci, genotyping errors are a major cause of non‐assignment. Using stochastic simulations, we tested the possibility to allow for mismatches at one or more allele as a way to recover assignment power. This was very efficient provided the set of loci used had a high assignment power (> 99%) and the error rate was not too high (below 3–4%). In these cases, most of the theoretical assignment power could be recovered. We also showed the efficiency of the method in a practical experiment with rainbow trout.