Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations

Using genome‐wide SNP data, we calculated genomic inbreeding coefficients (F_ROH > _1 Mb, F_{ROH > 2 Mb}, F_{ROH > 8 Mb} and F_{ROH > 16 Mb}) derived from runs of homozygosity (ROH) of different lengths (>1, >2, >8 and > 16 Mb) as well as from levels of homozygosity (F_HOM). We compared these values of inbreeding coefficients with those calculated from pedigrees (F_PED) of 1422 bulls comprising Brown Swiss (304), Fleckvieh (502), Norwegian Red (499) and Tyrol Grey (117) cattle breeds. For all four breeds, population inbreeding levels estimated by the genomic inbreeding coefficients F_{ROH > 8 Mb} and F_{ROH > 16 Mb} were similar to the levels estimated from pedigrees. The lowest values were obtained for Fleckvieh (F_PED = 0.014, F_{ROH > 8 Mb} = 0.019 and F_{ROH > 16 Mb} = 0.008); the highest, for Brown Swiss (F_PED = 0.048, F_{ROH > 8 Mb} = 0.074 and F_{ROH > 16 Mb} = 0.037). In contrast, inbreeding estimates based on the genomic coefficients F_{ROH > 1 Mb} and F_{ROH > 2 Mb} were considerably higher than pedigree‐derived estimates. Standard deviations of genomic inbreeding coefficients were, on average, 1.3–1.7‐fold higher than those obtained from pedigrees. Pearson correlations between genomic and pedigree inbreeding coefficients ranged from 0.50 to 0.62 in Norwegian Red (lowest correlations) and from 0.64 to 0.72 in Tyrol Grey (highest correlations). We conclude that the proportion of the genome present in ROH provides a good indication of inbreeding levels and that analysis based on ROH length can indicate the relative amounts of autozygosity due to recent and remote ancestors.     

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (F_PED), or runs of homozygosity (ROH; F_ROH) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of F_ROH. The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.     

Genetic diversity and genomic inbreeding in Japanese Black cows in the islands of Okinawa Prefecture evaluated using single-nucleotide polymorphism array

Maintaining genetic diversity and inbreeding control are important in Japanese Black cattle production, especially in remote areas such as the islands of Okinawa Prefecture. Using a single-nucleotide polymorphism (SNP) array, we evaluated the genetic diversity and genomic inbreeding in Japanese Black cows from the islands of Okinawa Prefecture and compared them to those from other locations across Japan. Linkage disequilibrium decay was slower in cows in the islands of Okinawa Prefecture. The estimated effective population size declined over time in both populations. The genomic inbreeding coefficient (F_ROH) was estimated using long stretches of consecutive homozygous SNPs (runs of homozygosity; ROH). F_ROH was higher in the cows on the islands of Okinawa Prefecture than on other locations. In total, 818 ROH fragments, including those containing NCAPG and PLAG1, which are major quantitative trait loci for carcass weight in Japanese Black cattle, were present at significantly higher frequencies in cows in the islands of Okinawa Prefecture. This suggests that the ROH fragments are under strong selection and that cows in the islands of Okinawa Prefecture have low genetic diversity and high genomic inbreeding relative to those at other locations. SNP arrays are useful tools for evaluating genetic diversity and genomic inbreeding in cattle.     

Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP‐derived metrics

Multilocus homozygosity, measured as the proportion of the autosomal genome in homozygous genotypes or in runs of homozygosity, was compared with the respective pedigree inbreeding coefficients in 64 Iberian pigs genotyped using the Porcine SNP60 Beadchip. Pigs were sampled from a set of experimental animals with a large inbreeding variation born in a closed strain with a completely recorded multi‐generation genealogy. Individual inbreeding coefficients calculated from pedigree were strongly correlated with the different SNP‐derived metrics of homozygosity (r = 0.814–0.919). However, unequal correlations between molecular and pedigree inbreeding were observed at chromosomal level being mainly dependent on the number of SNPs and on the correlation between heterozygosities measured across different loci. A panel of 192 SNPs of intermediate frequencies was selected for genotyping 322 piglets to test inbreeding depression on postweaning growth performance (daily gain and weight at 90 days). The negative effects on these traits of homozygosities calculated from the genotypes of 168 quality‐checked SNPs were similar to those of inbreeding coefficients. The results support that few hundreds of SNPs may be useful for measuring inbreeding and inbreeding depression, when the population structure or the mating system causes a large variance of inbreeding.     

Estimation of inbreeding and effective population size of full‐blood wagyu cattle registered with the American Wagyu Cattle Association

The objective of this research was to examine the population structure of full‐blood (100%) Wagyu cattle registered in the United States with the American Wagyu Association, with the aim of estimating and comparing the levels of inbreeding from both pedigree and genotypic data. A total of 4132 full‐blood Wagyu cattle pedigrees were assessed and used to compute the inbreeding coefficients (F_IT and F_ST) and the effective population size (N_e) from pedigree data for the period 1994 to 2011. In addition to pedigree analysis, 47 full‐blood Wagyu cattle representing eight prominent sire lines in the American Wagyu cattle population were genotyped using the Illumina BovineSNP50 BeadChip. Genotypic data were then used to estimate genomic inbreeding coefficients (F_ROH) by calculating runs of homozygosity. The mean inbreeding coefficient based on the pedigree data was estimated at 4.80%. The effective population size averaged 17 between the years 1994 and 2011 with an increase of 42.9 in 2000 and a drop of 1.8 in 2011. Examination of the runs of homozygosity revealed that the 47 Wagyu cattle from the eight prominent sire lines had a mean genomic inbreeding coefficient (F_ROH) estimated at 9.08% compared to a mean inbreeding coefficient based on pedigree data of 4.8%. These data suggest that the mean genotype inbreeding coefficient of full‐blood Wagyu cattle exceeds the inbreeding coefficient identified by pedigree. Inbreeding has increased slowly at a rate of 0.03% per year over the past 17 years. Wagyu breeders should continue to utilize many sires from divergent lines and consider outcrossing to other breeds to enhance genetic diversity and minimize the adverse effects of inbreeding in Wagyu.     

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.     

Assessing effective population size,coancestry and inbreeding effects on litter size using the pedigree and SNP data in closed lines of the Iberian pig breed

The complete pedigree of two closed Iberian pig lines (Gamito and Torbiscal), with 798 and 4077 reproducers, has been used to measure the evolution of coancestry (f) and inbreeding (F) for autosomal and X‐linked genes along 16 and 28 respective equivalent discrete generations. At the last generation, the mean values of each line were f = 0.41 and 0.22, F = 0.35 and 0.18, f_X = 0.46 and 0.22 and F_X = 0.47 and 0.19, respectively. Other calculated parameters were the effective number of founders (final values, 6.8 and 35.2) and non‐founders (1.5 and 2.4), founder genome equivalents (1.2 and 2.3) and effective population size (16.0 and 57.7). Measures of Torbiscal effective size based on rates of coancestry (66.1), inbreeding (65.0) and linkage disequilibrium (71.0) were estimated from whole‐genome SNP genotyping data. Values of new and old inbreeding and their respective rates by generation were computed to detect purging effects of natural selection. The analysis of 6854 Torbiscal litters showed significant negative impacts of new and fast inbreeding on litter size, as expected from the purging hypothesis: ?0.20 born piglets per litter by a 10% of new inbreeding, and ?0.03 and ?0.02 piglets by 1% of total and new inbreeding rates, respectively. The analysis performed on 1274 litters of the Gamito line failed to show purging effects. The only significant results were reductions in ?0.91 and ?0.17 piglets by a 10% of old and X‐linked genes inbreeding, respectively. These results may be useful for some practical issues in conservation programs of farm or captive wild animals.     

Hepatic steroid inactivating enzymes,hepatic portal blood flow and corpus luteum blood perfusion in cattle

Production from the corpus luteum (CL) and/or hepatic steroid inactivation impacts peripheral concentrations of P4, which can alter reproductive performance. Our primary objective was to examine hepatic steroid inactivating enzymes, portal blood flow, and luteal blood perfusion at 10 days post‐insemination in pregnant versus non‐pregnant beef and dairy cows. Twenty early lactation Holstein cows and 20 lactating commercial beef cows were utilized for this study. At day 10 post‐insemination, hepatic portal blood flow and CL blood perfusion were measured via Doppler ultrasonography. Liver biopsies were collected and frozen for later determination of cytochrome P450 1A (CYP1A), 2C (CYP2C), 3A (CYP3A), uridine diphosphate‐glucuronosyltransferase (UGT) and aldo‐keto reductase 1C (AKR1C) activities. Pregnancy was determined at day 30 post‐insemination and treatment groups were retrospectively assigned as pregnant or non‐pregnant. Data were analyzed using the mixed procedure of SAS. Steroid metabolizing enzyme activity was not different (p > .10) between pregnant versus non‐pregnant beef or dairy cows. Hepatic portal blood flow tended (p < .10) to be increased in pregnant versus non‐pregnant dairy cows. Luteal blood perfusion was increased (p < .05) in pregnant versus non‐pregnant dairy cows. Pregnant dairy cows appear to have an increased rate of hepatic clearance of P4 in combination with increased synthesis from the CL. This could account for the lack of difference in peripheral P4 concentrations between pregnant and non‐pregnant dairy cows. This study highlights the relevance of further investigation into steroid secretion and inactivation and their impact on the maintenance of pregnancy in cattle.     

Different prostaglandin F2α secretion in response to oxytocin injection between pregnant and non‐pregnant cows: effect of the day of oxytocin challenge test for determining the difference

The present study was conducted to determine the difference in plasma prostaglandin F₂α metabolite concentrations following oxytocin (OT) challenge between pregnant and non‐pregnant cows. Experiment 1: cows were subjected to the OT challenge test on days 12, 14 or 16 (day of estrus = day 0) with or without prior insemination and plasma 13,14‐dihydro‐15‐keto prostaglandin F₂α (PGFM) concentrations were measured from ?30 to 180 min after OT injection. On day 16, the increment of plasma PGFM concentrations in response to OT injection was significantly smaller in pregnant than that in cyclic cows. On days 12 and 14, there was little OT‐induced PGFM secretion and no difference in PGFM increase between the pregnant and cyclic cows. Experiment 2: cows were inseminated on day 0 and subjected to the OT challenge test on day 16. Cows were classified into non‐pregnant/early embryonic death (NP/EED), late embryonic death (LED) and pregnant (PREG) groups. The increment of PGFM concentrations in response to OT injection was less in both PREG and LED groups than that in the NP/EED group. In conclusion, plasma PGFM secretion induced by OT is suggested as the base of pregnancy diagnosis prior to returning estrus in cows.     

The contribution of dominance and inbreeding depression in estimating variance components for litter size in Pannon White rabbits

In a synthetic closed population of Pannon White rabbits, additive (V_A), dominance (V_D) and permanent environmental (V_Pe) variance components as well as doe (b_Fd) and litter (b_Fl) inbreeding depression were estimated for the number of kits born alive (NBA), number of kits born dead (NBD) and total number of kits born (TNB). The data set consisted of 18,398 kindling records of 3883 does collected from 1992 to 2009. Six models were used to estimate dominance and inbreeding effects. The most complete model estimated V_A and V_D to contribute 5.5 ± 1.1% and 4.8 ± 2.4%, respectively, to total phenotypic variance (V_P) for NBA; the corresponding values for NBD were 1.9 ± 0.6% and 5.3 ± 2.4%, for TNB, 6.2 ± 1.0% and 8.1 ± 3.2% respectively. These results indicate the presence of considerable V_D. Including dominance in the model generally reduced V_A and V_Pe estimates, and had only a very small effect on inbreeding depression estimates. Including inbreeding covariates did not affect estimates of any variance component. A 10% increase in doe inbreeding significantly increased NBD (b_Fd = 0.18 ± 0.07), while a 10% increase in litter inbreeding significantly reduced NBA (b_Fl = ?0.41 ± 0.11) and TNB (b_Fl = ?0.34 ± 0.10). These findings argue for including dominance effects in models of litter size traits in populations that exhibit significant dominance relationships.     

Implementation of genetic evaluation and mating designs for the endangered local pig breed ‘Bunte Bentheimer’

A pedigree including 1538 individuals of the endangered pig breed ‘Bunte Bentheimer’ and 3008 records of the fertility traits ‘number of piglets born alive’ (NBA) and ‘number of piglets weaned’ (NW) were used to i) characterize the population structure, ii) to estimate genetic (co)variance components and estimated breeding values (EBVs) and iii) to use EBVs for the application of the concept of optimal genetic contributions. The average coefficient of inbreeding increased from F = 0.103 to F = 0.121 within the two recent cohorts. Average rate of inbreeding amounted to 1.66%, which resulted in an effective population size of N_e = 30 animals in the recent cohort. Average generation interval was 3.07 years considering the whole pedigree, and in total, only 612 sows and boars generated offspring. Estimated heritabilities for both traits NBA and NW were 0.12, and the estimated genetic correlation between both traits was 0.96. The variance component due to the service sire was higher than in commercial pig breeds, presumably due to the widespread use of natural service boars. The EBVs for NBA from 333 selection candidates (63 boars and 270 sows) were used to determine optimal genetic contributions. Based on selected animals and their optimal genetic contributions, specific mating designs were evaluated to minimize inbreeding in the next generation. Best results were achieved when using a simulated annealing algorithm and allowing artificial insemination.     

Evaluation of three hormonal protocols for anovulatory lactating cows under regulations restricting the use of estrogenic compounds

When European Union regulations restricted the use of estrogenic compounds in food‐producing animals, refined hormonal protocols were no longer applicable for anovulatory cows. However, Ovsynch and its adaptations are routinely and uniformly applied to all cows regardless of ovarian function. To evaluate their efficacy on anovulatory cows, 143, 147 and 144 anovulatory cows received Ovsynch, Presynch and G6G protocols, respectively. In comparison, 150 cyclic cows were bred without using a synchronized protocol. Results showed that cows in the Presynch group had luteolysis responding to the last prostaglandin F_2α (PGF_2α) injection greater than the Ovsynch group. The serous progesterone levels at the first gonadotropin‐releasing hormone of Ovsych and the last PGF_2α injection was greater in the G6G group than the other two hormonal treatment groups. Concentrations of Ca²⁺ and total protein in cervical mucus in all three hormone‐treated groups before artificial insemination (AI) were significantly different from the controls. The G6G group obtained a greater pregnancy rate compared with Ovsynch and Presynch, but significantly less than the controls. For open cows in the Ovsynch group, estrus rate within 24 days after the first AI was significantly less than the controls. In conclusion, the G6G treatment resulted to better reproductive performance in anovulatory cows.     

Gonadotropin‐Releasing Hormone Administration to Dairy Cows without a Corpus Luteum 4 Weeks after Calving Increases Reproductive Performance

This field study investigated whether the administration of a single dose of gonadotropin‐releasing hormone (GnRH) to dairy cows without a corpus luteum (CL) 4 weeks after calving can improve reproductive performance. Holstein dairy cows underwent ultrasonography to assess the presence of ovarian structures at 29.2 ± 5.2 days post‐partum, and cows were divided into two main groups based on the presence (CL group, n = 230) or absence (non‐CL group, n = 460) of a CL. The non‐CL group was further randomly divided into two subgroups based on the administration of GnRH (non‐CL GnRH group, n = 230) or no GnRH (non‐CL control group, n = 230). Subsets of cows from non‐CL control (n = 166) and non‐CL GnRH (n = 175) groups received a second ultrasonography at 44.5 ± 5.4 days post‐partum to assess CL formation. The percentage of cows with CL at the second ultrasonography was greater in the non‐CL GnRH group (70.9%) than in the non‐CL control group (53.0%, p = 0.0006). The hazard of the first post‐partum insemination by 150 days in milk (DIM) was higher in the CL group than in the non‐CL control group (hazard ratio [HR]: 1.36, p = 0.001). The probability of a pregnancy to the first insemination was higher in non‐CL GnRH (odds ratio [OR]: 1.50, p = 0.04) and CL groups (OR: 1.55, p = 0.03) compared to the non‐CL control group. Furthermore, the hazard of pregnancy by 210 DIM was higher in non‐CL GnRH (HR: 1.30, p = 0.01) and CL (HR: 1.51, p = 0.0001) groups than in the non‐CL control group. In conclusion, administration of GnRH to dairy cows without a CL 4 weeks after calving was associated with an increase in ovulation and improved reproductive performance.     

Estimation of ancestral inbreeding effects on stillbirth,calving ease and birthweight in German Holstein dairy cattle

In this study, the effect of different measurements of ancestral inbreeding on birthweight, calving ease and stillbirth were analysed. Three models were used to estimate the effect of ancestral inbreeding, and the estimated regression coefficient of phenotypic data on different measurements of ancestral inbreeding was used to quantify the effect of ancestral inbreeding. The first model included only one measurement of inbreeding, whereas the second model included the classical inbreeding coefficients and one alternative inbreeding coefficient. The third model included the classical inbreeding coefficients, the interaction between classical inbreeding and ancestral inbreeding, and the classical inbreeding coefficients of the dam. Phenotypic data for this study were collected from February 1998 to December 2008 on three large commercial milk farms. During this time, 36 477 calving events were recorded. All calves were weighed after birth, and 8.08% of the calves died within 48 h after calving. Calving ease was recorded on a scale between 1 and 4 (1 = easy birth, 4 = surgery), and 69.95, 20.91, 8.92 and 0.21% of the calvings were scored with 1, 2, 3 and 4, respectively. The average inbreeding coefficient of inbred animals was 0.03, and average ancestral inbreeding coefficients were 0.08 and 0.01, depending on how ancestral inbreeding was calculated. Approximately 26% of classically non‐inbred animals showed ancestral inbreeding. Correlations between different inbreeding coefficients ranged between 0.46 and 0.99. No significant effect of ancestral inbreeding was found for calving ease, because the number of animals with reasonable high level of ancestral inbreeding was too low. Significant effects of ancestral inbreeding were estimated for birthweight and stillbirth. Unfavourable effects of ancestral inbreeding were observed for birthweight. However, favourable purging effects were estimated for stillbirth, indicating that purging could be partly beneficial for genetic improvement of stillbirth.     

Detection of APAF1 mutation in Holstein cows and mummified foetuses in Japanese dairy herds

Some of the highest genetic merit sires have been shown to harbour recessive mutations affecting fertility, which may spread rapidly in the population through AI. These disorders may result in abortion and decline in pregnancy per insemination in cows. This study was carried out on 240 Holstein‐Friesian cows and 15 mummified foetuses. Blood and tissue samples were collected from the cows and mummified foetuses, respectively, for DNA extraction. Allele‐specific PCR was designed for the detection of the cows and foetuses carrying the nonsense mutation (C/T) in apoptosis peptide activating factor 1 gene (APAF1). The mutant allele frequency of the APAF1 in carrier cows and mummified foetuses was calculated. Milk samples were taken from the carrier and non‐carrier cows for progesterone assay. The allele‐specific PCR reaction efficiently distinguished the C/T mutation in APAF1. Of 240 cows, seven cows (2.9%) were diagnosed to carry one copy of the mutant allele of APAF1. However, the carrier frequency was 33.3% in mummified foetuses (five of 15). The mutant allele frequency was 0.02 and 0.17 in the cows and mummified foetuses, respectively. Concentrations of progesterone did not differ between cows with APAF1 mutation and non‐carrier cows during 45 days post‐insemination. This study provided allele‐specific PCR for the detection of APAF1 mutation in cows. Moreover, it reports the carrier and mutant allele frequencies of APAF1 in dairy cows and mummified foetuses in Japan.     

Detection of runs of homozygosity in conserved and commercial pig breeds in Poland

Runs of homozygosity (ROH) are continuous segments of the genome that arose as a result of inbreeding, resulting in the inheritance of identical haplotypes from both parents who shared a common ancestor. In the present study, we performed a detailed characterization and comparison of ROH in four pig breeds, including intensively selected Polish Landrace as well as native unselected animals of Puławska and two Złotnicka breeds (White and Spotted). We used a medium-density PorcineSNP60 BeadChip assay (Illumina) and cgaTOH software to detect ROH covering a minimum of 30 adjacent SNPs and maintaining a size over 1 Mb. By analysing ROH distribution and frequency across the genome, we also identified genomic regions with high ROH frequency (so-called “ROH hotspots”). The obtained results showed that the analysed conserved breeds were characterized by a higher ROH span and higher ROH-based inbreeding coefficients (F_ROH), which likely result from past population bottlenecks, increasing the overall inbreeding level within these populations. The analysis of ROH distribution across the genomes revealed the presence of both shared and breed-specific ROH hotspots. These hotspots, presumably representing genome regions under selection, overlapped with a variety of genes associated with processes connected with immune system functioning, reproduction, glucose homeostasis and metabolism. The genome regions with ROH hotspots overlapping in all analysed populations, located on SSC4 (51.9–55.9 Mb) and 13 (92.6–97.8 Mb), covered thirty-one different genes, including MMP16, SLC7A13, ATP6V0D2, CNGB3, WWiP1, RiMDN1 and CPNE3. These genes are primarily associated with biological regulation and metabolism, processes that could be responsible for the variety of the selected production and functional features.     

Reducing inbreeding rates with a breeding circle: Theory and practice in Veluws Heideschaap

Breeding circles allow genetic management in closed populations without pedigrees. In a breeding circle, breeding is split over sub‐populations. Each sub‐population receives breeding males from a single sub‐population and supplies breeding males to one other sub‐population. Donor‐recipient combinations of sub‐populations remain the same over time. Here, we derive inbreeding levels both mathematically and by computer simulation and compare them to actual inbreeding rates derived from DNA information in a real sheep population. In Veluws Heideschaap, a breeding circle has been in operation for over 30 years. Mathematically, starting with inbreeding levels and kinships set to zero, inbreeding rates per generation (ΔF) initially were 0.29%–0.47% within flocks but later converged to 0.18% in all flocks. When, more realistically, inbreeding levels at the start were high and kinship between flocks low, inbreeding levels immediately dropped to the kinship levels between flocks and rates more gradually converged to 0.18%. In computer simulations with overlapping generations, inbreeding levels and rates followed the same pattern, but converged to a lower ΔF of 0.12%. ΔF was determined in the real population with a 12 K SNP chip in recent generations. ΔF in the real population was 0.29%, based on markers to 0.41% per generation based on heterozygosity levels. This is two to three times the theoretically derived values. These increased rates in the real population are probably due to selection and/or the presence of dominant rams siring a disproportionate number of offspring. When these were simulated, ΔF agreed better: 0.35% for selection, 0.38% for dominant rams and 0.67% for both together. The realized inbreeding rates are a warning that in a real population inbreeding rates in a breeding circle can be higher than theoretically expected due to selection and dominant rams. Without a breeding circle, however, inbreeding rates would have been even higher.     

Accuracy of genomic prediction using imputed whole‐genome sequence data in white layers

There is an increasing interest in using whole‐genome sequence data in genomic selection breeding programmes. Prediction of breeding values is expected to be more accurate when whole‐genome sequence is used, because the causal mutations are assumed to be in the data. We performed genomic prediction for the number of eggs in white layers using imputed whole‐genome resequence data including ~4.6 million SNPs. The prediction accuracies based on sequence data were compared with the accuracies from the 60 K SNP panel. Predictions were based on genomic best linear unbiased prediction (GBLUP) as well as a Bayesian variable selection model (BayesC). Moreover, the prediction accuracy from using different types of variants (synonymous, non‐synonymous and non‐coding SNPs) was evaluated. Genomic prediction using the 60 K SNP panel resulted in a prediction accuracy of 0.74 when GBLUP was applied. With sequence data, there was a small increase (~1%) in prediction accuracy over the 60 K genotypes. With both 60 K SNP panel and sequence data, GBLUP slightly outperformed BayesC in predicting the breeding values. Selection of SNPs more likely to affect the phenotype (i.e. non‐synonymous SNPs) did not improve the accuracy of genomic prediction. The fact that sequence data were based on imputation from a small number of sequenced animals may have limited the potential to improve the prediction accuracy. A small reference population (n = 1004) and possible exclusion of many causal SNPs during quality control can be other possible reasons for limited benefit of sequence data. We expect, however, that the limited improvement is because the 60 K SNP panel was already sufficiently dense to accurately determine the relationships between animals in our data.     

Effect of LEPR,ABCG2 and SCD1 Gene Polymorphisms on Reproductive Traits in the Iranian Holstein Cattle

During the last decades, genetic selection for milk production traits has led to increased fertility and health problems in dairy cattle. The aim of this study was to investigate the impact of three polymorphisms located in the ATP‐binding cassette superfamily G member 2 transporter (ABCG2), stearoyl‐CoA desaturase 1 (SCD1) and leptin receptor (LEPR) genes on reproductive traits and somatic cell count (SCC). The analysis was conducted on 408 randomly selected cows. The SNPs within the genes (LEPR, ABCG2 and SCD1) were genotyped using the PCR‐RFLP method. All three possible genotypes were observed for SCD1‐T878C and LEPR‐T945M SNPs, but not for ABCG2‐Y581S SNP. LEPR‐T945M and ABCG2‐Y581S SNPs had no statistically significant effect on the studied reproductive traits and SCC. However, SCD1‐T878C SNP were negatively and significantly related to pregnancy length, dry days and open days (p < 0.05), which lead to decreased profitability in dairy herds. The results suggest that the T878C SNP of SCD1 might be useful as a DNA marker to decrease reproductive problems and improve production traits in Iranian Holstein dairy cows.     

500 years of breeding in the Carthusian Strain of Pura Raza Español horse: An evolutional analysis using genealogical and genomic data

The Carthusian horse is a Pura Raza Español (PRE) strain (CS), bred as a closed population since its creation more than 500 years ago. The aim of this study was to analyse for the first time its population structure and situation of variability combining both genealogical (GEL) and genomic (GEN) data. The GEL data comprised 348,429 pedigree records (56,105 CS horses), while the GEN analysis included the high-density genotypes (670,804 SNPs) of 287 horses. Pedigree completeness demonstrated its accuracy, showing a good correlation of GEL (F) and GEN (F_ROH) inbreeding coefficient in the case of PRE subpopulations partially related and non-related to Carthusian strain (0.68) but a lower value in the 100% Carthusian horses (0.42), due to the high weight of founders not detected by GEL analysis. GEN (PCA, AMOVA, and Admixture) and GEL analysis showed a good differentiation of subpopulations, but also a high level of introgression of the CS in the breed during past decades. A recent change in this trend was noteworthy, with a considerable reduction in CS variability and a genetic bottleneck (effective population sizes of 31.57 and 30.20 in GEL and GEN analysis, respectively, in last generation). The PRE has maintained its variability, and a considerable difference in estimated N_e by GEL (60.77) and GEN (188.0) data was observed. Using two sources of complementary information, it was found the existence of an ancient PRE strain with a unique genetic landmark, practically free from the influence of other equine populations.