Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

Background: Chaka sheep, named after Chaka Salt Lake, are adapted to a harsh, highly saline environment. They are known for their high-grade meat quality and are a valuable genetic resource in China. Furthermore, the Chaka sheep breed has been designated a geographical symbol of agricultural products by the Chinese Ministry of Agriculture. Results: The genomes of 10 Chaka sheep were sequenced using next-generation sequencing, and compared to that of additional Chinese sheep breeds (Mongolian: Bayinbuluke and Tan; Tibetan: Oula sheep) to explore its population structure, genetic diversity and positive selection signatures. Principle component analysis and a neighbor-joining tree indicated that Chaka sheep significantly diverged from Bayinbuluke, Tan, and Oula sheep. Moreover, they were found to have descended from unique ancestors (K = 2 and K = 3) according to the structure analysis. The Chaka sheep genome demonstrated comparable genetic diversity from the other three breeds, as indicated by observed heterozygosity (Ho), expected heterozygosity (He), runs of homozygosity (ROH), linkage disequilibrium (LD) decay. The enrichment analysis revealed that in contrast to Mongolian or Tibetan lineage groups, the genes annotated by specific missense mutations of Chaka sheep were enriched with muscle structure development (GO:0061061) factors including insulin-like growth factor 1 (IGF1), growth differentiation factor 3 (GDF3), histone deacetylase 9 (HDAC9), transforming growth factor beta receptor 2 (TGFBR2), and calpain 3 (CAPN3), among others. A genome-wide scan using Fst and XP-CLR revealed a list of muscle-related genes, including neurofibromin 1 (NF1) and myomesin 1 (MYOM1), under potential selection in Chaka sheep compared with other breeds. Conclusions: The comprehensive genome-wide characterization provided the fundamental footprints for breeding and management of the Chaka sheep and confirmed that they harbor unique genetic resources.

Based on the previous GWAS research related to bovine heat tolerance trait, this study aimed to explore the effect of myosin-1a (MYO1A) gene on bovine heat tolerance trait, and find the molecular markers related to the heat tolerance of Chinese cattle. In our study, four novel candidate SNPs highly conserved in B. indicus breeds but barely existed in B. taurus were identified in MYO1A gene according to Bovine Genome Variation Database and Selective Signatures (BGVD). PCR and DNA sequencing were used to genotype 1072 individuals including 34 Chinese indigenous cattle breeds as well as Angus and Indian zebu. Two synonymous mutations (rs208210464 and rs110123931), one missense mutation (rs209999142; Phe172Ser), and one intron mutation (rs135771836) were detected. The frequencies of mutant alleles of the four SNPs gradually increased from northern groups to southern groups of Chinese cattle, which was consistent with the distribution of various climatic conditions of China. Additionally, four SNPs were significantly associated with four climatic conditions including annual mean temperature (T), relative humidity (H), temperature-humidity index (THI), and average annual sunshine hours (100-cloudiness) (SR). Among these, rs209999142 and Hap 1/1 had better performance than others. Our results suggested that rs209999142 was associated with heat-tolerance trait and rs208210464, rs110123931, and rs135771836 showed high phenotypic effect on heat-tolerance trait because of the strong linkage with rs209999142. These SNPs could be used as candidates for marker-assisted selection (MAS) in cattle breeding.

We report genome characterization of three native Chinese cattle breeds discovering ~34.3 M SNPs and ~3.8 M InDels using whole genome resequencing. On average, 10.4 M SNPs were shared amongst the three cattle breeds, whereas, 3.0 M, 4.9 M and 5.8 M were specific to LQ, WN and WS breeds, respectively. Gene ontology (GO)analysis revealed four immune response-related GO terms were over represented in all samples, while two immune signaling pathways were significantly over-represented in WS cattle. Altogether, we found immune related genes (PGLYRP2, ROMO1, FYB2, CD46, TSC1) in the three cattle breeds. Our study provides insights into the genetic basis of Chinese indicine adaptation to the tropic and subtropical environment, and provides a valuable resource for further investigations of genetic characteristics of the three breeds.

The Du’an goat is one of the most important farm animals in the Guangxi Autonomous Region of China, but the genetic basis underlying its adaptive traits has still not been investigated. Firstly, in this study, the genomes of 15 Du’an goats from a breeding farm were sequenced (mean depth: 9.50X) to analyze the patterns of genetic variation. A comparable diversity (17.3 million single nucleotide polymorphisms and 2.1 million indels) was observed to be associated with a lower runs of homozygosity-based inbreeding coefficient and smaller effective population size in comparison with other breeds. From selective sweep and gene set enrichment analyses, we revealed selective signals related to adaptive traits, including immune resistance (serpin cluster, INFGR1, TLR2, and immune-related pathways), body size (HMGA2, LCOR, ESR1, and cancer-related pathways) and heat tolerance (MTOR, ABCG2, PDE10A, and purine metabolism pathway). Our findings uncovered the unique diversity at the genomic level and will provide the opportunities for improvement of productivity in the Du’an goat.