Genomic patterns of selection in morphometric traits across diverse Indian cattle breeds.

This study seeks a comprehensive exploration of genome-wide selective processes impacting morphometric traits across diverse cattle breeds, utilizing an array of statistical methods. Morphometric traits, encompassing both qualitative and quantitative variables, play a pivotal role in characterizing and selecting livestock breeds based on their external appearance, size, and physical attributes. While qualitative traits, such as color, horn structure, and coat type, contribute to adaptive features and breed identification, quantitative traits like body weight and conformation measurements bear a closer correlation with production characteristics. This study employs advanced genotyping technologies, including the Illumina BovineSNP50 Bead Chip and next-generation sequencing methods like Reduced Representation sequencing, to identify genomic signatures associated with these traits. We applied four intra-population methods to find evidence of selection, such as Tajima's D, CLR, iHS, and ROH. We found a total of 40 genes under the selection signature, that were associated with morphometric traits in five cattle breeds (Kankrej, Tharparkar, Nelore, Sahiwal, and Gir). Crucial genes such as ADIPDQ, DPP6, INSIG1, SLC35D2 in Kankrej, LPL, ATP6V1B2, CDC14B in Tharparkar, HPSE2, PLAG1 in Nelore, PCSK1, PRKD1 in Sahiwal, and GNAQ, HPCAL1 in Gir were identified in our study. This approach provides valuable insights into the genetic basis of variations in body weight and conformation traits, facilitating informed selection processes and offering a deeper understanding of the evolutionary and domestication processes in diverse cattle breeds.

Relevance of pharmacogenetics and pharmacogenomics in veterinary clinical practice: A review.

The recent advances in high-throughput next-generation sequencing technologies have heralded the arrival of the Big Data era. As a result, the use of pharmacogenetics in drug discovery and individualized drug therapy has transformed the field of precision medicine. This paradigm shift in drug development programs has effectively reshaped the old drug development practices, which were primarily concerned with the physiological status of patients for drug development. Pharmacogenomics bridges the gap between pharmacodynamics and pharmacokinetics, advancing current diagnostic and treatment strategies and enabling personalized and targeted drug therapy. The primary goals of pharmacogenetic studies are to improve drug efficacy and minimize toxicities, to identify novel drug targets, to estimate drug dosage for personalized medicine, and to incorporate it as a routine diagnostic for disease susceptibility. Although pharmacogenetics has numerous applications in individualized drug therapy and drug development, it is in its infancy in veterinary medicine. The objective of this review is to present an overview of historical landmarks, current developments in various animal species, challenges and future perspectives of genomics in drug development and dosage optimization for individualized medicine in veterinary subjects.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses.

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data.

The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine 50K SNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milk production, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.

TXA2 mediates LPA1-stimulated uterine contraction in late pregnant mouse.

Lysophosphatidic acid (LPA) is known to increase uterine contraction in the estrus cycle and early pregnancy, however, the effect of LPA in late pregnant uterus and its mechanisms are not clear. In the present study, we show the LPA receptor subtypes expressed and the mechanism of LPA-induced contractions in late pregnant mouse uterus. We determined the relative mRNA expression of LPA receptor genes by quantitative PCR and elicited log concentration-response curves to oleoyl-L-α-LPA by performing tension experiments in the presence and absence of nonselective and selective receptor antagonists and inhibitors of the TXA2 pathway. LPA1 was the most highly expressed receptor subtype in the late pregnant mouse uterus and LPA1/2/3 agonist (Oleoyl-L-α LPA) elicited increased contractions in this tissue that had lesser efficacy compared to oxytocin. LPA1/3 antagonist, Ki-16425, and a potent LPA1 antagonist (AM-095) significantly inhibited the LPA-induced contractions. Further, the nonselective COX inhibitor, indomethacin, and potent thromboxane A2 synthase inhibitor, furegrelate significantly impaired LPA-induced contractions. Moreover, selective thromboxane receptor (TP) antagonist, SQ-29548, and Rho kinase inhibitor, Y-27632 almost eliminated LPA-induced uterine contractions. LPA1 stimulation elicits contractions in the late pregnant mouse uterus using the contractile prostanoid, TXA2 and may be targeted to induce labor in uterine dysfunctions/ dystocia.

Selection signatures for fiber production in commercial species: A review.

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Landmarks in the history of selective sweeps.

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

A whole genome scan reveals distinct features of selection in Zhaotong cattle of Yunnan province.

Over the years, indigenous cattle have not only played an essential role in securing primary food sources but have also been utilized for labor by humans, making them invaluable genetic resources. The Zhaotong cattle, a native Chinese breed from the Yunnan province, possess excellent meat quality and resistance to heat and humidity. Here we used whole genome sequencing data of 104 animals to delve into the population structure, genomic diversity and potential positive selection signals in Zhaotong cattle. The findings of this study demonstrate that the genetic composition of Zhaotong cattle was primarily derived from Chinese indicine cattle and East Asian cattle. The nucleotide diversity of Zhaotong cattle was only lower than that of Chinese indicine cattle, which was much higher than that of other taurine cattle. Genome-wide selection scans detected a series of positive candidate regions containing multiple key genes related to bone development and metabolism (CA10, GABRG3, GLDN and NOTUM), meat quality traits (ALG8, LINGO2, MYO5B, PRKG1 and GABRB1), immune response (ADA2, BMF, LEF1 and PAK6) and heat resistance (EIF2AK4 and LEF1). In summary, this study supplies essential genetic insights into the genome diversity within Zhaotong cattle and provides a foundational framework for comprehending the genetic basis of indigenous cattle breeds.

Molecular and phylogenetic analysis of MHC class I exons 7-8 in a variety of cattle and buffalo breeds.

The present study was conducted on the MHC class I (BoLA-A/BuLA-A) gene in Sahiwal, Jersey, Hariana, and Tharparkar breeds of cattle and Murrah, Mehsana, and Bhadawari breeds of buffalo to study the polymorphism. Exons 7-8 of the MHC class I gene was first characterized for polymorphism study in buffalo and the results reveal that this gene has a higher level of nucleotide changes than the cattle. Genes were investigated for polymorphisms in 285 animals of cattle and buffalo breeds. Molecular characterization of the MHC class I (BoLa-A/Bula-A) gene reveals a higher degree of polymorphism at the nucleotide level in cattle and buffalo. Results revealed this region has a higher level of polymorphisms in buffalo as campared to the cattle. Alul restriction patterns were monomorphic except for three different patterns but it was able to illustrate the differences in buffalo and cattle. SSCP analysis of exons 7-8 showed remarkable differences in cattle and buffalo. Sequence analysis revealed more closeness of Murrah breed with crossbred and indigenous cattle than Holstein Friesian. Exon 8 had more deletion and stop codon as compared to exon 7. The investigation confirmed that MHC class I BoLa-A/Bula-A exons 7-8 is highly polymorphic in buffalo as compared to cattle.

Characterization of BoLA class II DQA and DQB by PCR-RFLP, cloning, and sequencing reveals sequence diversity in crossbred cattle.

The BoLA class II DQA and DQB genes in crossbred cattle were studied using PCR-RFLP, cloning, and sequencing techniques. Seventy-two crossbred cattle (Vrindavani) were used in the current study. HaeIII and XbaI restriction enzymes digested DQA exon 2-3, revealing seven (HaeIII-A-G) and three (XbaI A-C) motifs, respectively. The BoLA-DQB gene was analyzed using PCR-RFLP with PstI and TaqI restriction enzymes, yielding five restriction motifs for each restriction enzyme (PstI-A-E and TaqI-A-E). In crossbred cattle, addition, deletion, and substitutions were observed in distinct sequences, resulting in variations in overall gene length. Changes in nucleotides at positions 64-80, 110-200, and 207-264 were largely responsible for polymorphism in DQA exon 2. The phylogenetic analysis predicted a high degree of nucleotide and amino acid changes in DQA exon 2-3 and DQB exon 2. DQA genes had a nucleotide dissimilarity of 0.3-25.4 percent, while DQB genes had a nucleotide dissimilarity of 1.5-14.3 percent. We cloned and sequenced 20 genotypes based on PCR-RFLP of the DQA and DQB genes. The current study observed variation in the DQA and DQB genes and will serve as a foundation for future research on the BoLA DQA and DQB genes.

Progress and future perspectives of livestock genomics in India: a mini review.

The field of genetics has evolved a lot after the emergence of molecular and advanced genomic technologies. The advent of Next Generation Sequencing, SNP genotyping platforms and simultaneous reduction in the cost of sequencing had opened the door to genomic research in farm animals. There are various applications of genomics in livestock, such as the use of genomic data: (i) to investigate genetic diversity and breed composition/population structure (ii) to identify genetic variants and QTLs related to economically important and ecological traits, genome-wide association studies (GWAS) and genomic signatures of selection; (iii) to enhance breeding programs by genomic selection. Compared to traditional methods, genomic selection is expected to improve selection response by increasing selection accuracy and reducing the generation interval due to early selection. Genomic selection (GS) in developed countries has led to rapid genetic gains, especially in dairy cattle, due to a well-established genetic evaluation system. Indian livestock system is still lagging behind developed nations in adopting these technologies. This review discusses the current status, challenges, and future perspectives of livestock genomics in India.

Evidence for selective sweeps in the MHC gene repertoire of various cattle breeds.

Major Histocompatibility Complex (MHC) genes are among the immune genes that have been extensively studied in vertebrates and are necessary for adaptive immunity. In the immunological response to infectious diseases, they play several significant roles. This research paper provides the selection signatures in the MHC region of the bovine genome as well as how certain genes related to innate immunity are undergoing a positive selective sweep. Here, we investigated signatures of historical selection on MHC genes in 15 different cattle populations and a total of 427 individuals. To identify the selection signatures, we have used three separate summary statistics. The findings show potential selection signatures in cattle from whom we isolated genes involved in the MHC. The most significant regions related to the bovine MHC are BOLA, non-classical MHC class I antigen (BOLA-NC1), Microneme protein 1 (MIC1) , Cluster of Differentiation 244 (CD244), Gap Junction Alpha-5 Protein (GJA5). It will be possible to gain new insight into immune system evolution by understanding the distinctive characteristics of MHC in cattle.

Genome-wide detection of copy number variations in Tharparkar cattle.

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

Comparative gene expression profile in circulating PBMCs of Bos indicus and crossbred cattle to understand disease tolerance mechanism.

The present investigation was performed to compare the global gene expression profile in peripheral blood mononuclear cells (PBMCs) of Bos indicus and crossbred (Bos taurus × B. indicus) cattle. Previously, several studies revealed the disease tolerance potential of B. indicus cattle but underlying genetic mechanism is still not fully explored. The PBMCs model was used for this investigation as it plays crucial role in the immune system regulation. Transcriptomic analysis revealed total 6767 significantly differentially expressed transcripts (fold change (absolute) >2.0, p < .05). In addition, 4149 transcripts were upregulated, 2618 transcripts were downregulated and fold change (absolute) of differentially expressed transcript varied from -223.32 to 213.63. Functional annotation analysis of differentially expressed genes confirmed their role in various molecular pathways viz. innate immune response, antigen processing and presentation, MHC protein complex, defense response to bacterium, regulation of immune response, positive regulation of JAK-STAT cascade, cytoskeletal protein binding, etc. Protein-protein interaction network analysis provided understanding of inter-relationship of immune genes with differentially expressed genes. In conclusion, this study could provide comprehensive information about the dysregulated genes and biological pathways in PBMCs which might be responsible for disease tolerance in B. indicus cattle.

Detection of genome-wide copy number variation in Murrah buffaloes.

Riverine Buffaloes, especially the Murrah breed because of their adaptability to harsh climatic conditions, is farmed in many countries to convert low-quality feed into valuable dairy products and meat. Here, we investigated the copy number variations (CNVs) in 296 Murrah buffalo using the Axiom® Buffalo Genotyping Array 90K (Affymetrix, Santa Clara, CA, USA). The CNVs were detected on the autosomes, using the Copy Number Analysis Module (CNAM) using the univariate analysis. 7937 CNVs were detected in 279 Buffaloes, the average length of the CNVs was 119,048.87 bp that ranged between 7800 and 4,561,030 bp. These CNVs were accounting for 10.33% of the buffalo genome, which was comparable to cattle, sheep, and goat CNV analyses. Further, CNVs were merged and 1541 CNVRs were detected using the Bedtools-mergeBed command. 485 genes were annotated within 196 CNVRs that were identified in at least 10 animals of Murrah population. Out of these, 40 CNVRs contained 59 different genes that were associated with 69 different traits. Overall, the study identified a significant number of CNVs and CNVRs in the Murrah breed of buffalo, with a wide range of lengths and frequencies across the autosomes. The identified CNVRs contained genes associated with important traits related to production and reproduction, making them potentially important targets for future breeding and genetic improvement efforts.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection.

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Genetic group and swim-up processing affect SYBR Green real-time PCR-based porcine semen sex ratio in Landrace and its crossbred boars.

Determination of factors affecting sex ratio is important while considering application of sex ratio enrichment approach. Present study aimed to design a SYBR Green qPCR-based method for measurement of primary sex ratio and to evaluate different factors (genetic group, sire, spermiogenic cycle and processing layer) affecting boar sperm sex ratio. The qPCR was based on relative copy number analysis of sex chromosome-specific single copy gene fragments with an autosomal gene as reference and was evaluated using DNA dilution series from pigs with numerically normal karyotype. The sex ratio was estimated from genomic DNA samples isolated from boar semen collected from different genetic groups at different time points and different processing layers. The X chromosome frequencies of semen samples revealed significant effect of genetic group. However, significant variation was observed neither within same genetic group nor between ejaculates of different spermatogenic cycles. Among the processing techniques studied, swim-up technique produced a significant X sperm enrichment in comparison to control whereas, Percoll density gradient failed to show any significant difference among layers. The lower layer in swim-up technique was found to contain higher proportion of X sperms. The designed qPCR is found to be an easy, less time-consuming method and does not require high end laboratory facilities or the specialized expertise. The lower layer of swim-up processing has a scope for X sperm enrichment in boar semen with proper validation.

Fat augments leptin-induced uterine contractions by decreasing JAK2 and BKCa channel expressions in late pregnant rats.

Altered lipid metabolism in obesity causes pregnancy complications in humans and animals. Leptin levels increase in pregnancy, as well as obesity. However, the effect of obesity on uterine leptin receptors and its distal signaling is not clear. The present study aimed to understand the effect of increased fat on leptin signaling in rat uterus. Wistar female rats were fed with an HF diet (40% Fat, 17% Sucrose, 1.25% Cholesterol, 0.75% Cholic acid) for 6 weeks before the mating and during pregnancy. HF diet significantly increased the fat depots, liver weight, serum, and tissue cholesterol levels. It produced fatty degeneration in the liver and caused infiltration of inflammatory cells, cystic endometrial glands, and sub endometrial fibrosis of the uterus. In isometric tension experiments, leptin caused a significant increase in uterine contractions in high fat-fed animals compared to control animals. Analysis of receptor expressions revealed no significant difference between the groups. However, a significant decrease in the JAK2 and BKCaα mRNA expression was observed in the uterus of high fat-fed rats. No change in the BKCaß, eNOS, iNOS, MLCP, and MLCK mRNA expressions was noticed in the HF group compared to the control. The findings of the present study suggest that the contractile response to leptin in the uterus of high fat-fed rats may be attributed to reduced signaling through JAK2 and, lowered expressions of BKCa channel α subunits.

Molecular characterization of CRBR2 fragment of TLR4 gene in association with mastitis in Vrindavani cattle.

The bovine TLR4 gene is an interesting candidate marker for mastitis resistance, since it is involved in neutrophil migration to and from the mammary gland during mastitis. TLR4 detects pathogen ligands, such as the Escherichia coli lipopolysaccharide (LPS) endotoxin and facilitates innate and adaptive immune responses. In the current study, a total of 130 crossbred cows (74 mastitis tolerant and 56 with clinical mastitis) kept at the Cattle and Buffalo Farm, IVRI, Izatnagar, were selected to explore the polymorphism in the co-receptor binding region 2 (CRBR2) fragment of the TLR4 gene. PCR-SSCP and sequence analysis showed two genotypes of the TLR4 gene's CRBR2 fragment, AA and AB, which were polymorphic in both the afflicted and tolerant groups. Sequencing revealed eight single nucleotide polymorphisms (SNPs) in allele A and ten SNPs in allele B. This genotype had no significant effect on the incidence of clinical mastitis according to the logistic regression model. Our study found insufficient evidence linking SNP variants in the CRBR2 region of the TLR4 gene to mastitis susceptibility in crossbred cattle.

Identification of important miRNAs in bubaline mammary gland at heifer stage - An in-silico approach.

The present study was undertaken to identify the microRNAs (miRNAs) expressed in the mammary tissue of a bubaline heifer. Small RNAs were isolated from the mammary gland tissue and enriched for miRNA fraction. The linker-ligated small RNAs were reverse transcribed to synthesize cDNA and amplified by PCR. The PCR products were ligated to the pGEM-T Easy vector; cloned into DH5 alpha cells and sequenced. Sequencing of 40 clones, randomly chosen from this library, produced 115 concatemerized short sequences. The short sequences were checked for their matches with the help of refseq_rna database, EST database (NCBI) and Ensembl. The analysis was performed for more than 90% identity with miRNAs across different species in miRBase. Alignment of putative small RNA sequences with the bovine genome was attempted in GenBank (NCBI) and Ensembl using BLAST. The small RNA sequences with a partial matches within the bovine genome and/or with flanking sequences (upstream or downstream) were analyzed for hairpin structures using the Mfold web server. Reverse complements were also assessed for the homology search. miRNA sequences showing only a partial match with already reported sequences were considered as a putative bubaline miRNAs. Six developmentally important putative miRNA precursors were identified from this study using cloning and sequencing followed by the Bioinformatics approach. This study will help in the elucidation of pathways involving miRNAs in bubaline species at the heifer stage.