The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

Anthropogenic events and responses to environmental stress are shaping the genomes of Ethiopian indigenous goats

Anthropological and biophysical processes have shaped livestock genomes over Millenia and can explain their current geographic distribution and genetic divergence. We analyzed 57 Ethiopian indigenous domestic goat genomes alongside 67 equivalents of east, west, and north-west African, European, South Asian, Middle East, and wild Bezoar goats. Cluster, ADMIXTURE (K = 4) and phylogenetic analysis revealed four genetic groups comprising African, European, South Asian, and wild Bezoar goats. The Middle Eastern goats had an admixed genome of these four genetic groups. At K = 5, the West African Dwarf and Moroccan goats were separated from East African goats demonstrating a likely historical legacy of goat arrival and dispersal into Africa via the coastal Mediterranean Sea and the Horn of Africa. FST, XP-EHH, and Hp analysis revealed signatures of selection in Ethiopian goats overlaying genes for thermo-sensitivity, oxidative stress response, high-altitude hypoxic adaptation, reproductive fitness, pathogen defence, immunity, pigmentation, DNA repair, modulation of renal function and integrated fluid and electrolyte homeostasis. Notable examples include TRPV1 (a nociception gene); PTPMT1 (a critical hypoxia survival gene); RETREG (a regulator of reticulophagy during starvation), and WNK4 (a molecular switch for osmoregulation). These results suggest that human-mediated translocations and adaptation to contrasting environments are shaping indigenous African goat genomes.

A Review on Indigenous Goats of East Africa: A Case for Conservation and Management

Indigenous goats are important in the livelihoods of rural households in East African countries. This is due to their ability to produce and reproduce in different environments and climatic conditions. Even though these indigenous goats are important, there is little available information on the genetic characterization of these breeds in Africa and at the international level. This paper reviews the status of indigenous goats, highlighting their production systems, phenotypic and genetic characteristics, and genetic diversity, and proposes potential ways for sustainable improvement and conservation in East African countries. Most households use traditional production systems with various uncharacterized goat breeds and ecotypes, which are hence named after the tribe or locality in which they are found. Most of these goats are classified as small East African breeds, with significant variability in morphological features. Some of the challenges to goat production in this region are indiscriminate crossbreeding, lack of pedigree records, parasites and disease incidences, low-quality pastures, and low levels of management. There is a need for a collaborative approach amongst the actors in goat breeding value chains as well as integrating modern genomic tools into breeding programs to enhance selection. This will ensure the resilience and sustainability of these unique indigenous goat populations in East Africa.

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.

Genome-Wide Association Study of Body Conformation Traits in Tashi Goats (Capra hircus)

Identifying genetic markers of economically valuable traits has practical benefits for the meat goat industry. To better understand the genomic variations influencing body conformation traits, a genome-wide association study was performed on Tashi goats, an indigenous Chinese goat breed. A total of 155 Tashi goats were phenotyped for eight body conformation traits: body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Then, 100 Tashi goats were randomly selected for whole-genome sequencing and genotyped. We obtained 1676.4 Gb of raw data with an average sequencing depth of 6.2X. Clean reads were aligned to the ARS1.2 reference genome, and 11,257,923 single nucleotide polymorphisms (SNPs) were identified. The structure analysis showed that these Tashi goats were almost not genetically related. The 109, 20, 52, 14, 62, 51, 70, and 7 SNPs were significantly associated with body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Within the ±500 kb region of significant SNPs, 183 genes were annotated. The most significantly enriched KEGG pathway was "olfactory transduction", and the most significantly enriched gene ontology (GO) terms were "cellular process", "cellular anatomical entity", and "molecular transducer activity". Interestingly, we found several SNPs on chromosomes 10 and 11 that have been identified multiple times for all eight body conformation traits located in two fragments (114 kb and 1.03 Mb). In chr.10:25988403-26102739, the six SNPs were tightly linked, the TACTAG genotype was the highest at 91.8%, and the FNTB (Farnesyltransferase, CAAX Box Beta) and CHURC1 (Churchill Domain Containing 1) genes were located. In chr.11:88216493-89250659, ten SNPs were identified with several dependent linkage disequilibrium (LD) blocks, and seven related genes were annotated, but no significant SNP was located in them. Our results provide valuable biological information for improving growth performance with practical applications for genomic selection in goats.

Coping with extremes: the rumen transcriptome and microbiome co-regulate plateau adaptability of Xizang goat

The interactions between the rumen microbiota and the host are crucial for the digestive and absorptive processes of ruminants, and they are heavily influenced by the climatic conditions of their habitat. Owing to the harsh conditions of the high-altitude habitat, little is known about how ruminants regulate the host transcriptome and the composition of their rumen microbiota. Using the model species of goats, we examined the variations in the rumen microbiota, transcriptome regulation, and climate of the environment between high altitude (Lhasa, Xizang; 3650 m) and low altitude (Chengdu, Sichuan, China; 500 m) goats. The results of 16 S rRNA sequencing revealed variations in the abundance, diversity, and composition of rumen microbiota. Papillibacter, Quinella, and Saccharofermentans were chosen as potential microbes for the adaptation of Xizang goats to the harsh climate of the plateau by the Spearman correlation study of climate and microbiota. Based on rumen transcriptome sequencing analysis, 244 genes were found to be differentially expressed between Xizang goats and low-altitude goats, with 127 genes showing up-regulation and 117 genes showing down-regulation. SLC26A9, GPX3, ARRDC4, and COX1 were identified as potential candidates for plateau adaptation in Xizang goats. Moreover, the metabolism of fatty acids, arachidonic acids, pathway involving cytokines and their receptors could be essential for adaptation to plateau hypoxia and cold endurance. The expression of GPX3, a gene linked to plateau acclimatization in Xizang goats, was linked to the abundance of Anaerovibrio, and the expression of SLC26A9 was linked to the quantity of Selenomonas, according to ruminal microbiota and host Spearman correlation analysis. Our findings imply that in order to adapt harsh plateau conditions, Xizang goats have evolved to maximize digestion and absorption as well as to have a rumen microbiota suitable for the composition of their diet.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Cranial form differences in goats by breed and domestic status

Domestic goats (Capra hircus) are globally represented by over 300 breeds, making them a useful model for investigating patterns of morphological change related to domestication. However, they have been little studied, likely due to their poor representation in museum collections and the difficulty in obtaining truly wild goat (Capra aegagrus, the bezoar) samples. Similar studies on other species reveal that domestication correlates with craniofacial alterations in domestics, which are non-uniform and often species-specific. Here, we use three-dimensional geometric morphometric methods (3DGMM) to describe and quantify cranial shape variation in wild (n = 21) versus domestic (n = 54) goats. We find that mean cranial shapes differ significantly between wild and domestic goats as well as between certain breeds. The detected differences are lower in magnitude than those reported for other domestic groups, possibly explained by the fewer directions of artificial selection in goat breeding, and their low global genetic diversity compared to other livestock. We also find tooth-row length reduction in the domestics, suggestive of rostral shortening-a prediction of the "domestication syndrome" (DS). The goat model thus expands the array-and combinations of-morphological changes observed under domestication, notably detecting alterations to the calvarium form which could be related to the ~ 15% brain size reduction previously reported for domestic compared to wild goats. The global success of domestic goats is due more to their ability to survive in a variety of harsh environments than to systematized human management. Nonetheless, their domestication has resulted in a clear disruption from the wild cranial form, suggesting that even low-intensity selection can lead to significant morphological changes under domestication.

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

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

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

Morphophysiological Responses of the Goat Mammary Gland to Water Scarcity in Arid and Semi-Arid Environments: Are They Enough to Generate Adaptation to New Climatic Challenges?

Due to climate change, diverse territories of the planet will suffer from water restrictions. Goats are perceived as the most resilient ruminants in this scenario. So, various studies have focused on describing how a lower water intake influences milk production, especially in breeds adapted to desert environments. In water-stress situations, goats lose up to 32% of their body weight (BW), the rate of passage is reduced, and the digestibility of the feed increases. When goats consume water again, the rumen prevents hemolysis and osmotic shock from occurring. Regarding milk production, the response varies depending on the breed and the level of water restriction, maintaining the milk volume or reducing it by up to 41%. Systemically, it decreases the urinary volume and glomerular filtration rate, increasing blood osmolality and the vasopressin (ADH) concentration. Studies are scarce regarding changes in blood flow to the mammary gland, but there would be a reduction in blood flow velocity of up to 40% without changing blood pressure. New studies must be undertaken to determine which breeds or crosses are the best adapted to changing environmental conditions and to improve our understanding of the changes that occur at the morphophysiological level of the caprine mammary gland.

Comparison of genetic diversity and phylogenetic structure of BRCA1 gene of some domestic and wild sheep breeds in different countries

BRCA1 gene plays an important role in DNA damage repair, cell cycle, and transcription process regulation; hence it's called gate keeper. The current research aims to perform bioinformatics analyzes of the BRCA1 gene of different breeds of domestic and wild sheep from 49 breeds in 14 countries using the NCBI genome database. The desired sequences were aligned using MEGA11 software and a phylogenetic tree was drawn by Neighbor-Joining method. The number of mutations, nucleotide diversity, and haploid diversity were also analyzed using Dnaspv5 software. The analyses showed 296 polymorphisms, which led to the creation of 45 different haplotypes with a haplotype diversity of 0.035. Nucleotide diversity and average nucleotide differences among breeds were estimated as 0.259 and 0.052, respectively. The average genetic distance within the population of countries was calculated as 0.052. The amount of sequence conservation in this research was 0.313 on average, which indicates the high polymorphism of this gene and the emergence of new proteins. Tajima's D value in Tajima's neutrality test was -2.421, which was significant (p < 0.05). One of the reasons for the high genetic diversity in Iran's wild sheep population is the existence of forests and open environments, which prevent genetic drift and reduce inbreeding.

Domestication and breeding objective did not shape the interpretation of physical and social cues in goats (Capra hircus)

Artificial selection by humans, either through domestication or subsequent selection for specific breeding objectives, drives changes in animal cognition and behaviour. However, most previous cognitive research comparing domestic and wild animals has focused on companion animals such as canids, limiting any general claims about the effects of artificial selection by humans. Using a cognitive test battery, we investigated the ability of wild goats (non-domestic, seven subjects), dwarf goats (domestic, not selected for milk production, 15 subjects) and dairy goats (domestic, selected for high milk yield, 18 subjects) to utilise physical and social cues in an object choice task. To increase the heterogeneity of our test samples, data for domestic goats were collected by two experimenters at two research stations (Agroscope; Research Institute for Farm Animal Biology). We did not find performance differences between the three groups in the cognitive test battery for either physical or social cues. This indicates that for a domestic non-companion animal species, domestication and selection for certain breeding objectives did not measurably shape the physical and cognitive skills of goats.

Harvesting cereals at Tappeh Sang-e Chakhmaq and the introduction of farming in Northeastern Iran during the Neolithic

Tappeh Sang-e Chakhmaq is the only Neolithic site in Northeastern Iran, characterised by aceramic and ceramic levels corresponding to an occupation of 1500 years from the eighth to the end of the sixth millennium BCE. The Western and Eastern Mounds represent the oldest and longest occupation among the sites identified East of the Zagros, providing a unique context to explore the origin and spread of farming outside the core area of the Eastern Fertile Crescent. We present data about the first harvesting activities in the Northeastern Iranian Central Plateau by applying usewear and microtexture analysis through confocal microscopy on sickle gloss blades. Our results indicate a community of pioneer farmers who settled down in the area carrying with them both domestic cereals as well as advanced techniques of cereal cultivation. We demonstrate that most of the tools were used for harvesting cereals in a fully ripened state collected near the ground, indicating a well-established cereal cultivation strategy. The use of straight shafts with parallel inserts in Tappeh Sang-e Chakhmaq, as known in some sites in the Zagros, suggests the dispersal of farming practices and technologies from the Eastern Fertile Crescent north-eastward across Iran. We observe an evolution in the degree of ripeness of harvested cereals along the first four levels of occupation of the Western Mound, where semi-ripe harvesting is relatively important, suggesting that domestic cereals to be harvested before full maturity were introduced into the village. From the topmost of the Western Mound and along the occupation of the Eastern Mound, ripe harvesting is dominant, showing a well-established cultivation strategy of fully mature cereal. This shift could indicate an in-situ evolution towards a better-established agricultural technology, including harvesting riper crops, that would have resulted in higher yields, as cereals were collected when the grain was fully formed.

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

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

Whole-genome sequencing and genomic analysis of Norduz goat (Capra hircus)

Artificial and natural selective breeding of goats has resulted in many different goat breeds all around the world. Norduz goat is one of these breeds, and it is a local goat breed of Turkey. The goats are favorable due to pre-weaning viability and reproduction values compared to the regional breeds. Development in sequencing technologies has let to understand huge genomic structures and complex phenotypes. Until now, such a comprehensive study has not been carried out to understand the genomic structure of the Norduz goats, yet. In the study, the next-generation sequencing was carried out to understand the genomic structure of Norduz goat. Real-time PCR was used to evaluate prominent CNVs in the Norduz goat individuals. Whole genome of the goat was constructed with an average of 33.1X coverage level. In the stringent filtering condition, 9,757,980 SNPs, 1,536,715 InDels, and 290 CNVs were detected in the Norduz goat genome. Functional analysis of high-impact SNP variations showed that the classical complement activation biological process was affected significantly in the goat. CNVs in the goat genome were found in genes related to defense against viruses, immune response, and cell membrane transporters. It was shown that GBP2, GBP5, and mammalian ortholog GBP1, which are INF-stimulated GTPases, were found to be high copy numbers in the goats. To conclude, genetic variations mainly in immunological response processes suggest that Norduz goat is an immunologically improved goat breed and natural selection could take an important role in the genetical improvements of the goats.

Peeping into Mitochondrial Diversity of Andaman Goats: Unveils Possibility of Maritime Transport with Diversified Geographic Signaling

Andaman and Nicobar Islands, a part of South-East Asia, is enriched with the presence of native breeds of livestock (cattle, pig, goat) and poultry. There are two native goat breeds, viz., Andaman local goat and Teressa goat in Andaman and Nicobar Islands. However, to date, the origin and genetic makeup of these two breeds have not been detailed. Therefore, the present study describes the genetic makeup of Andaman goats through analysis of mitochondrial D-loop sequence for sequence polymorphism, phylogeographical signaling and population expansion events. The genetic diversity of the Teressa goat was less compared to the Andaman local goat due to its sole presence on Teressa Island. Out of 38 well-defined haplotypes of Andaman goats, the majority of haplotypes belonged to haplogroup A followed by haplogroup B and haplogroup D. The result of mismatch distribution and neutrality tests indicated no population expansion event of haplogroup A and B. Finally, based on poor geographical signaling, we hypothesize that Andaman goats have been imported to these Islands either through multidirectional diffusion or unidirectional diffusion. We justify our hypothesis of multidirectional diffusion on the basis of observation of the haplotype and nucleotide diversity of Andaman goats. Simultaneously, the probability of unidirectional diffusion of goats in these islands from the Indian subcontinent in different spells of domestication events through maritime routes cannot be ignored.

Isolation, morphological identification, and xylanase characteristics of anaerobic gut fungi Neocallimastix from Anatolian wild goat

This study shows the morphological identification of anaerobic fungal strains isolated from fecal samples of goats inhabiting Turkey and the effects of various metal ions and chemicals on extracellular xylanase production. Three different anaerobic gut fungi isolated from wild goats in Turkey were identified as Neocallimastix spp. xylanase, cellulase, and lichenase production were tested in culture supernatants, and the maximum-specific activities were found as 560.42 ± 9.39, 159.70 ± 3.88, and 157.36 ± 3.83 (μmol/min/mg protein), respectively. While the optimum temperature range of exo-xylanases was found as 40-50°C, their optimum pH range was determined as 6.0-6.5. Xylanase activity decreased in metal ions and other chemical reactants based on dose. The metal ion that significantly inhibited xylanase activity was Fe⁺³ . It was found that the ferric ions inhibited xylanase activity in all three anaerobic gut fungi by 30%-90% depending on molarity. On the contrary, the 1 mM concentrations of the Mn⁺² , Ba⁺² , Co⁺² , Cu⁺² , Sn⁺² , and Mg⁺² metal ions and the ethylenediaminetetraacetic acid and β-mercaptoethanol reagents had a positive effect at rates in the range of 3%-92%. In conclusion, these findings demonstrate that anaerobic gut fungus has very stable fibrolytic enzymes that need to be separated, as well and the existence of a unique resource for industrial applications.

Human and Animal Fascioliasis: Origins and Worldwide Evolving Scenario

Fascioliasis is a plant- and waterborne zoonotic parasitic disease caused by two trematode species: (i) Fasciola hepatica in Europe, Asia, Africa, the Americas, and Oceania and (ii) F. gigantica, which is restricted to Africa and Asia. Fasciolid liver flukes infect mainly herbivores as ruminants, equids, and camelids but also omnivore mammals as humans and swine and are transmitted by freshwater Lymnaeidae snail vectors. Two phases may be distinguished in fasciolid evolution. The long predomestication period includes the F. gigantica origin in east-southern Africa around the mid-Miocene, the F. hepatica origin in the Near-Middle East of Asia around the latest Miocene to Early Pliocene, and their subsequent local spread. The short postdomestication period includes the worldwide spread by human-guided movements of animals in the last 12,000 years and the more recent transoceanic anthropogenic introductions of F. hepatica into the Americas and Oceania and of F. gigantica into several large islands of the Pacific with ships transporting livestock in the last 500 years. The routes and chronology of the spreading waves followed by both fasciolids into the five continents are redefined on the basis of recently generated knowledge of human-guided movements of domesticated hosts. No local, zonal, or regional situation showing disagreement with historical records was found, although in a few world zones the available knowledge is still insufficient. The anthropogenically accelerated evolution of fasciolids allows us to call them "peridomestic endoparasites." The multidisciplinary implications for crucial aspects of the disease should therefore lead the present baseline update to be taken into account in future research studies.

The mammalian brain under domestication: Discovering patterns after a century of old and new analyses

Comparisons of wild and domestic populations have established brain reduction as one of the most consistent patterns correlated with domestication. Over a century of scholarly work has been devoted to this subject, and yet, new data continue to foster its debate. Current arguments, both for and against the validity of brain reduction occurring in domestic taxa, have repeatedly cited a small set of reviews on this subject. The original works, their sampling, methodological details, and nuances of results that would be key to establishing validity, particularly in light of new data, have not been investigated. To facilitate and encourage a more informed discussion, we present a comprehensive review of original brain reduction literature for four mammalian clades: Artiodactyla, Perissodactyla, Carnivora, and Glires. Among these are studies that generated the most cited brain reduction values in modern domestication literature. In doing so, we provide a fairer stage for the critique of traits associated with domestication. We conclude that while brain reduction magnitudes may contain error, empirical data collectively support the reduction in brain size and cranial capacity for domestic forms.

Analysis of ancient and modern horse genomes reveals the critical impact of lncRNA-mediated epigenetic regulation on horse domestication

Background: The domestication of horses has played critical roles in human civilizations. The excavation of ancient horse DNA provides crucial data for studying horse domestication. Studies of horse domestication can shed light on the general mechanisms of animal domestication.

Objective: We wish to explore the gene transcription regulation by long noncoding RNAs (lncRNAs) that influence horse domestication.

Methods: First, we assembled the ancient DNA sequences of multiple horses at different times and the genomes of horses, donkeys, and Przewalski horses. Second, we extracted sequences of lncRNA genes shared in ancient horses and sequences of lncRNA genes and the promoter regions of domestication-critical genes shared in modern horses, modern donkeys, and Przewalski horses to form two sample groups. Third, we used the LongTarget program to predict potential regulatory interactions between these lncRNAs and these domestication-critical genes and analyzed the differences between the regulation in ancient/modern horses and between horses/donkeys/Przewalski horses. Fourth, we performed functional enrichment analyses of genes that exhibit differences in epigenetic regulation.

Results: First, genes associated with neural crest development and domestication syndrome are important targets of lncRNAs. Second, compared with undomesticated Przewalski horses, more lncRNAs participate in the epigenetic regulation in modern horses and donkeys, suggesting that domestication is linked to more epigenetic regulatory changes. Third, lncRNAs’ potential target genes in modern horses are mainly involved in two functional areas: 1) the nervous system, behavior, and cognition, and 2) muscle, body size, cardiac function, and metabolism.

Conclusion: Domestication is linked to substantial epigenetic regulatory changes. Genes associated with neural crest development and domestication syndrome underwent noticeable lncRNA-mediated epigenetic regulation changes during horse domestication.

A novel lineage of the Capra genus discovered in the Taurus Mountains of Turkey using ancient genomics

Direkli Cave, located in the Taurus Mountains of southern Turkey, was occupied by Late Epipaleolithic hunters-gatherers for the seasonal hunting and processing of game including large numbers of wild goats. We report genomic data from new and published Capra specimens from Direkli Cave and, supplemented with historic genomes from multiple Capra species, find a novel lineage best represented by a ~14,000 year old 2.59 X genome sequenced from specimen Direkli4. This newly discovered Capra lineage is a sister clade to the Caucasian tur species (Capra cylindricornis and Capra caucasica), both now limited to the Caucasus region. We identify genomic regions introgressed in domestic goats with high affinity to Direkli4, and find that West Eurasian domestic goats in the past, but not those today, appear enriched for Direkli4-specific alleles at a genome-wide level. This forgotten 'Taurasian tur' likely survived Late Pleistocene climatic change in a Taurus Mountain refuge and its genomic fate is unknown.

The genetic diversity and structure in the European polecat were not affected by the introduction of the American mink in Poland

The introduction and expansion of an invasive non-native species could have important consequences for the genetic patterns and processes of native species, moreover if the new arrival competes strongly for resources and space. This may result in the demographic decline of the native species. Knowing the effects on the levels of genetic diversity and structure in native species is key in terms of their conservation. We analysed temporal (over 50 years) genetic variation of the population of the European polecat (Mustela putorius), a species under threat in several European countries, in the Białowieża Primeval Forest (BPF), Poland, before and after the invasion of the American mink (Neovison vison). Using 11 microsatellite loci and a fragment of the mitochondrial control region we show that levels of diversity changed in the polecat population over 53 generations (over the period 1959–2012) and after the invasion of mink. When compared with other threatened European polecat populations, high levels of diversity are observed in the population in BPF in both periods, as well as in other areas in Poland. Our data shows that genetic structure was not present either before or after the mink invasion in BPF. This would suggest that the polecat population in Poland was not affected by invasive species and other negative factors and would be a potential good source of individuals for captive breeding or genetic rescue conservation management actions in areas where such actions are needed, for example the UK.

Paleobiogeographical origins of Fasciola hepatica and F. gigantica in light of new DNA sequence characteristics of F. nyanzae from hippopotamus

Fascioliasis is a highly pathogenic disease affecting humans and livestock worldwide. It is caused by the liver flukes Fasciola hepatica transmitted by Galba/Fossaria lymnaeid snails in Europe, Asia, Africa, the Americas and Oceania, and F. gigantica transmitted by Radix lymnaeids in Africa and Asia. An evident founder effect appears in genetic studies as the consequence of their spread by human-guided movements of domestic ruminants, equines and Old World camelids in the post-domestication period from the beginning of the Neolithic. Establishing the geographical origins of fasciolid expansion is multidisciplinary crucial for disease assessment. Sequencing of selected nuclear ribosomal and mitochondrial DNA markers of F. nyanzae infecting hippopotamuses (Hippopotamus amphibius) in South Africa and their comparative analyses with F. hepatica and F. gigantica, and the two Fascioloides species, Fs. jacksoni from Asian elephants and Fs. magna from Holarctic cervids, allow to draw a tuned-up evolutionary scenario during the pre-domestication period. Close sequence similarities indicate a direct derivation of F. hepatica and F. gigantica from F. nyanzae by speciation after host capture phenomena. Phylogenetic reconstruction, genetic distances and divergence estimates fully fit fossil knowledge, past interconnecting bridges between continents, present fasciolid infection in the wild fauna, and lymnaeid distribution. The paleobiogeographical analyses suggest an origin for F. gigantica by transfer from primitive hippopotamuses to grazing bovid ancestors of Reduncinae, Bovinae and Alcelaphinae, by keeping the same vector Radix natalensis in warm lowlands of southeastern Africa in the mid-Miocene, around 13.5 mya. The origin of F. hepatica should have occurred after capture from primitive, less amphibious Hexaprotodon hippopotamuses to mid-sized ovicaprines as the wild bezoar Capra aegagrus and the wild mouflon Ovis gmelini, and from R. natalensis to Galba truncatula in cooler areas and mountainous foothills of Asian Near East in the latest Miocene to Early Pliocene, around 6.0 to 4.0 mya and perhaps shortly afterwards.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

Analysis of World-Scale Mitochondrial DNA Reveals the Origin and Migration Route of East Asia Goats

Despite much attention on the history of goat evolution, information on origin, demographic history, and expansion route remains controversial. To address these questions, we collected 4,189 published goat DNA sequences including 1,228 sequences from 57 breeds in China and 2,961 sequences including 193 goat breeds from 71 other countries and carried out an integrated analysis. We found goat breeds from South China had the highest genetic diversity of lineage B, and subclades B2 only were found in Southwest China, suggesting that lineage B (particularly, subclade B2) probably originated from Southwest China and its surrounding areas. In addition, in this study, we found that lineage A from South China also presented higher genetic diversity and earlier expansion time (10, 606 years ago), even earlier than breeds from the Middle East. Hence, we speculated that South China and surrounding areas were the origin of lineage B and also the transportation hub for lineage A spreading to North China and Southwest Asia. Furthermore, according to the analysis of correlation between genetic differentiation value λ1 and λ2 and geographical distance, we further confirmed two phases of migration in goat breeds of North China. These results will contribute to a better understanding of the origin and migration history of domestic goat.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Dieudinné K, Vumilia Kizungu R, Ochieng JW, Pelle R. Haplotype analysis of the mitochondrial DNA d-loop region reveals the maternal origin and historical dynamics among the indigenous goat populations in east and west of the Democratic Republic of Congo

This study aimed at assessing haplotype diversity and population dynamics of three Congolese indigenous goat populations that included Kasai goat (KG), small goat (SG), and dwarf goat (DG) of the Democratic Republic of Congo (DRC). The 1169 bp d-loop region of mitochondrial DNA (mtDNA) was sequenced for 339 Congolese indigenous goats. The total length of sequences was used to generate the haplotypes and evaluate their diversities, whereas the hypervariable region (HVI, 453 bp) was analyzed to define the maternal variation and the demographic dynamic. A total of 568 segregating sites that generated 192 haplotypes were observed from the entire d-loop region (1169 bp d-loop). Phylogenetic analyses using reference haplotypes from the six globally defined goat mtDNA haplogroups showed that all the three Congolese indigenous goat populations studied clustered into the dominant haplogroup A, as revealed by the neighbor-joining (NJ) tree and median-joining (MJ) network. Nine haplotypes were shared between the studied goats and goat populations from Pakistan (1 haplotype), Kenya, Ethiopia and Algeria (1 haplotype), Zimbabwe (1 haplotype), Cameroon (3 haplotypes), and Mozambique (3 haplotypes). The population pairwise analysis (FST ) indicated a weak differentiation between the Congolese indigenous goat populations. Negative and significant (p-value <.05) values for Fu's Fs (-20.418) and Tajima's (-2.189) tests showed the expansion in the history of the three Congolese indigenous goat populations. These results suggest a weak differentiation and a single maternal origin for the studied goats. This information will contribute to the improvement of the management strategies and long-term conservation of indigenous goats in DRC.

Genome-Wide Association Study of Body Conformation Traits by Whole Genome Sequencing in Dazu Black Goats

Simple Summary

Body conformation traits are economically important in the goat meat industry. Good growth performance in goats, including an accelerated growth rate, can improve carcass weight and meat yield. The identification of genetic variants associated with these traits provides a basis for the genetic improvement of growth performance. In this study, we measured six body conformation traits, including body height, body length, cannon circumference, chest depth, chest width, and heart girth. By a genome-wide association study of a Chinese meat goat breed, 53 significant single nucleotide polymorphisms and 42 candidate genes associated with these traits were detected. These findings improve our understanding of the genetic basis of body conformation traits in goats.

Abstract

Identifying associations between genetic markers and economic traits has practical benefits for the meat goat industry. To better understand the genomic regions and biological pathways contributing to body conformation traits of meat goats, a genome-wide association study was performed using Dazu black goats (DBGs), a Chinese indigenous goat breed. In particular, 150 DBGs were genotyped by whole-genome sequencing, and six body conformation traits, including body height (BH), body length (BL), cannon circumference (CC), chest depth (CD), chest width (CW), and heart girth (HG), were examined. In total, 53 potential SNPs were associated with these body conformation traits. A bioinformatics analysis was performed to evaluate the genes located close to the significant SNPs. Finally, 42 candidate genes (e.g., PSTPIP2, C7orf57, CCL19, FGF9, SGCG, FIGN, and SIPA1L) were identified as components of the genetic architecture underlying body conformation traits. Our results provide useful biological information for the improvement of growth performance and have practical applications for genomic selection in goats.

Mitochondrial DNA D-loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro-ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D-loop. High haplotype (H d = 0.9619-0.9945) and nucleotide (π = 0.0120-0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro-ecological regions of the country.

Genome-wide scans for signatures of selection in Mangalarga Marchador horses using high-throughput SNP genotyping

BACKGROUND

The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system.

RESULTS

Three different methods were used to search for signals of selection: Tajima's D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array​ (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system.

CONCLUSIONS

Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima's D and iHS results point also to the presence of balancing selection in specific regions of the genome.

Detection of whole genome selection signatures of Pakistani Teddy goat

BACKGROUND

Natural and artificial selection tend to cause variability that contributes to shape the genome of livestock in a way that differentiates them among the animal kingdom. The particular aim here is to identify positive selection signatures with whole genome pooled-sequence data of Pakistani Teddy goat.

METHODS AND RESULTS

Paired-end alignment of 635,357,043 reads of Teddy goat with (ARS1) reference genome assembly was carried out. Pooled-Heterozygosity (Hp) and Tajima's D (TD) are applied for validation and getting better hits of selection signals, while pairwise F_ST statistics is conducted on Teddy vs. Bezoar (wild goat ancestor) for genomic differentiation, moreover annotation of regions under positive selection was also performed. Hp score with - ZHp > 5 detected six windows having highest hits on Chr. 29, 9, 25, 15 and 14 that harbor HRASLS5, LACE1 and AXIN1 genes which are candidate for embryonic development, lactation and body height. Secondly, - ZTD value of > 3.3 showed 4 windows with very strong hits on Chr.5 & 9 which harbor STIM1 and ADM genes related to body mass and weight. Lastly, - ZF_ST < - 5 generated four strong signals on Chr.5 & 12 harbor LOC102183233 gene. Other significant selection signatures encompass genes associated with wool production, prolificacy and coat colors traits in this breed.

CONCLUSIONS

In brief, this study identified the genes under selection in Pakistani Teddy goat that will be helpful to refining the marker-assisted breeding policies and converging required production traits within and across other goat breeds and to explore full genetic potential of this valued species of livestock.

Population validation of reproductive gene mutation loci and association with the litter size in Nubian goat

Litter size is an important component trait of doe reproduction. By improving it, production efficiency and economic benefits can be significantly provided. Genetic marker-assisted selection (MAS) based on proven molecular indicators could enhance the efficacy of goat selection, as well as litter size trait. Many molecular markers have been identified that they can be used to improve litter size in different goat breeds. However, the presence and value of these markers vary among goat breeds. In the present study, we used the reported loci on other breeds of goat as candidate loci to detect whether these loci appear in this Nubian goat population; then we proceed to genotype and detect surrounding loci (50 bp) by multiplex PCR and sequencing technology. As a result, 69 mutations (59 SNPs and 10 indels) were screened out from 23 candidate genes in Nubian goat population, 12 loci were significantly associated with the litter size of first-parity individuals; 5 loci were significantly associated with the litter size of second-parity individuals; 3 loci were significantly associated with the litter size of third-parity individuals. In addition, five loci were significantly associated with the average litter size. The additive effect value of KITLG: g.18047318 G>A in first parity, KITLG: g.18152042G>A in third parity, KISS-1: g.1341674 C>G in first parity, and GHR: g.32134187G>A in second parity exceed more than 0.40, and the preponderant alleles are G, C, A and G, respectively. Further, linkage disequilibrium analysis of 21 mutation loci shows that 3 haplotype blocks are formed, and the litter size of combination type AACC in KISS-1 gene and AAGG in KITLG gene are significantly lower than that of other combinations genotype in first parity (P<0.05). These findings can provide effective candidate DNA markers for selecting superior individuals in Nubian goat breeding.

Mitochondrial Sequence Variation, Haplotype Diversity, and Relationships Among Dromedary Camel-Types

Dromedary camels are outstanding livestock that developed efficient abilities to tolerate desert conditions. Many dromedary camel-types (i.e., named populations) exist but lack defined specific breed standards, registries, and breeders’ governing organizations. The breed status of dromedary camel-types can partly be assessed by exploring mitochondrial DNA (mtDNA) variation. Accordingly, this study aimed to examine the breed status and the inter-population relationships of dromedary camel-types by analyzing sequence variation in the mtDNA control region and in three coding genes [cytochrome b, threonine, and proline tRNA, and part of the displacement loop (D-loop)] (867 bp region). Tail hair samples (n = 119) that represent six camel-types from Kuwait were collected, extracted, sequenced, and compared to other publicly available sequences (n = 853). Within the sequenced mitochondrial region, 48 polymorphic sites were identified that contributed to 82 unique haplotypes across 37 camel-types. Haplotype names and identities were updated to avoid previous discrepancies. When all sequences were combined (n = 972), a nucleotide diversity of 0.0026 and a haplotype diversity of 0.725 was observed across the dromedary-types. Two major haplogroups (A and B) were identified and the B1 haplotype was predominant and found in almost all dromedary-types whereas the A haplotypes were more abundant in African regions. Non-metric multidimensional scaling revealed an increased similarity among Arabian Peninsula “Mezayen” camel-types, despite their defining coat colors. The relationships among dromedary camel-types can partly be explained by mtDNA. Future work aimed at a deeper understanding of camel-type breed status should focus on a high number of nuclear markers.

Ethiopian indigenous goats offer insights into past and recent demographic dynamics and local adaptation in sub-Saharan African goats

Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats.

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.

Distinguishing African bovids using Zooarchaeology by Mass Spectrometry (ZooMS): New peptide markers and insights into Iron Age economies in Zambia

Assessing past foodways, subsistence strategies, and environments depends on the accurate identification of animals in the archaeological record. The high rates of fragmentation and often poor preservation of animal bones at many archaeological sites across sub-Saharan Africa have rendered archaeofaunal specimens unidentifiable beyond broad categories, such as “large mammal” or “medium bovid”. Identification of archaeofaunal specimens through Zooarchaeology by Mass Spectrometry (ZooMS), or peptide mass fingerprinting of bone collagen, offers an avenue for identification of morphologically ambiguous or unidentifiable bone fragments from such assemblages. However, application of ZooMS analysis has been hindered by a lack of complete reference peptide markers for African taxa, particularly bovids. Here we present the complete set of confirmed ZooMS peptide markers for members of all African bovid tribes. We also identify two novel peptide markers that can be used to further distinguish between bovid groups. We demonstrate that nearly all African bovid subfamilies are distinguishable using ZooMS methods, and some differences exist between tribes or sub-tribes, as is the case for Bovina (cattle) vs. Bubalina (African buffalo) within the subfamily Bovinae. We use ZooMS analysis to identify specimens from extremely fragmented faunal assemblages from six Late Holocene archaeological sites in Zambia. ZooMS-based identifications reveal greater taxonomic richness than analyses based solely on morphology, and these new identifications illuminate Iron Age subsistence economies c. 2200–500 cal BP. While the Iron Age in Zambia is associated with the transition from hunting and foraging to the development of farming and herding, our results demonstrate the continued reliance on wild bovids among Iron Age communities in central and southwestern Zambia Iron Age and herding focused primarily on cattle. We also outline further potential applications of ZooMS in African archaeology.

Rediscovery and first genetic description of some poorly known tick species: Haemaphysalis kopetdaghica Kerbabaev, 1962 and Dermacentor raskemensis Pomerantzev, 1946

This study aimed to provide novel information for some poorly known/rare tick species collected from wild goats (Capra aegagrus) in the mountains of Eastern Anatolia, Turkey and to expand upon the available genetic data. The collected ticks were morphologically identified as Haemaphysalis kopetdaghica (all active stages, n = 140), Dermacentor raskemensis (adults, n = 7), Ixodes gibbosus (adults, n = 15), Rhipicephalus kohlsi (female, n = 1), and R. bursa (nymphs, n = 2). A total of 32 engorged ticks (6 larvae, 6 nymphs, and 20 females) collected were allowed to molt to the next stage or for egg laying and larval hatching, respectively. In addition, one R. kohlsi female (previously confirmed by SEM microscopy) collected from a wild goat in the neighboring province of Erzurum was included in this study for further genetic comparison. The partial mitochondrial 16S rDNA and cytochrome c oxidase subunit 1 (barcoding regions) genes of each tick species were sequenced. All DNA samples obtained from the ticks were checked by PCR for the presence of Anaplasma spp., Babesia spp., Borrelia burgdorferi sensu lato, spotted fever group rickettsiae, and Theileria spp., but were found to be negative. Phylogenetic analyses of the 16S rDNA and cox1 genes were performed using the ML method to determine their genetic relationship with related ticks. As a result, this study has: i) rediscovered and provided two new tick records (H. kopetdaghica and D. raskemensis) for Turkey, ii) provided the first genetic data for H. kopetdaghica and D. raskemensis and revealed their phylogenetic relationships, iii) characterized the cox1 region of I. gibbosus for the first time, and iv) revealed significant genetic diversity between R. kohlsi from Anatolia and R. kohlsi from Oman, suggesting that R. kohlsi could include a cryptic species.

Characterizing the Mitochondrial Diversity of Arbi Goats from Tunisia

Arbi is one of the main local goat breeds in Tunisia, representing an important economic resource in arid and hot areas where cattle and sheep cannot thrive successfully. In the current work, we have characterized the mitochondrial diversity of 26 Arbi goats by partially sequencing the mitochondrial D-loop region. These sequences plus 10 retrieved from GenBank were analyzed with the DnaSP v.5.10.1, evidencing the existence of 12 different haplotypes. Nucleotide and haplotype diversities were 0.02 and 0.96. Moreover, median-joining network analysis showed that all D-loop sequences from Arbi goats correspond to haplogroup A and that in general they do not cluster with sequences from other goat breeds. The high diversity that has been observed in North African goats is compatible with the maritime diffusion of the Neolithic package 10,000-7000 YBP. Moreover, there are evidences that local Tunisian breeds have been extensively crossed with highly productive transboundary breeds in order to improve meat and milk yields. These uncontrolled crossing practices may lead to the loss of alleles that play key roles in the adaptation of Tunisian local breeds to a harsh environment.

Keep Garfagnina alive. An integrated study on patterns of homozygosity, genomic inbreeding, admixture and breed traceability of the Italian Garfagnina goat breed

The objective of this study was to investigate the genetic diversity of the Garfagnina (GRF) goat, a breed that currently risks extinction. For this purpose, 48 goats were genotyped with the Illumina CaprineSNP50 BeadChip and analyzed together with 214 goats belonging to 9 other Italian breeds (~25 goats/breed), whose genotypes were available from the AdaptMap project [Argentata (ARG), Bionda dell'Adamello (BIO), Ciociara Grigia (CCG), Di Teramo (DIT), Garganica (GAR), Girgentana (GGT), Orobica (ORO), Valdostana (VAL) and Valpassiria (VSS)]. Comparative analyses were conducted on i) runs of homozygosity (ROH), ii) admixture ancestries and iii) the accuracy of breed traceability via discriminant analysis on principal components (DAPC) based on cross-validation. ROH analyses was used to assess the genetic diversity of GRF, while admixture and DAPC to evaluate its relationship to the other breeds. For GRF, common ROH (more than 45% in GRF samples) was detected on CHR 12 at, roughly 50.25-50.94Mbp (ARS1 assembly), which spans the CENPJ (centromere protein) and IL17D (interleukin 17D) genes. The same area of common ROH was also present in DIT, while a broader region (~49.25-51.94Mbp) was shared among the ARG, CCG, and GGT. Admixture analysis revealed a small region of common ancestry from GRF shared by BIO, VSS, ARG and CCG breeds. The DAPC model yielded 100% assignment success for GRF. Overall, our results support the identification of GRF as a distinct native Italian goat breed. This work can contribute to planning conservation programmes to save GRF from extinction and will improve the understanding of the socio-agro-economic factors related with the farming of GRF.

Evaluation of morphological traits and physiological variables of several Chinese goat breeds and their crosses

The current study was undertaken to evaluate some morphological traits of the goat breeds raised in Southwest China. The field experimentation and data collection were from 434 animals presenting seven breeds of the Dazu black goat (DBG; n = 203), Saanen milk goat (SMG; n = 50), Black Boer × Dazu black goat (BXC; n = 28), Hechuan white goat (HWG; n = 49), Inner Mongolia Cashmere goat (IMCG; n = 25), IMCG × DBG (F₁; n = 57) and F₁ × F₁ (F₂; n = 22). All studied animals were adult and selected to be at the same age (36.50 ± 0.75 months). After editing, more than 20 morphological and production traits like body condition score (BCS), testicle measurements, coat colour, fibre traits, skin colour, horn colour, horn shape, horn orientation, wattles, front hair, beard, ear shape, ear size, rump angle, hind leg angulation and physiological variables were analysed. BXC and DBG had dark coat colour, whilst SMG, HWG and IMCG had light colour, whilst F₁ and F₂ ranged from light to dark coat colour. Concerning BCS, the breeds BXC and DBG were characterized as fat goats, whilst SMG, HWG, F₁ and F₂ were average, whilst IMCG was thin. The maximum values for testis measurements were recorded for BXC, SMG and DBG. For fibre traits, IMCG, F₁ and F₂ were the most superior. BXC and DBG males have good fertility parameters. The highest values for rectal temperature, skin temperature and breath rate were recorded for SMG. These findings revealed the presence of a wide range of morphological differences among studied goat breeds. Such diversity in the performance of goat breeds raised in Southwest China is crucial to implement reliable selection strategies for breeding goats in this area.

Indonesian native goats (Capra hircus) reveal highest genetic frequency of mitochondrial DNA haplogroup B in the world

The objective of this study was to determine mtDNA sequences of the 481 bp HV1 region from two Indonesian native goat breeds, Kacang and Marica, to confirm the phylogeographic distribution of caprine haplogroup B in Southeast Asia. Based on these sequences, 12 haplotypes were observed and categorized into the predominant haplogroup B and minority haplogroup A, indicating that Indonesian native goats present the highest frequency (0.950) of the haplogroup B in the world. These results strongly emphasize previous observations of the haplogroup B frequencies tending to increase southeastward in Southeast Asia. Additionally, this suggests that goats primary bred into Southeast Asia might predominantly carry haplogroup B. This could essentially contribute to the understanding of the origin, propagation route and/or introgression history of Southeast Asian goats.

Assessing the Diversity and Population Substructure of Sarda Breed Bucks by Using Mtdna and Y-Chromosome Markers

A sample of 146 Sarda bucks from eight subregions of Sardinia, Italy (Nuorese, Barbagia, Baronia, Ogliastra, Sarrabus, Guspinese, Iglesiente, Sulcis) were characterized for Y-chromosome and mtDNA markers to assess the levels of population substructure. Five polymorphic loci (SRY, AMELY, ZFY, and DDX3Y) on the Y-chromosome were genotyped. The control region of mtDNA was sequenced as a source of complementary information. Analysis of Y-chromosome data revealed the segregation of 5 haplotypes: Y1A (66.43%), Y2 (28.57%), Y1C (3.57%), Y1B1 (0.71%), and Y1B2 (0.71%). High levels of Y-chromosome diversity were observed in populations from Southwest Sardinia. The FST values based on Y-chromosome and mtDNA data were low, although a paternal genetic differentiation was observed when comparing the Nuorese and Barbagia populations (Central Sardinia) with the Sulcis, Iglesiente, and Sarrabus populations (Southern Sardinia). AMOVA analysis supported the lack of population substructure. These results suggest the occurrence of a historical and extensive gene flow between Sarda goat populations from different locations of Sardinia, despite the fact that this island is covered by several large mountain ranges. Introgression with foreign caprine breeds in order to improve milk production might have also contributed to avoiding the genetic differentiation amongst Sarda populations.