Temporal changes in genetic variation of north European cattle breeds

Genetic diversity and recent ancestry based on whole-genome sequencing of endangered Swedish cattle breeds

Several indigenous cattle breeds in Sweden are endangered. Conservation of their genetic diversity and genomic characterization is a priority.Whole-genome sequences (WGS) with a mean coverage of 25X, ranging from 14 to 41X were obtained for 30 individuals of the breeds Fjällko, Fjällnära, Bohuskulla, Rödkulla, Ringamåla, and Väneko. WGS-based genotyping revealed 22,548,028 variants in total, comprising 18,876,115 single nucleotide polymorphisms (SNPs) and 3,671,913 indels. Out of these, 1,154,779 SNPs and 304,467 indels were novel. Population stratification based on roughly 19 million SNPs showed two major groups of the breeds that correspond to northern and southern breeds. Overall, a higher genetic diversity was observed in the southern breeds compared to the northern breeds. While the population stratification was consistent with previous genome-wide SNP array-based analyses, the genealogy of the individuals inferred from WGS based estimates turned out to be more complex than expected from previous SNP-array based estimates. Polymorphisms and their predicted phenotypic consequences were associated with differences in the coat color phenotypes between the northern and southern breeds. Notably, these high-consequence polymorphisms were not represented in SNP arrays, which are used routinely for genotyping of cattle breeds.This study is the first WGS-based population genetic analysis of Swedish native cattle breeds. The genetic diversity of native breeds was found to be high. High-consequence polymorphisms were linked with desirable phenotypes using whole-genome genotyping, which highlights the pressing need for intensifying WGS-based characterization of the native breeds.

Genetic diversity and population structure in divergent German cattle selection lines on the basis of milk protein polymorphisms

The aim of this study was to analyze the genetic structure of the casein cluster in eight selection lines of the Holstein Friesian (HF), German Simmental (GS) and German Black Pied cattle (“Deutsches Schwarzbuntes Niederungsrind”, DSN) breeds. A total of 2962 milk samples were typed at αs1-casein (αs1-CN), β-casein (β-CN), αs2-casein (αs2-CN) and κ-casein (κ-CN) loci using isoelectric focusing. The number of alleles per locus ranged from one (αs2-CN) to five (β-CN), and the average expected heterozygosity and polymorphic information content of all loci were 0.33 and 0.27, respectively. The unrooted dendrogram revealed that the selection lines of the endangered DSN breed were clearly separated from the HF and GS breeds due to their predominance of the β-CN A1 allele and the comprehensive haplotype BA1A (in the abbreviation of αs1-β-κ-CN). Temporal changes in allele distributions indicated decreasing genetic diversity at the casein loci, explaining the moderate level of genetic differentiation among selection lines (7.1 %). The variability of the casein should be exploited in future using breeding programs to select genetic lines for specific protein production in bovine milk but also to preserve biodiversity.

Utilization of farm animal genetic resources in a changing agro-ecological environment in the Nordic countries

Livestock production is the most important component of northern European agriculture and contributes to and will be affected by climate change. Nevertheless, the role of farm animal genetic resources in the adaptation to new agro-ecological conditions and mitigation of animal production’s effects on climate change has been inadequately discussed despite there being several important associations between animal genetic resources and climate change issues. The sustainability of animal production systems and future food security require access to a wide diversity of animal genetic resources. There are several genetic questions that should be considered in strategies promoting adaptation to climate change and mitigation of environmental effects of livestock production. For example, it may become important to choose among breeds and even among farm animal species according to their suitability to a future with altered production systems. Some animals with useful phenotypes and genotypes may be more useful than others in the changing environment. Robust animal breeds with the potential to adapt to new agro-ecological conditions and tolerate new diseases will be needed. The key issue in mitigation of harmful greenhouse gas effects induced by livestock production is the reduction of methane (CH₄) emissions from ruminants. There are differences in CH₄ emissions among breeds and among individual animals within breeds that suggest a potential for improvement in the trait through genetic selection. Characterization of breeds and individuals with modern genomic tools should be applied to identify breeds that have genetically adapted to marginal conditions and to get critical information for breeding and conservation programs for farm animal genetic resources. We conclude that phenotyping and genomic technologies and adoption of new breeding approaches, such as genomic selection introgression, will promote breeding for useful characters in livestock species.

New World cattle show ancestry from multiple independent domestication events

Previous archeological and genetic research has shown that modern cattle breeds are descended from multiple independent domestication events of the wild aurochs (Bos primigenius) ∼10,000 y ago. Two primary areas of domestication in the Middle East/Europe and the Indian subcontinent resulted in taurine and indicine lines of cattle, respectively. American descendants of cattle brought by European explorers to the New World beginning in 1493 generally have been considered to belong to the taurine lineage. Our analyses of 47,506 single nucleotide polymorphisms show that these New World cattle breeds, as well as many related breeds of cattle in southern Europe, actually exhibit ancestry from both the taurine and indicine lineages. In this study, we show that, although European cattle are largely descended from the taurine lineage, gene flow from African cattle (partially of indicine origin) contributed substantial genomic components to both southern European cattle breeds and their New World descendants. New World cattle breeds, such as Texas Longhorns, provide an opportunity to study global population structure and domestication in cattle. Following their introduction into the Americas in the late 1400s, semiferal herds of cattle underwent between 80 and 200 generations of predominantly natural selection, as opposed to the human-mediated artificial selection of Old World breeding programs. Our analyses of global cattle breed population history show that the hybrid ancestry of New World breeds contributed genetic variation that likely facilitated the adaptation of these breeds to a novel environment.

A review of some characteristics, socio-economic aspects and utilization of Zulu sheep: implications for conservation

Zulu sheep are Nguni sheep of Zululand and are adapted to the harsh conditions of KwaZulu-Natal. They are used by rural farmers for economic purposes. Their numbers are declining, indicating a potential extinction threat. Knowledge of their phenotypic and genotypic characteristics is essential for conservation planning. In this review, there is a focus on the utilization, socio-economic aspects, phenotypic and genotypic characteristics as well as a proposed breeding programme. A survey has shown that rural farmers in the areas of northern KwaZulu-Natal prefer to keep this breed for its adaptability, resistance to diseases and meat quality. Zulu sheep are small-framed multi-coloured animals. Mature males weigh up to 38 kg and females up to 32 kg. Based on four morphological traits and live weight, phenotypic diversity between three populations was estimated at 48%. A genetic diversity between these three populations was estimated at 22%. Live weight of Zulu sheep can be estimated using the heart girth and wither height measurements. Scrotum circumference of young rams (up to 22 months old) is reliable for estimating the live weight. Animals that were characterized in the studies were grazed extensively and no supplements were provided. There is therefore a potential of weight increase if these animals are reared in a semi-extensive environment. An open nucleus breeding scheme is thus recommended for a sustainable use and conservation of this breed. For more conclusive results, larger numbers of phenotypic and genetic characteristics, in larger numbers of Zulu sheep populations, should be investigated.

The phylogeographic system survey of native sheep breeds in the eastern and southern Central Asia

The genetic diversity and phylogenetic survey of native sheep breeds in the eastern and southern Central Asia were assessed in the present study. The clustering, principal components, structure and F statistics all demonstrate that the native sheep breeds in these regions be classified into two genetic groups: Mongolia-Tibetan sheep group and South-Southeast Asia sheep group. The Mongolia sheep group and the Tibetan sheep group had a certain degree of gene communication from the ancient times. In the present study we demonstrated that the Chinese native sheep populations belonged to Mongolia-Tibetan sheep group. However, the relationships among the sheep populations in Mongolia sheep group in China were not closely related to the geographical distance among sheep populations.

Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia

Background

Identification of global livestock diversity hotspots and their importance in diversity maintenance is essential for making global conservation efforts. We screened 52 sheep breeds from the Eurasian subcontinent with 20 microsatellite markers. By estimating and weighting differently within- and between-breed genetic variation our aims were to identify genetic diversity hotspots and prioritize the importance of each breed for conservation, respectively. In addition we estimated how important within-species diversity hotspots are in livestock conservation.

Results

Bayesian clustering analysis revealed three genetic clusters, termed Nordic, Composite and Fat-tailed. Southern breeds from close to the region of sheep domestication were more variable, but less genetically differentiated compared with more northern populations. Decreasing weight for within-breed diversity component led to very high representation of genetic clusters or regions containing more diverged breeds, but did not increase phenotypic diversity among the high ranked breeds. Sampling populations throughout 14 regional groups was suggested for maximized total genetic diversity.

Conclusions

During initial steps of establishing a livestock conservation program populations from the diversity hot-spot area are the most important ones, but for the full design our results suggested that approximately equal population presentation across environments should be considered. Even in this case, higher per population emphasis in areas of high diversity is appropriate. The analysis was based on neutral data, but we have no reason to think the general trend is limited to this type of data. However, a comprehensive valuation of populations should balance production systems, phenotypic traits and available genetic information, and include consideration of probability of success.

A genome map of divergent artificial selection between Bos taurus dairy cattle and Bos taurus beef cattle

A number of cattle breeds have become highly specialized for milk or beef production, following strong artificial selection for these traits. In this paper, we compare allele frequencies from 9323 single nucleotide polymorphism (SNP) markers genotyped in dairy and beef cattle breeds averaged in sliding windows across the genome, with the aim of identifying divergently selected regions of the genome between the production types. The value of the method for identifying selection signatures was validated by four sources of evidence. First, differences in allele frequencies between dairy and beef cattle at individual SNPs were correlated with the effects of those SNPs on production traits. Secondly, large differences in allele frequencies generally occurred in the same location for two independent data sets (correlation 0.45) between sliding window averages. Thirdly, the largest differences in sliding window average difference in allele frequencies were found on chromosome 20 in the region of the growth hormone receptor gene, which carries a mutation known to have an effect on milk production traits in a number of dairy populations. Finally, for the chromosome tested, the location of selection signatures between dairy and beef cattle was correlated with the location of selection signatures within dairy cattle.

Changes in inbreeding of U.S. Herefords during the twentieth century

Genetic diversity in the U.S. Hereford population was characterized by examining the level and rate of inbreeding and effective population size. Pedigree records for 20,624,418 animals were obtained from the American Hereford Association, of which 96.1% had both parents identified. Inbreeding coefficients were computed and mean inbreeding (Fx) calculated by year from 1900 to 2001. Inbreeding increased rapidly between 1900 and 1945. From 1946, inbreeding increased linearly to a maximum of 11.5% in 1966. Throughout the 1970s and 1980s, mean inbreeding decreased to mid-century levels. Several alternatives were investigated to explain this decline. The average relationship between prominent sires fell from 20 to 12% during the time that the level of inbreeding decreased, which reflects an increase in the popularity of certain less fashionable sire lines that would have temporarily decreased inbreeding. Pedigrees were constructed for animals born after 1990. This subsample of animals with no missing ancestors in at least 12 generations did not exhibit a decrease in inbreeding. Missing ancestral information therefore contributed to the apparent decline. One cause of missing ancestry results from outcrossing to imported animals. The effect of missing ancestry was investigated by simulating the missing ancestors. In 2001, Fx was 9.8%, and approximately 95% of individuals were inbred. The maximal inbreeding coefficient was 76%. The annual change in mean inbreeding (DeltaFx) was estimated for Herefords born during five time periods from 1946 to 2001, where inbreeding was changing at different linear rates. The DeltaFx for the most recent generation (1990 to 2001) was 0.12%/yr. Assuming a generation interval of 4.88 yr, the estimated effective population size was 85. This study provides a benchmark of current genetic diversity in the Hereford population. Results indicate that inbreeding is accumulating linearly and below critical levels. Increases in the adoption of reproductive technologies could decrease genetic diversity, and in the future, we may need to consider strategies to minimize inbreeding.

Prioritization for conservation of northern European cattle breeds based on analysis of microsatellite data

Northern European indigenous cattle breeds are currently endangered and at a risk of becoming extinct. We analyzed variation at 20 microsatellite loci in 23 indigenous, 3 old imported, and 9 modern commercial cattle breeds that are presently distributed in northern Europe. We measured the breeds' allelic richness and heterozygosity, and studied their genetic relationships with a neighbor-joining tree based on the Chord genetic distance matrix. We used the Weitzman approach and the core set diversity measure of Eding et al. (2002) to quantify the contribution of each breed to the maximum amount of genetic diversity and to identify breeds important for the conservation of genetic diversity. We defined 11 breeds as a "safe set" of breeds (not endangered) and estimated a reduction in genetic diversity if all nonsafe (endangered) breeds were lost. We then calculated the increase in genetic diversity by adding one by one each of the nonsafe breeds to the safe set (the safe-set-plus-one approach). The neighbor-joining tree grouped the northern European cattle breeds into Black-and-White type, Baltic Red, and Nordic cattle groups. Väne cattle, Bohus Poll, and Danish Jersey had the highest relative contribution to the maximum amount of genetic diversity when the diversity was quantified by the Weitzman diversity measure. These breeds not only showed phylogenetic distinctiveness but also low within-population variation. When the Eding et al. method was applied, Eastern Finncattle and Lithuanian White Backed cattle contributed most of the genetic variation. If the loss of the nonsafe set of breeds happens, the reduction in genetic diversity would be substantial (72%) based on the Weitzman approach, but relatively small (1.81%) based on the Eding et al. method. The safe set contained only 66% of the observed microsatellite alleles. The safe-set-plus-one approach indicated that Bohus Poll and Väne cattle contributed most to the Weitzman diversity, whereas the Eastern Finncattle contribution was the highest according to the Eding et al. method. Our results indicate that both methods of Weitzman and Eding et al. recognize the importance of local populations as a valuable resource of genetic variation.

Conservation value of non-native banteng in northern Australia

The global species extinction crisis has provided the impetus for elaborate translocation, captive breeding, and cloning programs, but more extreme actions may be necessary. We used mitochondrial DNA, Y-chromosome, and nuclear lactoferrin-encoding gene sequencing to identify a wild population of a pure-strain endangered bovid (Bos javanicus) introduced into northern Australia over 150 years ago. This places the Australian population in a different conservation category relative to its domesticated conspecific in Indonesia (i.e., Bali cattle) that has varying degrees of introgression from other domesticated Bos spp. The success of this endangered non-native species demonstrates that although risky, the deliberate introduction of threatened exotic species into non-native habitat may provide, under some circumstances, a biologically feasible option for conserving large herbivores otherwise imperiled in their native range.

Maximum-likelihood and markov chain monte carlo approaches to estimate inbreeding and effective size from allele frequency changes

Maximum-likelihood and Bayesian (MCMC algorithm) estimates of the increase of the Wright-Malécot inbreeding coefficient, F(t), between two temporally spaced samples, were developed from the Dirichlet approximation of allelic frequency distribution (model MD) and from the admixture of the Dirichlet approximation and the probabilities of fixation and loss of alleles (model MDL). Their accuracy was tested using computer simulations in which F(t) = 10% or less. The maximum-likelihood method based on the model MDL was found to be the best estimate of F(t) provided that initial frequencies are known exactly. When founder frequencies are estimated from a limited set of founder animals, only the estimates based on the model MD can be used for the moment. In this case no method was found to be the best in all situations investigated. The likelihood and Bayesian approaches give better results than the classical F-statistics when markers exhibiting a low polymorphism (such as the SNP markers) are used. Concerning the estimations of the effective population size all the new estimates presented here were found to be better than the F-statistics classically used.

Estimating effective population size and migration rates from genetic samples over space and time

In the past, moment and likelihood methods have been developed to estimate the effective population size (N(e)) on the basis of the observed changes of marker allele frequencies over time, and these have been applied to a large variety of species and populations. Such methods invariably make the critical assumption of a single isolated population receiving no immigrants over the study interval. For most populations in the real world, however, migration is not negligible and can substantially bias estimates of N(e) if it is not accounted for. Here we extend previous moment and maximum-likelihood methods to allow the joint estimation of N(e) and migration rate (m) using genetic samples over space and time. It is shown that, compared to genetic drift acting alone, migration results in changes in allele frequency that are greater in the short term and smaller in the long term, leading to under- and overestimation of N(e), respectively, if it is ignored. Extensive simulations are run to evaluate the newly developed moment and likelihood methods, which yield generally satisfactory estimates of both N(e) and m for populations with widely different effective sizes and migration rates and patterns, given a reasonably large sample size and number of markers.

Genetic distances within and across cattle breeds as indicated by biallelic AFLP markers

We tested the use of biallelic Amplified fragment length polymorphism (AFLP) polymorphisms for the estimation of relative genetic distances of cattle individuals within or across breeds. An allele permutation procedure was developed to estimate the stochastic variation of the genetic distance that is inherent to a given dataset. In a panel of 47 Holstein-Friesian cattle analysed with 248 polymorphic markers, the average genetic distance of bulls selected for breeding was slightly lower than the distance of the cows. The observed standard deviation (SD) of the distance indicated genetic subdivision, which for the bulls was explained by variation in the additive relationship derived from herdbook data. Animals from three different breeds, the highly selected Holstein-Friesian, the Italian Brown and the historic Maremmana, were compared on the basis of 106 polymorphic markers. No breed-specific fragments were observed. The mean pair-wise genetic distance within breeds was 85% of the value across breeds, but principal coordinates analysis clustered the animals according to their breed of origin. Calculation of distances between the breeds indicated a relatively divergent position of the Maremmana, relative to the two other breeds. However, biallelic markers indicate that the process of breed formation had only a limited effect on the diversity at marker loci.

Likelihood-based estimation of the effective population size using temporal changes in allele frequencies: a genealogical approach

A new genetic estimator of the effective population size (N(e)) is introduced. This likelihood-based (LB) estimator uses two temporally spaced genetic samples of individuals from a population. We compared its performance to that of the classical F-statistic-based N(e) estimator (N(eFk)) by using data from simulated populations with known N(e) and real populations. The new likelihood-based estimator (N(eLB)) showed narrower credible intervals and greater accuracy than (N(eFk)) when genetic drift was strong, but performed only slightly better when genetic drift was relatively weak. When drift was strong (e.g., N(e) = 20 for five generations), as few as approximately 10 loci (heterozygosity of 0.6; samples of 30 individuals) are sufficient to consistently achieve credible intervals with an upper limit <50 using the LB method. In contrast, approximately 20 loci are required for the same precision when using the classical F-statistic approach. The N(eLB) estimator is much improved over the classical method when there are many rare alleles. It will be especially useful in conservation biology because it less often overestimates N(e) than does N(eLB) and thus is less likely to erroneously suggest that a population is large and has a low extinction risk.