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.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

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

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

Genetic Architecture and Selection of Chinese Cattle Revealed by Whole Genome Resequencing

Genome-wide Assessment of Genetic Diversity and Population Structure Insights into Admixture and Introgression in Chinese Indigenous Cattle.

Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history

Genome-wide Scanning Reveals Genetic Diversity and Signatures of Selection in Chinese Indigenous Cattle Breeds

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

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

Species composition and environmental adaptation of indigenous Chinese cattle

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

Whole-genome Sequencing Reveals Mutational Landscape Underlying Phenotypic Differences Between Two Widespread Chinese Cattle Breeds

Genome-Wide SNPs and InDels Characteristics of Three Chinese Cattle Breeds

Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia

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

Matrilineal genetic inter-introgression of Bos taurus and Bos indicus in China

Genetic diversity and origin of Chinese cattle revealed by mtDNA D-loop sequence variation.

Population Genetic Diversity and Genetic Evaluation Models Reveal the Maternal Genetic Structure and Conservation Priority Characteristics of Indigenous Cattle in China

mtDNA Diversity and genetic lineages of eighteen cattle breeds from Bos taurus and Bos indicus in China

Genome-wide SNP discovery and evaluation of genetic diversity among six Chinese indigenous cattle breeds in Sichuan

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing.