Genome-wide runs of homozygosity signatures in diverse Indian goat breeds.

The present study analyzed ROH and consensus ROH regions in 102 animals of eleven diverse Indian goat (Capra hircus) breeds using whole genome sequencing. A total of 51,705 ROH and 21,271 consensus regions were identified. The mean number of ROH per animal was highest in the meat breed, Jharkhand Black (2693) and lowest in the pashmina breed, Changthangi (60). The average length of ROH (ALROH) was maximum in Kanniadu (974.11 Kb) and minimum in Tellicherry (146.98 Kb). Long ROH is typically associated with more recent inbreeding, whereas short ROH is connected to more ancient inbreeding. The overall ROH-based genomic inbreeding (FROH) was highest for Jharkhand Black (0.602) followed by Kanniadu (0.120) and Sangamneri (0.108) among all breeds. FROH of Jharkhand Black was higher than Kanniadu  up to 5 Mb ROH length category. However, in > 20 Mb ROH length category, Kanniadu (0.98)  exhibited significantly higher FROH than Jharkhand Black (0.46). This  implies that Kanniadu had higher levels of recent inbreeding than Jharkhand Black. Despite this, due to the presence of both recent and ancient inbreeding, Jharkhand Black  demonstrated higher overall FROH compared to Kanniadu. ROH patterns revealed dual purpose (meat and dairy) and pashmina breeds as less consanguineous while recent inbreeding was apparent in meat breeds. Analysis of ROH consensus regions identified selection sweeps in key genes governing intramuscular fat deposition, meat tenderisation, lean meat production and carcass weight (CDK4, ALOX15, CASP9, PRDM16, DVL1) in meat breeds; milk fat percentage and mammary gland development (POLD1, NOTCH2, ARHGAP35) in dual purpose (meat and dairy) breeds; while cold adaptation and hair follicle development (APOBEC1, DNAJC3, F2RL1, FGF9) in pashmina breed. MAPK, RAS, BMP and Wnt signaling pathways associated with hair follicle morphogenesis in Changthangi were also identified. PCA analysis based on ROH consensus regions revealed that meat breeds are more diverse than other goat breeds/populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-03921-y.

Exploring haplotype block structure, runs of homozygosity, and effective population size among dairy cattle breeds of India.

The present study aimed to explore haplotype structure, runs of homozygosity (ROH), effective population size and persistence of gametic phase among three indigenous dairy cattle breeds, viz., Sahiwal (n = 19), Tharparkar (n = 17), and Gir (n = 16) by using BovineHD single nucleotide polymorphism (SNP) genotyping assay. The filtered SNPs after quality control ranged from 44% in Sahiwal to 53% in Gir. The highest number of haplotype blocks was observed in Tharparkar (15,640) and the lowest in Sahiwal (8027) spanning 17.3% and 7.8% of genome, respectively. The average block length was found close to 26 kb which suggests that multiple recombination events fragmented the ancestral haplotypes into smaller sizes. Gir cattle had the largest number of runs of homozygosity (ROH) regions (1762) followed by Tharparkar (1528) and Sahiwal (1138). Without pedigree information, inbreeding coefficients estimated from ROH (FROH) revealed that Gir had the highest FROH (0.099) proposing more inbreeding rate in this population. Effective population size (Ne) decreased slowly over the last 60 generations and at 13 generations ago; Ne was estimated as 70 for all the three dairy breeds. The highest gametic phase correlation (r = 0.78) was observed for Sahiwal and Tharparkar breed pair suggesting formulation of multi-breed reference population for successful implementation of genomic selection among dairy breeds. The decline in effective population size among native Indian cattle breeds may help in formulating strategies for conservation and genetic improvement of native germplasm for future use.

Deciphering local adaptation of native Indian cattle (Bos indicus) breeds using landscape genomics and in-silico prediction of deleterious SNP effects on protein structure and function.

India has 50 registered breeds of native cattle (Bos indicus) which are locally adapted to diverse environmental conditions. This study aimed to investigate the genomic basis of adaptation of native Indian cattle and to predict the impact of key SNPs on the amino acid changes that affect protein function. The Illumina 777 K BovineHD BeadChip was used to genotype 178 native cattle belonging to contrasting landscapes and agro-climatic conditions. The genotype-environment association was investigated with R. SamBada, using 5,74,382 QC passed SNPs and 11 predictor variables (10 multi-collinearity controlled environmental variables and 1 variable as "score of PCA" on ancestry coefficients of individuals). In total, 1,12,780 models were selected as significant (q < 0.05) based on G score. The pathway ontology of the annotated genes revealed many important pathways and genes having a direct and indirect role in cold and hot adaptation. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Adaptive alleles identified in the present investigation might be responsible for the local adaptation of these cattle breeds. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03493-3.

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent.

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

Genome analyses revealed genetic admixture and selection signatures in Bos indicus.

The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (FST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.

Y-chromosome genetic diversity of Bos indicus cattle in close proximity to the centre of domestication.

Y-chromosome genetic diversity in and around its domestication origin and a better understanding of indicine-specific microsatellite alleles are imperative concerns but less -targeted. We analysed Y-chromosome markers in 301 bulls representing 19 native Indian cattle (Bos indicus) and identified new alleles and haplotypes. Compared to other indicine studies, the high Y-haplotype diversity found in Indian cattle supports the hypothesis of greater genetic variability across the centre of origin decreasing along migratory routes with increasing distance. Hence, a considerable paternal genetic diversity of Indian cattle appears to have been lost in transboundary commercial indicine breeds. The Khillar and Gir are the most diversified populations where the first tends to be the well-differentiated traditional breed carrying strikingly distinct Y-lineages with typical BM861-158 bp allele, characteristics of taurine cattle, while retaining standard indicine lineages for all other markers. Geographical distribution found to be an unreliable predictor of parental variation, and Y-lineages seemed closely related to Indian breed function/utility. The comprehensive Y-chromosome information will be useful to examine the demographic expansion/spread of Bos indicus lineages from close proximity to the domestication centre across different countries worldwide and such diversity should be preserved through effective management and conservation programs.

Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus.

Genome-wide runs of homozygosity (ROH) are suitable for understanding population history, calculating genomic inbreeding, deciphering genetic architecture of complex traits and diseases as well as identifying genes linked with agro-economic traits. Autozygosity and ROH islands, genomic regions with elevated ROH frequencies, were characterized in 112 animals of seven Indian native cattle breeds (B. indicus) using BovineHD BeadChip. In total, 4138 ROH were detected. The average number of ROH per animal was maximum in draft breed, Kangayam (63.62 ± 22.71) and minimum in dairy breed, Sahiwal (24.62 ± 11.03). The mean ROH length was maximum in Vechur (6.97 Mb) and minimum in Hariana (4.04 Mb). Kangayam revealed the highest ROH based inbreeding (FROH > 1Mb = 0.113 ± 0.059), whereas Hariana (FROH > 1Mb = 0.042 ± 0.031) and Sahiwal (FROH > 1Mb = 0.043 ± 0.048) showed the lowest. The high standard deviation observed in each breed highlights a considerable variability in autozygosity. Out of the total autozygous segments observed in each breed except Vechur, > 80% were of short length (< 8 Mb) and contributed almost 50% of the genome proportion under ROH. However, in Vechur cattle, long ROH contributed 75% of the genome proportion under ROH. ROH patterns revealed Hariana and Sahiwal breeds as less consanguineous, while recent inbreeding was apparent in Vechur. Maximum autozygosity observed in Kangayam is attributable to both recent and ancient inbreeding. The ROH islands were harbouring higher proportion of QTLs for production traits (20.68% vs. 14.64%; P≤ 0.05) but lower for reproductive traits (11.49% vs. 15.76%; P≤ 0.05) in dairy breeds compared to draft breed. In draft cattle, genes associated with resistant to diseases/higher immunity (LYZL1, SVIL, and GPX4) and stress tolerant (CCT4) were identified in ROH islands; while in dairy breeds, for milk production (PTGFR, CSN1S1, CSN2, CSN1S2, and CSN3). Significant difference in ROH islands among large and short statured breeds was observed at chromosome 3 and 5 involving genes like PTGFR and HMGA2 responsible for milk production and stature, respectively. PCA analysis on consensus ROH regions revealed distinct clustering of dairy, draft and short stature cattle breeds.

Genomic divergence reveals unique populations among Indian Yaks.

The present study focused upon identification of genome-wide SNPs through the reduced representation approach and to study the genomic divergence of the Indian yak populations. A total of 80 samples belonging to Arunachali yak (N = 20), Himachali yak (N = 20), Ladakhi yak (N = 20) and Sikkimi yak (N = 20) of India were used in the study. The results of the study revealed a total of 579575 high quality SNPs along with 50319 INDELs in the Indian yaks. The observed heterozygosity was found to be high in Himachali yak, followed by Arunachali yak, Ladakhi yak and Sikkimi yaks. The Sikkimi yaks was found to be genetically distant, followed by Ladakhi yaks which was observed to have some few individuals from Arunachali and Himachali yaks. Arunachali and Himachali yaks are found to get clustered together and are genetically similar. The study provides evidence about the genomic diversity in the Indian yak populations and information generated in the present study may help to formulate a suitable breeding plan for endangered Indian yaks. Moreover, the unique yak populations identified in the study will further help to focus attention for future characterization and prioritization of the animals for conservation purposes through the ddRAD approach.

Genome-wide discovery of SNPs in candidate genes related to production and fertility traits in Sahiwal cattle.

The present study was carried out to identify genome-wide genetic markers and variants in candidate genes for production and reproduction traits in Sahiwal cattle using a cost-effective reduced representation sequencing method. A total of 258,231 genome-wide SNPs were identified in Sahiwal cattle with reference to Bos indicus genome, of which 150,231 were novel SNPs. Among the high-confidence SNPs identified, 91.86% and 27.30% were genotyped in 50% and 100% of the samples. Mapping of the identified SNPs revealed 525 SNPs in candidate genes related to production traits while 333 SNPs were mapped to candidate genes related to reproduction traits. The SNPs identified in this study will facilitate further insights on tropical adaptation, domestication history and population structure of indigenous cattle. The variants in candidate genes identified in this study will serve as useful genetic tools, in the quest for phenotype modifying nucleotide change and help in designing appropriate genetic improvement programs.

Genomewide identification and annotation of SNPs in Bubalus bubalis.

The present study was carried out to identify and annotate the genome wide SNPs in Murrah buffalo genome. A total of 21.2 million raw reads from 4 pooled female Murrah buffalo samples were obtained using restriction enzyme digestion followed by sequencing with Illumina Hiseq 2000. After quality filtration, the reads were aligned to Murrah buffalo genome (ICAR-NBAGR) and Water buffalo genome (UMD_CASPUR_WB_2.0) which resulted in 99.37% and 99.67% of the reads aligning, respectively. A total of 130,688 high quality SNPs along with 35,110 indels were identified versus the Murrah bufffalo genome. Similarly 219,856 high quality SNPs along with 15,201 indels were identified versus the Water buffalo genome. We report 483 SNPs in 66 genes affecting Milk Production, 436 SNPs in 38 genes affecting fertility and 559 SNPs in 72 genes affecting other major traits. The average genome coverage was 13.4% and 14.8% versus the Murrah and Water buffalo genomes, respectively.

Evaluation of Bovine High-Density SNP Genotyping Array in Indigenous Dairy Cattle Breeds.

In total 52 samples of Sahiwal ( 19 ), Tharparkar ( 17 ), and Gir ( 16 ) were genotyped by using BovineHD SNP chip to analyze minor allele frequency (MAF), genetic diversity, and linkage disequilibrium among these cattle. The common SNPs of BovineHD and 54K SNP Chips were also extracted and evaluated for their performance. Only 40%-50% SNPs of these arrays was found informative for genetic analysis in these cattle breeds. The overall mean of MAF for SNPs of BovineHD SNPChip was 0.248 ± 0.006, 0.241 ± 0.007, and 0.242 ± 0.009 in Sahiwal, Tharparkar and Gir, respectively, while that for 54K SNPs was on lower side. The average Reynold's genetic distance between breeds ranged from 0.042 to 0.055 based on BovineHD Beadchip, and from 0.052 to 0.084 based on 54K SNP Chip. The estimates of genetic diversity based on HD and 54K chips were almost same and, hence, low density chip seems to be good enough to decipher genetic diversity of these cattle breeds. The linkage disequilibrium started decaying (r2 < 0.2) at 140 kb inter-marker distance and, hence, a 20K low density customized SNP array from HD chip could be designed for genomic selection in these cattle else the 54K Bead Chip as such will be useful.

Comparative screening of single nucleotide polymorphisms in ß-casein and κ-casein gene in different livestock breeds of India.

The most polymorphic milk protein gene is ß-casein; 13 protein variants are known in cattle. Milk protein genetic polymorphism has received considerable research interest in recent years because of possible associations between milk protein and economically important traits in livestock. The present study was undertaken to explore the genetic polymorphisms in exon 7 of ß-casein and exon 4 of κ-casein genes in Arunachali yaks (Bos grunniens), Sahiwal (Bos indicus) cattle, malpura sheep (Ovis aries) and Gaddi goat (Capra hircus). Results of the study revealed presence of 11 SNP variants in all livestock species. Four SNPs were observed in Bos indicus; two SNPs in Bos grunniens; three SNPs in Ovis aries and three SNPs in Capra hircus. These variations are found to be synonymous in nature as these variations do not result in their corresponding amino acids. A total of five polymorphic sites have been described at the κ-casein (CSN3) locus in the Indian domestic Gaddi goat (Capra hircus) when compared with exotic goat (X60763) while sequence analysis of κ-casein gene in sheep showed three novel nucleotide changes in malpura sheep when compared with the exotic sheep (AY237637). These results highlight the importance of taking into consideration the CSN3 SNPs when performing selection for milk composition in dairy livestock breeds.

Association of novel SNPs in the candidate genes affecting caprine milk fatty acids related to human health.

In the present investigation, 618 milk samples of Sirohi breed of goat were collected, and analyzed for conjugated linoleic acid (CLA, C18:2) and other fatty acids. The CLA in studied goat milk samples was 4.87 mg/g of milk fat and C18:2 cis-9, trans-11 contributes 2.9 mg/g of milk fat and trans10 cis12 contributes 0.82 mg/g of milk fat. The saturated fatty acids in the milk accounted for 69.55% and unsaturated fatty acid accounted for 28.50%. The unsaturated fatty acid was constituted by monounsaturated fatty acid (24.57%) and polyunsaturated fatty acids (3.96%.). The major contribution (45.56%) in total fatty acid was of C12:0, C14:0 and C16:0. C18:0 and short chain ones (C4:0, C6:0, C8:0, and C10:0) have a neutral or cholesterol-decreasing effect. The DNA sequence analysis of the genes (DGAT1, SCAP, PPARG, OLR, FABP3 and PRL) in a random panel of 8 Sirohi goats revealed 38 SNPs across the targeted regions. Out of the studied SNPs (38) across these genes, 22 SNPs had significant effect on one or a group of fatty acids including CLA. The genotypes at these loci showed significant differences in the least square means of a particular fatty acid or a group of fatty acids including CLA and its isomers.

Locus minimization in breed prediction using artificial neural network approach.

Molecular markers, viz. microsatellites and single nucleotide polymorphisms, have revolutionized breed identification through the use of small samples of biological tissue or germplasm, such as blood, carcass samples, embryos, ova and semen, that show no evident phenotype. Classical tools of molecular data analysis for breed identification have limitations, such as the unavailability of referral breed data, causing increased cost of collection each time, compromised computational accuracy and complexity of the methodology used. We report here the successful use of an artificial neural network (ANN) in background to decrease the cost of genotyping by locus minimization. The webserver is freely accessible (http://nabg.iasri.res.in/bisgoat) to the research community. We demonstrate that the machine learning (ANN) approach for breed identification is capable of multifold advantages such as locus minimization, leading to a drastic reduction in cost, and web availability of reference breed data, alleviating the need for repeated genotyping each time one investigates the identity of an unknown breed. To develop this model web implementation based on ANN, we used 51,850 samples of allelic data of microsatellite-marker-based DNA fingerprinting on 25 loci covering 22 registered goat breeds of India for training. Minimizing loci to up to nine loci through the use of a multilayer perceptron model, we achieved 96.63% training accuracy. This server can be an indispensable tool for identification of existing breeds and new synthetic commercial breeds, leading to protection of intellectual property in case of sovereignty and bio-piracy disputes. This server can be widely used as a model for cost reduction by locus minimization for various other flora and fauna in terms of variety, breed and/or line identification, especially in conservation and improvement programs.

Comparative sequence analysis in the exon 5 of growth hormone gene in the various livestock species of India.

The aim of the study was to identify genetic polymorphism in growth hormone (GH) gene locus of six different livestock species using PCR-Direct DNA sequencing method. In exon 5 of GH gene, 10 SNPs variants were identified in all livestock species studied, namely Bubalus bubalis, Bos indicus, Bos frontalis, Bos grunniens, Ovis aries, and Capra hircus. Four SNPs were observed in Bubalus bubalis, two SNPs in Bos indicus, one SNP in Ovis aries, and three SNPs in Capra hircus. No changes were observed in Bos grunniens and Bos frontalis when compared with the template sequence and the SNPs observed in the present investigation may be useful in the marker assisted selection.

Genetic structuring of nine Indian domestic goat breeds based on SNPs identified in IGF-1 gene.

The caprine Insulin like Growth Factor1 (IGF1) gene was analyzed for identification of single nucleotide polymorphisms (SNPs) and genetic structuring of Indian goat breeds. A panel of 80 samples belonging to nine Indian goat breeds (Capra hircus) including three large sized breeds (Jamunapari, Beetal and Jakhrana); three medium sized breeds (Sirohi, Barbari, and Osmanabadi) and three small sized breeds (Black Bengal, Changthangi, and Gaddi) were screened for SNP identification and diversity analysis. The comparative gene sequence analysis of all the nine goat breeds studied revealed a total of 18 SNPs in IGF1 gene. All the nucleotide changes were found to be synonymous. The mean observed heterozygosity was found to be maximum (0.074) in Sirohi, Beetal, Osmanabadi, and Gaddi breeds of goat, whereas it is found to be minimum (0.019) in Black Bengal breed of goat. The rest of the breeds were intermediate in terms of heterozygosity. The same has been confirmed by allele frequency distribution across the studied loci. Barbari and Gaddi were found to be more differentiated (0.0123), Changthangi and Jamunapari were least differentiated (0.00110) based on Nei's genetic distance.

Detection of polymorphism and sequence characterization of Toll-like receptor 7 gene of Indian goat revealing close relationship between ruminant species.

In this study, approximately 3.4 kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046 amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.

Management of attention deficit/hyperactivity disorder--use of an effective paradigm.

The management of ADHD poses a great problem before the psychologists, psychiatrists and different behavioural scientists. The multi-model approach, combining pharmacologic with different non-pharmacologic interventions, is more effective than any form of therapy. In the present study the theta feedback, a non-pharmacologic technique, has shown beneficial role among the low medicated ADHD cases. The present study also signifies the role of multi-model intervention in the management of ADHD.

Individuals at Risk of Coronary Heart Disease (CHD), its Prevention and Management by an Indigenous Compound.

A variety of rist factors have been suspected for causing the coronary heart disease. 406 cases of both sex groups with age range of 35 to 55 years were selected from three distinct localities of varanasi city. Individuals who reported single or more risk factors of CHD were isolated from the population of the particular areas. After a detailed preliminary screening of the subjects various physical, physiological, psychological and biochemical measurements were carried out. Other basis of initial findings, the cases who sowed abnormal lipid profile with dominant psychological involvement were given the organic extract of Inula racemosa (Pushkarmool), commiphora mukul (Guggulu), centella asiatica (Mandukaparni) and Hypericum perforatum (Basant in prescribed doses continuously for 6 moths. Correction in the lipid profile including triglycerides, blood pressure and the psychological factors like anxiety and depression to a significant level following least drugs treatment indicated the cardioprotective and therapeutic effects of the present formulation.Hence, by modifying the coronary risk factors, the incidence of CD can be minimized to a great extent as well as the test formulation may also be advocated in the prevention and management of CHD.