Search

It is of interest to evaluate crossbred pigs for hot carcass weight (HCW) and birth weight (BW); however, obtaining a HCW record is dependent on livability (LIV) and retained tag (RT). The purpose of this study is to analyze how HCW evaluations are affected when herd removal and missing identification are included in the model and examine if accounting for the reasons for missing traits improves the accuracy of predicting breeding values. Pedigree information was available for 1,965,077 purebred and crossbred animals. Records for 503,716 commercial three-way crossbred terminal animals from 2014 to 2019 were provided by Smithfield Premium Genetics. Two pedigree-based models were compared; model 1 (M1) was a threshold-linear model with all four traits (BW, HCW, RT, and LIV), and model 2 (M2) was a linear model including only BW and HCW. The fixed effects used in the model were contemporary group, sex, age at harvest (for HCW only), and dam parity. The random effects included direct additive genetic and random litter effects. Accuracy, dispersion, bias, and Pearson correlations were estimated using the linear regression method. The heritabilities were 0.11, 0.07, 0.02, and 0.04 for BW, HCW, RT, and LIV, respectively, with standard errors less than 0.01. No difference was observed in heritabilities or accuracies for BW and HCW between M1 and M2. Accuracies were 0.33, 0.37, 0.19, and 0.23 for BW, HCW, RT, and LIV, respectively. The genetic correlation between BW and RT was 0.34 ± 0.03, and between BW and LIV was 0.56 ± 0.03. Similarly, the genetic correlation between HCW and RT was 0.26 ± 0.04, and between HCW and LIV was 0.09 ± 0.05, respectively. The positive and moderate genetic correlations between BW and other traits imply a heavier BW resulted in a higher probability of surviving to harvest. Genetic correlations between HCW and other traits were lower due to the large quantity of missing records. Despite the heritable and correlated aspects of RT and LIV, results imply no major differences between M1 and M2; hence, it is unnecessary to include these traits in classical models for BW and HCW.

Inbreeding depression is a growing concern in livestock because it can detrimentally affect animal fitness, health, and production levels. Genomic information can be used to more effectively capture variance in Mendelian sampling, thereby enabling more accurate estimation of inbreeding, but further progress is still required. The calculation of inbreeding for herd management purposes is largely still done using pedigree information only, although inbreeding coefficients calculated in this manner have been shown to be less accurate than genomic inbreeding measures. Continuous stretches of homozygous genotypes, so called runs of homozygosity, have been shown to provide a better estimate of autozygosity at the genomic level than conventional measures based on inbreeding coefficients calculated through conventional pedigree information or even genomic relationship matrices. For improved and targeted management of genomic inbreeding at the population level, the development of methods that incorporate genomic information in mate selection programs may provide a more precise tool for reducing the detrimental effects of inbreeding in dairy herds. Additionally, a better understanding of the genomic architecture of inbreeding and incorporating that knowledge into breeding programs could significantly refine current practices. Opportunities to maintain high levels of genetic progress in traits of interest while managing homozygosity and sustaining acceptable levels of heterozygosity in highly selected dairy populations exist and should be examined more closely for continued sustainability of both the dairy cattle population as well as the dairy industry. The inclusion of precise genomic measures of inbreeding, such as runs of homozygosity, inbreeding, and mating programs, may provide a path forward. In this symposium review article, we describe traditional measures of inbreeding and the recent developments made toward more precise measures of homozygosity using genomic information. The effects of homozygosity resulting from inbreeding on phenotypes, the identification and mapping of detrimental homozygosity haplotypes, management of inbreeding with genomic data, and areas in need of further research are discussed.

Genomic information has a limited dimensionality (number of independent chromosome segments [Me]) related to the effective population size. Under the additive model, the persistence of genomic accuracies over generations should be high when the nongenomic information (pedigree and phenotypes) is equivalent to Me animals with high accuracy. The objective of this study was to evaluate the decay in accuracy over time and to compare the magnitude of decay with varying quantities of data and with traits of low and moderate heritability. The dataset included 161,897 phenotypic records for a growth trait (GT) and 27,669 phenotypic records for a fitness trait (FT) related to prolificacy in a population with dimensionality around 5,000. The pedigree included 404,979 animals from 2008 to 2020, of which 55,118 were genotyped. Two single-trait models were used with all ancestral data and sliding subsets of 3-, 2-, and 1-generation intervals. Single-step genomic best linear unbiased prediction (ssGBLUP) was used to compute genomic estimated breeding values (GEBV). Estimated accuracies were calculated by the linear regression (LR) method. The validation population consisted of single generations succeeding the training population and continued forward for all generations available. The average accuracy for the first generation after training with all ancestral data was 0.69 and 0.46 for GT and FT, respectively. The average decay in accuracy from the first generation after training to generation 9 was −0.13 and −0.19 for GT and FT, respectively. The persistence of accuracy improves with more data. Old data have a limited impact on the predictions for young animals for a trait with a large amount of information but a bigger impact for a trait with less information.

The science around the use of masks by the general public to impede COVID-19 transmission is advancing rapidly. Policymakers need guidance on how masks should be used by the general population to combat the COVID-19 pandemic. In this narrative review, we develop an analytical framework to examine mask usage, considering and synthesizing the relevant literature to inform multiple areas: population impact; transmission characteristics; source control; PPE; sociological considerations; and implementation considerations. A primary route of transmission of COVID-19 is via respiratory droplets, and is known to be transmissible from presymptomatic and asymptomatic individuals. Reducing disease spread requires two things: first, limit contacts of infected individuals via physical distancing and other measures, and second, reduce the transmission probability per contact. The preponderance of evidence indicates that mask wearing reduces the transmissibility per contact by reducing transmission of infected droplets in both laboratory and clinical contexts. Public mask wearing is most effective at reducing spread of the virus when compliance is high. The decreased transmissibility could substantially reduce the death toll and economic impact while the cost of the intervention is low. Given the current shortages of medical masks we recommend the adoption of public cloth mask wearing, as an effective form of source control, in conjunction with existing hygiene, distancing, and contact tracing strategies. Because many respiratory droplets become smaller due to evaporation, we recommend increasing focus on a previously overlooked aspect of mask usage: mask-wearing by infectious people ("source control") with benefits at the population-level, rather than mask-wearing by susceptible people, such as health-care workers, with focus on individual outcomes. We recommend that public officials and governments strongly encourage the use of widespread face masks in public, including the use of appropriate regulation.

Across the majority livestock species, routinely collected genomic and pedigree information has been incorporated into evaluations using single-step methods. As a result, strategies that reduce genotyping costs without reducing the response to selection are important as they could have substantial economic impacts on breeding programs. Therefore, the objective of the current study was to investigate the impact of selectively genotyping selection candidates on the selection response using simulation. Populations were simulated to mimic the genome and population structure of a swine and cattle population undergoing selection on an index comprised of the estimated breeding values (EBV) for 2 genetically correlated quantitative traits. Ten generations were generated and genotyping began generation 7. Two phenotyping scenarios were simulated that assumed the first trait was recorded early in life on all individuals and the second trait was recorded on all versus a random subset of the individuals. The EBV were generated from a bivariate animal model. Multiple genotyping scenarios were generated that ranged from not genotyping any selection candidates, a proportion of the selection candidates based on either their index value or chosen at random, and genotyping all selection candidates. An interim index value was utilized to decide who to genotype for the selective genotype strategy. The interim value assumed only the first trait was observed and the only genotypic information available was on animals in previous generations. Within each genotyping scenario 25 replicates were generated. Within each genotyping scenario the mean response per generation and the degree to which EBV were inflated/deflated was calculated. Across both species and phenotyping strategies, the plateau of diminishing returns was observed when 60% of the selection candidates with the largest index values were genotyped. When randomly genotyping selection candidates, either 80 or 100% of the selection candidates needed to be genotyped for there not to be a reduction in the index response. Across both populations, no differences in the degree that EBV were inflated/deflated for either trait 1 or 2 were observed between nongenotyped and genotyped animals. The current study has shown that animals can be selectively genotyped in order to optimize the response to selection as a function of the cost to conduct a breeding program using single-step genomic best linear unbiased prediction.

Vaginal and rectal prolapse incidence has been increasing in the U.S. commercial breeding sow population in recent years. This has resulted in increased sow breeding mortality rates in the U.S. due to both increased mortality and euthanasia. The objective of the present study was to estimate the potential role that genetics plays in the prolapse incidence of the U.S. commercial sow population. Data from 11,481 purebred Large White sows from two commercial farms were used in this study. Genetic parameters for prolapse incidence were estimated based on linear and threshold animal models utilizing AIREMLF90 and THRGIBBS1F90, respectively. Across both models, the fixed effects included contemporary group as year-month of removal and removal parity and a random effect was sow. A pedigree-based relationship matrix was included in the model. The highest prolapse incidence occurred in the 1st when compared to other parities. Heritability estimates for combined vaginal and rectal prolapse incidence were 0.03 ± 0.01 and 0.003 ± 0.003 for the linear and threshold models, respectively. The results from this study indicate that there is very little if any genetic impact on the prolapse incidence. Pork producers should focus on the environmental effects contributing to reducing prolapse incidence.

Traditionally, pedigree-based relationship coefficients have been used to manage the inbreeding and degree of inbreeding depression that exists within a population. The widespread incorporation of genomic information in dairy cattle genetic evaluations allows for the opportunity to develop and implement methods to manage populations at the genomic level. As a result, the realized proportion of the genome that 2 individuals share can be more accurately estimated instead of using pedigree information to estimate the expected proportion of shared alleles. Furthermore, genomic information allows genome-wide relationship or inbreeding estimates to be augmented to characterize relationships for specific regions of the genome. Region-specific stretches can be used to more effectively manage areas of low genetic diversity or areas that, when homozygous, result in reduced performance across economically important traits. The use of region-specific metrics should allow breeders to more precisely manage the trade-off between the genetic value of the progeny and undesirable side effects associated with inbreeding. Methods tailored toward more effectively identifying regions affected by inbreeding and their associated use to manage the genome at the herd level, however, still need to be developed. We have reviewed topics related to inbreeding, measures of relatedness, genetic diversity and methods to manage populations at the genomic level, and we discuss future challenges related to managing populations through implementing genomic methods at the herd and population levels.

The largest gains in accuracy in a genomic selection program come from genotyping young selection candidates who have not yet produced progeny and who might, or might not, have a phenotypic record recorded. To reduce genotyping costs and to allow for an increased amount of genomic data to be available in a population, young selection candidates may be genotyped with low-density (LD) panels and imputed to a higher density. However, to ensure that a reasonable imputation accuracy persists overtime, some parent animals originally genotyped at LD must be re-genotyped at a higher density. This study investigated the long-term impact of selectively re-genotyping parents with a medium-density (MD) SNP panel on the accuracy of imputation and on the genetic predictions using ssGBLUP in a simulated beef cattle population. Assuming a moderately heritable trait (0.25) and a population undergoing selection, the simulation generated sequence data for a founder population (100 male and 500 female individuals) and 9,000 neutral markers, considered as the MD panel. All selection candidates from generation 8 to 15 were genotyped with LD panels corresponding to a density of 0.5% (LD_0.5), 2% (LD_2), and 5% (LD_5) of the MD. Re-genotyping scenarios chose parents at random or based on EBV and ranged from 10% of male parents to re-genotyping all male and female parents with MD. Ranges in average imputation accuracy at generation 15 were 0.567 to 0.936, 0.795 to 0.985, and 0.931 to 0.995 for the LD_0.5, LD_2, and LD_5, respectively, and the average EBV accuracies ranged from 0.453 to 0.735, 0.631 to 0.784, and 0.748 to 0.807 for LD_0.5, LD_2, and LD_5, respectively. Re-genotyping parents based on their EBV resulted in higher imputation and EBV accuracies compared to selecting parents at random and these values increased with the size of LD panels. Differences between re-genotyping scenarios decreased when the density of the LD panel increased, suggesting fewer animals needed to be re-genotyped to achieve higher accuracies. In general, imputation and EBV accuracies were greater when more parents were re-genotyped, independent of the proportion of males and females. In practice, the relationship between the density of the LD panel used and the target panel must be considered to determine the number (proportion) of animals that would need to be re-genotyped to enable sufficient imputation accuracy.

For genomic predictors to be of use in genetic evaluation, their predicted accuracy must be a reliable indicator of their utility, and thus unbiased. The objective of this paper was to evaluate the accuracy of prediction of genomic breeding values (GBV) using different clustering strategies and response variables. Red Angus genotypes (n = 9,763) were imputed to a reference 50K panel. The influence of clustering method [k-means, k-medoids, principal component (PC) analysis on the numerator relationship matrix (A) and the identical-by-state genomic relationship matrix (G) as both data and covariance matrices, and random] and response variables [deregressed estimated breeding values (DEBV) and adjusted phenotypes] were evaluated for cross-validation. The GBV were estimated using a Bayes C model for all traits. Traits for DEBV included birth weight (BWT), marbling (MARB), rib-eye area (REA), and yearling weight (YWT). Adjusted phenotypes included BWT, YWT, and ultrasonically measured intramuscular fat percentage and REA. Prediction accuracies were estimated using the genetic correlation between GBV and associated response variable using a bivariate animal model. A simulation mimicking a cattle population, replicated 5 times, was conducted to quantify differences between true and estimated accuracies. The simulation used the same clustering methods and response variables, with the addition of 2 genotyping strategies (random and top 25% of individuals), and forward validation. The prediction accuracies were estimated similarly, and true accuracies were estimated as the correlation between the residuals of a bivariate model including true breeding value (TBV) and GBV. Using the adjusted Rand index, random clusters were clearly different from relationship-based clustering methods. In both real and simulated data, random clustering consistently led to the largest estimates of accuracy, while no method was consistently associated with more or less bias than other methods. In simulation, random genotyping led to higher estimated accuracies than selection of the top 25% of individuals. Interestingly, random genotyping seemed to overpredict true accuracy while selective genotyping tended to underpredict accuracy. When forward in time validation was used, DEBV led to less biased estimates of GBV accuracy. Results suggest the highest, least biased GBV accuracies are associated with random genotyping and DEBV.

In livestock, the regulation of drugs used to treat livestock has received increased attention and it is currently unknown how much of the phenotypic variation in drug metabolism is due to the genetics of an animal. Therefore, the objective of the study was to determine the amount of phenotypic variation in fenbendazole and flunixin meglumine drug metabolism due to genetics. The population consisted of crossbred female and castrated male nursery pigs (n = 198) that were sired by boars represented by four breeds. The animals were spread across nine batches. Drugs were administered intravenously and blood collected a minimum of 10 times over a 48 h period. Genetic parameters for the parent drug and metabolite concentration within each drug were estimated based on pharmacokinetics (PK) parameters or concentrations across time utilizing a random regression model. The PK parameters were estimated using a non-compartmental analysis. The PK model included fixed effects of sex and breed of sire along with random sire and batch effects. The random regression model utilized Legendre polynomials and included a fixed population concentration curve, sex, and breed of sire effects along with a random sire deviation from the population curve and batch effect. The sire effect included the intercept for all models except for the fenbendazole metabolite (i.e., intercept and slope). The mean heritability across PK parameters for the fenbendazole and flunixin meglumine parent drug (metabolite) was 0.15 (0.18) and 0.31 (0.40), respectively. For the parent drug (metabolite), the mean heritability across time was 0.27 (0.60) and 0.14 (0.44) for fenbendazole and flunixin meglumine, respectively. The errors surrounding the heritability estimates for the random regression model were smaller compared to estimates obtained from PK parameters. Across both the PK and plasma drug concentration across model, a moderate heritability was estimated. The model that utilized the plasma drug concentration across time resulted in estimates with a smaller standard error compared to models that utilized PK parameters. The current study found a low to moderate proportion of the phenotypic variation in metabolizing fenbendazole and flunixin meglumine that was explained by genetics in the current study.

While for the most part genome-wide metrics are currently employed in managing livestock inbreeding, genomic data offer, in principle, the ability to identify functional inbreeding. Here we present a heuristic method to identify haplotypes contained within a run of homozygosity (ROH) associated with reduced performance. Results are presented for simulated and swine data. The algorithm comprises 3 steps. Step 1 scans the genome based on marker windows of decreasing size and identifies ROH genotypes associated with an unfavorable phenotype. Within this stage, multiple aggregation steps reduce the haplotype to the smallest possible length. In step 2, the resulting regions are formally tested for significance with the use of a linear mixed model. Lastly, step 3 removes nested windows. The effect of the unfavorable haplotypes identified and their associated haplotype probabilities for a progeny of a given mating pair or an individual can be used to generate an inbreeding load matrix (ILM). Diagonals of ILM characterize the functional inbreeding load of individual (IIL). We estimated the accuracy of predicting the phenotype based on ILL. We further compared the significance of the regression coefficient for IIL on phenotypes to genome-wide inbreeding metrics. We tested the algorithm using simulated scenarios (n =12) combining different levels of linkage disequilibrium (LD) and number of loci impacting a quantitative trait. Additionally, we investigated 9 traits from two maternal purebred swine lines. In simulated data, as the LD in the population increased the algorithm identified a greater proportion of the true unfavorable ROH effects. For example, the proportion of highly unfavorable true ROH effects identified raised from 32 to 41 % for the low to the high LD scenario. In both simulated and real data the haplotypes identified were contained within a much larger ROH (9.12-12.1 Mb). The IIL prediction accuracy was greater than zero across all scenarios for simulated data (high LD scenario mean (95% confidence interval): 0.49 (0.47-0.52)) and for nearly all swine traits (mean ± SD: 0.17±0.10). On average across simulated and swine datasets the IIL regression coefficient was more closely related to progeny performance than any genome-wide inbreeding metric. A heuristic method was developed that identified ROH genotypes with reduced performance and characterized the combined effects of ROH genotypes within and across individuals.

Drug use in livestock has received increased attention due to welfare concerns and food safety. Characterizing heterogeneity in the way swine populations respond to drugs could allow for group-specific dose or drug recommendations. Our objective was to determine whether drug clearance differs across genetic backgrounds and sex for sulfamethazine, enrofloxacin, fenbendazole and flunixin meglumine. Two sires from each of four breeds were mated to a common sow population. The nursery pigs generated (n = 114) were utilized in a random crossover design. Drugs were administered intravenously and blood collected a minimum of 10 times over 48 h. A non-compartmental analysis of drug and metabolite plasma concentration vs. time profiles was performed. Within-drug and metabolite analysis of pharmacokinetic parameters included fixed effects of drug administration date, sex and breed of sire. Breed differences existed for flunixin meglumine (P-value

Background In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome. Methods We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated. Results We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (−0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions. Conclusions We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions. Electronic supplementary material The online version of this article (doi:10.1186/s12711-016-0269-y) contains supplementary material, which is available to authorized users.

Identifying individual genetic variation in drug metabolism pathways is of importance not only in livestock, but also in humans in order to provide the ultimate goal of giving the right drug at the right dose at the right time. Our objective was to identify individual genes and gene networks involved in metabolizing fenbendazole (FBZ) and flunixin meglumine (FLU) in swine liver. The population consisted of female and castrated male pigs that were sired by boars represented by 4 breeds. Progeny were randomly placed into groups: no drug (UNT), FLU or FBZ administered. Liver transcriptome profiles from 60 animals with extreme (i.e. fast or slow drug metabolism) pharmacokinetic (PK) profiles were generated from RNA sequencing. Multiple cytochrome P450 (CYP1A1, CYP2A19 and CYP2C36) genes displayed different transcript levels across treated versus UNT. Weighted gene co-expression network analysis identified 5 and 3 modules of genes correlated with PK parameters and a portion of these were enriched for biological processes relevant to drug metabolism for FBZ and FLU, respectively. Genes within identified modules were shown to have a higher transcript level relationship (i.e. connectivity) in treated versus UNT animals. Investigation into the identified genes would allow for greater insight into FBZ and FLU metabolism.

“Screening schoolhood” attempts to move on both the debate and material covered by the six lead articles in this special issue. It appreciates the ways the articles, and the films they examine, deal with educational issues while at the same time evaluating their range in terms of creative construct and place in the history of documentary film. Considerate of the lead articles’ probing of their selected films as “artefacts” ripe for usage as data sources and as analytical frameworks for understanding, this article posits an addendum that develops around the fabricated nature of the films, their coherence as a category, their overlap with fictional films and their potential for formal, semiological and interdisciplinary enquiry.

Characterizing the variability in transcript levels across breeds and sex in swine for genes that play a role in drug metabolism may shed light on breed and sex differences in drug metabolism. The objective of the study is to determine if there is heterogeneity between swine breeds and sex in transcript levels for genes previously shown to play a role in drug metabolism for animals administered flunixin meglumine or fenbendazole. Crossbred nursery female and castrated male pigs (n = 169) spread across 5 groups were utilized. Sires (n = 15) of the pigs were purebred Duroc, Landrace, Yorkshire or Hampshire boars mated to a common sow population. Animals were randomly placed into the following treatments: no drug (control), flunixin meglumine, or fenbendazole. One hour after the second dosing, animals were sacrificed and liver samples collected. Quantitative Real-Time PCR was used to measure liver gene expression of the following genes: SULT1A1, ABCB1, CYP1A2, CYP2E1, CYP3A22 and CYP3A29. The control animals were used to investigate baseline transcript level differences across breed and sex. Post drug administration transcript differences across breed and sex were investigated by comparing animals administered the drug to the controls. Contrasts to determine fold change were constructed from a model that included fixed and random effects within each drug. Significant (P-value

Background Variation in environment, management practices, nutrition or selection objectives has led to a variety of different choices being made in the use of genetic material between countries. Differences in genome-level homozygosity between countries may give rise to regions that result in inbreeding depression to differ. The objective of this study was to characterize regions that have an impact on a runs of homozygosity (ROH) metric and estimate their association with the additive genetic effect of milk (MY), fat (FY) and protein yield (PY) and calving interval (CI) using Australia (AU) and United States (US) Jersey cows. Methods Genotyped cows with phenotypes on MY, FY and PY (n = 6751 US; n = 3974 AU) and CI (n = 5816 US; n = 3905 AU) were used in a two-stage analysis. A ROH statistic (ROH4Mb), which counts the frequency of a SNP being in a ROH of at least 4 Mb was calculated across the genome. In the first stage, residuals were obtained from a model that accounted for the portion explained by the estimated breeding value. In the second stage, these residuals were regressed on ROH4Mb using a single marker regression model and a gradient boosted machine (GBM) algorithm. The relationship between the additive and ROH4Mb of a region was characterized based on the (co)variance of 500 kb estimated genomic breeding values derived from a Bayesian LASSO analysis. Phenotypes to determine ROH4Mb and additive effects were residuals from the two-stage approach and yield deviations, respectively. Results Associations between yield traits and ROH4Mb were found for regions on BTA13, BTA23 and BTA25 for the US population and BTA3, BTA7, BTA17 for the AU population. Only one association (BTA7) was found for CI and ROH4Mb for the US population. Multiple potential epistatic interactions were characterized based on the GBM analysis. Lastly, the covariance sign between ROH4Mb and additive SNP effect of a region was heterogeneous across the genome. Conclusion We identified multiple genomic regions associated with ROH4Mb in US and AU Jersey females. The covariance of regions impacting ROH4Mb and the additive genetic effect were positive and negative, which provides evidence that the homozygosity effect is location dependent. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2001-7) contains supplementary material, which is available to authorized users.