Skip to main content

The effect of mislabeled phenotypic status on the identification of mutation-carriers from SNP genotypes in dairy cattle



Statistical and machine learning applications are increasingly popular in animal breeding and genetics, especially to compute genomic predictions for phenotypes of interest. Noise (errors) in the data may have a negative impact on the accuracy of predictions. The effects of noisy data have been investigated in genome-wide association studies for case–control experiments, and in genomic predictions for binary traits in plants. No studies have been published yet on the impact of noisy data in animal genomics. In this work, the susceptibility to noise of five classification models (Lasso-penalised logistic regression—Lasso, K-nearest neighbours—KNN, random forest—RF, support vector machines with linear—SVML—or radial—SVMR—kernel) was tested. As illustration, the identification of carriers of a recessive mutation in cattle (Bos taurus) was used. A population of 3116 Fleckvieh animals with SNP genotypes on the same chromosome as the mutation locus (BTA 19) was available. The carrier status (0/1 phenotype) was randomly sampled to generate noise. Increasing proportions of noise—up to 20%— were introduced in the data.


SVMR and Lasso were relatively more robust to noise in the data, with total accuracy still above 0.975 and TPR (true positive rate; accuracy in the minority class) in the range 0.5–0.80 also with 17.5–20% mislabeled observations. The performance of SVML and RF decreased monotonically with increasing noise in the data, while KNN constantly failed to identify mutation carriers (observations in the minority class). The computation time increased with noise in the data, especially for the two support vector machines classifiers.


This work was the first to assess the impact of phenotyping errors on the accuracy of genomic predictions in animal genetics. The choice of the classification method can influence results in terms of higher or lower susceptibility to noise. In the presented problem, SVM with radial kernel performed relatively well even when the proportion of errors in the data reached 12.5%. Lasso was the second best method, while SVML, RF and KNN were very sensitive to noise. Taking into account both accuracy and computation time, Lasso provided the best combination.


In data science, statistical and machine learning approaches are used to identify patterns within data, with the primary objective of making predictions on future or unobserved data. Their popularity has increased with the size of available data: the advent of “big data” [1] has outdated many classical data analysis and statistical approaches. From a search on Google Scholar the number of publications related to statistical and machine learning increased from 10,690 in year 2000 to 1,211,400 in year 2016, with a peak rate between years 2011 and 2013, to then continue to increase at a slower pace (Fig. 1). Statistical and machine learning are nowadays applied to many different areas like Web Search, spam filters, recommender systems, ad placement, credit scoring, fraud detection [2,3,4,5], and to diverse biological disciplines, like drug development, DNA sequence analysis, cell biology and animal genetics [6,7,8,9].

A common learning task is classification (binomial or multinomial), where the objective is to build a classifier that can correctly predict the class of a new object given some training examples of known objects [10].

Machine learning methods may be susceptible to biases, especially if we consider that the training data can contain errors. Errors in the data are known as noise, and can arise because of different reasons (e.g. instrument errors, quantization errors, environmental noise, model mis-specification, human errors, inherent randomness in the physical processes): the consequence is that the classifier learns from a distorted version of the actual data and its predictive ability will be biased upwards or downwards, or randomly unreliable [11, 12].

In the field of genomics, genotypes are typically used together with phenotypes to either detect associations or make whole-genome predictions [13, 14]. Errors may be found in the genotypic and/or in the phenotypic data. The consequences of genotyping errors [15, 16], and of errors in the imputation of missing genotypes [17,18,19,20] on genome-wide association studies and genomic predictions have been addressed. Scientific literature on phenotypic errors in genomics is much scarcer. The effect of phenotype misclassifications on the statistical power of genome-wide association studies (GWAS) has been addressed in case–control studies in human medicine [21]. More recently, the influence of noisy data on the accuracy of whole-genome predictions has been examined in sugar beets [22]. No studies have been published yet on the impact of noisy phenotypes on genome-enabled predictions in human or animal genomics.

In this paper, the impact of randomly mislabeled observations on the accuracy of genomic predictions for binary traits is investigated. A cattle (Bos taurus) population with known carrier/non-carrier status for a harmful recessive genetic mutation was used for illustration. SNP genotypes were used to classify animals. Starting from a dataset with known mutation carrier status (no errors), increasing proportions of noisy labels were randomly generated, and the performance of different classification methods was measured.

Fig. 1

Cumulative number of machine and statistical learning-related publications over time. From Google Scholar queries for “machine learning” and “statistical learning” publications between 2000 and 2016 [machine learning (ML) solid gold line, statistical learning (SL) dotted red line, total dashed dark gray line]


Experimental data

SNP genotypes and mutation carrier status were available for a dairy cattle population of 3116 Fleckvieh animals. The mutation used for illustration is the TUBD1 recessive mutation [23, 24] at the beginning of BTA19 (Bos taurus autosome 19). The TUBD1 mutation (and the associated BH2 haplotype: [25]) were reported to be associated with stillbirth and calf survival rate [26]. Animals were labeled as carriers or not of the mutation (coded as 1 or 0). There were 126 carriers (4.04%) and 2990 non-carriers (95.96%). All animals were genotyped with the Bovine SNP50 v2 (54K) Illumina BeadChip. Only the 1512 SNPs on BTA19 were used for the analysis. No individual animal had a call-rate lower than 95%. SNPs with a call-rate lower than 95% (195 SNP) were removed from the analysis. Residual missing SNP genotypes were imputed based on linkage disequilibrium, using the localized haplotype clustering imputation method implemented in the computer package “Beagle” v.3 [27]. Data for the present study were provided by ZuchData EDV-Dienstleistungen GmbH (Austria).

Classification models

Five machine learning (ML) algorithms were used to identify mutation carriers from SNP genotypes: Lasso-penalised logistic regression (Lasso), Support Vector Machines using either a linear (SVML) or a radial (SVMR) kernel, K-nearest neighbours (KNN) and random forest (RF). In order to explore the effect of noisy labels on genomic classifications, 10 different scenarios were simulated. In each scenario an increasing proportion of noise was introduced by flipping independently the original carrier state. The following noise proportions were tested: 0% (original data with no errors), 1, 2.5, 5, 7.5, 10, 12.5, 15, 17.5 and 20%. For each proportion of mislabelled observations, the five classification models were tested.

Lasso-penalised logistic regression (Lasso)

The probability of carrying the mutation (\(P(Y = 1|X) = p(x)\)) was modeled as a linear combination of SNP genotypes in a logistic regression model:

$$\begin{aligned} logit(p(x_i))=\mu +\sum _{j=1}^m z_{ij}SNP_j \end{aligned}$$

where p(x) is the \(P(Y = 1|X)\) for individual i with vector of SNP genotypes \(x_i\); S\(NP_j\) is the effect of the \(j\)th marker; \(z_{ij}\) is the genotype of individual i at locus j (0, 1 or 2 for AA, AB and BB genotypes). The model in Eq. 1 was fitted by maximizing the corresponding Lasso-penalized log likelihood function [28]. The tuning parameter \(\lambda\) controls the degree of regularization, and was specified through cross-validation. Logistic regression returns the log-odds of p(x) which are back-transformed to \(P(Y = 1|X)\) through the cumulative distribution function of the logistic distribution (i.e. the logistic function). Individuals with \(p(x) > / < 0.5\) were classified as carriers or not of the mutation.

Support vector machines (SVM)

Two support vect//.oor machines (SVM) models were fitted for the classification of carriers and non-carriers of the mutation: with linear (SVML) and radial (SVMR) kernel functions. SVM maps the vector of SNP genotypes \(x \in \mathcal {R}\) into a higher dimensional feature space \(\phi (x) \in \mathcal {H}\) and constructs a decision boundary which is linear in \(\mathcal {H}\), and possibly non-linear in \(\mathcal {H}\). Animals are then classified into carriers and non-carriers of the mutation based on the width of the margin \(\mathcal {M}\) and the sign of the classifier:

$$\begin{aligned} f(x)=\beta _0 + \sum _{i=1}^n \alpha _i K(x,x_i) \end{aligned}$$

The kernel function K has the form \(K(x_i,x_{i'})=\sum _{j=1}^m x_{ij}x_{i'j}\) in SVML and \(K(x_i,x_{i'})=exp \left( -\gamma \sum _{j=1}^m (x_{ij}-x_{i'j})^2 \right)\) in SVMR. The hyperparameters C (which controls the width of the margin \(\mathcal {M}\)) and \(\gamma\) (which controls the degree of non-linearity in SVMR) were chosen so to minimize the classification error through cross-validation in the training set. A full description of SVM can be found in [29].

K-nearest neighbours (KNN)

The predicted carrier/non-carrier status for animal \(x_0\) was obtained by majority vote among the K closest neighbours. The neighbourhood was determined by Euclidean distances based on SNP genotypes, for each neighbour i over m SNP dimensions:

$$\begin{aligned} D_E=d(x_0,x_i)=\sqrt{\sum _{j=1}^m (x_{0j}-x_{ij})^2} \end{aligned}$$

The size of the neighbourhood K was determined through cross-validation in the training data.

Random forest (RF) classifier

A large number of classification trees was built on \(B=500\) bootstrapped samples of the data. Classification trees were decorrelated by using, at each node, a random subset s of the 1512 SNPs on BTA19. The size of the random feature subset s was optimized around \(\sqrt{1512} \approx 39\) SNPs The final classifier was obtained by majority vote over the B classification trees:

$$\begin{aligned} \hat{f}_{avg}(x_i)=\frac{1}{B}\sum _{b=1}^B I\left(\hat{f}_b(x_i)=[0/1]\right) \end{aligned}$$

where \(x_i\) is the vector of SNP genotypes for animal i, and \(\hat{f}_b(x_i)\) is the prediction (carrier/non-carrier) from the classification tree built on the \(b\)th bootstrapped data sample. More details on random forest can be found in [30].

Prediction accuracy

In order to compare the predictive ability of the five classifiers, the data were initially split in a training and a testing data set: \(70\%\) of the observations used for training, \(30\%\) of the observations used for testing. The training dataset (which contained increasing proportions of random noise) was used to tune the hyperparameters (\(\lambda\) for Lasso; C and \(\gamma\) for SVML and SVMR; K for KNN; s for RF) and train the classifier through a 10-fold cross-validation procedure: the hyperparameters that gave the lowest average balanced accuracy in the validation sets (the 10th fold, in turn) were selected. The final model was then applied to the testing set to predict the original carrier-non carrier status and measure the accuracy of classification. Prior to fitting the model, monomorphic and collinear (correlation >0.99) SNPs were edited out of the training set, to remove non-informative and redundant predictors and avoid problems due to linear dependencies. This procedure was repeated 10 times per each proportion of noise (0–20%), using different training and testing subsets each time. The following measures of classification accuracy were calculated in the testing data set: (1) accuracy (ACC): the proportion of the total number of correct predictions over the total test sample size; (2) true positive rate (TPR, sensitivity): the proportion of mutation carriers (positives) that were correctly identified, over the total number of carriers in the test set; and (3) true negative rate (TNR, specificity), the proportion of non-carriers (negatives) that were correctly identified over the total number of non-carriers in the test set. Results were averaged over replicates by noise proportion.


Data preparation and editing, and all statistical analysis were performed using the R programming environment v.3.2.3 [31], except missing genotype imputation, which was carried out with the computer package “Beagle” v.3.3.2 [27]. The R packages glmnet [32], e1071 [33], class and caret [34] were used to fit the Lasso logistic regression, SVM with linear and radial kernels, KNN and RF classification models. The analyses were run on the bioinformatics platform at PTP Science Park (, which includes a high performance computing cluster with 600 CPUs, 2.5 TB of RAM and 100 TB of storage space for archiving and back-up.


The total prediction accuracy for the five classification methods over the ten proportions of random errors introduced in the data is shown in Fig. 2. Total accuracy (ACC) was above 95% for all methods and proportions of errors. Lasso and SVML reached 100% accuracy with no errors in the data. When errors began to be introduced, the accuracy of SVML decreased, down to 95.02% with 12.5% errors. For Lasso and SVMR, ACC was above 99% from 0 to 12.5% errors in the data, dropping to 98% for 15 and 17.5% errors, and eventually relapsing back above 99% with 20% errors in the training set. With RF, ACC was 98.9% with no errors in the data and went down to 96.1% with 20% mislabelings. KNN gave a lower average ACC (95–96%) which remained fairly constant over different percentages of noise in the data.

Fig. 2

Total prediction accuracy (ACC) as a function of noise in the data. Proportion of observations (both carriers and non-carriers of the mutation) that were correctly identified by the five classification methods over the 10 proportions of errors introduced. Lasso black, SVML green, SVMR light blue, RF blue, KNN red

The TNR (specificity) and TPR (sensitivity) for the five classification methods over the 10 proportions of errors are shown in Fig. 3. All methods showed a power of detecting non-carriers of the mutation (TNR) above 98%, with very small variation with increasing amounts of errors. KNN always attained 100% TNR, except at 12.5% noise. SVML and SVMR had an opposite behaviour: the former showing some false positives when the noise proportion was below 10% and the latter when the proportion was above 5%. RF showed the largest variability of TNR (98.3–100%).

Fig. 3

TPR (sensitivity) and TNR (specificity) as a function of noise in the data. Proportion of true negatives—i.e. non-carriers of the mutation-and true positives—i.e. carriers of the mutation, below—correctly identified (TNR and TPR) by the five methods over 10 error proportions. Lasso black, SVML green, SVMR light blue, RF blue, KNN red

The TPR (sensitivity) for the five methods shows much larger proportions of errors in response to increasing noise in the data. SVML correctly identified all carriers of the mutation only when no errors were introduced in the training set. As the noise proportion increased, TPR approached 0 (i.e. no detection power). A similar trend was shown by RF, which started at TPR = 72% with no mislabelings, and plummeted to TPR = 2.5% when 17.5% mislabeled observations were introduced in the data. For SVMR, TPR ranged between a minimum of 88% and a maximum of 96% when the noise proportion was below 15%, dropping to 55% when the error rate was 17.5%. Lasso showed a very similar TPR pattern as SVMR, with false negatives in the range 0–17% up to 12.5% noise, thereafter jumping to 50% false negatives, and finally relapsing to about 25%. The TPR for KNN was constantly very low (0–19%), irrespective of the amount of noise in the data. Table 1 reports the reciprocals of ACC (total error rate: \(TER=1-ACC\)), TNR (false positive rate: \(FPR=1-TNR\)) and TPR (false negative rate: \(FNR=1-TPR\)).

The total computation time for the five classification methods as a function of error percentage can be seen in Fig. 4. The elapsed time to run 10 times a 10-fold cross-validation scheme ranged from a minimum of 45 min in the scenario with no errors in the training set using KNN to a maximum of 7 h and 24 min using RF with 2.5% errors in the training set. The computation time remained more or less stable for KNN and Lasso over noise thresholds, while it increased approximately linearly with noise both for SVML and SVMR. RF required large computation times at all noise thresholds. Overall, RF was the most computationally demanding algorithm, followed by the two SVM implementations. SVML and SVR took longer than RF only with >15% noise in the data. Only with 20% noise in the data SVML took longer than SVMR to run.

Fig. 4

Computation time as a function of noise in the data. Computation time for the five classification methods and 10 proportions of mislabeled observations. Results from a 10-fold cross-validation procedure, repeated 10 times. Lasso black, SVML green, SVMR light blue, RF blue, KNN red


In this paper, we presented how five classification methods responded to noise in the target variable. We selected two “global” linear methods (Lasso and SVML) and three “local” non-linear methods (KNN, RF and SVMR) in order to explore possible scenarios with state-of-the-art classification methods, each with specific properties.

The overall prediction accuracy in the base scenario (no noise) was close to \(100\%\) with all five classification methods (with KNN providing the lower bound at \(96\%\)). However, the overall accuracy is known to be biased upwards when data are unbalanced [35], as is the ratio between carriers and non-carriers of the TUBD1 mutation (\(4\%\)/\(96\%\)). In such cases, the proportion of errors in the two classes (carriers/non-carriers) gives a better representation of the relative performance of classifiers. In the analysed problem, the true positive and true negative ratios highlight the difficulty of correctly identifying carriers of the mutation (true positives), i.e. of predicting unobserved examples belonging to the minority class. All five methods identified non-carriers with virtually 100% accuracy (TNR = 100%), but they display different behaviour with respect to the prediction of carriers: KNN had a TPR very close to 0%, and never above 20%. With SVML, TPR was 100% with no noise in the data, and then rapidly decreased with increasing errors in the labels, eventually approaching TPR = 0 for noise >12.5%. RF followed a similar pattern, with a starting TPR of 72% that quickly decreased below 25% with minimum around 3%. SVMR and Lasso proved to be relatively more robust to noisy labels in the classification of mutation carriers: their TPR was larger than 80% up to 12.5% noise in the data, and only for larger proportions of errors in the data these two classifiers began to be unreliable. Standard classification algorithms have been shown to perform poorly with unbalanced data, and strategies to deal with unbalancedness have been proposed to improve the prediction accuracy [36, 37].

The worse relative performance of KNN, besides imbalance in the data, can be explained by the difficulty to handle large feature spaces: KNN is known to particularly suffer from the “curse of dimensionality” [38], especially when predictors are collinear, which can well be the case for SNP loci on the same chromosome, likely to be in moderate to high linkage disequilibrium (LD) with each other. The average pairwise LD between SNP loci on BTA19 was estimated as \(r^2=0.126\) [8]. Support vector machines, focussing chiefly on pivotal training observations that define the classification margin (support vectors), are much less affected by high dimensional data. The selection of an appropriate kernel function is however important, since it defines the transformed feature space in which the training set instances will be classified. At base scenario (no noise), SVML outperformed SVMR and had, together with Lasso, a TPR of 100%. This indicates that the decision boundary in this problem is very likely linear, and a method like SVMR—which is known to potentially produce highly non-linear decision boundaries—is expected to perform relatively worse. When errors were introduced in the data, though, the ability of SVMR to accommodate non-linear relationships appeared to be helpful in maintaining relatively high predictive ability.

When introducing incremental percentages of errors in the labels (mislabeled carriers and non-carriers of the mutation), the dataset becomes increasingly noisy, and the task of correctly identifying true carriers and true non-carriers gets more challenging. The overall accuracy decreased as more errors were introduced, but on the whole seemed quite robust to mislabeled observations. This was however true for the accuracy of classifying observations belonging to the majority class (non-carriers), which is trivial with unbalanced data: the TNR remained above 99% irrespective of the amount of noise introduced in the data. On the other hand, the classification of carriers (minority class) gave a very different picture: with TPR suffering much more from noise in the data.

The use of a radial rather than a linear kernel in SVM seemed to make the classification more robust to errors in the labels. With increasing noise, the TPA, TNR and TPR curves became more wiggly, and higher accuracy in the testing rather than the training set was sometimes observed (results not shown: see [8]). When data get noisier, it is more difficult for predicting algorithms to classify observations correctly, as shown also by the increased computation time (Fig. 4); after a certain proportion of errors in the data, predictive models may break down and yield unreliable results (garbage in, garbage out: [39]). In the present dataset, this appeared to happen after 12.5% mislabeled observations in the training data.

Table 1 Total error rate (TER), false positive (FPR) and false negative (FNR) rates for the five classification models over the ten thresholds of random noise introduced in the data

If computation time is also considered, Lasso provided the best combination in terms of classification accuracy and use of computer resources. SVMR showed comparable accuracy, but took much longer at base scenario and, especially, with noise in the data. RF was confirmed to be a demanding algorithm in terms of computing resources (see for instance Nazzicari et al. [40] for imputation of missing genotypes), unless computation strategies like parallelization are adopted; however, RF computation time seemed to be unaffected by noise in the data.

This paper focussed on the different behaviour of some standard machine/statistical learning methods for classification in response to mislabeled observations. When data are noisy, however, active strategies may be adopted to counteract—at least partially—the detrimental effect of noise on results from the statistical analysis: (a) data could be carefully cleaned before analysis [41]; (b) the loss functions by which the predictive equations are optimized can be modified to accommodate errors in the data e.g. by modelling explicitly or implicitly random and non-random errors [12, 42]; (c) locally adaptive approaches may be used to minimize the impact of errors in the data [43, 44].


Machine learning methods have many applications and are gaining increasing popularity also in animal genetics. Data coming from animal recording are not free from errors or inconsistencies. The advent of precision livestock farming and automated data collection can on one hand alleviate the problem of manual or clerical errors, but may on the other hand introduce new sources of noise e.g. random spurious errors, bias in the machine, lack of double checking for errors. When such data are used for predictions, aspects related to the presence of noise have to be taken into account.

This work was the first to assess the impact of phenotyping errors on the accuracy of genomic predictions in animal genetics. The choice of the method used for predictions can influence results, being more or less susceptible to noise. With the present problem of classifying mutation carriers from SNP genotypes, SVM with radial kernel performed relatively well even when the proportion of errors in the data reached 12.5%. Lasso was the second best method, while SVML, RF and KNN were very sensitive to noise (KNN also to data unbalancedness). Taking into account both accuracy and computation time, Lasso provided the best combination among the options considered here (Additional file 1).



single nucleotide polymorphism


K-nearest neighbors


Lasso-penalised logistic regression


support vector machines


SVM with linear kernel


SVM with radial basis function kernel


random forest


total accuracy


true positive rate


true negative rate


  1. 1.

    Lohr S. The age of big data. N Y Times. 2012;11.

  2. 2.

    Domingos P. A few useful things to know about machine learning. Commun ACM. 2012;55(10):78–87.

    Article  Google Scholar 

  3. 3.

    Benevenuto F, Magno G, Rodrigues T, Almeida V. Detecting spammers on twitter. In: Collaboration, electronic messaging, anti-abuse and spam Conference (CEAS), vol. 6; 2010. p. 12.

  4. 4.

    Sebastiani F. Machine learning in automated text categorization. ACM Comput Surv (CSUR). 2002;34(1):1–47.

    Article  Google Scholar 

  5. 5.

    Wang G, Hao J, Ma J, Jiang H. A comparative assessment of ensemble learning for credit scoring. Expert Syst Appl. 2011;38(1):223–30.

    Article  Google Scholar 

  6. 6.

    Sommer C, Gerlich DW. Machine learning in cell biology-teaching computers to recognize phenotypes. J Cell Sci. 2013;126(24):5529–39.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Castelo R, Guigó R. Splice site identification by idlBNs. Bioinformatics. 2004;20(suppl 1):69–76.

    Article  Google Scholar 

  8. 8.

    Biffani S, Dimauro C, Macciotta N, Rossoni A, Stella A, Biscarini F. Predicting haplotype carriers from snp genotypes in bos taurus through linear discriminant analysis. Genet Sel Evol. 2015;47(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Libbrecht MW, Noble WS. Machine learning applications in genetics and genomics. Nat Rev Genet. 2015;16(6):321–32.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Mitchell TM. Machine learning. Burr Ridge: McGraw Hill; 1997.

    Google Scholar 

  11. 11.

    Cesa-Bianchi N, Shalev-Shwartz S, Shamir O. Online learning of noisy data. IEEE Trans Inf Theory. 2011;57(12):7907–31.

    Article  Google Scholar 

  12. 12.

    Natarajan N, Dhillon IS, Ravikumar PK, Tewari A. Learning with noisy labels. In: Advances in neural information processing systems; 2013. p. 1196–204.

  13. 13.

    Bush WS, Moore JH. Genome-wide association studies. PLoS Comput Biol. 2012;8(12):1002822.

    Article  Google Scholar 

  14. 14.

    Meuwissen T, Goddard M. Accurate prediction of genetic values for complex traits by whole-genome resequencing. Genetics. 2010;185(2):623–31.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Kang SJ, Gordon D, Finch SJ. What snp genotyping errors are most costly for genetic association studies? Genet Epidemiol. 2004;26(2):132–41.

    Article  PubMed  Google Scholar 

  16. 16.

    Pompanon F, Bonin A, Bellemain E, Taberlet P. Genotyping errors: causes, consequences and solutions. Nat Rev Genet. 2005;6(11):846–7.

    Article  Google Scholar 

  17. 17.

    Dassonneville R, Brøndum RF, Druet T, Fritz S, Guillaume F, Guldbrandtsen B, Lund MS, Ducrocq V, Su G. Effect of imputing markers from a low-density chip on the reliability of genomic breeding values in holstein populations. J Dairy Sci. 2011;94(7):3679–86.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Poland J, Endelman J, Dawson J, Rutkoski J, Wu S, Manes Y, Dreisigacker S, Crossa J, Sánchez-Villeda H, Sorrells M, et al. Genomic selection in wheat breeding using genotyping-by-sequencing. Plant Genome. 2012;5(3):103–13.

    CAS  Article  Google Scholar 

  19. 19.

    Rutkoski JE, Poland J, Jannink J-L, Sorrells ME. Imputation of unordered markers and the impact on genomic selection accuracy. G3: Genes-Genomes-Genet. 2013;3(3):427–39.

    Article  Google Scholar 

  20. 20.

    Weigel K, de Los Campos G, Vazquez A, Rosa G, Gianola D, Van Tassell C. Accuracy of direct genomic values derived from imputed single nucleotide polymorphism genotypes in jersey cattle. J Dairy Sci. 2010;93(11):5423–35.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Buyske S, Yang G, Matise TC, Gordon D. When a case is not a case: effects of phenotype misclassification on power and sample size requirements for the transmission disequilibrium test with affected child trios. Hum Hered. 2009;67(4):287–92.

    Article  PubMed  Google Scholar 

  22. 22.

    Biscarini F, Nazzicari N, Broccanello C, Stevanato P, Marini S. “Noisy beets”: impact of phenotyping errors on genomic predictions for binary traits in beta vulgaris. Plant Methods. 2016;12(1):36.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Schwarzenbacher H, Burgstaller J, Seefried FR, Wurmser C, Hilbe M, Jung S, Fuerst C, Dinhopl N, Weissenboeck H, Fuerst-Waltl B, et al. A missense mutation in tubd1 is associated with high juvenile mortality in braunvieh and fleckvieh cattle. BMC Genom. 2016;17(1):1.

    Article  Google Scholar 

  24. 24.

    Biscarini F, Schwarzenbacher H, Pausch H, Nicolazzi EL, Pirola Y, Biffani S. Use of snp genotypes to identify carriers of harmful recessive mutations in cattle populations. BMC Genom. 2016;17(1):857.

    Article  Google Scholar 

  25. 25.

    BH2-a haplotype impacting stillbirth.

  26. 26.

    Schwarzenbacher H, Fürst C, Fürst-Waltl B, Dolezal M. A genome-wide search for harmful recessive haplotypes in brown swiss and fleckvieh cattle. In: EAAP (ed.) Book of abstracts of the 63rd EAAP Annual Meeting, Bratislava; 2012. p. 171.

  27. 27.

    Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81(5):1084–97.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Lee SI, Lee H, Abbeel P, Ng AY. Efficient \(\ell\)1 regularized logistic regression. In: Proceedings of the National Conference on artificial intelligence, vol. 21. Menlo Park: AAAI Press; MIT Press; 1999. p. 401.

  29. 29.

    Vapnik VN, Vapnik V. Statistical learning theory, vol. 1. New York: Wiley; 1998.

    Google Scholar 

  30. 30.

    Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.

    Article  Google Scholar 

  31. 31.

    R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.

  32. 32.

    Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Meyer D, Wien FT. Support vector machines. The Interface to libsvm in package; 2015. p. e1071.

  34. 34.

    Kuhn M. Caret package. J Stat Softw. 2008;28(5):1–26.

    Article  Google Scholar 

  35. 35.

    Chawla NV, Japkowicz N, Kotcz A. Editorial: special issue on learning from imbalanced data sets. ACM Sigkdd Explor Newsl. 2004;6(1):1–6.

    Article  Google Scholar 

  36. 36.

    Japkowicz N, et al. Learning from imbalanced data sets: a comparison of various strategies. In: AAAI Workshop on learning from imbalanced data sets, vol. 68, Menlo Park; 2000. p. 10–5.

  37. 37.

    Ertekin S, Huang J, Bottou L, Giles L. Learning on the border: active learning in imbalanced data classification. In: Proceedings of the Sixteenth ACM Conference on conference on information and knowledge management: ACM; 2007. p. 127–36.

  38. 38.

    James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning. New York: Springer; 2013.

    Google Scholar 

  39. 39.

    Tweedie R, Mengersen K, Eccleston J. Garbage in, garbage out: can statisticians quantify the effects of poor data. Chance. 1994;7(2):20–7.

    Google Scholar 

  40. 40.

    Nazzicari N, Biscarini F, Cozzi P, Brummer EC, Annicchiarico P. Marker imputation efficiency for genotyping-by-sequencing data in rice (Oryza sativa) and alfalfa (Medicago sativa). Mol Breed. 2016;36(6):1–16.

    CAS  Article  Google Scholar 

  41. 41.

    Rahm E, Do HH. Data cleaning: problems and current approaches. IEEE Data Eng Bull. 2000;23(4):3–13.

    Google Scholar 

  42. 42.

    Mnih V, Hinton GE. Learning to label aerial images from noisy data. In: Proceedings of the 29th International conference on machine learning (ICML-12); 2012. p. 567–74.

  43. 43.

    Wettschereck, D., Dietterich, T.G.: Locally adaptive nearest neighbor algorithms. Adv Neural Inf Proc Syst. 1994: 184.

  44. 44.

    Chu F, Wang Y, Zaniolo C. An adaptive learning approach for noisy data streams. In: Fourth IEEE International Conference on data mining, ICDM’04: IEEE; 2004. p. 351–54.

Download references

Authors' contributions

SB and FB analysed the data and prepared the manuscript. HS and HP provided the data and contributed to the conception of the study. All authors read and approved the final manuscript.


Animal samples and genotypes were provided by ZuchtData GmbH. The research leading to these results has received funding from the Italian national research project “GenHome”. FB was in part financially supported by the EU-FP7 project “FECUND” (No. 312097).

Competing interests

The authors declare that they have no competing interests.

Availability of supporting data

The data used in this study are a subset of the data from the recently published paper (“A missense mutation in TUBD1 is associated with high juvenile mortality in Braunvieh and Fleckvieh cattle”) by Schwarzenbacher et al. [23], BMC Genomics. Whole-genome sequencing data of BH2hom were deposited in the European Nucleotide Archive under accession number PRJEB12807 ( The sequencing data of 241 animals from the 1000 bull genomes project are publically available in the Sequence Read Archive of NCBI ( under Accession Number SRP039339.

Ethics approval and consent to participate

All data for the present study were provided by ZuchtData EDV-Dienstleistungen GmbH (Austria; No ethical approval was required for this study. The data used in this study are a subset of the data recently published by Schwarzenbacher et al. [23] in BMC Genomics.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Filippo Biscarini.

Additional file

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Biffani, S., Pausch, H., Schwarzenbacher, H. et al. The effect of mislabeled phenotypic status on the identification of mutation-carriers from SNP genotypes in dairy cattle. BMC Res Notes 10, 230 (2017).

Download citation


  • Noisy data
  • Genomic predictions
  • SNP genotypes
  • Bos taurus
  • Animal genomics
  • Classification
  • Phenotyping errors