A performance programmed method for computing inbreeding coefficients from large data sets for use in mixed-model analyses

Coefficients of inbreeding are commonly used in mixed-model methods for forming inverses of Wright's numerator relationship matrix and transformation matrices used in variance component estimation and national cattle evaluation. Computation of exact coefficients of inbreeding from very large data sets has been believed to be too expensive or too difficult a task to perform. Approximate methods have been used instead. The effects of using approximation methods for inbred data that appear in national cattle data sets are demonstrated. An algorithm is given for the computation of inbreeding coefficients for large data sets. The algorithm feasibly computes inbreeding coefficients for large data sets even on small computing architectures.     

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.     

Selection response in Tribolium to an index modified by introducing between-family or within-family or both restrictions in two steps

A fundamental strategy in selection programs is to combine maximum rate of response and minimum rate of inbreeding, these goals being in conflict with each other. Maximum selection response can be achieved at a cost of erosion in the effective number of breeding animals (a measure of the inbreeding level); reciprocally, the maximum effective number under selection can be preserved with a low response. The simultaneous consideration of both factors makes it difficult to decide on the use of individual (more effective in conserving effective number) or combined selection (maximizes response but yields low effective size). Q uinton et al. (1992) showed that comparing selection methods at the same level of inbreeding, rather than at the same selection intensity, changes the perspectives of current selection theory. If low to moderate inbreeding levels are considered, then phenotypic selection can yield higher response than selection on more accurate methods. Different methods have been proposed for maximizing selection response at the same level of inbreeding, i.e. to restrict the number of close relatives selected (N icholas and S mith 1983), to use false high heritability estimates in the genetic evaluation (G rundy and H ill 1993), to use assortative (S mith and H ammond 1986) or compensatory (G rundy et al. 1994) matings, to adjust estimated breeding values for the relationship with the already selected ones (G oddard and S mith 1990), to avoid matings of related individuals (T oro and P erez -E nciso 1990), or to use factorial rather than hierarchical matings (W oolliams 1989; L eitch et al. 1994). Q uinton and S mith (1995) compared the merits of these methods using stochastic simulation; they concluded that none of the methods was best over all conditions, and that the use of false high heritabilities, or adjusted estimated breeding values with the relationships, does not seem to be recommended; besides, mating together those individuals with the lowest relationship has little effect on the accumulated inbreeding. W ray and G oddard (1994), and B risbane and G ibson (1995) indicated that if G_n is the genetic mean after n generations of selection and F_n is the mean inbreeding coefficient, a reasonable selection objective is G_n ? DF_n, where D is the value of a unit of inbreeding relative to a unit of genetic gain. M euwissen (1997) showed that these methods do not guarantee maximum genetic gains at some level of inbreeding and presented a rule for maximizing the genetic response with a predefined rate of inbreeding. His algorithm can be used to put a constraint on the variance of the selection response by replacing the additive relationship matrix by the prediction error variance (W oolliams and M euwissen 1993). W ei (1995a) developed a restricted phenotypic selection by considering limits on the number of individuals that will be selected from a family or on the family number selected. This less sophisticated method balances response and inbreeding. A restriction on the family number may lead to an increased response (but a decreased effective size), whereas restricting the proportion of selected individuals from a family is an efficient way to control the inbreeding (decreased response). W ei (1995b) generalized the method by introducing both restrictions. In this study, rates of response were compared under between-family, within-family, or both restrictions for a two-trait selection index in a short-term experiment with Tribolium.     

Case study: The effect of inbreeding on the production and reproduction traits in the Elsenburg Dormer sheep stud

Data of the Elsenburg Dormer sheep stud, which was kept closed since inception, were collected over a period of 62 years (1941–2002). The breed is a composite, resulting from a cross of Dorset Horn rams with South African Mutton Merino ewes. These data were analysed to quantify the increase in actual level of inbreeding and to investigate the effect of inbreeding on phenotypic values, genetic parameters and estimated breeding values. After editing 11954 pedigree, 11721 birth weight (BW) and survival, 9205 weaning weight (WW) and 7504 reproduction records were available for analysis. The mean level of inbreeding (F) of all animals over all years was 16%; 14% for dams and 16% for sires. Mean, minimum and maximum F for the lambs in 1997 (when 3 rams from outside were introduced) were 22%, 21% and 24% respectively. Estimates of inbreeding depression for individual inbreeding of 1% were − 0.006 kg for birth and − 0.093 kg for weaning weight respectively. These were the only estimates that were significantly (P < 0.01) different from zero. No significant effects of inbreeding on the other traits were found. There were virtually no differences in the genetic parameters estimated when fitting the two models (inclusion or exclusion of inbreeding coefficients as covariates). Estimates of the phenotypic variance differed slightly for WW between the two models. Ranking of animals were studied for weaning weight when the two models were considered. The high correlation coefficients (0.990) indicate that the use of inbreeding coefficients did not cause important changes in ranking of animals and sires for WW. It was concluded that slow inbreeding (rate of inbreeding of approximately 1.53% per generation over 19 generations) allows natural selection to operate and to remove the less fit animals. At any given mean level of F, less inbreeding depression would then be expected among the individuals who accumulated the inbreeding over a larger number of generations. Nevertheless, inbreeding coefficients should be considered when mating decisions are made, to limit the possible deleterious effects of inbreeding on productive and reproductive traits and to detect animals “resilient to” higher levels of inbreeding.     

Additional considerations to the use of single‐step genomic predictions in a dominance setting

Recent publications indicate that single‐step models are suitable to estimate breeding values, dominance deviations and total genetic values with acceptable quality. Additive single‐step methods implicitly extend known number of allele information from genotyped to non‐genotyped animals. This theory is well derived in an additive setting. It was recently shown, at least empirically, that this basic strategy can be extended to dominance with reasonable prediction quality. Our study addressed two additional issues. It illustrated the theoretical basis for extension and validated genomic predictions to dominance based on single‐step genomic best linear unbiased prediction theory. This development was then extended to include inbreeding into dominance relationships, which is a currently not yet solved issue. Different parametrizations of dominance relationship matrices were proposed. Five dominance single‐step inverse matrices were tested and described as C¹ , C² , C³ , C⁴ and C⁵ . Genotypes were simulated for a real pig population (n = 11,943 animals). In order to avoid any confounding issues with additive effects, pseudo‐records including only dominance deviations and residuals were simulated. SNP effects of heterozygous genotypes were summed up to generate true dominance deviations. We added random noise to those values and used them as phenotypes. Accuracy was defined as correlation between true and predicted dominance deviations. We conducted five replicates and estimated accuracies in three sets: between all ( S₁ ), non‐genotyped ( S₂ ) and inbred non‐genotyped ( S₃ ) animals. Potential bias was assessed by regressing true dominance deviations on predicted values. Matrices accounting for inbreeding ( C³ , C⁴ and C⁵ ) best fit. Accuracies were on average 0.77, 0.40 and 0.46 in S₁ , S₂ and S₃ , respectively. In addition, C³ , C⁴ and C⁵ scenarios have shown better accuracies than C¹ and C² , and dominance deviations were less biased. Better matrix compatibility (accuracy and bias) was observed by re‐scaling diagonal elements to 1 minus the inbreeding coefficient ( C⁵ ).     

Evaluation of inbreeding and genetic diversity in Japanese Shorthorn cattle by pedigree analysis

The Japanese Shorthorn is a Japanese Wagyu breed maintained at a small population size. We assessed the degree of inbreeding and genetic diversity among Japanese Shorthorn cattle using pedigree analysis. We analyzed the pedigree records of registered Japanese Shorthorn born between 1980 and 2018, after evaluating the pedigree completeness. The average of the actual inbreeding coefficients increased at the same rates annually from approximately 1.5% in 1980 to 4.2% in 2018 and was higher than the expected inbreeding coefficients over time. The effective population size based on the individual coancestry rate largely decreased from 127.8 in 1980 to 82.6 in 1999, and then remained almost constant at approximately 90. Three effective numbers of ancestors decreased over time until 1995, then remained almost constant. In particular, the effective number of founder genomes (N_ge) decreased from 43.8 in 1980 to 11.9 in 2018. The index of genetic diversity based on N_ge decreased from 0.99 in 1980 to 0.96 in 2018 due to genetic drift in non-founder generations. Changes in inbreeding and genetic diversity parameters were similar between Japanese Shorthorn and other Japanese Wagyu breeds, but the magnitude of the changes was lower in the Japanese Shorthorn.     

Assessment of inbreeding depression for body measurements in Spanish Purebred (Andalusian) horses

Our aim was to ascertain inbreeding depression in the Spanish Purebred horses for eight body measurements. A total of 16,472 individuals were measured for height at withers, height at chest, leg length, body length, width of chest, heart girth circumference, knee perimeter and cannon bone circumference. Three different multivariate animal models including, respectively, no measure of inbreeding, individual inbreeding coefficients (F_i) or individual increase in inbreeding coefficients (ΔF_i) as linear covariates were used. Significant inbreeding depression was assessed. Even though the models including measures of inbreeding fitted better with data, no effect on estimates of genetic parameters was assessed. However, the inclusion of inbreeding measures affected the ranking order according to the Expected Breeding Values (EBV). Due to the better fit with data and nice properties (the adjustment of individual inbreeding coefficients with the pedigree depth and linear behaviour) the use of ΔF_i in the evaluation models can be recommended for morphological traits in horses.     

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.     

An algorithm to compute optimal genetic contributions in selection programs with large numbers of candidates

A novel algorithm, OCSELECT, is presented for the calculation of optimal genetic contributions with a restricted rate of inbreeding when the number of selection candidates is very large. The calculation of optimal genetic contributions requires the relationship matrix between the candidates and its inverse. The relationship matrix was written as: A = ZA(p)Z' + D, where A(p) is the relationship matrix of the parents, D is a diagonal matrix of Mendelian sampling variances, and Z contains genetic contributions from parents to offspring. Therefore, A(-1) = d(-1) - d(-1)Z(Z'd(-1)Z + A(P)(-1))(-1) Z'd(-1), requires only inversion of matrices of the size of the number of parents instead of the number of offspring. The new algorithm was compared with the software package GENCONT on a salmon data set containing 39,214 selection candidates and 45,846 pedigreed fish in total. Because GENCONT could not handle such a large data set, this data set was split into 19 smaller data sets. Both algorithms gave the same solution with respect to the genetic gain and very similar solutions with respect to the number of selected animals. The OCSELECT algorithm was able to calculate the optimal contributions for the complete data set of 39,214, and therefore no preselection of the 39,214 fish was necessary before entering the fish into the new optimal contribution selection procedure.     

Microsatellite marker analysis of the genetic variability in Hanoverian Hounds

Genetic variability of the dog breed Hanoverian Hound was analysed using a set of 16 microsatellites. The sample of 92 dogs was representative for the total current population [n = 334, inbreeding coefficient 9.2%, relationship coefficient 11.2%] with respect to the level and distribution of the inbreeding and relationship coefficients. All microsatellites used were in Hardy–Weinberg equilibrium. The average number of alleles was 6.4. The average observed heterozygosity (H_O) was slightly higher than the expected heterozygosity (H_E). Dinucleotide microsatellites exhibited lower polymorphism information content (PIC) than tetranucleotide microsatellites (0.52 versus 0.66). The average PIC was 0.61. The individual inbreeding coefficient was negatively related to the average H_O of all microsatellites, whereas the proportion of genes from introducing of Hanoverian Hounds from abroad showed no relationships to H_O. We found that the genetic variability in the Hanoverian Hounds analysed here was unexpectedly higher than that previously published for dog breeds of similar population size. Even in dog breeds of larger population size heterogyzosity was seldom higher than that observed here. The rather high genetic variability as quantified by polymorphic microsatellites in Hanoverian Hounds may be due to a large genetic variation in the founder animals of this breed and to the fact that this genetic diversity could be maintained despite genetic bottlenecks experienced by this breed in the 1920s and 1950s and despite the presence of high inbreeding and relationship coefficients for more than 50 years.     

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.     

Simple preconditioners for the conjugate gradient method: experience with test day models

Preconditioned conjugate gradient method can be used to solve large mixed model equations quickly. Convergence of the method depends on the quality of the preconditioner. Here, the effect of simple preconditioners on the number of iterations until convergence was studied by solving breeding values for several test day models. The test day records were from a field data set, and several simulated data sets with low and high correlations among regression coefficients. The preconditioner matrices had diagonal or block diagonal parts. Transformation of the mixed model equations by diagonalization of the genetic covariance matrix was studied as well. Preconditioner having the whole block of the fixed effects was found to be advantageous. A block diagonal preconditioner for the animal effects reduced the number of iterations the higher the correlations among animal effects, but increased memory usage of the preconditioner. Diagonalization of the animal genetic covariance matrix often reduced the number of iterations considerably without increased memory usage.     

不同分蜂方法的繁蜂效果研究

为了弄清西方蜜蜂在不同方法下进行分蜂后的繁蜂效果和生产性能,对2011年石榴花期前1个月左右已达11脾的蜂群,采取异地和本地分蜂两种方法进行人工分蜂,并对两种分蜂方法的蜂群增殖情况和蜂蜜产量进行了分析。结果表明,第1种分蜂方法优于第2种分蜂方法。第1种分蜂方法分出的蜂繁殖快、进蜜好,是一种理想的分蜂方法。     

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.     

Optimized management of genetic variability in selected pig populations

Controlling the increase of coancestry and inbreeding coefficients in selected populations is made possible through calculation of the optimal contributions allowed to breeding animals, given the current situation with regard to genetic diversity, and further, through optimal design of matings. The potential of such an approach for pig breeding was tested by retrospective optimization on the French Landrace population in reference to the matings actually carried out during a 21-week test period. The major constraint was that the average overall estimated breeding value (EBV) should be the same as the observed one, for not decreasing short-term genetic gain. Optimizing breeding allocations to boars would have led one to decrease coancestry and inbreeding coefficients by approximately 20%. This decrease would have even increased to approximately 30%, would have replacements and disposals been optimized after accounting for genetic variability, keeping the same constraint of genetic level identical to the observed one. These results showed the potential value, in the future, of completing each periodical calculation of EBVs by optimizations considering genetic variability and of releasing corresponding information to breeders, in order to enhance maintenance of genetic variability.     

广西巴马小型猪近交群体的遗传结构分析

本研究以33头广西巴马小型猪为研究对象,利用19个微卫星标记,通过估算其等位基因数(N)、多态信息含量(PIC )和近交系数(F)等参数,对广西巴马小型猪近交群F10~F15世代进行了遗传学检测。试验结果显示,19个微卫星座位上共检测到31个等位基因,其中F10~F15分别具有27、29、28、27、26和25个等位基因,平均每个位点分别仅有1.42、1.53、1.47、1.42、1.37和1.32个等位基因;6个世代群体平均多态信息含量为0.1044。广西巴马小型猪F10~F15世代的平均近交系数分别为0.8070、0.8263、0.8491、0.8710、0.8904、0.9118。结果表明,广西巴马小型猪近交群体的基因多态性低、近交程度高,是一个遗传性稳定的群体。     

Inbreeding depression in Zebu cattle traits

The productivity of herds may be negatively affected by inbreeding depression, and it is important to know how intense is this effect on the livestock performance. We performed a comprehensive analysis involving five Zebu breeds reared in Brazil to estimate inbreeding depression in productive and reproductive traits. Inbreeding depression was estimated for 13 traits by including the individual inbreeding rate as a linear covariate in the standard genetic evaluation models. For all breeds and for almost all traits (no effect was observed on gestation length), the performance of the animals was compromised by an increase in inbreeding. The average inbreeding depression was ?0.222% and ?0.859% per 1% of inbreeding for linear regression coefficients scaled on the percentage of mean (β_m) and standard deviation (β_σ), respectively. The means for β_m (and β_σ) were ?0.269% (?1.202%) for weight/growth traits and ?0.174% (?0.546%) for reproductive traits. Hence, inbreeding depression is more pronounced in weight/growth traits than in reproductive traits. These findings highlight the need for the management of inbreeding in the respective breeding programmes of the breeds studied here.     

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.     

Variance components and breeding values for growth traits from different statistical models.

Estimates of genetic parameters resulting from various analytical models for birth weight (BWT, n = 4,155), 205-d weight (WWT, n = 3,884), and 365-d weight (YWT, n = 3,476) were compared. Data consisted of records for Line 1 Hereford cattle selected for postweaning growth from 1934 to 1989 at ARS-USDA, Miles City, MT. Twelve models were compared. Model 1 included fixed effects of year, sex, age of dam; covariates for birth day and inbreeding coefficients of animal and of dam; and random animal genetic and residual effects. Model 2 was the same as Model 1 but ignored inbreeding coefficients. Model 3 was the same as Model 1 and included random maternal genetic effects with covariance between direct and maternal genetic effects, and maternal permanent environmental effects. Model 4 was the same as Model 3 but ignored inbreeding. Model 5 was the same as Model 1 but with a random sire effect instead of animal genetic effect. Model 6 was the same as Model 5 but ignored inbreeding. Model 7 was a sire model that considered relationships among males. Model 8 was a sire model, assuming sires to be unrelated, but with dam effects as uncorrelated random effects to account for maternal effects. Model 9 was a sire and dam model but with relationships to account for direct and maternal genetic effects; dams also were included as uncorrelated random effects to account for maternal permanent environmental effects. Model 10 was a sire model with maternal grandsire and dam effects all as uncorrelated random effects. Model 11 was a sire and maternal grandsire model, with dams as uncorrelated random effects but with sires and maternal grandsires assumed to be related using male relationships. Model 12 was the same as Model 11 but with all pedigree relationships from the full animal model for sires and maternal grandsires. Rankings on predictions of breeding values were the same regardless of whether inbreeding coefficients for animal and dam were included in the models. Heritability estimates were similar regardless of whether inbreeding effects were in the model. Models 3 and 9 best fit the data for estimation of variances and covariances for direct, maternal genetic, and permanent environmental effects. Other models resulted in changes in ranking for predicted breeding values and for estimates of direct and maternal heritability. Heritability estimates of direct effects were smallest with sire and sire-maternal grandsire models.     

基于简化基因组测序的黄色波尔山羊公羊亲缘关系及近交系数分析

为提高云南黄山羊新品种培育的育种效率和遗传进展，准确度量育种素材波尔山羊黄色群体的近交程度及亲缘关系是一个有效途径。本研究采用简化基因组测序（GBS）技术对来自云南省种羊推广中心的37只黄色波尔山羊公羊进行测序，通过质控获取高质量高密度SNP变异信息，并使用Gmatrix v2、Plink v1.90、MegaX v10.0等软件进行主成分分析（PCA）、状态同源距离计算（identity by state，IBS）、基因型亲缘关系G矩阵构建、亲缘系数计算、NJ聚类分析和基因组近交系数计算。结果显示，在37只黄色公羊中检测出位于29条常染色体上的高质量SNPs位点88 393个，共检测出长纯合片段（ROH）1 537条，大小在1 000.582~18 400.12 kb之间，平均每条ROH长2 576.34 kb，平均含有93.15个SNPs；37只黄色公羊被分为11个家系，其中3个家系仅各有1只黄色公羊，家系A与K亲缘关系最远；基于ROH的群体平均基因组近交系数为0.043，其中3只黄色公羊的基因组近交系数>0.125，存在较多的近交积累。本研究结果为波尔山羊黄色群体在云南黄山羊新品种培育中的合理利用提供了科学依据，也为评估山羊个体近交水平、防止近交衰退、优化选种选配方案提供了有力的技术手段。