Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima’s D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Pasteurella multocida causes acute septicemic and respiratory diseases, including haemorrhagic septicaemia, in cattle and buffalo with case fatality of 100%. In the present study, mice were infected with P. multocida (1.6 × 10³ cfu, intraperitoneal) to evaluate host gene expression profile at early and late stages of infection using high throughput microarray transcriptome analyses. Several differentially expressed genes (DEGs) at both the time points were identified in P.multocida infected spleen, liver and lungs. Functional annotation of these DEGs showed enrichment of key pathways such as TLR, NF-κB, MAPK, TNF, JAK-STAT and NOD like receptor signaling pathways. Several DEGs overlapped across different KEGG pathways indicating a crosstalk between them. The predicted protein—protein interaction among these DEGs suggested, that the recognition of P. multocida LPS or outer membrane components by TLR4 and CD14, results in intracellular signaling via MyD88, IRAKs and/or TRAF6 leading to activation of NFκB and MAPK pathways and associated cytokines.

Highly pathogenic Avian influenza (HPAI) is a notifiable viral disease caused by avian influenza type A viruses of the Orthomyxoviridae family. Type A influenza genome consists of eight segments of negative-sense RNA. RNA segment 2 encodes three proteins, PB1, PB1-F2, and N40, which are translated from the same mRNA by ribosomal leaky scanning and reinitiation. Since these proteins are critical for viral replication and pathogenesis, targeting their expression can be one of the approaches to control and resist HPAI. MicroRNAs are short noncoding RNAs that regulate a variety of biological processes such as cell growth, tissue differentiation, apoptosis, and viral infection. In this study, a set of 300 miRNAs expressed in chicken lungs were screened against the HPAI virus (H5N1) segment 2 with different screening parameter like thermodynamic stability of heteroduplex, seed sequence complementarity, conserved target sequence, and target-site accessibility for identifying miRNAs that can potentially target the transcript of segment 2 of H5N1. Chicken miRNAs gga-mir-133c, gga-mir-1710, and gga-mir-146c* are predicted to target the expression of PB1, PB1-F2, and N40 proteins. This indicates that chicken has genetic potential to resist/tolerate H5N1 infection and these can be suitably exploited in designing strategies for control of avian influenza in chicken.

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.