Abstract: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.

Whole-genome comparative analysis reveals genetic mechanisms of disease resistance and heat tolerance of tropical Bos indicus cattle breeds

Comparative Analyses of Copy Number Variations Between Bos Taurus and Bos Indicus

Genetic basis of immunity in Indian cattle as revealed by comparative analysis of Bos genome

Genome-wide expression analysis reveals different heat shock responses in indigenous (Bos indicus) and crossbred (Bos indicus X Bos taurus) cattle

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

Transcriptome Analysis Reveals Potential Regulatory Genes Related To Heat Tolerance In Holstein Dairy Cattle

Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data

Genome-wide DNA methylation profiles regulate distinct heat stress response in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) cattle

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

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)

Global Genetic Diversity, Introgression, and Evolutionary Adaptation of Indicine Cattle Revealed by Whole Genome Sequencing.

Whole-Genome Resequencing Reveals Selection Signatures of Abigar Cattle for Local Adaptation

Transcriptomic profiling identifies breed-specific immune signatures of Tuberculosis susceptibility in cattle

Functionally prioritised whole-genome sequence variants improve the accuracy of genomic prediction for heat tolerance

Functional transcriptome analysis revealed upregulation of MAPK-SMAD signalling pathways in chronic heat stress in crossbred cattle

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

Whole-Genome Analysis Deciphers Population Structure and Genetic Introgression among Bovine Species.

Molecular markers for resistance against infectious diseases of economic importance

Advancing the Indian Cattle Pangenome: Characterizing Non-Reference Sequences in Bos indicus

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions.

Unraveling the Genomic Association for Milk Production Traits and Signatures of Selection of Cattle in a Harsh Tropical Environment