Identification and prevalence of a genetic defect that causes leukocyte adhesion deficiency in Holstein cattle

Multiplex fluorescent amplification-refractory mutation system PCR method for the detection of 10 genetic defects in Holstein cattle and its comparison with the KASP genotyping assay

The common deleterious genetic defects in Holstein cattle include haplotypes 1-6 (HH1-HH6), haplotypes for cholesterol deficiency (HCD), bovine leukocyte adhesion deficiency (BLAD), complex vertebral malformation (CVM) and brachyspina syndrome (BS). Recessive inheritance patterns of these genetic defects permit the carriers to function normally, but homozygous recessive genotypes cause embryo loss or neonatal death. Therefore, rapid detection of the carriers is essential to manage these genetic defects. This study was conducted to develop a single-tube multiplex fluorescent amplification-refractory mutation system (mf-ARMS) PCR method for efficient genotyping of these 10 genetic defects and to compare its efficiency with the kompetitive allele specific PCR (KASP) genotyping assay. The mf-ARMS PCR method introduced 10 sets of tri-primers optimized with additional mismatches in the 3' end of wild and mutant-specific primers, size differentiation between wild and mutant-specific primers, fluorescent labeling of universal primers, adjustment of annealing temperatures and optimization of primer concentrations. The genotyping of 484 Holstein cows resulted in 16.12% carriers with at least one genetic defect, while no homozygous recessive genotype was detected. This study found carrier frequencies ranging from 0.0% (HH6) to 3.72% (HH3) for individual defects. The mf-ARMS PCR method demonstrated improved detection, time and cost efficiency compared with the KASP method for these defects. Therefore, the application of mf-ARMS PCR for genotyping Holstein cattle is anticipated to decrease the frequency of lethal alleles and limit the transmission of these genetic defects.

Genome-wide in silico analysis leads to identification of deleterious L290V mutation in RBBP5 gene in Bos indicus

Genome-wide deleterious mutations were identified in zebu cattle (Bos indicus) using in silico approach. The ddRAD sequence data of Sahiwal cattle were annotated and aligned with the cattle reference genome (ARS-UCD1.2). A total of 279,383 SNPs were identified at Read Depth10, which were further filtered to 692 missense SNPs. These SNPs were further analyzed, for functional consequences, by using Variant Effect Predictor, PolyPhen, PROVEAN, and PANTHER tools. A total of 18 SNPs, were finally identified as deleterious, and among these, 12 SNPs were mapped on nine different genes. ERRAT, ProSA-web, Project HOPE, TM-Align, and YASSARA tools, further confirmed the protein malfunctioning of one missense (L290V) mutation of Retinoblastoma binding protein-5 (RBBP5) gene, transcribing a cell cycle regulatory protein and associated with Retinoblastoma in human. This derived bioinformatics pipeline may be useful for preliminarily identifying the deleterious DNA mutations in livestock, specifically in absence of any genetic disease records.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects

BACKGROUND

Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs.

RESULTS

The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings.

CONCLUSIONS

Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Mating allocations in Holstein combining genomic information and linear programming optimization at the herd level

In this study, we explored mating allocation in Holstein using genomic information for 24,333 Holstein females born in Denmark, Finland, and Sweden. We used 2 data sets of bulls: the top 50 genotyped bulls and the top 25 polled genotyped bulls on the Nordic total merit scale. We used linear programming to optimize economic scores within each herd, considering genetic level, genetic relationship, semen cost, the economic impact of genetic defects, polledness, and β-casein. We found that it was possible to reduce genetic relationships and eliminate expression of genetic defects with minimal effect on the genetic level in total merit index. Compared with maximizing only Nordic total merit index, the relative frequency of polled offspring increased from 13.5 to 22.5%, and that of offspring homozygous for β-casein (A2A2) from 66.7 to 75.0% in one generation, without any substantial negative impact on other comparison criteria. Using only semen from polled bulls, which might become necessary if dehorning is banned, considerably reduced the genetic level. We also found that animals carrying the polled allele were less likely to be homozygous for β-casein (A2A2) and more likely to be carriers of the genetic defect HH1. Hence, adding economic value to a monogenic trait in the economic score used for mating allocation sometimes negatively affected another monogenetic trait. We recommend that the comparison criteria used in this study be monitored in a modern genomic mating program.

Identification of a putative haplotype associated with recumbency in Holstein calves

Thirty-four Holstein calves from multiple farms were found recumbent during the neonatal period with no detectable neurologic, infectious, or metabolic abnormalities. Most calves did not survive beyond 6 wk of age. The objective of this study was to conduct a genome-wide association and pedigree analysis to determine if a genetic origin was plausible. There were 101,917 DNA markers for 18 affected calves and 26 unaffected family controls available for analysis. Genome-wide association, homozygosity screening, and a parental based transmission disequilibrium test were conducted in PLINK. A genomic region on the end of chromosome 16 that contained 78 markers based on a recessive inheritance model and that spanned 5.1 million bp was considered the most probable region for a genetic defect; the region was narrowed to 2.1 million bp following homozygosity screening and the transmission disequilibrium test with all affected calves homozygous in the candidate region and 1 homozygous control. A genotyped sire and 2 dams with imputed genotypes were heterozygous in the candidate region. A common sire born in 2008 was identified that was present for both paternal and maternal lineages of all affected calves; nearly all lineages traced through a prolific son born in 2010 who was genotyped and was heterozygous for the candidate region. Therefore, a possible genetic defect with incomplete penetrance on chromosome 16 that results in recumbency has been identified. Further efforts with an increase in families represented are needed to confirm a genetic basis, and identify the mutation and mode of inheritance.

Inbreeding and pedigree analysis of the European red dairy cattle

BACKGROUND

Red dairy cattle breeds have an important role in the European dairy sector because of their functional characteristics and good health. Extensive pedigree information is available for these breeds and provides a unique opportunity to examine their population structure, such as effective population size, depth of the pedigree, and effective number of founders and ancestors, and inbreeding levels. Animals with the highest genetic contributions were identified. Pedigree data included 9,073,403 animals that were born between 1900 and 2019 from Denmark, Finland, Germany, Latvia, Lithuania, the Netherlands, Norway, Poland, and Sweden, and covered 32 breeds. The numerically largest breeds were Red Dairy Cattle and Meuse-Rhine-Yssel.

RESULTS

The deepest average complete generation equivalent (9.39) was found for Red Dairy Cattle in 2017. Mean pedigree completeness ranged from 0.6 for Finncattle to 7.51 for Red Dairy Cattle. An effective population size of 166 animals was estimated for the total pedigree and ranged from 35 (Rotes Höhenvieh) to 226 (Red Dairy Cattle). Average generation intervals were between 5 and 7 years. The mean inbreeding coefficient for animals born between 1960 and 2018 was 1.5%, with the highest inbreeding coefficients observed for Traditional Angler (4.2%) and Rotes Höhenvieh (4.1%). The most influential animal was a Dutch Meuse-Rhine-Yssel bull born in 1960. The mean inbreeding level for animals born between 2016 and 2018 was 2% and highest for the Meuse-Rhine-Yssel (4.64%) and Rotes Hohenvieh breeds (3.80%).

CONCLUSIONS

We provide the first detailed analysis of the genetic diversity and inbreeding levels of the European red dairy cattle breeds. Rotes Höhenvieh and Traditional Angler have high inbreeding levels and are either close to or below the minimal recommended effective population size, thus it is necessary to implement tools to monitor the selection process in order to control inbreeding in these breeds. Red Dairy Cattle, Vorderwälder, Swedish Polled and Hinterwälder hold more genetic diversity. Regarding the Meuse-Rhine-Yssel breed, given its decreased population size, increased inbreeding and low effective population size, we recommend implementation of a breeding program to prevent further loss in its genetic diversity.

Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle

Mendelian variants can determine both insemination success and neonatal survival and thus influence fertility and rearing success of cattle. We present 24 deficient homozygous haplotype regions in the Holstein population of Switzerland and provide an overview of the previously identified haplotypes in the global Holstein breed. This study encompasses massive genotyping, whole-genome sequencing (WGS) and phenotype association analyses. We performed haplotype screenings on almost 53 thousand genotyped animals including 114 k SNP data with two different approaches. We revealed significant haplotype associations to several survival, birth and fertility traits. Within haplotype regions, we mined WGS data of hundreds of bovine genomes for candidate causal variants, which were subsequently evaluated by using a custom genotyping array in several thousand breeding animals. With this approach, we confirmed the known deleterious SMC2:p.Phe1135Ser missense variant associated with Holstein haplotype (HH) 3. For two previously reported deficient homozygous haplotypes that show negative associations to female fertility traits, we propose candidate causative loss-of-function variants: the HH13-related KIR2DS1:p.Gln159* nonsense variant and the HH21-related NOTCH3:p.Cys44del deletion. In addition, we propose the RIOX1:p.Ala133_Glu142del deletion as well as the PCDH15:p.Leu867Val missense variant to explain the unexpected low number of homozygous haplotype carriers for HH25 and HH35, respectively. In conclusion, we demonstrate that with mining massive SNP data in combination with WGS data, we can map several haplotype regions and unravel novel recessive protein-changing variants segregating at frequencies of 1 to 5%. Our findings both confirm previously identified loci and expand the spectrum of undesired alleles impairing reproduction success in Holstein cattle, the world's most important dairy breed.

Pre-selection against a lethal recessive allele in breeding schemes with optimum-contribution selection or truncation selection

BACKGROUND

We tested the hypothesis that breeding schemes with a pre-selection step, in which carriers of a lethal recessive allele (LRA) were culled, and with optimum-contribution selection (OCS) reduce the frequency of a LRA, control rate of inbreeding, and realise as much genetic gain as breeding schemes without a pre-selection step.

METHODS

We used stochastic simulation to estimate true genetic gain realised at a 0.01 rate of true inbreeding (ΔF_true) by breeding schemes that combined one of four pre-selection strategies with one of three selection strategies. The four pre-selection strategies were: (1) no carriers culled, (2) male carriers culled, (3) female carriers culled, and (4) all carriers culled. Carrier-status was known prior to selection. The three selection strategies were: (1) OCS in which [Formula: see text] was predicted and controlled using pedigree relationships (POCS), (2) OCS in which [Formula: see text] was predicted and controlled using genomic relationships (GOCS), and (3) truncation selection of parents. All combinations of pre-selection strategies and selection strategies were tested for three starting frequencies of the LRA (0.05, 0.10, and 0.15) and two linkage statuses with the locus that has the LRA being on a chromosome with or without loci affecting the breeding goal trait. The breeding schemes were simulated for 10 discrete generations (t = 1, …, 10). In all breeding schemes, ΔF_true was calibrated to be 0.01 per generation in generations t = 4, …, 10. Each breeding scheme was replicated 100 times.

RESULTS

We found no significant difference in true genetic gain from generations t = 4, …, 10 between breeding schemes with or without pre-selection within selection strategy. POCS and GOCS schemes realised similar true genetic gains from generations t = 4, …, 10. POCS and GOCS schemes realised 12% more true genetic gain from generations t = 4, …, 10 than truncation selection schemes.

CONCLUSIONS

We advocate for OCS schemes with pre-selection against the LRA that cause animal suffering and high costs. At LRA frequencies of 0.10 or lower, OCS schemes in which male carriers are culled reduce the frequency of LRA, control rate of inbreeding, and realise no significant reduction in true genetic gain compared to OCS schemes without pre-selection against LRA.

Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins

Background

The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes.

Methods

In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits.

Results

In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a − 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles.

Conclusions

ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00660-z.

Prevalence of nine genetic defects in Chinese Holstein cattle

Worldwide use of elite sires has caused inbreeding accumulation and high frequencies of genetic defects in dairy cattle populations. In recent years, several genetic defect genes or haplotypes have been identified in Holstein cattle. A rapid and reliable microfluidic chip with Kompetitive allele-specific PCR (KASP) assay was developed in our previous study for the detection of heterozygotes at eight genetic defect loci of bovine leukocyte adhesion deficiency (BLAD), Brachyspina syndrome (BS), complex vertebral malformation (CVM), Holstein haplotype 1 (HH1), Holstein haplotype 3 (HH3), Holstein haplotype 4 (HH4), Holstein haplotype 5 (HH5) and haplotype for cholesterol deficiency (HCD). This study aimed to extend that assay to include a newly identified genetic defect of Holstein haplotype 6 (HH6) and to estimate the frequencies of carriers for each of the nine genetic defects in six Chinese Holstein herds. Of the 1633 cows, carrier frequencies of the genetic defects were 6.92%, 5.76%, 4.46%, 4.30%, 3.62%, 2.94%, 1.86% and 0.37% for HH1, HH3, CVM, HH5, HCD, BS, HH6 and BLAD, respectively. No carrier was found for HH4. Notably, 27.43% of cows carried at least one genetic defect, while 2.27% and 0.12% of cows carried double and triple genetic defect alleles, respectively. The existence of genetic defects calls for routine molecular testing and effective management of genetic defects by avoiding carrier-to-carrier mating in production herds and eliminating or at least reducing the frequency of the defective alleles through marker-assisted selection in breeding herds.

Induced pluripotent stem cells from farm animals

The development of the induced pluripotent stem cells (iPSCs) technology has revolutionized the world on the establishment of pluripotent stem cells (PSCs) across a great variety of animal species. Generation of iPSCs from domesticated animals would provide unrestricted cell resources for the study of embryonic development and cell differentiation of these species, for screening and establishing desired traits for sustainable agricultural production, and as veterinary and preclinical therapeutic tools for animal and human diseases. Induced PSCs from domesticated animals thus harbor enormous scientific, economical, and societal values. Although much progress has been made toward the generation of PSCs from these species, major obstacles remain precluding the exclamation of the establishment of bona fide iPSCs. The most prominent of them remain the inability of these cells to silence exogenous reprogramming factors, the obvious reliance on exogenous factors for their self-renewal, and the restricted development potential in vivo. In this review, we summarize the history and current progress in domestic farm animal iPSC generation, with a focus on swine, ruminants (cattle, ovine, and caprine), horses, and avian species (quails and chickens). We also discuss the problems associated with the farm animal iPSCs and potential future directions toward the complete reprogramming of somatic cells from farm animals.

Novel approach to incorporate information about recessive lethal genes increases the accuracy of genomic prediction for mortality traits

The genetic underpinnings of calf mortality can be partly polygenic and partly due to deleterious effects of recessive lethal alleles. Prediction of the genetic merits of selection candidates should thus take into account both genetic components contributing to calf mortality. However, simultaneously modeling polygenic risk and recessive lethal allele effects in genomic prediction is challenging due to effects that behave differently. In this study, we present a novel approach where mortality risk probabilities from polygenic and lethal allele components are predicted separately to compute the total risk probability of an individual for its future offspring as a basis for selection. We present methods for transforming genomic estimated breeding values of polygenic effect into risk probabilities using normal density and cumulative distribution functions and show computations of risk probability from recessive lethal alleles given sire genotypes and population recessive allele frequencies. Simulated data were used to test the novel approach as implemented in probit, logit, and linear models. In the simulation study, the accuracy of predicted risk probabilities was computed as the correlation between predicted mortality probabilities and observed calf mortality for validation sires. The results indicate that our novel approach can greatly increase the accuracy of selection for mortality traits compared with the accuracy of predictions obtained without distinguishing polygenic and lethal gene effects.

An IL17RA frameshift variant in a Holstein cattle family with psoriasis-like skin alterations and immunodeficiency

BACKGROUND

Skin lesions and dermatoses in cattle are often associated with infections due to bacteria, fungi or environmental risk factors. Dermatoses with genetic etiology have been described in cattle. Among these rare disorders, there are primary congenital dermatoses that are associated with inherited nutritional deficiencies, such as bovine hereditary zinc deficiency or zinc deficiency-like syndrome. This study presents three cases of Holstein cattle with congenital skin lesions observed on a single farm that resemble zinc deficiency-like syndrome. Close clinical and pathological examinations took place in two cases. Pedigree analysis indicated autosomal recessive inheritance and whole-genome sequencing of both affected calves was performed.

RESULTS

The two calves showed retarded growth and suffered from severe ulcerative dermatitis with hyperkeratosis, alopecia furunculosis and subcutaneous abscess formation. Blood analysis showed correspondent leukocytosis with neutrophilia whereas minerals, macro- and micronutrients were within the reference ranges. Variant calling and filtering against the 1000 Bull Genomes variant catalogue resulted in the detection of a single homozygous protein-changing variant exclusively present in both sequenced genomes. This single-nucleotide deletion in exon 3 of IL17RA on bovine chromosome 5 was predicted to have a deleterious impact on the encoded protein due to a frameshift leading to a truncated gene product. Genotyping of the affected cattle family confirmed recessive inheritance.

CONCLUSIONS

A loss-of-function mutation of the IL17RA transmembrane protein could be identified as most likely pathogenic variant for the psoriasis-like skin alterations observed in the two affected Holstein calves. In man, rare recessive diseases associated with IL17RA include immunodeficiency 51 and chronic mucocutaneous candidiasis. This supports the observed immunodeficiency of the presented cases. This study reports the first naturally occurring IL17RA-associated animal model.

Technical note: Development and application of KASP assays for rapid screening of 8 genetic defects in Holstein cattle

Specific DNA mutations underlying several genetic defects associated with embryo loss or reduced calf survivability have been identified in dairy cattle, and a convenient and cost-effective platform is required for their routine screening. We developed Kompetitive allele-specific PCR (KASP) assays for discrimination of the wild-type alleles from the associated defective alleles at each of 8 common genetic defects in Holstein cattle, involving 5 SNP [HH1, HH3, HH4, bovine leukocyte adhesion deficiency (BLAD), and complex vertebral malformation (CVM)] and 3 insertion or deletion mutations [HH5, haplotype for cholesterol deficiency (HCD), and brachyspina (BS)]. A total of 390 cows from a Chinese Holstein herd were genotyped and the carriers identified at 7 of these 8 loci (except HH4), with the highest carrier frequencies found for CVM (10.5%) and HH1 (10.0%), followed by HH3 (2.6%), BS (2.1%), HCD (1.3%), HH5 (0.8%), and BLAD (0.5%). Surprisingly, 102 cows (26.2%) carried at least 1 of the 7 defective alleles. Our results demonstrate that these KASP assays are simple, rapid, and reliable for the detection of multiple genetic defects. The high carrier frequency of these genetic defects indicates an urgent need for routine molecular testing to eliminate the deleterious alleles from Chinese Holstein cattle.

Genotypes and allele frequencies of buried SNPs in a bovine single-nucleotide polymorphism array in Japanese Black cattle

Single nucleotide polymorphism (SNP) arrays are widely used for genetic and genomic analyses in cattle breeding; thus, data derived from SNP arrays have accumulated on a large scale nationwide. Commercial SNP arrays contain a considerable number of unassigned SNPs on the chromosome/position on the genome; these SNPs are excluded in subsequent analyses. Notably, the position-unassigned SNPs, or "buried SNPs" include some of the markers associated with genetic disease. In this study, we identified the position of buried SNPs using the Basic Local Alignment Search Tool against the surrounding sequences and characterized the relationship between SNPs and genetic diseases in Online Mendelian Inheritance in Animals based on the genomic position. We determined the position of 285 buried SNPs on the genome and surveyed the genotype and allele frequencies of these SNPs in 5,955 individual Japanese Black cattle. Eleven SNPs associated with genetic disease, which contained five buried SNPs, were found in the population with the risk allele frequency ranging from 0.00008396 to 0.46. These results indicate that buried SNPs in the bovine SNP array can be utilized to identify associations with genetic disorders from large scale accumulated SNP genotype data in Japanese Black cattle.

Management of lethal recessive alleles in beef cattle through the use of mate selection software

BACKGROUND

Recessive loss-of-function (LOF) alleles at genes which are essential for life, can result in early embryonic mortality. Cattle producers can use the LOF carrier status of individual animals to make selection and mate allocation decisions.

METHODS

Two beef cattle breeding strategies i.e. (1) selection against LOF carriers as parents and (2) simultaneous selection and mate allocation to avoid the occurrence of homozygous offspring in three scenarios, which differed in number and frequency of LOF alleles were evaluated using the mate selection program, MateSel. Scenarios included (a) seven loci with high-frequency LOF alleles, (b) 76 loci with low-frequency LOF alleles, and (c) 50 loci with random high- and low-frequency LOF alleles. In addition, any savings resulting from the information obtained by varying the percentage (0-100%) of the herd genotyped, together with segregation analysis to cover ungenotyped animals, were calculated to determine (1) which percentage optimized net profit for a fixed cost of genotyping ($30/test), and (2) the breakeven cost for genotyping.

RESULTS

With full knowledge of the LOF alleles carried by selection candidates, the most profitable breeding strategy was always simultaneous selection and mate allocation to avoid homozygous affected offspring (aa) as compared to indiscriminate selection against carrier parents (Aa). The breakeven value of genotyping depended on the number of loci modeled, the LOF allele frequencies, and the mating/selection strategies used. Genotyping was most valuable when it was used to avoid otherwise high levels of embryonic mortalities. As the number of essential loci with LOF alleles increased, especially when some were present at relatively high minor allele frequencies, embryonic losses increased, and profit was maximized by genotyping 10 to 20% of a herd and using that information to reduce these losses.

CONCLUSIONS

Genotyping 100% of the herd was never the most profitable outcome in any scenario; however, genotyping some proportion of the herd, together with segregation analysis to cover ungenotyped animals, maximized overall profit in scenarios with large numbers of loci with LOF alleles. As more LOF alleles are identified, such a mate selection software will likely be required to optimally select and allocate matings to balance the rate of genetic gain, embryonic losses, and inbreeding.

The role of molecular genetics in livestock production

Genetic variations that lead to easy-to-identify phenotypic changes have always been of interest to livestock breeders since domestication. Molecular genetics has opened up possibilities for identifying these variations and understanding their biological and population effects. Moreover, molecular genetics is part of the most diverse approaches and applications in animal production nowadays, including paternity testing, selection based on genetic variants, diagnostic of genetic diseases, reproductive biotechniques, fraud identification, differentiation of hybrids, parasite identification, genetic evaluation, diversity studies, and genome editing, among others. Therefore, the objective of this review was to describe the different applications of molecular genetics in livestock production, contextualising them with examples and highlighting the importance of the study of these topics and their applications.

A bovine CD18 signal peptide variant with increased binding activity to Mannheimia hemolytica leukotoxin

Background: Mannheimia haemolytica is the major bacterial infectious agent of bovine respiratory disease complex and causes severe morbidity and mortality during lung infections. M. haemolytica secretes a protein leukotoxin (Lkt) that binds to the CD18 receptor on leukocytes, initiates lysis, induces inflammation, and causes acute fibrinous bronchopneumonia. Lkt binds the 22-amino acid CD18 signal peptide domain, which remains uncleaved in ruminant species. Our aim was to identify missense variation in the bovine CD18 signal peptide and measure the effects on Lkt binding.

Methods: Missense variants in the integrin beta 2 gene ( ITGB2) encoding CD18 were identified by whole genome sequencing of 96 cattle from 19 breeds, and targeted Sanger sequencing of 1238 cattle from 46 breeds. The ability of different CD18 signal peptide variants to bind Lkt was evaluated by preincubating the toxin with synthetic peptides and applying the mixture to susceptible bovine cell cultures in cytotoxicity-blocking assays.

Results: We identified 14 missense variants encoded on 15 predicted haplotypes, including a rare signal peptide variant with a cysteine at position 5 (C ₅) instead of arginine (R ₅). Preincubating Lkt with synthetic signal peptides with C ₅ blocked cytotoxicity significantly better than those with R ₅. The most potent synthetic peptide (C ₅PQLLLLAGLLA) had 30-fold more binding activity compared to that with R ₅.

Conclusions: The results suggest that missense variants in the CD18 signal peptide affect Lkt binding, and animals carrying the C ₅ allele may be more susceptible to the effects of Lkt. The results also identify a potent class of non-antibiotic Lkt inhibitors that could potentially protect cattle from cytotoxic effects during acute lung infections.

Invited review: Reproductive and genomic technologies to optimize breeding strategies for genetic progress in dairy cattle

Abstract. Dairy cattle breeders have exploited technological advances that have emerged in the past in regards to reproduction and genomics. The implementation of such technologies in routine breeding programs has permitted genetic gains in traditional milk production traits as well as, more recently, in low-heritability traits like health and fertility. As demand for dairy products increases, it is important for dairy breeders to optimize the use of available technologies and to consider the many emerging technologies that are currently being investigated in various fields. Here we review a number of technologies that have helped shape dairy breeding programs in the past and present, along with those potentially forthcoming. These tools have materialized in the areas of reproduction, genotyping and sequencing, genetic modification, and epigenetics. Although many of these technologies bring encouraging opportunities for genetic improvement of dairy cattle populations, their applications and benefits need to be weighed with their impacts on economics, genetic diversity, and society.

Feline leukocyte adhesion (CD18) deficiency caused by a deletion in the integrin β2 (ITGB2) gene

BACKGROUND

Leukocyte adhesion deficiency (LAD) or CD18 deficiency is an autosomal recessive immunodeficiency which has been described in people, cattle, dogs, and knockout mice.

OBJECTIVES

The study goals were to characterize the clinicopathologic, immunologic, and molecular genetic features of feline LAD (FLAD) in a neutered male adult Domestic Longhair cat with severe leukocytosis and recurrent infections.

METHODS

Flow cytometry evaluated surface expression of CD18 on neutrophils. In vitro functional assays assessed CD18-dependent neutrophil adhesion and T-cell proliferation. Genomic DNA and cDNA were used to identify a causative mutation in the coding sequence of the integrin β₂ subunit (ITGB2) gene.

RESULTS

The affected cat developed periodontitis during the first months of life followed by recurrent infections poorly responsive to antibiotic therapy, accompanied by extreme neutrophilia. Neutrophils from the proband, compared to feline controls, did not express any CD18 on the cell surface. Adhesion of affected neutrophils was severely impaired with and without phorbol-myristate-acetate activation. The proband's T-cells proliferated weakly to 1 pg but normally to 100 pg staphylococcal enterotoxin A, suggesting a CD18-independent T-cell response at higher doses. Molecular genetic analysis of the ITGB2 gene revealed a 24 bp deletion at the exon 2 to intron 2 boundary (c.46_58 + 11del), predicting premature translational termination due to abnormal splicing of exon 1 to exon 3 or 4.

CONCLUSIONS

Feline LAD exhibits features similar to LAD in other species. However, clinical episodes in FLAD appeared milder allowing for an extended life expectancy under long-term antimicrobial therapy, possibly due to an alternative, CD18-independent T-cell proliferation pathway.

Use of haplotypes to identify regions harbouring lethal recessive variants in pigs

Background

Lethal recessive genetic variants are maintained at relatively low frequencies in a population in the heterozygous state, but by definition are fatal and therefore unobserved in the homozygous state. Since haplotypes allow the tagging of rare and untyped genetic variants, they have potential for studying lethal recessive variants. In this study, we used a large commercial population to identify putative lethal recessive haplotypes that impact either the total number born (TNB) or the number born alive (NBA) as a proportion of the total number born (NBA/TNB). We also compared the use of haplotypes with a single nucleotide polymorphism (SNP)-by-SNP approach and examined the benefits of using additional haplotypes imputed from low-density genotype data for the detection of lethal recessive variants. Candidate haplotypes were identified using population-wide haplotype frequencies and within-family analyses. These candidate haplotypes were subsequently assessed for putative lethal recessive effects on TNB and NBA/TNB by comparing carrier-to-carrier matings with carrier-to-non-carrier matings.

Results

Using both medium-density and imputed low-density genotype data six regions were identified as containing putative lethal recessive haplotypes that had an effect on TNB. It is likely that these regions were related to at least four putative lethal recessive variants, each located on a different chromosome. Evidence for putative lethal recessive effects on TNB was found on chromosomes 1, 6, 10 and 14 using haplotypes. Using haplotypes from individuals genotyped only at medium-density or a SNP-by-SNP approach did not detect any lethal recessive effects. No lethal recessive haplotypes or SNPs were detected that had an effect on NBA/TNB.

Conclusions

We show that the use of haplotypes from combining medium-density and imputed low-density genotype data is superior for the identification of lethal recessive variants compared to both a SNP-by-SNP approach and to the use of only medium-density data. We developed a formal statistical framework that provided sufficient power to detect lethal recessive variants in species, which produce large full-sib families, while reducing false positive or type I errors. Applying this framework results in improvements in reproductive performance by purging lethal recessive alleles from a population in a timely and cost-effective manner.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0332-3) contains supplementary material, which is available to authorized users.

A Review of Selected Genes with Known Effects on Performance and Health of Cattle

There are genetic conditions that influence production in dairy and beef cattle. The objective of this review was to describe relevant genetic conditions that have been associated with productivity and health in cattle. Genes or genomic regions that have been identified as a candidate for the condition will be included, and the genetic basis of the condition will be defined. Genes and genetic conditions included in this review are bovine leukocyte adhesion deficiency, deficiency of the uridine monophosphate synthase, bovine chronic interstitial nephritis, horn development, myostatin, complex vertebral malformation, leptin, osteopetrosis, apoptosis peptide activating factor 1, chondrodysplastic dwarfism, caseins, calpastatin, umbilical hernia, lactoglobulin, citrullinemia, cholesterol deficiency, prions, thyroglobulin, diacylglycerol acyltransferase, syndactyly, maple syrup urine disease, slick hair, Factor XI deficiency, and μ-Calpain. This review is not meant to be comprehensive, and relevant information is provided to ascertain genetic markers associated with the conditions.

Use of SNP genotypes to identify carriers of harmful recessive mutations in cattle populations

BACKGROUND

SNP (single nucleotide polymorphisms) genotype data are increasingly available in cattle populations and, among other things, can be used to predict carriers of specific mutations. It is therefore convenient to have a practical statistical method for the accurate classification of individuals into carriers and non-carriers. In this paper, we compared - through cross-validation- five classification models (Lasso-penalized logistic regression -Lasso, Support Vector Machines with either linear or radial kernel -SVML and SVMR, k-nearest neighbors -KNN, and multi-allelic gene prediction -MAG), for the identification of carriers of the TUBD1 recessive mutation on BTA19 (Bos taurus autosome 19), known to be associated with high calf mortality. A population of 3116 Fleckvieh and 392 Brown Swiss animals genotyped with the 54K SNP-chip was available for the analysis.

RESULTS

In general, the use of SNP genotypes proved to be very effective for the identification of mutation carriers. The best predictive models were Lasso, SVML and MAG, with an average error rate, respectively, of 0.2 %, 0.4 % and 0.6 % in Fleckvieh, and 1.2 %, 0.9 % and 1.7 % in Brown Swiss. For the three models, the false positive rate was, respectively, 0.1 %, 0.1 % and 0.2 % in Fleckvieh, and 3.0 %, 2.4 % and 1.6 % in Brown Swiss; the false negative rate was 4.4 %, 7.6 %1.0 % in Fleckvieh, and 0.0 %, 0.1% and 0.8 % in Brown Swiss. MAG appeared to be more robust to sample size reduction: with 25 % of the data, the average error rate was 0.7 % and 2.2 % in Fleckvieh and Brown Swiss, compared to 2.1 % and 5.5 % with Lasso, and 2.6 % and 12.0 % with SVML.

CONCLUSIONS

The use of SNP genotypes is a very effective and efficient technique for the identification of mutation carriers in cattle populations. Very few misclassifications were observed, overall and both in the carriers and non-carriers classes. This indicates that this is a very reliable approach for potential applications in cattle breeding.

Development of a fast and economical genotyping protocol for bovine leukocyte adhesion deficiency (BLAD) in cattle

Fast and economical means of assaying SNP's are important in diagnostic assays, especially when a large number of animals have to be screened for a genetic disease. This study was aimed at the development of a fast and economical screening assay for bovine leukocyte adhesion deficiency (BLAD) which is an important genetic disease of cattle industry. Four primers were designed where the outer primers amplify a 354 bp amplicon of CD18 gene carrying the polymorphism responsible for BLAD. The specifically designed inner primers in conjunction with the modified reaction mixture and cyclic conditions ensured amplification of either of wild or mutated alleles. Together with outer primers, the inner primers generated typical banding pattern in agarose gel which discriminated the normal animal against the carrier. We successfully used this protocol in 200 bulls for genotyping the BLAD allele which confirmed by sequencing, showing a cent percentage concordance. With the developed assay the need for restriction digestion or use of costly equipment viz. real time PCR was eliminated. This genotyping assay ensured fast and economical genotyping and could be adopted in every laboratory with a minimum equipment requirement of thermocycler and gel documentation system.

Predicting population gene frequency from sample data

With an increase in the number of candidate genes for important traits in livestock, effective strategies for incorporating such genes into selection programmes are increasingly important. Those strategies in part depend on the frequency of a favoured allele in a population. Since comprehensive genotyping of a population is seldom possible, we investigate the consequences of sampling strategies on the reliability of the gene frequency estimate for a bi-allelic locus. Even within a subpopulation or line, often only a proportion of individuals will be genotype tested. However, through segregation analysis, probable genotypes can be assigned to individuals that themselves were not tested, using known genotypes on relatives and a starting (presumed) gene frequency. The value of these probable genotypes in estimation of gene frequency was considered. A subpopulation or line was stochastically simulated and sampled at random, over a cluster of years or by favouring a particular genotype. Line was simulated (replicated) 1000 times. The reliability of gene frequency estimates depended on the sampling strategy used. With random sampling, even when a small proportion of a line was genotyped (0·10), the gene frequency of the population was well estimated from the across-line mean. When information on probable genotypes on untested individuals was combined with known genotypes, the between-line variance in gene frequency was estimated well; including probable genotypes overcame problems of statistical sampling. When the sampling strategy favoured a particular genotype, unsurprisingly the estimate of gene frequency was biased towards the allele favoured. In using probable genotypes the bias was lessened but the estimate of gene frequency still reflected the sampling strategy rather than the true population frequency. When sampling was confined to a few clustered years, the estimation of gene frequency was biased for those generations preceding the sampling event, particularly when the presumed starting gene frequency differed from the true population gene frequency. The potential risks of basing inferences about a population from a potentially biased sample are discussed.

Molecular characterization and combined genotype association study of bovine cluster of differentiation 14 gene with clinical mastitis in crossbred dairy cattle

Aim:

The present study was undertaken with the objectives to characterize and to analyze combined genotypes of cluster of differentiation 14 (CD14) gene to explore its association with clinical mastitis in Karan Fries (KF) cows maintained in the National Dairy Research Institute herd, Karnal.

Materials and Methods:

Genomic DNA was extracted using blood of randomly selected 94 KF lactating cattle by phenol-chloroform method. After checking its quality and quantity, polymerase chain reaction (PCR) was carried out using six sets of reported gene-specific primers to amplify complete KF CD14 gene. The forward and reverse sequences for each PCR fragments were assembled to form complete sequence for the respective region of KF CD14 gene. The multiple sequence alignments of the edited sequence with the corresponding reference with reported Bos taurus sequence (EU148610.1) were performed with ClustalW software to identify single nucleotide polymorphisms (SNPs). Basic Local Alignment Search Tool analysis was performed to compare the sequence identity of KF CD14 gene with other species. The restriction fragment length polymorphism (RFLP) analysis was carried out in all KF cows using Helicobacter pylori 188I (Hpy188I) (contig 2) and Haemophilus influenzae I (HinfI) (contig 4) restriction enzyme (RE). Cows were assigned genotypes obtained by PCR-RFLP analysis, and association study was done using Chi-square (χ²) test. The genotypes of both contigs (loci) number 2 and 4 were combined with respect to each animal to construct combined genotype patterns.

Results:

Two types of sequences of KF were obtained: One with 2630 bp having one insertion at 616 nucleotide (nt) position and one deletion at 1117 nt position, and the another sequence was of 2629 bp having only one deletion at 615 nt position. ClustalW, multiple alignments of KF CD14 gene sequence with B. taurus cattle sequence (EU148610.1), revealed 24 nt changes (SNPs). Cows were also screened using PCR-RFLP with Hpy188I (contig 2) and HinfI (contig 4) RE, which revealed three genotypes each that differed significantly regarding mastitis incidence. The maximum possible combination of these two loci shown nine combined genotype patterns and it was observed only eight combined genotypes out of nine: AACC, AACD, AADD, ABCD, ABDD, BBCC, BBCD, and BBDD. The combined genotype ABCC was not observed in the studied population of KF cows. Out of 94 animals, AACD combined genotype animals (10.63%) were found to be not affected with mastitis, and ABDD combined genotyped animals was observed having the highest mastitis incidence of 15.96%.

Conclusion:

AACD typed cows were found to be least susceptible to mastitis incidence as compared to other combined genotypes.

Determination of Sheep Prion Gene Polymorphisms from Paraffin-Embedded Tissue

Amino acid polymorphisms of the prion protein (PrP) greatly influence the susceptibility of sheep to scrapie. Selective breeding to increase the prevalence of PrP gene alleles associated with scrapie resistance is a flock management practice that is important for scrapie control programs. Determination of sheep PrP alleles typically has required extraction of DNA from host tissues that are freshly derived or stored frozen. We describe application of a DNA extraction procedure for formalin-fixed, paraffin-embedded tissues (PET) for the purpose of PCR amplification and nucleotide sequencing of relevant codons (136–171) of the sheep PrP gene. Tissues derived from 96 sheep were studied. The DNA sequence identity was confirmed in 87 of 94 matched samples of PET and frozen tissue specimens. DNA from brainstem PET of 2 sheep, from which fresh tissue was not available, was amplified and sequenced after formalin fixation for 7–70 days. This method will allow retrospective analysis of PrP genetics of sheep subsequent to postmortem diagnosis of scrapie when nonfixed tissue is unavailable for DNA extraction; however, it is not recommended that submission of fixed tissue supplant collection of fresh tissues for the purpose of determining PrP gene polymorphisms.

A 0.5-Mbp deletion on bovine chromosome 23 is a strong candidate for stillbirth in Nordic Red cattle

Background

A whole-genome association study of 4631 progeny-tested Nordic Red dairy cattle bulls using imputed next-generation sequencing data revealed a major quantitative trait locus (QTL) that affects birth index (BI) on Bos taurus autosome (BTA) 23. We analyzed this QTL to identify which of the component traits of BI are affected and understand its molecular basis.

Results

A genome-wide scan of BI in Nordic Red dairy cattle detected major QTL on BTA6, 14 and 23. The strongest associated single nucleotide polymorphism (SNP) on BTA23 was located at 13,313,896 bp with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- \log_{10} ({\text{p}}) = 50.63$$\end{document}-log10(p)=50.63. Analyses of component traits showed that the QTL had a large effect on stillbirth. Based on the 10 most strongly associated SNPs with stillbirth, we constructed a haplotype. Among this haplotype’s alleles, HAP_QTL had a large negative effect on stillbirth. No animals were found to be homozygous for HAP_QTL. Analysis of stillbirth records that were categorized by carrier status for HAP_QTL of the sire and maternal grandsire suggested that this haplotype had a recessive mode of inheritance. Illumina BovineHD BeadChip genotypes and genotype intensity data indicated a chromosomal deletion between 12.28 and 12.81 Mbp on BTA23. An independent set of Illumina Bovine50k BeadChip genotypes identified a recessive lethal haplotype that spanned the deleted region.

Conclusions

A deleted region of approximately 500 kb that spans three genes on BTA23 was identified and is a strong candidate QTL with a large effect on BI by increasing stillbirth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0215-z) contains supplementary material, which is available to authorized users.

Novel identification of Factor XI deficiency in Indian Sahiwal (Bos indicus) cattle

Factor-XI deficiency (FXID) is inherited as autosomal lethal recessive disorder of carrier Holstein-Friesian bulls. A 76 base pair segment insertion into exon 12 in Factor-XI gene causes FXID in cattle. Keeping this in view the present study was conducted to screen breeding bulls of both indigenous and exotic breeds for mutation in Factor-XI gene and to find out the frequency of FXID carrier animals in breeding bulls. A total of 120 bulls of different age group maintained at Frozen Semen Bull Station, India were randomly selected from different cattle breeds to screen presence of FXID syndrome in breeding sires. Genomic DNA was isolated from blood of the selected bulls. PCR parameters were standardized to obtain 244 and 320 bp amplicons. The results showed that 2 Sahiwal bulls out of 120 animals were carrier for FXID. Amplicons of the carrier animals were sequenced and annoted, which confirms a 76 bp insertion in the exon 12. Bleeding and clotting time showed considerable discrepancy in the carrier animals as compared to the normal animals. The findings of relative mRNA expression of Factor XI transcript revealed identical tendency in the carrier. The frequency of carrier animals and mutant allele was 2.5 % and 0.025 respectively. This study recommends for screening of breeding at AI bull centers in the country for FXID. The study also stands a merit for identification of FXID carrier in Bos indicus for the first time.

Immune response to bovine viral diarrhea virus--looking at newly defined targets

Bovine viral diarrhea virus (BVDV) has long been associated with a wide variety of clinical syndromes and immune dysregulation, many which result in secondary bacterial infections. Current understanding of immune cell interactions that result in activation and tolerance are explored in light of BVDV infection including: depletion of lymphocytes, effects on neutrophils, natural killer cells, and the role of receptors and cytokines. In addition, we review some new information on the effect of BVDV on immune development in the fetal liver, the role of resident macrophages, and greater implications for persistent infection.

A simple strategy for managing many recessive disorders in a dairy cattle breeding program

Background

High-density single nucleotide polymorphism genotypes have recently been used to identify a number of novel recessive mutations that adversely affect fertility in dairy cattle, as well as to track other conditions such as red coat color and polled. Most current methods for mate allocation fail to consider this information, and it will become increasingly difficult to manage matings as the number of recessive mutations to be accounted for increases.

Methods

A modified version of a mating strategy that constrains inbreeding based on genomics (the Pryce method) was developed that also accounts for the economic effects of Mendelian disorders on overall economic merit (modified Pryce method) and compared with random mating, truncation selection, and the Pryce scheme. Several scenarios were considered, including scenarios with six hypothetical recessive alleles and 12 recessive alleles that are currently segregating in the US Holstein population.

Results

The Pryce method and the modified Pryce method showed similar ability to reduce frequencies of recessive alleles, particularly for loci with frequencies greater than 0.30. The modified Pryce method outperformed the Pryce method for low-frequency alleles with small economic value. Cumulative genetic gain for the selection objective was slightly greater when using the Pryce method, but rates of inbreeding were similar across methods.

Conclusions

The proposed method reduces allele frequencies faster than other methods, and also can be used to maintain or increase the frequency of desirable recessives. It can be easily implemented in software for mate allocation, and the code used in this study is freely available as a reference implementation.

Considering genetic characteristics in German Holstein breeding programs

Recently, several research groups have demonstrated that several haplotypes may cause embryonic loss in the homozygous state. Up to now, carriers of genetic disorders were often excluded from mating, resulting in a decrease of genetic gain and a reduced number of sires available for the breeding program. Ongoing research is very likely to identify additional genetic defects causing embryonic loss and calf mortality by genotyping a large proportion of the female cattle population and sequencing key ancestors. Hence, a clear demand is present to develop a method combining selection against recessive defects (e.g., Holstein haplotypes HH1-HH5) with selection for economically beneficial traits (e.g., polled) for mating decisions. Our proposed method is a genetic index that accounts for the allele frequencies in the population and the economic value of the genetic characteristic without excluding carriers from breeding schemes. Fertility phenotypes from routine genetic evaluations were used to determine the economic value per embryo lost. Previous research has shown that embryo loss caused by HH1 and HH2 occurs later than the loss for HH3, HH4, and HH5. Therefore, an economic value of € 97 was used against HH1 and HH2 and € 70 against HH3, HH4, and HH5. For polled, € 7 per polled calf was considered. Minor allele frequencies of the defects ranged between 0.8 and 3.3%. The polled allele has a frequency of 4.1% in the German Holstein population. A genomic breeding program was simulated to study the effect of changing the selection criteria from assortative mating based on breeding values to selecting the females using the genetic index. Selection for a genetic index on the female path is a useful method to control the allele frequencies by reducing undesirable alleles and simultaneously increasing economical beneficial characteristics maintaining most of the genetic gain in production and functional traits. Additionally, we applied the genetic index to real data and found a decrease of the genetic trend for the birth years 1990 to 2006. Since 2010 the genetic index has increased due to a strong increase in the polled frequency. However, further investigation is needed to better understand the biology to determine the correct time of embryo loss and the economic value of fertility disorders.

Genetic tools to improve reproduction traits in dairy cattle

Fertility is a major concern in the dairy cattle industry and has been the subject of numerous studies over the past 20 years. Surprisingly, most of these studies focused on rough female phenotypes and, despite their important role in reproductive success, male- and embryo-related traits have been poorly investigated. In recent years, the rapid and important evolution of technologies in genetic research has led to the development of genomic selection. The generalisation of this method in combination with the achievements of the AI industry have led to the constitution of large databases of genotyping and sequencing data, as well as refined phenotypes and pedigree records. These resources offer unprecedented opportunities in terms of fundamental and applied research. Here we present five such examples with a focus on reproduction-related traits: (1) detection of quantitative trait loci (QTL) for male fertility and semen quality traits; (2) detection of QTL for refined phenotypes associated with female fertility; (3) identification of recessive embryonic lethal mutations by depletion of homozygous haplotypes; (4) identification of recessive embryonic lethal mutations by mining whole-genome sequencing data; and (5) the contribution of high-density single nucleotide polymorphism chips, whole-genome sequencing and imputation to increasing the power of QTL detection methods and to the identification of causal variants.

Single nucleotide polymorphisms in candidate genes and their relation with somatic cell scores in Argentinean dairy cattle

The prevention and control of bovine mastitis by enhancing natural defenses in animals is important to improve the quality of dairy products. Mastitis resistance is a complex trait which depends on genetic components, as well as environmental and physiological factors. The limitations of classical control measures have led to the search for alternative approaches to minimize the use of antibiotics by selecting naturally resistant animals. Polymorphisms in genes associated with the innate immune system are strong candidates to be evaluated as genetic markers. In this work, we evaluated a set of single nucleotide polymorphisms (SNPs) in candidate genes for health and production traits, and determined their association with the somatic cell score (SCS) as an indicator of mastitis in Argentinean dairy cattle. We evaluated 941 cows: Holstein (n = 677) and Holstein × Jersey (n = 264) crossbred, daughters from 22 bulls from 14 dairy farms located in the central dairy area of Argentina. Two of the 21 successfully genotyped markers were found to be significantly associated (p < 0.05) with the SCS: GHR_140 and OPN_8514C-T. The heterozygote genotype for GHR_140 showed a favorable effect in reducing the SCS. On the other hand, heterozygote genotypes for OPN8514C-T caused an increase in the SCS; moreover, combined genotypes for OPN SNPs showed an even larger effect. These findings can contribute to the design of effective marker-assisted selection programs.

A limited number of Y chromosome lineages is present in North American Holsteins

Holsteins are the most numerous dairy cattle breed in North America and the breed has undergone intensive selection for improving milk production and conformation. Theoretically, this intensive selection could lead to a reduction of the effective population size and reduced genetic diversity. The objective of this study was to investigate the effective population size of the Holstein Y chromosome and the effects of limited Y chromosome lineages on male reproduction and the future of the breed. Paternal pedigree information of 62,897 Holstein bulls born between 1950 and 2013 in North America and 220,872 bulls evaluated by multiple-trait across-country genetic evaluations of Interbull (Uppsala, Sweden) were collected and analyzed. The results indicated that the number of Y chromosome lineages in Holsteins has undergone a dramatic decrease during the past 50 years because of artificial selection and the application of artificial insemination (AI) technology. All current Holstein AI bulls in North America are the descendants of only 2 ancestors (Hulleman and Neptune H) born in 1880. These 2 ancestral Y-lineages are continued through 3 dominant pedigrees from the 1960s; namely, Pawnee Farm Arlinda Chief, Round Oak Rag Apple Elevation, and Penstate Ivanhoe Star, with a contribution of 48.78, 51.06, and 0.16% to the Holstein bull population in the 2010s, respectively. The Y-lineage of Penstate Ivanhoe Star is almost eliminated from the breed. The genetic variations in the 2 ancestral Y-lineages were evaluated among 257 bulls by determining the copy number variations (CNV) of 3 Y-linked gene families: PRAMEY, HSFY, and ZNF280BY, which are spread along the majority (95%) of the bovine Y chromosome male-specific region (MSY). No significant difference was found between the 2 ancestral Y-lineages, although large CNV were observed within each lineage. This study suggests minimal genetic diversity on the Y chromosome in Holsteins and provides a starting point for investigating the effect of the extremely limited number of Y-lineages on male reproduction and other traits important for the future of the Holstein breed.

Prediction of expected genetic variation within groups of offspring for innovative mating schemes

Background

Experience from progeny-testing indicates that the mating of popular bull sires that have high estimated breeding values with excellent dams does not guarantee the production of offspring with superior breeding values. This is explained partly by differences in the standard deviation of gamete breeding values (SDGBV) between animals at the haplotype level. The SDGBV depends on the variance of the true effects of single nucleotide polymorphisms (SNPs) and the degree of heterozygosity. Haplotypes of 58 035 Holstein animals were used to predict and investigate expected SDGBV for fat yield, protein yield, somatic cell score and the direct genetic effect for stillbirth.

Results

Differences in SDGBV between animals were detected, which means that the groups of offspring of parents with low SDGBV will be more homogeneous than those of parents with high SDGBV, although the expected mean breeding values of the progeny will be the same. SDGBV was negatively correlated with genomic and pedigree inbreeding coefficients and a small loss of SDGBV over time was observed. Sires that had relatively low mean gamete breeding values but high SDGBV had a higher probability of producing extremely positive offspring than sires that had a high mean gamete breeding value and low SDGBV.

Conclusions

An animal’s SDGBV can be estimated based on genomic information and used to design specific genomic mating plans. Estimated SDGBV are an additional tool for mating programs, which allows breeders to identify and match mating partners using specific haplotype information.