Genome-wide mapping and estimation of inbreeding depression of semen quality traits in a cattle population

Genomic analysis of conservation status, population structure and admixture in local Czech and Slovak dairy goat breeds

While dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 years due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair and Czech Landrace) and one Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, ROH-based inbreeding and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters equal to 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds.

Inbreeding depression for producer-recorded udder, metabolic, and reproductive diseases in US dairy cattle

This study leveraged a growing dataset of producer-recorded phenotypes for mastitis, reproductive diseases (metritis and retained placenta), and metabolic diseases (ketosis, milk fever, and displaced abomasum) to investigate the potential presence of inbreeding depression for these disease traits. Phenotypic, pedigree, and genomic information were obtained for 354,043 and 68,292 US Holstein and Jersey cows, respectively. Total inbreeding coefficients were calculated using both pedigree and genomic information; the latter included inbreeding estimates obtained using a genomic relationship matrix and runs of homozygosity. We also generated inbreeding coefficients based on the generational inbreeding for recent and old pedigree inbreeding, for different run-of-homozygosity length classes, and for recent and old homozygous-by-descent segment-based inbreeding. Estimates on the liability scale revealed significant evidence of inbreeding depression for reproductive-disease traits, with an increase in total pedigree and genomic inbreeding showing a notable effect for recent inbreeding. However, we found inconsistent evidence for inbreeding depression for mastitis or any metabolic diseases. Notably, in Holsteins, the probability of developing displaced abomasum decreased with inbreeding, particularly for older inbreeding. Estimates of disease probability for cows with low, average, and high inbreeding levels did not significantly differ across any inbreeding coefficient and trait combination, indicating that although inbreeding may affect disease incidence, it likely plays a smaller role compared with management and environmental factors.

Impact of inbreeding on production, fertility, and health traits in German Holstein dairy cattle utilizing various inbreeding estimators

In dairy cattle production, it is important to understand how inbreeding affects production, fertility, and health traits. However, there is still limited use of genomic information to estimate inbreeding, despite advancements in genotyping technologies. To address this gap, we investigated the impact of inbreeding on German Holstein dairy cattle using both pedigree-based and genomic-based inbreeding estimators. We employed one method based on pedigree information (Fped) together with 6 genomic-based methods, including 3 GCTA estimators (Fhat1, Fhat2, Fhat3), VanRaden's first method (FVR1, with observed allele frequencies and FVR0.5 when allele frequencies are set to 0.5), and one based on runs of homozygosity (Froh). Data from 24,489 cows with both phenotypes and genotypes were used, with a pedigree including 232,780 animals born between 1970 and 2018. We analyzed the effects of inbreeding depression on production, fertility, and health traits separately, using single-trait linear animal models as well as threshold models to account for the binary nature of the health traits. For the health traits, we transformed solutions from the liability scale to a probability scale for easier interpretation. Our results showed that the mean inbreeding coefficients from all estimators ranged from -0.003 to 0.243, with negative values observed for most genomic-based methods. We found out that a 1% increase in inbreeding caused a depression ranging from 25.94 kg (Fhat1) to 40.62 kg (Fhat3), 1.18 kg (Fhat2) to 1.70 kg (Fhat3), 0.90 kg (Fhat2) to 1.45 kg (Froh and Fhat3), 0.19 (Fped) to 0.34 d (Fhat3) for 305 d milk yield, fat, protein, and calving interval respectively. The health traits showed slight changes when inbreeding was gradually increased by 5% with digital dermatitis showing rather a contrasting trend to that of mastitis which reduced the more an animal was inbred. Overall, our study highlights the importance of considering both pedigree-based and genomic-based inbreeding estimators when assessing the impact on inbreeding, emphasizing that not all inbreeding is harmful.

Identification of runs of homozygosity associated with male fertility in Italian Brown Swiss cattle

Intensive selection for improved productivity has been accompanied by an increase in inbreeding rates and a reduction in genetic diversity. The increase in inbreeding tends to impact performance, especially fitness-related traits such as male fertility. Inbreeding can be monitored using runs of homozygosity (ROH), defined as contiguous lengths of homozygous genotypes observed in an individual's chromosome. The goal of this study was to evaluate the presence of ROH in Italian Brown Swiss cattle and assess its association with bull fertility. First, we evaluated the association between ROH and male fertility using 1,102 Italian Brown Swiss bulls with sire conception rate records and 572 K SNPs spanning the entire genome. Second, we split the entire population into 100 high-fertility and 100 low-fertility bulls to investigate the potential enrichment of ROH segments in the low-fertility group. Finally, we mapped the significant ROH regions to the bovine genome to identify candidate genes associated with sperm biology and male fertility. Notably, there was a negative association between bull fertility and the amount of homozygosity. Four different ROH regions located in chromosomes 6, 10, 11, and 24 were significantly overrepresented in low-fertility bulls (Fisher's exact test, p-value <0.01). Remarkably, these four genomic regions harbor many genes such as WDR19, RPL9, LIAS, UBE2K, DPF3, 5S-rRNA, 7SK, U6, and WDR7 that are related to sperm biology and male fertility. Overall, our findings suggest that inbreeding and increased homozygosity have a negative impact on male fertility in Italian Brown Swiss cattle. The quantification of ROH can contribute to minimizing the inbreeding rate and avoid its negative effect on fitness-related traits, such as male fertility.

Comparing pedigree and genomic inbreeding coefficients, and inbreeding depression of reproductive traits in Japanese Black cattle

BACKGROUND

Pedigree-based inbreeding coefficients have been generally included in statistical models for genetic evaluation of Japanese Black cattle. The use of genomic data is expected to provide precise assessment of inbreeding level and depression. Recently, many measures have been used for genome-based inbreeding coefficients; however, with no consensus on which is the most appropriate. Therefore, we compared the pedigree- ([Formula: see text]) and multiple genome-based inbreeding coefficients, which were calculated from the genomic relationship matrix with observed allele frequencies ([Formula: see text]), correlation between uniting gametes ([Formula: see text]), the observed vs expected number of homozygous genotypes ([Formula: see text]), runs of homozygosity (ROH) segments ([Formula: see text]) and heterozygosity by descent segments ([Formula: see text]). We quantified inbreeding depression from estimating regression coefficients of inbreeding coefficients on three reproductive traits: age at first calving (AFC), calving difficulty (CD) and gestation length (GL) in Japanese Black cattle.

RESULTS

The highest correlations with [Formula: see text] were for [Formula: see text] (0.86) and [Formula: see text] (0.85) whereas [Formula: see text] and [Formula: see text] provided weak correlations with [Formula: see text], with range 0.33-0.55. Except for [Formula: see text] and [Formula: see text], there were strong correlations among genome-based inbreeding coefficients ([Formula: see text] 0.94). The estimates of regression coefficients of inbreeding depression for [Formula: see text] was 2.1 for AFC, 0.63 for CD and -1.21 for GL, respectively, but [Formula: see text] had no significant effects on all traits. Genome-based inbreeding coefficients provided larger effects on all reproductive traits than [Formula: see text]. In particular, for CD, all estimated regression coefficients for genome-based inbreeding coefficients were significant, and for GL, that for [Formula: see text] had a significant.. Although there were no significant effects when using overall genome-level inbreeding coefficients for AFC and GL, [Formula: see text] provided significant effects at chromosomal level in four chromosomes for AFC, three chromosomes for CD, and two chromosomes for GL. In addition, similar results were obtained for [Formula: see text].

CONCLUSIONS

Genome-based inbreeding coefficients can capture more phenotypic variation than [Formula: see text]. In particular, [Formula: see text] and [Formula: see text] can be considered good estimators for quantifying inbreeding level and identifying inbreeding depression at the chromosome level. These findings might improve the quantification of inbreeding and breeding programs using genome-based inbreeding coefficients.

Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls

BACKGROUND

The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., F_ROH > 1 Mb) and excess of SNP homozygosity (F_SNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants.

RESULTS

Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in F_ROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting F_ROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility.

CONCLUSIONS

Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

Review: Genetic mutations affecting bull fertility

Cattle are a well-suited "model organism" to study the genetic underpinnings of variation in male reproductive performance. The adoption of artificial insemination and genomic prediction in many cattle breeds provide access to microarray-derived genotypes and repeated measurements for semen quality and insemination success in several thousand bulls. Similar-sized mapping cohorts with phenotypes for male fertility are not available for most other species precluding powerful association testing. The repeated measurements of the artificial insemination bulls' semen quality enable the differentiation between transient and biologically relevant trait fluctuations, and thus, are an ideal source of phenotypes for variance components estimation and genome-wide association testing. Genome-wide case-control association testing involving bulls with either aberrant sperm quality or low insemination success revealed several causal recessive loss-of-function alleles underpinning monogenic reproductive disorders. These variants are routinely monitored with customised genotyping arrays in the male selection candidates to avoid the use of subfertile or infertile bulls for artificial insemination and natural service. Genome-wide association studies with quantitative measurements of semen quality and insemination success revealed quantitative trait loci for male fertility, but the underlying causal variants remain largely unknown. Moreover, these loci explain only a small part of the heritability of male fertility. Integrating genome-wide association studies with gene expression and other omics data from male reproductive tissues is required for the fine-mapping of candidate causal variants underlying variation in male reproductive performance in cattle.

A comparison of marker-based estimators of inbreeding and inbreeding depression

BACKGROUND

The availability of genome-wide marker data allows estimation of inbreeding coefficients (F, the probability of identity-by-descent, IBD) and, in turn, estimation of the rate of inbreeding depression (ΔID). We investigated, by computer simulations, the accuracy of the most popular estimators of inbreeding based on molecular markers when computing F and ΔID in populations under random mating, equalization of parental contributions, and artificially selected populations. We assessed estimators described by Li and Horvitz (F_LH1 and F_LH2), VanRaden (F_VR1 and F_VR2), Yang and colleagues (F_YA1 and F_YA2), marker homozygosity (F_HOM), runs of homozygosity (F_ROH) and estimates based on pedigree (F_PED) in comparison with estimates obtained from IBD measures (F_IBD).

RESULTS

If the allele frequencies of a base population taken as a reference for the computation of inbreeding are known, all estimators based on marker allele frequencies are highly correlated with F_IBD and provide accurate estimates of the mean ΔID. If base population allele frequencies are unknown and current frequencies are used in the estimations, the largest correlation with F_IBD is generally obtained by F_LH1 and the best estimator of ΔID is F_YA2. The estimators F_VR2 and F_LH2 have the poorest performance in most scenarios. The assumption that base population allele frequencies are equal to 0.5 results in very biased estimates of the average inbreeding coefficient but they are highly correlated with F_IBD and give relatively good estimates of ΔID. Estimates obtained directly from marker homozygosity (F_HOM) substantially overestimated ΔID. Estimates based on runs of homozygosity (F_ROH) provide accurate estimates of inbreeding and ΔID. Finally, estimates based on pedigree (F_PED) show a lower correlation with F_IBD than molecular estimators but provide rather accurate estimates of ΔID. An analysis of data from a pig population supports the main findings of the simulations.

CONCLUSIONS

When base population allele frequencies are known, all marker-allele frequency-based estimators of inbreeding coefficients generally show a high correlation with F_IBD and provide good estimates of ΔID. When base population allele frequencies are unknown, F_LH1 is the marker frequency-based estimator that is most correlated with F_IBD, and F_YA2 provides the most accurate estimates of ΔID. Estimates from F_ROH are also very precise in most scenarios. The estimators F_VR2 and F_LH2 have the poorest performances.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Inbreeding depression and the probability of racing in the Thoroughbred horse

Small effective population sizes and active inbreeding can lead to inbreeding depression due to deleterious recessive mutations exposed in the homozygous state. The Thoroughbred racehorse has low levels of population genetic diversity, but the effects of genomic inbreeding in the population are unknown. Here, we quantified inbreeding based on runs of homozygosity (ROH) using 297 K SNP genotypes from 6128 horses born in Europe and Australia, of which 13.2% were unraced. We show that a 10% increase in inbreeding (F_ROH) is associated with a 7% lower probability of ever racing. Moreover, a ROH-based genome-wide association study identified a haplotype on ECA14 which, in its homozygous state, is linked to a 32.1% lower predicted probability of ever racing, independent of F_ROH. The haplotype overlaps a candidate gene, EFNA5, that is highly expressed in cartilage tissue, which when damaged is one of the most common causes of catastrophic musculoskeletal injury in racehorses. Genomics-informed breeding aiming to reduce inbreeding depression and avoid damaging haplotype carrier matings will improve population health and racehorse welfare.

Runs of Homozygosity and Quantitative Trait Locus/Association for Semen Parameters in Selected Chinese and South African Beef Cattle

In this study, runs of homozygosity (ROH) and quantitative trait locus/association (QTL) for semen parameters in selected Chinese and South African beef cattle breed were estimated. The computed results showed 7516 ROH were observed between classes 0−5 Mb with no ROH observed in classes >40 Mb. Distribution of ROH showed high level of genomic coverage for ANG, NGU, CSI, and BEL breeds. Approximately 13 genomic regions with QTL were controlling sperm motility, sperm concentration, semen volume, sperm count, sperm head abnormalities, sperm tail abnormalities, sperm integrity, and percentage of abnormal sperm traits. Nine candidate genes, CDF9, MARCH1, WDR19, SLOICI, ST7, DOP1B, CFAF9, INHBA, and ADAMTS1, were suggested to be associated with above mentioned QTL traits. The results for inbreeding coefficient showed moderate correlation between FROH vs FHOM at 0.603 and high correlation between FROH 0−5 Mb 0.929, and lowest correlation for 0−>40 Mb 0.400. This study suggested recent inbreeding in CSI, BEL, ANG, BON, SIM, and NGU breeds. Furthermore, it highlighted varied inbreeding levels and identified QTL for semen traits and genes of association. These results can assist in implementation of genetic improvement strategies for bulls and provide awareness and proper guidelines in developing breeding programs.

Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data

Background

Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F_ROH) and homozygosity (F_HOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein.

Results

The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F_ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F_ROH and F_HOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits.

Conclusion

Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08583-9.

Candidate Genes in Bull Semen Production Traits: An Information Approach Review

Semen quality plays a crucial role in the successful implementation of breeding programs, especially where artificial insemination (AI) is practiced. Bulls with good semen traits have good fertility and can produce a volume of high semen per ejaculation. The aim of this review is to use an information approach to highlight candidate genes and their relation to bull semen production traits. The use of genome-wide association studies (GWAS) has been demonstrated to be successful in identifying genomic regions and individual variations associated with production traits. Studies have reported over 40 genes associated with semen traits using Illumina BeadChip single-nucleotide polymorphism (SNPs).

A 1-bp deletion in bovine QRICH2 causes low sperm count and immotile sperm with multiple morphological abnormalities

Background

Semen quality and insemination success are monitored in artificial insemination bulls to ensure high male fertility rates. Only ejaculates that fulfill minimum quality requirements are processed and eventually used for artificial inseminations. We examined 70,990 ejaculates from 1343 Brown Swiss bulls to identify bulls from which all ejaculates were rejected due to low semen quality. This procedure identified a bull that produced 12 ejaculates with an aberrantly small number of sperm (0.2 ± 0.2 × 109 sperm per mL) which were mostly immotile due to multiple morphological abnormalities.

Results

The genome of this bull was sequenced at a 12× coverage to investigate a possible genetic cause. Comparing the sequence variant genotypes of this bull with those from 397 fertile bulls revealed a 1-bp deletion in the coding sequence of the QRICH2 gene which encodes the glutamine rich 2 protein, as a compelling candidate causal variant. This 1-bp deletion causes a frameshift in translation and a premature termination codon (ENSBTAP00000018337.1:p.Cys1644AlafsTer52). The analysis of testis transcriptomes from 76 bulls showed that the transcript with the premature termination codon is subject to nonsense-mediated mRNA decay. The 1-bp deletion resides in a 675-kb haplotype that includes 181 single nucleotide polymorphisms (SNPs) from the Illumina BovineHD Bead chip. This haplotype segregates at a frequency of 5% in the Brown Swiss cattle population. Our analysis also identified another bull that carried the 1-bp deletion in the homozygous state. Semen analyses from the second bull confirmed low sperm concentration and immotile sperm with multiple morphological abnormalities that primarily affect the sperm flagellum and, to a lesser extent, the sperm head.

Conclusions

A recessive loss-of-function allele of the bovine QRICH2 gene likely causes low sperm concentration and immotile sperm with multiple morphological abnormalities. Routine sperm analyses unambiguously identify homozygous bulls for this allele. A direct gene test can be implemented to monitor the frequency of the undesired allele in cattle populations.

A bioinformatic search for correspondence between differentially expressed genes of domestic versus wild animals and orthologous human genes altering reproductive potential

One of the greatest achievements of genetics in the 20th century is D.K. Belyaev’s discovery of destabilizing selection during the domestication of animals and that this selection affects only gene expression regulation (not gene structure) and inf luences systems of neuroendocrine control of ontogenesis in a stressful environment. Among the experimental data generalized by Belyaev’s discovery, there are also f indings about accelerated extinc tion of testes’ hormonal function and disrupted seasonality of reproduction of domesticated foxes in comparison with their wild congeners. To date, Belyaev’s discovery has already been repeatedly conf irmed, for example, by independent observations during deer domestication, during the use of rats as laboratory animals, after the reintroduction of endangered species such as Przewalski’s horse, and during the creation of a Siberian reserve population of the Siberian grouse when it had reached an endangered status in natural habitats. A genome-wide comparison among humans, several domestic animals, and some of their wild congeners has given rise to the concept of self-domestication syndrome, which includes autism spectrum disorders. In our previous study, we created a bioinformatic model of human self-domestication syndrome using differentially expressed genes (DEGs; of domestic animals versus their wild congeners) orthologous to the human genes (mainly, nervous-system genes) whose changes in expression affect reproductive potential, i. e., growth of the number of humans in the absence of restrictions caused by limiting factors. Here, we applied this model to 68 human genes whose changes in expression alter the reproductive health of women and men and to 3080 DEGs of domestic versus wild animals. As a result, in domestic animals, we identif ied 16 and 4 DEGs, the expression changes of which are codirected with changes in the expression of the human orthologous genes decreasing and increasing human reproductive potential, respectively. The wild animals had 9 and 11 such DEGs, respectively. This difference between domestic and wild animals was signif icant according to Pearson’s χ2 test (p < 0.05) and Fisher’s exact test (p < 0.05). We discuss the results from the standpoint of restoration of endangered animal species whose natural habitats are subject to an anthropogenic impact.

Genome-wide detection of non-additive quantitative trait loci for semen production traits in beef and dairy bulls

Semen production traits are important aspects of bull fertility, because semen quantity leads to direct profits for artificial insemination centres, and semen quality is associated with the probability of achieving a pregnancy. Most genome-wide association studies (GWASs) for semen production traits have assumed that each quantitative trait locus (QTL) has an additive effect. However, GWASs that account for non-additive effects are also important in fitness traits, such as bull fertility. Here, we performed a GWAS using models that accounted for additive and non-additive effects to evaluate the importance of non-additive effects on five semen production traits in beef and dairy bulls. A total of 65 463 records for 615 Japanese Black bulls (JB) and 50 734 records for 873 Holstein bulls (HOL), which were previously genotyped using the Illumina BovineSNP50 BeadChip, were used to estimate genetic parameters and perform GWAS. The heritability estimates were low (ranged from 0.11 to 0.23), and the repeatability estimates were low to moderate (ranged from 0.28 to 0.45) in both breeds. The estimated repeatability was approximately twice as high as the estimated heritability for all traits. In this study, only one significant region with an additive effect was detected in each breed, but multiple significant regions with non-additive effects were detected for each breed. In particular, the region at approximately 64 Mbp on Bos taurus autosome 17 had the highest significant non-additive effect on four semen production traits in HOL. The rs41843851 single nucleotide polymorphism (SNP) in the region had a much lower P-value for the non-additive effect (P-value = 1.1 × 10^-31) than for the additive effect (P-value = 1.1 × 10^-8) in sperm motility. The AA and AB genotypes on the SNP had a higher phenotype than the BB genotype in HOL, and there was no bull with the BB genotype in JB. Our results showed that non-additive QTLs affect semen production traits, and a novel QTL accounting for non-additive effects could be detected by GWAS. This study provides new insights into non-additive QTLs that affect fitness traits, such as semen production traits in beef and dairy bulls.

Genetic Diversity, Admixture and Analysis of Homozygous-by-Descent (HBD) Segments of Russian Wild Boar

The wild boar is the wild ancestor of the domestic pig and one of the most common species of ungulates. At the beginning of the 20th century, the wild boar was practically exterminated in the European part of Russia. In the period 1935-1988, 7705 boars were caught in various regions of the European part of Russia, the Far East, Ukraine, Belarus, Kyrgyzstan, Kazakhstan, Latvia, Lithuania, Estonia, Tajikistan and resettled in the territory of Russia. Asian and European wild boars dwell the territory of Russia. The aim of our research was to study the genetic diversity and structure of wild boar populations in different regions of Russia using genome-wide genotyping. We have determined the genetic distances, population structure, parameters of genetic diversity and significantly expanded our understanding of the genetic state of the Russian wild boar. For the first time, we calculated autozygosity of the wild boar of the European and Asian subspecies using Homozygous-by-Descent (HBD) Segments analysis, which is important in terms of population recovery. We also found evidence of hybridization between Russian wild boar and domestic pigs. A group of European wild boars showed introgression of the Asian boar into population. The mean level of the inbreeding coefficient in European wild boar was higher than in Asian wild boar, and combined groups of the European boar had higher inbreeding coefficient than Russian wild boars. These results obtained can be used in population management.

Evolution of inbreeding: a gaze into five Italian beef cattle breeds history

In the last decades, intensive selection programs have led to sustained increases of inbreeding in dairy cattle, a feature that might have adverse consequences on the viability and phenotypic performance of their offspring. This study aimed to determine the evolution of inbreeding of five Italian beef cattle breeds (Marchigiana, Chianina, Romagnola, Maremmana, and Podolica) during a period of almost 20 years (2002–2019). The estimates of Ho, He, F_hat2, and F_ped averaged across years (2002–2019) in the studied breeds fluctuated between 0.340–0.401, 0.348–0.392, –0.121–0.072, and 0.000–0.068, respectively. Moreover, annual rates of increase of the estimated inbreeding coefficients have been very low (F_hat2 = 0.01–0.02%; F_ped = 0.003–0.004%). The use of a high number of bulls combined with strategies implemented by the Association of Italian Beef Cattle Breeders ANABIC to minimize inbreeding might explain these results. Despite the fact that diversity and inbreeding have remained quite stable during the last two decades, we have detected a sustained decrease of the population effective size of these five breeds. Such results should be interpreted with caution due to the inherent difficulty of estimating N_e from SNPs data in a reliable manner.

Runs of homozygosity analysis reveals consensus homozygous regions affecting production traits in Chinese Simmental beef cattle

BACKGROUND

Genomic regions with a high frequency of runs of homozygosity (ROH) are related to important traits in farm animals. We carried out a comprehensive analysis of ROH and evaluated their association with production traits using the BovineHD (770 K) SNP array in Chinese Simmental beef cattle.

RESULTS

We detected a total of 116,953 homozygous segments with 2.47Gb across the genome in the studied population. The average number of ROH per individual was 99.03 and the average length was 117.29 Mb. Notably, we detected 42 regions with a frequency of more than 0.2. We obtained 17 candidate genes related to body size, meat quality, and reproductive traits. Furthermore, using Fisher's exact test, we found 101 regions were associated with production traits by comparing high groups with low groups in terms of production traits. Of those, we identified several significant regions for production traits (P < 0.05) by association analysis, within which candidate genes including ECT2, GABRA4, and GABRB1 have been previously reported for those traits in beef cattle.

CONCLUSIONS

Our study explored ROH patterns and their potential associations with production traits in beef cattle. These results may help to better understand the association between production traits and genome homozygosity and offer valuable insights into managing inbreeding by designing reasonable breeding programs in farm 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.

Genomic partitioning of inbreeding depression in humans

Across species, offspring of related individuals often exhibit significant reduction in fitness-related traits, known as inbreeding depression (ID), yet the genetic and molecular basis for ID remains elusive. Here, we develop a method to quantify enrichment of ID within specific genomic annotations and apply it to human data. We analyzed the phenomes and genomes of ∼350,000 unrelated participants of the UK Biobank and found, on average of over 11 traits, significant enrichment of ID within genomic regions with high recombination rates (>21-fold; p < 10^-5), with conserved function across species (>19-fold; p < 10^-4), and within regulatory elements such as DNase I hypersensitive sites (∼5-fold; p = 8.9 × 10^-7). We also quantified enrichment of ID within trait-associated regions and found suggestive evidence that genomic regions contributing to additive genetic variance in the population are enriched for ID signal. We find strong correlations between functional enrichment of SNP-based heritability and that of ID (r = 0.8, standard error: 0.1). These findings provide empirical evidence that ID is most likely due to many partially recessive deleterious alleles in low linkage disequilibrium regions of the genome. Our study suggests that functional characterization of ID may further elucidate the genetic architectures and biological mechanisms underlying complex traits and diseases.

Grid search approach to discriminate between old and recent inbreeding using phenotypic, pedigree and genomic information

Background

Although inbreeding caused by the mating of animals related through a recent common ancestor is expected to have more harmful effects on phenotypes than ancient inbreeding (old inbreeding), estimating these effects requires a clear definition of recent (new) and ancient (old) inbreeding. Several methods have been proposed to classify inbreeding using pedigree and genomic data. Unfortunately, these methods are largely based on heuristic criteria such as the number of generations from a common ancestor or length of runs of homozygosity (ROH) segments. To mitigate these deficiencies, this study aimed to develop a method to classify pedigree and genomic inbreeding into recent and ancient classes based on a grid search algorithm driven by the assumption that new inbreeding tends to have a more pronounced detrimental effect on traits. The proposed method was tested using a cattle population characterized by a deep pedigree.

Results

Effects of recent and ancient inbreeding were assessed on four growth traits (birth, weaning and yearling weights and average daily gain). Thresholds to classify inbreeding into recent and ancient classes were trait-specific and varied across traits and sources of information. Using pedigree information, inbreeding generated in the last 10 to 11 generations was considered as recent. When genomic information (ROH) was used, thresholds ranged between four to seven generations, indicating, in part, the ability of ROH segments to characterize the harmful effects of inbreeding in shorter periods of time. Nevertheless, using the proposed classification method, the discrimination between new and old inbreeding was less robust when ROH segments were used compared to pedigree. Using several model comparison criteria, the proposed approach was generally better than existing methods. Recent inbreeding appeared to be more harmful across the growth traits analyzed. However, both new and old inbreeding were found to be associated with decreased yearling weight and average daily gain.

Conclusions

The proposed method provided a more objective quantitative approach for the classification of inbreeding. The proposed method detected a clear divergence in the effects of old and recent inbreeding using pedigree data and it was superior to existing methods for all analyzed traits. Using ROH data, the discrimination between old and recent inbreeding was less clear and the proposed method was superior to existing approaches for two out of the four analyzed traits. Deleterious effects of recent inbreeding were detected sooner (fewer generations) using genomic information than pedigree. Difference in the results using genomic and pedigree information could be due to the dissimilarity in the number of generations to a common ancestor. Additionally, the uncertainty associated with the identification of ROH segments and associated inbreeding could have an effect on the results. Potential biases in the estimation of inbreeding effects may occur when new and old inbreeding are discriminated based on arbitrary thresholds. To minimize the impact of inbreeding, mating designs should take the different inbreeding origins into consideration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07872-z.

How Depressing Is Inbreeding? A Meta-Analysis of 30 Years of Research on the Effects of Inbreeding in Livestock

Inbreeding depression has been widely documented for livestock and other animal and plant populations. Inbreeding is generally expected to have a stronger unfavorable effect on fitness traits than on other traits. Traditionally, the degree of inbreeding depression in livestock has been estimated as the slope of the linear regression of phenotypic values on pedigree-based inbreeding coefficients. With the increasing availability of SNP-data, pedigree inbreeding can now be replaced by SNP-based measures. We performed a meta-analysis of 154 studies, published from 1990 to 2020 on seven livestock species, and compared the degree of inbreeding depression (1) across different trait groups, and (2) across different pedigree-based and SNP-based measures of inbreeding. Across all studies and traits, a 1% increase in pedigree inbreeding was associated with a median decrease in phenotypic value of 0.13% of a trait’s mean, or 0.59% of a trait’s standard deviation. Inbreeding had an unfavorable effect on all sorts of traits and there was no evidence for a stronger effect on primary fitness traits (e.g., reproduction/survival traits) than on other traits (e.g., production traits or morphological traits). p-values of inbreeding depression estimates were smaller for SNP-based inbreeding measures than for pedigree inbreeding, suggesting more power for SNP-based measures. There were no consistent differences in p-values for percentage of homozygous SNPs, inbreeding based on runs of homozygosity (ROH) or inbreeding based on a genomic relationship matrix. The number of studies that directly compares these different measures, however, is limited and comparisons are furthermore complicated by differences in scale and arbitrary definitions of particularly ROH-based inbreeding. To facilitate comparisons across studies in future, we provide the dataset with inbreeding depression estimates of 154 studies and stress the importance of always reporting detailed information (on traits, inbreeding coefficients, and models used) along with inbreeding depression estimates.

Genome-wide association study and inbreeding depression on body size traits in Qira black sheep (Ovis aries)

Qira black sheep is a famous indigenous sheep breed in China. The objectives of this study are to identify candidate genes related to body size, and to estimate the level of inbreeding depression on body size based on runs of homozygosity in Qira black sheep. Here, 188 adult Qira black sheep were genotyped with a high density (630 K) SNP chip and genome-wide association study for body weight and body size traits (including withers height, body slanting length, tail length, chest girth, chest width, and chest depth) were performed using an additive linear model. In consequence, 12 genome- and chromosome-wide significant SNPs and, accordingly, six candidate genes involved in muscle differentiation, metabolism and cell processes were identified. Of them, ZNF704 (zinc finger protein 704) was identified for body weight; AK2 (adenylate kinase 2) and PARK2 (parkin RBR E3 ubiquitin protein ligase) for tail length; MOCOS (molybdenum cofactor sulfurase) and ELP2 (elongator acetyltransferase complex subunit 2) for chest width; and MFAP1 (microfibril associated protein 1) for chest girth. Additionally, inbreeding depressions on body size were observed in the current herd. These results will provide insightful understandings into the genetic mechanisms of adult body size, and into the conservation and utilization of Qira black sheep.

Mutation load decreases with haplotype age in wild Soay sheep

Runs of homozygosity (ROH) are pervasive in diploid genomes and expose the effects of deleterious recessive mutations, but how exactly these regions contribute to variation in fitness remains unclear. Here, we combined empirical analyses and simulations to explore the deleterious effects of ROH with varying genetic map lengths in wild Soay sheep. Using a long-term dataset of 4879 individuals genotyped at 417K SNPs, we found that inbreeding depression increases with ROH length. A 1% genomic increase in long ROH (>12.5 cM) reduced the odds of first-year survival by 12.4% compared to only 7.7% for medium ROH (1.56-12.5 cM), whereas short ROH (<1.56 cM) had no effect on survival. We show by forward genetic simulations that this is predicted: compared to shorter ROH, long ROH will have higher densities of deleterious alleles, with larger average effects on fitness and lower population frequencies. Taken together, our results are consistent with the idea that the mutation load decreases in older haplotypes underlying shorter ROH, where purifying selection has had more time to purge deleterious mutations. Finally, our study demonstrates that strong inbreeding depression can persist despite ongoing purging in a historically small population.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

Autosomal recessive loci contribute significantly to quantitative variation of male fertility in a dairy cattle population

BACKGROUND

Cattle are ideally suited to investigate the genetics of male fertility. Semen from individual bulls is used for thousands of artificial inseminations for which the fertilization success is monitored. Results from the breeding soundness examination and repeated observations of semen quality complement the fertility evaluation for each bull.

RESULTS

In a cohort of 3881 Brown Swiss bulls that had genotypes at 683,609 SNPs, we reveal four novel recessive QTL for male fertility on BTA1, 18, 25, and 26 using haplotype-based association testing. A QTL for bull fertility on BTA1 is also associated with sperm head shape anomalies. All other QTL are not associated with any of the semen quality traits investigated. We perform complementary fine-mapping approaches using publicly available transcriptomes as well as whole-genome sequencing data of 125 Brown Swiss bulls to reveal candidate causal variants. We show that missense or nonsense variants in SPATA16, VWA3A, ENSBTAG00000006717 and ENSBTAG00000019919 are in linkage disequilibrium with the QTL. Using whole-genome sequence data, we detect strong association (P = 4.83 × 10^- 12) of a missense variant (p.Ile193Met) in SPATA16 with male fertility. However, non-coding variants exhibit stronger association at all QTL suggesting that variants in regulatory regions contribute to variation in bull fertility.

CONCLUSION

Our findings in a dairy cattle population provide evidence that recessive variants may contribute substantially to quantitative variation in male fertility in mammals. Detecting causal variants that underpin variation in male fertility remains difficult because the most strongly associated variants reside in poorly annotated non-coding regions.

Runs of homozygosity and analysis of inbreeding depression

Pedigree information was traditionally used to assess inbreeding. The availability of high-density marker panels provides an alternative to assess inbreeding, particularly in the presence of incomplete and error-prone pedigrees. Assessment of autozygosity across chromosomal segments using runs of homozygosity (ROH) has emerged as a valuable tool to estimate inbreeding due to its general flexibility and ability to quantify the chromosomal contribution to genome-wide inbreeding. Unfortunately, the identification of ROH segments is sensitive to the parameters used during the search process. These parameters are heuristically set, leading to significant variation in the results. The minimum length required to identify an ROH segment has major effects on the estimation of inbreeding and inbreeding depression, yet it is arbitrarily set. To overcome this limitation, a search algorithm to approximate mutation enrichment was developed to determine the minimum length of ROH segments. It consists of finding genome segments with significant effect differences in trait means between animals with high and low burdens of autozygous intervals with a specific length. The minimum length could be determined heuristically as the smallest interval at which a significant signal is detected. The proposed method was tested in an inbred Hereford cattle population genotyped for 30,220 SNPs. Phenotypes recorded for six traits were used for the approximation of mutation loads. The estimated minimum length was around 1 Mb for yearling weight (YW) and average daily gain (ADG) and 4 Mb for birth weight and weaning weight. These trait-specific thresholds estimated using the proposed method could be attributed to a trait-dependent effect of homozygosity. The detection of significant inbreeding effects was well aligned with the estimated thresholds, especially for YW and ADG. Although highly deleterious alleles are expected to be more frequent in recent inbreeding (long ROH), short ROH segments (<5 Mb) could contain a large number of less deleterious mutations with substantial joint effects on some traits (YW and ADG). Our results highlight the importance of accurate estimation of the ROH-based inbreeding and the necessity to consider a trait-specific minimum length threshold for the identification of ROH segments in inbreeding depression analyses. These thresholds could be determined using the proposed method provided the availability of phenotypic information.

On the estimation of inbreeding depression using different measures of inbreeding from molecular markers

The inbreeding coefficient (F) of individuals can be estimated from molecular marker data, such as SNPs, using measures of homozygosity of individual markers or runs of homozygosity (ROH) across the genome. These different measures of F can then be used to estimate the rate of inbreeding depression (ID) for quantitative traits. Some recent simulation studies have investigated the accuracy of this estimation with contradictory results. Whereas some studies suggest that estimates of inbreeding from ROH account more accurately for ID, others suggest that inbreeding measures from SNP-by-SNP homozygosity giving a large weight to rare alleles are more accurate. Here, we try to give more light on this issue by carrying out a set of computer simulations considering a range of population genetic parameters and population sizes. Our results show that the previous studies are indeed not contradictory. In populations with low effective size, where relationships are more tight and selection is relatively less intense, F measures based on ROH provide very accurate estimates of ID whereas SNP-by-SNP-based F measures with high weight to rare alleles can show substantial upwardly biased estimates of ID. However, in populations of large effective size, with more intense selection and trait allele frequencies expected to be low if they are deleterious for fitness because of purifying selection, average estimates of ID from SNP-by-SNP-based F values become unbiased or slightly downwardly biased and those from ROH-based F values become slightly downwardly biased. The noise attached to all these estimates, nevertheless, can be very high in large-sized populations. We also investigate the relationship between the different F measures and the homozygous mutation load, which has been suggested as a proxy of inbreeding depression.

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

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

An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree

The estimation of the inbreeding coefficient (F) is essential for the study of inbreeding depression (ID) or for the management of populations under conservation. Several methods have been proposed to estimate the realized F using genetic markers, but it remains unclear which one should be used. Here we used whole-genome sequence data for 245 individuals from a Holstein cattle pedigree to empirically evaluate which estimators best capture homozygosity at variants causing ID, such as rare deleterious alleles or loci presenting heterozygote advantage and segregating at intermediate frequency. Estimators relying on the correlation between uniting gametes (F_UNI) or on the genomic relationships (F_GRM) presented the highest correlations with these variants. However, homozygosity at rare alleles remained poorly captured. A second group of estimators relying on excess homozygosity (F_HOM), homozygous-by-descent segments (F_HBD), runs-of-homozygosity (F_ROH) or on the known genealogy (F_PED) was better at capturing whole-genome homozygosity, reflecting the consequences of inbreeding on all variants, and for young alleles with low to moderate frequencies (0.10 < . < 0.25). The results indicate that F_UNI and F_GRM might present a stronger association with ID. However, the situation might be different when recessive deleterious alleles reach higher frequencies, such as in populations with a small effective population size. For locus-specific inbreeding measures or at low marker density, the ranking of the methods can also change as F_HBD makes better use of the information from neighboring markers. Finally, we confirmed that genomic measures are in general superior to pedigree-based estimates. In particular, F_PED was uncorrelated with locus-specific homozygosity.

Inbreeding depression across the genome of Dutch Holstein Friesian dairy cattle

BACKGROUND

Inbreeding depression refers to the decrease in mean performance due to inbreeding. Inbreeding depression is caused by an increase in homozygosity and reduced expression of (on average) favourable dominance effects. Dominance effects and allele frequencies differ across loci, and consequently inbreeding depression is expected to differ along the genome. In this study, we investigated differences in inbreeding depression across the genome of Dutch Holstein Friesian cattle, by estimating dominance effects and effects of regions of homozygosity (ROH).

METHODS

Genotype (75 k) and phenotype data of 38,792 cows were used. For nine yield, fertility and udder health traits, GREML models were run to estimate genome-wide inbreeding depression and estimate additive, dominance and ROH variance components. For this purpose, we introduced a ROH-based relationship matrix. Additive, dominance and ROH effects per SNP were obtained through back-solving. In addition, a single SNP GWAS was performed to identify significant additive, dominance or ROH associations.

RESULTS

Genome-wide inbreeding depression was observed for all yield, fertility and udder health traits. For example, a 1% increase in genome-wide homozygosity was associated with a decrease in 305-d milk yield of approximately 99 kg. For yield traits only, including dominance and ROH effects in the GREML model resulted in a better fit (P < 0.05) than a model with only additive effects. After correcting for the effect of genome-wide homozygosity, dominance and ROH variance explained less than 1% of the phenotypic variance for all traits. Furthermore, dominance and ROH effects were distributed evenly along the genome. The most notable region with a favourable dominance effect for yield traits was on chromosome 5, but overall few regions with large favourable dominance effects and significant dominance associations were detected. No significant ROH-associations were found.

CONCLUSIONS

Inbreeding depression was distributed quite equally along the genome and was well captured by genome-wide homozygosity. These findings suggest that, based on 75 k SNP data, there is little benefit of accounting for region-specific inbreeding depression in selection schemes.

Effect of recent and ancient inbreeding on production and fertility traits in Canadian Holsteins

BACKGROUND

Phenotypic performances of livestock animals decline with increasing levels of inbreeding, however, the noticeable decline known as inbreeding depression, may not be due only to the total level of inbreeding, but rather could be distinctly associated with more recent or more ancient inbreeding. Therefore, splitting inbreeding into different age classes could help in assessing detrimental effects of different ages of inbreeding. Hence, this study sought to investigate the effect of recent and ancient inbreeding on production and fertility traits in Canadian Holstein cattle with both pedigree and genomic records. Furthermore, inbreeding coefficients were estimated using traditional pedigree measure (F_PED) and genomic measures using segment based (F_ROH) and marker-by-marker (F_GRM) based approaches.

RESULTS

Inbreeding depression was found for all production and most fertility traits, for example, every 1% increase in F_PED, F_ROH and F_GRM was observed to cause a - 44.71, - 40.48 and - 48.72 kg reduction in 305-day milk yield (MY), respectively. Similarly, an extension in first service to conception (FSTC) of 0.29, 0.24 and 0.31 day in heifers was found for every 1% increase in F_PED, F_ROH and F_GRM, respectively. Fertility traits that did not show significant depression were observed to move in an unfavorable direction over time. Splitting both pedigree and genomic inbreeding into age classes resulted in recent age classes showing more detrimental inbreeding effects, while more distant age classes caused more favorable effects. For example, a - 1.56 kg loss in 305-day protein yield (PY) was observed for every 1% increase in the most recent pedigree age class, whereas a 1.33 kg gain was found per 1% increase in the most distant pedigree age class.

CONCLUSIONS

Inbreeding depression was observed for production and fertility traits. In general, recent inbreeding had unfavorable effects, while ancestral inbreeding had favorable effects. Given that more negative effects were estimated from recent inbreeding when compared to ancient inbreeding suggests that recent inbreeding should be the primary focus of selection programs. Also, further work to identify specific recent homozygous regions negatively associated with phenotypic traits could be investigated.

Genome-wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH-based coefficient of inbreeding (F_ROH ) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome-wide ROH-traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.

Genome-Wide Assessment of Runs of Homozygosity in Chinese Wagyu Beef Cattle

Runs of homozygosity (ROH) are continuous homozygous regions that generally exist in the DNA sequence of diploid organisms. Identifications of ROH leading to reduction in performance can provide valuable insight into the genetic architecture of complex traits. Here, we evaluated genome-wide patterns of homozygosity and their association with important traits in Chinese Wagyu beef cattle. We identified a total of 29,271 ROH segments from 462 animals. Within each animal, an average number of ROH was 63.36 while an average length was 62.19 Mb. To evaluate the enrichment of ROH across genomes, we initially identified 280 ROH regions by merging ROH events across all individuals. Of these, nine regions containing 154 candidate genes, were significantly associated with six traits (body height, chest circumference, fat coverage, backfat thickness, ribeye area, and carcass length; p < 0.01). Moreover, we found 26 consensus ROH regions with frequencies exceeding 10%, and several regions overlapped with QTLs, which are associated with body weight, calving ease, and stillbirth. Among them, we observed 41 candidate genes, including BCKDHB, MAB21L1, SLC2A13, FGFR3, FGFRL1, CPLX1, CTNNA1, CORT, CTNNBIP1, and NMNAT1, which have been previously reported to be related to body conformation, meat quality, susceptibility, and reproductive traits. In summary, we assessed genome-wide autozygosity patterns and inbreeding levels in Chinese Wagyu beef cattle. Our study identified many candidate regions and genes overlapped with ROH for several important traits, which could be unitized to assist the design of a selection mating strategy in beef cattle.

Ancient Homozygosity Segments in West African Djallonké Sheep Inform on the Genomic Impact of Livestock Adaptation to the Environment

A sample of Burkina Faso Djallonké (West African Dwarf) sheep was analyzed to identify stretches of homozygous segments (runs of homozygosity; ROH) overlapping with ancient homozygosity-by-descent (HBD) segments. HBD segments were considered ancient if they were likely to be inherited from ancestors living from 1024 to 2048 generations ago, roughly coinciding with the time in which sheep entered into West Africa. It is hypothesized that such homozygous segments can inform on the effect of the sheep genome of human-mediated selection for adaptation to this harsh environment. PLINK analyses allowed to identify a total of 510 ROH segments in 127 different individuals that could be summarized into 124 different ROH. A total of 32,968 HBD segments were identified on 119 individuals using the software ZooRoH. HBD segments inherited from ancestors living 1024 and 2048 generations ago were identified on 61 individuals. The overlap between consensus ROH identified using PLINK and HBD fragments putatively assigned to generations 1024 and 2048 gave 108 genomic areas located on 17 different ovine chromosomes which were considered candidate regions for gene-annotation enrichment analyses. Functional annotation allowed to identify six statistically significant functional clusters involving 50 candidate genes. Cluster 1 was involved in homeostasis and coagulation; functional clusters 2, 3, and 6 were associated to innate immunity, defense against infections, and white blood cells proliferation and migration, respectively; cluster 4 was involved in parasite resistance; and functional cluster 5, formed by 20 genes, was involved in response to stress. The current analysis confirms the importance of genomic areas associated to immunity, disease resistance, and response to stress for adaptation of sheep to the challenging environment of humid Sub-Saharan West Africa.

Activation of cryptic splicing in bovine WDR19 is associated with reduced semen quality and male fertility

Cattle are ideally suited to investigate the genetics of male reproduction, because semen quality and fertility are recorded for all ejaculates of artificial insemination bulls. We analysed 26,090 ejaculates of 794 Brown Swiss bulls to assess ejaculate volume, sperm concentration, sperm motility, sperm head and tail anomalies and insemination success. The heritability of the six semen traits was between 0 and 0.26. Genome-wide association testing on 607,511 SNPs revealed a QTL on bovine chromosome 6 that was associated with sperm motility (P = 2.5 x 10⁻²⁷), head (P = 2.0 x 10⁻⁴⁴) and tail anomalies (P = 7.2 x 10⁻⁴⁹) and insemination success (P = 9.9 x 10⁻¹³). The QTL harbors a recessive allele that compromises semen quality and male fertility. We replicated the effect of the QTL on fertility (P = 7.1 x 10⁻³²) in an independent cohort of 2481 Brown Swiss bulls. The analysis of whole-genome sequencing data revealed that a synonymous variant (BTA6:58373887C>T, rs474302732) in WDR19 encoding WD repeat-containing protein 19 was in linkage disequilibrium with the fertility-associated haplotype. WD repeat-containing protein 19 is a constituent of the intraflagellar transport complex that is essential for the physiological function of motile cilia and flagella. Bioinformatic and transcription analyses revealed that the BTA6:58373887 T-allele activates a cryptic exonic splice site that eliminates three evolutionarily conserved amino acids from WDR19. Western blot analysis demonstrated that the BTA6:58373887 T-allele decreases protein expression. We make the remarkable observation that, in spite of negative effects on semen quality and bull fertility, the BTA6:58373887 T-allele has a frequency of 24% in the Brown Swiss population. Our findings are the first to uncover a variant that is associated with quantitative variation in semen quality and male fertility in cattle.

In cattle farming, artificial insemination is the most common method of breeding. To ensure high fertilization rates, ejaculate quality and insemination success are closely monitored in artificial insemination bulls. We analyse semen quality, insemination success and microarray-called genotypes at more than 600,000 genome-wide SNP markers of 794 bulls to identify a recessive allele that compromises semen quality. We take advantage of whole-genome sequencing to pinpoint a variant in the coding sequence of WDR19 encoding WD repeat-containing protein 19 that activates a novel exonic splice site. Our results indicate that cryptic splicing in WDR19 is associated with reduced male reproductive performance. This is the first report of a variant that contributes to quantitative variation in bovine semen quality.

Whole-genome homozygosity mapping reveals candidate regions affecting bull fertility in US Holstein cattle

Background

Achieving rapid genetic progress while maintaining adequate genetic diversity is one of the main challenges facing the dairy industry. The increase in inbreeding can be used to monitor the loss of genetic diversity. Inbreeding tends to increase the proportion of homozygous loci, some of which cause homozygosity of recessive alleles that results in reduced performance. This phenomenon is known as inbreeding depression and tends to be most prominent on fitness-related traits, such as male fertility. Traditionally, inbreeding has been monitored using pedigree information, or more recently, genomic data. Alternatively, it can be quantified using runs of homozygosity (ROH), i.e., contiguous lengths of homozygous genotypes observed in an individual’s chromosome.

Results

The objective of this study was to evaluate the association between ROH and sire conception rate. ROH were evaluated using 268 k genetic markers in 11,790 US Holstein bulls. Interestingly, either the sum, mean, or maximum length of ROH were negatively associated with bull fertility. The association analysis between ROH and sire fertility was performed comparing 300 high-fertility vs. 300 low-fertility bulls. Both the average and sum of ROH length were higher in the low-fertility group. The enrichment of ROH regions in bulls with low fertility was assessed using a Fisher’s exact test. Nine regions were significantly enriched in low-fertility compared to high-fertility bulls. Notably, these regions harbor genes that are closely related to sperm biology and male fertility, including genes exclusively or highly expressed in testis.

Conclusions

The results of this study can help not only to manage inbreeding in genomic selection programs by designing custom mating schemes, but also to better understand the mechanisms underlying male fertility in dairy cattle.

Genomic Signatures After Five Generations of Intensive Selective Breeding: Runs of Homozygosity and Genetic Diversity in Representative Domestic and Wild Populations of Turbot (Scophthalmus maximus)

Massive genotyping of single nucleotide polymorphisms (SNP) has opened opportunities for analyzing the way in which selection shapes genomes. Artificial or natural selection usually leaves genomic signatures associated with selective sweeps around the responsible locus. Strong selective sweeps are most often identified either by lower genetic diversity than the genomic average and/or islands of runs of homozygosity (ROHi). Here, we conducted an analysis of selective sweeps in turbot (Scophthalmus maximus) using two SNP datasets from a Northeastern Atlantic population (36 individuals) and a domestic broodstock (46 individuals). Twenty-six families (∼ 40 offspring per family) from this broodstock and three SNP datasets applying differing filtering criteria were used to adjust ROH calling parameters. The best-fitted genomic inbreeding estimate (F_ROH) was obtained by the sum of ROH longer than 1 Mb, called using a 21,615 SNP panel, a sliding window of 37 SNPs and one heterozygous SNP per window allowed. These parameters were used to obtain the ROHi distribution in the domestic and wild populations (49 and 0 ROHi, respectively). Regions with higher and lower genetic diversity within each population were obtained using sliding windows of 37 SNPs. Furthermore, those regions were mapped in the turbot genome against previously reported genetic markers associated with QTL (Quantitative Trait Loci) and outlier loci for domestic or natural selection to identify putative selective sweeps. Out of the 319 and 278 windows surpassing the suggestive pooled heterozygosity thresholds (ZHp) in the wild and domestic population, respectively, 78 and 54 were retained under more restrictive ZHp criteria. A total of 116 suggestive windows (representing 19 genomic regions) were linked to either QTL for production traits, or outliers for divergent or balancing selection. Twenty-four of them (representing 3 genomic regions) were retained under stricter ZHp thresholds. Eleven QTL/outlier markers were exclusively found in suggestive regions of the domestic broodstock, 7 in the wild population and one in both populations; one (broodstock) and two (wild) of those were found in significant regions retained under more restrictive ZHp criteria in the broodstock and the wild population, respectively. Genome mining and functional enrichment within regions associated with selective sweeps disclosed relevant genes and pathways related to aquaculture target traits, including growth and immune-related pathways, metabolism and response to hypoxia, which showcases how this genome atlas of genetic diversity can be a valuable resource to look for candidate genes related to natural or artificial selection in turbot populations.

A multivariate analysis with direct additive and inbreeding depression load effects

Background

Inbreeding is caused by mating between related individuals and its most common consequence is inbreeding depression. Several studies have detected heterogeneity in inbreeding depression among founder individuals, and recently a procedure for predicting hidden inbreeding depression loads associated with founders and the Mendelian sampling of non-founders has been developed. The objectives of our study were to expand this model to predict the inbreeding loads for all individuals in the pedigree and to estimate the covariance between the inbreeding loads and the additive genetic effects for the trait of interest. We tested the proposed approach with simulated data and with two datasets of records on weaning weight from the Spanish Pirenaica and Rubia Gallega beef cattle breeds.

Results

The posterior estimates of the variance components with the simulated datasets did not differ significantly from the simulation parameters. In addition, the correlation between the predicted and simulated inbreeding loads were always positive and ranged from 0.27 to 0.82. The beef cattle datasets comprised 35,126 and 75,194 records on weights between 170 and 250 days of age, and pedigrees of 308,836 and 384,434 individual-sire-dam entries for the Pirenaica and Rubia Gallega breeds, respectively. The posterior mean estimates of the variance of inbreeding depression loads were 29,967.8 and 28,222.4 for the Pirenaica and Rubia Gallega breeds, respectively. They were larger than those of the additive variance (695.0 and 439.8 for Pirenaica and Rubia Gallega, respectively), because they should be understood as the variance of the inbreeding depression achieved by a fully inbred (100%) descendant. Therefore, the inbreeding loads have to be rescaled for smaller inbreeding coefficients. In addition, a strong negative correlation (− 0.43 ± 0.10) between additive effects and inbreeding loads was detected in the Pirenaica, but not in the Rubia Gallega breed.

Conclusions

The results of the simulation study confirmed the ability of the proposed procedure to predict inbreeding depression loads for all individuals in the populations. Furthermore, the results obtained from the two real datasets confirmed the variability in the inbreeding depression loads in both breeds and suggested a negative correlation of the inbreeding loads with the additive genetic effects in the Pirenaica breed.

Selection for bull fertility: a review

Fertility is a critically important factor in cattle production because it directly relates to the ability to produce the offspring necessary to offset costs in production systems. Female fertility has received much attention and has been enhanced through assisted reproductive technologies, as well as genetic selection; however, improving bull fertility has been largely ignored. Improvements in bull reproductive performance are necessary to optimize the efficiency of cattle production. Selection and management to improve bull fertility not only have the potential to increase conception rates but also have the capacity to improve other economically relevant production traits. Bull fertility has reportedly been genetically correlated with traits such as average daily gain, heifer pregnancy, and calving interval. Published studies show that bull fertility traits are low to moderately heritable, indicating that improvements in bull fertility can be realized through selection. Although female fertility has continued to progress according to increasing conception rates, the reported correlation between male and female fertility is low, indicating that male fertility cannot be improved by selection for female fertility. Correlations between several bull fertility traits, such as concentration, number of spermatozoa, motility, and number of spermatozoa abnormalities, vary among studies. Using male fertility traits in selection indices would provide producers with more advanced selection tools. The objective of this review was to discuss current beef bull fertility measurements and to discuss the future of genetic evaluation of beef bull fertility and potential genetic improvement strategies.

Inbreeding depression due to recent and ancient inbreeding in Dutch Holstein-Friesian dairy cattle

Background

Inbreeding decreases animal performance (inbreeding depression), but not all inbreeding is expected to be equally harmful. Recent inbreeding is expected to be more harmful than ancient inbreeding, because selection decreases the frequency of deleterious alleles over time. Selection efficiency is increased by inbreeding, a process called purging. Our objective was to investigate effects of recent and ancient inbreeding on yield, fertility and udder health traits in Dutch Holstein–Friesian cows.

Methods

In total, 38,792 first-parity cows were included. Pedigree inbreeding (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED) was computed and 75 k genotype data were used to compute genomic inbreeding, among others based on regions of homozygosity (ROH) in the genome (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ROH}$$\end{document}FROH).

Results

Inbreeding depression was observed, e.g. a 1% increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ROH}$$\end{document}FROH was associated with a 36.3 kg (SE = 2.4) decrease in 305-day milk yield, a 0.48 day (SE = 0.15) increase in calving interval and a 0.86 unit (SE = 0.28) increase in somatic cell score for day 150 through to 400. These effects equalled − 0.45, 0.12 and 0.05% of the trait means, respectively. When \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED was split into generation-based components, inbreeding on recent generations was more harmful than inbreeding on more distant generations for yield traits. When \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{PED}$$\end{document}FPED was split into new and ancestral components, based on whether alleles were identical-by-descent for the first time or not, new inbreeding was more harmful than ancestral inbreeding, especially for yield traits. For example, a 1% increase in new inbreeding was associated with a 2.42 kg (SE = 0.41) decrease in 305-day fat yield, compared to a 0.03 kg (SE = 0.71) increase for ancestral inbreeding. There were no clear differences between effects of long ROH (recent inbreeding) and short ROH (ancient inbreeding).

Conclusions

Inbreeding depression was observed for yield, fertility and udder health traits. For yield traits and based on pedigree, inbreeding on recent generations was more harmful than inbreeding on distant generations and there was evidence of purging. Across all traits, long and short ROH contributed to inbreeding depression. In future work, inbreeding depression and purging should be assessed in more detail at the genomic level, using higher density information and genomic time series.

Genome-wide mapping of the dominance effects based on breed ancestry for semen traits in admixed Swiss Fleckvieh bulls

Heterosis is the beneficial deviation of crossbred progeny from the average of parental lines for a particular trait. Heterosis is due to nonadditive genetic effects with dominance and epistatic components. Recent advances in genotyping technology have encouraged researchers to estimate and scan heterosis components for a range of traits in crossbred populations, applying various definitions of such components. In this study, we defined the intralocus (dominance) component of heterosis using local genetic ancestry and performed genome-wide association analysis for admixed Swiss Fleckvieh bulls and their parental populations, Red Holstein Friesian and Swiss Simmental, for semen traits. A linear mixed model for 41,824 SNP, including SNP additive genetic, breed additive, and breed dominance effects on 1,178 bulls (148 Red Holstein Friesian, 213 Swiss Simmental, and 817 Swiss Fleckvieh) with a total of 43,782 measurements was performed. In total, 19 significant regions for breed dominance were identified for volume (2 regions on Bos taurus autosome 10 and 22) and percentage of live spermatozoa (17 regions on Bos taurus autosome 3, 4, 5, 7, 13, 14, and 17), and genes associated with spermatogenesis, sperm motility, and male fertility traits were located there. No significant region for breed dominance was detected for total number of spermatozoa. The signals for breed dominance were relatively wide, most likely due to limited numbers of recombination events in a small number of generations (10-15 generations) of crossbreeding in the recent Swiss Fleckvieh composite.

A splice donor variant in CCDC189 is associated with asthenospermia in Nordic Red dairy cattle

Background

Cattle populations are highly amenable to the genetic mapping of male reproductive traits because longitudinal data on ejaculate quality and dense microarray-derived genotypes are available for thousands of artificial insemination bulls. Two young Nordic Red bulls delivered sperm with low progressive motility (i.e., asthenospermia) during a semen collection period of more than four months. The bulls were related through a common ancestor on both their paternal and maternal ancestry. Thus, a recessive mode of inheritance of asthenospermia was suspected.

Results

Both bulls were genotyped at 54,001 SNPs using the Illumina BovineSNP50 Bead chip. A scan for autozygosity revealed that they were identical by descent for a 2.98 Mb segment located on bovine chromosome 25. This haplotype was not found in the homozygous state in 8557 fertile bulls although five homozygous haplotype carriers were expected (P = 0.018). Whole genome-sequencing uncovered that both asthenospermic bulls were homozygous for a mutation that disrupts a canonical 5′ splice donor site of CCDC189 encoding the coiled-coil domain containing protein 189. Transcription analysis showed that the derived allele activates a cryptic splice site resulting in a frameshift and premature termination of translation. The mutated CCDC189 protein is truncated by more than 40%, thus lacking the flagellar C1a complex subunit C1a-32 that is supposed to modulate the physiological movement of the sperm flagella. The mutant allele occurs at a frequency of 2.5% in Nordic Red cattle.

Conclusions

Our study in cattle uncovered that CCDC189 is required for physiological movement of sperm flagella thus enabling active progression of spermatozoa and fertilization. A direct gene test may be implemented to monitor the asthenospermia-associated allele and prevent the birth of homozygous bulls that are infertile. Our results have been integrated in the Online Mendelian Inheritance in Animals (OMIA) database (https://omia.org/OMIA002167/9913/).

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5628-y) contains supplementary material, which is available to authorized users.

AUTALASSO: an automatic adaptive LASSO for genome-wide prediction

Background

Genome-wide prediction has become the method of choice in animal and plant breeding. Prediction of breeding values and phenotypes are routinely performed using large genomic data sets with number of markers on the order of several thousands to millions. The number of evaluated individuals is usually smaller which results in problems where model sparsity is of major concern. The LASSO technique has proven to be very well-suited for sparse problems often providing excellent prediction accuracy. Several computationally efficient LASSO algorithms have been developed, but optimization of hyper-parameters can be demanding.

Results

We have developed a novel automatic adaptive LASSO (AUTALASSO) based on the alternating direction method of multipliers (ADMM) optimization algorithm. The two major hyper-parameters of ADMM are the learning rate and the regularization factor. The learning rate is automatically tuned with line search and the regularization factor optimized using Golden section search. Results show that AUTALASSO provides superior prediction accuracy when evaluated on simulated and real bull data compared to the adaptive LASSO, LASSO and ridge regression implemented in the popular glmnet software.

Conclusions

The AUTALASSO provides a very flexible and computationally efficient approach to GWP, especially when it is important to obtain high prediction accuracy and genetic gain. The AUTALASSO also has the capability to perform GWAS of both additive and dominance effects with smaller prediction error than the ordinary LASSO.

Inbreeding depression in line 1 Hereford cattle population using pedigree and genomic information

This study aimed at assessing inbreeding and its effect on growth and fertility traits using the longtime closed line 1 Hereford cattle population. Inbreeding was estimated based on pedigree (FPED) and genomic information. For the latter, three estimates were derived based on the diagonal elements of the genomic relationship matrix using estimated (FGRM) or fixed (FGRM0.5) minor allele frequencies or runs of homozygosity (ROH) (FROH). A pedigree containing 10,186 animals was used to calculate FPED. Genomic inbreeding was evaluated using 785 animals genotyped for 30,810 SNP. Traits analyzed were birth weight (BWT), weaning weight (WWT), yearling weight (YWT), ADG, and age at first calving (AFC). The number of ROH per animal ranged between 6 and 119 segments with an average of 83. The shortest and longest segments were 1.36 and 64.86 Mb long, respectively, reflecting both ancient and recent inbreeding occurring in the last 30 to 40 generations. The average inbreeding was 29.2%, 16.1%, 30.2%, and 22.9% for FPED, FGRM, FGRM0.5, and FROH, respectively. FROH provided the highest correlations with FPED (r = 0.66). Across paternal half-sib families, with minimal variation in FPED, there were substantial variations in their genomic inbreeding. Inbreeding depression analyses showed that a 1% increase in an animal's FPED resulted in a decrease of 1.20 kg, 2.03 kg, and 0.004 kg/d in WWT, YWT, and ADG, respectively. Maternal inbreeding showed significantly negative effects on progeny growth performance. AFC increased by 1.4 and 0.8 d for each 1% increase in FPED of the cow and her dam, respectively. Using genomic inbreeding, similar impact on growth traits was observed although the magnitude of the effect varied between methods. Across all genomic measures, WWT, YWT, and ADG decreased by 0.21 to 0.53 kg, 0.46 to 1.13 kg, and 0.002 to 0.006 kg/d for each 1% increase in genomic inbreeding, respectively. Four chromosomes (9, 12, 17, and 27) were identified to have a significant association between their homozygosity (FROH-CHR) and growth traits. Variability in genomic inbreeding could be useful when deciding between full and half-sib selection candidates. Despite the high level of inbreeding in this study, its negative impact on growth performance was not as severe as expected, which may be attributed to the purging of the deleterious alleles due to natural or artificial selection over time.