Using genomic tools to maintain diversity and fitness in conservation programmes

Genomic analysis of inbreeding and coancestry in Nordic Jersey and Holstein dairy cattle populations

In recent years, genomic selection (GS) has accelerated genetic gain in dairy cattle breeds worldwide. Despite the evident genetic progress, several dairy populations have also encountered challenges such as heightened inbreeding rates and reduced effective population sizes. The challenge has been to find a balance between achieving substantial genetic gain while managing genetic diversity within the population, thereby mitigating the negative effects of inbreeding depression. This study aims to elucidate the impact of GS on pedigree and genomic rates of inbreeding (ΔF) and coancestry (ΔC) in Nordic Jersey (NJ) and Holstein (NH) cattle populations. Furthermore, key genetic metrics, including the generation interval (L), effective population size (Ne), and future effective population size (FNe) were assessed between 2 time periods, before and after GS, and across distinct animal cohorts in both breeds: females, bulls, and approved semen-producing bulls (AI-sires). Analysis of ΔF and ΔC revealed distinct trends across the studied periods and animal groups. Notably, there was a consistent increase in yearly ΔF for most animal groups in both breeds. An exception was observed in NH AI-sires, which demonstrated a slight decrease in yearly ΔF. Moreover, NJ displayed minimal changes in yearly ΔC between the periods, whereas NH exhibited elevated ΔC values across all animal groups. Particularly striking was the substantial increase in yearly ΔC within the NH female population, surging from 0.02% to 0.39% between the periods. Implementation of GS resulted in a reduction of the generation interval across all animal cohorts in both NJ and NH breeds. However, the extent of reduction was more pronounced in males compared with females. This reduction in generation interval influenced generational changes in ΔF and ΔC. Bulls and AI-sires of both breeds exhibited reduced generational ΔF between periods, in contrast to females that demonstrated an opposing pattern. Between the periods, NJ maintained a relatively stable Ne (29.4 before and 30.3 after GS), whereas NH experienced a notable decline from 54.3 to 42.8. Female groups in both breeds displayed a negative Ne trend, whereas males demonstrated either neutral or positive Ne developments. Regarding FNe, NJ exhibited positive FNe development with an increase from 40.7 to 57.2. The opposite was observed in NH, where FNe decreased from 198.8 to 42.7. In summary, it was evident that the genomic methods could detect differences between the populations and changes in ΔF and ΔC more efficiently than pedigree methods. Implementation of GS yielded positive outcomes within the NJ population regarding the rate of coancestry but the opposite was observed with NH. Moreover, analysis of ΔC data hints at the potential to decrease future ΔF through informed mating strategies. Conversely, NH faces more pressing concerns, even though ΔF remains comparatively modest in contrast to what has been observed in other Holstein populations. These findings underscore the necessity of genomic control of inbreeding and coancestry with strategic changes in the Nordic breeding schemes for dairy to ensure long-term sustainability in the forthcoming years.

Mating allocations in Holstein combining genomic information and linear programming optimization at the herd level

In this study, we explored mating allocation in Holstein using genomic information for 24,333 Holstein females born in Denmark, Finland, and Sweden. We used 2 data sets of bulls: the top 50 genotyped bulls and the top 25 polled genotyped bulls on the Nordic total merit scale. We used linear programming to optimize economic scores within each herd, considering genetic level, genetic relationship, semen cost, the economic impact of genetic defects, polledness, and β-casein. We found that it was possible to reduce genetic relationships and eliminate expression of genetic defects with minimal effect on the genetic level in total merit index. Compared with maximizing only Nordic total merit index, the relative frequency of polled offspring increased from 13.5 to 22.5%, and that of offspring homozygous for β-casein (A2A2) from 66.7 to 75.0% in one generation, without any substantial negative impact on other comparison criteria. Using only semen from polled bulls, which might become necessary if dehorning is banned, considerably reduced the genetic level. We also found that animals carrying the polled allele were less likely to be homozygous for β-casein (A2A2) and more likely to be carriers of the genetic defect HH1. Hence, adding economic value to a monogenic trait in the economic score used for mating allocation sometimes negatively affected another monogenetic trait. We recommend that the comparison criteria used in this study be monitored in a modern genomic mating program.

Maintenance of genetic diversity in subdivided populations using genomic coancestry matrices

For both undivided and subdivided populations, the consensus method to maintain genetic diversity is the Optimal Contribution (OC) method. For subdivided populations, this method determines the optimal contribution of each candidate to each subpopulation to maximize global genetic diversity (which implicitly optimizes migration between subpopulations) while balancing the relative levels of coancestry between and within subpopulations. Inbreeding can be controlled by increasing the weight given to within-subpopulation coancestry (λ). Here we extend the original OC method for subdivided populations that used pedigree-based coancestry matrices, to the use of more accurate genomic matrices. Global levels of genetic diversity, measured as expected heterozygosity and allelic diversity, their distributions within and between subpopulations, and the migration pattern between subpopulations, were evaluated via stochastic simulations. The temporal trajectory of allele frequencies was also investigated. The genomic matrices investigated were (i) the matrix based on deviations of the observed number of alleles shared by two individuals from the expected number under Hardy-Weinberg equilibrium; and (ii) a matrix based on a genomic relationship matrix. The matrix based on deviations led to higher global and within-subpopulation expected heterozygosities, lower inbreeding and similar allelic diversity than the second genomic and pedigree-based matrices when a relatively high weight was given to the within-subpopulation coancestries (λ ≥ 5). Under this scenario, allele frequencies moved only slightly away from the initial frequencies. Therefore, the recommended strategy is to use the former matrix in the OC methodology giving a high weight to the within-subpopulation coancestry.

Genomic rotational crossbreeding with advanced optimum contribution selection methods applied to simulated German Angler and German Holstein dairy cattle populations

Many local dairy cattle breeds are facing genetic extinction due to a large proportion of foreign genes, which have been introgressed in the past. In addition, the performance gap to popular high-yielding breeds is increasing, resulting in a risk of numeric extinction. In the present simulation study, a genomic rotational crossbreeding scheme with the high-yielding German Holstein breed and the numerically small German Angler breed was analysed with the aim to utilize heterosis effects in the crossbred animals. Simultaneously inbreeding was controlled, and the amount of Holstein introgression observed in the Angler breed was reduced. Different scenarios of implementing OCS methods for Angler individuals were evaluated, which differed in their restrictions regarding kinship, native kinship, as well as the amount of genetic contributions from German Holstein. The results showed that rotational crossbreeding can result in superior crossbred offspring compared to the purebred parental lines, whereby OCS methods can simultaneously restrict the increase in inbreeding and keep the Holstein contributions at their current level. However, reducing the amount of migrant contributions while restricting the increase in the native kinship in Angler turned out to be a costly restriction. The reason was that Angler with low genetic contributions from Holsteins tended to have similar Angler ancestors. Consequently, reducing Holstein contributions would considerably increase the native kinship in Angler if it were not constrained. The constraint on the native kinship made a constraint on the conventional kinship superfluous and caused it to increase at a much lower rate than envisaged. This led to both, a high genetic diversity and a low genetic gain. The high genetic diversity in Angler also resulted in lower and oscillating heterosis effects in the crossbred animals. Thus, the reduction of migrant contribution did not increase heterosis effects in the crossbred offspring, and did not result in superior crossbred offspring in general.

De novo transcriptomes for two species of snappers from the Western Atlantic (Lutjanus griseus and Lutjanus synagris)

Transcriptomes are appropriate resources for studying species that lack sequenced genomes, as they can serve as references for a broad suite of genetic applications, including: phylogenetic assessments, population genomics, and evaluate responses to environmental fluctuations. Here, we present the transcriptomes of two species of marine fishes of commercial and ecological relevance in the Western Atlantic: Lutjanus griseus and L. synagris. This project represents a step forward on developing genomic resources for important species of the Atlantic Ocean.

The effect of new and ancestral inbreeding on milk production traits in Iranian Holstein cattle

Inbreeding depression, the reduction of fitness and performance, is due to an increase in the mating of related individuals. Based on the purge hypothesis, inbreeding and breeding over generations reduce the effect of deleterious alleles responsible for inbreeding depression. Thus, recent inbreeding is assumed to be more harmful than ancestral inbreeding. This study aimed at evaluating the effects of new and ancestral inbreeding on milk, fat and protein production in Iranian Holstein cattle. The secondary objective was to examine the changes in predicted breeding values when the inbreeding effect was included in the model's analysis. To this end, inbreeding coefficients were calculated using the pedigree of 2,394,517 Holstein cattle to achieve these goals. In addition, 419,132 records of milk, fat and protein yields of first parity cows were collected to assess inbreeding depression and breeding values. The average inbreeding coefficients were 0.83% and 1.68% for the whole population and the inbred animals, respectively. A 1% increase in classical pedigree-based inbreeding coefficient was associated with a decrease of 11.99 kg in milk, 0.39 kg in fat and 0.29 kg in protein. The effect of ancestral inbreeding was more detrimental to performance traits than the effect of new inbreeding. This result contradicted the hypothesis of purging. By including the inbreeding coefficient in the model, the rank of animals remained unchanged, but the average predicted breeding values increased. In general, inbreeding depression was observed in Iranian Holstein cows; however, no evidence of purging was observed. The average of inbreeding coefficients was not high in this population, although accounting for inbreeding coefficients in the analytical model did significantly increase the predicted breeding values. It is recommended that the analytical model incorporate the inbreeding coefficient to improve the accuracy of genetic evaluation. In future studies, inbreeding depression should be assessed using genomic data for performance and reproduction traits.

Genomics-Enabled Management of Genetic Resources in Radiata Pine

Traditional tree improvement is cumbersome and costly. Our main objective was to assess the extent to which genomic data can currently accelerate and improve decision making in this field. We used diameter at breast height (DBH) and wood density (WD) data for 4430 tree genotypes and single-nucleotide polymorphism (SNP) data for 2446 tree genotypes. Pedigree reconstruction was performed using a combination of maximum likelihood parentage assignment and matching based on identity-by-state (IBS) similarity. In addition, we used best linear unbiased prediction (BLUP) methods to predict phenotypes using SNP markers (GBLUP), recorded pedigree information (ABLUP), and single-step “blended” BLUP (HBLUP) combining SNP and pedigree information. We substantially improved the accuracy of pedigree records, resolving the inconsistent parental information of 506 tree genotypes. This led to substantially increased predictive ability (i.e., by up to 87%) in HBLUP analyses compared to a baseline from ABLUP. Genomic prediction was possible across populations and within previously untested families with moderately large training populations (N = 800–1200 tree genotypes) and using as few as 2000–5000 SNP markers. HBLUP was generally more effective than traditional ABLUP approaches, particularly after dealing appropriately with pedigree uncertainties. Our study provides evidence that genome-wide marker data can significantly enhance tree improvement. The operational implementation of genomic selection has started in radiata pine breeding in New Zealand, but further reductions in DNA extraction and genotyping costs may be required to realise the full potential of this approach.

Development and Validation of a 36K SNP Array for Radiata Pine (Pinus radiata D.Don)

Radiata pine (Pinus radiata D.Don) is one of the world’s most domesticated pines and a key economic species in New Zealand. Thus, the development of genomic resources for radiata pine has been a high priority for both research and commercial breeding. Leveraging off a previously developed exome capture panel, we tested the performance of 438,744 single nucleotide polymorphisms (SNPs) on a screening array (NZPRAD01) and then selected 36,285 SNPs for a final genotyping array (NZPRAD02). These SNPs aligned to 15,372 scaffolds from the Pinus taeda L. v. 1.01e assembly, and 20,039 contigs from the radiata pine transcriptome assembly. The genotyping array was tested on more than 8000 samples, including material from archival progenitors, current breeding trials, nursery material, clonal lines, and material from Australia. Our analyses indicate that the array is performing well, with sample call rates greater than 98% and a sample reproducibility of 99.9%. Genotyping in two linkage mapping families indicated that the SNPs are well distributed across the 12 linkage groups. Using genotypic data from this array, we were also able to differentiate representatives of the five recognized provenances of radiata pine, Año Nuevo, Monterey, Cambria, Cedros and Guadalupe. Furthermore, principal component analysis of genotyped trees revealed clear patterns of population structure, with the primary axis of variation driven by provenance ancestry and the secondary axis reflecting breeding activities. This represents the first commercial use of genomics in a radiata pine breeding program.

Mating allocations in Nordic Red Dairy Cattle using genomic information

In this study, we compared mating allocations in Nordic Red Dairy Cattle using genomic information. We used linear programming to optimize different economic scores within each herd, considering genetic level, semen cost, the economic impact of recessive genetic defects, and genetic relationships. We selected 9,841 genotyped females born in Denmark, Finland, or Sweden in 2019 for mating allocations. We used 2 different pedigree relationship coefficients, the first tracing the pedigree 3 generations back from the parents of the potential mating and the second based on all available pedigree information. We used 3 different genomic relationship coefficients, 1 SNP-by-SNP genomic relationship and 2 based on shared genomic segments. We found high correlations (≥0.83) between the pedigree and genomic relationship measures. The mating results showed that it was possible to reduce the different genetic relationships between parents with minimal effect on genetic level. Including the cost of known recessive genetic defects eliminated expression of genetic defects. It was possible to reduce genomic relationships between parents with pedigree measures, but it was best done with genomic measures. Linear programming maximized the economic score for all herds studied within seconds, which means that it is suitable for implementation in mating software to be used by advisors and farmers.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Changes in Allele Frequencies When Different Genomic Coancestry Matrices Are Used for Maintaining Genetic Diversity

A main objective in conservation programs is to maintain genetic variability. This can be achieved using the Optimal Contributions (OC) method that optimizes the contributions of candidates to the next generation by minimizing the global coancestry. However, it has been argued that maintaining allele frequencies is also important. Different genomic coancestry matrices can be used on OC and the choice of the matrix will have an impact not only on the genetic variability maintained, but also on the change in allele frequencies. The objective of this study was to evaluate, through stochastic simulations, the genetic variability maintained and the trajectory of allele frequencies when using two different genomic coancestry matrices in OC to minimize the loss of diversity: (i) the matrix based on deviations of the observed number of alleles shared between two individuals from the expected numbers under Hardy–Weinberg equilibrium (θ_LH); and (ii) the matrix based on VanRaden’s genomic relationship matrix (θ_VR). The results indicate that the use of θ_LH resulted in a higher genetic variability than the use of θ_VR. However, the use of θ_VR maintained allele frequencies closer to those in the base population than the use of θ_LH.

Characterization of Autozygosity in Pigs in Three-Way Crossbreeding

Crossbreeding in livestock can be used to increase genetic diversity. The resulting increase in variability is related to the heterozygosity of the crossbred animal. The evolution of diversity during crossbreeding can be assessed using genomic data. The objective of this study was to describe patterns of runs of homozygosity (ROH) in animals resulting from three-way crossbreeding, from parental pure lines, and in their crossbred offspring. The crossbreeding scheme consisted of a first crossbreeding Pietrain boars and Large White sows, after which the offspring of the Pietrain × Large White were crossed with Duroc boars. The offspring of the second crossbreeding are called G0, the offspring of G0 boars and G0 sows are called G1. All the animals were genotyped using the Illumina SNP60 porcine chip. After filtering, analyses were performed with 2,336 animals and 48,579 autosomal single nucleotide polymorphism (SNP). The mean ROH-based inbreeding coefficients were shown to be 0.27 ± 0.05, 0.23 ± 0.04, and 0.26 ± 0.04 for Duroc, Large White, and Pietrain, respectively. ROH were detected in the Pietrain × Large White crossbred but the homozygous segments were fewer and smaller than in their parents. Similar results were obtained in the G0 crossbred. However, in the G1 crossbreds the number and the size of ROH were higher than in G0 parents. Similar ROH hotspots were detected on SSC1, SSC4, SSC7, SSC9, SSC13, SSC14, and SSC15 in both G0 and G1 animals. Long ROH (>16 Mb) were observed in G1 animals, suggesting regions with low recombination rates. The conservation of these homozygous segments in the three crossbred populations means that some haplotypes were shared between parental breeds. Gene annotation in ROH hotspots in G0 animals identified genes related to production traits including carcass composition and reproduction. These findings advance our understanding of how to manage genetic diversity in crossbred populations.

Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (F_ROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that F_ROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.

Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments

Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.

Evaluating wildlife translocations using genomics: A bighorn sheep case study

Wildlife restoration often involves translocation efforts to reintroduce species and supplement small, fragmented populations. We examined the genomic consequences of bighorn sheep (Ovis canadensis) translocations and population isolation to enhance understanding of evolutionary processes that affect population genetics and inform future restoration strategies. We conducted a population genomic analysis of 511 bighorn sheep from 17 areas, including native and reintroduced populations that received 0-10 translocations. Using the Illumina High Density Ovine array, we generated datasets of 6,155 to 33,289 single nucleotide polymorphisms and completed clustering, population tree, and kinship analyses. Our analyses determined that natural gene flow did not occur between most populations, including two pairs of native herds that had past connectivity. We synthesized genomic evidence across analyses to evaluate 24 different translocation events and detected eight successful reintroductions (i.e., lack of signal for recolonization from nearby populations) and five successful augmentations (i.e., reproductive success of translocated individuals) based on genetic similarity with the source populations. A single native population founded six of the reintroduced herds, suggesting that environmental conditions did not need to match for populations to persist following reintroduction. Augmentations consisting of 18-57 animals including males and females succeeded, whereas augmentations of two males did not result in a detectable genetic signature. Our results provide insight on genomic distinctiveness of native and reintroduced herds, information on the relative success of reintroduction and augmentation efforts and their associated attributes, and guidance to enhance genetic contribution of augmentations and reintroductions to aid in bighorn sheep restoration.

Inbreeding depression across the genome of Dutch Holstein Friesian dairy cattle

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Intensified Use of Reproductive Technologies and Reduced Dimensions of Breeding Schemes Put Genetic Diversity at Risk in Dairy Cattle Breeds

In the management of dairy cattle breeds, two recent trends have arisen that pose potential threats to genetic diversity: the use of reproductive technologies (RT) and a reduction in the number of bulls in breeding schemes. The expected outcome of these changes, in terms of both genetic gain and genetic diversity, is not trivial to predict. Here, we simulated 15 breeding schemes similar to those carried out in large French dairy cattle breeds; breeding schemes differed with respect to their dimensions, the intensity of RT use, and the type of RT involved. We found that intensive use of RT resulted in improved genetic gain, but deteriorated genetic diversity. Specifically, a reduction in the interval between generations through the use of ovum pick-up and in vitro fertilization (OPU-IVF) resulted in a large increase in the inbreeding rate both per year and per generation, suggesting that OPU-IVF could have severe adverse effects on genetic diversity. To achieve a given level of genetic gain, the scenarios that best maintained genetic diversity were those with a higher number of sires/bulls and a medium intensity of RT use or those with a higher number of female donors to compensate for the increased intensity of RT.

RG, Bauck S. Genomic mating as sustainable breeding for Chinese indigenous Ningxiang pigs

An important economic reason for the loss of local breeds is that they tend to be less productive, and hence having less market value than commercial breeds. Nevertheless, local breeds often have irreplaceable values, genetically and sociologically. In the breeding programs with local breeds, it is crucial to balance the selection for genetic gain and the maintaining of genetic diversity. These two objectives are often conflicting, and finding the optimal point of the trade-off has been a challenge for breeders. Genomic selection (GS) provides a revolutionary tool for the genetic improvement of farm animals. At the same time, it can increase inbreeding and produce a more rapid depletion of genetic variability of the selected traits in future generations. Optimum-contribution selection (OCS) represents an approach to maximize genetic gain while constraining inbreeding within a targeted range. In the present study, 515 Ningxiang pigs were genotyped with the Illumina Porcine SNP60 array or the GeneSeek Genomic Profiler Porcine 50K array. The Ningxiang pigs were found to be highly inbred at the genomic level. Average locus-wise inbreeding coefficients were 0.41 and 0.37 for the two SNP arrays used, whereas genomic inbreeding coefficients based on runs of homozygosity were 0.24 and 0.25, respectively. Simulated phenotypic data were used to assess the utility of genomic OCS (GOCS) in comparison with GS without inbreeding control. GOCS was conducted under two scenarios, selecting sires only (GOCS_S) or selecting sires and dams (GOCS_SD), while kinships were constrained on selected parents. The genetic gain for average daily body weight gain (ADG) per generation was between 18.99 and 20.55 g with GOCS_S, and between 23.20 and 28.92 with GOCS_SD, and it varied from 25.38 to 48.38 g under GS without controlling inbreeding. While the rate of genetic gain per generation obtained using GS was substantially larger than that obtained by the two scenarios of genomic OCS in the beginning generations of selection, the difference in the genetic gain of ADG between GS and GOCS reduced quickly in latter generations. At generation ten, the difference in the realized rates of genetic gain between GS and GOCS_SD diminished and ended up with even a slightly higher genetic gain with GOCS_SD, due to the rapid loss of genetic variance with GS and fixation of causative genes. The rate of inbreeding was mostly maintained below 5% per generation with genomic OCS, whereas it increased to between 10.5% and 15.3% per generation with GS. Therefore, genomic OCS appears to be a sustainable strategy for the genetic improvement of local breeds such as Ningxiang pigs, but keeping mind that a variety of GOCS methods exist and the optimal forms remain to be exploited further.

Management of Genetic Diversity in the Era of Genomics

Management of genetic diversity aims to (i) maintain heterozygosity, which ameliorates inbreeding depression and loss of genetic variation at loci that may become of importance in the future; and (ii) avoid genetic drift, which prevents deleterious recessives (e.g., rare disease alleles) from drifting to high frequency, and prevents random drift of (functional) traits. In the genomics era, genomics data allow for many alternative measures of inbreeding and genomic relationships. Genomic relationships/inbreeding can be classified into (i) homozygosity/heterozygosity based (e.g., molecular kinship matrix); (ii) genetic drift-based, i.e., changes of allele frequencies; or (iii) IBD-based, i.e., SNPs are used in linkage analyses to identify IBD segments. Here, alternative measures of inbreeding/relationship were used to manage genetic diversity in genomic optimal contribution (GOC) selection schemes. Contrary to classic inbreeding theory, it was found that drift and homozygosity-based inbreeding could differ substantially in GOC schemes unless diversity management was based upon IBD. When using a homozygosity-based measure of relationship, the inbreeding management resulted in allele frequency changes toward 0.5 giving a low rate of increase in homozygosity for the panel used for management, but not for unmanaged neutral loci, at the expense of a high genetic drift. When genomic relationship matrices were based on drift, following VanRaden and as in GCTA, drift was low at the expense of a high rate of increase in homozygosity. The use of IBD-based relationship matrices for inbreeding management limited both drift and the homozygosity-based rate of inbreeding to their target values. Genetic improvement per percent of inbreeding was highest when GOC used IBD-based relationships irrespective of the inbreeding measure used. Genomic relationships based on runs of homozygosity resulted in very high initial improvement per percent of inbreeding, but also in substantial discrepancies between drift and homozygosity-based rates of inbreeding, and resulted in a drift that exceeded its target value. The discrepancy between drift and homozygosity-based rates of inbreeding was caused by a covariance between initial allele frequency and the subsequent change in frequency, which becomes stronger when using data from whole genome sequence.

Assessing the power of principal components and wright's fixation index analyzes applied to reveal the genome-wide genetic differences between herds of Holstein cows

Background

Due to the advent of SNP array technology, a genome-wide analysis of genetic differences between populations and breeds has become possible at a previously unattainable level. The Wright’s fixation index (F_st) and the principal component analysis (PCA) are widely used methods in animal genetics studies. In paper we compared the power of these methods, their complementing each other and which of them is the most powerful.

Results

Comparative analysis of the power Principal Components Analysis (PCA) and F_st were carried out to reveal genetic differences between herds of Holsteinized cows. Totally, 803 BovineSNP50 genotypes of cows from 13 herds were used in current study. Obtained F_st values were in the range of 0.002–0.012 (mean 0.0049) while for rare SNPs with MAF 0.0001–0.005 they were even smaller in the range of 0.001–0.01 (mean 0.0027). Genetic relatedness of the cows in the herds was the cause of such small F_st values. The contribution of rare alleles with MAF 0.0001–0.01 to the F_st values was much less than common alleles and this effect depends on linkage disequilibrium (LD). Despite of substantial change in the MAF spectrum and the number of SNPs we observed small effect size of LD - based pruning on F_st data. PCA analysis confirmed the mutual admixture and small genetic difference between herds. Moreover, PCA analysis of the herds based on the visualization the results of a single eigenvector cannot be used to significantly differentiate herds. Only summed eigenvectors should be used to realize full power of PCA to differentiate small between herds genetic difference. Finally, we presented evidences that the significance of F_st data far exceeds the significance of PCA data when these methods are used to reveal genetic differences between herds.

Conclusions

LD - based pruning had a small effect on findings of F_st and PCA analyzes. Therefore, for weakly structured populations the LD - based pruning is not effective. In addition, our results show that the significance of genetic differences between herds obtained by F_st analysis exceeds the values of PCA. Proposed, to differentiate herds or low structured populations we recommend primarily using the F_st approach and only then PCA.

Effect of genomic selection on rate of inbreeding and coancestry and effective population size of Holstein and Jersey cattle populations

Genetic diversity in livestock populations is a significant contributor to the sustainability of animal production. Also, genetic diversity allows animal production to become more responsive to environmental changes and market demands. The loss of genetic diversity can result in a plateau in production and may also result in loss of fitness or viability in animal production. In this study, we investigated the rate of inbreeding (ΔF), rate of coancestry (Δf), and effective population size (N_e) as important quantitative indicators of genetic diversity and evaluated the effect of the recent implementation of genomic selection on the loss of genetic diversity in North American Holstein and Jersey dairy cattle. To estimate the rate of inbreeding and coancestry, inbreeding and coancestry coefficients were calculated using the traditional pedigree method and genomic methods estimated from segment- and marker-based approaches. Furthermore, we estimated N_e from the rate of inbreeding and coancestry and extent of linkage disequilibrium. A total of 205,755 and 89,238 pedigreed and genotyped animals born between 1990 and 2018 inclusively were available for Holsteins and Jerseys, respectively. The estimated average pedigree inbreeding coefficients were 7.74 and 7.20% for Holsteins and Jerseys, respectively. The corresponding values for the segment and marker-by-marker genomic inbreeding coefficients were 13.61, 15.64, and 31.40% for Holsteins and 21.16, 22.54, and 42.62% for Jerseys, respectively. The average coancestry coefficients were 8.33 and 15.84% for Holsteins and 9.23 and 23.46% for Jerseys with pedigree and genomic measures, respectively. Generation interval for the whole 29-yr time period averaged approximately 5 yr for all selection pathways combined. The ΔF per generation based on pedigree, segment, and marker-by-marker genomic measures for the entire 29-yr period was estimated to be 0.75, 1.10, 1.16, and 1.02% for Holstein animals and 0.67, 0.62, 0.63, and 0.59% for Jersey animals, respectively. The Δf was estimated to be 0.98 and 0.98% for Holsteins and 0.73 and 0.78% for Jerseys with pedigree and genomic measures, respectively. These ΔF and Δf translated to an N_e that ranged from 43 to 66 animals for Holsteins and 64 to 85 animals for Jerseys. In addition, the N_e based on linkage disequilibrium was 58 and 120 for Holsteins and Jerseys, respectively. The 10-yr period that involved the application of genomic selection resulted in an increased ΔF per generation with ranges from 1.19 to 2.06% for pedigree and genomic measures in Holsteins. Given the rate at which inbreeding is increasing after the implementation of genomic selection, there is a need to implement measures and means for controlling the rate of inbreeding per year, which will help to manage and maintain farm animal genetic resources.

Genome-Wide Runs of Homozygosity, Effective Population Size, and Detection of Positive Selection Signatures in Six Chinese Goat Breeds

Detection of selection footprints provides insight into the evolution process and the underlying mechanisms controlling the phenotypic diversity of traits that have been exposed to selection. Selection focused on certain characters, mapping certain genomic regions often shows a loss of genetic diversity with an increased level of homozygosity. Therefore, the runs of homozygosity (ROHs), homozygosity by descent (HBD), and effective population size (N_e) are effective tools for exploring the genetic diversity, understanding the demographic history, foretelling the signature of directional selection, and improving the breeding strategies to use and conserve genetic resources. We characterized the ROH, HBD, N_e, and signature of selection of six Chinese goat populations using single nucleotide polymorphism (SNP) 50K Illumina beadchips. Our results show an inverse relationship between the length and frequency of ROH. A long ROH length, higher level of inbreeding, long HBD segment, and smaller N_e in Guangfeng (GF) goats suggested intensive selection pressure and recent inbreeding in this breed. We identified six reproduction-related genes within the genomic regions with a high ROH frequency, of which two genes overlapped with a putative selection signature. The estimated pair-wise genetic differentiation (F_ST) among the populations is 9.60% and the inter- and intra-population molecular variations are 9.68% and 89.6%, respectively, indicating low to moderate genetic differentiation. Our selection signatures analysis revealed 54 loci harboring 86 putative candidate genes, with a strong signature of selection. Further analysis showed that several candidate genes, including MARF1, SYCP2, TMEM200C, SF1, ADCY1, and BMP5, are involved in goat fecundity. We identified 11 candidate genes by using cross-population extended haplotype homozygosity (XP-EHH) estimates, of which MARF1 and SF1 are under strong positive selection, as they are differentiated in high and low reproduction groups according to the three approaches used. Gene ontology enrichment analysis revealed that different biological pathways could be involved in the variation of fecundity in female goats. This study provides a new insight into the ROHs patterns for maintenance of within breed diversity and suggests a role of positive selection for genetic variation influencing fecundity in Chinese goat.

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.

Key Genetic Parameters for Population Management

Population management has the primary task of maximizing the long-term competitiveness of a breed. Breeds compete with each other for being able to supply consumer demands at low costs and also for funds from conservation programs. The competition for consumer preference is won by breeds with high genetic gain for total merit who maintained a sufficiently high genetic diversity, whereas the competition for funds is won by breeds with high conservation value. The conservation value of a breed could be improved by increasing its contribution to the gene pool of the species. This may include the recovery of its original genetic background and the maintenance of a high genetic diversity at native haplotype segments. The primary objective of a breeding program depends on the genetic state of the population and its intended usage. In this paper, we review the key genetic parameters that are relevant for population management, compare the methods for estimating them, derive the formulas for predicting their value at a future time, and clarify their usage in various types of breeding programs that differ in their main objectives. These key parameters are kinships, native kinships, breeding values, Mendelian sampling variances, native contributions, and mutational effects. Population management currently experiences a transition from using pedigree-based estimates to marker-based estimates, which improves the accuracies of these estimates and thereby increases response to selection. In addition, improved measures of the factors that determine the competitiveness of a breed and utilize auxiliary parameters, such as Mendelian sampling variances, mutational effects, and native kinships, enable to improve further upon historic recommendations for genetic population management.

A double-labeling marker-based method for estimating inbreeding and parental genomic components in a population under conservation

Objective

The objective of a conservation program is to maintain maximum genetic diversity and preserve the viability of a breed. However, the efficiency of a program is influenced by the ability to accurately measure and predict genetic diversity.

Methods

To examine this question, we conducted a simulation in which common measures (i.e. heterozygosity) and novel measures (identity-by-descent probabilities and parental genomic components) were used to estimate genetic diversity within a conserved population using double-labeled single nucleotide polymorphism markers.

Results

The results showed that the accuracy and sensitivity of identity-by-state probabilities and heterozygosity were close to identity by descent (IBD) probabilities, which reflect the true genetic diversity. Expected heterozygosity most closely aligned with IBD. All common measures suggested that practices used in the current Chinese pig conservation program result in a ~5% loss in genetic diversity every 10 generations. Parental genomic components were also analyzed to monitor real-time changes in genomic components for each male and female ancestor. The analysis showed that ~7.5% of male families and ~30% of female families were lost every 5 generations. After 50 generations of simulated conservation, 4 male families lost ~50% of their initial genomic components, and the genomic components for 24.8% of the female families were lost entirely.

Conclusion

In summary, compared with the true genetic diversity value obtained using double-labeled markers, expected heterozygosity appears to be the optimal indicator. Parental genomic components analysis provides a more detailed picture of genetic diversity and can be used to guide conservation management practices.

Managing genomes of conserved livestock breeds with historical introgression to decrease genetic overlap with other breeds

Recovering the native genetic background of a breed and increasing the founder genome equivalent (FGE) that is contributed by a breed to the gene pool of the species can increase its value for conservation. The suitability of several strategies was compared, whereby a hypothetical multi-breed population, the core set, was used to approximate the genetic diversity of the species. Twenty-five generations of management were simulated based on genotypes of German Angler cattle. The scenarios were compared when the kinship reached 0.10. The native contribution (NC) increased in a population with 400 births per generation from 0.317 to 0.706, whereas 1,000 births enabled to reach 0.894. This scenario maximized the NC, constrained the native kinship, and the kinship of the core set so that its genetic diversity could not decrease. It increased the proportions of mainstream breeds because their genes were removed from the target breed. A substantial increase of the FGE was achieved in some other scenarios, which arose from reduced genetic overlap and from increased diversity within the breed. The latter factor is especially important for breeds with high contributions to the core set.

Deleterious alleles in the context of domestication, inbreeding, and selection

Each individual has a certain number of harmful mutations in its genome. These mutations can lower the fitness of the individual carrying them, dependent on their dominance and selection coefficient. Effective population size, selection, and admixture are known to affect the occurrence of such mutations in a population. The relative roles of demography and selection are a key in understanding the process of adaptation. These are factors that are potentially influenced and confounded in domestic animals. Here, we hypothesize that the series of events of bottlenecks, introgression, and strong artificial selection associated with domestication increased mutational load in domestic species. Yet, mutational load is hard to quantify, so there are very few studies available revealing the relevance of evolutionary processes. The precise role of artificial selection, bottlenecks, and introgression in further increasing the load of deleterious variants in animals in breeding and conservation programmes remains unclear. In this paper, we review the effects of domestication and selection on mutational load in domestic species. Moreover, we test some hypotheses on higher mutational load due to domestication and selective sweeps using sequence data from commercial pig and chicken lines. Overall, we argue that domestication by itself is not a prerequisite for genetic erosion, indicating that fitness potential does not need to decline. Rather, mutational load in domestic species can be influenced by many factors, but consistent or strong trends are not yet clear. However, methods emerging from molecular genetics allow discrimination of hypotheses about the determinants of mutational load, such as effective population size, inbreeding, and selection, in domestic systems. These findings make us rethink the effect of our current breeding schemes on fitness of populations.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits

Background

The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (F_ROH), genomic relationship matrix (F_GRM), and pedigree-based coefficient (F_PED).

Results

The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of F_PED and F_ROH, and their correlations (− 0.05 to 0.26) were low. Estimates of correlation between F_GRM-F_PED was zero, while the correlation (− 0.01 to − 0.07) between F_GRM-F_ROH decreased as a function of ROH length, except for F_{ROH > 8Mb} (− 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands.

Conclusions

Low F_PED-F_ROH correlation estimates indicate that F_PED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and F_ROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5060-8) contains supplementary material, which is available to authorized users.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Evaluating sample size to estimate genetic management metrics in the genomics era

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free-ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross-species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra- and interpopulation relationships using kinship and identity by state metrics, as well as FST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.

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.

Population genetic analysis of the recently rediscovered Hula painted frog (Latonia nigriventer) reveals high genetic diversity and low inbreeding

After its recent rediscovery, the Hula painted frog (Latonia nigriventer) has remained one of the world's rarest and least understood amphibian species. Together with its apparently low dispersal capability and highly disturbed niche, the low abundance of this living fossil calls for urgent conservation measures. We used 18 newly developed microsatellite loci and four different models to calculate the effective population size (N_e) of a total of 125 Hula painted frog individuals sampled at a single location. We compare the N_e estimates to the estimates of potentially reproducing adults in this population (N_ad) determined through a capture-recapture study on 118 adult Hula painted frogs captured at the same site. Surprisingly, our data suggests that, despite N_ad estimates of only ~234-244 and N_e estimates of ~16.6-35.8, the species appears to maintain a very high genetic diversity (H_O = 0.771) and low inbreeding coefficient (F_IS = -0.018). This puzzling outcome could perhaps be explained by the hypotheses of either genetic rescue from one or more unknown Hula painted frog populations nearby or by recent admixture of genetically divergent subpopulations. Independent of which scenario is correct, the original locations of these populations still remain to be determined.

Genomic Characterisation of the Indigenous Irish Kerry Cattle Breed

Kerry cattle are an endangered landrace heritage breed of cultural importance to Ireland. In the present study we have used genome-wide SNP array data to evaluate genomic diversity within the Kerry population and between Kerry cattle and other European breeds. Patterns of genetic differentiation and gene flow among breeds using phylogenetic trees with ancestry graphs highlighted historical gene flow from the British Shorthorn breed into the ancestral population of modern Kerry cattle. Principal component analysis (PCA) and genetic clustering emphasised the genetic distinctiveness of Kerry cattle relative to comparator British and European cattle breeds. Modelling of genetic effective population size (N_e) revealed a demographic trend of diminishing N_e over time and that recent estimated N_e values for the Kerry breed may be less than the threshold for sustainable genetic conservation. In addition, analysis of genome-wide autozygosity (F_ROH) showed that genomic inbreeding has increased significantly during the 20 years between 1992 and 2012. Finally, signatures of selection revealed genomic regions subject to natural and artificial selection as Kerry cattle adapted to the climate, physical geography and agro-ecology of southwest Ireland.

Which Individuals To Choose To Update the Reference Population? Minimizing the Loss of Genetic Diversity in Animal Genomic Selection Programs

Genomic selection (GS) is commonly used in livestock and increasingly in plant breeding. Relying on phenotypes and genotypes of a reference population, GS allows performance prediction for young individuals having only genotypes. This is expected to achieve fast high genetic gain but with a potential loss of genetic diversity. Existing methods to conserve genetic diversity depend mostly on the choice of the breeding individuals. In this study, we propose a modification of the reference population composition to mitigate diversity loss. Since the high cost of phenotyping is the limiting factor for GS, our findings are of major economic interest. This study aims to answer the following questions: how would decisions on the reference population affect the breeding population, and how to best select individuals to update the reference population and balance maximizing genetic gain and minimizing loss of genetic diversity? We investigated three updating strategies for the reference population: random, truncation, and optimal contribution (OC) strategies. OC maximizes genetic merit for a fixed loss of genetic diversity. A French Montbéliarde dairy cattle population with 50K SNP chip genotypes and simulations over 10 generations were used to compare these different strategies using milk production as the trait of interest. Candidates were selected to update the reference population. Prediction bias and both genetic merit and diversity were measured. Changes in the reference population composition slightly affected the breeding population. Optimal contribution strategy appeared to be an acceptable compromise to maintain both genetic gain and diversity in the reference and the breeding populations.

Long-Term Impact of Optimum Contribution Selection Strategies on Local Livestock Breeds with Historical Introgression Using the Example of German Angler Cattle

The long-term performance of different selection strategies was evaluated via simulation using the example of a local cattle breed, German Angler cattle. Different optimum contribution selection (OCS) approaches to maximize genetic gain were compared to a reference scenario without selection and truncation selection. The kinships and migrant contribution (MC) were estimated from genomic data. Truncation selection achieved the highest genetic gain but decreased diversity considerably at native alleles. It also caused the highest increase in MCs. Traditional OCS, which only constrains kinship, achieved almost the same genetic gain but also caused a small increase of MC and remarkably reduced the diversity of native alleles. When MC was required not to increase and the increase of kinship at native alleles was restricted, the MC levels and the diversity at native alleles were well managed, and the genetic gain was only slightly reduced. However, genetic progress was substantially lower in the scenario that aimed to recover the original genetic background. Truncation selection and traditional OCS selection both reduce the genetic originality of breeds with historical introgression. The inclusion of MC and kinship at native alleles as additional constraints in OCS showed great potential for conservation. Recovery of the original genetic background is possible but requires many generations of selection and reduces the genetic progress in performance traits. Hence, constraining MCs at their current values can be recommended to avoid further reduction of genetic originality.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

Applying SNP-Derived Molecular Coancestry Estimates to Captive Breeding Programs

Captive breeding programs for wildlife species typically rely on pedigrees to inform genetic management. Although pedigree-based breeding strategies are quite effective at retaining long-term genetic variation, management of zoo-based breeding programs continues to be hampered when pedigrees are poorly known. The objective of this study was to evaluate 2 options for generating single nucleotide polymorphism (SNP) data to resolve unknown relationships within captive breeding programs. We generated SNP data for a zoo-based population of addax (Addax nasomasculatus) using both the Illumina BovineHD BeadChip and double digest restriction site-associated DNA (ddRAD) sequencing. Our results demonstrated that estimates of allele sharing (AS) between pairs of individuals exhibited low variances. Average AS variances were highest when using 50 loci (SNPchipall = 0.00159; ddRADall = 0.0249), but fell below 0.0003 for the SNP chip dataset when sampling ≥250 loci and below 0.0025 for the ddRAD dataset when sampling ≥500 loci. Furthermore, the correlation between the SNPchipall and ddRADall AS datasets was 0.88 (95%CI = 0.84-0.91) when subsampling 500 loci. Collectively, our results indicated that both SNP genotyping methods produced sufficient data for accurately estimating relationships, even within an extremely bottlenecked population. Our results also suggested that analytic assumptions historically integrated into the addax pedigree are not adversely impacting long-term pedigree-based management; kinships calculated from the analytic pedigree were significantly correlated (P << 0.001) with AS estimates. Overall, our conclusions are intended to serve as both a proof of concept and a model for applying molecular data to the genetic management of captive breeding programs.

The Use of Genomics in Conservation Management of the Endangered Visayan Warty Pig (Sus cebifrons)

The list of threatened and endangered species is growing rapidly, due to various anthropogenic causes. Many endangered species are present in captivity and actively managed in breeding programs in which often little is known about the founder individuals. Recent developments in genetic research techniques have made it possible to sequence and study whole genomes. In this study we used the critically endangered Visayan warty pig (Sus cebifrons) as a case study to test the use of genomic information as a tool in conservation management. Two captive populations of S. cebifrons exist, which originated from two different Philippine islands. We found some evidence for a recent split between the two island populations; however all individuals that were sequenced show a similar demographic history. Evidence for both past and recent inbreeding indicated that the founders were at least to some extent related. Together with this, the low level of nucleotide diversity compared to other Sus species potentially poses a threat to the viability of the captive populations. In conclusion, genomic techniques answered some important questions about this critically endangered mammal and can be a valuable toolset to inform future conservation management in other species as well.

The use of coancestry based on shared segments for maintaining genetic diversity

We have evaluated the use of genomic coancestry coefficients based on shared segments for the maintenance of genetic diversity through optimal contributions methodology for populations of three different Austrian cattle breeds. This coancestry measure has been compared with the genomic coancestry coefficient calculated on a SNP-by-SNP basis and with pedigree-based coancestry. The regressions of the shared segments coancestry on the other two coefficients suggest that the former mainly reflect Identity By Descent but with the advantage over pedigree-based coancestry of providing the realized Identity By Descent rather than an expectation. The effective population size estimated from the rate of coancestry based on shared segments was very similar to those obtained with the other coefficients and of small magnitude (from 26.24 to 111.90). This result highlights the importance of implementing active management strategies to control the increase of inbreeding and the loss of genetic diversity in livestock breeds, even when the population size is reasonably large. One problem for the implementation of coancestry based on shared segments is the need of estimating the gametic phases of the SNPs which, given the techniques used to obtain the genotypes, are a priori unknown. This study shows, through computer simulations, that using estimates of gametic phases for computing coancestry based on shared segments does not lead to a significant loss in the diversity maintained. This has been shown to be true even when the size of the population is very small as it is usually the case in populations subjected to conservation programmes.

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F ROH) in three local dairy cattle breeds. F ROH values were compared with F estimated from the genomic relationship matrix (F GRM), based on the difference between observed v. expected number of homozygous genotypes (F HOM) and the genomic homozygosity of individual i (F MOL i ). The molecular coancestry coefficient (f MOL ij ) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F HOM and F ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F ROH and other estimated F coefficients were generally lower ranged from 0.45 (F MOL i -F ROH) in Cinisara to 0.17 (F GRM-F ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.

Very Low Population Structure in a Highly Mobile and Wide-Ranging Endangered Bird Species

The loss of biodiversity following fragmentation and degradation of habitat is a major issue in conservation biology. As competition for resources increases following habitat loss and fragmentation, severe population declines may occur even in common, highly mobile species; such demographic decline may cause changes within the population structure of the species. The regent honeyeater, Anthochaera phrygia, is a highly nomadic woodland bird once common in its native southeast Australia. It has experienced a sharp decline in abundance since the late 1970s, following clearing of large areas of its preferred habitat, box-ironbark woodland, within the last 200 years. A captive breeding program has been established as part of efforts to restore this species. This study used genetic data to examine the range-wide population structure of regent honeyeaters, including spatial structure, its change through time, sex differences in philopatry and mobility, and genetic differences between the captive and wild populations. There was low genetic differentiation between birds captured in different geographic areas. Despite the recent demographic decline, low spatial structure appears to have some temporal consistency. Both sexes appear to be highly mobile, and there does not seem to be significant genetic differentiation between the captive and wild populations. We conclude that management efforts for survival of this species, including habitat protection, restoration, and release of captive-bred birds into the wild, can treat the species as effectively a single genetic population.

Using genome-wide measures of coancestry to maintain diversity and fitness in endangered and domestic pig populations

Conservation and breeding programs aim at maintaining the most diversity, thereby avoiding deleterious effects of inbreeding while maintaining enough variation from which traits of interest can be selected. Theoretically, the most diversity is maintained using optimal contributions based on many markers to calculate coancestries, but this can decrease fitness by maintaining linked deleterious variants. The heterogeneous patterns of coancestry displayed in pigs make them an excellent model to test these predictions. We propose methods to measure coancestry and fitness from resequencing data and use them in population management. We analyzed the resequencing data of Sus cebifrons, a highly endangered porcine species from the Philippines, and genotype data from the Pietrain domestic breed. By analyzing the demographic history of Sus cebifrons, we inferred two past bottlenecks that resulted in some inbreeding load. In Pietrain, we analyzed signatures of selection possibly associated with commercial traits. We also simulated the management of each population to assess the performance of different optimal contribution methods to maintain diversity, fitness, and selection signatures. Maximum genetic diversity was maintained using marker-by-marker coancestry, and least using genealogical coancestry. Using a measure of coancestry based on shared segments of the genome achieved the best results in terms of diversity and fitness. However, this segment-based management eliminated signatures of selection. We demonstrate that maintaining both diversity and fitness depends on the genomic distribution of deleterious variants, which is shaped by demographic and selection histories. Our findings show the importance of genomic and next-generation sequencing information in the optimal design of breeding or conservation programs.

Genome-wide estimates of coancestry, inbreeding and effective population size in the Spanish Holstein population

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.

Artificial selection with traditional or genomic relationships: consequences in coancestry and genetic diversity

Estimated breeding values (EBVs) are traditionally obtained from pedigree information. However, EBVs from high-density genotypes can have higher accuracy than EBVs from pedigree information. At the same time, it has been shown that EBVs from genomic data lead to lower increases in inbreeding compared with traditional selection based on genealogies. Here we evaluate the performance with BLUP selection based on genealogical coancestry with three different genome-based coancestry estimates: (1) an estimate based on shared segments of homozygosity, (2) an approach based on SNP-by-SNP count corrected by allelic frequencies, and (3) the identity by state methodology. We evaluate the effect of different population sizes, different number of genomic markers, and several heritability values for a quantitative trait. The performance of the different measures of coancestry in BLUP is evaluated in the true breeding values after truncation selection and also in terms of coancestry and diversity maintained. Accordingly, cross-performances were also carried out, that is, how prediction based on genealogical records impacts the three other measures of coancestry and inbreeding, and viceversa. Our results show that the genetic gains are very similar for all four coancestries, but the genomic-based methods are superior to using genealogical coancestries in terms of maintaining diversity measured as observed heterozygosity. Furthermore, the measure of coancestry based on shared segments of the genome seems to provide slightly better results on some scenarios, and the increase in inbreeding and loss in diversity is only slightly larger than the other genomic selection methods in those scenarios. Our results shed light on genomic selection vs. traditional genealogical-based BLUP and make the case to manage the population variability using genomic information to preserve the future success of selection programmes.

Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data

BACKGROUND

The current availability of genotypes for very large numbers of single nucleotide polymorphisms (SNPs) is leading to more accurate estimates of inbreeding coefficients and more detailed approaches for detecting inbreeding depression. In the present study, genome-wide information was used to detect inbreeding depression for two reproductive traits (total number of piglets born and number of piglets born alive) in an ancient strain of Iberian pigs (the Guadyerbas strain) that is currently under serious danger of extinction.

METHODS

A total of 109 sows with phenotypic records were genotyped with the PorcineSNP60 BeadChip v1. Inbreeding depression was estimated using a bivariate animal model in which the inbreeding coefficient was included as a covariate. We used two different measures of genomic inbreeding to perform the analyses: inbreeding estimated on a SNP-by-SNP basis and inbreeding estimated from runs of homozygosity. We also performed the analyses using pedigree-based inbreeding.

RESULTS

Significant inbreeding depression was detected for both traits using all three measures of inbreeding. Genome-wide information allowed us to identify one region on chromosome 13 associated with inbreeding depression. This region spans from 27 to 54 Mb and overlaps with a previously detected quantitative trait locus and includes the inter-alpha-trypsin inhibitor gene cluster that is involved with embryo implantation.

CONCLUSIONS

Our results highlight the value of high-density SNP genotyping for providing new insights on where genes causing inbreeding depression are located in the genome. Genomic measures of inbreeding obtained on a SNP-by-SNP basis or those based on the presence/absence of runs of homozygosity represent a suitable alternative to pedigree-based measures to detect inbreeding depression, and a useful tool for mapping studies. To our knowledge, this is the first study in domesticated animals using the SNP-by-SNP inbreeding coefficient to map specific regions within chromosomes associated with inbreeding depression.