Ancient goat genomes reveal mosaic domestication in the Fertile Crescent

Microbial dynamics and Pseudomonas natural product production in milk and dairy products

Covering: 2000 up to the first half of 2024Milk and its derived dairy products have long been integral to the human diet, with evidence of consumption dating back over 9000 years. Milk's high nutritional value renders dairy products an important element of human diet while also offering a fertile environment for microbial growth. Beneficial microorganisms in dairy products are often associated with biogenic and probiotic effects, whereas spoilage or pathogenic microorganisms can pose health risks. Fermentation is a key method to preserve milk. Whereas dairying practices in most parts of the world have been highly altered by industrialization over the past century, nomadic pastoralists in Mongolia notably retain a rich tradition of household-level dairy fermentation that has been practiced since 3000 BC. Milk-associated microorganisms produce a vast number of low molecular weight natural products that can mediate beneficial and detrimental interactions. Bacteria of the genus Pseudomonas are found in traditional Mongolian dairy products and are common contaminants in commercial dairy products, and they can strongly impact the quality and shelf-life of dairy products. These bacteria are well known for their ability to produce a variety of secondary metabolites, including nonribosomal (lipo)peptides, which are both structurally and functionally diverse. Lipopeptides can have antimicrobial properties, act as quorum sensing molecules, and contribute to biofilm formation due to their amphiphilic nature. Although often associated with spoilage, some of these natural products can also exhibit positive effects with potential beneficial applications in the dairy industry. This review aims to provide a comprehensive overview of the interplay between culinary fermentation and the production and activities of microbial-derived natural products.

Two species, two stories: divergent morphometric evolution of sheep and goats over the last 8000 years

Since their domestication about 10 500 years ago, sheep and goats have been central to human subsistence economies. Despite shared traits and frequent co-herding, these species are biologically distinct, with their uses reflecting diverse cultural practices. Advances in geometric morphometrics enable detailed analyses of their phenotypic evolution, shaped by natural and artificial selection, genetic drift and gene flow. These evolutionary mechanisms are often challenging to disentangle in archaeozoology, but patterns of morphometric variation can reveal evolutionary parallels and divergences. This study investigates the morphometric evolution of sheep and goat teeth in the northwestern Mediterranean basin over 8000 years. Dental remains, which preserve well in archaeological contexts, provide rich data for evolutionary studies, including species-level identifications. Using two-dimensional geometric morphometrics, third lower molars from 1602 sheep and 635 goats, spanning archaeological and modern specimens, were analysed. Molar size, shape and variance were compared across six chrono-cultural phases, tracing diachronic morphological changes. Distinct patterns emerge: sheep exhibit greater variability likely reflecting selective breeding for diverse purposes. Goats, in contrast, show greater uniformity. These findings underscore unique evolutionary trajectories for sheep and goats, offering new perspectives on their biocultural evolution within the dynamic environmental and anthropic contexts that shaped their current diversity.This article is part of the theme issue 'Unravelling domestication: multi-disciplinary perspectives on human and non-human relationships in the past, present and future'.

Genomic analyses reveal a lack of widespread strong selection in indigenous chickens

The study of domestication has been revolutionized with the advent of molecular genetics. Chickens, with their clear domestication history, emerge as an excellent model for study into the paths of evolution in domestication and improvement. Here we used genomic data from wild, indigenous, and commercial chickens to better understand how genetic drift and selection translate into their differentiations. Our investigation into the patterns of allelic change and divergence reveals a polygenic architecture governing genetic differentiation during domestication and improvement. We uncover distinctive population-specific differentiations in terms of genes and functions among wild, indigenous, and commercial chickens. Using Runs Of Homozygosity (ROH) based mixed model approach developed in this study, we identified only directional selection signatures occurring in wild and commercial chickens. Notably, our findings suggest that indigenous chickens serve as reservoirs of genetic diversity, necessary for rapid adaptation to new environments or subsequent modern breeding. This work provides unprecedented insights into the chicken domestication and improvement, and it illuminates our understanding of the domestication of other animal species.

Identification and expression profile analysis of circRNAs associated with goat uterus with different fecundity during estrous cycle

BACKGROUND

The Yunshang Black Goat, a distinguished meat goat breed native to China, is renowned for its superior reproductive capabilities. Despite this, there is considerable phenotypic variability within the breed. During the reproductive cycle, the uterus plays a pivotal role, with its functions evolving in line with the different stages of the cycle. This study focuses on the uterine tissues, including both the endometrium and myometrium, of Yunshang Black Goats with high fecundity (HF) and low fecundity (LF) during the proliferative (FP) and secretory (LP) phases of the estrous cycle. By examining these tissues, we aim to elucidate the underlying molecular and physiological mechanisms of the observed differences in reproductive success.

RESULTS

High-throughput sequencing was conducted, followed by bioinformatics analysis to identify the expression profiles of circRNAs. A total of 7,445 circRNAs were identified through the integration of findings from find_circ and CIRI2 software. Comparative analyses between the FPLF vs. FPHF and LPLF vs. LPHF revealed 149 differentially expressed (DE) circRNAs (94 up-regulated and 55 down-regulated) and 276 DE circRNAs (56 up-regulated and 220 down-regulated), respectively. The enrichment analysis indicated that the primary pathways involved were the Sphingolipid signaling pathway, MAPK signaling pathway, and GnRH signaling pathway, all of which are closely associated with cellular growth and development. Additionally, several key candidate genes were identified, such as FGF2 and MBTPS1. We also predicted a total of 281 miRNA-circRNA binding pairs, encompassing 263 circRNAs and 60 miRNAs, and simultaneously, 14 coding circRNAs were anticipated.

CONCLUSION

Based on the analysis, we have established the expression profiles of circRNAs during the follicular and luteal phases, respectively. Furthermore, using various analytical methods and data from high- and low-yield experimental control groups over different periods, we have identified multiple circRNAs that affect the high reproductive capacity of goats. Through enrichment analysis of the host genes of these circRNAs, we have discovered several key candidate genes. These findings provide fundamental data for the study of the molecular mechanisms underlying the fecundity of goats and pave the way for future genetic improvement strategies.

Genetic advancements and future directions in ruminant livestock breeding: from reference genomes to multiomics innovations

Ruminant livestock provide a rich source of products, such as meat, milk, and wool, and play a critical role in global food security and nutrition. Over the past few decades, genomic studies of ruminant livestock have provided valuable insights into their domestication and the genetic basis of economically important traits, facilitating the breeding of elite varieties. In this review, we summarize the main advancements for domestic ruminants in reference genome assemblies, population genomics, and the identification of functional genes or variants for phenotypic traits. These traits include meat and carcass quality, reproduction, milk production, feed efficiency, wool and cashmere yield, horn development, tail type, coat color, environmental adaptation, and disease resistance. Functional genomic research is entering a new era with the advancements of graphical pangenomics and telomere-to-telomere (T2T) gap-free genome assembly. These advancements promise to improve our understanding of domestication and the molecular mechanisms underlying economically important traits in ruminant livestock. Finally, we provide new perspectives and future directions for genomic research on ruminant genomes. We suggest how ever-increasing multiomics datasets will facilitate future studies and molecular breeding in livestock, including the potential to uncover novel genetic mechanisms underlying phenotypic traits, to enable more accurate genomic prediction models, and to accelerate genetic improvement programs.

Cost-effectiveness and other considerations for different research techniques applied in ancient DNA analysis

Ancient DNA (aDNA) analysis has developed rapidly since it first emerged in the 1980s, becoming an almost indispensable tool in anthropological and archaeological sciences. Earlier aDNA study was based on the polymerase chain reaction (PCR) technique, with which, unfortunately, modern DNA contamination and other authenticity issues were often incurred. These technical hurdles were soon overcome by application of advancements in the forms of the next generation sequencing (NGS) technique and others. However, since NGS requires money, time, and, in the case of large projects, manpower as well, genetic analysis of some ancient samples considered to be insignificant is commonly delayed or, in the worst cases, neglected entirely. We acknowledge that as a diagnostic tool in aDNA analysis, PCR is less accurate than NGS and more easily affected by modern DNA contamination; but it also has advantages, such as simplicity, time-saving, and greater ease of interpretation, among others. The role of PCR in aDNA analysis, then, should be reconsidered.

Whole-genome variants resource of 298 Saanen dairy goats

The Saanen breeds are often used as terminal sires for hybridization and play an important role in the global dairy food industry. However, there is still a lack of genomics information on the Saanen dairy goats. Whole-genome sequencing offers a promising approach to identify genetic markers associated with economic traits and discover new candidate genes. This can effectively utilize genetic resources to accelerate breeding processes and improve lactation performance in Saanen dairy goats. In this study, we present the genomes of 298 Saanen dairy goats. Through rigorous sequencing and quality control, we achieved an average sequencing depth of 14.6X, with 92.3% of high-quality (Q30 > 90%) data and an average mapping ratio of 99.9%, indicating reliable results. By comparing our data to a reference genome of Saanen dairy goats, we identified14.59 million single nucleotide polymorphisms (SNPs) and 1.34 million insertions-deletions (InDels). This dataset significantly contributes to enriching public databases in dairy goats and provides valuable resources for studying genetic diversity, improving breeds, and developing new varieties.

Archaeogenomic insights into commensalism and regional variation in pig management in Neolithic northwest Europe

The relationship between humans and pigs has changed dramatically since their domestication in southwest Asia and subsequent human-induced introduction into Europe. Introgression between incoming southwest Asian pigs and European boar resulted in the gradual replacement of southwest Asian ancestry in European pigs. However, we currently lack genomic data required to explore the regional trajectories, nature, and extent of contact between European boar and pigs that led to this turnover, and how this process was facilitated by human activity. We addressed this deficit by sequencing four Mesolithic boar and seven Neolithic pig samples from six archaeological sites in the Netherlands and Britain ranging from the Mesolithic (5500 BCE) to Neolithic (2500 BCE). Our data show that despite continuous gene flow with European boar, Neolithic European pigs show varying levels of southwest Asian ancestry. The low and varying southwest Asian ancestry in pigs from the Early Neolithic Dutch settlement Swifterbant indicates a high contribution of wild ancestry. The genetic profile, enriched δ15N values, on-site presence, and wide size distribution of Swifterbant Sus scrofa suggest a commensal relationship. Runs of homozygosity (ROH) imply that both closed-breeding and free-ranging management occurred in Neolithic communities, where the former showed an extreme burden of long segments of ROH. We further show selection signatures, associated with coat color and behavior, in Neolithic herds despite recurrent wild gene flow. Altogether, our results show distinct husbandry practices through space and time in Neolithic Europe, with heavy reliance on boar recruitment via the commensal pathway in northwest Europe.

The genomic legacy of aurochs hybridisation in ancient and modern Iberian cattle

Cattle (Bos taurus) play an important role in the life of humans in the Iberian Peninsula not just as a food source but also in cultural events. When domestic cattle were first introduced to Iberia, wild aurochs (Bos primigenius) were still present, leaving ample opportunity for mating (whether intended by farmers or not). Using a temporal bioarchaeological dataset covering eight millennia, we trace gene flow between the two groups. Our results show frequent hybridisation during the Neolithic and Chalcolithic, likely reflecting a mix of hunting and herding or relatively unmanaged herds, with mostly male aurochs and female domestic cattle involved. This is supported by isotopic evidence consistent with ecological niche sharing, with only a few domestic cattle possibly being managed. The proportion of aurochs ancestry in domestic cattle remains relatively constant from about 4000 years ago, probably due to herd management and selection against first generation hybrids, coinciding with other cultural transitions. The constant level of wild ancestry (~20%) continues into modern Western European breeds including Iberian cattle selected for aggressiveness and fighting ability. This study illuminates the genomic impact of human actions and wild introgression in the establishment of cattle as one of the most important domestic species today.

Maternal and paternal lineage analysis of Island Southeast Asian goats reveals continental propagation routes and introgression through the Indian ocean

The Philippines and Indonesia, located at the crossroads of historical human migrations in the Asia-Pacific region, lie in the southeasternmost part of Mainland Southeast Asia. Both countries are recognized as mega-diverse nations, hosting a wide variety of species. However, the history of goat propagation in Island Southeast Asia remains unclear due to limited molecular studies and the absence of archaeological evidence in this region. This study conducted phylogeographic analyses of mtDNA hypervariable region 1 and SRY 3'UTR sequences of goats from the Philippines and Indonesia, comparing them with Old World domestic goats. Philippine goat mtDNA haplotypes were classified into haplogroups B (152/206, 73.8%) and A (54/206, 26.2%), whereas all Indonesian goat mtDNA haplotypes belonged to haplogroup B (72/72, 100%). The distribution of mtDNA haplogroup B and regional differences suggest two propagation routes: one from China/Taiwan to the Philippines and another from Mainland Southeast Asia to Indonesia. The distribution of five SRY haplotypes demonstrated that haplotype Y2A, absent in Mainland Southeast Asia, was present in Island Southeast Asia, indicating gene introgression through the Indian Ocean from the African continent and/or southern Europe. This introgression may have been driven by human activities, such as European migrations between the 16th and 19th centuries.

Whole-Genome Selective Scans Detect Genes Associated with Cashmere Traits and Climatic Adaptation in Cashmere Goats (Capra hircus) in China

Background: Cashmere, valued for its exceptional softness and warmth, is a major focus in goat breeding due to its high economic importance. However, the molecular mechanisms underlying cashmere production remain largely unknown, hindering efforts to optimize yield and quality. Additionally, domestic goats exhibit remarkable adaptability to diverse climates, ranging from arid northern regions to humid southern areas, yet the genetic basis for these adaptations is poorly understood. Exploring the genetic factors driving cashmere production and climatic adaptation could provide crucial insights into trait evolution and support the development of breeding strategies for improved productivity and resilience. Methods: We utilized whole-genome resequencing data from 157 samples representing 14 goat populations to analyze the genetic diversity between cashmere and non-cashmere breeds. Additionally, we conducted the tests of selective sweeps (i.e., pairwise FST, θπ and XP-EHH) for cashmere traits and genome-environment association analysis (i.e., XtX statistic), respectively. Results: We identified strong selective signatures in previous reports (e.g., AKT3, FOXP1, FGF5, TGFBR3) and novel genes (e.g., ZEB1, ZNRF3, MAPK8IP3, MAPK8IP2, AXIN1) associated with cashmere traits. Further gene annotation and KEGG analyses showed that these genes were identified to be the most probable genes accounting for the cashmere traits. Also, we detected some genes such as PDGFRB, PRDM8, SLC26A2, SCAMP1, EPHX1, CDC25A, and POLK that played critical roles in the adaptation of goats to local climate variation. Conclusions: Collectively, our results provide novel insights into the genetic mechanisms underlying the cashmere traits and climatic adaptation, and also identified new genetic markers for genetic improvement in goats.

Ancient genomics and the origin, dispersal, and development of domestic sheep

The origins and prehistory of domestic sheep (Ovis aries) are incompletely understood; to address this, we generated data from 118 ancient genomes spanning 12,000 years sampled from across Eurasia. Genomes from Central Türkiye ~8000 BCE are genetically proximal to the domestic origins of sheep but do not fully explain the ancestry of later populations, suggesting a mosaic of wild ancestries. Genomic signatures indicate selection by ancient herders for pigmentation patterns, hornedness, and growth rate. Although the first European sheep flocks derive from Türkiye, in a notable parallel with ancient human genome discoveries, we detected a major influx of Western steppe-related ancestry in the Bronze Age.

Does strain-level persistence of lactobacilli in long-term back-slopped sourdoughs inform on domestication of food-fermenting lactic acid bacteria?

Sourdoughs are maintained by back-slopping over long time periods. To determine strain-level persistence of bacteria, we characterized four sourdoughs from bakeries over a period of 3.3, 11.0, 18.0, and 19.0 years. One sourdough included isolates of Levilactobacillus spp. and Fructilactobacillus spp. that differed by fewer than 10 single-nucleotide polymorphisms (SNPs) from the isolates obtained 3.3 years earlier and thus likely represent the same strain. Isolates of Levilactobacillus parabrevis differed by 200-300 SNPs; their genomes were under positive selection, indicating transmission from an external source. In two other sourdoughs, isolates of Fructilactobacillus sanfranciscensis that were obtained 11 and 18 years apart differed by 19 and 29 SNPs, respectively, again indicating repeated isolation of the same strain. The isolate of Fl. sanfranciscensis from the fourth sourdough differed by 45 SNPs from the isolate obtained 19 years previously. We thus identified strain-level persistence in three out of four long-term back-slopped sourdoughs, making it possible that strains persisted over periods that are long enough to allow bacterial speciation and domestication.IMPORTANCEThe assembly of microbial communities in sourdough is shaped by dispersal and selection. Speciation and domestication of fermentation microbes in back-slopped food fermentations have been documented for food-fermenting fungi including sourdough yeasts but not for bacteria, which evolve at a slower rate. Bacterial speciation in food fermentations requires strain-level persistence of fermentation microbes over hundreds or thousands of years. By documenting strain-level persistence in three out of four sourdoughs over a period of up to 18 years, we demonstrate that persistence over hundreds or thousands of years is possible, if not likely. We thus not only open a new perspective on fermentation control in bakeries but also support the possibility that all humans, despite their cultural diversity, share the same fermentation microbes.

Telomere-to-telomere genome assembly of a male goat reveals variants associated with cashmere traits

A complete goat (Capra hircus) reference genome enhances analyses of genetic variation, thus providing insights into domestication and selection in goats and related species. Here, we assemble a telomere-to-telomere (T2T) gap-free genome (2.86 Gb) from a cashmere goat (T2T-goat1.0), including a Y chromosome of 20.96 Mb. With a base accuracy of >99.999%, T2T-goat1.0 corrects numerous genome-wide structural and base errors in previous assemblies and adds 288.5 Mb of previously unresolved regions and 446 newly assembled genes to the reference genome. We sequence the genomes of five representative goat breeds for PacBio reads, and use T2T-goat1.0 as a reference to identify a total of 63,417 structural variations (SVs) with up to 4711 (7.42%) in the previously unresolved regions. T2T-goat1.0 was applied in population analyses of global wild and domestic goats, which revealed 32,419 SVs and 25,397,794 SNPs, including 870 SVs and 545,026 SNPs in the previously unresolved regions. Also, our analyses reveal a set of selective variants and genes associated with domestication (e.g., NKG2D and ABCC4) and cashmere traits (e.g., ABCC4 and ASIP).

Evolution and legacy of East Asian aurochs

Aurochs (Bos primigenius), once widely distributed in Afro-Eurasia, became extinct in the early 1600 s. However, their phylogeography and relative contributions to domestic cattle remain unknown. In this study, we analyzed 16 genomes of ancient aurochs and three mitogenomes of ancient bison (Bison priscus) excavated in East Asia, dating from 43,000 to 3,590 years ago. These newly generated data with previously published genomic information on aurochs as well as ancient/extant domestic cattle worldwide through genome analysis. Our findings revealed significant genetic divergence between East Asian aurochs and their European, Near Eastern, and African counterparts on the basis of both mitochondrial and nuclear genomic data. Furthermore, we identified evidence of gene flow from East Asian aurochs into ancient and present-day taurine cattle, suggesting their potential role in facilitating the environmental adaptation of domestic cattle.

Bronze Age cheese reveals human-Lactobacillus interactions over evolutionary history

Despite the long history of consumption of fermented dairy, little is known about how the fermented microbes were utilized and evolved over human history. Here, by retrieving ancient DNA of Bronze Age kefir cheese (∼3,500 years ago) from the Xiaohe cemetery, we explored past human-microbial interactions. Although it was previously suggested that kefir was spread from the Northern Caucasus to Europe and other regions, we found an additional spreading route of kefir from Xinjiang to inland East Asia. Over evolutionary history, the East Asian strains gained multiple gene clusters with defensive roles against environmental stressors, which can be a result of the adaptation of Lactobacillus strains to various environmental niches and human selection. Overall, our results highlight the role of past human activities in shaping the evolution of human-related microbes, and such insights can, in turn, provide a better understanding of past human behaviors.

Chromosome-level genome assembly of the cashmere goat

The goat, an early domesticated ruminant, is a reliable source of cashmere, meat and milk in global agricultural production. Despite this, the genome of cashmere-rich goats has yet to be characterized. Here, we assembled the nearly complete genome of a cashmere goat from a highly economically valuable Inner Mongolian Cashmere buck, utilizing a combination of PacBio HiFi, ONT ultra-long reads, and Hi-C technologies. The size of this genome is 2.76 Gb, with a contig N50 of 95.22 Mb. All assembled sequences were anchored onto 29 autosomes and both sex chromosomes, with only two gaps present on the X chromosome. We identified 1,333.29 Mb (48.26%) of repetitive sequences and predicted 22,480 protein-coding genes. Assembly quality assessment of the genome demonstrated that our assembled cashmere goat genome surpasses the continuity, completeness, and accuracy of other published goat genomes. Taken together, we provided the first cashmere goat assembly, bridging the gap in the genome of important economic breeds of domestic goats, and providing a valuable reference resource for goat genetics and genome research.

The goat pan-genome reveals patterns of gene loss during domestication

BACKGROUND

Unveiling genetic diversity features and understanding the genetic mechanisms of diverse goat phenotypes are pivotal in facilitating the preservation and utilization of these genetic resources. However, the total genetic diversity within a species can't be captured by the reference genome of a single individual. The pan-genome is a collection of all the DNA sequences that occur in a species, and it is expected to capture the total genomic diversity of the specific species.

RESULTS

We constructed a goat pan-genome using map-to-pan assemble based on 813 individuals, including 723 domestic goats and 90 samples from their wild relatives, which presented a broad regional and global representation. In total, 146 Mb sequences and 974 genes were identified as absent from the reference genome (ARS1.2; GCF_001704415.2). We identified 3,190 novel single nucleotide polymorphisms (SNPs) using the pan-genome analysis. These novel SNPs could properly reveal the population structure of domestic goats and their wild relatives. Presence/absence variation (PAV) analysis revealed gene loss and intense negative selection during domestication and improvement.

CONCLUSIONS

Our research highlights the importance of the goat pan-genome in capturing the missing genetic variations. It reveals the changes in genomic architecture during goat domestication and improvement, such as gene loss. This improves our understanding of the evolutionary and breeding history of goats.

The Population History of Domestic Sheep Revealed by Paleogenomes

Sheep was one of the first domesticated animals in Neolithic West Eurasia. The zooarchaeological record suggests that domestication first took place in Southwest Asia, although much remains unresolved about the precise location(s) and timing(s) of earliest domestication, or the post-domestication history of sheep. Here, we present 24 new partial sheep paleogenomes, including a 13,000-year-old Epipaleolithic Central Anatolian wild sheep, as well as 14 domestic sheep from Neolithic Anatolia, two from Neolithic Iran, two from Neolithic Iberia, three from Neolithic France, and one each from Late Neolithic/Bronze Age Baltic and South Russia, in addition to five present-day Central Anatolian Mouflons and two present-day Cyprian Mouflons. We find that Neolithic European, as well as domestic sheep breeds, are genetically closer to the Anatolian Epipaleolithic sheep and the present-day Anatolian and Cyprian Mouflon than to the Iranian Mouflon. This supports a Central Anatolian source for domestication, presenting strong evidence for a domestication event in SW Asia outside the Fertile Crescent, although we cannot rule out multiple domestication events also within the Neolithic Fertile Crescent. We further find evidence for multiple admixture and replacement events, including one that parallels the Pontic Steppe-related ancestry expansion in Europe, as well as a post-Bronze Age event that appears to have further spread Asia-related alleles across global sheep breeds. Our findings mark the dynamism of past domestic sheep populations in their potential for dispersal and admixture, sometimes being paralleled by their shepherds and in other cases not.

A label-free quantification method for assessing sex from modern and ancient bovine tooth enamel

Identification of the sex of modern, fossil and archaeological animal remains offers many insights into their demography, mortality profiles and domestication pathways. However, due to many-factors, sex determination of osteological remains is often problematic. To overcome this, we have developed an innovative protocol to determine an animal's sex from tooth enamel, by applying label-free quantification (LFQ) of two unique AmelY peptides 'LRYPYP' (AmelY;[M+2]2 + 404.7212 m/z) and 'LRYPYPSY' (AmelY;[M+2]2 + 529.7689 m/z) that are only present in the enamel of males. We applied this method to eight modern cattle (Bos taurus) of known sex, and correctly assigned them to sex. We then applied the same protocol to twelve archaeological Bos teeth from the Neolithic site of Beisamoun, Israel (8-th-7-th millennium BC) and determined the sex of the archaeological samples. Since teeth are usually better preserved than bones, this innovative protocol has potential to facilitate sex determination in ancient and modern bovine remains that currently cannot be sexed.

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.

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.

Facilitated introgression from domestic goat into Alpine ibex at immune loci

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1-2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

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.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

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.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats

BACKGROUND

Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate.

RESULTS

To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds.

CONCLUSIONS

The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

Selection signatures and landscape genomics analysis to reveal climate adaptation of goat breeds

Goats have achieved global prominence as essential livestock since their initial domestication, primarily owing to their remarkable adaptability to diverse environmental and production systems. Differential selection pressures influenced by climate have led to variations in their physical attributes, leaving genetic imprints within the genomes of goat breeds raised in diverse agroecological settings. In light of this, our study pursued a comprehensive analysis, merging environmental data with single nucleotide polymorphism (SNP) variations, to unearth indications of selection shaped by climate-mediated forces in goats. Through the examination of 43,300 SNPs from 51 indigenous goat breeds adapting to different climatic conditions using four analytical methods: latent factor mixed models (LFMM), F-statistics (Fst), Extended haplotype homozygosity across populations (XPEHH), and spatial analysis method (SAM), A total of 74 genes were revealed to display clear signs of selection, which are believed to be influenced by climatic conditions. Among these genes, 32 were consistently identified by at least two of the applied methods, and three genes (DENND1A, PLCB1, and ITPR2) were confirmed by all four approaches. Moreover, our investigation yielded 148 Gene Ontology (GO) terms based on these 74 genes, underlining pivotal biological pathways crucial for environmental adaptation. These pathways encompass functions like vascular smooth muscle contraction, cellular response to heat, GTPase regulator activity, rhythmic processes, and responses to temperature stimuli. Of significance, GO terms about endocrine regulation and energy metabolic responses, key for local adaptation were also uncovered, including biological processes, such as cell differentiation, regulation of peptide hormone secretion, and lipid metabolism. These findings contribute to our knowledge of the genetic structure of climate-triggered adaptation across the goat genome and have practical implications for marker-assisted breeding in goats.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries)

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Aşıklı Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Aşıklı Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Aşıklı Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Genomic and pangenomic analyses provide insights into the population history and genomic diversification of bottle gourd

Bottle gourd (Lagenaria siceraria (Mol.) Strandl.) is an economically important vegetable crop and one of the earliest domesticated crops. However, the population history and genomic diversification of bottle gourd have not been extensively studied. We generated a comprehensive bottle gourd genome variation map from genome sequences of 197 world-wide representative accessions, which enables a genome-wide association study for identifying genomic loci associated with resistance to zucchini yellow mosaic virus, and constructed a bottle gourd pangenome that harbors 1534 protein-coding genes absent in the reference genome. Demographic analyses uncover that domesticated bottle gourd originated in Southern Africa c. 12 000 yr ago, and subsequently radiated to the New World via the Atlantic drift and to Eurasia through the efforts of early farmers in the initial Holocene. The identified highly differentiated genomic regions among different bottle gourd populations harbor many genes contributing to their local adaptations such as those related to disease resistance and stress tolerance. Presence/absence variation analysis of genes in the pangenome reveals numerous genes including those involved in abiotic/biotic stress responses that have been under selection during the world-wide expansion of bottle gourds. The bottle gourd variation map and pangenome provide valuable resources for future functional studies and genomics-assisted breeding.

Runs of Homozygosity Detection and Selection Signature Analysis for Local Goat Breeds in Yunnan, China

Runs of Homozygosity (ROH) are continuous homozygous DNA segments in diploid genomes, which have been used to estimate the genetic diversity, inbreeding levels, and genes associated with specific traits in livestock. In this study, we analyzed the resequencing data from 10 local goat breeds in Yunnan province of China and five additional goat populations obtained from a public database. The ROH analysis revealed 21,029 ROH segments across the 15 populations, with an average length of 1.27 Mb, a pattern of ROH, and the assessment of the inbreeding coefficient indicating genetic diversity and varying levels of inbreeding. iHS (integrated haplotype score) was used to analyze high-frequency Single-Nucleotide Polymorphisms (SNPs) in ROH regions, specific genes related to economic traits such as coat color and weight variation. These candidate genes include OCA2 (OCA2 melanosomal transmembrane protein) and MLPH (melanophilin) associated with coat color, EPHA6 (EPH receptor A6) involved in litter size, CDKAL1 (CDK5 regulatory subunit associated protein 1 like 1) and POMC (proopiomelanocortin) linked to weight variation and some putative genes associated with high-altitude adaptability and immune. This study uncovers genetic diversity and inbreeding levels within local goat breeds in Yunnan province, China. The identification of specific genes associated with economic traits and adaptability provides actionable insights for utilization and conservation efforts.

Limited historical admixture between European wildcats and domestic cats

Domestic cats were derived from the Near Eastern wildcat (Felis lybica), after which they dispersed with people into Europe. As they did so, it is possible that they interbred with the indigenous population of European wildcats (Felis silvestris). Gene flow between incoming domestic animals and closely related indigenous wild species has been previously demonstrated in other taxa, including pigs, sheep, goats, bees, chickens, and cattle. In the case of cats, a lack of nuclear, genome-wide data, particularly from Near Eastern wildcats, has made it difficult to either detect or quantify this possibility. To address these issues, we generated 75 ancient mitochondrial genomes, 14 ancient nuclear genomes, and 31 modern nuclear genomes from European and Near Eastern wildcats. Our results demonstrate that despite cohabitating for at least 2,000 years on the European mainland and in Britain, most modern domestic cats possessed less than 10% of their ancestry from European wildcats, and ancient European wildcats possessed little to no ancestry from domestic cats. The antiquity and strength of this reproductive isolation between introduced domestic cats and local wildcats was likely the result of behavioral and ecological differences. Intriguingly, this long-lasting reproductive isolation is currently being eroded in parts of the species' distribution as a result of anthropogenic activities.

Demographic reconstruction of the Western sheep expansion from whole-genome sequences

As one of the earliest livestock, sheep (Ovis aries) were domesticated in the Fertile Crescent about 12,000-10,000 years ago and have a nearly worldwide distribution today. Most of our knowledge about the timing of their expansions stems from archaeological data but it is unclear how the genetic diversity of modern sheep fits with these dates. We used whole-genome sequencing data of 63 domestic breeds and their wild relatives, the Asiatic mouflon (O. gmelini, previously known as O. orientalis), to explore the demographic history of sheep. On the global scale, our analysis revealed geographic structuring among breeds with unidirectional recent gene flow from domestics into Asiatic mouflons. We then selected 4 representative breeds from Spain, Morocco, the United Kingdom, and Iran to build a comprehensive demographic model of the Western sheep expansion. We inferred a single domestication event around 11,000 years ago. The subsequent westward expansion is dated to approximately 7,000 years ago, later than the original Neolithic expansion of sheep and slightly predating the Secondary Product Revolution associated with wooly sheep. We see some signals of recent gene flow from an ancestral population into Southern European breeds which could reflect admixture with feral European mouflon. Furthermore, our results indicate that many breeds experienced a reduction of their effective population size during the last centuries, probably associated with modern breed development. Our study provides insights into the complex demographic history of Western Eurasian sheep, highlighting interactions between breeds and their wild counterparts.

Genome-Wide Association and Pathway Analysis of Carcass and Meat Quality Traits in Karachai Young Goats

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. The work aimed to search for candidate genes related to body measurements and body weight of Karachai goats and develop an experimental PCR-RV test system for genotyping significant SNPs. Comparison of GWAS results for ages 4 and 8 months revealed 58 common SNPs for significant genotypes. 11 common SNPs were identified for body weight, 4 SNPs-for group of traits withers height, rump height, body length, 2 SNPs-for withers height and rump height, 1 SNP-for body length and chest depth. Structural annotation of genomic regions covering a window of ±0.20 Mb showed the presence of 288 genes; 52 of them had the described functions in accordance with gene ontology. The main molecular functions of proteins encoded by these genes are the regulation of transcription, cell proliferation, angiogenesis, body growth, fatty acid and lipid metabolism, nervous system development, and spermatogenesis. SNPs common to body weight and localized within a window of ±200 kb from the structural genes CRADD, HMGA2, MSRB3, FUT8, MAX, and RAB15 were selected to create a test system. The study of meat productivity after slaughter and chemical analysis of muscle tissue in Karachai goats at the age of 8 months of different genotypes according to the identified SNPs revealed that rs268269710 is the most promising for further research and use in breeding. The GG genotype is associated with a larger live weight of animals, a larger carcass yield, the content of the boneless part in it, and the ratio of protein and adipose tissue in meat preferred for dietary nutrition. These results will contribute to the genetic improvement of Karachai goats.

MTaxi: A comparative tool for taxon identification of ultra low coverage ancient genomes

A major challenge in zooarchaeology is to morphologically distinguish closely related species' remains, especially using small bone fragments. Shotgun sequencing aDNA from archeological remains and comparative alignment to the candidate species' reference genomes will only apply when reference nuclear genomes of comparable quality are available, and may still fail when coverages are low. Here, we propose an alternative method, MTaxi, that uses highly accessible mitochondrial DNA (mtDNA) to distinguish between pairs of closely related species from ancient DNA sequences. MTaxi utilises mtDNA transversion-type substitutions between pairs of candidate species, assigns reads to either species, and performs a binomial test to determine the sample taxon. We tested MTaxi on sheep/goat and horse/donkey data, between which zooarchaeological classification can be challenging in ways that epitomise our case. The method performed efficiently on simulated ancient genomes down to 0.3x mitochondrial coverage for both sheep/goat and horse/donkey, with no false positives. Trials on n=18 ancient sheep/goat samples and n=10 horse/donkey samples of known species identity also yielded 100% accuracy. Overall, MTaxi provides a straightforward approach to classify closely related species that are difficult to distinguish through zooarchaeological methods using low coverage aDNA data, especially when similar quality reference genomes are unavailable. MTaxi is freely available at https://github.com/goztag/MTaxi.

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.

Ancient Mongolian aurochs genomes reveal sustained introgression and management in East Asia

Societies in East Asia have utilized domesticated cattle for over 5000 years, but the genetic history of cattle in East Asia remains understudied. Genome-wide analyses of 23 ancient Mongolian cattle reveal that East Asian aurochs and ancient East Asian taurine cattle are closely related, but neither are closely related to any modern East Asian breeds. We observe binary variation in aurochs diet throughout the early Neolithic, and genomic evidence shows millennia of sustained male-dominated introgression. We identify a unique connection between ancient Mongolian aurochs and the European Hereford breed. These results point to the likelihood of human management of aurochs in Northeast Asia prior to and during the initial adoption of taurine cattle pastoralism.

A missense mutation in RRM1 contributes to animal tameness

The increased tameness to reduce avoidance of human in wild animals has been long proposed as the key step of animal domestication. The tameness is a complex behavior trait and largely determined by genetic factors. However, the underlying genetic mutations remain vague and how they influence the animal behaviors is yet to be explored. Behavior tests of a wild-domestic hybrid goat population indicate the locus under strongest artificial selection during domestication may exert a huge effect on the flight distance. Within this locus, only one missense mutation RRM1I241V which was present in the early domestic goat ~6500 years ago. Genome editing of RRM1I241V in mice showed increased tameness and sociability and reduced anxiety. These behavioral changes induced by RRM1I241V were modulated by the alternation of activity of glutamatergic synapse and some other synapse-related pathways. This study established a link between RRM1I241V and tameness, demonstrating that the complex behavioral change can be achieved by mutations under strong selection during animal domestication.

Selection signatures for fiber production in commercial species: A review

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

Genomic analysis of indigenous goats in Southwest Asia reveals evidence of ancient adaptive introgression related to desert climate

Understanding how evolutionary pressures related to climate change have shaped the current genetic background of domestic animals is a fundamental pursuit of biology. Here, we generated whole-genome sequencing data from native goat populations in Iraq and Pakistan. Combined with previously published data on modern, ancient (Late Neolithic to Medieval periods), and wild Capra species worldwide, we explored the genetic population structure, ancestry components, and signatures of natural positive selection in native goat populations in Southwest Asia (SWA). Results revealed that the genetic structure of SWA goats was deeply influenced by gene flow from the eastern Mediterranean during the Chalcolithic period, which may reflect adaptation to gradual warming and aridity in the region. Furthermore, comparative genomic analysis revealed adaptive introgression of the KITLG locus from the Nubian ibex ( C. nubiana) into African and SWA goats. The frequency of the selected allele at this locus was significantly higher among goat populations located near northeastern Africa. These results provide new insights into the genetic composition and history of goat populations in the SWA region.

Neolithic hydroclimatic change and water resources exploitation in the Fertile Crescent

In the first millennia of the Holocene, human communities in the Fertile Crescent experienced drastic cultural and technological transformations that modified social and human-environments interactions, ultimately leading to the rise of complex societies. The potential influence of climate on this "Neolithic Revolution" has long been debated. Here we present a speleothem record from the Kurdistan Region of Iraq, covering from Early Neolithic to Early Chalcolithic periods (~ 11 to 7.3 ka, 9000-5300 BCE). The record reveals the influence of the Siberian High on regional precipitation, and shows large hydroclimatic variability at the multicentennial scale. In particular, it highlights wetter conditions between 9.7 and 9.0 ka, followed by an abrupt reduction of precipitation between 9.0 and 8.5 ka, and a wetter interval between 8.5 and 8.0 ka. A comparison with regional and local archaeological data demonstrates an influence of recorded hydroclimatic changes on settlement patterns (size, distribution, permanent vs. seasonal occupation) and on the exploitation of water resources by Neolithic to Chalcolithic populations. Our record does not show prominent hydroclimatic changes at 9.3 and 8.2 ka, thus not supporting direct influence of such rapid and widespread events on the process of Neolithization and its cultural dispersal.

Whole-genome selective scans detect genes associated with important phenotypic traits in goat (Capra hircus)

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. In this study, we focused on the worldwide goat breeds with outstanding traits and used whole-genome resequencing data in 361 samples from 68 breeds to detect genomic selection sweep regions. We identified 210-531 genomic regions with six phenotypic traits, respectively. Further gene annotation analysis revealed 332, 203, 164, 300, 205, and 145 candidate genes corresponding with dairy, wool, high prolificacy, poll, big ear, and white coat color traits. Some of these genes have been reported previously (e.g., KIT, KITLG, NBEA, RELL1, AHCY, and EDNRA), while we also discovered novel genes, such as STIM1, NRXN1, LEP, that may be associated with agronomic traits like poll and big ear morphology. Our study found a set of new genetic markers for genetic improvement in goats and provided novel insights into the genetic mechanisms of complex traits.

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.

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.

Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.

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.