Whole genome population genetics analysis of Sudanese goats identifies regions harboring genes associated with major traits

Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats

BACKGROUND

Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated.

RESULTS

Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh.

CONCLUSIONS

The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

Genetic Diversity, Population Structure and Selection Signature in Begait Goats Revealed by Whole-Genome Sequencing

Goats belong to a group of animals called small ruminants and are critical sources of livelihood for rural people. Genomic sequencing can provide information ranging from basic knowledge about goat diversity and evolutionary processes that shape genomes to functional information about genes/genomic regions. In this study, we exploited a whole-genome sequencing data set to analyze the genetic diversity, population structure and selection signatures of 44 individuals belonging to 5 Ethiopian goat populations: 12 Aberegalle (AB), 5 Afar (AF), 11 Begait (BG), 12 Central highlands (CH) and 5 Meafure (MR) goats. Our results revealed the highest genetic diversity in the BG goat population compared to the other goat populations. The pairwise genetic differentiation (FST) among the populations varied and ranged from 0.011 to 0.182, with the closest pairwise value (0.003) observed between the AB and CH goats and a distant correlation (FST = 0.182) between the BG and AB goats, indicating low to moderate genetic differentiation. Phylogenetic tree, ADMIXTURE and principal component analyses revealed a classification of the five Ethiopian goat breeds in accordance with their geographic distribution. We also found three top genomic regions that were detected under selection on chromosomes 2, 5 and 13. Moreover, this study identified different candidate genes related to milk characteristics (GLYCAM1 and SRC), carcass (ZNF385B, BMP-7, PDE1B, PPP1R1A, FTO and MYOT) and adaptive and immune response genes (MAPK13, MAPK14, SCN7A, IL12A, EST1 DEFB116 and DEFB119). In conclusion, this information could be helpful for understanding the genetic diversity and population structure and selection scanning of these important indigenous goats for future genetic improvement and/or as an intervention mechanism.

Diverse WGBS profiles of longissimus dorsi muscle in Hainan black goats and hybrid goats

BACKGROUND

Goat products have played a crucial role in meeting the dietary demands of people since the Neolithic era, giving rise to a multitude of goat breeds globally with varying characteristics and meat qualities. The primary objective of this study is to pinpoint the pivotal genes and their functions responsible for regulating muscle fiber growth in the longissimus dorsi muscle (LDM) through DNA methylation modifications in Hainan black goats and hybrid goats.

METHODS

Whole-genome bisulfite sequencing (WGBS) was employed to scrutinize the impact of methylation on LDM growth. This was accomplished by comparing methylation differences, gene expression, and their associations with growth-related traits.

RESULTS

In this study, we identified a total of 3,269 genes from differentially methylated regions (DMR), and detected 189 differentially expressed genes (DEGs) through RNA-seq analysis. Hypo DMR genes were primarily enriched in KEGG terms associated with muscle development, such as MAPK and PI3K-Akt signaling pathways. We selected 11 hub genes from the network that intersected the gene sets within DMR and DEGs, and nine genes exhibited significant correlation with one or more of the three LDM growth traits, namely area, height, and weight of loin eye muscle. Particularly, PRKG1 demonstrated a negative correlation with all three traits. The top five most crucial genes played vital roles in muscle fiber growth: FOXO3 safeguarded the myofiber's immune environment, FOXO6 was involved in myotube development and differentiation, and PRKG1 facilitated vasodilatation to release more glucose. This, in turn, accelerated the transfer of glucose from blood vessels to myofibers, regulated by ADCY5 and AKT2, ultimately ensuring glycogen storage and energy provision in muscle fibers.

CONCLUSION

This study delved into the diverse methylation modifications affecting critical genes, which collectively contribute to the maintenance of glycogen storage around myofibers, ultimately supporting muscle fiber growth.

Population structure of wild soybean (Glycine soja) based on SLAF-seq have implications for its conservation

Background

Glycine soja Sieb. & Zucc. is the wild ancestor from which the important crop plant soybean was bred. G. soja provides important germplasm resources for the breeding and improvement of cultivated soybean crops, however the species is threatened by habitat loss and fragmentation, and is experiencing population declines across its natural range. Understanding the patterns of genetic diversity in G. soja populations can help to inform conservation practices.

Methods

In this study, we analyzed the genetic diversity and differentiation of G. soja at different sites and investigated the gene flow within the species. We obtained 147 G. soja accessions collected from 16 locations across the natural range of the species from China, Korea and Japan. Samples were analyzed using SLAF-seq (Specific-Locus Amplified Fragment Sequencing).

Results

We obtained a total of 56,489 highly consistent SNPs. Our results suggested that G. soja harbors relatively high diversity and that populations of this species are highly differentiated. The populations harboring high genetic diversity, especially KR, should be considered first when devising conservation plans for the protection of G. soja, and in situ protection should be adopted in KR. G. soja populations from the Yangtze River, the Korean peninsula and northeastern China have a close relationship, although these areas are geographically disconnected. Other populations from north China clustered together. Analysis of gene flow suggested that historical migrations of G. soja may have occurred from the south northwards across the East-Asia land-bridge, but not across north China. All G. soja populations could be divided into one of two lineages, and these two lineages should be treated separately when formulating protection policies.

High-density genotyping reveals candidate genomic regions for chicken body size in breeds of Asian origin

Body size is one of the main selection indices in chicken breeding. Although often investigated, knowledge of the underlying genetic mechanisms is incomplete. The aim of the current study was to identify genomic regions associated with body size differences between Asian Game and Asian Bantam type chickens. In this study, 94 and 107 chickens from 4 Asian Game and 5 Asian Bantam type breeds, respectively, were genotyped using the chicken 580K single nucleotide polymorphism (SNP) array. A genome-wide association study (GWAS) and principal component analyses (PCA) were performed to identify genomic regions associated with body size related-traits such as wing length, shank length, shank thickness, keel length, and body weight. Hierarchical clustering of genotype data showed a clear genetic difference between the investigated Asian Game and Asian Bantam chicken types. GWAS identified 16 genomic regions associated with wing length (2, FDR ≤ 0.018), shank thickness (6, FDR ≤ 0.008), keel length (5, FDR ≤ 0.023), and body weight (3, FDR ≤ 0.041). PCA showed that the first principal component (PC1) separated the 2 chicken types and significantly correlated with the measured body size related-traits (P ≤ 2.24e-40). SNPs contributing significantly to PC1 were subjected to a more detailed investigation. This analysis identified 11 regions potentially associated with differences in body size related-traits. A region on chromosome 4 (GGA4) (17.3-21.3 Mb) was detected in both analyses GWAS and PCA. This region harbors 60 genes. Among them are myotubularin 1 (MTM1) and secreted frizzled-related protein 2 (SFPR2) which can be considered as potential candidate genes for body size related-traits. Our results clearly show that the investigated Asian Game type chicken breeds are genetically different from the Asian Bantam breeds. A region on GGA4 between 17.3 and 21.3 Mb was identified which contributes to the phenotypic difference, though further validation of candidate genes is necessary.

Genome-Wide Population Structure and Selection Signatures of Yunling Goat Based on RAD-seq

Simple Summary

Goats are important domestic animals that provide meat, milk, fur, and other products for humans. The demand for these products has increased in recent years. Disease resistance among goat breeds is different, but the genetic basis of the differences in resistance to diseases is still unclear and needs to be further studied. In this study, many genes and pathways related to immunity and diseases were identified to be under positive selection between Yunling and Nubian goats using RAD-seq technology. This study on the selection signatures of Yunling goats provides the scientific basis and technical support for the breeding of domestic goats for disease resistance, which has important social and economic significance.

Abstract

Animal diseases impose a huge burden on the countries where diseases are endemic. Conventional control strategies of vaccines and veterinary drugs are to control diseases from a pharmaceutical perspective. Another alternative approach is using pre-existing genetic disease resistance or tolerance. We know that the Yunling goat is an excellent local breed from Yunnan, southwestern China, which has characteristics of strong disease resistance and remarkable adaptability. However, genetic information about the selection signatures of Yunling goats is limited. We reasoned that the genes underlying the observed difference in disease resistance might be identified by investigating selection signatures between two different goat breeds. Herein, we selected the Nubian goat as the reference group to perform the population structure and selection signature analysis by using RAD-seq technology. The results showed that two goat breeds were divided into two clusters, but there also existed gene flow. We used Fst (F-statistics) and π (pi/θπ) methods to carry out selection signature analysis. Eight selected regions and 91 candidate genes were identified, in which some genes such as DOK2, TIMM17A, MAVS, and DOCK8 related to disease and immunity and some genes such as SPEFI, CDC25B, and MIR103 were associated with reproduction. Four GO (Gene Ontology) terms (GO:0010591, GO:001601, GO:0038023, and GO:0017166) were associated with cell migration, signal transduction, and immune responses. The KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways were mainly associated with immune responses, inflammatory responses, and stress reactions. This study preliminarily revealed the genetic basis of strong disease resistance and adaptability of Yunling goats. It provides a theoretical basis for the subsequent genetic breeding of disease resistance of goats.

Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection

Background

The domestication of goat (Capra hircus) started 11,000 years ago in the fertile crescent. Breed formation in the nineteenth century, establishment of herd books, and selection for specific traits resulted in 10 modern goat breeds in Switzerland. We analyzed whole-genome sequencing (WGS) data from 217 modern goats and nine wild Bezoar goats (Capra aegagrus). After quality control, 27,728,288 biallelic single nucleotide variants (SNVs) were used for the identification of runs of homozygosity (ROH) and the detection of ROH islands.

Results

Across the 226 caprine genomes from 11 populations, we detected 344 ROH islands that harbor 1220 annotated genes. We compared the ROH islands between the modern breeds and the Bezoar goats. As a proof of principle, we confirmed a signature of selection, which contains the ASIP gene that controls several breed-specific coat color patterns. In two other ROH islands, we identified two missense variants, STC1:p.Lys139Arg and TSHR:p.Ala239Thr, which might represent causative functional variants for domestication signatures.

Conclusions

We have shown that the information from ROH islands using WGS data is suitable for the analysis of signatures of selection and allowed the detection of protein coding variants that may have conferred beneficial phenotypes during goat domestication. We hypothesize that the TSHR:p.Ala239Thr variant may have played a role in changing the seasonality of reproduction in modern domesticated goats. The exact functional significance of the STC1:p.Lys139Arg variant remains unclear and requires further investigation. Nonetheless, STC1 might represent a new domestication gene affecting relevant traits such as body size and/or milk yield in goats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00695-w.

SNP-Based Genotyping Provides Insight Into the West Asian Origin of Russian Local Goats

Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (<2 Mb) were predominant in all breeds, while the longest ROH class (>16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.

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

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

Capture Sequencing to Explore and Map Rare Casein Variants in Goats

Genetic variations in the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3 have obtained substantial attention since they affect the milk protein yield, milk composition, cheese processing properties, and digestibility as well as tolerance in human nutrition. Furthermore, milk protein variants are used for breed characterization, biodiversity, and phylogenetic studies. The current study aimed at the identification of casein protein variants in five domestic goat breeds from Sudan (Nubian, Desert, Nilotic, Taggar, and Saanen) and three wild goat species [Capra aegagrus aegagrus (Bezoar ibex), Capra nubiana (Nubian ibex), and Capra ibex (Alpine ibex)]. High-density capture sequencing of 33 goats identified in total 22 non-synonymous and 13 synonymous single nucleotide polymorphisms (SNPs), of which nine non-synonymous and seven synonymous SNPs are new. In the CSN1S1 gene, the new non-synonymous SNP ss7213522403 segregated in Alpine ibex. In the CSN2 gene, the new non-synonymous SNPs ss7213522526, ss7213522558, and ss7213522487 were found exclusively in Nubian and Alpine ibex. In the CSN1S2 gene, the new non-synonymous SNPs ss7213522477, ss7213522549, and ss7213522575 were found in Nubian ibex only. In the CSN3 gene, the non-synonymous SNPs ss7213522604 and ss7213522610 were found in Alpine ibex. The identified DNA sequence variants led to the detection of nine new casein protein variants. New variants were detected for alpha S1 casein in Saanen goats (CSN1S1^∗C1), Bezoar ibex (CSN1S1^∗J), and Alpine ibex (CSN1S1^∗K), for beta and kappa caseins in Alpine ibex (CSN2^∗F and CSN3^∗X), and for alpha S2 casein in all domesticated and wild goats (CSN1S2^∗H), in Nubian and Desert goats (CSN1S2^∗I), or in Nubian ibex only (CSN1S2^∗J and CSN1S2^∗K). The results show that most novel SNPs and protein variants occur in the critically endangered Nubian ibex. This highlights the importance of the preservation of this endangered breed. Furthermore, we suggest validating and further characterizing the new casein protein variants.

Whole-genome sequencing to identify candidate genes for litter size and to uncover the variant function in goats (Capra hircus)

To select candidate genes for goat prolificacy, we managed six multi- and six single-kid female goats at the same feeding level and in the same management mode over a 4-year period. These goats showed stable differences in litter size over five continuous parturition records. Whole-genome re-sequencing was used in all 12 to select candidate genes, namely, AURKA, ENDOG, SOX2, RORA, GJA10, RXFP2, CDC25C, and NANOS3, by the strength of their differentiation signals. Most of the selected genes were enriched in the coiled coil process and ovarian development, which suggests that the coiled coil process has a potential regulatory effect on fecundity. Detection of the distribution of variants and association analyses with litter size in 400 goats showed that NANOS3 exon mutations may lead to a transformation of the protein structure. The variation in CDC25C, ENDOG, and NANOS3 showed a significant association with litter size. These results can contribute to the improvement of reproduction traits in the artificial breeding of goats.

Skin transcriptome reveals the periodic changes in genes underlying cashmere (ground hair) follicle transition in cashmere goats

Background

Cashmere goats make an outstanding contribution to the livestock textile industry and their cashmere is famous for its slenderness and softness and has been extensively studied. However, there are few reports on the molecular regulatory mechanisms of the secondary hair follicle growth cycle in cashmere goats. In order to explore the regular transition through the follicle cycle and the role of key genes in this cycle, we used a transcriptome sequencing technique to sequence the skin of Inner Mongolian cashmere goats during different months. We analyzed the variation and difference in genes throughout the whole hair follicle cycle. We then verified the regulatory mechanism of the cashmere goat secondary hair follicle growth cycle using fluorescence quantitative PCR.

Results

The growth cycle of cashmere hair could be divided into three distinct periods: a growth period (March–September), a regression period (September–December), and a resting period (December–March). The results of differential gene analyses showed that March was the most significant month. Cluster analysis of gene expression throughout the whole growth cycle further supported the key nodes of the three periods of cashmere growth, and the differential gene expression of keratin corresponding to the ground haircashmere growth cycle further supported the results from tissue slices. Quantitative fluorescence analysis showed that KAP3–1, KRTAP 8–1, and KRTAP 24–1 genes had close positive correlation with the cashmere growth cycle, and their regulation was consistent with the growth cycle of cashmere.

Conclusion

The growth cycle of cashmere cashmere could be divided into three distinct periods: a growth period (March–September), a regression period (September–December) and a resting period (December–March). March was considered to be the beginning of the cycle. KAP and KRTAP showed close positive correlation with the growth cycle of secondary hair follicle cashmere growth, and their regulation was consistent with the cashmere growth cycle. But hair follicle development-related genes are expressed earlier than cashmere growth, indicating that cycle regulation could alter the temporal growth of cashmere. This study laid a theoretical foundation for the study of the cashmere development cycle and provided evidence for key genes during transition through the cashmere cycle. Our study provides a theoretical basis for cashmere goat breeding.

Population structure of indigenous southern African goats based on the Illumina Goat50K SNP panel

In this study, the genetic structure of indigenous Tswana and Swazi goats using the Illumina Goat50K SNP array was investigated. Two South African commercial goat breeds were included to investigate admixture with the indigenous populations in southern Africa. A total of 144 DNA samples including Boer goats (n = 24), Kalahari Red (n = 24), Swazi (n = 48), and Tswana goats (n = 48) were genotyped. Statistical analysis was performed using PLINK version 1.07. Genetic diversity, measured as expected heterozygosity, was estimated at 0.390, 0.398, 0.413, and 0.387 for Boer, Kalahari Red, Tswana, and Swazi goats, respectively. The individual inbreeding coefficient varied from 0.019 ± 0.05 to 0.011 ± 0.06 for the Tswana and Swazi goats, respectively. The Principal component analysis clustered the populations according to geographical origin and breed type. Linkage disequilibrium (LD) for shorter intervals (0-10 kb) ranged from 0.44 to 0.56 and commercial breeds had higher values. Effective population sizes decreased with generations and at the 13th generation ranged between 87 for Boer to 266 for Tswana goats. The Tswana population exhibited the highest level of genetic variation and effective population size, which holds potential for improved production in marginal regions. A national strategy is required to maintain genetic diversity in communal goat production systems through well-structured breeding and conservation programs.

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 F_ST 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.

Whey protein polymorphisms in Sudanese goat breeds

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

Genome-wide diversity and demographic dynamics of Cameroon goats and their divergence from east African, north African, and Asian conspecifics

Indigenous goats make significant contributions to Cameroon’s national and local economy, but little effort has been devoted to identifying the populations. Here, we assessed the genetic diversity and demographic dynamics of Cameroon goat populations using mitochondrial DNA (two populations) and autosomal markers (four populations) generated with the Caprine 50K SNP chip. To infer genetic relationships at continental and global level, genotype data on six goat populations from Ethiopia and one population each from Egypt, Morocco, Iran, and China were included in the analysis. The mtDNA analysis revealed 83 haplotypes, all belonging to haplogroup A, in Cameroon goats. Four haplotypes were shared between goats found in Cameroon, Mozambique, Namibia, Zimbabwe, Kenya, and Ethiopia. Analysis of autosomal SNPs in Cameroon goats revealed the lowest H_O (0.335±0.13) and H_E (0.352±0.15) in the North-west Highland and Central Highland populations, respectively. Overall, the highest H_O (0.401±0.12) and H_E (0.422±0.12) were found for Barki and Iranian goats, respectively. Barki goats had the highest average MAF, while Central Highland Cameroon goats had the lowest. Overall, Cameroon goats demonstrated high F_IS. AMOVA revealed that 13.29% of the variation was explained by genetic differences between the six population groups. Low average F_ST (0.01) suggests intermixing among Cameroon goats. All measures indicated that Cameroon goats are closer to Moroccan goats than to other goat populations. PCA and STRUCTURE analyses poorly differentiated the Cameroon goats, as did genetic distance, Neighbor-Net network, and neighbor-joining tree analyses. The haplotype analysis of mtDNA showed the initial dispersion of goats to Cameroon and central Africa from north-east Africa following the Nile Delta. Whereas, the approximate Bayesian computation indicated Cameroon goats were separated from Moroccan goats after 506 generations in later times (~1518 YA), as supported by the phylogenetic net-work and admixture outputs. Overall, indigenous goats in Cameroon show weak phylogenetic structure, suggesting either extensive intermixing.

The Origin of a Coastal Indigenous Horse Breed in China Revealed by Genome-Wide SNP Data

The Jinjiang horse is a unique Chinese indigenous horse breed distributed in the southern coastal areas, but the ancestry of Jinjiang horses is not well understood. Here, we used Equine SNP70 Bead Array technology to genotype 301 horses representing 10 Chinese indigenous horse breeds, and we integrated the published genotyped data of 352 individuals from 14 foreign horse breeds to study the relationships between Jinjiang horses and horse breeds from around the world. Principal component analysis (PCA), linkage disequilibrium (LD), runs of homozygosity (ROH) analysis, and ancestry estimating methods were conducted to study the population relationships and the ancestral sources and genetic structure of Jinjiang horses. The results showed that there is no close relationship between foreign horse breeds and Jinjiang horses, and Jinjiang horses shared a similar genetic background with Baise horses. TreeMix analysis revealed that there was gene flow from Chakouyi horses to Jinjiang horses. The ancestry analysis showed that Baise horses and Chakouyi horses are the most closely related ancestors of Jinjiang horses. In conclusion, our results showed that Jinjiang horses have a native origin and that Baise horses and Chakouyi horses were key ancestral sources of Jinjiang horses. The study also suggested that ancient trade activities and the migration of human beings had important effects on indigenous horse breeds in China.

Genetic Basis of Phenotypic Differences Between Chinese Yunling Black Goats and Nubian Goats Revealed by Allele-Specific Expression in Their F1 Hybrids

Chinese Yunling black goats and African Nubian goats are divergent breeds showing significant differences in body size, milk production, and environmental adaptation. However, the genetic mechanisms underlying these phenotypic differences remain to be elucidated. In this report, we provide a detailed portrait of allele-specific expression (ASE) from 54 RNA-Seq analyses across six tissues from nine F1 hybrid offspring generated by crossing the two breeds combined with 13 genomes of the two breeds. We identified a total of 524 genes with ASE, which are involved in bone development, muscle cell differentiation, and the regulation of lipid metabolic processes. We further found that 38 genes with ASE were also under directional selection by comparing 13 genomes of the two breeds; these 38 genes play important roles in metabolism, immune responses, and the adaptation to hot and humid environments. In conclusion, our study shows that the exploration of genes with ASE in F1 hybrids provides an efficient way to understand the genetic basis underlying the phenotypic differences of two diverse goat breeds.

Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds

Both natural and artificial selection are among the main driving forces shaping genetic variation across the genome of livestock species. Selection typically leaves signatures in the genome, which are often characterized by high genetic differentiation across breeds and/or a strong reduction in genetic diversity in regions associated with traits under intense selection pressure. In this study, we evaluated selection signatures and genomic inbreeding coefficients, FROH, based on runs of homozygosity (ROH), in six Ugandan goat breeds: Boer (n = 13), and the indigenous breeds Karamojong (n = 15), Kigezi (n = 29), Mubende (n = 29), Small East African (n = 29), and Sebei (n = 29). After genotyping quality control, 45,294 autosomal single nucleotide polymorphisms (SNPs) remained for further analyses. A total of 394 and 6 breed-specific putative selection signatures were identified across all breeds, based on marker-specific fixation index (FST-values) and haplotype differentiation (hapFLK), respectively. These regions were enriched with genes involved in signaling pathways associated directly or indirectly with environmental adaptation, such as immune response (e.g., IL10RB and IL23A), growth and fatty acid composition (e.g., FGF9 and IGF1), and thermo-tolerance (e.g., MTOR and MAPK3). The study revealed little overlap between breeds in genomic regions under selection and generally did not display the typical classic selection signatures as expected due to the complex nature of the traits. In the Boer breed, candidate genes associated with production traits, such as body size and growth (e.g., GJB2 and GJA3) were also identified. Furthermore, analysis of ROH in indigenous goat breeds showed very low levels of genomic inbreeding (with the mean FROH per breed ranging from 0.8% to 2.4%), as compared to higher inbreeding in Boer (mean FROH = 13.8%). Short ROH were more frequent than long ROH, except in Karamojong, providing insight in the developmental history of these goat breeds. This study provides insights into the effects of long-term selection in Boer and indigenous Ugandan goat breeds, which are relevant for implementation of breeding programs and conservation of genetic resources, as well as their sustainable use and management.

Genome-wide association study of body morphological traits in Sudanese goats

Long-term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome-wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185-scaffold1365-620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482-scaffold89-467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1-Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes.