Population dynamics of potentially harmful haplotypes: a pedigree analysis

BACKGROUND

The identification of low-frequency haplotypes, never observed in homozygous state in a population, is considered informative on the presence of potentially harmful alleles (candidate alleles), putatively involved in inbreeding depression. Although identification of candidate alleles is challenging, studies analyzing the dynamics of potentially harmful alleles are lacking. A pedigree of the highly endangered Gochu Asturcelta pig breed, including 471 individuals belonging to 51 different families with at least 5 offspring each, was genotyped using the Axiom PigHDv1 Array (658,692 SNPs). Analyses were carried out on four different cohorts defined according to pedigree depth and at the whole population (WP) level.

RESULTS

The 4,470 Linkage Blocks (LB) identified in the Base Population (10 individuals), gathered a total of 16,981 alleles in the WP. Up to 5,466 (32%) haplotypes were statistically considered candidate alleles, 3,995 of them (73%) having one copy only. The number of alleles and candidate alleles varied across cohorts according to sample size. Up to 4,610 of the alleles identified in the WP (27% of the total) were present in one cohort only. Parentage analysis identified a total of 67,742 parent-offspring incompatibilities. The number of mismatches varied according to family size. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the LB in which they were located. Segregation analyses informed that most potential candidate alleles appeared de novo in the pedigree. Only 17 candidate alleles were identified in the boar, sow, and paternal and maternal grandparents and were considered segregants.

CONCLUSIONS

Our results suggest that neither mutation nor recombination are the major forces causing the apparition of candidate alleles. Their occurrence is more likely caused by Allele-Drop-In events due to SNP calling errors. New alleles appear when wrongly called SNPs are used to construct haplotypes. The presence of candidate alleles in either parents or grandparents of the carrier individuals does not ensure that they are true alleles. Minimum Allele Frequency thresholds may remove informative alleles. Only fully segregant candidate alleles should be considered potentially harmful alleles. A set of 16 candidate genes, potentially involved in inbreeding depression, is described.

The effect of genetic defects on pregnancy loss in Swedish dairy cattle

The effect of carrier status of 10 lethal recessive genetic defects on pregnancy maintenance in Swedish dairy cattle was examined. The genetic defects were Ayrshire Haplotype 1, Ayrshire Haplotype 2, BTA12, BTA23, and Brown Swiss Haplotype 2 in Red Dairy Cattle (RDC), and Holstein Haplotype 1, 3, 4, 6, and 7 (HH1-HH7) in Holstein. Effects of carrier status of BTA12 and HH3 on conception rate (CR), interval from first to last service (FLS), and milk production were also examined. Data were obtained for 1,429 herds in the Swedish milk recording system, while information on carrier status of genetic defects was obtained from the Nordic Cattle Genetic Evaluation. In total, data on 158,795 inseminations in 28,432 RDC and 22,018 Holstein females were available. Data permitted separate analyses of BTA12 and HH3, but carrier frequencies of other defects were too low to enable further analysis. Pregnancy loss was defined as failure to maintain pregnancy, where pregnancy status was confirmed with manual and chemical pregnancy diagnosis, insemination, calving, sales and culling data. Odds ratios (OR) and probabilities of pregnancy loss and CR were estimated using generalized linear mixed models, while pregnancy loss, CR, FLS, milk, protein, and fat yields were analyzed using linear mixed models. Pregnancy losses were reported on average within the first month post-AI. At-risk matings were more prone to suffer pregnancy loss in BTA12 (OR = 1.79) and HH3 carriers (OR = 1.77) than not-at-risk matings. At-risk matings also had lower CR (OR = 0.62 and 0.63 for BTA12 and HH3, respectively) than not-at-risk matings. Carrier females of BTA12 had longer FLS and higher milk production than noncarriers. Conception rate and pregnancy maintenance could be improved by avoiding at-risk matings. This finding could help reduce pregnancy loss due to genetic defects in the breeding program for improved fertility.

Role of genetic introgression in introducing mutant alleles in Bos indicus cattle and prevalence of lethal genetic disorders in Bos taurus × Bos indicus and Bos indicus cattle in India

Fertility is an important trait associated with reproductive performance and animal welfare concern. Lethal alleles affect fertility through early embryonic death, abortions, and stillbirth depending on the genetic expression of the allele. Holstein Friesian and Jersey are two major Bos taurus breeds used widely for increasing milk yield along with purebreds of Bos indicus breeds like Gir, Kankrej, Sahiwal, and Tharparkar. In the present study, prevalence of lethal mutants in crossbred Holstein Friesian (CBHF, n = 2435), crossbred Jersey (CBJY, n = 2874), Gir (n = 3288), Kankrej (n = 593), Sahiwal (n = 965), and Tharparkar (n = 18) were studied. Heterozygous carrier animals were identified for bovine leukocyte adhesion deficiency (BLAD), Citrullinemia, complex vertebral malformation (CVM), Brachyspina, Holstein Haplotype 1 (HH1), Holstein Haplotype 3 (HH3),Holstein Haplotype 4 (HH4) and Jersey Haplotype 1 (JH1). Breed purity analysis confirmed inheritance of Bos taurus genes contributing to the presence of lethal mutant alleles like BLAD, Citrullinemia, HH1, and JH1 in apparently phenotypic Bos indicus animals. Screening and elimination of heterozygous carrier bulls/cows is essential to control fertility loss associated with lethal alleles.

Genomic evaluation of late-term abortion in cows recorded through Dairy Herd Improvement test plans

Late-term abortions cause significant economic loss and are of great concern for dairy herds. Late-term abortions >152 d and <251 d of gestation that terminate a lactation or initiate a new lactation have long been recorded by Dairy Herd Improvement (DHI). For 24.8 million DHI lactations, the average recorded incidence of late-term abortions across all years (2001-2018) was 1.2%. However, the 1.3% incidence of abortions reported in 2012 has declined to <1.0% incidence since 2015. Small adjustments were applied to the 82 million daughter pregnancy rate (DPR), 29 million cow conception rate (CCR), and 9 million heifer conception rate (HCR) records to account for late-term abortions more accurately. Fertility credits for CCR and HCR were changed to treat the last breeding as a failure instead of success if the next calving was coded as a late-term abortion. Similarly, when computing DPR, days open is now set to the maximum value of 250 instead of the reported days open if the next reported calving is an abortion. The test of these changes showed very small changes in standard deviation and high correlations (0.997) of adjusted predicted transmitting abilities (PTA) with official PTA from about 20,000 Holstein bulls born since 2000 with >50% reliability. For late-term abortion as a trait, estimated heritability was only 0.001 and PTA had a standard deviation of only 0.1% for recent sires with high reliability (>75%). Young animal genomic PTA have near 50% reliability but range only from -0.5 to +0.4 because of the low incidence and heritability. Genetic trend was slightly favorable and late-term abortion PTA were correlated favorably by 0.27 with net merit, 0.49 with productive life, 0.33 with livability, 0.23 with CCR, 0.20 with HCR, 0.26 with DPR, -0.31 with somatic cell score, -0.24 with daughter stillbirth, and -0.26 with daughter dystocia. Thus, PTA for late-term abortions should not be needed as a separate fertility trait and instead these minor edit changes should suffice. The PTA for late-term abortions would add little value because national evaluations for current fertility traits already account for those economic losses.

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.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

Truncation of IFT80 causes early embryonic loss in Holstein cattle associated with Holstein haplotype 2

Recessive alleles represent genetic risk in populations that have undergone bottleneck events. We present a comprehensive framework for identification and validation of these genetic defects, including haplotype-based detection, variant selection from sequence data, and validation using knockout embryos. Holstein haplotype 2 (HH2), which causes embryonic death, was used to demonstrate the approach. Holstein haplotype 2 was identified using a deficiency-of-homozygotes approach and confirmed to negatively affect conception rate and stillbirths. Five carriers were present in a group of 183 sequenced Holstein bulls selected to maximize the coverage of unique haplotypes. Three variants concordant with haplotype calls were found in HH2: a high-priority frameshift mutation resulting, and 2 low-priority variants (1 synonymous variant, 1 premature stop codon). The frameshift in intraflagellar 80 (IFT80) was confirmed in a separate group of Holsteins from the 1000 Bull Genomes Project that shared no animals with the discovery set. IFT80-null embryos were generated by truncating the IFT80 transcript at exon 2 or 11 using a CRISPR-Cas9 system. Abattoir-derived oocytes were fertilized in vitro, and zygotes were injected at the one-cell stage either with a guide RNA and CAS9 mRNA complex (n = 100) or Cas9 mRNA (control, n = 100) before return to culture, and replicated 3 times. IFT80 is activated at the 8-cell stage, and IFT80-null embryos arrested at this stage of development, which is consistent with data from mouse hypomorphs and HH2 carrier-to-carrier matings. This frameshift in IFT80 on chromosome 1 at 107,172,615 bp (p.Leu381fs) disrupts WNT and hedgehog signaling, and is responsible for the death of homozygous embryos.

Genomic Prediction for Abortion in Lactating Holstein Dairy Cows

Abortion in dairy cattle causes great economic losses due to reduced animal health, increase in culling rates, reduction in calf production, and milk yield, among others. Although the etiology of abortions can be of various origins, previous research has shown a genetic component. The objectives of this study were to (1) describe the development of the genomic prediction for cow abortions in lactating Holstein dairy cattle based on producer-recorded data and ssGBLUP methodology and (2) evaluate the efficacy of genomic predictions for cow abortions in commercial herds of US Holstein cows using data from herds that do not contribute phenotypic information to the evaluation. We hypothesized that cows with greater genomic predictions for cow abortions (Z_Abort STA) would have a reduced incidence of abortion. Phenotypic data on abortions, pedigree, and genotypes were collected directly from commercial dairy producers upon obtaining their permission. Abortion was defined as the loss of a confirmed pregnancy after 42 and prior to 260 days of gestation, treated as a binary outcome (0, 1), and analyzed using a threshold model. Data from a different subset of animals were used to test the efficacy of the prediction. The additive genetic variance for the cow abortion trait (Z_Abort) was 0.1235 and heritability was 0.0773. For all animals with genotypes (n = 1,662,251), mean reliability was 42%, and genomic predicted transmitting abilities (gPTAs) ranged from −8.8 to 12.4. Z_Abort had a positive correlation with cow and calf health traits and reproductive traits, and a negative correlation with production traits. Z_Abort effectively identified cows with a greater or lesser risk of abortion (16.6% vs. 11.0% for the worst and best genomics groups, respectively; p < 0.0001). The inclusion of cow abortion genomic predictions in a multi-trait selection index would allow dairy producers and consultants to reduce the incidence of abortion and to select high-producing, healthier, and more profitable cows.

Allele-biased expression of the bovine APOB gene associated with the cholesterol deficiency defect suggests cis-regulatory enhancer effects of the LTR retrotransposon insertion

The insertion of an endogenous retroviral long terminal repeat (LTR) sequence into the bovine apolipoprotein B (APOB) gene is causal to the inherited genetic defect cholesterol deficiency (CD) observed in neonatal and young calves. Affected calves suffer from developmental abnormalities, symptoms of incurable diarrhoea and often die within weeks to a few months after birth. Neither the detailed effects of the LTR insertion on APOB expression profile nor the specific mode of inheritance nor detailed phenotypic consequences of the mutation are undisputed. In our study, we analysed German Holstein dairy heifers at the peak of hepatic metabolic load and exposed to an additional pathogen challenge for clinical, metabolic and hepatic transcriptome differences between wild type (CDF) and heterozygote carriers of the mutation (CDC). Our data revealed that a divergent allele-biased expression pattern of the APOB gene in heterozygous CDC animals leads to a tenfold higher expression of exons upstream and a decreased expression of exons downstream of the LTR insertion compared to expression levels of CDF animals. This expression pattern could be a result of enhancer activity induced by the LTR insertion, in addition to a previously reported artificial polyadenylation signal. Thus, our data support a regulatory potential of mobile element insertions. With regard to the phenotype generated by the LTR insertion, heterozygote CDC carriers display significantly differential hepatic expression of genes involved in cholesterol biosynthesis and lipid metabolism. Phenotypically, CDC carriers show a significantly affected lipomobilization compared to wild type animals. These results reject a completely recessive mode of inheritance for the CD defect, which should be considered for selection decisions in the affected population. Exemplarily, our results illustrate the regulatory impact of mobile element insertions not only on specific host target gene expression but also on global transcriptome profiles with subsequent biological, functional and phenotypic consequences in a natural in-vivo model of a non-model mammalian organism.

A review of inbreeding depression in dairy cattle: current status, emerging control strategies, and future prospects

Dairy cattle breeding has historically focused on relatively small numbers of elite bulls as sires of sons. In recent years, even if generation intervals were reduced and more diverse sires of sons could have been selected, genomic selection has not fundamentally changed the fact that a large number of individuals are being analyzed. However, a relatively small number of elite bulls are still siring those animals. Therefore inbreeding-derived negative consequences in the gene pool have brought concern. The detrimental effects of non-additive genetic changes such as inbreeding depression and dominance have been widely disseminated while seriously affecting bioeconomically important parameters because of an antagonistic relationship between dairy production and reproductive traits. Therefore, the estimation of benefits and limitations of inbreeding and variance of the selection response deserves to be evaluated and discussed to preserve genetic variability, a significant concern in the selection of individuals for reproduction and production. Short-term strategies for genetic merit improvement through modern breeding programs have severely lowered high-producing dairy cattle fertility potential. Since the current selection programs potentially increase long-term costs, genetic diversity has decreased globally as a consequence. Therefore, a greater understanding of the potential that selection programs have for supporting long-term genetic sustainability and genetic diversity among dairy cattle populations should be prioritized in managing farm profitability. The present review provides a broad approach to current inbreeding-derived problems, identifying critical points to be solved and possible alternative strategies to control selection against homozygous haplotypes while maintaining sustained selection pressure. Moreover, this manuscript explores future perspectives, emphasizing theoretical applications and critical points, and strategies to avoid the adverse effects of inbreeding in dairy cattle. Finally, this review provides an overview of challenges that will soon require multidisciplinary approaches to managing dairy cattle populations, intending to combine increases in productive trait phenotypes with improvements in reproductive, health, welfare, linear conformation, and adaptability traits into the foreseeable future.

Performance of Swedish Warmblood fragile foal syndrome carriers and breeding prospects

Background

Warmblood fragile foal syndrome (WFFS) is a monogenetic defect caused by a recessive lethal missense point mutation in the procollagen‐lysine, 2‐oxoglutarate 5‐dioxygenase 1 gene (PLOD1, c.2032G>A). The majority of homozygous WFFS horses are aborted during gestation. Clinical signs of affected horses include fragile skin, skin and mucosa lacerations, hyperextension of the articulations, and hematomas. In spite of its harmful effect, a relatively high frequency of WFFS carriers has been found in Warmblood horses, suggesting a heterozygote advantage. Thus, in this study our aims were to: (1) estimate the frequency of WFFS carriers in the Swedish Warmblood breed (SWB), (2) estimate the effect of WFFS carrier genotype on performance traits in two SWB subpopulations bred for different disciplines, and (3) simulate the potential effects of balancing selection and different selection strategies on the frequency of carriers.

Methods

In total, 2288 SWB sport horses born between 1971 and 2020 were tested for the WFFS mutation and had estimated breeding values (EBV) for ten traditional evaluating and 50 linear descriptive traits.

Results

The frequency of WFFS carriers calculated from a pool of 511 randomly selected SWB horses born in 2017 was equal to 7.4% and ranged from 0.0 to 12.0% among the whole set of tested SWB horses, starting from 1971 till 2020. The effect of the WFFS carrier genotype was significant for several EBV mainly related to movements and dressage traits and especially for horses not bred for the show jumping discipline. Using simulation, we showed that balancing selection can maintain a recessive lethal allele in populations such as the SWB breed over generations and that the frequency is expected to slowly decrease in absence of balancing selection. Finally, we showed that selection against carrier sires can result in a more rapid decrease of the frequency of the mutant allele over time.

Conclusion

Further research is needed to confirm the apparent association between equine performance and the WFFS carrier genotype. Identification of such associations or new causative mutations for horse performance traits can serve as new tools in horse breeding to select for healthy, sustainable, and better performing horses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00693-4.

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

Integrated analysis of whole genome and transcriptome sequencing reveals a frameshift mutation associated with recessive embryonic lethality in Holstein cattle

Embryo loss is an important factor affecting fertility in dairy production. HH2 was identified as a haplotype on chromosome 1 associated with embryonic lethality in Holstein cattle. In the current study, both short- and long-read WGS was performed on four carriers and four non-carriers of HH2 to screen for variants in concordance with HH2 haplotype status. Sequence variation analysis revealed five putative functional variants of protein-coding genes, including a frameshift mutation (g.107172616delT) in intraflagellar transport protein 80 (IFT80) gene. Transcriptome analysis of whole blood indicated that no gene exhibited significantly differential expression or allele-specific expression between carriers and non-carriers in the candidate region. This evidence points to g.107172616delT as the highest priority causative mutation for HH2. Protein prediction reveals that the frameshift mutation results in a premature stop codon to reduce the peptide chain from 760 to 383 amino acids and greatly alters the structure and function of IFT80 protein. Our results demonstrate that the use of a combination of multiple high-throughput sequencing technologies is an efficient strategy to screen for the candidate causative mutations responsible for Mendelian traits, including genetic disorders.

Frequency of genotypic markers for genetic disorders, colour, polledness, and major genes in Blanco Orejinegro cattle

The purpose of this work was to establish DNA marker frequencies for genetic disorders, colour, horned/polled trait, and major genes of importance for productive and reproductive traits in Blanco Orejinegro (BON) cattle. The Blanco Orejinegro breed is a Colombian creole breed characterized by a white hair coat on black skin with black ears, black hair on the middle part of the legs, and absence of horns. We genotyped 420 animals of Colombia, 70 with the GGP-Bovine 150 K chip and 350 with the GGP-Bovine 50 K chip. The markers were associated with 50 genetic diseases, 52 major gene variants related to productive traits, and 12 variants related to coat coloration, presence of horns, and adaptation, selected from the information contained in the chips. Genotype frequencies were estimated using the R statistical program. Genetic disorder annotations were derived using the Online Mendelian Inheritance in Animals tool (OMIA) and the average inbreeding coefficient (F) (n = 7799) using the MTDFNRM program. Carriers were found for 16 of the genetic disorders evaluated but with low frequencies (0.24 to 2.46%); no homozygous animals were found for the disorders. Carriers were associated with disorders such as bovine leukocyte adhesion deficiency (BLAD), deficiency of uridine monophosphate synthase (DUMPS), syndactyly, and epidermolysis bullosa (EB). The F was 4.41%. Concerning the genes associated with colour (TYR, MC1R, and PMEL), alleles related to black pigmentation, the absence of horns (polledness), and slick coat (an adaptive trait) were highly frequent (> 81.90%). Markers associated with milk production and quality, yellow fat, and fertility showed variable frequencies, indicating selection potential. Allele frequency of genetic disorders in BON cattle was low, suggesting few genetic disorder problems, with syndactyly being the most frequent condition. The markers associated with colour and polledness were almost fixed, with a frequency at or near 100%. Production and reproduction markers showed variability for selection.

Identification of candidate genes on the basis of SNP by time-lagged heat stress interactions for milk production traits in German Holstein cattle

The aim of this study was to estimate genotype by time-lagged heat stress (HS) variance components as well as main and interaction SNP-marker effects for maternal HS during the last eight weeks of cow pregnancy, considering milk production traits recorded in the offspring generation. The HS indicator was the temperature humidity index (THI) for each week. A dummy variable with the code = 1 for the respective week for THI ≥ 60 indicated HS, otherwise, for no HS, the code = 0 was assigned. The dataset included test-day and lactation production traits from 14,188 genotyped first parity Holstein cows. After genotype quality control, 41,139 SNP markers remained for the genomic analyses. Genomic animal models without (model VC_nHS) and with in-utero HS effects (model VC_wHS) were applied to estimate variance components. Accordingly, for genome-wide associations, models GWA_nHS and GWA_wHS, respectively, were applied to estimate main and interaction SNP effects. Common genomic and residual variances for the same traits were very similar from models VC_nHS and VC_wHS. Genotype by HS interaction variances varied, depending on the week with in-utero HS. Among all traits, lactation milk yield with HS from week 5 displayed the largest proportion for interaction variances (0.07). For main effects from model GWA_wHS, 380 SNPs were suggestively associated with all production traits. For the SNP interaction effects from model GWA_wHS, we identified 31 suggestive SNPs, which were located in close distance to 62 potential candidate genes. The inferred candidate genes have various biological functions, including mechanisms of immune response, growth processes and disease resistance. Two biological processes excessively represented in the overrepresentation tests addressed lymphocyte and monocyte chemotaxis, ultimately affecting immune response. The modelling approach considering time-lagged genotype by HS interactions for production traits inferred physiological mechanisms being associated with health and immunity, enabling improvements in selection of robust animals.

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.

First report of cholesterol deficiency associated APOB mutation causing calf mortality in Indian Holstein Friesian population

Cholesterol deficiency is a new monogenic autosomal recessive inherited defect in Holstein cattle and is associated with chronic and recurrent diarrhea followed by death of calves during the first six months of life. We have successfully screened the causative mutation for cholesterol deficiency in indigenous, exotic and crossbred cattle of India. This is the first report about the presence of APOB mutation causing calf mortality in the Indian Holstein Friesian population. The mutation associated with cholesterol deficiency was not detected in any of indigenous breeds of cattle suggesting that the mutation could be specific to Holstein Friesian breed only. The occurrence of causal mutation for cholesterol deficiency in one of the Holstein Friesian bull under study (with carrier frequency of 1.67%) clearly indicates that mutation has already been transmitted into Indian Holstein-Friesian population. In our opinion it is the right time to take stringent actions in order to check its further spread. Further, this study accentuates the need for regular screening of existing as well as importing breeding bulls, especially of Holstein-Friesian lineage, for CDH genetic defect before introduction of bulls into genetic improvement programs. It is also highly recommended that the screening for CDH mutation must be made mandatory for all the semen producing/importing agencies under the minimum standards for production of bovine semen.

Evaluation of bull fertility in Italian Brown Swiss dairy cattle using cow field data

Dairy bull fertility is traditionally evaluated using semen production and quality traits; however, these attributes explain only part of the differences observed in fertility among bulls. Alternatively, bull fertility can be directly evaluated using cow field data. The main objective of this study was to investigate bull fertility in the Italian Brown Swiss dairy cattle population using confirmed pregnancy records. The data set included a total of 397,926 breeding records from 1,228 bulls and 129,858 lactating cows between first and fifth lactation from 2000 to 2019. We first evaluated cow pregnancy success, including factors related to the bull under evaluation, such as bull age, bull inbreeding, and AI organization, and factors associated with the cow that receives the dose of semen, including herd-year-season, cow age, parity, and milk yield. We then estimated sire conception rate using only factors related to the bull. Model predictive ability was evaluated using 10-fold cross-validation with 10 replicates. Interestingly, our analyses revealed that there is a substantial variation in conception rate among Brown Swiss bulls, with more than 20% conception rate difference between high-fertility and low-fertility bulls. We also showed that the prediction of bull fertility is feasible as our cross-validation analyses achieved predictive correlations equal to 0.30 for sire conception rate. Improving reproduction performance is one of the major challenges of the dairy industry worldwide, and for this, it is essential to have accurate predictions of service sire fertility. This study represents the foundation for the development of novel tools that will allow dairy producers, breeders, and artificial insemination companies to make enhanced management and selection decisions on Brown Swiss male fertility.

A variance component estimation approach to infer associations between Mendelian polledness and quantitative production and female fertility traits in German Simmental cattle

Background

Managing beneficial Mendelian characteristics in dairy cattle breeding programs implies that the correlated genetic effects are considered to avoid possible adverse effects in selection processes. The Mendelian trait polledness in cattle is traditionally associated with the belief that the polled locus has unfavorable effects on breeding goal traits. This may be due to the inferior breeding values of former polled bulls and cows in cattle breeds, such as German Simmental, or to pleiotropic or linkage effects of the polled locus.

Methods

We focused on a variance component estimation approach that uses a marker-based numerator relationship matrix reflecting gametic relationships at the polled locus to test for direct pleiotropic or linked quantitative trait loci (QTL) effects of the polled locus on relevant traits. We applied the approach to performance, health, and female fertility traits in German Simmental cattle.

Results

Our results showed no evidence for any pleiotropic QTL effects of the polled locus on test-day production traits milk yield and fat percentage, on the mastitis indicator ‘somatic cell score’, and on several female fertility traits, i.e. 56 days non return rate, days open and days to first service. We detected a significant and unfavorable QTL effect accounting for 6.6% of the genetic variance for protein percentage only.

Conclusions

Pleiotropy does not explain the lower breeding values and phenotypic inferiority of polled German Simmental sires and cows relative to the horned population in the breed. Thus, intensified selection in the polled population will contribute to increased selection response in breeding goal traits and genetic merit and will narrow the deficit in breeding values for production traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00652-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.

Identification of homozygous haplotypes carrying putative recessive lethal mutations that compromise fertility traits in French Lacaune dairy sheep

Background

Homozygous recessive deleterious mutations can cause embryo/fetal or neonatal lethality, or genetic defects that affect female fertility and animal welfare. In livestock populations under selection, the frequency of such lethal mutations may increase due to inbreeding, genetic drift, and/or the positive pleiotropic effects of heterozygous carriers on selected traits.

Results

By scanning the genome of 19,102 Lacaune sheep using 50 k single nucleotide polymorphism (SNP) phased genotypes and pedigree data, we identified 11 Lacaune deficient homozygous haplotypes (LDHH1 to LDHH11) showing a highly significant deficit of homozygous animals ranging from 79 to 100%. These haplotypes located on chromosomes 3, 4, 13, 17 and 18, spanned regions from 1.2 to 3.0 Mb long with a frequency of heterozygous carriers between 3.7 and 12.1%. When we compared at-risk matings (between carrier rams and daughters of carrier rams) and safe matings, seven of the 11 haplotypes were associated with a significant alteration of two fertility traits, a reduced success of artificial insemination (LDHH1, 2, 8 and 9), and/or an increased stillbirth rate (LDHH3, 6, 8, 9, and 10). The 11 haplotypes were also tested for a putative selective advantage of heterozygous carrier rams based on their daughter yield deviation for six dairy traits (milk, fat and protein yields, fat and protein contents and lactation somatic cell score). LDHH1, 3, 4, 5, 7, 9 and 11 were associated with positive effects on at least one selected dairy trait, in particular milk yield. For each haplotype, the most probable candidate genes were identified based on their roles in lethality of mouse knock-out models and in mammalian genetic disorders.

Conclusions

Based on a reverse genetic strategy, we identified at least 11 haplotypes with homozygous deficiency segregating in French Lacaune dairy sheep. This strategy represents a first tool to limit at-risk matings in the Lacaune dairy selection scheme. We assume that most of the identified LDHH are in strong linkage disequilibrium with a recessive lethal mutation that affects embryonic or juvenile survival in sheep but is yet to be identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00634-1.

Alternate PCR assays for screening of JH1 mutation associated with embryonic death in Jersey cattle

Jersey haplotype (JH) 1, a stop-gain lethal mutation in the CWC15 gene, causes embryonic losses in Jersey cattle. Two PCR based assays using Amplification Refractory Mutation System (T-ARMS-PCR) and restriction fragment length polymorphism (PCR-RFLP) were developed for screening of the JH1 in cattle. During the screening, seven among 30 Indian Jersey bulls were identified as carriers of the mutant JH1 allele, the first time in the country. These PCR assays are economical, rapid and accurate; and can be used separately or in combination for screening and cross-validation of the JH1 carriers in Jersey cattle.

Cholesterol deficiency haplotype frequency and its impact on milk production and milk cholesterol content in Canadian Holstein cows

An autosomal recessive defect caused by a loss of function mutation in the bovine apolipoprotein B (APOB) gene causes cholesterol deficiency with strong effects on calf survival and development. This study examined the distribution of cholesterol deficiency haplotype (CDH) in the APOB gene in Canadian Holstein cows and its impacts on milk production and milk cholesterol content. The CDH in 917 cows from 29 herds was determined via polymerase chain reaction (PCR). A longer insertion of approximately 7 kb potentially characterizing the mutation was amplified. The average carrier frequency of herds in Quebec based on PCR results was 14.6%. Cholesterol deficiency haplotype status had non-significant (P > 0.05) effects on milk cholesterol content, other milk components, and milk yield.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach

Deviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific- and unspecific-parent TRD) and genotypic (additive- and dominance-TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent-specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.

Cholesterol deficiency carriers have lowered serum cholesterol and perform well at an elite cattle show

•

Cholesterol deficiency carriers had lower serum cholesterol than noncarriers.

•

Cholesterol deficiency carriers had faster milking speeds than noncarriers.

•

Cholesterol deficiency carriers had lower body weights than noncarriers.

•

Cows that carry cholesterol deficiency tend to have an advantage in the show ring.

Cholesterol deficiency (CD) in Holstein cattle results from an insertion in the APOB gene, which traces to the bull Maughlin Storm and was spread to elite show cattle largely through his grandson Braedale Goldwyn. Calves that are homozygous for CD can neither transport dietary fat from the gut nor synthesize low-density lipoproteins, and they perish before 6 mo of age. The mutation may have a partially dominant effect, with lowered cholesterol in heterozygous animals. Our aims were to evaluate serum cholesterol in heterozygous animals, evaluate daily phenotypic records of lactating CD carriers, and determine whether CD carrier status is associated with show ring success. Blood was collected on a single date in May 2018 from 12 carriers and 14 controls and 1 yr later from 17 carriers (7 repeated for 2018) and 27 controls (6 repeated from 2018); animals ranged in age from 14 to 84 mo and varied in stage of lactation and pregnancy status. Serum samples were analyzed for cholesterol concentration (mg/dL) and results were analyzed using mixed models that included linear and quadratic effects of days in milk and days carried calf. Daily phenotypic records of milk yield, milking time, milk conductivity, activity, and body weight were compared between carriers and controls. Additionally, the CD carrier status of cows that were sired by known CD carriers and had available genotypes was recorded. Cows placing in the top 5 and top 10 of their classes at World Dairy Expo between the years of 2006 and 2019 were included in our analysis. Chi-squared tests for equal proportions were evaluated to determine whether the following were associated with CD status: placing in the top 10 with repeated placings from the same cow included, placing in the top 5 with repeated placings, top 10 cows considering individuals only once, and top 5 cows considering individuals only once. Control animals (noncarriers) had significantly higher serum cholesterol (116.21 ± 6.45 mg/dL) than CD carriers (90.15 ± 7.49 mg/dL). There was no difference between carriers and controls in milk yield, milk conductivity, or activity. However, CD carriers had shorter milking times and lower body weights than their herdmates. Further, there were more CD carriers in the top 5 and more carriers tended to place in the top 10 at World Dairy Expo when repeated placings were included. These data support observations that CD is partially dominant and indicate that CD carriers display other phenotypic differences from their herdmates that may give them a small advantage at elite cattle shows.

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.

Symposium review: Exploiting homozygosity in the era of genomics-Selection, inbreeding, and mating programs

The advent of genomic selection paved the way for an unprecedented acceleration in genetic progress. The increased ability to select superior individuals has been coupled with a drastic reduction in the generation interval for most dairy populations, representing both an opportunity and a challenge. Homozygosity is now rapidly accumulating in dairy populations. Currently, inbreeding depression is managed mostly by culling at the farm level and by controlling the overall accumulation of homozygosity at the population level. A better understanding of how homozygosity and recessive load are related will guarantee continued genetic improvement while curtailing the accumulation of harmful recessives and maintaining enough genetic variability to ensure the possibility of selection in the face of changing environmental conditions. In this review, we present a snapshot of the current dairy selection structure as it relates to response to selection and accumulation of homozygosity, briefly outline the main approaches currently used to manage inbreeding and overall variability, and present some approaches that can be used in the short term to control accumulation of harmful recessives while maintaining sustained selection pressure.

Reconstitution and modernization of lost Holstein male lineages using samples from a gene bank

More than 99% of all known Holstein artificial insemination (AI) bulls in the United States can be traced through their male lineage to just 2 bulls born in the 1950s, and all Holstein bulls can be traced back to 2 bulls born in the late 1800s. As the Y chromosome is passed exclusively from sire to son, this suggests that variation is limited for much of the Y chromosome. Two additional male lineages that are separate from modern lineages before 1890 were present at the start of the AI era and had semen available from the USDA National Animal Germplasm Program (Fort Collins, CO). Semen from representatives of those lineages were used for in vitro embryo production by mating to elite modern genetic females, resulting in the birth of 7 bulls and 8 heifers. Genomic evaluation of the bulls suggested that lineages from the beginning of the AI era could be reconstituted to breed average for total economic merit in 1 generation when mated to elite females due to high genetic merit for fertility, near-average genetic merit for fat and protein yield, and below-average genetic merit for udder and physical conformation. Semen from the bulls is commercially available to facilitate Y chromosome research and efforts to restore lost genetic diversity.

Genetic and genomic analysis of long insemination interval in Israeli dairy cattle as an indicator of early abortions

One of the causes of observed low fertility is embryo loss after fertilization. Previous findings suggested that more than half of fertilizations result in embryo loss before pregnancy is detected. We proposed reinsemination between 49 and 100 d after the first insemination as an indicator trait for early abortion (EA) in dairy cattle based on the mean estrus interval of 21 d. This trait was compared with conception rate from first insemination and conception status, computed as the inverse of the number of inseminations to conception. Animal model variance components were estimated by REML, including parents and grandparents of cows with records. First-parity heritability for first insemination conception rate was 3%. In the multitrait analysis of parities 1 to 3 for putative EA, heritabilities ranged from 8.9% for first parity to 10.4% for second parity. All genetic correlations were >0.9, whereas all environmental correlations were <0.12. The variance component for the service sire effect for putative EA rate was less than half the variance component for conception rate. Thus, genetic control of the 2 traits is clearly different, and analysis of EA rate by a single-trait animal model is justified. Genetic evaluation for putative EA was computed using this model, including all first- through third-parity cows with freshening dates from January 1, 1985, through December 31, 2016, that either became pregnant on first insemination or were reinseminated between 49 and 100 d after the first insemination. All known parents and grandparents of cows with records were included in the analysis. The regression of the breeding value for non-abortion rate on the cows' birth year was 0.083%/yr. The genetic correlation between first-parity EA and conception status was 0.995. The genetic correlations between first-parity EA and milk, fat, and protein production were all negative, whereas the genetic correlation between EA and herd life was 0.33. Inclusion of putative EA in the selection index instead of conception status resulted in 10 to 20% greater genetic gain for both fertility traits. In a genome-wide association study based on 1,200 dairy bulls with reliabilities >50% for abortion rate genotyped for 41,000 markers, 6 markers were found with nominal probabilities of <10^-12 to reject the null hypothesis of no effect on EA rate. The markers with the lowest probabilities for EA rate were also included among the markers with the lowest probabilities for female fertility, but not vice versa. The marker explaining the most variance for abortion rate is located within the ABCA9 gene, which is found within an ATP-binding cassette (ABC) genes cluster. The ABC family is the major class of primary active transporters in the placenta.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires

Background

Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers.

Methods

The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries.

Results

A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles’ QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG “Circadian rhythm” and “Neurotrophin signaling pathway” were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones.

Conclusions

The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

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.

The APOB loss-of-function mutation of Holstein dairy cattle does not cause a deficiency of cholesterol but decreases the capacity for cholesterol transport in circulation

The loss-of-function mutation of the apolipoprotein (APO) B gene (APOB) in Holstein cattle accounts for increased losses in calves that are homozygous for this mutation. Heterozygous carriers of the APOB mutation are clinically healthy but show decreased concentrations of plasma cholesterol and lipoproteins. So far, the metabolic effects of the mutation have only been investigated in heterozygous calves, bulls, and nonlactating females. In high-yielding dairy cows, a marked decrease in cholesterol concentration in plasma during early lactation is part of the usual metabolic changes. Given the essential role of cholesterol in fatty acid and lipid metabolism, a specific effect of the APOB mutation on metabolism and performance in dairy cows is expected. Therefore, the aim of the present study was to investigate the effects of different APOB genotypes on metabolic parameters, hepatic metabolism, and lactation and reproductive performance. Twenty pairs of full siblings with similar age, performance, and calving were investigated. Both animals of each pair were kept on the same farm and consisted of a heterozygous carrier (CDC) and a noncarrier (CDF) of the APOB mutation associated with cholesterol deficiency. Blood samples were taken in early (25.5 ± 4.7 d in milk) and mid lactation (158.2 ± 11.1 d in milk; mean ± SD), and analyzed for nonesterified fatty acids, β-hydroxybutyrate, glucose, insulin-like growth factor-1, aspartate aminotransferase and gamma-glutamyltransferase activity, total cholesterol, free cholesterol, triacylglycerols, high density lipoprotein-cholesterol, and phospholipids. The evaluation of milk production, milk gross composition, and lactation persistency was based on official Dairy Herd Improvement Association recordings. Cholesterol and lipoprotein concentrations were lower in CDC cows than in CDF cows in early and mid lactation. Metabolic parameters, triacylglycerol concentration in plasma, and lactation and reproductive performance did not differ between CDC cows and CDF cows. The low cholesterol concentrations associated with the APOB mutation in heterozygous carriers are not because of a primary deficiency of cholesterol at a cellular level, as the term "cholesterol deficiency" suggests, but rather a consequence of reduced capacity for its transport in circulation. Overall, the data of the present study suggest that, despite the presence of the APOB mutation, cholesterol is not limiting for animals' metabolic adaptation and performance in heterozygous Holstein cows.

Association of lipid-related genes implicated in conceptus elongation with female fertility traits in dairy cattle

Elongation of the preimplantation conceptus is a requirement for pregnancy success in ruminants, and failures in this process are highly associated with subfertility in dairy cattle. Identifying genetic markers that are related to early conceptus development and survival and utilizing these markers in selective breeding can improve the reproductive efficiency of dairy herds. Here, we evaluated the association of 1,679 SNP markers within or close to 183 candidate genes involved in lipid metabolism of the elongating conceptus with different fertility traits in US Holstein cattle. A total of 27,371 bulls with predicted transmitting ability records for daughter pregnancy rate, cow conception rate, and heifer conception rate were used as the discovery population. The associations found in the discovery population were validated using 2 female populations (1,122 heifers and 2,138 lactating cows) each with 4 fertility traits, including success to first insemination, number of services per conception, age at first conception for heifers, or days open for cows. Marker effects were estimated using a linear mixed model with SNP genotype as a linear covariate and a random polygenic effect. After multiple testing correction, 39 SNP flagging 27 candidate genes were associated with at least one fertility trait in the discovery population. Of these 39 markers, 3 SNP were validated in the heifer population and 4 SNP were validated in the cow population. The 3 SNP validated in heifers are located within or near genes CAT, MYOF, and RBP4, and the 4 SNP validated in lactating cows are located within or close to genes CHKA, GNAI1, and HMOX2. These validated genes seem to be relevant for reducing pregnancy losses, and the SNP within these genes are excellent candidates for inclusion in genomic tests to improve reproductive performance in dairy 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.

Short communication: Phenotypic and genetic effects of the polled haplotype on yield, longevity, and fertility in US Brown Swiss, Holstein, and Jersey cattle

Phenotypes from the December 2018 US national genetic evaluations were used to compute effects of the polled haplotype in US Brown Swiss (BS), Holstein (HO), and Jersey (JE) cattle on milk, fat, and protein yields, somatic cell score, single-trait productive life, daughter pregnancy rate, heifer conception rate, and cow conception rate. Lactation records pre-adjusted for nongenetic factors and direct genomic values were used to estimate phenotypic and genetic effects of the polled haplotype, respectively. No phenotypic or direct genomic values effects were different from zero for any trait in any breed. Genomic PTA (gPTA) for the lifetime net merit (NM$) selection index of bulls born since January 1, 2012, that received a marketing code from the National Association of Animal Breeders (Madison, WI), and cows born on or after January 1, 2015, were compared to determine whether there was a systematic benefit to polled or horned genetics. Horned bulls had the highest average gPTA for NM$ in all 3 breeds, but that difference was significant only in HO and JE (HO: 615.4 ± 1.9, JE: 402.3 ± 3.4). Homozygous polled BS cows had significantly higher average gPTA for NM$ than their heterozygous polled or horned contemporaries (PP = 261.4 ± 43.5, Pp = 166.1 ± 13.7, pp = 174.1 ± 1.8), but the sample size was very small (n = 9). In HO and JE, horned cows had higher gPTA for NM$ (HO = 378.3 ± 0.2, JE = 283.3 ± 0.3). Selection for polled cattle should not have a detrimental effect on yield, fertility, or longevity, but these differences show that, in the short term, selection for polled over horned cattle will result in lower rates of genetic gain.

The cholesterol-deficiency associated mutation in APOB segregates at low frequency in Chinese Holstein cattle

Cholesterol deficiency (CD) is a recessive genetic defect in Holstein cattle, which affects the rearing success of calves. We investigated whether CD has been spread into the People’s Republic of China. A total of 138 Holstein bulls and 90 cows were screened for CD using polymerase chain reaction amplification and then seven bulls (5.07%) and one cow (1.11%) were identified as CD carriers. Pedigree analyses revealed that all seven CD carrier bulls were traced back to MAUGHLIN STORM and CD may have spread into the People’s Republic of China through his daughter, Braedale Baler Twine. Our results demonstrated that the CD defect exists in the Chinese Holstein population.

Identification of quantitative trait loci for milk protein percentage in Murrah buffaloes

Milk protein is an important constituent of milk in buffaloes and is moderately heritable. The milk protein percentage varies significantly between breeds/herds/species. Buffaloes can be selected for higher milk protein percentage and this paper provides QTLs for marker assisted selection in buffaloes. The milk protein percentage records on 2,028 daughters belonging to 12 half sib families were analyzed for the identification of QTLs on 8 chromosomes in buffaloes using chromosome scans. The single marker analysis revealed 74 markers to be associated with milk protein percentage in 12 sire families. When common markers were removed from the analysis, 51 markers remained. The Interval mapping using R/qtl identified 69 QTLs in 12 half sib families on 8 chromosomes of buffalo. The meta QTL analysis defined 25 consensus QTL regions in buffaloes for milk protein percentage. Most of the QTLs identified have been reported for cattle however few new chromosomal locations were also identified to be associated with milk protein percentage in buffaloes. Comparative genomics revealed 1117 genes underlying the QTL regions associated with milk protein percentage. Among these, 109 genes were directly associated with protein metabolism. The protein-protein interaction among the genes and gene ontology analysis and pathways have been identified. These 109 genes have potential to be candidate genes for milk protein percentage in buffaloes.

Removal of alleles by genome editing (RAGE) against deleterious load

BACKGROUND

In this paper, we simulate deleterious load in an animal breeding program, and compare the efficiency of genome editing and selection for decreasing it. Deleterious variants can be identified by bioinformatics screening methods that use sequence conservation and biological prior information about protein function. However, once deleterious variants have been identified, how can they be used in breeding?

RESULTS

We simulated a closed animal breeding population that is subject to both natural selection against deleterious load and artificial selection for a quantitative trait representing the breeding goal. Deleterious load was polygenic and was due to either codominant or recessive variants. We compared strategies for removal of deleterious alleles by genome editing (RAGE) to selection against carriers. When deleterious variants were codominant, the best strategy for prioritizing variants was to prioritize low-frequency variants. When deleterious variants were recessive, the best strategy was to prioritize variants with an intermediate frequency. Selection against carriers was inefficient when variants were codominant, but comparable to editing one variant per sire when variants were recessive.

CONCLUSIONS

Genome editing of deleterious alleles reduces deleterious load, but requires the simultaneous editing of multiple deleterious variants in the same sire to be effective when deleterious variants are recessive. In the short term, selection against carriers is a possible alternative to genome editing when variants are recessive. Our results suggest that, in the future, there is the potential to use RAGE against deleterious load in animal breeding.

Analysis of a large dataset reveals haplotypes carrying putatively recessive lethal and semi-lethal alleles with pleiotropic effects on economically important traits in beef cattle

BACKGROUND

In livestock, deleterious recessive alleles can result in reduced economic performance of homozygous individuals in multiple ways, e.g. early embryonic death, death soon after birth, or semi-lethality with incomplete penetrance causing reduced viability. While death is an easy phenotype to score, reduced viability is not as easy to identify. However, it can sometimes be observed as reduced conception rates, longer calving intervals, or lower survival for live born animals.

METHODS

In this paper, we searched for haplotypes that carry putatively recessive lethal or semi-lethal alleles in 132,725 genotyped Irish beef cattle from five breeds: Aberdeen Angus, Charolais, Hereford, Limousin, and Simmental. We phased the genotypes in sliding windows along the genome and used five tests to identify haplotypes with absence of or reduced homozygosity. Then, we associated the identified haplotypes with 44,351 insemination records that indicated early embryonic death, and postnatal survival records. Finally, we assessed haplotype pleiotropy by estimating substitution effects on estimates of breeding value for 15 economically important traits in beef production.

RESULTS

We found support for one haplotype that carries a putatively recessive lethal (chromosome 16 in Simmental) and two haplotypes that carry semi-lethal alleles (chromosome 14 in Aberdeen Angus and chromosome 19 in Charolais), with population frequencies of 8.8, 15.2, and 14.4%, respectively. These three haplotypes showed pleiotropic effects on economically important traits for beef production. Their allele substitution effects are €2.30, €3.42, and €1.47 for the terminal index and €1.03, - €3.11, and - €0.88 for the replacement index, where the standard deviations for the terminal index are €22.52, €18.65, and €22.70 and for the replacement index they are €31.35, €29.82, and €35.79. We identified ZFAT as the candidate gene for semi-lethality in Aberdeen Angus, several candidate genes for the lethal Simmental haplotype, and no candidate genes for the semi-lethal Charolais haplotype.

CONCLUSIONS

We analysed genotype, reproduction, survival, and production data to detect haplotypes that carry putatively recessive lethal or semi-lethal alleles in Irish beef cattle and identified one lethal and two semi-lethal haplotypes, which have pleiotropic effects on economically important traits in beef production.

Diversity of copy number variation in the worldwide goat population

Goats (Capra hircus) are an important farm animal species. Copy number variation (CNV) represents a major source of genomic structural variation. We investigated the diversity of CNV distribution in goats using CaprineSNP50 genotyping data generated by the ADAPTmap Project. We identified 6286 putative CNVs in 1023 samples from 50 goat breeds using PennCNV. These CNVs were merged into 978 CNV regions, spanning ~262 Mb of total length and corresponding to ~8.96% of the goat genome. We then divided the samples into six subgroups per geographic distribution and constructed a comparative CNV map. Our results revealed a population differentiation in CNV across different geographical areas, including Western Asia, Eastern Mediterranean, Alpine & Northern Europe, Madagascar, Northwestern Africa, and Southeastern Africa groups. The results of a cluster heatmap analysis based on the CNV count per individual across different groups was generally consistent with the one generated from the SNP data, likely reflecting the population history of different goat breeds. We sought to determine the gene content of these CNV events and found several important CNV-overlapping genes (e.g. EDNRA, ADAMTS20, ASIP, KDM5B, ADAM8, DGAT1, CHRNB1, CLCN7, and EXOSC4), which are involved in local adaptations such as coat color, muscle development, metabolic processes, osteopetrosis, and embryonic development. Therefore, this research generated an extensive CNV map in the worldwide population of goat, which offers novel insight into the goat genome and its functional annotation.

Screening genetic diseases prevalence in Braunvieh cattle

Heritable abnormalities can cause a reduction in productive performance, structural defects, or death of the animal. There are reports of hereditary abnormalities in Braunvieh cattle from several countries, but no evidence was found on their existence in Mexico. In this study, 28 genes associated with hereditary diseases were screened with the GGP-LD 30K array (GeneSeek^®) in 300 Mexican registered Braunvieh animals. Allelic frequencies of the markers associated with illness were obtained for the following: citrullinaemia, spinal dysmyelination, spinal muscular atrophy, Brows Swiss fertility haplotype 2, congenital muscular dystonia, epidermolysis bullosa, Pompes, maple syrup urine, syndactyly, Weaver syndrome, crooked tail, deficiency of uridine monophosphate synthase, hypotrichosis, Marfan syndrome, and weak calf syndrome. The allelic frequency values were low for all the analysed loci (from 0.0015 to 0.0110), with exception of syndactyly (0.4145). Although homozygous animals for these genetic conditions were detected, no physical or physiological abnormalities associated with the clinical form of the diseases were observed in the sampled animals. Markers associated with a crooked tail, deficiency of uridine monophosphate synthase, hypotrichosis, Marfan syndrome, and weak calf syndrome were absent. The studied Mexican Braunvieh population does not present clinical or subclinical effects for ten diseases in homozygous animals. However, since the assessed animals are considered as breeding stock, the monitoring of carrier animals might be periodically necessary.

Review: Genomics of bull fertility

Fertility is one of the most economically important traits in both beef and dairy cattle production; however, only female fertility is typically subjected to selection. Male and female fertility have only a small positive genetic correlation which is likely due to the existence of a relatively small number of genetic variants within each breed that cause embryonic and developmental losses. Genomic tools have been developed that allow the identification of lethal recessive loci based upon marker haplotypes. Selection against haplotypes harbouring lethal alleles in conjunction with selection to improve female fertility will result in an improvement in male fertility. Genomic selection has resulted in a two to fourfold increase in the rate of genetic improvement of most dairy traits in US Holstein cattle, including female fertility. Considering the rapidly increasing rate of adoption of high-throughput single nucleotide polymorphism genotyping in both the US dairy and beef industries, genomic selection should be the most effective of all currently available approaches to improve male fertility. However, male fertility phenotypes are not routinely recorded in natural service mating systems and when artificial insemination is used, semen doses may be titrated to lower post-thaw progressively motile sperm numbers for high-merit and high-demand bulls. Standardization of sperm dosages across bull studs for semen distributed from young bulls would allow the capture of sire conception rate phenotypes for young bulls that could be used to generate predictions of genetic merit for male fertility in both males and females. These data would allow genomic selection to be implemented for male fertility in addition to female fertility within the US dairy industry. While the rate of use of artificial insemination is much lower within the US beef industry, the adoption of sexed semen in the dairy industry has allowed dairy herds to select cows from which heifer replacements are produced and cows that are used to produce terminal crossbred bull calves sired by beef breed bulls. Capture of sire conception rate phenotypes in dairy herds utilizing sexed semen will contribute data enabling genomic selection for male fertility in beef cattle breeds. As the commercial sector of the beef industry increasingly adopts fixed-time artificial insemination, sire conception rate phenotypes can be captured to facilitate the development of estimates of genetic merit for male fertility within US beef breeds.