High altitude represents some of the most extreme environments worldwide. The genetic changes underlying adaptation to such environments have been recently identified in multiple animals but remain unknown in horses. Here, we sequence the complete genome of 138 domestic horses encompassing a whole altitudinal range across China to uncover the genetic basis for adaptation to high-altitude hypoxia. Our genome data set includes 65 lowland animals across ten Chinese native breeds, 61 horses living at least 3,300 m above sea level across seven locations along Qinghai-Tibetan Plateau, as well as 7 Thoroughbred and 5 Przewalski’s horses added for comparison. We find that Tibetan horses do not descend from Przewalski’s horses but were most likely introduced from a distinct horse lineage, following the emergence of pastoral nomadism in Northwestern China ∼3,700 years ago. We identify that the endothelial PAS domain protein 1 gene (EPAS1, also HIF2A) shows the strongest signature for positive selection in the Tibetan horse genome. Two missense mutations at this locus appear strongly associated with blood physiological parameters facilitating blood circulation as well as oxygen transportation and consumption in hypoxic conditions. Functional validation through protein mutagenesis shows that these mutations increase EPAS1 stability and its hetero dimerization affinity to ARNT (HIF1B). Our study demonstrates that missense mutations in the EPAS1 gene provided key evolutionary molecular adaptation to Tibetan horses living in high-altitude hypoxic environments. It reveals possible targets for genomic selection programs aimed at increasing hypoxia tolerance in livestock and provides a textbook example of evolutionary convergence across independent mammal lineages.

Detection of selection footprints provides insight into the evolution process and the underlying mechanisms controlling the phenotypic diversity of traits that have been exposed to selection. Selection focused on certain characters, mapping certain genomic regions often shows a loss of genetic diversity with an increased level of homozygosity. Therefore, the runs of homozygosity (ROHs), homozygosity by descent (HBD), and effective population size (N_e) are effective tools for exploring the genetic diversity, understanding the demographic history, foretelling the signature of directional selection, and improving the breeding strategies to use and conserve genetic resources. We characterized the ROH, HBD, N_e, and signature of selection of six Chinese goat populations using single nucleotide polymorphism (SNP) 50K Illumina beadchips. Our results show an inverse relationship between the length and frequency of ROH. A long ROH length, higher level of inbreeding, long HBD segment, and smaller N_e in Guangfeng (GF) goats suggested intensive selection pressure and recent inbreeding in this breed. We identified six reproduction-related genes within the genomic regions with a high ROH frequency, of which two genes overlapped with a putative selection signature. The estimated pair-wise genetic differentiation (F_ST) among the populations is 9.60% and the inter- and intra-population molecular variations are 9.68% and 89.6%, respectively, indicating low to moderate genetic differentiation. Our selection signatures analysis revealed 54 loci harboring 86 putative candidate genes, with a strong signature of selection. Further analysis showed that several candidate genes, including MARF1, SYCP2, TMEM200C, SF1, ADCY1, and BMP5, are involved in goat fecundity. We identified 11 candidate genes by using cross-population extended haplotype homozygosity (XP-EHH) estimates, of which MARF1 and SF1 are under strong positive selection, as they are differentiated in high and low reproduction groups according to the three approaches used. Gene ontology enrichment analysis revealed that different biological pathways could be involved in the variation of fecundity in female goats. This study provides a new insight into the ROHs patterns for maintenance of within breed diversity and suggests a role of positive selection for genetic variation influencing fecundity in Chinese goat.

Long intergenic noncoding RNAs (lincRNAs) have been recognized in recent years as key regulators of biological processes. However, lincRNAs in goat remain poorly characterized both across various tissues and during different developmental stages in goat (Capra hircus). We performed the genome-wide discovery of the lincRNAs in goat by combining the RNA-seq dataset that were generated from 28 cashmere goat skin samples and the 12 datasets of goat tissues downloaded from the NCBI database. We identified a total of 5546 potential lincRNA transcripts that overlapped 3641 lincRNA genes. These lincRNAs exhibited a tissue-specific pattern. Specifically, there are 584 lincRNAs expressed exclusively in only one tissue, and 91 were highly expressed in hair follicle (HF). In addition, 2350 protein-coding genes and 492 lincRNAs were differentially expressed in the skin of goat. The majority exhibited the remarkable differential expression during the transition of the goat skin from the May–June to August–October time point, which covered the different seasons. Fundamental biological processes, such as skin development, were significantly enriched in these genes. Furthermore, we identified several lincRNAs highly expressed in the HF, which exhibited not only the co-expression pattern with the key factors to the HF development but also the activated expression in the August to October time point. Intriguingly, one of spatiotemporal lincRNAs, linc-chig1598 could be a potential regulator of distal-less homeobox 3 expression during the secondary hair follicle growth. This study will facilitate future studies aimed at unravelling the function of lincRNAs in hair follicle development.

Conservation of genetic resources is of great concern globally to maintain genetic diversity for sustainable food security. Comprehensive identification of the breed composition, estimation of inbreeding and effective population size are essential for the effective management of farm animal genetic resources and to prevent the animals from genetic erosion. The Zhongwei male (ZWM), Arbas Cashmere male (ACM) and Jining Grey male (JGM) goats are conserved in three different state goat farms in China but their family information, level of inbreeding and effective population size are unknown. We investigated the genomic relationship, inbreeding coefficient and effective population size in these three breeds from three state goat farms using the Illumina goat SNP50 BeadChip. Genomic relationships and phylogenetic analysis revealed that the breeds are clearly separated and formed separate clusters based on their genetic relationship. We obtained a high proportion of informative SNPs, ranging from 91.8% in the Arbas Cashmere male to 96.2% in the Jining Grey male goat breeds with an average mean of 96.8%. Inbreeding, as measured by F_ROH, ranged from 1.79% in ZWM to 8.62% in ACM goat populations. High F_ROH values, elevated genomic coverage of very long ROH (>30 Mb) and severe decline in effective population size were recorded in ACM goat farm. The existence of a high correlation between F_HOM and F_ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records. The N_e estimates 13 generations ago were 166, 69 and 79 for ZWM, ACM and JGM goat farm, respectively indicating that these goat breeds were strongly affected by selection pressure or genetic drift. This study provides insight into the genomic relationship, levels of inbreeding and effective population size in the studied goat populations conserved in the state goat farms which will be valuable in prioritizing populations for conservation and for developing suitable management practices for further genetic improvement of these Chinese male goats.

Tan sheep is a special breed of locally protected sheep in China, one of the best quality meat sheep in the world. Due to the unclear pedigree of the rams on the Ningxia Tan sheep breeding farm, we investigated 74 rams in the field and explored a new method for family division. Genomic DNA was extracted from the blood of breeding rams. Using Plink software, GCTA tools and R language, we analyzed the genetic structure, kinship, and inbreeding coefficient of the breeding sheep, which revealed the genetic relationship between the individuals. The results showed that there was no obvious clustering phenomenon in the PCA, and the genetic background of the samples was similar. The G matrix and IBS distance matrix indicated that most individuals were far away from each other. Paternity testing identified 24 pairs of unknown parent–child pairs, and all the Tan sheep could be divided into 12 families, which provided a reference for sheep breeding. The average inbreeding coefficient based on the ROH of this population was 0.049, so there was a low degree of inbreeding and the rams in the field were able to maintain high genetic diversity. Overall, we explored a more accurate method through paternity and kinship analysis; it provides a scientific basis for pedigree construction, which has an important application value for Tan sheep breeding.