Are cattle, sheep, and goats endangered species?

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1-5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Insights into the genomic homogeneity of Moroccan indigenous sheep breeds though the lens of runs of homozygosity

Numerous studies have indicated that Morocco's indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D'man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24-48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1-6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds' genomic landscapes and adaptive mechanisms.

Circulation of small ruminant lentivirus in endangered goat and sheep breeds of Southern Italy

According to the Domestic Animal Diversity Information System (DAD-IS) of the FAO, Italy has one of the largest numbers of local small ruminant breeds among European countries. In Southern Italy, namely the Campania Region, Bagnolese and Laticauda sheep breeds and Cilentana goat breeds are considered endangered according to the DAD-IS. Conservation of endangered animal breeds is a goal of the European Union (EU). However, the role of infectious diseases as risk factors for endangered breeds has rarely been considered. Small ruminant lentiviruses (SRLV) infect sheep and goats, causing slow-progressive, persistent, and debilitating diseases that can lead to animal death and productivity loss. In this study, we investigated the presence of SRLV in Bagnolese, Laticauda, and Cilentana breeds using a commercial ELISA in parallel with an in-house ELISA. The results of the two tests were in good agreement (Cohen Kappa 0.84, 95 % CI = 0.76-0.93). Discrepancies between the two tests were resolved using western blotting. In total, 430 samples were tested (248 Bagnolese, 125 Laticauda, and 57 Cilentana). The apparent prevalence rates were 12.5 %, 6.4 %, and 1.7 % in Bagnolese, Laticauda, and Cilentana, respectively. In the molecular analysis of 11 proviral partial sequences, subtypes B2 and A24 were identified in two Bagnolese herds. Owing to the beneficial role of sheep and goat breeding in marginal areas, it is important to screen the entire population and implement control/eradication of SRLV infections in conjunction with each conservation program.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs

BACKGROUND

Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear.

RESULTS

Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments.

CONCLUSIONS

In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

The effects of repeated freezing and thawing on bovine sperm morphometry and function

Refreezing the remaining genetic resources after in vitro fertilization (IVF) can conserve genetic materials. However, the precise damage inflicted by repeated freezing and thawing on bovine sperm and its underlying mechanism remain largely unexplored. Thus, this study investigates the impact of repeated freeze-thaw cycles on sperm. Our findings indicate that such cycles significantly reduce sperm viability and motility. Furthermore, the integrity of the sperm plasma membrane and acrosome is compromised during this process, exacerbating the advanced apoptosis triggered by oxidative stress. Additionally, transmission electron microscopy exposed severe damage to the plasma membranes of both the sperm head and tail. Notably, the "9 + 2" structure of the tail was disrupted, along with a significant decrease in the level of the axonemal protein DNAH10, leading to reduced sperm motility. IVF outcomes revealed that repeated freeze-thaw cycles considerably impair sperm fertilization capability, ultimately reducing the blastocyst rate. In summary, our research demonstrates that repeated freeze-thaw cycles lead to a decline in sperm viability and motility, attributed to oxidative stress-induced apoptosis and DNAH10-related dynamic deficiency. As a result, the utility of semen is compromised after repeated freezing.

Detection of selective sweep in European wild sheep breeds

In wild animal populations, there is a differentiation between populations due to natural selection. The direction and pressure of natural selection in the wild sheep are different in the various geographic areas. Linkage disequilibrium studies showed that regions of the genome in whole wild sheep are under natural selection and that natural selection can affect immune or reproductive or metabolic traits. The study aimed to identify genomic regions under natural selection in wild sheep. For this purpose, the genetic information of 24 European wild sheep and 24 Sardinian wild sheep was used. The genotypes were determined using Illumina 50 K SNPChip arrays based on Oar_4.0 version of the sheep genome. After quality control steps, finally, 31,560 SNP markers were analyzed. The value of LD was calculated by calculating the r2 statistic between all pairs of locations through PLINK software. To identify signs of selection based on linkage disequilibrium methods, an extended haplotype homozygosity test of XP-EHH crossing population and iHS intrapopulation was used. The results of iHS studies showed that in European and Sardinian wild sheep, the highest iHS coefficient under natural selection was observed on 3 and 2 chromosome numbers, respectively. Also, the results of XP-EHH studies showed that the largest XP-EHH coefficients under natural selection in European wild sheep compared to Sardinian and vice versa in Sardinian wild sheep compared to European wild sheep were observed on 3 and 16 chromosome numbers, respectively. In addition, the results of gene cycle studies showed that COPB1, SEC24D, ZDHHC17, BBS4, RFX3, SLC26A8, CAMK2D, GRIA1, GRM1, GRID2, PPP2R1A, CPEB4, PLEKHA5 and KIF13A, VPS39, VPS53, DTNBP1, DYNC1I1, FAM91A genes are under natural selection in Sardinian and European wild sheeps, respectively. The direction and selection pressure of natural selection in the two breeds of wild sheep is different due to different geographic conditions.

Late History of Cattle Breeds in Central Europe in Light of Genetic and Archaeogenetic Sources-Overview, Thoughts, and Perspectives

Although Europe was not a primary centre of cattle domestication, its expansion from the Middle East and subsequent development created a complex pattern of cattle breed diversity. Many isolated populations of local historical breeds still carry the message about the physical and genetic traits of ancient populations. Since the way of life of human communities starting from the eleventh millennium BP was strongly determined by livestock husbandry, the knowledge of cattle diversity through the ages is helpful in the interpretation of many archaeological findings. Historical cattle diversity is currently at the intersection of two leading directions of genetic research. Firstly, it is archaeogenetics attempting to recover and interpret the preserved genetic information directly from archaeological finds. The advanced archaeogenetic approaches meet with the population genomics of extant cattle populations. The immense amount of genetic information collected from living cattle, due to its key economic role, allows for reconstructing the genetic profiles of the ancient populations backwards. The present paper aims to place selected archaeogenetic, genetic, and genomic findings in the picture of cattle history in Central Europe, as suggested by archaeozoological and historical records. Perspectives of the methodical connection between the genetic approaches and the approaches of traditional archaeozoology, such as osteomorphology and osteometry, are discussed. The importance, actuality, and effectiveness of combining different approaches to each archaeological find, such as morphological characterization, interpretation of the historical context, and molecular data, are stressed.

The intestinal microflora diversity of aboriginal chickens in Jiangxi province, China

Intestinal microbiota can coevolve with host to form symbiotic relationship and be participated in the regulation of host physiological function. At present, there is no clear explanation on the effect of intestinal microflora in Jiangxi aboriginal chickens. Here, we investigated the association between gut microbiota and host genome of Jiangxi local chickens using 16S rRNA sequencing and genome-wide association studies (GWAS). The results showed that the breeds and genders had important effects on the intestinal microbiota of chickens. A total of 28 SNPs in 14 regions of the chicken genome were related to the relative abundance of microorganisms in 5 genera: Clostridium_sensu_stricto_1, Enterococcus, Gallibacterium, Turicibacter, and Rikenellaceae_RC9_gut_group. A total of 17 candidate genes were identified composition of chicken microbiome and show an association between the host genome and the chicken intestinal microbiota, which also unveiled the diversity of intestinal microbes in Jiangxi chickens. Given the correlation between chicken genome and intestinal microbe found in the present study, a new idea for the protection of aboriginal chicken genetic resources in China could be provided.

Phenotype transition from wild mouflon to domestic sheep

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

IBD sharing patterns as intra-breed admixture indicators in small ruminants

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Parasite exchange and hybridisation at a wild-feral-domestic interface

Interactions between wild, feral, and domestic animals are of economic and conservation significance. The pigeon Columba livia is a synanthropic species in a feral form, but it also includes the rare Rock Dove. Columba livia is an important player at the wild-domestic interface, acting as a carrier of avian diseases, and the feral form threatens Rock Doves with extinction via hybridisation. Despite its abundance, little is known about drivers of disease prevalence in C. livia, or how disease and hybridisation represent synergistic threats to Rock Doves. We focused on infection by the parasite Trichomonas, first collating prevalence estimates in domestic and free-living populations from relevant studies of C. livia. Second, we characterised variation in the diversity and prevalence of Trichomonas among three C. livia populations in the United Kingdom: a feral, a Rock Dove, and a feral-wild hybrid population. Across multiple continents, free-living pigeons had lower Trichomonas infection than captive conspecifics, but the effect was weak. Environmental factors which could impact Trichomonas infection status did not explain variation in infection among populations. Among the British populations, strain diversity varied, and there was lower parasite prevalence in Rock Doves than feral pigeons. Individual infection status was not explained by the available covariates, including hybrid score and site. The drivers of Trichomonas prevalence are unclear, perhaps due to idiosyncratic local-scale drivers. However, given the population-level variation in both infection prevalence and introgressive hybridisation, the potential combined effects could accelerate the extinction of the Rock Dove. Further study of the synergistic effects of multiple types of biotic interactions at the wild-feral-domestic interface is warranted, especially where vagile, globally distributed and superabundant animals are involved.

In vitro and in vivo embryo production efficiency in Flemish and Holstein donor females

The aim of this study was to compare embryo production efficiency in Flemish and Holstein donor females using ovum pick-up and in vitro fertilization (OPU-IVF) or in vivo production (superovulation; SOV) procedures. The study was conducted using a split-plot design, with eight Flemish and eight Holstein non-lactating cycling females. Females were subjected to ten weekly OPU/IVF sessions and/or two SOV/embryo collections sessions at a 63-day interval, for a total of 160 OPU-IVF and 32 SOV sessions. Mean numbers of follicles and corpora lutea, and cumulus-oocyte complex (COC) recovery rates were similar between breeds after the OPU and SOV sessions. However, Flemish donors yielded better quality grade II COCs (301, 41.9%) than Holstein females (609, and 202, 33.1%). Also, cleavage and blastocyst rates, and the total number and the mean number of viable embryos obtained after OPU-IVF were higher in Flemish (49.6% and 11.8%, and 63 and 11.8 per donor, respectively) than in Holstein (32.8% and 7.2%, and 34 and 7.2 per donor, respectively) females. Flemish females were also more efficient in yielding viable embryos after SOV (111, 7.3 per donor) than Holstein (48, 3.3 per donor) females. Overall, Flemish donor females had better responses to OPU-IVF or SOV procedures than Holstein counterparts. Irrespective of the breeds, SOV procedures were more efficient than OPU-IVF in yielding more viable embryos, under the conditions of this study. Both reproductive procedures were useful tools for the genetic conservation of the Flemish cattle breed in Southern Brazil.

ZooMS confirms geometric morphometrics species identification of ancient sheep and goat

Geometric morphometrics can effectively distinguish isolated third lower molars of present-day sheep and goat, but its applicability to archaeological specimens has yet to be established. Using a modern reference collection of 743 sheep and goats and a two-dimensional landmark-based geometric morphometric (GMM) protocol, this study aimed to morphometrically identify 109 archaeological specimens, used as case studies, dating from the Late Neolithic to the modern period/era. These morphometric identifications were then compared to molecular identifications via collagen peptide mass fingerprinting, known as Zooarcheology by Mass Spectrometry (ZooMS). ZooMS confirmed the morphometric identifications for 104 specimens, with the five misidentified specimens all morphometrically identified as goat. Modern sheep and goats have larger teeth and distinct shapes compared to their archaeological counterparts, suggesting strong differences between archaeological and modern specimens potentially linked with recent breed improvement or geographical origin of the specimens. In addition, for both species, some of the archaeological dental morphologies do not match with any of our modern references. This study validates the applicability of geometric morphometrics for identifying isolated archaeological sheep and goat teeth. It represents a stepping stone for future, non-destructive, bioarchaeological studies of the two species.

Comparative genomic analysis uncovers candidate genes related with milk production and adaptive traits in goat breeds

During the process of animal domestication, both natural and artificial selection cause variation in allele frequencies among populations. Identifying genomic areas of selection in domestic animals may aid in the detection of genomic areas linked to ecological and economic traits. We studied genomic variation in 140 worldwide goat individuals, including 75 Asian, 30 African and 35 European goats. We further carried out comparative population genomics to detect genomic regions under selection for adaptability to harsh conditions in local Asian ecotypes and also milk production traits in European commercial breeds. In addition, we estimated the genetic distances among 140 goat individuals. The results showed that among all studied goat groups, local breeds from West and South Asia emerged as an independent group. Our search for selection signatures in local goats from West and South Asia revealed candidate genes related to adaptation to hot climate (HSPB6, HSF4, VPS13A and NBEA genes) and immune response (IL7, IL5, IL23A and LRFN5) traits. Furthermore, selection signatures in European commercial goats involved several milk production related genes, such as VPS13C, NCAM2, TMPRSS15, CSN3 and ABCG2. The identified candidate genes could be the fundamental genetic resource for enhancement of goat production and environmental-adaptive traits, and as such they should be used in goat breeding programs to select more efficient breeds.

Characterization of the genetic pool of the Canadienne dairy cattle breed

BACKGROUND

Canadienne cattle are the oldest breed of dairy cattle in North America. The Canadienne breed originates from cattle that were brought to America by the mid-seventeenth century French settlers. The herd book was established in 1886 and the current breed characteristics include dark coat color, small size compared to the modern Holstein breed, and overall rusticity shaped by the harsh environmental conditions that were prevalent during the settlement of North America. The Canadienne breed is an invaluable genetic resource due to its high resilience, longevity and fertility. However, it is heavily threatened with a current herd limited to an estimated 1200 registered animals, of which less than 300 are fullblood. To date, no effort has been made to document the genetic pool of this heritage breed in order to preserve it.

RESULTS

In this project, we used genomic data, which allow a precise description of the genetic makeup of a population, to provide valuable information on the genetic diversity of this heritage breed and suggest management options for its long-term viability. Using a panel that includes 640,000 single nucleotide polymorphisms (SNPs), we genotyped 190 animals grouped into six purity ranges. Unsupervised clustering analyses revealed three genetically distinct groups among those with the higher levels of purity. The observed heterozygosity was higher than expected even in the 100% purebreds. Comparison with Holstein genotypes showed significantly shorter runs of homozygosity for the Canadienne breed, which was unexpected due to the high inbreeding value calculated from pedigree data.

CONCLUSIONS

Overall, our data indicate that the fullblood gene pool of the Canadienne breed is more diversified than expected and that bloodline management could promote breed sustainability. In its current state, the Canadienne is not a dead-end breed but remains highly vulnerable due to its small population size.

Genetic diversity evaluation of Sojat goat population of India

Sojat goat forms the part and parcel of the lives of the farmers in the Barmer region of Rajasthan. These are large-sized goats, white in colour, and have dual utility. These goats fetch premium price during the Eid festival. Assessment of diversity is essential for germplasm characterization and management. Genomic microsatellite markers being a valuable tool for estimating genetic diversity were selected for exploring existing genetic variability in the Sojat goat population. The standard metrics of genomic diversity detected moderate variability with a total of 162 alleles across 22 loci in this lesser-known population. The expected number of alleles had a mean value of 3.40±0.39. Similarly, a moderate magnitude of diversity was recorded in the Sojat population as the mean observed heterozygosity was 0.54±0.05. Expected heterozygosity was higher than the observed (0.60±0.06), indicating a deviation from Hardy-Weinberg Equilibrium (HWE) and the possibility of inbreeding due to the non-random mating in the population. Accordingly, significant heterozygote deficiency was noticed (F=0.08±0.03). The population did not suffer a reduction in effective population size in the last few generations. Mutation drift equilibrium did not reveal significant heterozygosity excess under different models of microsatellite evolution and no shift was recorded in the frequency distribution of alleles. To conclude, the results provided the first insights into the genetic diversity of Sojat goats. A moderate genetic variability with heterozygote deficiency within the population warrants immediate attention for scientific management of this unique goat population to conserve the existing genetic variation and to avoid any escalation of inbreeding.

Legacies of domestication, Neolithic diffusion and trade between Indian subcontinent and Island Southeast Asia shape maternal genetic diversity of Andaman cattle

Andaman cattle is a precious indigenous livestock species endemic to Andaman and Nicobar Islands, India. Till date, origin and genetic makeup of the breed which is warranted for breed conservation is not known. Moreover, the spread of zebu cattle from Indus valley to different parts of Island Southeast Asia (ISEA) is not properly understood. Here, we report the genetic diversity, population structure of Andaman cattle and their evolution in the context of epicentre of zebu domestication and ISEA. High genetic diversity in complete mitochondrial D-loop sequences indicated the ability of the breed to withstand impending climate change. Total 81 haplotypes were detected and all of them except three belonged to Bos indicus. The presence of taurine haplotypes in Andaman cattle indicate introgression by European-derived cattle. A poor phylogenetic signal of Andaman cattle with genetic affinities with cattle of Indian subcontinent and ISEA was observed. The poor phylogenetic structure may be due to multidirectional gene flow from Indian subcontinent and ISEA, with which Andaman shares a close cultural and trade relationship from Neolithic age. We hypothesize that Andaman cattle is the outcome of Neolithic diffusion from centre of zebu domestication along with multidirectional commercial exchange between Indian subcontinent and ISEA.

Animal board invited review: Genomic-based improvement of cattle in response to climate change

Climate change brings challenges to cattle production, such as the need to adapt to new climates and pressure to reduce greenhouse emissions (GHG). In general, the improvement of traits in current breeding goals is favourably correlated with the reduction of GHG. Current breeding goals and tools for increasing cattle production efficiency have reduced GHG. The same amount of production can be achieved by a much smaller number of animals. Genomic selection (GS) may offer a cost-effective way of using an efficient breeding approach, even in low- and middle-income countries. As climate change increases the intensity of heatwaves, adaptation to heat stress leads to lower efficiency of production and, thus, is unfavourable to the goal of reducing GHG. Furthermore, there is evidence that heat stress during cow pregnancy can have many generation-long lowering effects on milk production. Both adaptation and reduction of GHG are among the difficult-to-measure traits for which GS is more efficient and suitable than the traditional non-genomic breeding evaluation approach. Nevertheless, the commonly used within-breed selection may be insufficient to meet the new challenges; thus, cross-breeding based on selecting highly efficient and highly adaptive breeds may be needed. Genomic introgression offers an efficient approach for cross-breeding that is expected to provide high genetic progress with a low rate of inbreeding. However, well-adapted breeds may have a small number of animals, which is a source of concern from a genetic biodiversity point of view. Furthermore, low animal numbers also limit the efficiency of genomic introgression. Sustainable cattle production in countries that have already intensified production is likely to emphasise better health, reproduction, feed efficiency, heat stress and other adaptation traits instead of higher production. This may require the application of innovative technologies for phenotyping and further use of new big data techniques to extract information for breeding.

Inbreeding and pedigree analysis of the European red dairy cattle

BACKGROUND

Red dairy cattle breeds have an important role in the European dairy sector because of their functional characteristics and good health. Extensive pedigree information is available for these breeds and provides a unique opportunity to examine their population structure, such as effective population size, depth of the pedigree, and effective number of founders and ancestors, and inbreeding levels. Animals with the highest genetic contributions were identified. Pedigree data included 9,073,403 animals that were born between 1900 and 2019 from Denmark, Finland, Germany, Latvia, Lithuania, the Netherlands, Norway, Poland, and Sweden, and covered 32 breeds. The numerically largest breeds were Red Dairy Cattle and Meuse-Rhine-Yssel.

RESULTS

The deepest average complete generation equivalent (9.39) was found for Red Dairy Cattle in 2017. Mean pedigree completeness ranged from 0.6 for Finncattle to 7.51 for Red Dairy Cattle. An effective population size of 166 animals was estimated for the total pedigree and ranged from 35 (Rotes Höhenvieh) to 226 (Red Dairy Cattle). Average generation intervals were between 5 and 7 years. The mean inbreeding coefficient for animals born between 1960 and 2018 was 1.5%, with the highest inbreeding coefficients observed for Traditional Angler (4.2%) and Rotes Höhenvieh (4.1%). The most influential animal was a Dutch Meuse-Rhine-Yssel bull born in 1960. The mean inbreeding level for animals born between 2016 and 2018 was 2% and highest for the Meuse-Rhine-Yssel (4.64%) and Rotes Hohenvieh breeds (3.80%).

CONCLUSIONS

We provide the first detailed analysis of the genetic diversity and inbreeding levels of the European red dairy cattle breeds. Rotes Höhenvieh and Traditional Angler have high inbreeding levels and are either close to or below the minimal recommended effective population size, thus it is necessary to implement tools to monitor the selection process in order to control inbreeding in these breeds. Red Dairy Cattle, Vorderwälder, Swedish Polled and Hinterwälder hold more genetic diversity. Regarding the Meuse-Rhine-Yssel breed, given its decreased population size, increased inbreeding and low effective population size, we recommend implementation of a breeding program to prevent further loss in its genetic diversity.

High-Density Genomic Characterization of Native Croatian Sheep Breeds

A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH>2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.

Comparison of Female Verzaschese and Camosciata delle Alpi Goats' Hematological Parameters in The Context of Adaptation to Local Environmental Conditions in Semi-Extensive Systems in Italy

Local livestock breeds are pivotal to ensure sustainable agriculture and represent a real genetic resource in the context of genetic variability reduction. This study aimed at broadening the knowledge about hematological values in female Verzaschese and Camosciata delle Alpi goats (an Italian local and a cosmopolitan goat breed, respectively) and investigating factors affecting them (i.e., breed, age, season). Blood samples were collected monthly from 34 Verzaschese and 37 Camosciata delle Alpi female goats kept under the same semi-extensive farming system for a whole year. The main hematological parameters were evaluated, and descriptive as well as inferential statistical analyses were performed. Reference intervals for complete blood cell count were established and several inter-breed differences were found. In particular, most of the red blood cell parameters were higher in Verzaschese than Camosciata delle Alpi (p < 0.05) suggesting a greater gastrointestinal parasites’ resilience of the local breed. The age effect (p < 0.05) was consistent with the literature while the season effect (p < 0.05) could be explained by the melatonin immunostimulant action and gastrointestinal parasitism influence. Overall, differences in blood values could be attributed to physiological changes and adaptive strategies developed through centuries highlighting the remarkable rusticity and adaptation to the environment and farming system of the local breed.

The wild ancestors of domestic animals as a neglected and threatened component of biodiversity

Domestic animals have immense economic, cultural, and practical value and have played pivotal roles in the development of human civilization. Many domesticates have, among their wild relatives, undomesticated forms representative of their ancestors. Resurgent interest in these ancestral forms has highlighted the unclear genetic status of many, and some are threatened with extinction by hybridization with domestic conspecifics. We considered the contemporary status of these ancestral forms relative to their scientific, practical, and ecological importance; the varied impacts of wild-domestic hybridization; and the challenges and potential resolutions involved in conservation efforts. Identifying and conserving ancestral forms, particularly with respect to disentangling patterns of gene flow from domesticates, is complex because of the lack of available genomic and phenotypic baselines. Comparative behavioral, ecological, and genetic studies of ancestral-type, feral, and domestic animals should be prioritized to establish the contemporary status of the former. Such baseline information will be fundamental in ensuring successful conservation efforts.

Capturing Genetic Diversity and Selection Signatures of the Endangered Kosovar Balusha Sheep Breed

There is a growing concern about the loss of animal genetic resources. The aim of this study was to analyze the genetic diversity and potential peculiarity of the endangered Kosovar sheep breed Balusha. For this purpose, a dataset consisting of medium-density SNP chip genotypes (39,879 SNPs) from 45 Balusha sheep was generated and compared with SNP chip genotypes from 29 individuals of a second Kosovar breed, Bardhoka. Publicly available SNP genotypes from 39 individuals of the relatively closely located sheep breeds Istrian Pramenka and Ruda were additionally included in the analyses. Analysis of heterozygosity, allelic richness and effective population size was used to assess the genetic diversity. Inbreeding was evaluated using two different methods (F_IS, F_ROH). The standardized F_ST (d_i) and cross-population extended haplotype homozygosity (XPEHH) methods were used to detect signatures of selection. We observed the lowest heterozygosity (H_O = 0.351) and effective population size (Ne₅ = 25, Ne₅₀ = 228) for the Balusha breed. The mean allelic richness levels (1.780–1.876) across all analyzed breeds were similar and also comparable with those in worldwide breeds. F_ROH estimates (0.023–0.077) were highest for the Balusha population, although evidence of decreased inbreeding was observed in F_IS results for the Balusha breed. Two Gene Ontology (GO) TERMs were strongly enriched for Balusha, and involved genes belonging to the melanogenesis and T cell receptor signaling pathways, respectively. This could result from selection for the special coat color pattern of Balusha (black head) and resistance to certain infectious diseases. The analyzed diversity parameters highlight the urgency to preserve the local Kosovar Balusha sheep as it is clearly distinguished from other sheep of Southeastern Europe, has the lowest diversity level and may harbor valuable genetic variants, e.g., for resistance to infectious diseases.

Selection Signature Analyses Revealed Genes Associated With Adaptation, Production, and Reproduction in Selected Goat Breeds in Kenya

Artificial and natural selection in livestock is expected to leave unique footprints on their genomes. Goat breeds in Kenya have evolved for survival, breeding, and production in various harsh ecological areas, and their genomes are likely to have acquired unique alleles for adaptation to such diverse production environments and other traits of economic importance. To investigate signals of selection for some selected goat breeds in Kenya, Alpine (n = 29), Galla (n = 12), Saanen (n = 24), and Toggenburg (n = 31) were considered. A total of 53,347 single-nucleotide polymorphisms (SNPs) generated using the Illumina GoatSNP50 BeadChip were analyzed. After quality control, 47,663 autosomal single-nucleotide polymorphisms remained for downstream analyses. Several complementary approaches were applied for the following analyses: integrated Haplotype Score (iHS), cross-population-extended haplotype homozygosity (XP-EHH), hapFLK, and FLK. A total of 404 top genomic regions were identified across all the four breeds, based on the four complementary analyses. Out of the 16 identified putative selection signature regions by the intersection of multiple-selective signal analyses, most of the putative regions were found to overlap significantly with the iHS and XP-EHH analyses on chromosomes 3, 4, 10, 15, 22, and 26. These regions were enriched with some genes involved in pathways associated directly or indirectly with environmental adaptation regulating immune responses (e.g., HYAL1 and HYAL3), milk production (e.g., LEPR and PDE4B), and adaptability (e.g., MST1 and PCK). The results revealed few intersect between breeds in genomic selection signature regions. In general, this did not present the typical classic selection signatures as predicted due to the complex nature of the traits. The results support that some various selection pressures (e.g., environmental challenges, artificial selection, and genome admixture challenges) have molded the genome of goat breeds in Kenya. Therefore, the research provides new knowledge on the conservation and utilization of these goat genetic resources in Kenya. In-depth research is needed to detect precise genes connected with adaptation and production in goat breeds in Kenya.

Added Value of Local Sheep Breeds in Alpine Agroecosystems

This study is part of a project (Sheep Al.L. Chain, RDP Veneto Region) aiming to improve the competitiveness of local sheep breed farms through valorization of their links with mountain agroecosystems. We considered two local sheep breeds of the eastern Italian Alps, “Alpagota” and “Lamon”, which have a population of 400 and 3000 heads, respectively, and are used to produce lambs for typical products. A total of 35 farms (17 for Alpagota, with a total of 1652 heads; 18 for Lamon, with a total of 337 heads) were surveyed to collect data on farm organization, flock structure and management (farm questionnaire), land use management (GIS approach), and value chain organization (participatory processes). The link between the two local sheep breeds and mountain agroecosystems is very strong: land use is characterized by a large number of small patches of grasslands in marginal areas. Moving from the results of this study, a set of strategies aiming at improving the competitiveness of these systems have been proposed. Communication to the consumers and to the relevant stakeholders of the added value of local sheep breeds in marginal mountain agroecosystems can contribute to favor the resilience of small ruminant farms and the conservation of Alpine sheep breeds.

Single nucleotide polymorphisms and metabolic biochemical profile of productive markers characterize three European breeds of dairy cattle

The objective of this study was to investigate polymorphisms of DGAT1, FABP, OLR1 and ATP1A1 genes using PCR-DNA sequencing, and to associate these genetic structures to changes in metabolic biochemical markers and milk composition indicators in a total of 90 dairy cows of the Holstein, Simmental, and Brown Swiss breeds (30 cows each). PCR was carried out for amplification of 411-bp of DGAT1, 525-bp of FABP, 582-bp of OLR1, and 300-bp of ATP1A1 genes. Three breeds’ nucleotide sequence variations in the form of single nucleotide polymorphisms (SNPs) were detailed by DNA sequencing analysis. Chisquare analysis showed that the distribution of all discovered SNPs varied significantly (P < 0.001). Biochemical indices in cow’s serum revealed no significant difference in serum total protein, albumin, and total cholesterol among the three breeds. However, triglyceride showed a significant increase in Simmental compared to either Holsteins or Brown Swiss, while the highest mean value of triiodothyronine (T3) and tetraiodothyronine (T4) was detected in Holstein dairy cows The milk composition indicators analysis revealed that milk protein, sugar, and density were significantly higher in Holsteins than both Simmental and Brown Swiss. Meanwhile, milk fat and total solids revealed a significantly higher increase in Simmental than both brown Swiss and Holstein. As a result, the metabolic biochemical markers profile along with the identified SNPs could be used as a candidate and a reference guide for effective characterization of the Holstein, Simmental, and Brown Swiss breeds, leading to the creation of a marker-assisted selection system for production traits in dairy cattle breeds.

Genomic Uniqueness of Local Sheep Breeds From Morocco

Sheep farming is a major source of meat in Morocco and plays a key role in the country's agriculture. This study aims at characterizing the whole-genome diversity and demographic history of the main Moroccan sheep breeds, as well as to identify selection signatures within and between breeds. Whole genome data from 87 individuals representing the five predominant local breeds were used to estimate their level of neutral genetic diversity and to infer the variation of their effective population size over time. In addition, we used two methods to detect selection signatures: either for detecting selective sweeps within each breed separately or by detecting differentially selected regions by contrasting different breeds. We identified hundreds of genomic regions putatively under selection, which related to several biological terms involved in local adaptation or the expression of zootechnical performances such as Growth, UV protection, Cell maturation or Feeding behavior. The results of this study revealed selection signatures in genes that have an important role in traits of interest and increased our understanding of how genetic diversity is distributed in these local breeds. Thus, Moroccan local sheep breeds exhibit both a high genetic diversity and a large set of adaptive variations, and therefore, represent a valuable genetic resource for the conservation of sheep in the context of climate change.

Local adaptations of Mediterranean sheep and goats through an integrative approach

Small ruminants are suited to a wide variety of habitats and thus represent promising study models for identifying genes underlying adaptations. Here, we considered local Mediterranean breeds of goats (n = 17) and sheep (n = 25) from Italy, France and Spain. Based on historical archives, we selected the breeds potentially most linked to a territory and defined their original cradle (i.e., the geographical area in which the breed has emerged), including transhumant pastoral areas. We then used the programs PCAdapt and LFMM to identify signatures of artificial and environmental selection. Considering cradles instead of current GPS coordinates resulted in a greater number of signatures identified by the LFMM analysis. The results, combined with a systematic literature review, revealed a set of genes with potentially key adaptive roles in relation to the gradient of aridity and altitude. Some of these genes have been previously implicated in lipid metabolism (SUCLG2, BMP2), hypoxia stress/lung function (BMPR2), seasonal patterns (SOX2, DPH6) or neuronal function (TRPC4, TRPC6). Selection signatures involving the PCDH9 and KLH1 genes, as well as NBEA/NBEAL1, were identified in both species and thus could play an important adaptive role.

Genetic Variations and Differential DNA Methylation to Face Contrasted Climates in Small Ruminants: An Analysis on Traditionally-Managed Sheep and Goats

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (F_ST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

Naturally Occurring Glycosidases in Milk from Native Cattle Breeds: Activity and Consequences on Free and Protein Bound-Glycans

Little is known about the extent of variation and activity of naturally occurring milk glycosidases and their potential to degrade milk glycans. A multi-omics approach was used to investigate the relationship between glycosidases and important bioactive compounds such as free oligosaccharides and O-linked glycans in bovine milk. Using 4-methylumbelliferone (4-MU) assays activities of eight indigenous glycosidases were determined, and by mass spectrometry and ¹H NMR spectroscopy various substrates and metabolite products were quantified in a subset of milk samples from eight native North European cattle breeds. The results showed a clear variation in glycosidase activities among the native breeds. Interestingly, negative correlations between some glycosidases including β-galactosidase, N-acetyl-β-d-glucosaminidase, certain oligosaccharide isomers as well as O-linked glycans of κ-casein were revealed. Further, a positive correlation was found for free fucose content and α-fucosidase activity (r = 0.37, p-value < 0.001) indicating cleavage of fucosylated glycans in milk at room temperature. The results obtained suggest that milk glycosidases might partially degrade valuable glycans, which would result in lower recovery of glycans and thus represent a loss for the dairy ingredients industry if these activities are pronounced.

Disentangling the Multidimensional Relationship between Livestock Breeds and Ecosystem Services

Simple Summary

Livestock breeds represent the diversity of livestock animals. They participate in the delivery of ecosystem services (ES), i.e., the benefits to humans provided by nature. In recent years, the contribution of livestock breeds to ES has received attention in livestock research. Additionally, there is increasing interest in integrating this knowledge into policies to make agriculture more sustainable. In this work, we elaborate on livestock breed characteristics that are key to the study of livestock breed contributions to ES. Thus, we explore the natural and human factors that have produced livestock breeds as ecologically and culturally mediated entities. In addition, we review the different roles of livestock breeds as biodiversity components. Finally, we examine how livestock breeds participate in livestock system heterogeneity. By integrating these aspects, we might better understand how livestock breeds provide and modulate ES provision and, therefore, how to improve breed conservation and livestock policies toward more sustainable farming.

Abstract

There is an increasing interest in assessing livestock breed contributions to ecosystem services (ES) and including this knowledge in decision making. However, this task has been limited due to the complexity of the multidimensional relationship between livestock diversity and ecosystem services. In this work, we elaborate on the livestock breed characteristics central to developing a comprehensive approach to livestock breed inclusion in the ecosystem services framework. Thus, we explore the multidimensional nature of livestock breeds, i.e., as eco-cultural entities, biodiversity components, and drivers of livestock system heterogeneity and functioning. First, anthropogenic and natural factors have acted jointly to develop breeds as eco-cultural entities. This fact represents an opportunity to move toward farming system sustainability by Nature-Based Solutions and Nature’s Contribution to People paradigms. Second, livestock breeds are components of biodiversity, and as such, can be framed as goods, as final ecosystem services, and as regulators of ecosystem processes. Third, livestock breeds contribute to livestock system heterogeneity and resilience. By integrating these aspects, we might better understand how livestock breeds provide and modulate ecosystem service provision and, therefore, how to improve breed conservation and livestock policies toward farming system sustainability.

Morphometric and microsatellite-based comparative genetic diversity analysis in Bubalus bubalis from North India

To understand the similarities and dissimilarities of a breed structure among different buffalo breeds of North India, it is essential to capture their morphometric variation, genetic diversity, and effective population size. In the present study, diversity among three important breeds, namely, Murrah, Nili-Ravi and Gojri were studied using a parallel approach of morphometric characterization and molecular diversity. Morphology was characterized using 13 biometric traits, and molecular diversity through a panel of 22 microsatellite DNA markers recommended by FAO, Advisory Group on Animal Genetic Diversity, for diversity studies in buffaloes. Canonical discriminate analysis of biometric traits revealed different clusters suggesting distinct genetic entities among the three studied populations. Analysis of molecular variance revealed 81.8% of genetic variance was found within breeds, while 18.2% of the genetic variation was found between breeds. Effective population sizes estimated based on linkage disequilibrium were 142, 75 and 556 in Gojri, Nili-Ravi and Murrah populations, respectively, indicated the presence of sufficient genetic variation and absence of intense selection among three breeds. The Bayesian approach of STRUCTURE analysis (at K = 3) assigned all populations into three clusters with a degree of genetic admixture in the Murrah and Nili-Ravi buffalo populations. Admixture analysis reveals introgression among Murrah and Nili-Ravi breeds while identified the Gojri as unique buffalo germplasm, indicating that there might be a common origin of Murrah and Nili-Ravi buffaloes. The study provides important insights on buffalo breeds of North India that could be utilized in designing an effective breeding strategy, with an appropriate choice of breeds for upgrading local non-descript buffaloes along with conservation of unique germplasm.

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

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

Genome-wide DNA arrays profiling unravels the genetic structure of Iranian sheep and pattern of admixture with worldwide coarse-wool sheep breeds

Archaeological and genetic evidence show that sheep were originally domesticated in area around the North of Zagros mountains, North-west of Iran. The Persian plateau exhibits a variety of native sheep breeds with a common characteristic of coarse-wool production. Therefore, knowledge about the genetic structure and diversity of Iranian sheep and genetic connections with other sheep breeds is of great interest. To this end, we genotyped 154 samples from 11 sheep breeds distributed across Iran with the Ovine Infinium HD SNP 600 K BeadChip array, and analyzed this dataset combined with the retrieved data of 558 samples from 19 worldwide coarse-wool sheep breeds. The average genetic diversity ranged from 0.315 to 0.354, while the F_ST values ranged from 0.016 to 0.177 indicating a low differentiation of Iranian sheep. Analysis of molecular variance showed that 90.21 and 9.79% of the source of variation were related to differences within and between populations, respectively. Our results indicated that the coarse-wool sheep from Europe were clearly different from those of the Asia. Accordingly, the Asiatic mouflon was positioned between Asian and European countries. In addition, we found that the genetic background of Iranian sheep is present in sheep from China and Kyrgyzstan, as well as India. The revealed admixture patterns of the Iranian sheep and other coarse-wool sheep breeds probably resulted from the expansion of nomads and through the Silk Road trade network.

Genome Divergence and Dynamics in the Thin-Tailed Desert Sheep From Sudan

With climate change bound to affect food and feed production, emphasis will shift to resilient and adapted indigenous livestock to sustain animal production. However, indigenous livestock comprise several varieties, strains and ecotypes whose genomes are poorly characterized. Here, we investigated genomic variation in an African thin-tailed Desert Sheep sampled in Sudan, using 600K genotype data generated from 92 individuals representing five ecotypes. We included data from 18 fat-tailed and 45 thin-tailed sheep from China, to investigate shared ancestry and perform comparative genomic analysis. We observed a clear genomic differentiation between the African thin-tailed Desert Sheep and the Chinese thin-tailed and fat-tailed sheep, suggesting a broad genetic structure between the fat-tailed and thin-tailed sheep in general, and that at least two autosomal gene pools comprise the genome profile of the thin-tailed sheep. Further analysis detected two distinct genetic clusters in both the African thin-tailed Desert Sheep and the Chinese thin-tailed sheep, suggesting a fine-scale and complex genome architecture in thin-tailed sheep. Selection signature analysis suggested differences in adaptation, production, reproduction and morphology likely underly the fine-scale genetic structure in the African thin-tailed Desert Sheep. This may need to be considered in designing breeding programs and genome-wide association studies.

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y-linked scaffolds using existing data. A total of 11.5 Mb Y-linked sequences were clustered into 33 scaffolds, and 187 protein-coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole-genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post-Neolithic period.

Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population

The Brazilian Guzerá population originated from a few founders introduced from India. These animals adapted well to the harsh environments in Brazil, were selected for beef, milk, or dual-purpose (beef and milk), and were extensively used to produce crossbred animals. Here, the impact of these historical events with regard to the population structure and genetic diversity in a Guzerá meta-population was evaluated. DNA samples of 744 animals (one dairy, nine dual-purpose, and five beef herds) were genotyped for 21 microsatellite loci. Ho, He, PIC, Fis, Fit, and Fst estimates were obtained considering either farms or lineages as subpopulations. Mean Ho (0.73) and PIC (0.75) suggest that genetic diversity was efficiently conserved. Fit, Fis and Fst values (95% CI) pointed to a low fixation index, and large genetic diversity: Fit (Farms = 0.021-0.100; lineages = 0.021-0.100), Fis (Farms = -0.007-0.076; lineages = -0.014-0.070), and Fst (Farms = 0.0237-0.032; lineages = 0.029-0.038). The dual-purpose herds/selection lines are the most uniform subpopulation, while the beef one preserved larger amounts of genetic diversity among herds. In addition, the dairy herd showed to be genetically distant from other herds. Taken together, these results suggest that this Guzerá meta-population has high genetic diversity, a low degree of population subdivision, and a low inbreeding level.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Genetic characterisation of non-descript cattle populations in communal areas of South Africa

Context Indigenous cattle breeds represent an important genetic resource for livelihood of communal-area inhabitants. Indigenous breeds have the ability to withstand harsh climatic conditions, can adapt genetically to poor-quality forages and are resistant to parasites and diseases. These unique traits possessed by indigenous breeds are under threat because of unrestrained crossing with exotic commercial breeds, and this can lead to total loss of a breed. Aims The study was conducted to assess the genetic diversity and population structure of South African non-descript communal beef cattle populations by using 25 microsatellite markers. Methods Unrelated and non-descript animals (n = 150) were sampled from communal areas from five (5) provinces of South Africa, namely, Eastern Cape, KwaZulu–Natal, Limpopo, Mpumalanga and the North West, with 30 samples per breed taken. Six (6) known cattle breeds (n = 180) were used as a reference population. This included Angus, Afrikaner, Bonsmara, Brahman, Drakensberger and the Nguni, with 30 samples per breed. Key results High level of genetic diversity was found across the five non-descript populations, with an average heterozygosity of 75%. The Limpopo population was found to be the most diverse population, with the highest average number of alleles (8.5) and heterozygosity (ranging between observed heterozygosity of 70% and expected heterozygosity of 79%). STRUCTURE software assigned populations (2 ≤ K ≤ 20), with the most probable cluster being at K = 7. The Eastern Cape, KwaZulu–Natal and Limpopo populations had genetic material similar to those possessed by the Nguni and Bonsmara reference populations. Conclusions Results from the study showed that most genetic differentiation occurred within populations rather than among populations, and this might be due to the fact that there is no selection for or against any specific production trait expressed in the populations. Implications The obtained information will serve as a baseline for the development and implementation of sound breeding programs that will assist in controlling the gene flow, so as to lower the possible genetic dilution of the currently available genetic material.

Recent land use and management changes decouple the adaptation of livestock diversity to the environment

Native livestock breeds, i.e. those autochthonous to a specific region, are locally adapted domesticated animals that conserve genetic resources, guaranty food security and provide agroecosystem services. Native breeds are largely threatened worldwide by agricultural intensification and rural areas abandonment processes related to recent changes in production schemes and planning. Yet, our gap of knowledge regarding livestock breed-environment relationships may prevent the design of successful conservation measures. In this work, we analyse the links between livestock diversity -i.e. richness of native breeds- and a selection of environmental factors that express at broad scales, with a temporal perspective. We compare native breeds distributional patterns before and after the agricultural intensification, in the context of land-use change in mainland Spain. Our results confirm the existence of strong associations between the distribution of native livestock breeds and environmental factors. These links, however, weaken for contemporary distributions. In fact, changes in breed distribution reflect a shift towards more productive environments. Finally, we found that the areas having higher breed richness are undergoing land abandonment processes. Succeeding in the conservation of threatened native breeds will require going beyond merely genetic and production-oriented views. Ecological and sociocultural perspectives should also be accounted for as global change processes are determinant for livestock agrobiodiversity.

Genetic material from buffalo and cattle: crucial importance in the formalization of bilateral trade between India and Brazil

The trade in live animals between India and Brazil dates from the late nineteenth century when European travellers traded animals of Indian origin for display in zoos. Considering the origin of coffee and sugar cane, as well as the expertise related to mineral evaluation, we need to consider that India was involved in important economic cycles of Brazil, even indirectly. This virtuous flow of trade has been maintained and intensified throughout modern history, especially after these two nations gained political independence from their colonisers, thereby becoming independent in mercantile affairs. This paper addresses the main points related to the use of animals of Indian origin in Brazil. We revisit some of the historical aspects of the process of colonisation of Brazil, as well as the importation of animals from India. The restrictions imposed on this process due to the occurrence of diseases in cattle and buffalo in India will be examined. At the end of the text, emphasis will be given to the risks of introducing exotic diseases into Brazil.

Population genetic structure of two herds of Aberdeen Angus cattle breed in Colombia

Background: Two biotypes of Aberdeen Angus cattle breed, known as Old Type and New Type, that differ in their origin and beef production are formally recognized. In Colombia, this breed has been commercialized for approximately 80 years. Studies on the origin, kinship and levels of genetic diversity of this breed in Colombian herds are scarce, yet important for planning crossing and management strategies. Objective: To measure the genetic diversity and structure of two Colombian herds of Old Type and New Type biotypes of Aberdeen Angus from Huila and Cundinamarca provinces and assess mitochondrial introgression with other breeds. Methods: A set of ten microsatellites and sequences of the Mitochondrial Control Region were characterized. Estimators of genetic diversity and population differentiation along with tests of population assignment were applied. Results: Nuclear loci were highly polymorphic as shown by the Polymorphic Information Content (0.599) and the Probability of Identity (1.896 10-08). Both populations were highly diverse and clearly differentiated into two groups corresponding to the Old Type and New Type phenotypes. In contrast, mitochondrial data failed to distinguish these two groups and showed extensive admixture. Conclusions: This study optimized a set of ten highly polymorphic nuclear markers that may be used for parentage and population genetic studies of Aberdeen Angus. Genetic differentiation in these loci agreed with phenotypic differences of the Old and New Types. However, mitochondrial data indicated ancestry of multiple European breeds in the origin of Colombian Aberdeen Angus.

Genetic signature of hybridization between Chinese spot-billed ducks and domesticated ducks

In this study, we analyzed 93 whole genomes from Chinese spot-billed ducks (CSB), meat-type ducks (MET), and egg and dual purpose-type ducks (EDT) to characterize the genetic material flowing between the CSB and modern ducks. Using a frequency of shared identical-by-descent method, approximately 10.9 Mb introgression segments containing 140 genes were identified showing the signatures of introgression between CSB and EDT. Meanwhile, nearly 10.6 M introgression regions containing 149 genes were identified between CSB and MET. Based on the haplotypes tree of each segment, we found that the introgression between CSB and domesticated ducks was asymmetric with a high level of gene flow from domestic to CSB and a low level of migration in the opposite direction. Moreover, we identified several genes that were introgressions from CSB and showed the signature of positive selection, which may contribute to the breeding of modern ducks. Our results provide new insight into the evolution and breeding history of domestic ducks and may be useful for the future management of wild and domestic duck populations.

Genetic diversity status of only registered cattle breed of Chhattisgarh-Kosali

Assessment of diversity is essential for germplasm characterization and management. Kosali is the first and only registered breed of Chhattisgarh state which immensely contributes to the local economy of 70% farmers. Genomic microsatellite markers being valuable tool for estimating genetic diversity were selected for exploring existing genetic variability in Kosali cattle population. The standard metrics of genomic diversity detected high variability in this breed. All the loci were polymorphic resulting in 297 alleles. Mean values of observed and expected number of alleles were 11.423±0.877 and 4.989±0.372, respectively. Similarly, higher mean values of observed heterozygosity (0.693±0.031) also corroborated with the allelic diversity. Mean expected heterozygosity (0.765±0.02) under Hardy-Weinberg equilibrium (HWE) was higher than the observed values indicating some deviations from assumptions of the model. It can be attributed to the forces such as inbreeding. In fact a positive FIS value of 0.088±0.038 indicated some heterozygote deficiency in the population. Bottleneck analysis indicated that Kosali cattle have not suffered any population bottleneck event during evolutionary trajectory. This study is first to report the genetic diversity status of Kosali cattle based on microsatellite markers.The results imply the necessity of management programs in order to conserve the existing genetic variation and to avoid any escalation of inbreeding.

Genomic diversity and population structure of the indigenous Greek and Cypriot cattle populations

Background

The indigenous cattle populations from Greece and Cyprus have decreased to small numbers and are currently at risk of extinction due to socio-economic reasons, geographic isolation and crossbreeding with commercial breeds. This study represents the first comprehensive genome-wide analysis of 10 indigenous cattle populations from continental Greece and the Greek islands, and one from Cyprus, and compares them with 104 international breeds using more than 46,000 single nucleotide polymorphisms (SNPs).

Results

We estimated several parameters of genetic diversity (e.g. heterozygosity and allelic diversity) that indicated a severe loss of genetic diversity for the island populations compared to the mainland populations, which is mainly due to the declining size of their population in recent years and subsequent inbreeding. This high inbreeding status also resulted in higher genetic differentiation within the Greek and Cyprus cattle group compared to the remaining geographical breed groups. Supervised and unsupervised cluster analyses revealed that the phylogenetic patterns in the indigenous Greek breeds were consistent with their geographical origin and historical information regarding crosses with breeds of Anatolian or Balkan origin. Cyprus cattle showed a relatively high indicine ancestry. Greek island populations are placed close to the root of the tree as defined by Gir and the outgroup Yak, whereas the mainland breeds share a common historical origin with Buša. Unsupervised clustering and D-statistics analyses provided strong support for Bos indicus introgression in almost all the investigated local cattle breeds along the route from Anatolia up to the southern foothills of the Alps, as well as in most cattle breeds along the Apennine peninsula to the southern foothills of the Alps.

Conclusions

All investigated Cyprus and Greek breeds present complex mosaic genomes as a result of historical and recent admixture events between neighbor and well-separated breeds. While the contribution of some mainland breeds to the genetic diversity pool seems important, some island and fragmented mainland breeds suffer from a severe decline of population size and loss of alleles due to genetic drift. Conservation programs that are a compromise between what is feasible and what is desirable should focus not only on the still highly diverse mainland breeds but also promote and explore the conservation possibilities for island breeds.