An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

Unveiling population structure and selection signatures of riverine and genetically improved rohu, Labeo rohita using genome wide SNPs

BACKGROUND

Captive breeding, along with artificial selection can significantly impact population structure by influencing allele frequencies and driving populations towards specific adaptation. Selective sweeps are powerful forces in shaping genetic variation within populations and can drive rapid spread of beneficial alleles while simultaneously reducing genetic diversity in localized regions of the genome. The present work was undertaken to assess the genetic structure and consequences of artificial selection in 10th generation of genetically improved rohu by comparing with wild populations.

METHODS AND RESULTS

The present study used 11,022 high-quality genome wide SNPs to compare the population genetic structure and signatures of selection between Jayanti rohu population and its wild counterpart. Outlier analysis revealed presence of 14 adaptive SNPs, out of which 5 were classified to be under decisive selection pressure. Notably, Jayanti rohu (JR) displayed 297 private alleles exclusive to its population. Chromosomes 7 and 16 emerged as potential hotspots containing a majority of the identified SNPs. Structure and principal component analysis revealed two distinct clusters, effectively distinguishing the JR and wild rohu populations. Phylogenetic analysis indicated a separate cluster of JR population distant from wild groups.

CONCLUSION

The results of present study shall help in elucidating patterns of genetic variation and characterizing selection signatures associated with captive bred and natural populations of rohu. The genomic resources generated through this work shall be helpful in improving the traceability of selectively bred germplasm for developing future strategies of genetic management.

Signatures of selection in Angus and Hanwoo beef cattle using imputed whole genome sequence data

In this study, we detected signatures of selection in Hanwoo and Angus beef cattle using allele frequency and haplotype-based methods based on imputed whole genome sequence variants. Our dataset included 13,202 Angus animals with 10,057,633 imputed SNPs and 10,437 Hanwoo animals with 13,241,550 imputed SNPs. The dataset was subset down to 6,873,624 SNPs in common between the two populations to identify within population (runs of homozygosity, extended haplotype homozygosity) and between population signals of selection (allele fixation index, extended haplotype homozygosity). Assuming these selection signals were complementary to each other, they were combined into a decorrelated composite of multiple signals to identify regions under selection for each of the breeds. 27 genomic regions spanning 25.15 Mb and harboring 360 genes were identified in Angus on chromosomes 1,3, 4, 5, 6, 7, 8, 12, 13, 14, 16, 20, 21 and 28. Similarly, in Hanwoo, 59 genes and 17 genomic regions spanning 5.21 Mb on chromosomes 2, 4, 5, 6, 7, 8, 9, 10, 13, 17, 20 and 24 were identified. Apart from a small region on chromosome 13, there was no major overlap of selection signals between the two breeds reflecting their largely different selection histories, environmental challenges, breeding objectives and breed characteristics. Positional candidate genes identified in selected genomic regions in Angus have been previously associated with growth, immunity, reproductive development, feed efficiency and adaptation to environment while the candidate genes identified in Hanwoo included important genes regulating meat quality, fat deposition, cholesterol metabolism, lipid synthesis, neuronal development, and olfactory reception.

Beef Cattle Genome Project: Advances in Genome Sequencing, Assembly, and Functional Genes Discovery

Beef is a major global source of protein, playing an essential role in the human diet. The worldwide production and consumption of beef continue to rise, reflecting a significant trend. However, despite the critical importance of beef cattle resources in agriculture, the diversity of cattle breeds faces severe challenges, with many breeds at risk of extinction. The initiation of the Beef Cattle Genome Project is crucial. By constructing a high-precision functional annotation map of their genome, it becomes possible to analyze the genetic mechanisms underlying important traits in beef cattle, laying a solid foundation for breeding more efficient and productive cattle breeds. This review details advances in genome sequencing and assembly technologies, iterative upgrades of the beef cattle reference genome, and its application in pan-genome research. Additionally, it summarizes relevant studies on the discovery of functional genes associated with key traits in beef cattle, such as growth, meat quality, reproduction, polled traits, disease resistance, and environmental adaptability. Finally, the review explores the potential of telomere-to-telomere (T2T) genome assembly, structural variations (SVs), and multi-omics techniques in future beef cattle genetic breeding. These advancements collectively offer promising avenues for enhancing beef cattle breeding and improving genetic traits.

Uncovering the architecture of selection in two Bos taurus cattle breeds

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle.

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 differentiation following recent domestication events: A study of farmed Nile tilapia (Oreochromis niloticus) populations

Nile tilapia (Oreochromis niloticus) is among the most farmed finfish worldwide, distributed across different environmental conditions. Its wide distribution has mainly been facilitated by several breeding programs and widespread dissemination of genetically improved strains. In the first Nile tilapia study exploiting a whole-genome pooled sequencing (Poolseq) approach, we identified the genetic structure and signatures of selection in diverse, farmed Nile tilapia populations, with a particular focus on the GIFT strain, developed in the 1980s, and currently managed by WorldFish (GIFTw). We also investigated important farmed strains from The Philippines and Africa. Using both SNP array data and Poolseq SNPs, we characterized the population structure of these samples. We observed the greatest separation between the Asian and African populations and greater admixture in the Asian populations than in the African ones. We also established that the SNP array data were able to successfully resolve relationships between these diverse Nile tilapia populations. The Poolseq data identified genomic regions with high levels of differentiation (F _ST) between GIFTw and the other populations. Gene ontology terms associated with mesoderm development were significantly enriched in the genes located in these regions. A region on chromosome Oni06 was genetically differentiated in pairwise comparisons between GIFTw and all other populations. This region contains genes associated with muscle-related traits and overlaps with a previously published QTL for fillet yield, suggesting that these traits may have been direct targets for selection on GIFT. A nearby region was also identified using XP-EHH to detect genomic differentiation using the SNP array data. Genomic regions with high or extended homozygosity within each population were also identified. This study provides putative genomic landmarks associated with the recent domestication process in several Nile tilapia populations, which could help to inform their genetic management and improvement.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Genome-wide scan for selective sweeps identifies novel loci associated with resistance to mastitis in German Holstein cattle

Domestication and selection significantly changed phenotypic and behavioural traits in modern domestic animals. In this study, to identify the genomic regions associated with mastitis, genomic data of German Holstein dairy cattle were analysed. The samples were genotyped using the Bovine 50 K SNP chip. For each defined healthy and sick group, 133 samples from 13,276 genotyped dairy cows were selected based on mastitis random residual effects. Grouping was done to infer selection signatures based on XP-EHH statistic. The results revealed that for the top 0.01 percentile of the obtained XP-EHH values, five genomic regions on chromosomes 8, 11, 12, 14 and 26 of the control group, and four regions on chromosomes 3, 4 (two regions) and 22 of the case group, have been under selection. Also, consideration of the top 0.1 percentile of the XP-EHH values, clarified 21 and 15 selective sweeps in the control and case group, respectively. This study identified some genomic regions containing potential candidate genes associated with resistance and susceptibility to mastitis, immune system and inflammation, milk traits, udder morphology and different types of cancers. In addition, these regions overlap with some quantitative trait loci linked to clinical mastitis, immunoglobulin levels, somatic cell score, udder traits, milk fat and protein, milk yield, milking speed and veterinary treatments. It is noteworthy that we found two regions in the healthy group (on chromosomes 12 and 14) with strong signals, which were not described previously. It is likely that future research could link these identified genomic regions to mastitis. The results of the current study contribute to the identification of causal mutations, genomic regions and genes affecting mastitis incidence in dairy cows.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Integrative QTL mapping and selection signatures in Groningen White Headed cattle inferred from whole-genome sequences

Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2-71.3 Mb), BTA5 (10.0-19.7 Mb), BTA20 (10.0-19.9 and 20.0-22.7 Mb), and BTA25 (0.5-9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., CWC27 and CLUAP1), ultraviolet protection (e.g., ERCC8), and pigmentation (e.g. PDE4D) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.

Genome-wide scan for signatures of selection in the Brangus cattle genome

Brangus is a composite cattle breed developed with the objective of combining the advantages of Angus and Zebuine breeds (Brahman, mainly) in tropical climates. The aim of this work was to estimate breed composition both genome-wide and locally, at the chromosome level, and to uncover genomic regions evidencing positive selection in the Argentinean Brangus population/nucleus. To do so, we analysed marker data from 478 animals, including Brangus, Angus and Brahman. Average breed composition was 35.0% ± 9.6% of Brahman, lower than expected according to the theoretical fractions deduced by the usual cross-breeding practice in this breed. Local ancestry analysis evidenced that breed composition varies between chromosomes, ranging from 19.6% for BTA26 to 56.1% for BTA5. Using approaches based on allelic frequencies and linkage disequilibrium, genomic regions with putative selection signatures were identified in several chromosomes (BTA1, BTA5, BTA6 and BTA14). These regions harbour genes involved in horn development, growth, lipid metabolism, reproduction and immune response. We argue that the overlapping of a chromosome segment originated in one of the parental breeds and over-represented in the sample with the location of a signature of selection constitutes evidence of a selection process that has occurred in the breed since its take off in the 1950s. In this regard, our results could contribute to the understanding of the genetic mechanisms involved in cross-bred cattle adaptation and productivity in tropical environments.

Analysis on the desert adaptability of indigenous sheep in the southern edge of Taklimakan Desert

The southern margin of the Taklimakan Desert is characterized by low rainfall, heavy sandstorms, sparse vegetation and harsh ecological environment. The indigenous sheep in this area are rich in resources, with the advantages of perennial estrus and good resistance to stress in most sheep. Exploring the molecular markers of livestock adaptability in this environment will provide the molecular basis for breeding research to cope with extreme future changes in the desert environment. In this study, we analyzed the population genetic structure and linkage imbalance of five sheep breeds with three different agricultural geographic characteristics using four complementary genomic selection signals: fixation index (FST), cross-population extended haplotype homozygosity (xp-EHH), Rsb (extended haplotype homozygosity between-populations) and iHS (integrated haplotype homozygosity score). We used Illumina Ovine SNP 50K Genotyping BeadChip Array, and gene annotation and enrichment analysis were performed on selected regions of the obtained genome. The ovary of Qira Black sheep (Follicular phase, Luteal phase, 30th day of pregnancy, 45th day of pregnancy) was collected, and the differentially expressed genes were screened by transcriptomic sequencing. Genome-wide selective sweep results and transcriptome data were combined for association analysis to obtain candidate genes associated with perennial estrus and stable reproduction. In order to verify the significance of the results, 15 resulting genes were randomly selected for fluorescence quantitative analysis. The results showed that Dolang sheep and Qira Black sheep evolved from Kazak sheep. Linkage disequilibrium analysis showed that the decay rate of sheep breeds in the Taklimakan Desert was higher than that in Yili grassland. The signals of FST, xp-EHH, Rsb and iHS detected 526, 332, 308 and 408 genes, respectively, under the threshold of 1% and 17 overlapping genes under the threshold of 5%. A total of 29 genes were detected in association analysis of whole-genome and transcriptome data. This study reveals the genetic mechanism of perennial estrus and environmental adaptability of indigenous sheep breeds in the Taklimakan Desert. It provides a theoretical basis for the conservation and exploitation of genetic resources of indigenous sheep breeds in extreme desert environment. This provides a new perspective for the quick adaptation of sheep and other mammals to extreme environments and future climate changes.

Along the Bos taurus genome, uncover candidate imprinting control regions

Background

In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin.

Results

Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression.

Conclusion

Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08694-3.

Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data

BACKGROUND

Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide.  RESULTS: Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, F_ST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2).

CONCLUSION

In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak

Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima's D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Powerful detection of polygenic selection and evidence of environmental adaptation in US beef cattle

Selection on complex traits can rapidly drive evolution, especially in stressful environments. This polygenic selection does not leave intense sweep signatures on the genome, rather many loci experience small allele frequency shifts, resulting in large cumulative phenotypic changes. Directional selection and local adaptation are changing populations; but, identifying loci underlying polygenic or environmental selection has been difficult. We use genomic data on tens of thousands of cattle from three populations, distributed over time and landscapes, in linear mixed models with novel dependent variables to map signatures of selection on complex traits and local adaptation. We identify 207 genomic loci associated with an animal's birth date, representing ongoing selection for monogenic and polygenic traits. Additionally, hundreds of additional loci are associated with continuous and discrete environments, providing evidence for historical local adaptation. These candidate loci highlight the nervous system's possible role in local adaptation. While advanced technologies have increased the rate of directional selection in cattle, it has likely been at the expense of historically generated local adaptation, which is especially problematic in changing climates. When applied to large, diverse cattle datasets, these selection mapping methods provide an insight into how selection on complex traits continually shapes the genome. Further, understanding the genomic loci implicated in adaptation may help us breed more adapted and efficient cattle, and begin to understand the basis for mammalian adaptation, especially in changing climates. These selection mapping approaches help clarify selective forces and loci in evolutionary, model, and agricultural contexts.

Detection of Genomic Regions with Pleiotropic Effects for Growth and Carcass Quality Traits in the Rubia Gallega Cattle Breed

Simple Summary

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. We have developed a ssGWAS by backsolving the SNP effects after implementing a ssGBLUP. The results showed an apparent heterogeneity of the additive genetic variance across the genome. Some of the genomic regions explaining the most of this additive variance were shared across traits, indicating the presence of pleiotropic effects, which were reflected in their genetic correlations.

Abstract

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. In recent years, genomic evaluation has been implemented by using a ssGBLUP (single-step Genomic Best Linear Unbiased Prediction). This procedure can reparameterized to perform ssGWAS (single-step Genome Wide Association Studies) by backsolving the SNP (single nucleotide polymorphisms) effects. Therefore, the objective of this study was to identify genomic regions associated with the genetic variability in growth and carcass quality traits. We implemented a ssGBLUP by using a database that included records for Birth Weight (BW-327,350 records-), Weaning Weight (WW-83,818-), Cold Carcass Weight (CCW-91,621-), Fatness (FAT-91,475-) and Conformation (CON-91,609-). The pedigree included 464,373 individuals, 2449 of which were genotyped. After a process of filtering, we ended up using 43,211 SNP markers. We used the GBLUP and SNPBLUP model equivalences to obtain the effects of the SNPs and then calculated the percentage of variance explained by the regions of the genome between 1 Mb. We identified 7 regions of the genome for CCW; 8 regions for BW, WW, FAT and 9 regions for CON, which explained the percentage of variance above 0.5%. Furthermore, a number of the genome regions had pleiotropic effects, located at: BTA1 (131–132 Mb), BTA2 (1–11 Mb), BTA3 (32–33 Mb), BTA6 (36–38 Mb), BTA16 (24–26 Mb), and BTA 21 (56–57 Mb). These regions contain, amongst others, the following candidate genes: NCK1, MSTN, KCNA3, LCORL, NCAPG, and RIN3.

Whole-Genome Sequence Data Suggest Environmental Adaptation of Ethiopian Sheep Populations

Great progress has been made over recent years in the identification of selection signatures in the genomes of livestock species. This work has primarily been carried out in commercial breeds for which the dominant selection pressures are associated with artificial selection. As agriculture and food security are likely to be strongly affected by climate change, a better understanding of environment-imposed selection on agricultural species is warranted. Ethiopia is an ideal setting to investigate environmental adaptation in livestock due to its wide variation in geo-climatic characteristics and the extensive genetic and phenotypic variation of its livestock. Here, we identified over three million single nucleotide variants across 12 Ethiopian sheep populations and applied landscape genomics approaches to investigate the association between these variants and environmental variables. Our results suggest that environmental adaptation for precipitation-related variables is stronger than that related to altitude or temperature, consistent with large-scale meta-analyses of selection pressure across species. The set of genes showing association with environmental variables was enriched for genes highly expressed in human blood and nerve tissues. There was also evidence of enrichment for genes associated with high-altitude adaptation although no strong association was identified with hypoxia-inducible-factor (HIF) genes. One of the strongest altitude-related signals was for a collagen gene, consistent with previous studies of high-altitude adaptation. Several altitude-associated genes also showed evidence of adaptation with temperature, suggesting a relationship between responses to these environmental factors. These results provide a foundation to investigate further the effects of climatic variables on small ruminant populations.

Insights into adaption and growth evolution: a comparative genomics study on two distinct cattle breeds from Northern and Southern China

Mongolian cattle (MG, Bos taurus) and Minnan cattle (MN, Bos indicus) are two different breeds of Chinese indigenous cattle, representing North type and South type, respectively. However, their value and potential have not yet been discovered at the genomic level. In this study, 26 individuals of MN and MG were sequenced for the first time at an average of 13.9- and 12.8-fold, respectively. Large numbers of different variations were identified. In addition, the analyses of phylogenetic and population structure showed that these two cattle breeds are distinct from each other, and results of linkage disequilibrium analysis revealed that these two cattle breeds have undergone various degrees of intense natural or artificial selection. Subsequently, 496 and 306 potential selected genes (PSRs) were obtained in MN and MG, containing 1,096 and 529 potential selected genes (PSGs), respectively. These PSGs, together with the analyzed copy number variation (CNV)-related genes, showed potential relations with their phenotypic characteristics, including environmental adaptability (e.g., DVL2, HSPA4, CDHR4), feed efficiency (e.g., R3HDM1, PLAG1, XKR4), and meat/milk production (e.g., PDHB, LEMD3, APOF). The results of this study help to gain new insights into the genetic characteristics of two distinct cattle breeds and will contribute to future cattle breeding.

The Genomic Variation in the Aosta Cattle Breeds Raised in an Extensive Alpine Farming System

The Aosta Red Pied (Valdostana Pezzata Rossa (VRP)), the Aosta Black Pied (Valdostana Pezzata Nera (VBP)) and the Aosta Chestnut (Valdostana Castana (CAS)) are dual-purpose cattle breeds (meat and milk), very well adapted to the harsh environmental conditions of alpine territories: their farming is in fact characterized by summer pasture at very high altitude. A total of 728 individuals were genotyped with the GeenSeek Genomic Profiler® (GGP) Bovine 150K Illumina SNP chip as a part of the DUALBREEDING-PSRN Italian-funded research project. The genetic diversity among populations showed that the three breeds are distinct populations based on the F_ST values, ADMIXTURE and Principal Component Analysis (PCA) results. Runs of Homozygosity (ROH) were obtained for the three populations to disclose recent autozygosity. The genomic inbreeding based on the ROH was calculated and coupled with information derived from the F (inbreeding coefficient) and F_ST parameters. The mean F_ROH values were low: CAS = 0.06, VBP = 0.05 and VRP = 0.07, while the average F values were -0.003, -0.01 and -0.003, respectively. The annotation and enrichment analysis, performed in the identified most frequent ROH (TOP_ROH), showed genes that can be linked to the resilience capacity of these populations to harsh environmental farming conditions, and to the peculiar characteristics searched for by farmers in each breed.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds

Background

Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions.

Results

Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C).

Conclusions

We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Signatures of positive selection underlying beef production traits in Korean cattle breeds

The difference in the breeding programs and population history may have diversely shaped the genomes of Korean native cattle breeds. In the absence of phenotypic data, comparisons of breeds that have been subjected to different selective pressures can aid to identify genomic regions and genes controlling qualitative and complex traits. In this study to decipher genetic variation and identify evidence of divergent selection, 3 Korean cattle breeds were genotyped using the recently developed high-density GeneSeek Genomic Profiler F250 (GGP-F250) array. The three Korean cattle breeds clustered according to their coat color phenotypes and breeding programs. The Heugu breed reliably showed smaller effective population size at all generations considered. Across the autosomal chromosomes, 113 and 83 annotated genes were identified from Hanwoo-Chikso and Hanwoo-Heugu comparisons, respectively of which 16 genes were shared between the two pairwise comparisons. The most important signals of selection were detected on bovine chromosomes 14 (24.39–25.13 Mb) and 18 (13.34–15.07 Mb), containing genes related to body size, and coat color (XKR4, LYN, PLAG1, SDR16C5, TMEM68, CDH15, MC1R, and GALNS). Some of the candidate genes are also associated with meat quality traits (ACSF3, EIF2B1, BANP, APCDD1, and GALM) and harbor quantitative trait locus (QTL) for beef production traits. Further functional analysis revealed that the candidate genes (DBI, ACSF3, HINT2, GBA2, AGPAT5, SCAP, ELP6, APOB, and RBL1) were involved in gene ontology (GO) terms relevant to meat quality including fatty acid oxidation, biosynthesis, and lipid storage. Candidate genes previously known to affect beef production and quality traits could be used in the beef cattle selection strategies.

High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data

Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.

Genomic Signatures After Five Generations of Intensive Selective Breeding: Runs of Homozygosity and Genetic Diversity in Representative Domestic and Wild Populations of Turbot (Scophthalmus maximus)

Massive genotyping of single nucleotide polymorphisms (SNP) has opened opportunities for analyzing the way in which selection shapes genomes. Artificial or natural selection usually leaves genomic signatures associated with selective sweeps around the responsible locus. Strong selective sweeps are most often identified either by lower genetic diversity than the genomic average and/or islands of runs of homozygosity (ROHi). Here, we conducted an analysis of selective sweeps in turbot (Scophthalmus maximus) using two SNP datasets from a Northeastern Atlantic population (36 individuals) and a domestic broodstock (46 individuals). Twenty-six families (∼ 40 offspring per family) from this broodstock and three SNP datasets applying differing filtering criteria were used to adjust ROH calling parameters. The best-fitted genomic inbreeding estimate (F_ROH) was obtained by the sum of ROH longer than 1 Mb, called using a 21,615 SNP panel, a sliding window of 37 SNPs and one heterozygous SNP per window allowed. These parameters were used to obtain the ROHi distribution in the domestic and wild populations (49 and 0 ROHi, respectively). Regions with higher and lower genetic diversity within each population were obtained using sliding windows of 37 SNPs. Furthermore, those regions were mapped in the turbot genome against previously reported genetic markers associated with QTL (Quantitative Trait Loci) and outlier loci for domestic or natural selection to identify putative selective sweeps. Out of the 319 and 278 windows surpassing the suggestive pooled heterozygosity thresholds (ZHp) in the wild and domestic population, respectively, 78 and 54 were retained under more restrictive ZHp criteria. A total of 116 suggestive windows (representing 19 genomic regions) were linked to either QTL for production traits, or outliers for divergent or balancing selection. Twenty-four of them (representing 3 genomic regions) were retained under stricter ZHp thresholds. Eleven QTL/outlier markers were exclusively found in suggestive regions of the domestic broodstock, 7 in the wild population and one in both populations; one (broodstock) and two (wild) of those were found in significant regions retained under more restrictive ZHp criteria in the broodstock and the wild population, respectively. Genome mining and functional enrichment within regions associated with selective sweeps disclosed relevant genes and pathways related to aquaculture target traits, including growth and immune-related pathways, metabolism and response to hypoxia, which showcases how this genome atlas of genetic diversity can be a valuable resource to look for candidate genes related to natural or artificial selection in turbot populations.

Linkage Disequilibrium, Effective Population Size and Genomic Inbreeding Rates in American Mink Using Genotyping-by-Sequencing Data

Knowledge of linkage disequilibrium (LD) patterns is necessary to determine the minimum density of markers required for genomic studies and to infer historical changes as well as inbreeding events in the populations. In this study, we used genotyping-by-sequencing (GBS) approach to detect single nucleotide polymorphisms (SNPs) across American mink genome and further to estimate LD, effective population size (Ne), and inbreeding rates based on excess of homozygosity (F_HOM) and runs of homozygosity (ROH). A GBS assay was constructed based on the sequencing of ApeKI-digested libraries from 285 American mink using Illumina HiSeq Sequencer. Data of 13,321 SNPs located on 46 scaffolds was used to perform LD analysis. The average LD (r² ± SD) between adjacent SNPs was 0.30 ± 0.35 over all scaffolds with an average distance of 51 kb between markers. The average r² < 0.2 was observed at inter-marker distances of >40 kb, suggesting that at least 60,000 informative SNPs would be required for genomic selection in American mink. The Ne was estimated to be 116 at five generations ago. In addition, the most rapid decline of population size was observed between 100 and 200 generations ago. Our results showed that short extensions of homozygous genotypes (500 kb to 1 Mb) were abundant across the genome and accounted for 33% of all ROH identified. The average inbreeding coefficient based on ROH longer than 1 Mb was 0.132 ± 0.042. The estimations of F_HOM ranged from −0.44 to 0.34 among different samples with an average of 0.15 over all individuals. This study provided useful insights to determine the density of SNP panel providing enough statistical power and accuracy in genomic studies of American mink. Moreover, these results confirmed that GBS approach can be considered as a useful tool for genomic studies in American mink.

Assessing genomic diversity and signatures of selection in Original Braunvieh cattle using whole-genome sequencing data

Background

Autochthonous cattle breeds are an important source of genetic variation because they might carry alleles that enable them to adapt to local environment and food conditions. Original Braunvieh (OB) is a local cattle breed of Switzerland used for beef and milk production in alpine areas. Using whole-genome sequencing (WGS) data of 49 key ancestors, we characterize genomic diversity, genomic inbreeding, and signatures of selection in Swiss OB cattle at nucleotide resolution.

Results

We annotated 15,722,811 SNPs and 1,580,878 Indels including 10,738 and 2763 missense deleterious and high impact variants, respectively, that were discovered in 49 OB key ancestors. Six Mendelian trait-associated variants that were previously detected in breeds other than OB, segregated in the sequenced key ancestors including variants causal for recessive xanthinuria and albinism. The average nucleotide diversity (1.6  × 10^− 3) was higher in OB than many mainstream European cattle breeds. Accordingly, the average genomic inbreeding derived from runs of homozygosity (ROH) was relatively low (F_ROH = 0.14) in the 49 OB key ancestor animals. However, genomic inbreeding was higher in OB cattle of more recent generations (F_ROH = 0.16) due to a higher number of long (> 1 Mb) runs of homozygosity. Using two complementary approaches, composite likelihood ratio test and integrated haplotype score, we identified 95 and 162 genomic regions encompassing 136 and 157 protein-coding genes, respectively, that showed evidence (P < 0.005) of past and ongoing selection. These selection signals were enriched for quantitative trait loci related to beef traits including meat quality, feed efficiency and body weight and pathways related to blood coagulation, nervous and sensory stimulus.

Conclusions

We provide a comprehensive overview of sequence variation in Swiss OB cattle genomes. With WGS data, we observe higher genomic diversity and less inbreeding in OB than many European mainstream cattle breeds. Footprints of selection were detected in genomic regions that are possibly relevant for meat quality and adaptation to local environmental conditions. Considering that the population size is low and genomic inbreeding increased in the past generations, the implementation of optimal mating strategies seems warranted to maintain genetic diversity in the Swiss OB cattle population.

Genetic and genomic analyses underpin the feasibility of concomitant genetic improvement of milk yield and mastitis resistance in dairy sheep

Milk yield is the most important dairy sheep trait and constitutes the key genetic improvement goal via selective breeding. Mastitis is one of the most prevalent diseases, significantly impacting on animal welfare, milk yield and quality, while incurring substantial costs. Our objectives were to determine the feasibility of a concomitant genetic improvement programme for enhanced milk production and resistance to mastitis. Individual records for milk yield, and four mastitis-related traits (milk somatic cell count, California Mastitis Test score, total viable bacterial count in milk and clinical mastitis presence) were collected monthly throughout lactation for 609 ewes of the Chios breed. All ewes were genotyped with a mastitis specific custom-made 960 single nucleotide polymorphism (SNP) array. We performed targeted genomic association studies, (co)variance component estimation and pathway enrichment analysis, and characterised gene expression levels and the extent of allelic expression imbalance. Presence of heritable variation for milk yield was confirmed. There was no significant genetic correlation between milk yield and mastitis traits. Environmental factors appeared to favour both milk production and udder health. There were no overlapping of SNPs associated with mastitis resistance and milk yield in Chios sheep. Furthermore, four distinct Quantitative Trait Loci (QTLs) affecting milk yield were detected on chromosomes 2, 12, 16 and 19, in locations other than those previously identified to affect mastitis resistance. Five genes (DNAJA1, GHR, LYPLA1, NUP35 and OXCT1) located within the QTL regions were highly expressed in both the mammary gland and milk transcriptome, suggesting involvement in milk synthesis and production. Furthermore, the expression of two of these genes (NUP35 and OXCT1) was enriched in immune tissues implying a potentially pleiotropic effect or likely role in milk production during udder infection, which needs to be further elucidated in future studies. In conclusion, the absence of genetic antagonism between milk yield and mastitis resistance suggests that simultaneous genetic improvement of both traits be achievable.

Genome-wide association study of milk and reproductive traits in dual-purpose Xinjiang Brown cattle

Background

Dual-purpose cattle are more adaptive to environmental challenges than single-purpose dairy or beef cattle. Balance among milk, reproductive, and mastitis resistance traits in breeding programs is therefore more critical for dual-purpose cattle to increase net income and maintain well-being. With dual-purpose Xinjiang Brown cattle adapted to the Xinjiang Region in northwestern China, we conducted genome-wide association studies (GWAS) to dissect the genetic architecture related to milk, reproductive, and mastitis resistance traits. Phenotypic data were collected for 2410 individuals measured during 1995–2017. By adding another 445 ancestors, a total of 2855 related individuals were used to derive estimated breeding values for all individuals, including the 2410 individuals with phenotypes. Among phenotyped individuals, we genotyped 403 cows with the Illumina 150 K Bovine BeadChip.

Results

GWAS were conducted with the FarmCPU (Fixed and random model circulating probability unification) method. We identified 12 markers significantly associated with six of the 10 traits under the threshold of 5% after a Bonferroni multiple test correction. Seven of these SNPs were in QTL regions previously identified to be associated with related traits. One identified SNP, BovineHD1600006691, was significantly associated with both age at first service and age at first calving. This SNP directly overlapped a QTL previously reported to be associated with calving ease. Within 160 Kb upstream and downstream of each significant SNP identified, we speculated candidate genes based on functionality. Four of the SNPs were located within four candidate genes, including CDH2, which is linked to milk fat percentage, and GABRG2, which is associated with milk protein yield.

Conclusions

These findings are beneficial not only for breeding through marker-assisted selection, but also for genome editing underlying the related traits to enhance the overall performance of dual-purpose cattle.

Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows

Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.

A Population Genomics Analysis of the Native Irish Galway Sheep Breed

The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.

Genome-wide scan reveals genetic divergence and diverse adaptive selection in Chinese local cattle

Background

Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs.

Results

We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions.

Conclusions

Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5822-y) contains supplementary material, which is available to authorized users.

Genome-Wide SNP Discovery in Indigenous Cattle Breeds of South Africa

Single nucleotide polymorphism arrays have created new possibilities for performing genome-wide studies to detect genomic regions harboring sequence variants that affect complex traits. However, the majority of validated SNPs for which allele frequencies have been estimated are limited primarily to European breeds. The objective of this study was to perform SNP discovery in three South African indigenous breeds (Afrikaner, Drakensberger, and Nguni) using whole genome sequencing. DNA was extracted from blood and hair samples, quantified and prepared at 50 ng/μl concentration for sequencing at the Agricultural Research Council Biotechnology Platform using an Illumina HiSeq 2500. The fastq files were used to call the variants using the Genome Analysis Tool Kit. A total of 1,678,360 were identified as novel using Run 6 of 1000 Bull Genomes Project. Annotation of the identified variants classified them into functional categories. Within the coding regions, about 30% of the SNPs were non-synonymous substitutions that encode for alternate amino acids. The study of distribution of SNP across the genome identified regions showing notable differences in the densities of SNPs among the breeds and highlighted many regions of functional significance. Gene ontology terms identified genes such as MLANA, SYT10, and CDC42EP5 that have been associated with coat color in mouse, and ADAMS3, DNAJC3, and PAG5 genes have been associated with fertility in cattle. Further analysis of the variants detected 688 candidate selective sweeps (ZHp Z-scores ≤ -4) across all three breeds, of which 223 regions were assigned as being putative selective sweeps (ZHp scores ≤-5). We also identified 96 regions with extremely low ZHp Z-scores (≤-6) in Afrikaner and Nguni. Genes such as KIT and MITF that have been associated with skin pigmentation in cattle and CACNA1C, which has been associated with biopolar disorder in human, were identified in these regions. This study provides the first analysis of sequence data to discover SNPs in indigenous South African cattle breeds. The information will play an important role in our efforts to understand the genetic history of our cattle and in designing appropriate breed improvement programmes.

Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions

Northern Fennoscandia and the Sakha Republic in the Russian Federation represent the northernmost regions on Earth where cattle farming has been traditionally practiced. In this study, we performed whole-genome sequencing to genetically characterize three rare native breeds Eastern Finncattle, Western Finncattle and Yakutian cattle adapted to these northern Eurasian regions. We examined the demographic history, genetic diversity and unfolded loci under natural or artificial selection. On average, we achieved 13.01-fold genome coverage after mapping the sequencing reads on the bovine reference genome (UMD 3.1) and detected a total of 17.45 million single nucleotide polymorphisms (SNPs) and 1.95 million insertions-deletions (indels). We observed that the ancestral species (Bos primigenius) of Eurasian taurine cattle experienced two notable prehistorical declines in effective population size associated with dramatic climate changes. The modern Yakutian cattle exhibited a higher level of within-population variation in terms of number of SNPs and nucleotide diversity than the contemporary European taurine breeds. This result is in contrast to the results of marker-based cattle breed diversity studies, indicating assortment bias in previous analyses. Our results suggest that the effective population size of the ancestral Asiatic taurine cattle may have been higher than that of the European cattle. Alternatively, our findings could indicate the hybrid origins of the Yakutian cattle ancestries and possibly the lack of intensive artificial selection. We identified a number of genomic regions under selection that may have contributed to the adaptation to the northern and subarctic environments, including genes involved in disease resistance, sensory perception, cold adaptation and growth. By characterizing the native breeds, we were able to obtain new information on cattle genomes and on the value of the adapted breeds for the conservation of cattle genetic resources.

Genetic mechanisms underlying spermatic and testicular traits within and among cattle breeds: systematic review and prioritization of GWAS results

Reduced bull fertility imposes economic losses in bovine herds. Specifically, testicular and spermatic traits are important indicators of reproductive efficiency. Several genome-wide association studies (GWAS) have identified genomic regions associated with these fertility traits. The aims of this study were as follows: 1) to perform a systematic review of GWAS results for spermatic and testicular traits in cattle and 2) to identify key functional candidate genes for these traits. The identification of functional candidate genes was performed using a systems biology approach, where genes shared between traits and studies were evaluated by a guilt by association gene prioritization (GUILDify and ToppGene software) in order to identify the best functional candidates. These candidate genes were integrated and analyzed in order to identify overlapping patterns among traits and breeds. Results showed that GWAS for testicular-related traits have been developed for beef breeds only, whereas the majority of GWAS for spermatic-related traits were conducted using dairy breeds. When comparing traits measured within the same study, the highest number of genes shared between different traits was observed, indicating a high impact of the population genetic structure and environmental effects. Several chromosomal regions were enriched for functional candidate genes associated with fertility traits. Moreover, multiple functional candidate genes were enriched for markers in a species-specific basis, taurine (Bos taurus) or indicine (Bos indicus). For the different candidate regions identified in the GWAS in the literature, functional candidate genes were detected as follows: B. Taurus chromosome X (BTX) (TEX11, IRAK, CDK16, ATP7A, ATRX, HDAC6, FMR1, L1CAM, MECP2, etc.), BTA17 (TRPV4 and DYNLL1), and BTA14 (MOS, FABP5, ZFPM2). These genes are responsible for regulating important metabolic pathways or biological processes associated with fertility, such as progression of spermatogenesis, control of ciliary activity, development of Sertoli cells, DNA integrity in spermatozoa, and homeostasis of testicular cells. This study represents the first systematic review on male fertility traits in cattle using a system biology approach to identify key candidate genes for these traits.

Combining multi-OMICs information to identify key-regulator genes for pleiotropic effect on fertility and production traits in beef cattle

The identification of biological processes related to the regulation of complex traits is a difficult task. Commonly, complex traits are regulated through a multitude of genes contributing each to a small part of the total genetic variance. Additionally, some loci can simultaneously regulate several complex traits, a phenomenon defined as pleiotropy. The lack of understanding on the biological processes responsible for the regulation of these traits results in the decrease of selection efficiency and the selection of undesirable hitchhiking effects. The identification of pleiotropic key-regulator genes can assist in developing important tools for investigating biological processes underlying complex traits. A multi-breed and multi-OMICs approach was applied to study the pleiotropic effects of key-regulator genes using three independent beef cattle populations evaluated for fertility traits. A pleiotropic map for 32 traits related to growth, feed efficiency, carcass and meat quality, and reproduction was used to identify genes shared among the different populations and breeds in pleiotropic regions. Furthermore, data-mining analyses were performed using the Cattle QTL database (CattleQTLdb) to identify the QTL category annotated in the regions around the genes shared among breeds. This approach allowed the identification of a main gene network (composed of 38 genes) shared among breeds. This gene network was significantly associated with thyroid activity, among other biological processes, and displayed a high regulatory potential. In addition, it was possible to identify genes with pleiotropic effects related to crucial biological processes that regulate economically relevant traits associated with fertility, production and health, such as MYC, PPARG, GSK3B, TG and IYD genes. These genes will be further investigated to better understand the biological processes involved in the expression of complex traits and assist in the identification of functional variants associated with undesirable phenotypes, such as decreased fertility, poor feed efficiency and negative energetic balance.

Genetic diversity and selection signatures of the beef 'Charolais de Cuba' breed

In this study, we used BovineSNP50 Genotyping BeadChip data to estimate the structure, putative ancestral origin as well as to identify regions with selective sweeps that may have had an important role in the adaptation to tropical conditions of the 'Charolais de Cuba' (CHCU) breed. According to a principal component analysis, CHCU samples cluster together with taurine breeds with an estimated 93% of taurus ancestral alleles. Despite the short period since importation, we detected differentiation (Fst = 0.049) between the French Charolaise (CHA) and CHCU. However, CHA breed was the closest breed to CHCU followed by other hybrids breed with a clear CHA origin. Linkage disequilibrium (r²) decay tends to be lower in CHCU compared to CHA probably due to a less intense artificial selection programs of CHCU. Signals of recent adaptation to tropical conditions between CHCU and CHA were identified. Genes mapping within those regions reflect different functions related to immunity, metabolic changes and heat tolerance (CHCU) and muscle development and meat quality (CHA) that may have had an important role in the phenotypic differentiation of these breeds. Further studies will expand our knowledge on the molecular basis of adaptation of cattle to tropical conditions and molecular process associated with meat quality traits.

Genome-Wide Signatures of Selection Reveal Genes Associated With Performance in American Quarter Horse Subpopulations

Selective breeding for athletic performance in various disciplines has resulted in population stratification within the American Quarter Horse (QH) breed. The goals of this study were to utilize high density genotype data to: (1) identify genomic regions undergoing positive selection within and among QH subpopulations; (2) investigate haplotype structure within each QH subpopulation; and (3) identify candidate genes within genomic regions of interest (ROI), as well as biological pathways, predicted to play a role in elite performance in each group. For that, 65K SNP genotyping data on 143 elite individuals from 6 QH subpopulations (cutting, halter, racing, reining, western pleasure, and working cow) were imputed to 2M SNPs. Signatures of selection were identified using F_ST-based (d_i ) and haplotype-based (hapFLK) analyses, accompanied by identification of local haplotype structure and sharing within subpopulations (hapQTL). Regions undergoing positive selection were identified on all 31 autosomes, and ROI on 2 chromosomes were identified by all 3 methods combined. Genes within each ROI were retrieved and used to identify pathways and genes that might contribute to performance in each subpopulation. These included, among others, candidate genes associated with skeletal muscle development, metabolism, and central nervous system development. This work improves our understanding of equine breed development, and provides breeders with a better understanding of how selective breeding impacts the performance of QH populations.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing

BACKGROUND

Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids.

RESULTS

Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts.

CONCLUSIONS

This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

Species composition and environmental adaptation of indigenous Chinese cattle

Indigenous Chinese cattle combine taurine and indicine origins and occupy a broad range of different environments. By 50 K SNP genotyping we found a discontinuous distribution of taurine and indicine cattle ancestries with extremes of less than 10% indicine cattle in the north and more than 90% in the far south and southwest China. Model-based clustering and f4-statistics indicate introgression of both banteng and gayal into southern Chinese cattle while the sporadic yak influence in cattle in or near Tibetan area validate earlier findings of mitochondrial DNA analysis. Geographic patterns of taurine and indicine mitochondrial and Y-chromosomal DNA diversity largely agree with the autosomal cline. The geographic distribution of the genomic admixture of different bovine species is proposed to be the combined effect of prehistoric immigrations, gene flow, major rivers acting as genetic barriers, local breeding objectives and environmental adaptation. Whole-genome scan for genetic differentiation and association analyses with both environmental and morphological covariables are remarkably consistent with previous studies and identify a number of genes implicated in adaptation, which include TNFRSF19, RFX4, SP4 and several coat color genes. We propose indigenous Chinese cattle as a unique and informative resource for gene-level studies of climate adaptation in mammals.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations

BACKGROUND

Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF).

RESULTS

The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26.

CONCLUSIONS

The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.

High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds

BACKGROUND

With the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals.

RESULTS

Based on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included.

CONCLUSIONS

Here, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.

Genomic Comparison of Indigenous African and Northern European Chickens Reveals Putative Mechanisms of Stress Tolerance Related to Environmental Selection Pressure

Global climate change is increasing the magnitude of environmental stressors, such as temperature, pathogens, and drought, that limit the survivability and sustainability of livestock production. Poultry production and its expansion is dependent upon robust animals that are able to cope with stressors in multiple environments. Understanding the genetic strategies that indigenous, noncommercial breeds have evolved to survive in their environment could help to elucidate molecular mechanisms underlying biological traits of environmental adaptation. We examined poultry from diverse breeds and climates of Africa and Northern Europe for selection signatures that have allowed them to adapt to their indigenous environments. Selection signatures were studied using a combination of population genomic methods that employed FST , integrated haplotype score (iHS), and runs of homozygosity (ROH) procedures. All the analyses indicated differences in environment as a driver of selective pressure in both groups of populations. The analyses revealed unique differences in the genomic regions under selection pressure from the environment for each population. The African chickens showed stronger selection toward stress signaling and angiogenesis, while the Northern European chickens showed more selection pressure toward processes related to energy homeostasis. The results suggest that chromosomes 2 and 27 are the most diverged between populations and the most selected upon within the African (chromosome 27) and Northern European (chromosome 2) birds. Examination of the divergent populations has provided new insight into genes under possible selection related to tolerance of a population's indigenous environment that may be baselines for examining the genomic contribution to tolerance adaptions.

Use of high-density SNP data to identify patterns of diversity and signatures of selection in broiler chickens

The development of broiler chickens over the last 70 years has been accompanied by large phenotypic changes, so that the resulting genomic signatures of selection should be detectable by current statistical techniques with sufficiently dense genetic markers. Using two approaches, this study analysed high-density SNP data from a broiler chicken line to detect low-diversity genomic regions characteristic of past selection. Seven regions with zero diversity were identified across the genome. Most of these were very small and did not contain many genes. In addition, fifteen regions were identified with diversity increasing asymptotically from a low level. These regions were larger and thus generally included more genes. Several candidate genes for broiler traits were found within these 'regression regions', including IGF1, GPD2 and MTNR1AI. The results suggest that the identification of zero-diversity regions is too restrictive for characterizing regions under selection, but that regions showing patterns of diversity along the chromosome that are consistent with selective sweeps contain a number of genes that are functional candidates for involvement in broiler development. Many regions identified in this study overlap or are close to regions identified in layer chicken populations, possibly due to their shared precommercialization history or to shared selection pressures between broilers and layers.

Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers

Background

The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation.

Results

The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean H_O and H_E was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures F_EH, F_VR, F_LEUT, F_ROH and F_PED was 0.129, −0.012, −0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis.

Conclusions

In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds’ development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed F_ST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3610-0) contains supplementary material, which is available to authorized users.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.