Effects of pedigree errors on the efficiency of conservation decisions

Analysing the pedigree to identify undesirable losses of genetic diversity and to prioritize management decisions in captive breeding: a case study

When prevention of species extinction is the priority, captive breeding is a key component in conservation programmes, allowing the recording of pedigree information in studbooks. The genealogical information registered in Cuvier's gazelle studbook between 1975 and 2023 was analysed to (a) assess if the implemented mating policy was successful in preserving the genetic background of the founders (1 male:3 females) in the present population, and b) improve future management and breeding decisions. Although the maternal contribution of one founder female was lost and the mean inbreeding of the total live population was high (0.305 ± 0.095), the breeding policy applied produced better results than expected from a population starting from four founders. It was successful in keeping the individual increase in inbreeding low (0.047 ± 0.021), and, notably, the inbreeding tended to decrease during the last three decades of the breeding programme, ensuring the viability of this highly inbred population. Historical dissemination of individuals among the zoos of Europe and North America caused population structuring and genetic differentiation of the live North American population. However, it did not risk the viability of the captive population. The average relatedness coefficients allowed the identification of individuals with underrepresented genotypes, which is relevant to plan future mating guidelines to keep the founders' representation balanced in the next generations. This study highlights the importance of keeping long-term pedigree information to monitor changes in the genetic diversity of captive populations, which is crucial to implement optimal mating decisions and assuring their long-term viability within an ex situ conservation programme.

Pedigree reconstruction based on genotype data in chickens

A reliable pedigree serves as the backbone of genetic evolution in domesticated animals, providing guidance for daily management and breeding strategies. However, in commercial chicken breeding, pedigree errors and omissions are common. The large-scale application of genomic selection provides an opportunity to reconstruct chicken pedigrees using SNP markers. Here, to reconstruct pedigrees in chickens, we detected high-quality SNPs from 2866 parent-offspring pairs and calculated their genomic relationship and identity by descent (IBD). The results showed that the IBD values for parent-offspring pairs ranged from 0.48 to 0.58, clearly distinguishing them from nonparent-offspring pairs and demonstrating robustness in parentage assignment. In contrast, the genomic relatedness coefficients varied from 0.32 to 0.65. The accuracy of pedigree reconstruction significantly improved as the SNP number and minor allele frequency (MAF) increased. When the number of SNPs exceeded 200, better inference power was exhibited with IBD than with genomic relatedness. Upon reaching an effective SNP quantity of 350, despite a MAF of 0.01, the accuracy of the pedigrees inferred reached a remarkable level of 99%. Furthermore, with a doubled SNP quantity of 700 and a MAF of 0.05, the accuracy increased to a perfect 100%. This study demonstrated the feasibility of accurately constructing pedigrees in chickens using low-density SNP markers and emphasized the importance of considering the number and MAFs of these markers to achieve optimal outcomes. The adoption of the IBD as a suitable metric for pedigree inference is promising for improving the efficiency and accuracy of genetic breeding programs. These findings are paramount for the development of cost-effective yet accurate parentage verification systems.

APIS: an updated parentage assignment software managing triploids induced from diploid parents

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96-192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

Construction of relatedness matrices in autopolyploid populations using low-depth high-throughput sequencing data

KEY MESSAGE

An improved estimator of genomic relatedness using low-depth high-throughput sequencing data for autopolyploids is developed. Its outputs strongly correlate with SNP array-based estimates and are available in the package GUSrelate. High-throughput sequencing (HTS) methods have reduced sequencing costs and resources compared to array-based tools, facilitating the investigation of many non-model polyploid species. One important quantity that can be computed from HTS data is the genetic relatedness between all individuals in a population. However, HTS data are often messy, with multiple sources of errors (i.e. sequencing errors or missing parental alleles) which, if not accounted for, can lead to bias in genomic relatedness estimates. We derive a new estimator for constructing a genomic relationship matrix (GRM) from HTS data for autopolyploid species that accounts for errors associated with low sequencing depths, implemented in the R package GUSrelate. Simulations revealed that GUSrelate performed similarly to existing GRM methods at high depth but reduced bias in self-relatedness estimates when the sequencing depth was low. Using a panel consisting of 351 tetraploid potato genotypes, we found that GUSrelate produced GRMs from genotyping-by-sequencing (GBS) data that were highly correlated with a GRM computed from SNP array data, and less biased than existing methods when benchmarking against the array-based GRM estimates. GUSrelate provides researchers with a tool to reliably construct GRMs from low-depth HTS data.

Evaluation of the Impact of Population Management on the Genetic Parameters of Selected Spiral-Horned Antelopes

The rapid loss of biodiversity and the associated reduction and fragmentation of habitats means that ex situ populations have become an important part of species conservation. These populations, which are often established from a small number of founders, require careful management to avoid the negative effects of genetic drift and inbreeding. Although the inclusion of molecular data is recommended, their availability for captive breeding management remains limited. The aim of this study was to evaluate the relationship between the levels of genetic diversity in six spiral-horned antelope taxa bred under human care and their respective management strategies, conservation status, demography, and geographic origin, using 10 nuclear DNA microsatellite loci and mitochondrial control region DNA sequences. Our findings include associations between genetic diversity and management intensity but also with the diversity and contribution of wild populations to captive founders, with some populations apparently composed of animals from divergent wild lineages elevating captive genetic diversity. When population sizes are large, the potential advantages of maximizing genetic diversity in widely outcrossed populations may need careful consideration with respect to the potential disruption of adaptive diversity. Genetic data serve as a robust tool for managing captive populations, yet their interpretation necessitates a comprehensive understanding of species biology and history.

Genome-Wide Assessment of Runs of Homozygosity by Whole-Genome Sequencing in Diverse Horse Breeds Worldwide

In the genomes of diploid organisms, runs of homozygosity (ROH), consecutive segments of homozygosity, are extended. ROH can be applied to evaluate the inbreeding situation of individuals without pedigree data and to detect selective signatures via ROH islands. We sequenced and analyzed data derived from the whole-genome sequencing of 97 horses, investigated the distribution of genome-wide ROH patterns, and calculated ROH-based inbreeding coefficients for 16 representative horse varieties from around the world. Our findings indicated that both ancient and recent inbreeding occurrences had varying degrees of impact on various horse breeds. However, recent inbreeding events were uncommon, particularly among indigenous horse breeds. Consequently, the ROH-based genomic inbreeding coefficient could aid in monitoring the level of inbreeding. Using the Thoroughbred population as a case study, we discovered 24 ROH islands containing 72 candidate genes associated with artificial selection traits. We found that the candidate genes in Thoroughbreds were involved in neurotransmission (CHRNA6, PRKN, and GRM1), muscle development (ADAMTS15 and QKI), positive regulation of heart rate and heart contraction (HEY2 and TRDN), regulation of insulin secretion (CACNA1S, KCNMB2, and KCNMB3), and spermatogenesis (JAM3, PACRG, and SPATA6L). Our findings provide insight into horse breed characteristics and future breeding strategies.

Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection

Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (F_ROH). The average pedigree inbreeding was 0.07, while the average F_ROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for F_ped and F_ROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.

Comparative analysis of inbreeding parameters and runs of homozygosity islands in 2 Italian autochthonous cattle breeds mainly raised in the Parmigiano-Reggiano cheese production region

Reggiana and Modenese are autochthonous cattle breeds, reared in the North of Italy, that can be mainly distinguished for their standard coat color (Reggiana is red, whereas Modenese is white with some pale gray shades). Almost all milk produced by these breeds is transformed into 2 mono-breed branded Parmigiano-Reggiano cheeses, from which farmers receive the economic incomes needed for the sustainable conservation of these animal genetic resources. After the setting up of their herd books in 1960s, these breeds experienced a strong reduction in the population size that was subsequently reverted starting in the 1990s (Reggiana) or more recently (Modenese) reaching at present a total of about 2,800 and 500 registered cows, respectively. Due to the small population size of these breeds, inbreeding is a very important cause of concern for their conservation programs. Inbreeding is traditionally estimated using pedigree data, which are summarized in an inbreeding coefficient calculated at the individual level (F_PED). However, incompleteness of pedigree information and registration errors can affect the effectiveness of conservation strategies. High-throughput SNP genotyping platforms allow investigation of inbreeding using genome information that can overcome the limits of pedigree data. Several approaches have been proposed to estimate genomic inbreeding, with the use of runs of homozygosity (ROH) considered to be the more appropriate. In this study, several pedigree and genomic inbreeding parameters, calculated using the whole herd book populations or considering genotyping information (GeneSeek GGP Bovine 150K) from 1,684 Reggiana cattle and 323 Modenese cattle, were compared. Average inbreeding values per year were used to calculate effective population size. Reggiana breed had generally lower genomic inbreeding values than Modenese breed. The low correlation between pedigree-based and genomic-based parameters (ranging from 0.187 to 0.195 and 0.319 to 0.323 in the Reggiana and Modenese breeds, respectively) reflected the common problems of local populations in which pedigree records are not complete. The high proportion of short ROH over the total number of ROH indicates no major recent inbreeding events in both breeds. ROH islands spread over the genome of the 2 breeds (15 in Reggiana and 14 in Modenese) identified several signatures of selection. Some of these included genes affecting milk production traits, stature, body conformation traits (with a main ROH island in both breeds on BTA6 containing the ABCG2, NCAPG, and LCORL genes) and coat color (on BTA13 in Modenese containing the ASIP gene). In conclusion, this work provides an extensive comparative analysis of pedigree and genomic inbreeding parameters and relevant genomic information that will be useful in the conservation strategies of these 2 iconic local cattle breeds.

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.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray

Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

Assessment of genomic prediction accuracy using different selection and evaluation approaches in a simulated Korean beef cattle population

OBJECTIVE

This study assessed genomic prediction accuracies based on different selection methods, evaluation procedures, training population (TP) sizes, heritability (h2) levels, marker densities and pedigree error (PE) rates in a simulated Korean beef cattle population.

METHODS

A simulation was performed using two different selection methods, phenotypic and estimated breeding value (EBV), with an h2 of 0.1, 0.3, or 0.5 and marker densities of 10, 50, or 777K. A total of 275 males and 2,475 females were randomly selected from the last generation to simulate ten recent generations. The simulation of the PE dataset was modified using only the EBV method of selection with a marker density of 50K and a heritability of 0.3. The proportions of errors substituted were 10%, 20%, 30%, and 40%, respectively. Genetic evaluations were performed using genomic best linear unbiased prediction (GBLUP) and single-step GBLUP (ssGBLUP) with different weighted values. The accuracies of the predictions were determined.

RESULTS

Compared with phenotypic selection, the results revealed that the prediction accuracies obtained using GBLUP and ssGBLUP increased across heritability levels and TP sizes during EBV selection. However, an increase in the marker density did not yield higher accuracy in either method except when the h2 was 0.3 under the EBV selection method. Based on EBV selection with a heritability of 0.1 and a marker density of 10K, GBLUP and ssGBLUP_0.95 prediction accuracy was higher than that obtained by phenotypic selection. The prediction accuracies from ssGBLUP_0.95 outperformed those from the GBLUP method across all scenarios. When errors were introduced into the pedigree dataset, the prediction accuracies were only minimally influenced across all scenarios.

CONCLUSION

Our study suggests that the use of ssGBLUP_0.95, EBV selection, and low marker density could help improve genetic gains in beef cattle.

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.

Inbreeding, effective population size, and coancestry in the Latxa dairy sheep breed

Traditionally, breeding programs have estimated and managed inbreeding based on pedigree information. The availability of genomic marker panels has made possible new alternatives to achieve more precise estimates, for example in case of missing pedigree. The objective of the present study was to assess and compare, different estimation methods (pedigree-based methodologies, single SNP-based approach (homozygosity) and runs of homozygosity-based method) to analyze the evolution of genetic diversity measured as inbreeding or as coancestry of 3 selected populations of Latxa dairy sheep (Latxa Cara Rubia and Latxa Cara Negra from Euskadi and Navarre). Genomic data came from 972 artificial insemination rams genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) whose genealogy consisted of 4,484 animals. Inbreeding estimates based on molecular data were more similar between them than compared with those based on pedigree information. However, the SNP-based approach estimations of effective population size differed more, reflecting the sensitivity of effective population size to small changes in the evolution of inbreeding. The 2 Latxa Cara Negra populations showed increases of inbreeding rates with time and effective population sizes between 64 and 103 animals, depending on breed and methodology used. The Latxa Cara Rubia population did not show an increase in inbreeding rate, mainly due to semen importation from the related French population of Manech Tête Rousse. The effective size estimates based on coancestry increase show a higher variability and they are more sensitive to the source of information and the data structure considered. Realized effective population size based on individual increase in inbreeding were in agreement with the previous estimates. Coancestry evolution analysis based on DNA information showed an increase on coancestry during the last 10 yr in all breeds, as a consequence of the selection process. Moreover, the increase on coancestry between Latxa Cara Rubia and Manech Tête Rousse was more noticeable between than within each of those breeds.

Assessment of relationships between pigs based on pedigree and genomic information

Relationships play a very important role in studies on quantitative genetics. In traditional breeding, pedigree records are used to establish relationships between animals; while this kind of relationship actually represents one kind of relatedness, it cannot distinguish individual specificity, capture the variation between individuals or determine the actual genetic superiority of an animal. However, with the popularization of high-throughput genotypes, assessments of relationships among animals based on genomic information could be a better option. In this study, we compared the relationships between animals based on pedigree and genomic information from two pig breeding herds with different genetic backgrounds and a simulated dataset. Two different methods were implemented to calculate genomic relationship coefficients and genomic kinship coefficients, respectively. Our results show that, for the same kind of relative, the average genomic relationship coefficients (G matrix) were very close to the pedigree relationship coefficients (A matrix), and on average, the corresponding values were halved in genomic kinship coefficients (K matrix). However, the genomic relationship yielded a larger variation than the pedigree relationship, and the latter was similar to that expected for one relative with no or little variation. Two genomic relationship coefficients were highly correlated, for farm1, farm2 and simulated data, and the correlations for the parent-offspring, full-sib and half-sib were 0.95, 0.90 and 0.85; 0.93, 0.96 and 0.89; and 0.52, 0.85 and 0.77, respectively. When the inbreeding coefficient was measured, the genomic information also yielded a higher inbreeding coefficient and a larger variation than that yielded by the pedigree information. For the two genetically divergent Large White populations, the pedigree relationship coefficients between the individuals were 0, and 62 310 and 175 271 animal pairs in the G matrix and K matrix were greater than 0. Our results demonstrated that genomic information outperformed the pedigree information; it can more accurately reflect the relationships and capture the variation that is not detected by pedigree. This information is very helpful in the estimation of genomic breeding values or gene mapping. In addition, genomic information is useful for pedigree correction. Further, our findings also indicate that genomic information can establish the genetic connection between different groups with different genetic background. In addition, it can be used to provide a more accurate measurement of the inbreeding of an animal, which is very important for the assessment of a population structure and breeding plan. However, the approaches for measuring genomic relationships need further investigation.

Genome-Wide Analysis of SNPs Is Consistent with No Domestic Dog Ancestry in the Endangered Mexican Wolf (Canis lupus baileyi)

The Mexican gray wolf (Canis lupus baileyi) was historically distributed throughout the southwestern United States and northern Mexico. Extensive predator removal campaigns during the early 20th century, however, resulted in its eventual extirpation by the mid 1980s. At this time, the Mexican wolf existed only in 3 separate captive lineages (McBride, Ghost Ranch, and Aragón) descended from 3, 2, and 2 founders, respectively. These lineages were merged in 1995 to increase the available genetic variation, and Mexican wolves were reintroduced into Arizona and New Mexico in 1998. Despite the ongoing management of the Mexican wolf population, it has been suggested that a proportion of the Mexican wolf ancestry may be recently derived from hybridization with domestic dogs. In this study, we genotyped 87 Mexican wolves, including individuals from all 3 captive lineages and cross-lineage wolves, for more than 172000 single nucleotide polymorphisms. We identified levels of genetic variation consistent with the pedigree record and effects of genetic rescue. To identify the potential to detect hybridization with domestic dogs, we compared our Mexican wolf genotypes with those from studies of domestic dogs and other gray wolves. The proportion of Mexican wolf ancestry assigned to domestic dogs was only between 0.06% (SD 0.23%) and 7.8% (SD 1.0%) for global and local ancestry estimates, respectively; and was consistent with simulated levels of incomplete lineage sorting. Overall, our results suggested that Mexican wolves lack biologically significant ancestry with dogs and have useful implications for the conservation and management of this endangered wolf subspecies.

Effect of errors in pedigree on the accuracy of estimated breeding value for carcass traits in Korean Hanwoo cattle

OBJECTIVE

This study evaluated the effect of pedigree errors (PEs) on the accuracy of estimated breeding value (EBV) and genetic gain for carcass traits in Korean Hanwoo cattle.

METHODS

The raw data set was based on the pedigree records of Korean Hanwoo cattle. The animals' information was obtained using Hanwoo registration records from Korean animal improvement association database. The record comprised of 46,704 animals, where the number of the sires used was 1,298 and the dams were 38,366 animals. The traits considered were carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS). Errors were introduced in the pedigree dataset through randomly assigning sires to all progenies. The error rates substituted were 5%, 10%, 20%, 30%, 40%, 50%, 60%, 70%, and 80%, respectively. A simulation was performed to produce a population of 1,650 animals from the pedigree data. A restricted maximum likelihood based animal model was applied to estimate the EBV, accuracy of the EBV, expected genetic gain, variance components, and heritability (h2) estimates for carcass traits. Correlation of the simulated data under PEs was also estimated using Pearson's method.

RESULTS

The results showed that the carcass traits per slaughter year were not consistent. The average CWT, EMA, BFT, and MS were 342.60 kg, 78.76 cm2, 8.63 mm, and 3.31, respectively. When errors were introduced in the pedigree, the accuracy of EBV, genetic gain and h2 of carcass traits was reduced in this study. In addition, the correlation of the simulation was slightly affected under PEs.

CONCLUSION

This study reveals the effect of PEs on the accuracy of EBV and genetic parameters for carcass traits, which provides valuable information for further study in Korean Hanwoo cattle.

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.

The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds

Background

In France, implementation of genomic evaluations in dairy cattle breeds started in 2009 and this has modified the breeding schemes drastically. In this context, the goal of our study was to understand the impact of genomic selection on the genetic diversity of bulls from three French dairy cattle breeds born between 2005 and 2015 (Montbéliarde, Normande and Holstein) and the factors that are involved.

Methods

We compared annual genetic gains, inbreeding rates based on runs of homozygosity (ROH) and pedigree data, and mean ROH length within breeds, before and after the implementation of genomic selection.

Results

Genomic selection induced an increase in mean annual genetic gains of 50, 71 and 33% for Montbéliarde, Normande and Holstein bulls, respectively, and in parallel, the generation intervals were reduced by a factor of 1.7, 1.9 and 2, respectively. We found no significant change in inbreeding rate for the two national breeds, Montbéliarde and Normande, and a significant increase in inbreeding rate for the Holstein international breed, which is now as high as 0.55% per year based on ROH and 0.49% per year based on pedigree data (equivalent to a rate of 1.36 and 1.39% per generation, respectively). The mean ROH length was longer for bulls from the Holstein breed than for those from the other two breeds.

Conclusions

With the implementation of genomic selection, the annual genetic gain increased for bulls from the three major French dairy cattle breeds. At the same time, the annual loss of genetic diversity increased for Holstein bulls, possibly because of the massive use of a few elite bulls in this breed, but not for Montbéliarde and Normande bulls. The increase in mean ROH length in Holstein may reflect the occurrence of recent inbreeding. New strategies in breeding schemes, such as female donor stations and embryo transfer, and recent implementation of genomic evaluations in small regional breeds should be studied carefully in order to ensure the sustainability of breeding schemes in the future.

Advanced optimum contribution selection as a tool to improve regional cattle breeds: a feasibility study for Vorderwald cattle

In the middle of the 20th century, increasing inbreeding rates were identified as a threat to livestock breeding. Consequences include reduced fertility, fitness and phenotypic expression of lethal alleles. An important step in mitigating this inbreeding was the introduction of optimum contribution selection (OCS). OCS facilitates the simultaneous management of genetic gain and inbreeding rates. However, using a standard OCS methodology for regional breeds with historical introgression for upgrading reasons could lead to reinforced selection on introgressed genetic material since those alleles improve the rate of genetic gain and reduce the average kinship in the population. Consequently, regional breeds may become genetically extinct if a standard OCS approach is used. Thus, the advanced OCS (aOCS) approach takes introgressed genetic material into account. The major goals of this study were to (i) gather key information on the feasibility of aOCS under practical conditions of the actual breeding scheme of Vorderwald cattle, (ii) identify superior strategies for implementing the actual scheme and (iii) examine whether historical breeding decisions to increase genetic gain by introgression from commercial breeds could have been avoided by using aOCS. Stochastic simulations were designed in this study to create populations from the historical gene pool by using aOCS. Simultaneously, all practical constraints of a breeding scheme were met. Thus, the simulated populations were comparable with real data. The annual genetic gain was higher in reality (1.56) than in the simulation scenarios (1.12-1.40). The introgressed genetic material increased to 61.3% in reality but was conserved at a final value of 15.3% (±0.78) across simulations. The classical rate of inbreeding and rate of native inbreeding were constrained to 0.092% on an annual basis. This value is equal to an effective population size of 100. The observed values for rates of inbreeding were 0.082-0.087% and 0.087-0.088% for classical and native kinship, respectively. The corresponding figures in reality were 0.067% and 0.184%, respectively. This study suggests that aOCS is feasible for Vorderwald cattle. Strategies for implementation are identified. Finally, we conclude that historical breeding decisions could have been avoided by using aOCS. The genetic gain would have been reduced by at least 12.2%, but the introgressed genetic material, genetic diversity and native genetic diversity would have been more desirable for a breed under conservation.

Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNP-by-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNP-by-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.

Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein-Friesian breeding program from 1986 to 2015

Background

In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015.

Methods

Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time.

Results

Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods.

Conclusions

Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0385-y) contains supplementary material, which is available to authorized users.

Comparative Analysis of Genome Diversity in Bullmastiff Dogs

Management and preservation of genomic diversity in dog breeds is a major objective for maintaining health. The present study was undertaken to characterise genomic diversity in Bullmastiff dogs using both genealogical and molecular analysis. Genealogical analysis of diversity was conducted using a database consisting of 16,378 Bullmastiff pedigrees from year 1980 to 2013. Additionally, a total of 188 Bullmastiff dogs were genotyped using the 170,000 SNP Illumina CanineHD Beadchip. Genealogical parameters revealed a mean inbreeding coefficient of 0.047; 142 total founders (f); an effective number of founders (fe) of 79; an effective number of ancestors (fa) of 62; and an effective population size of the reference population of 41. Genetic diversity and the degree of genome-wide homogeneity within the breed were also investigated using molecular data. Multiple-locus heterozygosity (MLH) was equal to 0.206; runs of homozygosity (ROH) as proportion of the genome, averaged 16.44%; effective population size was 29.1, with an average inbreeding coefficient of 0.035, all estimated using SNP Data. Fine-scale population structure was analysed using NETVIEW, a population analysis pipeline. Visualisation of the high definition network captured relationships among individuals within and between subpopulations. Effects of unequal founder use, and ancestral inbreeding and selection, were evident. While current levels of Bullmastiff heterozygosity, inbreeding and homozygosity are not unusual, a relatively small effective population size indicates that a breeding strategy to reduce the inbreeding rate may be beneficial.

Rates of inbreeding and genetic adaptation for populations managed as herds in zoos with a rotational mating system or with optimized contribution of parents

This study compares two genetic management scenarios for species kept in herds, such as deer. The simulations were designed so that their results can be extended to a wide range of zoo populations. In the first scenario, the simulated populations of size 3 × 20, 6 × 40 or 20 × 60 (herds × animals in herd) were managed with a rotational mating (RM) scheme in which 10%, 20% or 50% of males were selected for breeding and moved between herds in a circular fashion. The second scenario was based on optimal contribution theory (OC). OC requires an accurate pedigree to calculate kinship; males were selected and assigned numbers of offspring to minimize kinship in the next generation. RM was efficient in restriction of inbreeding and produced results comparable with OC. However, RM can result in genetic adaptation of the population to the zoo environment, in particular when 20% or less males are selected for rotation and selection of animals is not random. Lowest rates of inbreeding were obtained by combining OC with rotation of males as in the RM scheme. RM is easy to implement in practice and does not require pedigree data. When full pedigree is available, OC management is preferable.

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 (FIT and FST ) and the effective population size (Ne ) 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 (FROH ) 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 (FROH ) 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.

Interspecific germline transmission of cultured primordial germ cells

In birds, the primordial germ cell (PGC) lineage separates from the soma within 24 h following fertilization. Here we show that the endogenous population of about 200 PGCs from a single chicken embryo can be expanded one million fold in culture. When cultured PGCs are injected into a xenogeneic embryo at an equivalent stage of development, they colonize the testis. At sexual maturity, these donor PGCs undergo spermatogenesis in the xenogeneic host and become functional sperm. Insemination of semen from the xenogeneic host into females from the donor species produces normal offspring from the donor species. In our model system, the donor species is chicken (Gallus domesticus) and the recipient species is guinea fowl (Numida meleagris), a member of a different avian family, suggesting that the mechanisms controlling proliferation of the germline are highly conserved within birds. From a pragmatic perspective, these data are the basis of a novel strategy to produce endangered species of birds using domesticated hosts that are both tractable and fecund.

There is room for selection in a small local pig breed when using optimum contribution selection: a simulation study

Selection progress must be carefully balanced against the conservation of genetic variation in small populations of local breeds. Well-defined breeding programs with specified selection traits are rare in local pig breeds. Given the small population size, the focus is often on the management of genetic diversity. However, in local breeds, optimum contribution selection can be applied to control the rate of inbreeding and to avoid reduced performance in traits with high market value. The aim of this study was to assess the extent to which a breeding program aiming for improved product quality in a small local breed would be feasible. We used stochastic simulations to compare 25 scenarios. The scenarios differed in size of population, selection intensity of boars, type of selection (random selection, truncation selection based on BLUP breeding values, or optimum contribution selection based on BLUP breeding values), and heritability of the selection trait. It was assumed that the local breed is used in an extensive system for a high-meat-quality market. The simulations showed that in the smallest population (300 female reproducers), inbreeding increased by 0.8% when selection was performed at random. With optimum contribution selection, genetic progress can be achieved that is almost as great as that with truncation selection based on BLUP breeding values (0.2 to 0.5 vs. 0.3 to 0.5 genetic SD, P < 0.05), but at a considerably decreased rate of inbreeding (0.7 to 1.2 vs. 2.3 to 5.7%, P < 0.01). This confirmation of the potential utilization of OCS even in small populations is important in the context of sustainable management and the use of animal genetic resources.

MIM: an indirect method to assess inbreeding and coancestry in large incomplete pedigrees of selected dairy cattle

In real data, inbreeding is usually underestimated because of missing pedigree information. A method adapted to the dairy cattle situation is presented to approximate inbreeding when the stored population pedigree is incomplete. Missing parents in incomplete pedigrees were given a dummy identification and assigned to groups (up to nine for a given birth date of progeny). These groups were linked to contemporary reference groups with known parents. An explicit model considered that polygenic breeding values in a censored group were centred on a function of the average breeding value in the corresponding reference group and deviated independently. Inbreeding coefficients were obtained progressively over birth dates starting from founders. For each date considered, the parameters pertaining to its groups were computed using the parameters already obtained from groups belonging to the previous dates. The updating algorithms were given in detail. An indirect method was implemented to expedite mass computations of the relationship coefficients involved (MIM). MIM was compared to Van Raden's (VR) method using simulated populations with 20 overlapping generations and different rates of missing sires and dams. In the situation of random matings, the average inbreeding coefficients by date obtained by MIM were close to true values, whereas they were strongly underestimated by VR. In the situation of assortative matings, MIM gave average inbreeding coefficients moderately underestimated, whereas those of VR's method were still strongly underestimated. The main conclusion of this study adapted to the situation of dairy cattle with incomplete pedigrees was that corrections for inbreeding and coancestry coefficients are more efficient with an explicit appropriate genetic model than without.

History and structure of the closed pedigreed population of Icelandic Sheepdogs

BACKGROUND

Dog breeds lose genetic diversity because of high selection pressure. Breeding policies aim to minimize kinship and therefore maintain genetic diversity. However, policies like mean kinship and optimal contributions might be impractical. Cluster analysis of kinship can elucidate the population structure, since this method divides the population in clusters of related individuals. Kinship-based analyses have been carried out on the entire Icelandic Sheepdog population, a sheep-herding breed.

RESULTS

Analyses showed that despite increasing population size and deliberately transferring dogs, considerable genetic diversity has been lost. When cluster analysis was based on kinships calculated seven generation backwards, as performed in previous studies, results differ markedly from those based on calculations going back to the founder-population, and thus invalidate recommendations based on previous research. When calculated back to the founder-population, kinship-based clustering reveals the distribution of genetic diversity, similarly to strategies using mean kinship.

CONCLUSION

Although the base population consisted of 36 Icelandic Sheepdog founders, the current diversity is equivalent to that of only 2.2 equally contributing founders with no loss of founder alleles in descendants. The maximum attainable diversity is 4.7, unlikely achievable in a non-supervised breeding population like the Icelandic Sheepdog. Cluster analysis of kinship coefficients can provide a supporting tool to assess the distribution of available genetic diversity for captive population management.