ABSTRACT
Improving the low productivity levels of native cattle breeds in smallholder farming systems is a pressing concern in Pakistan. Crossbreeding high milk-yielding holstein friesian (HF) breed with the adaptability and heat tolerance of Sahiwal cattle has resulted in offspring that are well-suited to local conditions and exhibit improved milk yield. The exploration of how desirable traits in crossbred dairy cattle are selected has not yet been investigated. This study aims to provide the first overview of the selective pressures on the genome of crossbred dairy cattle in Pakistan. A total of eighty-one crossbred, thirty-two HF and twenty-four Sahiwal cattle were genotyped, and additional SNP genotype data for HF and Sahiwal were collected from a public database to equate the sample size in each group. Within-breed selection signatures in crossbreds were investigated using the integrated haplotype score. Crossbreds were also compared to each of their parental breeds to discover between-population signatures of selection using two approaches: cross-population extended haplotype homozygosity and fixation index. We identified several overlapping genes associated with production, immunity, and adaptation traits, including U6, TMEM41B, B4GALT7, 5S_rRNA, RBM27, POU4F3, NSD1, PRELID1, RGS14, SLC34A1, TMED9, B4GALT7, OR2AK3, OR2T16, OR2T60, OR2L3, and CTNNA1. Our results suggest that regions responsible for milk traits have generally experienced stronger selective pressure than others.
ABSTRACT
Linkage disequilibrium (LD) affects genomic studies accuracy. High-density genotyping platforms identify SNPs across animal genomes, increasing LD evaluation resolution for accurate analysis. This study aimed to evaluate the decay and magnitude of LD in a cohort of 81 crossbred dairy cattle using the GGP_HDv3_C Bead Chip. After quality control, 116,710 Single Nucleotide Polymorphisms (SNPs) across 2520.241 Mb of autosomes were retained. LD extent was assessed between autosomal SNPs within a 10 Mb range using the r2 statistics. LD value declined as inter-marker distance increased. The average r2 value was 0.24 for SNP pairs < 10 kb apart, decreasing to 0.13 for 50-100 kb distances. Minor allele frequency (MAF) and sample size significantly impact LD. Lower MAF thresholds result in smaller r2 values, while higher thresholds show increased r2 values. Additionally, smaller sample sizes exhibit higher average r2 values, especially for larger physical distance intervals (> 50 kb) between SNP pairs. Effective population size and inbreeding coefficient were 150 and 0.028 for the present generation, indicating a decrease in genetic diversity over time. These findings imply that the utilization of high-density SNP panels and customized/breed-specific SNP panels represent a highly favorable approach for conducting genome-wide association studies (GWAS) and implementing genomic selection (GS) in the Bos indicus cattle breeds, whose genomes are still largely unexplored. Furthermore, it is imperative to devise a meticulous breeding strategy tailored to each herd, aiming to enhance desired traits while simultaneously preserving genetic diversity.