Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNP-by-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNP-by-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.     

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.     

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (F_PED), marker-based using a genomic relationship matrix (F_GRM), and segment-based using ROH (F_ROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using F_PED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using F_GRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using F_ROH. In addition, the effective population size at past generations (5–100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000–2009, 2010–2014, and 2015–2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.     

Genomic and pedigree estimation of inbreeding depression for semen traits in the Basco-Béarnaise dairy sheep breed

Inbreeding depression is associated with a decrease in performance and fitness of the animals. The goal of this study was to evaluate pedigree-based and genomic methods to estimate the level of inbreeding and inbreeding depression for 3 semen traits (volume, concentration, and motility score) in the Basco-Béarnaise sheep breed. Data comprised 16,196 (or 15,071) phenotypic records from 620 rams (of which 533 rams had genotypes of 36,464 SNPs). The pedigree included 8,266 animals, composed of the 620 rams and their ancestors. The number of equivalent complete generations for the 620 rams was 7.04. Inbreeding coefficients were estimated using genomic and pedigree-based information. Genomic inbreeding coefficients were estimated from individual SNP and using segments of homozygous SNP (runs of homozygosity, ROH). Short ROH are of old origin, whereas long ROH are due to recent inbreeding. Considering that the equivalent number of generations in Basco-Béarnaise was 6, inbreeding coefficients for ROH with a length >4 Mb refer to all (recent + old) inbreeding, those with a length >17 Mb correspond to recent inbreeding, and the difference between them indicates old inbreeding. Pedigree-based inbreeding coefficients were also estimated classically, or accounting for nonzero relationships for unknown parents, or including metafounder relationships (estimated using markers) to account for missing pedigree information. Finally, inbreeding coefficients combining genotyped and nongenotyped animal information were computed from matrix H of the single-step approach, also including metafounders. Inbreeding depression was estimated differently depending on the approach used to compute inbreeding coefficients. These 8 estimators of inbreeding coefficients were included as covariates in different animal models. No inbreeding depression was detected for sperm volume or sperm concentration. Inbreeding depression was significant for the motility of spermatozoa. The effect of old and recent inbreeding on motility was null and negative, respectively, demonstrating the existence of purging by selection of deleterious recessive alleles affecting motility. A 10% increase in inbreeding would result in a reduction in mean motility ranging between 0.09 and 0.22 points in the score (from 0 to 5). Motility is unfavorably affected by increasing recent inbreeding but the impact is very small. Runs of homozygosity and metafounders allow us to accurately estimate inbreeding depression and detect recent inbreeding.     

Genomic inbreeding and relationships among Holsteins, Jerseys, and Brown Swiss

Genomic measures of relationship and inbreeding within and across breeds were compared with pedigree measures using genotypes for 43,385 loci of 25,219 Holsteins, 3,068 Jerseys, and 872 Brown Swiss. Adjustment factors allow genomic and pedigree relationships to match more closely within breeds and in multibreed populations and were estimated using means and regressions of genomic on pedigree relationships and allele frequencies in base populations. Correlations of genomic relationships with pedigree inbreeding were higher within each breed when an allele frequency of 0.5, rather than base population frequencies, was used, whereas correlations of average genomic relationships with average pedigree relationships and also reliabilities of genomic evaluations were higher using base population frequencies. Allele frequencies differed in the 3 breeds and were correlated by 0.65 to 0.67 when estimated from genotyped animals compared with 0.72 to 0.74 when estimated from breed base populations. The largest difference in allele frequency was between Holstein and the other breeds on chromosome Bos taurus autosome 4 near a gene affecting appearance of white skin patches (vitiligo) in humans. Each animal's breed composition was predicted very accurately with a standard deviation of <3% using regressions on genotypes at all loci or less accurately with a standard deviation of <6% using subsets of loci. Genomic future inbreeding (half an animal's mean genomic relationship to current animals of the same breed) was correlated by 0.75 to 0.94 with expected future inbreeding (half the average pedigree relationship). Correlations of both were slightly higher with parent averages than with genomic evaluations for net merit of young Holstein bulls. Thus, rates of increase in genomic and pedigree inbreeding per generation should be slightly reduced with genomic selection, in agreement with previous simulations. Genomic inbreeding and future inbreeding have been provided with individual genomic predictions since 2008. New methods to adjust pedigree and genomic relationship matrices so that they match may provide an improved basis for multibreed genomic evaluation. Positive definite matrices can be obtained by adjusting pedigree relationships for covariances among base animals within breed, whereas adjusting genomic relationships to match pedigree relationships can introduce negative eigenvalues. Pedigree relationship matrices ignore common ancestry shared by base animals within breed and may not approximate genomic relationships well in multibreed populations.     

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1–4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.     

Genomic inbreeding coefficients using imputed genotypes: Assessing different estimators in Holstein-Friesian dairy cows

The objective of this study was to estimate inbreeding coefficients in Holstein dairy cattle using imputed SNPs data. A data set of 95,540 Italian Holstein dairy cows from the routine genomic evaluations of the Italian National Association of Holstein, Brown, and Jersey Breeders were analyzed, with 84,445 imputed SNP. Ten widely used genomic inbreeding estimators were tested, including 4 PLINK v1.9 estimators (F, F_HAT1, F_HAT2, F_HAT3), 3 genomic relationship matrix (GRM)-based methods [VanRaden's first method with observed allele frequencies (F_GRM) or with fixed frequencies at 0.5 (F_GRM05), VanRaden's third method, allelic frequency free and pedigree regressed (F_GRM2)], runs of homozygosity (ROH)-based estimators in a complete (F_ROH) and simplified version (F_ROH2), and proportion of homozygous SNP (F_PH). Pairwise comparisons among them were made, including the comparison with traditional pedigree-based inbreeding coefficients (F_PED). Our results showed variability among the genomic inbreeding estimators. Coefficients of F_GRM and F_HAT3 were >1, meaning that more variability has been lost than the variability that existed in the base population. Regarding the remaining ones, F_GRM05, F_ROH, F_ROH2, and F_PH provided coefficients within the [0,1] space and are considered comparable to F_PED. Not comparable to F_PED, yet with an interpretable value, can be considered the coefficients of F, F_HAT2, and F_GRM2. Estimators based on ROH had the highest correlation with pedigree-based coefficients (0.59–0.66), among all estimators tested. In this study, Spearman correlations were shown to possibly provide a clearer estimation of the strength of the relationship between estimators. We hypothesize that imputation might cause extreme genomic inbreeding values that deserves further investigation.     

Invited review: Inbreeding in the genomics era: Inbreeding,inbreeding depression,and management of genomic variability

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.     

Symposium review: The genomic architecture of inbreeding: How homozygosity affects health and performance

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.     

Genome-wide association study using high-density single nucleotide polymorphism arrays and whole-genome sequences for clinical mastitis traits in dairy cattle

Mastitis is a mammary disease that frequently affects dairy cattle. Despite considerable research on the development of effective prevention and treatment strategies, mastitis continues to be a significant issue in bovine veterinary medicine. To identify major genes that affect mastitis in dairy cattle, 6 chromosomal regions on Bos taurus autosome (BTA) 6, 13, 16, 19, and 20 were selected from a genome scan for 9 mastitis phenotypes using imputed high-density single nucleotide polymorphism arrays. Association analyses using sequence-level variants for the 6 targeted regions were carried out to map causal variants using whole-genome sequence data from 3 breeds. The quantitative trait loci (QTL) discovery population comprised 4,992 progeny-tested Holstein bulls, and QTL were confirmed in 4,442 Nordic Red and 1,126 Jersey cattle. The targeted regions were imputed to the sequence level. The highest association signal for clinical mastitis was observed on BTA 6 at 88.97 Mb in Holstein cattle and was confirmed in Nordic Red cattle. The peak association region on BTA 6 contained 2 genes: vitamin D-binding protein precursor (GC) and neuropeptide FF receptor 2 (NPFFR2), which, based on known biological functions, are good candidates for affecting mastitis. However, strong linkage disequilibrium in this region prevented conclusive determination of the causal gene. A different QTL on BTA 6 located at 88.32 Mb in Holstein cattle affected mastitis. In addition, QTL on BTA 13 and 19 were confirmed to segregate in Nordic Red cattle and QTL on BTA 16 and 20 were confirmed in Jersey cattle. Although several candidate genes were identified in these targeted regions, it was not possible to identify a gene or polymorphism as the causal factor for any of these regions.     

A mating advice system in dairy cattle incorporating genomic information

This study investigated the effects of alternative mating programs that incorporate genomic information on expected progeny herd performance and inbreeding, as well as methods to include un-genotyped animals in such mating programs. A total of 54,535 Holstein-Friesian cattle with imputed high-density genotypes (547,650 SNP after edits) were available. First, to quantify the accuracy of imputing un-genotyped animals (often an issue in populations), a sub-population of 729 genotyped animals had their genotypes masked, and their allele dosages were imputed, using linear regression exploiting information on genotyped relatives. The reference population for imputation included all genotyped animals, excluding the 729 selected animals and their sires, dams, and grandsires, and had either (1) their sires' genotypes, (2) their dams' genotypes (3) both their sires' and their dams' genotypes, or (4) both their sires' and maternal grandsires' genotypes introduced into the reference population. The correlations between true genotypes and the imputed allele dosages ranged from 0.58 (sire only) to 0.68 (both sire and dam). A herd of 100 cows was then simulated (1,000 replicates) from the sub-population of 729 imputed animals. The top 10 bulls from the genotyped population, based on their total genetic merit index (TMI) were selected to be used as sires. Three mating allotment methods were investigated: (1) random mating, (2) sequential mating based on maximizing only the expected TMI of the progeny, and (3) linear programming to maximize a generated index constructed to maximize genetic merit and minimize expected progeny inbreeding as well as intra- and inter-progeny variability in genetic merit. Relationships among candidate parents were calculated using either the pedigree relationship matrix or the genomic relationship matrix; the latter was constructed using either the true genotypes of both parents or the true genotypes of the sire plus the imputed allele dosages of the dam. Using the genomic co-ancestry estimates resulted in lower average herd expected genomic inbreeding levels compared with using the pedigree-based co-ancestry estimates. Additionally, if the dams were not genotyped, using their imputed allele dosages also resulted in lower average herd expected inbreeding levels compared with using the pedigree co-ancestry estimates. The inter-progeny coefficient of variation for selected traits, milk and fertility, estimated breeding values were reduced by 12 to 65% using the linear programing method compared with sequential mating.     

Controlling bias in genomic breeding values for young genotyped bulls

The objectives of this study were to investigate bias in genomic predictions for dairy cattle and to find a practical approach to reduce the bias. The simulated data included phenotypes, pedigrees, and genotypes, mimicking a dairy cattle population (i.e., cows with phenotypes and bulls with no phenotypes) and assuming selection by breeding values or no selection. With the simulated data, genomic estimated breeding values (GEBV) were calculated with a single-step genomic BLUP and compared with true breeding values. Phenotypes and genotypes were simulated in 10 generations and in the last 4 generations, respectively. Phenotypes in the last generation were removed to predict breeding values for those individuals using only genomic and pedigree information. Complete pedigrees and incomplete pedigrees with 50% missing dams were created to construct the pedigree-based relationship matrix with and without inbreeding. With missing dams, unknown parent groups (UPG) were assigned in relationship matrices. Regression coefficients (b₁) and coefficients of determination (R²) of true breeding values on (G)EBV were calculated to investigate inflation and accuracy in GEBV for genotyped animals, respectively. In addition to the simulation study, 18 linear type traits of US Holsteins were examined. For the 18 type traits, b₁ and R² of GEBV with full data sets on GEBV with partial data sets for young genotyped bulls were calculated. The results from the simulation study indicated inflation in GEBV for genotyped males that were evaluated with only pedigree and genomic information under BLUP selection. However, when UPG for only pedigree-based relationships were included, the inflation was reduced, accuracy was highest, and genetic trends had no bias. For the linear type traits, when UPG for only pedigree-based relationships were included, the results were generally in agreement with those from the simulation study, implying less bias in genetic trends. However, when including no UPG, UPG in pedigree-based relationships, or UPG in genomic relationships, inflation and accuracy in GEBV were similar. The results from the simulation and type traits suggest that UPG must be defined accurately to be estimable and inbreeding should be included in pedigree-based relationships. In dairy cattle, known pedigree information with inbreeding and estimable UPG plays an important role in improving compatibility between pedigree-based and genomic relationship matrices, resulting in more reliable genomic predictions.     

Milk coagulation ability of five dairy cattle breeds

Samples of herd milk (506) were analyzed to assess sources of variation for milk coagulation properties (MCP) for 5 different dairy cattle breeds. Data were recorded in 55 single-breed dairy herds in the Trento province, a mountain area in northeast Italy. The 5 cattle breeds were Holstein-Friesian (8 herds), Brown Swiss (16 herds), Simmental (10 herds), Rendena (13 herds), and Alpine Gray (8 herds). Herd milk samples were analyzed for the MCP traits, milk rennet coagulation time (RCT), curd-firming time, and curd firmness (a₃₀), as well as protein and fat percentages, somatic cell count, Soxhlet-Henkel acidity, and bacterial count. An ANOVA was performed to study the effect of breed, herd within breed, DIM, month of lactation, protein and fat percentages, somatic cell score, titratable acidity, and log bacterial count within breed on MCP. Breed was the most important source of variation. In particular, the Rendena breed showed the best MCP traits at 13.5 min and 27.0 mm for RCT and a₃₀, respectively. The Holstein-Friesian breed had the worst coagulation properties at 18.0 min and 17.5 mm for RCT and a₃₀, respectively. The other 3 breeds showed intermediate coagulation properties. The RCT values were better at the beginning of lactation, whereas RCT and a₃₀ values were better in September and October (14.3 min and 25.7 mm, respectively). Among the composition traits, only the titratable acidity affected MCP traits of herd milk positively.     

A genome-wide association study for left-sided displacement of the abomasum using a high-density single nucleotide polymorphism array

Left-sided displacement of the abomasum (LDA) is a frequent disease in dairy cattle causing significant financial losses for dairy farmers. Heritability (h²) of this complex disease was estimated at up to 0.5 in German Holstein (GH) cattle. Using the Bovine High Density BeadChip (Illumina Inc., San Diego, CA) comprising 588,753 single nucleotide polymorphisms (SNP) after quality control for 126 LDA cases and 280 population-based controls, we used a mixed linear model analysis in a genome-wide association study (GWAS). We identified 6 genomic regions for LDA on bovine chromosomes 2, 8, 13, 20, 24, and X that were significantly associated with LDA. Each of these regions was covered by 4 to 12 LDA-associated SNP. Single SNP within these regions explained up to 7.3% of the phenotypic variance. An independent sample of 1,554 GH cows, including 539 controls and 1,015 cases, were genotyped for 8 SNP highly associated with LDA on Bos taurus autosomes (BTA) 2, 8, 13, and 24, as well as 6 SNP located in previously identified LDA regions on BTA1, 5, 11, and 27 using competitive allele-specific PCR genotyping technology (KASP). The analysis using the KASP genotypes confirmed LDA-associated loci on BTA2, 8, 13, and 27. These genomic regions may contribute to the susceptibility to LDA in Holstein cows and may harbor functional variants for LDA.     

Genome-wide association study identifies loci influencing natural antibody titers in milk of Dutch Holstein-Friesian cattle

Natural antibodies (NAb) are produced without any antigenic stimulation as a part of the innate immune system and provide a first line of defense against pathogens. Hence, they may be a useful trait when estimating an animal's potential immune competence and in selection for disease resistance. The aim of this study was to identify genomic regions associated with different NAb traits in milk and potentially describe candidate genes. Milk samples from 1,695 first-lactation Holstein Friesian cows with titer measurements for keyhole limpet hemocyanin, lipopolysaccharide, lipoteichoic acid, and peptidoglycan-binding total NAb and isotypes IgG₁, IgM, and IgA were used. Genome-wide association study analyses were performed using imputed 777K SNP genotypes, accounting for relationships using pedigree information. Functional enrichment analysis was performed on the significantly associated genomic regions to look for candidate genes. For IgM NAb, significant associations (false discovery rate <0.05) were found on Bos taurus autosome (BTA) 17, 18, and 21 with candidate genes related to immunoglobulin structure and early B cell development. For IgG₁, associations were found on BTA3, and we confirmed a quantitative trait loci on BTA21 previously reported for IgG NAb in serum. Our results provide new insights into the regulation of milk NAb that will help unravel the complex relationship between milk immunoglobulins and disease resistance in dairy cattle.     

Novel strategies to minimize progeny inbreeding while maximizing genetic gain using genomic information

In this study, 3 strategies for controlling progeny inbreeding in mating plans were compared. The strategies used information from pedigree inbreeding coefficients, genomic relationships, or shared runs of homozygosity. The strategies were compared for the reduction in genetic gain and progeny inbreeding that would be expected from selected matings, and for the decrease of homozygosity of deleterious recessive alleles. Using real pedigree, genotype [43,115 single nucleotide polymorphism (SNP) markers], and estimated breeding value data from Holstein cattle, mating plans were derived for herds of 300 cows with 20 sires available for mating, replicated 50 times. Each of the 300 individuals allocated as dams were matched to 1 of 20 sires to maximize genetic merit minus the penalty for estimated progeny inbreeding, and given the restriction that the sire could not be mated to more than 10% of the cows. The strategy that used a genomic relationship matrix (GRM) was the most effective in reducing average progeny inbreeding; this strategy also resulted in fewer homozygous SNP out of 1,000 low-frequency SNP compared with the strategy using pedigree information. In the future, large numbers of cattle may be genotyped for low-density SNP panels. A GRM constructed using 3,123 SNP produced results similar to a GRM constructed using the full 43,115 SNP. These results demonstrate that using GRM information, a 1% reduction in progeny inbreeding (valued at around $5 per cow) can be made with very little compromise in the overall breeding objective. These results and the availability of low-cost, low-density genotyping make it attractive to apply mating plans that use genomic information in commercial dairy herds.     

The genetic architecture of milk ELISA scores as an indicator of Johne's disease (paratuberculosis) in dairy cattle

Johne's disease (or paratuberculosis), caused by Mycobacterium avium ssp. paratuberculosis (MAP) infection, is a globally prevalent disease with severe economic and welfare implications. With no effective treatment available, understanding the role of genetics influencing host infection status is essential to develop selection strategies to breed for increased resistance to MAP infection. The main objectives of this study were to estimate genetic parameters for the MAP-specific antibody response using milk ELISA scores in Canadian Holstein cattle as an indicator of resistance to Johne's disease, and to unravel genomic regions and candidate genes significantly associated with MAP infection. After data editing, 168,987 milk ELISA records from 2,306 herds, obtained from CanWest Dairy Herd Improvement, were used for further analyses. Variance and heritability estimates for MAP infection status were determined using univariate linear animal models under 3 scenarios: (a) SCEN1: the complete data set (all herds); (b) SCEN2: herds with at least one suspect or test-positive animal (ELISA optical density ≥0.07); and (c) SCEN3: herds with at least one test-positive animal (ELISA optical density ≥0.11). Heritability estimates were calculated as 0.066, 0.064, and 0.063 for SCEN1, SCEN2, and SCEN3, respectively. The correlations between estimated breeding values for resistance to MAP infection and other economically important traits, when significant, were favorable and of low magnitude. Genome-wide association analyses identified important genomic regions on Bos taurus autosome (BTA)1, BTA7, BTA9, BTA14, BTA15, BTA17, BTA19, and BTA25 showing significant association with MAP infection status. These regions included 2 single nucleotide polymorphisms located 2 kb upstream of positional candidate genes CD86 and WNT9B, which play key roles in host immune response and tissue homeostasis. This study revealed the genetic architecture of MAP infection in Canadian Holstein cattle as measured by milk ELISA scores by estimating genetic parameters along with the identification of genomic regions potentially influencing MAP infection status. These findings will be of significant value toward implementing genetic and genomic evaluations for resistance to MAP infection in Holstein cattle.     

Mating programs including genomic relationships and dominance effects

Computerized mating programs using genomic information are needed by breed associations, artificial-insemination organizations, and on-farm software providers, but such software is already challenged by the size of the relationship matrix. As of October 2012, over 230,000 Holsteins obtained genomic predictions in North America. Efficient methods of storing, computing, and transferring genomic relationships from a central database to customers via a web query were developed for approximately 165,000 genotyped cows and the subset of 1,518 bulls whose semen was available for purchase at that time. This study, utilizing 3 breeds, investigated differences in sire selection, methods of assigning mates, the use of genomic or pedigree relationships, and the effect of including dominance effects in a mating program. For both Jerseys and Holsteins, selection and mating programs were tested using the top 50 marketed bulls for genomic and traditional lifetime net merit as well as 50 randomly selected bulls. The 500 youngest genotyped cows in the largest herd in each breed were assigned mates of the same breed with limits of 10 cows per bull and 1 bull per cow (only 79 cows and 8 bulls for Brown Swiss). A dominance variance of 4.1 and 3.7% was estimated for Holsteins and Jerseys using 45,187 markers and management group deviation for milk yield. Sire selection was identified as the most important component of improving expected progeny value, followed by managing inbreeding and then inclusion of dominance. The respective percentage gains for milk yield in this study were 64, 27, and 9, for Holsteins and 73, 20, and 7 for Jerseys. The linear programming method of assigning a mate outperformed sequential selection by reducing genomic or pedigree inbreeding by 0.86 to 1.06 and 0.93 to 1.41, respectively. Use of genomic over pedigree relationship information provided a larger decrease in expected progeny inbreeding and thus greater expected progeny value. Based on lifetime net merit, the economic value of using genomic relationships was >$3 million per year for Holsteins when applied to all genotyped females, assuming that each will provide 1 replacement. Previous mating programs required transferring only a pedigree file to customers, but better service is possible by incorporating genomic relationships, more precise mate allocation, and dominance effects. Economic benefits will continue to grow as more females are genotyped.     

Genome-wide association study for milk production traits in a Brazilian Holstein population

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.     

Association analysis for young stock survival index with imputed whole-genome sequence variants in Nordic Holstein cattle

Identification of the genetic variants associated with calf survival in dairy cattle will aid in the elimination of harmful mutations from the cattle population and the reduction of calf and young stock mortality rates. We used de-regressed estimated breeding values for the young stock survival (YSS) index as response variables in a genome-wide association study with imputed whole-genome sequence variants. A total of 4,610 bulls with estimated breeding values were genotyped with the Illumina BovineSNP50 (Illumina, San Diego, CA) single nucleotide polymorphism (SNP) genotyping array. Genotypes were imputed to whole-genome sequence variants. After quality control, 15,419,550 SNP on 29 Bos taurus autosomes (BTA) were used for association analysis. A modified mixed-model association analysis was used for a genome scan, followed by a linear mixed-model analysis for selected genetic variants. We identified 498 SNP on BTA5 and BTA18 that were associated with the YSS index in Nordic Holstein. The SNP rs440345507 (Chr5:94721790) on BTA5 was the putative causal mutation affecting YSS. Two haplotype-based models were used to identify haplotypes with the largest detrimental effects on YSS index. For each association signal, 1 haplotype region with harmful effects and the lead associated SNP were identified. Detected haplotypes on BTA5 and BTA18 explained 1.16 and 1.20%, respectively, of genetic variance for the YSS index. We examined whether YSS quantitative trait loci (QTL) on BTA5 and BTA18 were associated with stillbirth. YSS QTL on BTA18 overlapped a QTL region for stillbirth, but most likely 2 different causal variants were responsible for these 2 QTL. Four component traits of the YSS index, defined by sex and age, were analyzed separately by the modified mixed-model approach. The same genomic regions were associated with both bull and heifer calf mortality. Several genes (EPS8, LOC100138951, and KLK family genes) contained a lead associated SNP or were included in haplotypes with large detrimental effects on YSS in Nordic Holstein cattle.