Ancient Genomes Reveal the Evolutionary History and Origin of Cashmere-Producing Goats in China

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.

Identification of genomic characteristics and selective signals in Guizhou black goat

BACKGROUND

Guizhou black goat is one of the indigenous black goat breeds in the southwest region of Guizhou, China, which is an ordinary goat for mutton production. They are characterized by moderate body size, black coat, favorite meat quality with tender meat and lower odor, and tolerance for cold and crude feed. However, little is known about the genetic characteristics or variations underlying their important economic traits.

RESULTS

Here, we resequenced the whole genome of Guizhou black goat from 30 unrelated individuals breeding in the five core farms. A total of 9,835,610 SNPs were detected, and 2,178,818 SNPs were identified specifically in this breed. The population structure analysis revealed that Guizhou black goat shared a common ancestry with Shaanbei white cashmere goat (0.146), Yunshang black goat (0.103), Iran indigenous goat (0.054), and Moroccan goat (0.002). However, Guizhou black goat showed relatively higher genetic diversity and a lower level of linkage disequilibrium than the other seven goat breeds by the analysis of the nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity. Based on FST and θπ values, we identified 645, 813, and 804 selected regions between Guizhou black goat and Yunshang black goat, Iran indigenous goat, and cashmere goats. Combined with the results of XP-EHH, there were 286, 322, and 359 candidate genes, respectively. Functional annotation analysis revealed that these genes are potentially responsible for the immune response (e.g., CD28, CD274, IL1A, TLR2, and SLC25A31), humility-cold resistance (e.g., HBEGF, SOSTDC1, ARNT, COL4A1/2, and EP300), meat quality traits (e.g., CHUK, GAB2, PLAAT3, and EP300), growth (e.g., GAB2, DPYD, and CSF1), fertility (e.g., METTL15 and MEI1), and visual function (e.g., PANK2 and NMNAT2) in Guizhou black goat.

CONCLUSION

Our results indicated that Guizhou black goat had a high level of genomic diversity and a low level of linkage disequilibrium in the whole genome. Selection signatures were detected in the genomic regions that were mainly related to growth and development, meat quality, reproduction, disease resistance, and humidity-cold resistance in Guizhou black goat. These results would provide a basis for further resource protection and breeding improvement of this very local breed.

Animal Use Strategies in the Longshan Mountain Region of Northern China during the First Millennium BC: A Zooarchaeological Analysis of Yucun

The first millennium BC saw the expansion of the Western Zhou dynasty in its northwestern frontier, alongside the rise and development of the Qin State in the Longshan Mountain region of northern China. Exploring the subsistence practices of these communities is crucial to gaining a better understanding of the social, cultural, and political landscape in this region at the time. While much of the research to date has focused on the Qin people, the subsistence practices of the Zhou people remain poorly understood. In this study, we analyzed animal remains from Yucun, a large settlement site associated with the Zhou people, located to the east of the Longshan Mountain. These animal remains were recovered in the excavation seasons of 2018-2020. Our results show that pigs, dogs, cattle, caprines, and horses, which were the major domestic animals at Yucun, accounted for over 90.8% of the animal remains examined in terms of the number of identified specimens (NISP) and 72.8% in terms of the minimum number of individuals (MNI), with cattle and caprines playing dominant roles. In terms of the taxonomic composition and the mortality profiles of pigs, caprines, and cattle, Yucun shared similarities with Maojiaping and Xishan, two contemporaneous Qin cultural sites located to the west of the Longshan Mountain, and differ from other farming societies in the middle and lower reaches of the Yellow River valley. Considering the cultural attributes and topographic conditions of these various sites, these findings imply that environmental conditions may have played a more significant role than cultural factors in shaping the animal-related subsistence practices in northern China during the first millennium BC.

The multi-kingdom microbiome of the goat gastrointestinal tract

BACKGROUND

Goat is an important livestock worldwide, which plays an indispensable role in human life by providing meat, milk, fiber, and pelts. Despite recent significant advances in microbiome studies, a comprehensive survey on the goat microbiomes covering gastrointestinal tract (GIT) sites, developmental stages, feeding styles, and geographical factors is still unavailable. Here, we surveyed its multi-kingdom microbial communities using 497 samples from ten sites along the goat GIT.

RESULTS

We reconstructed a goat multi-kingdom microbiome catalog (GMMC) including 4004 bacterial, 71 archaeal, and 7204 viral genomes and annotated over 4,817,256 non-redundant protein-coding genes. We revealed patterns of feeding-driven microbial community dynamics along the goat GIT sites which were likely associated with gastrointestinal food digestion and absorption capabilities and disease risks, and identified an abundance of large intestine-enriched genera involved in plant fiber digestion. We quantified the effects of various factors affecting the distribution and abundance of methane-producing microbes including the GIT site, age, feeding style, and geography, and identified 68 virulent viruses targeting the methane producers via a comprehensive virus-bacterium/archaea interaction network.

CONCLUSIONS

Together, our GMMC catalog provides functional insights of the goat GIT microbiota through microbiome-host interactions and paves the way to microbial interventions for better goat and eco-environmental qualities. Video Abstract.

Ecogroups and maternal haplogroups reveal the ancestral origin of native Chinese goat populations based on the variation of mtDNA D-loop sequences

China is rich in goat breeding resources. Officially recognized local goat breeds are mainly distributed in agro-ecological regions. The population structure and matrilineal origin of native Chinese goats can be used to formulate protection and utilization strategies for these genetic resources. In this study, the genetic structure and maternal origin of native Chinese goats were investigated using mtDNA D-loop sequences. A total of 329 goat samples from 25 Chinese indigenous goat populations and five introduced goat breeds from abroad were collected; these populations were distributed in four ecogroups designated as Southwest, South-central, the North China Plain, and Foreign-ecogroup. A larger average number of nucleotide differences and richer nucleotide diversity were observed in South-central and Foreign-ecogroup, whereas these were lower in Southwest. The 216 haplotypes divided into several haplogroups, of which HapA contained 99 haplotypes distributed in Southwest, the North China Plain, and Foreign-ecogroup with high frequency (0.53-0.77), whereas the frequency of HapA in South-central was <0.09. HapB was mostly found in South-central (0.5538) and was distributed to the North China Plain (0.2667), while it was rare in Southwest (<0.08) and Foreign-ecogroup (<0.07). According to the estimation of kinship and ancestry, HapA had five ancestors (A2, A3, A5, A10, and A12), HapB had a single maternal ancestor (A8), and HapC had two maternal ancestors (A1 and A4). This study showed that native Chinese goat breeds were mainly divided into three haplogroups (HapA, HapB, and HapC) and goat populations have expanded in the ecological regions.

Population Genetic Analysis of Six Chinese Indigenous Pig Meta-Populations Based on Geographically Isolated Regions

The diversification of indigenous pig breeds in China has resulted from multiple climate, topographic, and human cultural influences. The numerous indigenous pig breeds can be geographically divided into six meta-populations; however, their genetic relationships, contributions to genetic diversity, and genetic signatures remain unclear. Whole-genome SNP data for 613 indigenous pigs from the six Chinese meta-populations were obtained and analyzed. Population genetic analyses confirmed significant genetic differentiation and a moderate mixture among the Chinese indigenous pig meta-populations. The North China (NC) meta-population had the largest contribution to genetic and allelic diversity. Evidence from selective sweep signatures revealed that genes related to fat deposition and heat stress response (EPAS1, NFE2L2, VPS13A, SPRY1, PLA2G4A, and UBE3D) were potentially involved in adaptations to cold and heat. These findings from population genetic analyses provide a better understanding of indigenous pig characteristics in different environments and a theoretical basis for future work on the conservation and breeding of Chinese indigenous pigs.

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.

Convergent Genomic Signatures of Cashmere Traits: Evidence for Natural and Artificial Selection

Convergent evolution provides powerful opportunities to investigate the genetic basis of complex traits. The Tibetan antelope (Pantholops hodgsonii) and Siberian ibex (Capra sibirica) belong to different subfamilies in Bovidae, but both have evolved similar superfine cashmere characteristics to meet the cold temperature in plateau environments. The cashmere traits of cashmere goats underwent strong artificial selection, and some traces of domestication also remained in the genome. Hence, we investigated the convergent genomic signatures of cashmere traits between natural and artificial selection. We compared the patterns of convergent molecular evolution between Tibetan antelope and Siberian ibex by testing positively selected genes, rapidly evolving genes and convergent amino acid substitutions. In addition, we analyzed the selected genomic features of cashmere goats under artificial selection using whole-genome resequencing data, and skin transcriptome data of cashmere goats were also used to focus on the genes involved in regulating cashmere traits. We found that molecular convergent events were very rare, but natural and artificial selection genes were convergent enriched in similar functional pathways (e.g., ECM-receptor interaction, focal adhesion, PI3K-Akt signaling pathway) in a variety of gene sets. Type IV collagen family genes (COL4A2, COL4A4, COL4A5, COL6A5, COL6A6) and integrin family genes (ITGA2, ITGA4, ITGA9, ITGB8) may be important candidate genes for cashmere formation and development. Our results provide a comprehensive approach and perspective for exploring cashmere traits and offer a valuable reference for subsequent in-depth research on the molecular mechanisms regulating cashmere development and fineness.

Factors driving the halophyte rhizosphere bacterial communities in coastal salt marshes

Introduction

Root-associated microorganisms promote plant growth and provide protection from stresses. Halophytes are the fundamental components maintaining ecosystem functions of coastal salt marshes; however, it is not clear how their microbiome are structured across large spatial scales. Here, we investigated the rhizosphere bacterial communities of typical coastal halophyte species (Phragmites australis and Suaeda salsa) in temperate and subtropical salt marshes across 1,100 km in eastern China.

Methods

The sampling sites were located from 30.33 to 40.90°N and 119.24 to 121.79°E across east China. A total of 36 plots were investigated in the Liaohe River Estuary, the Yellow River Estuary, Yancheng, and Hangzhou Bay in August 2020. We collected shoot, root, and rhizosphere soil samples. the number of pakchoi leaves, total fresh and dry weight of the seedlings was counted. The soil properties, plant functional traits, the genome sequencing, and metabolomics assay were detected.

Results

The results showed that soil nutrients (total organic carbon, dissolved organic carbon, total nitrogen, soluble sugars, and organic acids) are high in the temperate marsh, while root exudates (measured by metabolite expressions) are significantly higher in the subtropical marsh. We observed higher bacterial alpha diversity, more complex network structure, and more negative connections in the temperate salt marsh, which suggested intense competition among bacterial groups. Variation partitioning analysis showed that climatic, edaphic, and root exudates had the greatest effects on the bacteria in the salt marsh, especially for abundant and moderate subcommunities. Random forest modeling further confirmed this but showed that plant species had a limited effect.

Conclutions

Taken together, the results of this study revealed soil properties (chemical properties) and root exudates (metabolites) had the greatest influence on the bacterial community of salt marsh, especially for abundant and moderate taxa. Our results provided novel insights into the biogeography of halophyte microbiome in coastal wetlands and can be beneficial for policymakers in decision-making on the management of coastal wetlands.

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.

Genomic insights into the conservation and population genetics of two Chinese native goat breeds

Chinese goats are an important group of goats worldwide. However, there are few studies on the conservation priority, genetic relationship, and potential gene flow between Chinese and global goat breeds. Here, we genotyped 239 goats from conservation populations of the Chinese Guangfeng and Ganxi breeds using the GoatSNP50 BeadChip. The conservation priority, population structure, selection signatures and introgression of these goats were analyzed in the context of 36 global goat breeds. First, we showed that Guangfeng and Ganxi goats had the largest effective population sizes across the global breeds 13 generations ago. Nevertheless, Ganxi goats have recently experienced a high degree of inbreeding, resulting in their conservation priority based on total gene and allelic diversities being lower than that of most other Chinese breeds (including Guangfeng goats). Population structure and admixture analyses showed that an average of 18% of Guangfeng genomic components were introgressed from Boer goats approximately 18-yr ago. Next, we reconstructed the subfamily structure of the core populations of Guangfeng and Ganxi goats, and proposed reasonable conservation strategies for inbreeding management. Moreover, a list of candidate genes under selection for fertility, immunity, growth, and meat quality were detected in Guangfeng and Ganxi goats. Finally, we identified some genes related to body development and reproduction, which were introgressed from Boer goats and may be beneficial for improving performance and productivity of Guangfeng goats. In conclusion, this study not only provides new insights into the conservation and utilization of Guangfeng and Ganxi goats but also enriches our understanding of artificial introgression from exotic goats into Chinese local goats.

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

Simple Summary

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

Abstract

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

A deletion variant within the FGF5 gene in goats is associated with gene expression levels and cashmere growth

The FGF5 gene has been associated with the regulation of fibre length in mammals, including cashmere goats. A deletion variant at ~14 kb downstream of the FGF5 gene showed significant divergence between cashmere and non-cashmere goats in previous studies. In this study, we designed specific primers to genotype the deletion variant. The results of gel electrophoresis and Sanger sequencing revealed that a 507-bp deletion mutation is located at 95 454 685-95 455 191 of chromosome 6 in goats. Genotyping data from a large panel of 288 goats showed that the deletion at the FGF5 gene locus appeared to be associated with cashmere length. The deletion variant was close to fixation (frequency 0.97) in cashmere goats. Furthermore, electrophoretic mobility shift assays for evaluating DNA-protein interaction and mRNA expression levels of FGF5 suggested that the deletion variant may serve as a cis-acting element by specifically binding transcription factors to mediate quantitative changes in FGF5 mRNA expression. Our study illustrates how a structural mutation of the FGF5 gene has contributed to the cashmere growth phenotype in domestic goats. The deletion mutation within the FGF5 gene could potentially serve as a molecular marker of cashmere growth in cashmere goat breeding.

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.

Selection of Cashmere Fineness Functional Genes by Translatomics

Cashmere fineness is an important index to evaluate cashmere quality. Liaoning Cashmere Goat (LCG) has a large cashmere production and long cashmere fiber, but its fineness is not ideal. Therefore, it is important to find genes involved in cashmere fineness that can be used in future endeavors aiming to improve this phenotype. With the continuous advancement of research, the regulation of cashmere fineness has made new developments through high-throughput sequencing and genome-wide association analysis. It has been found that translatomics can identify genes associated with phenotypic traits. Through translatomic analysis, the skin tissue of LCG sample groups differing in cashmere fineness was sequenced by Ribo-seq. With these data, we identified 529 differentially expressed genes between the sample groups among the 27197 expressed genes. From these, 343 genes were upregulated in the fine LCG group in relation to the coarse LCG group, and 186 were downregulated in the same relationship. Through GO enrichment analysis and KEGG enrichment analysis of differential genes, the biological functions and pathways of differential genes can be found. In the GO enrichment analysis, 491 genes were significantly enriched, and the functional region was mainly in the extracellular region. In the KEGG enrichment analysis, the enrichment of the human papillomavirus infection pathway was seen the most. We found that the COL6A5 gene may affect cashmere fineness.

A single-nucleotide mutation within the TBX3 enhancer increased body size in Chinese horses

Chinese ponies are endemic to the mountainous areas of southwestern China and were first reported in the archaeological record at the Royal Tomb of Zhongshan King, Mancheng, dated to approximately ∼2,100 YBP.¹ Previous work has started uncovering the genetic basis of size variation in western ponies and horses, revealing a limited number of loci, including HMGA2,²LCORL/NCAPG,³ZFAT, and LASP1.^4,5 Whether the same genetic pathways also drive the small body size of Chinese ponies, which show striking anatomical differences to Shetland ponies,⁶ remains unclear.^2,7 To test this, we combined whole-genome sequences of 187 horses across China. Statistical analyses revealed top association between genetic variation at the T-box transcription factor 3 (TBX3) and the body size. Fine-scale analysis across an extended population of 189 ponies and 574 horses narrowed down the association to one A/G SNP at an enhancer region upstream of the TBX3 (ECA8:20,644,555, p = 2.34e−39). Luciferase assays confirmed the single-nucleotide G mutation upregulating TBX3 expression, and enhancer-knockout mice exhibited shorter limbs than wild-type littermates (p < 0.01). Re-analysis of ancient DNA data showed that the G allele, which is most frequent in modern horses, first occurred some ∼2,300 years ago and rose in frequency since. This supports selection for larger size in Asia from approximately the beginning of the Chinese Empire. Overall, this study characterized the causal regulatory mutation underlying small body size in Chinese ponies and revealed size as one of the main selection targets of past Chinese breeders.

•

One single A/G SNP in TBX3 enhancer region drives size variation in Chinese horses

•

The frequency of the G variant correlates positively with size in 763 horses

•

Cellular and mice models confirm it affects TBX3 transcription and the limb length

•

The G variant first occurred ∼2,300 years ago and rose in frequency since

RGD v2.0: a major update of the ruminant functional and evolutionary genomics database

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.

Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains

The Aceramic Neolithic (∼9600 to 7000 cal BC) period in the Zagros Mountains, western Iran, provides some of the earliest archaeological evidence of goat (Capra hircus) management and husbandry by circa 8200 cal BC, with detectable morphological change appearing ∼1,000 y later. To examine the genomic imprint of initial management and its implications for the goat domestication process, we analyzed 14 novel nuclear genomes (mean coverage 1.13X) and 32 mitochondrial (mtDNA) genomes (mean coverage 143X) from two such sites, Ganj Dareh and Tepe Abdul Hosein. These genomes show two distinct clusters: those with domestic affinity and a minority group with stronger wild affinity, indicating that managed goats were genetically distinct from wild goats at this early horizon. This genetic duality, the presence of long runs of homozygosity, shared ancestry with later Neolithic populations, a sex bias in archaeozoological remains, and demographic profiles from across all layers of Ganj Dareh support management of genetically domestic goat by circa 8200 cal BC, and represent the oldest to-this-date reported livestock genomes. In these sites a combination of high autosomal and mtDNA diversity, contrasting limited Y chromosomal lineage diversity, an absence of reported selection signatures for pigmentation, and the wild morphology of bone remains illustrates domestication as an extended process lacking a strong initial bottleneck, beginning with spatial control, demographic manipulation via biased male culling, captive breeding, and subsequently phenotypic and genomic selection.

GGVD: A goat genome variation database for tracking the dynamic evolutionary process of selective signatures and ancient introgressions

Understanding the evolutionary history and adaptive process depends on the knowledge that we can acquire from both ancient and modern genomic data. With the availability of a deluge of whole-genome sequencing data from ancient and modern goat samples, a user-friendly database making efficient reuse of these important resources is needed. Here, we use the genomes of 208 modern domestic goats, 24 bezoars, 46 wild ibexes, and 82 ancient goats to present a comprehensive goat genome variation database (GGVD). GGVD hosts a total of ∼41.44 million SNPs, ∼5.14 million indels, 6,193 selected loci, and 112 introgression regions. Users can freely visualize the frequency of genomic variations in geographical maps, selective sweeps in interactive tables, Manhattan plots, or line charts, as well as the heatmap patterns of the SNP genotype. Ancient data can be shown in haplotypes to track the state of genetic variants of selection and introgression events in the early, middle, and late stages. For facilitating access to sequence features, the UCSC Genome Browser, BLAT, BLAST, LiftOver, and pcadapt are also integrated into GGVD. GGVD will be a convenient tool for population genetic studies and molecular marker designing in goat breeding programs, and it is publicly available at http://animal.nwsuaf.edu.cn/GoatVar.

Deletions in the KAP6-1 gene are associated with fiber traits in cashmere-producing goats

Keratin-associated proteins (KAPs) are important structural components of fibers that predominantly present in the ortho-cortex. These proteins form a cross-linked network with keratin intermediate filaments (KIFs), thus producing a strong hair shaft. The keratin-associated protein 6-1 gene (KAP6-1) is a member of the KAPs family that has a potential correlation with fiber traits. In this study, we investigated the influence of KAP6-1 sequence polymorphisms on the fiber characteristics of a Chinese cashmere-producing goat breed (n = 844). Two main variants were found, including a three base pair (bp) deletion (namely B) and a 36-bp deletion (namely C), while the reference genotype of KAP6-1 was named A. Among them, the B variant was first reported on cashmere goats. This study then correlated these genotypes with the collected fiber data to investigate the potential association of these variants. The results showed that variant A is associated with decreased fiber diameter (p < 0.01), while variant C is associated with deceased fiber length (p < 0.01). These two related variants of the KAP6-1 gene have potential applications as gene-makers to improve the fiber diameter and length in cashmere-producing goats.

Whole-genome resequencing reveals diversity and selective signals in Longlin goat

The Longlin goat is one of the most valuable livestock species in Guangxi Autonomous Region of China, but its genomic diversity and selective signals are not clearly elucidated. Here we compared 20 genomes of Longlin goat to 66 genomes of other seven goat breeds worldwide to analyze patterns of Longlin goat genetic variation. We found the lowest linkage disequilibrium at the large distances between SNPs associated with the highest effective population size in the recent generations ago in Longlin goat. The eight goat breeds could be divided into Euro-African and East Asian goat population. Interestingly, like East Asian taurine, the same two migration phases might have occurred in the history of East Asian goat. More importantly, we identified selective signals implicated in immune resistance to disease, especially for skin disease, in Longlin goat. Our findings will not only help understand the evolutionary history and breed characteristic but can provide valuable resources for conservation of germplasm resources and implementation of crossbreeding programs.

Identification and population genetic analyses of copy number variations in six domestic goat breeds and Bezoar ibexes using next-generation sequencing

Background

Copy number variations (CNVs) are a major form of genetic variations and are involved in animal domestication and genetic adaptation to local environments. We investigated CNVs in the domestic goat (Capra hircus) using Illumina short-read sequencing data, by comparing our lab data for 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) to public data for 26 individuals from three other breeds (two Moroccan and one Chinese) and 21samples from Bezoar ibexes.

Results

We obtained a total of 2394 CNV regions (CNVRs) by merging 208,649 high-confidence CNVs, which spanned ~ 267 Mb of total length and accounted for 10.80% of the goat autosomal genome. Functional analyses showed that 2322 genes overlapping with the CNVRs were significantly enriched in 57 functional GO terms and KEGG pathways, most related to the nervous system, metabolic process, and reproduction system. Clustering patterns of all 85 samples generated separately from duplications and deletions were generally consistent with the results from SNPs, agreeing with the geographical origins of these goats. Based on genome-wide F_ST at each CNV locus, some genes overlapping with the highly divergent CNVs between domestic and wild goats were mainly enriched for several immunity-related pathways, whereas the genes overlapping with the highly differentiated CNVs between highland and lowland goats were mainly related to vitamin and lipid metabolism. Remarkably, a 507-bp deletion at ~ 14 kb downstream of FGF5 on chromosome 6 showed highly divergent (F_ST = 0.973) between the highland and lowland goats. Together with an enhancer activity of this sequence shown previously, the function of this duplication in regulating fiber growth deserved to be further investigated in detail.

Conclusion

We generated a comprehensive map of CNVs in goats. Many genetically differentiated CNVs among various goat populations might be associated with the population characteristics of domestic goat breeds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07267-6.