ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure.

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

Selection of species specific panel of reference genes in peripheral blood mononuclear cells of native livestock species adapted to trans-Himalayan region of Leh-Ladakh.

The identification of appropriate references genes is an integral component of any gene expression-based study for getting accuracy and reliability in data interpretation. In this study, we evaluated the expression stability of 10 candidate reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) in peripheral blood mononuclear cells of livestock species that are adapted to high altitude hypoxia conditions of Leh-Ladakh. A total of 37 PBMCs samples from six native livestock species of Leh-Ladakh region such as Ladakhi cattle, Ladakhi yak, Ladakhi donkey, Chanthangi goat, Double hump cattle and Zanskar ponies were included in this study. The commonly used statistical algorithms such as geNorm, Normfinder, BestKeeper and RefFinder were employed to assess the stability of these RGs in all the livestock species. Our study has identified different panel of reference genes in each species; for example, EEF1A1, RPL4 in Ladakhi cattle; GAPDH, RPS9, ACTB in Ladakhi yak; HPRT1, B2M, ACTB in Ladakhi donkey; HPRT1, B2M, ACTB in Double hump camel, RPS9, HPRT1 in Changthangi goat, HPRT1 and ACTB in Zanskar ponies. To the best of our knowledge, this is the first systematic attempt to identify panel of RGs across different livestock species types adapted to high altitude hypoxia conditions. In future, the findings of the present study would be quite helpful in conducting any transcriptional studies to understand the molecular basis of high altitude adaptation of native livestock population of Leh-Ladakh.

Demographic pattern of A1/A2 beta casein variants indicates conservation of A2 type haplotype across native cattle breeds (Bos indicus) of India.

Genetic variations of the beta casein gene hold importance because of their probable association with human health. Comparative sequence analysis of ß-casein gene across Indian native, crossbred and exotic breeds in India revealed 15 SNPs and 4 INDELs corresponding to 14 haplotypes. The frequency of A2 type haplotype was maximum (0.941) across all Indian native breeds. Among the 15 variants reported for taurine breeds, only three (A1, A2 and B) were observed in analysed populations. Allelic profiling of A1/A2 ß-casein variants in ~ 4000 animals belonging to three cattle types and breeding bulls also revealed the predominance of A2 allele (0.95) in Indian cattle. The high proportion of A2 allele/haplotype indicates that Indian native cattle are the best suited to meet the demands for A2 milk globally. However, a higher percentage of heterozygous genotype (A1A2) in breeding bulls warrants the need to screen sire lines so as to drift the herd towards A2. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-022-03232-0.

Characterization of thermo-physiological, hematological, and molecular changes in response to seasonal variations in two tropically adapted native cattle breeds of Bos indicus lineage in hot arid ambience of Thar Desert.

The selection of climate resilient animal is necessary to secure the future of sustainable animal production. The present investigation therefore was an effort to unravel answers to the adaptation at physiological, hematological, and molecular levels in cows of hot arid region that helps them to survive harsh environment, to continue production and reproduction. This investigation was carried out in indicine cows over a period of one year, encompassing four seasons, wherein physiological data of 50 animals, hematological data of 15 animals, and gene expression profile of 5 animals from each of Sahiwal and Kankrej breeds per season was generated. In total, 5600 physiological observations, 1344 hematological observations, and 480 molecular samples were processed. The meteorological data revealed a high diurnal variation of temperature across seasons, with THI exceeding 80 during the months of summer and hot-humid seasons, indicating significant heat stress (HS). The physiological parameters showed an increasing trend with the incremental THI, with significantly (p < 0.05) higher values of rectal temperature (RT), respiration rate (RR), pulse rate (PR), and body surface temperature (BST) at ventral (VT), lateral (LT), dorsal (DT), and frontal (FT), in both breeds recorded during HS. The hematological pictures also revealed significant (p < 0.05) seasonal perturbations in erythrocytic and leucocytic parameters. Moreover, the molecular response was driven by a significant (p < 0.05) upregulation of all the key HSPs, HSP70, HSP90, HSP60, and HSP40, except HSP27 during the hotter months of summer and hot-humid seasons. The expression of HSF1, an important transcriptional regulator of  HSP70 was also significantly (p < 0.05) upregulated during summer season in both breeds. All the molecular chaperones revealed a significant upregulation during the summer season, followed by a decreasing trend by hot-humid season. The study indicated a well-developed thermotolerance mechanism in animals of both breeds, with Kankrej cows exhibiting better thermotolerance compared to Sahiwal cows.

Identification of Internal Reference Genes in Peripheral Blood Mononuclear Cells of Cattle Populations Adapted to Hot Arid Normoxia and Cold Arid Hypoxia Environments.

To estimate gene expression in a reliable manner, quantitative real-time polymerase chain reaction data require normalisation using a panel of stably expressed reference genes (RGs). To date, information on an appropriate panel of RGs in cattle populations reared at cold arid high-altitude hypoxia and hot arid tropical normoxia environments is not available. Therefore, the present study was carried out to identify a panel of stably expressed RGs from 10 candidate genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, HMBS, B2M, RPS15, and ACTB) in peripheral blood mononuclear cells (PBMCs) of cattle populations reared at cold arid high-altitude hypoxia and hot arid normoxia environments. Four different statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder were used to assess the stability of these genes. A total of 30 blood samples were collected: six adult heifers each of Ladakhi (LAC) and Holstein Frisian crosses (HFX) and 4 Jersey (JYC) cows from cold arid high-altitude hypoxia environments (group I) and five adult heifers each of Sahiwal (SAC), Karan Fries (KFC), and Holstein Friesian (HFC) cows from hot arid normoxia environments (group II). Combined analysis of group I and group II resulted in identification of a panel of RGs like RPS9, RPS15, and GAPDH that could act as a useful resource to unravel the accurate transcriptional profile of PBMCs from diverse cattle populations adapted to distinct altitudes.