Analysis of genome-wide DNA arrays reveals the genomic population structure and diversity in autochthonous Greek goat breeds

Goats play an important role in the livestock sector in Greece. The national herd consists mainly of two indigenous breeds, the Eghoria and Skopelos. Here, we report the population structure and genomic profiles of these two native goat breeds using Illumina's Goat SNP50 BeadChip. Moreover, we present a panel of candidate markers acquired using different genetic models for breed discrimination. Quality control on the initial dataset resulted in 48,841 SNPs kept for downstream analysis. Principal component and admixture analyses were applied to assess population structure. The rate of inbreeding within breed was evaluated based on the distribution of runs of homozygosity in the genome and respective coefficients, the genomic relationship matrix, the patterns of linkage disequilibrium, and the historic effective population size. Results showed that both breeds exhibit high levels of genetic diversity. Level of inbreeding between the two breeds estimated by the Wright's fixation index FST was low (Fst = 0.04362), indicating the existence of a weak genetic differentiation between them. In addition, grouping of farms according to their geographical locations was observed. This study presents for the first time a genome-based analysis on the genetic structure of the two indigenous Greek goat breeds and identifies markers that can be potentially exploited in future selective breeding programs for traceability purposes, targeted genetic improvement schemes and conservation strategies

Validation of the 50k Illumina goat SNP chip in the South African Angora goat

Tools for the genomic evaluation of goats have generally lagged behind those for other species. However, the recent availability of the goat SNP50 consortium bead chip has marked a positive change for this small ruminant species. Polymorphic loci can differ greatly between breeds of the same species. Exclusion of fibre-producing breeds, such as the Angora goat, during the development of this genotyping array necessitates the validation of SNPs included on the chip to allow for genomic applications that would accelerate genetic progress in mohair yield and quality. Forty eight unrelated Angora goats, displaying phenotypic variation in two important price-determining traits, namely fibre diameter and fleece weight, were genotyped with the goat SNP50 consortium bead chip. Results revealed that 46 983 SNP (88.1%) of the 53 347 called SNPs were polymorphic (MAF>0.05). After quality control, 3 960 SNP were filtered from further analysis for violating Hardy-Weinberg and call-rate parameter thresholds, leaving 43 759 (82%) of the 53 347 SNPs to be validated for downstream analysis. Observed and expected heterozygosity values of 0.365 and 0.370, respectively, were obtained for polymorphic SNPs. A total of 30 357 SNPs in linkage disequilibrium (LD) were removed to obtain a set of independent markers, resulting in a final SNP density of 1 SNP/ ~226 kb. Results indicate that the goat SNP50 bead chip was informative in the Angora goats that were studied, and should be useful in examining the underlying genetic variationKeywords: Call rate, Hardy-Weinberg equilibrium, linkage disequilibrium, minor allele frequenc

Data from: Landscape genomics and pathway analysis to understand genetic adaptation of South African indigenous goat populations

Extensively raised livestock populations in most smallholder farming communities are exposed to harsh and heterogeneous climates and disease pathogen that they adapt to in order to survive. Majority of these livestock species including goats are of non-descript breeds and their response to natural selection presented by heterogeneous environments is still unresolved. This study investigated genetic diversity and its association with environmental and geographic conditions in 194 South African indigenous goats from different geographic using Illumina goat SNP50K genotypes. The Tankwa goats formed a homogeneous genetic cluster restricted to the Northern Cape province. Overall, the Boer, Kalahari Red, and Savanna showed a wide geographic spread of shared genetic components whereas the village ecotypes revealed a longitudinal distribution. The relative importance of environmental factors on genetic variation of goat populations was assessed using Redundancy Analysis (RDA). Variance partitioning using Model I with both climatic and geographic variables explained 22% of the total variation while model II with climatic variables accounted for 17%. Model III, which accounted for geographic variables, explained 1% of the total variation. The first axis (Model I) of the RDA analysis, revealed 329 outlier SNPs. Landscape genomic approaches of Spatial Analysis Method identified a total of 843 (1.75%) SNPs, while Latent Factor Mixed Models identified 714 (1.48%) SNPs significantly associated with environmental variables. Significant markers were within genes involved in biological functions potentially important for environmental adaptation. Overall, the study suggested environmental factors to have some effect in shaping the genetic variation of South African indigenous goat populations

Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment

Goats and sheep are versatile domesticates that have been integrated into diverse environments and production systems. Natural and artificial selection have shaped the variation in the two species, but natural selection has played the major role among indigenous flocks. To investigate signals of natural selection, we analyzed genotype data generated using the caprine and ovine 50K SNP BeadChips from Barki goats and sheep that are indigenous to a hot arid environment in Egypt|[rsquo]|s Coastal Zone of the Western Desert. We identify several candidate regions under selection that spanned 119 genes. A majority of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. In particular, selection signatures spanning several genes that directly or indirectly influenced traits for adaptation to hot arid environments, such as thermo-tolerance (melanogenesis) (FGF2, GNAI3, PLCB1), body size and development (BMP2, BMP4, GJA3, GJB2), energy and digestive metabolism (MYH, TRHDE, ALDH1A3), and nervous and autoimmune response (GRIA1, IL2, IL7, IL21, IL1R1) were identified. We also identified eight common candidate genes under selection in the two species and a shared selection signature that spanned a conserved syntenic segment to bovine chromosome 12 on caprine and ovine chromosomes 12 and 10, respectively, providing, most likely, the evidence for selection in a common environment in two different but closely related species. Our study highlights the importance of indigenous livestock as model organisms for investigating selection sweeps and genome-wide association mapping