Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing.
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.
Chen N
，Xia X
，Hanif Q
，Zhang F
，Dang R
，Huang B
，Lyu Y
，Luo X
，Zhang H
，Yan H
，Wang S
，Wang F
，Chen J
，Guan X
，Liu Y
，Li S
，Jin L
，Wang P
，Sun L
，Zhang J
，Liu J
，Qu K
，Cao Y
，Sun J
，Liao Y
，Xiao Z
，Cai M
，Mu L
，Siddiki AZ
，Asif M
，Mansoor S
，Babar ME
，Hussain T
，Silva GLLP
，Gorkhali NA
，Terefe E
，Belay G
，Tijjani A
，Zegeye T
，Gebre MG
，Ma Y
，Wang Y
，Huang Y
，Lan X
，Chen H
，Migliore NR
，Colombo G
，Semino O
，Achilli A
，Sinding MS
，Lenstra JA
，Cheng H
，Lu W
，Hanotte O
，Han J
，Jiang Y
，Lei C
... -
《Nature Communications》
Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan.
The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.
Iqbal N
，Liu X
，Yang T
，Huang Z
，Hanif Q
，Asif M
，Khan QM
，Mansoor S
... -
《PLoS One》
Novel polymorphisms in UTR and coding region of inducible heat shock protein 70.1 gene in tropically adapted Indian zebu cattle (Bos indicus) and riverine buffalo (Bubalus bubalis).
Due to evolutionary divergence, cattle (taurine, and indicine) and buffalo are speculated to have different responses to heat stress condition. Variation in candidate genes associated with a heat-shock response may provide an insight into the dissimilarity and suggest targets for intervention. The present work was undertaken to characterize one of the inducible heat shock protein genes promoter and coding regions in diverse breeds of Indian zebu cattle and buffaloes. The genomic DNA from a panel of 117 unrelated animals representing 14 diversified native cattle breeds and 6 buffalo breeds were utilized to determine the complete sequence and gene diversity of HSP70.1 gene. The coding region of HSP70.1 gene in Indian zebu cattle, Bos taurus and buffalo was similar in length (1,926 bp) encoding a HSP70 protein of 641 amino acids with a calculated molecular weight (Mw) of 70.26 kDa. However buffalo had a longer 5' and 3' untranslated region (UTR) of 204 and 293 nucleotides respectively, in comparison to Indian zebu cattle and Bos taurus wherein length of 5' and 3'-UTR was 172 and 286 nucleotides, respectively. The increased length of buffalo HSP70.1 gene compared to indicine and taurine gene was due to two insertions each in 5' and 3'-UTR. Comparative sequence analysis of cattle (taurine and indicine) and buffalo HSP70.1 gene revealed a total of 54 gene variations (50 SNPs and 4 INDELs) among the three species in the HSP70.1 gene. The minor allele frequencies of these nucleotide variations varied from 0.03 to 0.5 with an average of 0.26. Among the 14 B. indicus cattle breeds studied, a total of 19 polymorphic sites were identified: 4 in the 5'-UTR and 15 in the coding region (of these 2 were non-synonymous). Analysis among buffalo breeds revealed 15 SNPs throughout the gene: 6 at the 5' flanking region and 9 in the coding region. In bubaline 5'-UTR, 2 additional putative transcription factor binding sites (Elk-1 and C-Re1) were identified, other than three common sites (CP2, HSE and Pax-4) observed across all the analyzed animals. No polymorphism was found within the 3'-UTR of Indian cattle or buffalo as it was found to be monomorphic. The promoter sequences generated in 117 individuals showed a rich array of sequence elements known to be involved in transcription regulation. A total of 11 nucleotide changes were observed in the promoter sequence across the analyzed species, 3 of these changes were located within the potential transcription factor binding domains. We also identified 4 microsatellite markers within the buffalo HSP70.1 gene and 3 microsatellites within bovine HSP70.1. The present study identified several distinct changes across indicine, taurine and bubaline HSP70.1 genes that could further be evaluated as molecular markers for thermotolerance.
Sodhi M
，Mukesh M
，Kishore A
，Mishra BP
，Kataria RS
，Joshi BK
... -
《-》
Whole-Genome Analysis Deciphers Population Structure and Genetic Introgression Among Bovine Species.
It is known that throughout history and presently, taurine () and indicine/zebu () cattle were crossed with other bovine species (e.g., gayal, gaur, banteng, yak, wisent, and bison). Information on the role of interspecific hybridization to facilitate faster adaptation of the newly arrived domestic species to new environments is poorly known. Herein, we collected 266 samples of bovine species of the taurine, zebu, yak, and gaur from West Europe, Southwest China, Indian subcontinent, and Southeast Asia to conduct the principal component analysis (PCA), admixture, gene flow, and selection signature analyses by using SNPs distributed across the bovine autosomes. The results showed that the genetic relationships between the zebu, yak, and gaur mirrored their geographical origins. Three ancestral components of the European taurine, East Asian taurine, and Indian zebu were found in domestic cattle, and the bidirectional genetic introgression between the Diqing cattle and Zhongdian yak was also detected. Simultaneously, the introgressed genes from the Zhongdian yak to the Diqing cattle were mainly enriched with immune-related pathways, and the , , and genes related to the adaptation of high-altitude hypoxia were detected. Additionally, we found the genetic components of the Zhongdian yak had introgressed into Tibetan cattle. The 30 selected genes were detected in Tibetan cattle, which were significantly enriched in the chemokine signaling pathway. Interestingly, some genes (, , and ) associated with hypoxia response were discovered, in which and played important roles in angiogenesis and erythropoiesis, and heart function, respectively. This result showed that genetic introgression was one of the important ways for the environmental adaptation of domestic cattle.
Li R
，Chen S
，Li C
，Xiao H
，Costa V
，Bhuiyan MSA
，Baig M
，Beja-Pereira A
... -
《Frontiers in Genetics》
ResearchGate
(全网免费下载)
钛学术
(全网免费下载)
