Associations between the breeding values of Holstein-Friesian bulls and longevity and culling reasons of their daughters

Taking into account functional traits in the breeding practice should lead to a longer productive life of cows. However, despite the increased contribution of these traits in bull selection indices, their daughters are frequently culled as early as the 2nd or 3rd lactation. The problem is whether and to what extent the genetic potential of animals is realized in the production practice. Therefore, the purpose of this study was to determine the associations between the breeding value (BV) of bulls and their daughters for cow longevity and culling reasons in the Holstein-Friesian cattle population in Poland. Data for 532 062 cows culled in 2012, 2015, and 2018 were analyzed. A majority of 5 045 cow sires originated from Poland, Germany, France, the Netherlands, and the United States. The highest variation in the contribution of culling reasons was for the cows culled at the age of 2–4 years. The contribution of the culling reasons, analyzed in relation to the cow culling age, remained similar and the only exception was culling because of old age, for which a significant increase was observed only for the culling age of at least 9 years (13.8%), which was reached by only 7.3% of the cows. The sires were characterized by generally high BV for conformation and reproductive traits. However, they had, at most, the average genetic potential for functional longevity. There were a number of beneficial associations found between the BV of bulls and the distribution of culling reasons in their daughters. For example, it concerns relations between the somatic cell score in milk and culling due to udder diseases and low milk yield, between the interval from calving to first insemination and low milk yield, between the protein yield and old age, or between the BV for certain conformation traits (size, udder) and cow culling due to age. In these cases, as the BV increased for a given trait, the contribution of the corresponding cow culling reason tended to decrease. Our study showed that it seems reasonable to consider Holstein-Friesian cows aged at least 9 years at culling to be long-living animals. This is primarily evidenced by the rapid increase in the culling due to old age in relation to younger cows. Nowadays the above age limit can be suggested as a criterion of longevity for Holstein-Friesian cows but the criterion should be updated to the relation genotype-environment-economy that tends to change over time.     

Converting multi-trait breeding objectives into operative selection indexes to ensure genetic gains in low-input sheep and goat breeding programmes

Small ruminant breeding programmes in low-input production systems are best organised at the community level. Participant farmers have to agree on goal traits and their relative importance. When BLUP breeding values of goal traits are not available in time, appropriate selection indexes can be used to aid visual selection. Taking Ethiopian Abergelle goat and Bonga sheep community-based breeding programmes (CBBPs) as an example, breeding objective functions were defined and selection indexes were constructed and evaluated. Breeding goals for Abergelle goats included early sale weight, survival and milk production. Breeding goals for Bonga included the number of offspring born, sale weight and survival. Economic weights of objective traits can be used in several ways depending on measured traits and the reliability of their genetic parameters. Selection indexes included combinations of objective traits measured on candidates and their dams and situations when Abergelle communities prefer to restrict genetic changes in number of offspring born or adult weight and when Bonga communities prefer to restrict changes in adult weight. Genetic and economic gains were evaluated as well as sensitivity to feed cost assumptions and to repeated dam records. After independent culling on preponderant traits such as coat colour and horn/tail type, sires in Abergelle goat community breeding programmes should be selected on indexes including at least own early live weight and their dams average milk production records. Sires for Bonga sheep programmes should be selected on own early live weight and desirably also on their dam’s number of offspring born. Sensitivity to feed cost assumptions was negligible but repeated measurements of dam records improved index accuracies considerably. Restricting genetic changes in number of offspring born or adult weight is not recommended.     

A genetic model involving fetal effects on traits of the dam

The genetic fetal effects model shows that the usual sire effect is composed of one-half the direct additive genetic value and one-fourth of the fetal additive genetic value of the sire. The usual sire component of variance is actually the variance of that function. Genetic covariances between records of relatives influenced by fetuses of related sires can easily be written. If the magnitude of fetal sire effects is such that nonrandom use of fetal sires on daughters of sires being evaluated on daughter performance results in bias, the bias can be eliminated (Henderson 1975) by considering the fetal sire effects to be fixed effects. Some reduction in prediction error variance is likely by including fetal sire in the sire evaluation model.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Estimation by simulation of the efficiency of the French marker-assisted selection program in dairy cattle (Open Access publication)

The efficiency of the French marker-assisted selection (MAS) was estimated by a simulation study. The data files of two different time periods were used: April 2004 and 2006. The simulation method used the structure of the existing French MAS: same pedigree, same marker genotypes and same animals with records. The program simulated breeding values and new records based on this existing structure and knowledge on the QTL used in MAS (variance and frequency). Reliabilities of genetic values of young animals (less than one year old) obtained with and without marker information were compared to assess the efficiency of MAS for evaluation of milk, fat and protein yields and fat and protein contents. Mean gains of reliability ranged from 0.015 to 0.094 and from 0.038 to 0.114 in 2004 and 2006, respectively. The larger number of animals genotyped and the use of a new set of genetic markers can explain the improvement of MAS reliability from 2004 to 2006. This improvement was also observed by analysis of information content for young candidates. The gain of MAS reliability with respect to classical selection was larger for sons of sires with genotyped progeny daughters with records. Finally, it was shown that when superiority of MAS over classical selection was estimated with daughter yield deviations obtained after progeny test instead of true breeding values, the gain was underestimated.     

Impact of bovine somatotropin on ranking for genetic value of dairy sires for milk yield traits and somatic cell score

Records of Holstein cows were used to examine how different models account for the effect of bovine somatotropin (bST) treatment on genetic evaluation of dairy sires for yield traits and somatic cell score. Data set 1 included 65,720 first-lactation records. Set 2 included 50,644 second-lactation records. Set 3 included 45,505 records for lactations three, four and five. Estimated breeding values (EBV) of sires were with three different animal models. With Model 1, bST administration was ignored. With Model 2, bST administration was used as a fixed effect. With Model 3, administration of bST was used to define the contemporary group (herd-year-month of calving-bST). Correlations for EBV of 1,366 sires with treated daughters between pairs of the three models were calculated for milk, fat and protein yields and somatic cell score for the three data sets. Correlations for EBV of sires between pairs of models for all traits ranged from 0.971 to 0.999. Fractions of sires with bST-treated progeny selected in common (top 10 to 15%) were 0.94 and usually greater for all pairs of models for all traits and parities. For this study, the method of statistical adjustment for bST treatment resulted in a negligible effect on genetic evaluations of sires when some daughters were treated with bST and suggests that selection of sires to produce the next generation of sires and cows might not be significantly affected by how the effect of bST is modeled for prediction of breeding values for milk, fat and protein yields and somatic cell score.     

Reproductive fitness and artificial selection in animal breeding: culling on fitness prevents a decline in reproductive fitness in lines of Drosophila melanogaster selected for increased inebriation time

Summary The maintenance of reproductive fitness in lines subjected to artificial selection is one of the major problems in animal breeding. The decline in reproductive performance has neither been predictable from heritabilities and genetic correlations, nor have conventional selection indices been adequate to avoid the problem. Gowe (1983) has suggested that the heritabilities of reproductive traits are non-linear, with heritabilities being higher on the lower fitness side. Consequently, he has predicted that culling on reproductive fitness in artificial selection lines will be effective in preventing the usual declines in fitness. An experimental evaluation of Gowe's prediction has been carried out by comparing fitnesses of replicated lines of three treatments: selection for increased inebriation time without culling on fitness (HO), selection for inebriation time with culling of 20% (4/20) of selected females on reproductive fitness (HS), and unselected controls (C). Response to selection for inebriation time in the two selection treatments was similar. After 25 generations, the competitive index, a measure of reproductive fitness, was significantly lower in the HO treatment than the HS treatment, while the HS treatment did not differ from the control lines or the base population. These results demonstrate for the first time that culling on reproductive fitness in selection lines can be used to prevent the usual decline in reproductive performance.     

The value of cows in reference populations for genomic selection of new functional traits

Today, almost all reference populations consist of progeny tested bulls. However, older progeny tested bulls do not have reliable estimated breeding values (EBV) for new traits. Thus, to be able to select for these new traits, it is necessary to build a reference population. We used a deterministic prediction model to test the hypothesis that the value of cows in reference populations depends on the availability of phenotypic records. To test the hypothesis, we investigated different strategies of building a reference population for a new functional trait over a 10-year period. The trait was either recorded on a large scale (30 000 cows per year) or on a small scale (2000 cows per year). For large-scale recording, we compared four scenarios where the reference population consisted of 30 sires; 30 sires and 170 test bulls; 30 sires and 2000 cows; or 30 sires, 2000 cows and 170 test bulls in the first year with measurements of the new functional trait. In addition to varying the make-up of the reference population, we also varied the heritability of the trait (h2 = 0.05 v. 0.15). The results showed that a reference population of test bulls, cows and sires results in the highest accuracy of the direct genomic values (DGV) for a new functional trait, regardless of its heritability. For small-scale recording, we compared two scenarios where the reference population consisted of the 2000 cows with phenotypic records or the 30 sires of these cows in the first year with measurements of the new functional trait. The results showed that a reference population of cows results in the highest accuracy of the DGV whether the heritability is 0.05 or 0.15, because variation is lost when phenotypic data on cows are summarized in EBV of their sires. The main conclusions from this study are: (i) the fewer phenotypic records, the larger effect of including cows in the reference population; (ii) for small-scale recording, the accuracy of the DGV will continue to increase for several years, whereas the increases in the accuracy of the DGV quickly decrease with large-scale recording; (iii) it is possible to achieve accuracies of the DGV that enable selection for new functional traits recorded on a large scale within 3 years from commencement of recording; and (iv) a higher heritability benefits a reference population of cows more than a reference population of bulls.     

Genetic parameters for feed efficiency in Romane rams and responses to single-generation selection

Feeding costs represent one of the highest expenditures in animal production systems. Breeding efficient animals that express their growth potential while eating less is therefore a major objective for breeders. We estimated the genetic parameters for feed intake, feed efficiency traits (residual feed intake (RFI) and feed conversion ratio (FCR)), growth and body composition traits in the Romane meat sheep breed. In these traits, selection responses to single-generation divergent selection on RFI were evaluated. From 2009 to 2016, a total of 951 male lambs were tested for 8 weeks starting from 3 months of age. They were weighed at the beginning and at the end of the testing period. Backfat thickness and muscle depth were recorded at the end of the testing period through ultrasound measurements. Feed intake was continuously recorded over the testing period using the automatic concentrate feeders. The heritability of RFI was estimated at 0.45 ± 0.08, which was higher than the heritability of FCR (0.30 ± 0.08). No significant genetic correlations were observed between RFI and growth traits. A favourable low negative genetic correlation was estimated between RFI and muscle depth (−0.30 ± 0.15), though additional data are needed to confirm these results. The selection of low RFI sires based on their breeding values led to the production of lambs eating significantly less concentrate (3% decrease in the average daily feed intake), but with the same growth as lambs from sires selected based on high RFI breeding values. We concluded that in meat sheep, RFI is a heritable trait that is genetically independent of post-weaning growth and body composition traits. A one-generation divergent selection based on RFI breeding values highlighted that substantial gains in feeding costs can be expected in selection schemes for meat sheep breeds.     

Simulation study on the efficiencies of MOET nucleus breeding schemes applying marker assisted selection in dairy cattle

Advantages of breeding schemes using genetic marker information and/or multiple ovulation and embryo transfer (MOET) technology over the traditional approach were extensively evaluated through simulation. Milk yield was the trait of interest and QTL was the genetic marker utilized. Eight dairy cattle breeding scenarios were considered, i.e., traditional progeny testing breeding scheme (denoted as STANPT), GASPT scheme including a pre-selection of young bulls entering progeny testing based on their own QTL information, MOETPT scheme using MOET technology to generate young bulls and a selection of young bulls limited within the full-sib family, GAMOPT scheme adopting both QTL pre-selection and MOET technology, COMBPT scheme using a mixed linear model which considered QTL genotype instead of the BLUP model in GAMOPT, and three non-progeny testing schemes, i.e. the MOET, GAMO and COMB schemes, corresponding to MOETPT, GAMOPT and COMBPT with progeny testing being part of the system. Animals were selected based on their breeding value which was estimated under an animal model framework. Sequential selection over 17 years was performed in the simulations and 30 replicates were designed for each scenario. The influences of using QTL information and MOET technology on favorable QTL allele frequency, true breeding values, polygenetic breeding values and the accumulated genetic superiority were extensively evaluated, for five different populations including active sires, lactating cows, bull dams, bull sires, and young bulls. The results showed that the combined schemes significantly outperformed other approaches wherein accumulated true breeding value progressed. The difference between schemes exclusively using QTL information or MOET technology was not significant. The STANPT scheme was the least efficient among the 8 schemes. The schemes using MOET technology had a higher polygenetic response than others in the 17th year. The increases of frequency of the favorable QTL allele varied more greatly across the 3 male groups than in the lactating cows group. The accumulated genetic superiorities of the GASPT scheme, MOETPT scheme, GAMOPT scheme, COMBPT scheme, MOET scheme, GAMO scheme and COMB scheme over the STANPT scheme were 8.42%, 3.59%, 14.58%, 18.54%, 4.12%, 14.12%, 16.50% in active sires and 2.70%, 5.00%, 11.05%, 12.78%, 7.51%, 17.12%, 25.38% in lactating cows. Supported by Key Project for Introducing Advanced International Agriculture Science & Technologies (Grant No. 2006-G48), the National Key Basic Research and Development Program of China (Grant No. 2006CB102107) and National Key Technology Research and Development Program of China (Grant No. 2006BAD04A01).     

An approximate method for optimum independent culling level selection for n stages of selection with explicit solutions

Summary An approximate method with explicit solutions to apply independent culling levels for multiple traits in n-stages of selection was developed. An approximate solution was found for sequentially selected traits. Two assumptions were necessary. The first was to assume that subsequent selection would not appreciably change the mean of traits already selected, and the second was to approximate the variance of a correlated trait in a selected population with an upward biased projection. The procedure was shown to give near optimal results regardless of selection intensity or genetic correlations if phenotypic correlations among traits were low. The procedure gave poor results only for certain sequences of selection when phenotypic correlations were high. However, in those cases good results were obtained using a different sequence of selection. With high correlations, the procedure is recommended only after comparing solutions and expected genetic gain for all sequences of selection. If the expected aggregate gain for the sequence of selection desired is less than that of another order, culling points associated with the optimal ordering must be determined. Genetic gain from use of culling points is independent of order of selection. The procedure is recommended for use with computer programs that attempt to find optimal culling points to reduce computational time and to check results.Journal Paper No. 12448 of the Purdue University Agricultural Experiment Station     

Turning science on robust cattle into improved genetic selection decisions

More robust cattle have the potential to increase farm profitability, improve animal welfare, reduce the contribution of ruminant livestock to greenhouse gas emissions and decrease the risk of food shortages in the face of increased variability in the farm environment. Breeding is a powerful tool for changing the robustness of cattle; however, insufficient recording of breeding goal traits and selection of animals at younger ages tend to favour genetic change in productivity traits relative to robustness traits. This paper has extended a previously proposed theory of artificial evolution to demonstrate, using deterministic simulation, how choice of breeding scheme design can be used as a tool to manipulate the direction of genetic progress, whereas the breeding goal remains focussed on the factors motivating individual farm decision makers. Particular focus was placed on the transition from progeny testing or mass selection to genomic selection breeding strategies. Transition to genomic selection from a breeding strategy where candidates are selected before records from progeny being available was shown to be highly likely to favour genetic progress in robustness traits relative to productivity traits. This was shown even with modest numbers of animals available for training and when heritability for robustness traits was only slightly lower than that for productivity traits. When transitioning from progeny testing to a genomic selection strategy without progeny testing, it was shown that there is a significant risk that robustness traits could become less influential in selection relative to productivity traits. Augmentations of training populations using genotyped cows and support for industry-wide improvements in phenotypic recording of robustness traits were put forward as investment opportunities for stakeholders wishing to facilitate the application of science on robust cattle into improved genetic selection schemes.     

Studies on changes of estimated breeding values of U.S. Holstein bulls for final score from the first to second crop of daughters

The purpose of this study was to find ways of reducing changes of sire predicted transmitting ability for type’s final scores (PTATs) from the first to second crop of daughters. The PTATs were estimated from two datasets: D01 (scores recorded up to 2001) and D05 (scores recorded up to 2005). The PTAT changes were calculated as the difference between the evaluations based on D01 and D05. The PTATs were adjusted to a common genetic base of all evaluated cows born in 1995. The single-trait (ST) animal model included the fixed effects of the herd–year–season–classifier, age by year group at classification, stage of lactation at classification, registry status of animals, and additive genetic and permanent environment random effects. Unknown parent groups (UPGs) were defined based on every other birth year starting from 1972. Modifications to the ST model included the usage of a single record per cow, separate UPGs for first and second crop daughters, separate UPGs for sires and dams, and deepened pedigrees for dams with missing phenotypic records. Also, the multiple-trait (MT) model treated records of registered and grade cows as correlated traits. The mean PTAT change, for all of the sires, was close to zero in all of the models analyzed. The estimated mean PTAT change for 145 sires with 40 to 100 first crop and ≥200 second crop daughters was −0.33, −0.20, −0.13, −0.28, and −0.12 with ST, only first records, only last records, updated pedigrees, and allowing separate parent groups (PGs) for sires and dams after updating the pedigrees, respectively. The percentages of sires showing PTAT decline were reduced from 74.5 (with ST) to 57.3 by using only the last records of cows, and to 56.4 by allowing separate UPGs for sires and dams after updating the pedigrees. Though updating of the pedigrees alone was not effective, separate UPGs for sires together with additional pedigree was helpful in reducing the bias.     

Genetic correlation estimations between artificial insemination sire performances and their progeny beef traits both measured in test stations

In France, beef traits of artificial insemination (AI) beef bulls are improved through the sequential selection for their own performances and for their male progeny performances, both being recorded in test stations. The efficiency of such programmes mainly depends on the genetic correlations between sire performances and progeny beef traits. Such correlations were independently estimated, using the multivariate REML (restricted maximum likelihood) method in a Limousin and a Charolais programme. In both breeds, high genetic correlations were observed between sires and progeny analogous morphology scores (from 0.64 to 0.82). Genetic correlations estimated between sires and progeny growth (from 0.41 to 0.70) were lower probably due to the difference of diet in central and progeny stations. Correlations between sire muscling scores and progeny skeletal frames (and vice-versa) were negative (from -0.05 to -0.58). The genetic correlations of sire traits with progeny dressing percentage (DP_p) and carcass fatness score (CF_p) were only low to moderate. These results show that the selection of bulls at the end of performance testing in test stations may be efficient in improving progeny growth and morphology. However, such a selection is insufficient in improving their dressing percentage and carcass composition.     

Physical examination of the reproductive organs of range beef bulls in Mozambique

Physical examinations were performed on 3991 bulls reared in subtropical (75%) and tropical (25%) Mozambique. A culling rate of 16.3% (651 sires) was found. The range of culling rate between farms varied from 3 to 44%. The main culling reasons were a) epididymitis, b) atrophy and/or hypoplasia, and c) testicular fibrosis (contributing to 30.4, 20.9, and 15.2%, respectively, of the culled bulls). Epididymitis-vaginitis (epivag) syndrome was considered the main reason for the high incidence of epididymitis and testicular fibrosis. A highly significant difference (P<0.001) in culling rate for Simmental (37.7%) and Brahman (13.1%) bulls was found. A more comprehensive investigation of culling rate of different breeds used in the country is needed, as well as a program for controlling epivag in the more affected areas. Examining the reproductive organs of the beef bulls in Mozambique before the breeding season is very important to improve fertility in the beef herds.     

Genetic improvement of laying hens viability using survival analysis

The survival of about eight generations of a large strain of laying hens was analysed separating the rearing period (RP) from the production period (PP), after hens were housed. For RP (respectively PP), 97.8% (resp., 94.1% ) of the 109 160 (resp., 100 665) female records were censored after 106 days (resp., 313 days) on the average. A Cox proportional hazards model stratified by flock (= season) and including a hatch-within-flock (HWF) fixed effect seemed to reasonably fit the RP data. For PP, this model could be further simplified to a non-stratified Weibull model. The extension of these models to sire-dam frailty (mixed) models permitted the estimation of the sire genetic variances at 0.261 ± 0.026 and 0.088 ± 0.010 for RP and PP, respectively. Heritabilities on the log scale were equal to 0.48 and 0.19. Non-additive genetic effects could not be detected. Selection was simulated by evaluating all sires and dams, after excluding all records from the last generation. Then, actual parents of this last generation were distributed into four groups according to their own pedigree index. Raw survivor curves of the progeny of extreme parental groups substantially differed (e.g., by 1.7% at 300 days for PP), suggesting that selection based on solutions from the frailty models could be efficient, despite the very large proportion of censored records.     

Genetic evaluation of ovulatory disorders in Austrian Fleckvieh cows: a comparison between linear models and survival analysis

The objective of this study was to compare linear models and survival analysis for genetic evaluation of ovulatory disorders, which included veterinary treatments of silent heat/anestrus and cystic ovaries. Data of 23 450 daughters of 274 Austrian Fleckvieh sires were analyzed. For linear model analyses, ovulatory disorders were defined as a binary response (presence or absence) in the time periods from calving to 150 days after calving and from calving to 300 days after calving. For survival analysis, ovulatory disorders were defined either as the number of days from calving to the day of the first treatment for an ovulatory disorder (uncensored record) or from calving to the day of culling, or the last day of the period under investigation (until 150 or 300 days after calving; censored record). Estimates of heritability were very similar (0.016 to 0.020) across methods and periods. Correlations between sire estimated breeding value from linear model and survival analysis were 0.98, whereas correlations between different time periods were somewhat lower (0.95 and 0.96). The results showed that the length of time period had a larger effect on genetic evaluation than methodology.     

An Epidemiological and Genetic Study on Registered Diseases in Finnish Ayrshire Cattle

The epidemiology and genetic variability of the most common dairy cow diseases were examined. This paper describes the data set, lactation incidence rates and culling during lactation. The data set consisted of the lactation records of 73,368 Finnish Ayrshire dairy cows. Each cow was under observation for 2 days before and 305 days after calving. Lactational incidence rates (%) for the most common diseases were: ovulatory dysfunction 7.0, ketosis 6.0, acute mastitis 5.4, an oestrus and suboestrus 5.2, retained placenta 4.5, parturient paresis 3.8 and teat injury 2.6. Multiple logistic regression was utilized to investigate the possible effects of certain factors on culling. The model predicted the log odds for culling as an additive function of the explanatory factors. Using the estimated odds and forming the odds ratios it was possible to investigate, relative risks between any combination of groups of the explanatory factors. The risk of culling increased with parity after the second parturition, and with increasing herd milk yield. Mastitis and parturient paresis had positive associations with culling, while ketosis and infertility had negative associations. Heritability estimates for culling in various parity groups were from 2 % to 9 % on the binomial scale corresponding from 5 % to 14 % on the normal scale. There was a neagtive genetic correlation between culling and previous milk production.     

The Wilhelmine E. Key 1986 invitational lecture. Fifty years of genetic load

The author's involvement with and his successive reactions to the genetic load concept [whose beginning is identified with Haldane's paper. The effect of variation on fitness] is presented in the form of a personal odyssey. The major change in attitude involved the realization that the density- and frequency-independent selection discussed by most population geneticists has little bearing on events transpiring within natural populations; instead, natural selection should be viewed primarily as a density- and frequency-dependent phenomenon. Under this view, the culling of a large number of young zygotes to the considerably smaller number of adults that can be sustained by the environment is an essential process enabling any population's continued existence; to the extent that genetic variation facilitates culling, a genetic load (in direct opposition to the early view) can enhance a population's persistence through time.     

The importance of culling in Johne's disease control

Johne's disease is caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection and results in economic losses in the dairy industry. To control MAP transmission in herds, test-based culling has been recommended and immediate culling of high shedding animals is typically implemented. In this study, we quantified the effects of MAP control in US dairy herds, using the basic reproduction ratio R₀. The effectiveness of culling strategies was evaluated for good and poor herd management (low- and high-transmission rates, respectively) by a phase diagram approach. To establish a quantitative relationship between culling rates and test properties, we defined the average detection times for low and high shedding animals. The effects of various culling strategies and test characteristics, such as test sensitivity, test turnaround time, and testing interval, were analyzed. To understand the overall effect of model parameters on R₀, we performed global uncertainty and sensitivity analyses. We also evaluated the effectiveness of culling only high shedding animals by comparing three test methods (fecal culture, fecal polymerase chain reaction, PCR, and enzyme-linked immunosorbent assay, ELISA). Our study shows that, in the case of good herd management, culling of only high shedding animals may be effective in controlling MAP transmission. However, in the case of poor management, in addition to immediate culling of high shedding animals, culling of low shedding animals (based on the fecal culture test) will be necessary. Culling of low shedding animals may be delayed 6-12 months, however, if a shorter testing interval is applied. This study suggests that if farmers prefer culling only high shedding animals, faster MAP detection tests (such as the fecal PCR and ELISA) of higher sensitivity should be applied with high testing frequency, particularly on farms with poor management. Culling of infectious animals with a longer testing interval is generally not effective to control MAP.