Evaluation of non-synonym mutation in DGAT1 K232A as a marker for milk production traits in Ongole cattle and Murrah buffalo from Southern India

Various candidate genes have been reported to affect milk yield and composition in dairy cattle. A non-synonymous mutation in the DGAT1 gene, i.e., K232A was reported to have a strong association with milk yield and milk composition of Bos taurus. A study has been undertaken on 502 unrelated individuals belonging to indigenous Ongole cattle, crossbred cattle, and Murrah buffaloes from the Indian sub-continent with the objective to determine the polymorphism of the K232A locus and their association with milk yield and composition. Typing DGAT1 K232A allelic variation by PCR-RFLP using CfrI restriction enzyme revealed three genotypes in crossbred cattle. Genotype KK was more prevalent (0.60) in Jersey crossbred, whereas in Holstein Friesian crossbred it was KA genotype (0.48). In Ongole cattle and Murrah buffaloes, the locus did not exhibit polymorphism. The least-square mean of milk yields pooled over lactations across the DGAT1 variants was significantly (P < 0.05) higher among the homozygous (AA) genotypes, both in Jersey crossbred and HF crossbred cattle after adjusting for the effects of farm, parity, and season. The fat, SNF, and protein content values of AA genotypes were less than the KK genotypes in both the genetic groups (P > 0.05). The fixation of the DGAT1^K allele at the locus in Bos indicus cattle and Bubalus bubalis in the present study did not support its use as a reliable universal marker for milk production and composition traits.

Novel mutation in UMPS gene leads to false positive result of DUMPS (genetic disorder) in buffaloes: need for gene sequencing before confirming results of RFLP in new species

Deficiency of uridine monophosphate synthase (DUMPS) is a lethal genetic disorder associated with early embryonic mortality. Murrah and Mehsana male buffaloes (n = 594) were screened for DUMPS by PCR-RFLP technique. A few Murrah buffalo male calves were found to be carriers of DUMPS in RFLP, which has not been reported earlier. On the Sanger sequencing, a novel A to G substitution mutation was identified in AvaI restriction recognition site of UMPS gene in buffaloes. This mutation hinders digestion of DNA by AvaI which leds to false positive results for DUMPS carrier in RFLP. The results indicated that genome sequencing must be performed before confirming results of RFLP in any new species. All the buffaloes that were tested had only wild-type genotype in exon 5 for DUMPS specific allele.

Association of GH polymorphisms with growth traits in buffaloes

Members of the somatotrophic axis, especially GH and IGF1, are essential for growth. The association between GH polymorphisms and growth traits was numerously studied in cattle; however, no data are available for such association studies in buffalo. Therefore, this study was conducted to screen for polymorphisms in the GH gene and to study their putative association with growth traits in 200 Egyptian buffaloes. Polymerase chain reaction single-strand conformation polymorphism and sequencing were applied to look for polymorphisms in 3 loci spanning all exons and introns of buffalo GH. The C (MspI+) >T (MspI-) SNP in intron3, which is well known in cattle, was not detected in the examined buffaloes. However, 2 missense mutations were detected in exon5: one previously detected p.Leu153Val SNP, with very low frequencies for the mutant (Val) allele and one novel p.Asn174His SNP. At weaning age, the p.Leu153Val SNP was significantly associated with weaning body weight and gain with the positive effect of the wild allele (Leu) and higher GH serum, mRNA, and protein levels in animals with Leu/Asn and Leu/His haplotypes. At yearling age, the 2 SNPs associated with yearling weight and gain with positive effect for the mutant (Val and His) alleles with increased GH, and IGF1 serum, mRNA, and protein and GHR mRNA and protein levels in animals with Val/Asn haplotype. Therefore, the selection of Egyptian buffaloes with the Val/Asn haplotype could improve the growth traits of Egyptian buffaloes at yearling age which is the target age for perfect growing.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

Genomic Identification, Evolution and Sequence Analysis of the Heat-Shock Protein Gene Family in Buffalo

Heat-shock proteins (HSP) are conserved chaperones crucial for protein degradation, maturation, and refolding. These adenosine triphosphate dependent chaperones were classified based on their molecular mass that ranges between 10-100 kDA, including; HSP10, HSP40, HSP70, HSP90, HSPB1, HSPD, and HSPH1 family. HSPs are essential for cellular responses and imperative for protein homeostasis and survival under stress conditions. This study performed a computational analysis of the HSP protein family to better understand these proteins at the molecular level. Physiochemical properties, multiple sequence alignment, and phylogenetic analysis were performed for 64 HSP genes in the Bubalus bubalis genome. Four genes were identified as belonging to the HSP90 family, 10 to HSP70, 39 to HSP40, 8 to HSPB, one for each HSPD, HSPH1, and HSP10, respectively. The aliphatic index was higher for HSP90 and HSP70 as compared to the HSP40 family, indicating their greater thermostability. Grand Average of hydropathicity Index values indicated the hydrophilic nature of HSP90, HSP70, and HSP40. Multiple sequence alignment indicated the presence of highly conserved consensus sequences that are plausibly significant for the preservation of structural integrity of proteins. In addition, this study has expanded our current knowledge concerning the genetic diversity and phylogenetic relatedness of HSPs of buffalo with other mammalian species. The phylogenetic tree revealed that buffalo is more closely related to Capra hircus and distantly associated with Danio rerio. Our findings provide an understanding of HSPs in buffalo at the molecular level for the first time. This study highlights functionally important HSPs and indicates the need for further investigations to better understand the role and mechanism of HSPs.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Sequencing and assembly of the Egyptian buffalo genome

Water buffalo (Bubalus bubalis) is an important source of meat and milk in countries with relatively warm weather. Compared to the cattle genome, a little has been done to reveal its genome structure and genomic traits. This is due to the complications stemming from the large genome size, the complexity of the genome, and the high repetitive content. In this paper, we introduce a high-quality draft assembly of the Egyptian water buffalo genome. The Egyptian breed is used as a dual purpose animal (milk/meat). It is distinguished by its adaptability to the local environment, quality of feed changes, as well as its high resistance to diseases. The genome assembly of the Egyptian water buffalo has been achieved using a reference-based assembly workflow. Our workflow significantly reduced the computational complexity of the assembly process, and improved the assembly quality by integrating different public resources. We also compared our assembly to the currently available draft assemblies of water buffalo breeds. A total of 21,128 genes were identified in the produced assembly. A list of milk virgin-related genes; milk pregnancy-related genes; milk lactation-related genes; milk involution-related genes; and milk mastitis-related genes were identified in the assembly. Our results will significantly contribute to a better understanding of the genetics of the Egyptian water buffalo which will eventually support the ongoing breeding efforts and facilitate the future discovery of genes responsible for complex processes of dairy, meat production and disease resistance among other significant traits.

Genetic diversity, relatedness and inbreeding of ranched and fragmented Cape buffalo populations in southern Africa

Wildlife ranching, although not considered a conventional conservation system, provides a sustainable model for wildlife utilization and could be a source of valuable genetic material. However, increased fragmentation and intensive management may threaten the evolutionary potential and conservation value of species. Disease-free Cape buffalo (Syncerus caffer caffer) in southern Africa exist in populations with a variety of histories and management practices. We compared the genetic diversity of buffalo in national parks to private ranches and found that, except for Addo Elephant National Park, genetic diversity was high and statistically equivalent. We found that relatedness and inbreeding levels were not substantially different between ranched populations and those in national parks, indicating that breeding practices likely did not yet influence genetic diversity of buffalo on private ranches in this study. High genetic differentiation between South African protected areas highlighted their fragmented nature. Structure analysis revealed private ranches comprised three gene pools, with origins from Addo Elephant National Park, Kruger National Park and a third, unsampled gene pool. Based on these results, we recommend the Addo population be supplemented with disease-free Graspan and Mokala buffalo (of Kruger origin). We highlight the need for more research to characterize the genetic diversity and composition of ranched wildlife species, in conjunction with wildlife ranchers and conservation authorities, in order to evaluate the implications for management and conservation of these species across different systems.

MetaRNAseq analysis of surti buffalo rumen content reveals that transcriptionally active microorganisms need not be abundant

The present study describes rumen microbiota composition and their functional profiles in Indian Surti buffaloes by metagenomic (MG) and metatranscriptomic (MT) approaches. The study compares samples from buffaloes fed three different proportion of roughages; green and dry type of roughage; and different rumen liquor fractions. Irrespective of sample, Bacteroidetes and Firmicutes were the most predominant bacterial phyla, followed by Proteobacteria, Fibrobacteres and Actinobacteria while, Prevotella, Bacteroides, Ruminococcus and Clostridium were the most abundant genera. Different proportions of taxa were observed in both MG and MT approaches indicating the differences in organisms present and organisms active in the rumen. Higher proportions of fungal taxa were observed in MT while important organisms like Fibrobacter and Butyrivibrio and abundant organisms like Bacteroides and Prevotella were underrepresented in MT data. Functionally, higher proportions of genes involved in Carbohydrate metabolism, Amino acid metabolism and Translation were observed in both data. Genes involved in Metabolism were observed to be underrepresented in MT data while, those involved in Genetic information processing were overrepresented in MT data. Further, genes involved in Carbohydrate metabolism were overexpressed compared to genes involved in Amino acid metabolism in MT data compared to MG data which had higher proportion of genes involved in Amino acid metabolism than Carbohydrate metabolism. In all significant differences were observed between both approaches, different fractions of rumen liquor (liquid and solid) and different proportions of roughage in diet.

Analysis of meiotic segregation by triple-color fish on both total and motile sperm fractions in a t(1p;18) river buffalo bull

Chromosomal aberrations are relatively frequent pathologies in both humans and animals. Among them, translocations present a specific meiotic segregation pattern able to give a higher percentage of unbalanced gametes that can induce fertility problems. In this study, the meiotic segregation patterns of 1p, 1q and 18 Bubalus bubalis chromosomes were analyzed in both total sperm fraction and motile sperm fraction of a t(1p;18) carrier and a control bulls by triple-color FISH analysis with a pool of specific BAC probes. The frequencies of each total sperm fraction products in the carrier resulting from alternate, adjacent I, adjacent II and 3:1 segregation were 39%, 20%, 1% and 38%, respectively. On the other hand, the frequencies of each motile sperm fraction products in the carrier resulting from alternate, adjacent I, adjacent II and 3:1 segregation were 93%, 5%, 0% and 2%, respectively. The frequencies of normal sperms in the carrier were 27% and 69% in total sperm fraction and motile sperm fraction, respectively. The frequencies detected in motile sperm fraction were also validated by comparison with bull’s progeny. To our knowledge, this is the first report on the meiotic segregation patterns in motile sperm fractions of B. bubalis bull carrying a chromosomal translocation. These data suggest that translocation has a very limited effect on aneuploidy in the gametes, and therefore, on the reproductive abilities of the bull.

miR-326 down-regulate CYP19A1 expression and estradiol-17b production in buffalo granulosa cells through CREB and C/EBP-β

Ovarian granulosa cells, known to be endocrine cells, have well active TLR4-/NFKB signalling mediated innate immune capabilities. We have previously shown that endotoxin not only transiently regulates proinflammatory cytokines but cells become tolerant on repeated exposure to endotoxin and impaired granulosa cells functions, which includes downregulation of CYP19A1 gene. To understand further endotoxin tolerance and impaired granulosa cells function, genome-wide transcriptomic profiling in endotoxin tolerant buffalo granulosa cells (bGCs) identified miR-326 as upregulated amongst top 5 DE miRNAs [unpublished data] and qPCR validation confirmed its upregulation during endotoxin tolerance. In silico analyses showed that miR-326 targets CYP19A1 gene. Therefore, in the present study, we elucidated the role of miR-326 in buffalo granulosa cells (bGCs). We first validated its expression vis-à-vis CYP19A1 gene expression in bGCs, both in vivo and in vitro. Results showed an inverse relationship between miR-326 and CYP19A1 expression. Similarly, transcription factors, known to be involved in CYP19A1 gene regulation, CREB and C/EBP-β expression was also found to be decreased in granulosa cells mimicking pre-ovulatory follicular stage. Further, miR-326 mimic was transfected to bGCs in culture and expression of CYP19A1 and CREB & C/EBP-β and genes encoding other enzymes of steroidogenesis pathway were also analyzed. The present study results showed that miR-326 significantly inhibits the expression of CYP19A1 gene while expression of transcription factors CREB and C/EBP-β was found to be upregulated. The expression of STAR and CYP11A1 was found to be unaffected. To elucidate the molecular mechanism of miR-326 mediated downregulation of CYP19A1, binding analyses of RNA polymerase II and CEBP-β to CYP19A1 gene promoter II was analyzed. The result also showed decreased binding of RNA polymerase II with increased binding of CEBP-β to CYP19A1 gene promoter II in bGCs, transfected with miR-326 as compared to control. In summary, our results suggest that miR-326 upregulate CREB and CREB may activate C/EBP-β and later inhibited the transcription of CYP19A1 and decreased estradiol-17b production. The miR-326 mediated down-regulation of the CYP19A1 gene involving CREB-C/EBP-β can be exploited in developing strategies to attenuate endotoxin-mediated tolerance induced impaired granulosa cells function to ensure proper fertility in females.

Potential authentication of various meat-based products using simple and efficient DNA extraction method

BACKGROUND

The growth of halal food consumption worldwide has resulted in an increase in the request for halal authentication. DNA-based detection using powerful real-time polymerase chain reaction (PCR) technique has been shown to be highly specific and sensitive authentication tool. The efficient DNA extraction method in terms of quality and quantity is a backbone step to obtain successful real-time PCR assays. In this study, different DNA extraction methods using three lysis buffers were evaluated and developed to recommend a much more efficient method as well as achieve a successful detection using real-time PCR.

RESULTS

The lysis buffer 2 (LB2) has been shown to be the best lysis buffer for DNA extraction from both raw and processed meat samples comparing to other lysis buffers tested. Hence, the LB2 has been found to be ideal to detect meat and porcine DNAs by real-time PCR using pairs of porcine specific primers and universal primers which amplified at 119 bp fragment and 93 bp fragment, respectively. This assay allows detection as low as 0.0001 ng of DNA. Higher efficiency and sensitivity of real-time PCR via a simplified DNA extraction method using LB2 have been observed, as well as a reproducible and high correlation coefficient (R²  = 0.9979) based on the regression analysis of the standard curve have been obtained.

CONCLUSION

This study has established a fast, simple, inexpensive and efficient DNA extraction method that is feasible for raw and processed meat products. This extraction technique allows an accurate DNA detection by real-time PCR and can also be implemented to assist the halal authentication of various meat-based products available in the market. © 2019 Society of Chemical Industry.

Investigating genetic variability in Hsp70 gene-5'-fragment and its association with thermotolerance in Murrah buffalo (Bubalus bubalis) under sub-tropical climate of India

The present study was undertaken to investigate genetic variability in a fragment comprising 5'UTR along with partial coding sequence of Hsp70 gene and its association with thermotolerance traits in Murrah buffalo at ICAR-Research Complex for Eastern Region, Patna (India). The allelic variants were identified from genomic DNA samples using SSCP technique. The PCR products were sequenced and analyzed. Data on different thermotolerance traits recorded in three seasons were analyzed by least squares ANOVA taking the SSCP genotypes as fixed effect. Two allelic variants (A and B), each of 503-bp in size, were documented with frequency of 0.59 and 0.41, respectively, and three genotypes (AA, AB and BB) with corresponding frequency of 0.30, 0.58 and 0.12. The allelic variants were due to single nucleotide substitution at 55th base position leading to a change of threonine (A) to methionine (B) in amino acid sequence. Both the allelic variants had 99.8% similarity in nucleotide sequence. In phylogenetic tree, allele A was in a cluster while allele B and Gangatiri cattle sequence formed a different cluster. The SSCP genotypes had significant effect on different thermotolerance traits in summer with thermo-humidity index of ≥ 84. Buffaloes with AA genotype had the highest (P ˂ 0.05) summer evening rectal temperature, respiration rate and pulse rate, inferring that the buffaloes carrying AA genotype had more stress in summer than those with AB and BB genotype. These SSCP genotypes might have differential role in heat shock protein response to induce thermotolerance of Murrah buffaloes in Gangetic plains.

Genomic analyses reveal distinct genetic architectures and selective pressures in buffaloes

BACKGROUND

The domestic buffalo (Bubalus bubalis) is an essential farm animal in tropical and subtropical regions, whose genomic diversity is yet to be fully discovered.

RESULTS

In this study, we describe the demographic events and selective pressures of buffalo by analyzing 121 whole genomes (98 newly reported) from 25 swamp and river buffalo breeds. Both uniparental and biparental markers were investigated to provide the final scenario. The ancestors of swamp and river buffalo diverged ∼0.23 million years ago and then experienced independent demographic histories. They were domesticated in different regions, the swamp buffalo at the border between southwest China and southeast Asia, while the river buffalo in south Asia. The domestic stocks migrated to other regions and further differentiated, as testified by (at least) 2 ancestral components identified in each subspecies. Different signals of selective pressures were also detected in these 2 types of buffalo. The swamp buffalo, historically used as a draft animal, shows selection signatures in genes associated with the nervous system, while in river dairy breeds, genes under selection are related to heat stress and immunity.

CONCLUSIONS

Our findings substantially expand the catalogue of genetic variants in buffalo and reveal new insights into the evolutionary history and distinct selective pressures in river and swamp buffalo.

Expression and functional role of bone morphogenetic proteins (BMPs) in placenta during different stages of pregnancy in water buffalo (Bubalus bubalis)

The objective of this study was to document the expression and functional role of BMPs in the placental (caruncle; CAR, cotyledon; COT) during different stages of pregnancy in water buffalo. Samples collected from Early pregnancy 1 (EP1); Early pregnancy 2 (EP2), Mid pregnancy (MP), Late pregnancy (LP) while the third stage of oestrus cycle (NP) was taken as control. Also, the synergistic role of BMP4/BMP7 or combination on mRNA expression of vWF, PCNA, StAR, CYP11A1, 3βHSD, and BAX were studied in trophoblast cells cultured (TCC) during an early stage. The qPCR and immunoblotting studies revealed that BMP2, BMPR1A, BMPR1B, and BMPR2 mRNA level was significantly (p < 0.05) upregulated during early pregnancy in COTs while in CARs it was significantly upregulated (p < 0.05) during all the stages of pregnancy.BMP4 mRNA level was significantly upregulated (p < 0.05) during early pregnancy in COTs as well as in CARs. BMP6 expression was significantly upregulated (p < 0.05) during early and late stages of pregnancy. BMP7 mRNA level was upregulated (p < 0.05) during the late stage of pregnancy in COTs. At 100 ng/ml, the BMP4 maximally stimulated the transcripts of StAR, CYP11A1, and 3βHSD while BMP7 maximally stimulated the transcripts of 3βHSD that paralleled with P4 accretion in the media (P < 0.05). BMP4 as well as BMP7 upregulated the transcripts of PCNA, vWF, and downregulated BAX in the TCC (P < 0.05). In conclusion, BMPs are expressed in a regulated manner with stage-specific differences in the placenta and promotes the angiogenesis, proliferation, cell survivability, and steroidogenesis thereby regulating placental function in an autocrine/paracrine manner in water buffalo.

Unusual interplay of contrasting selective pressures on β-defensin genes implicated in male fertility of the Buffalo (Bubalus bubalis)

BACKGROUND

The buffalo, despite its superior milk-producing ability, suffers from reproductive limitations that constrain its lifetime productivity. Male sub-fertility, manifested as low conception rates (CRs), is a major concern in buffaloes. The epididymal sperm surface-binding proteins which participate in the sperm surface remodelling (SSR) events affect the survival and performance of the spermatozoa in the female reproductive tract (FRT). A mutation in an epididymal secreted protein, beta-defensin 126 (DEFB-126/BD-126), a class-A beta-defensin (CA-BD), resulted in decreased CRs in human cohorts across the globe. To better understand the role of CA-BDs in buffalo reproduction, this study aimed to identify the BD genes for characterization of the selection pressure(s) acting on them, and to identify the most abundant CA-BD transcript in the buffalo male reproductive tract (MRT) for predicting its reproductive functional significance.

RESULTS

Despite the low protein sequence homology with their orthologs, the CA-BDs have maintained the molecular framework and the structural core vital to their biological functions. Their coding-sequences in ruminants revealed evidence of pervasive purifying and episodic diversifying selection pressures. The buffalo CA-BD genes were expressed in the major reproductive and non-reproductive tissues exhibiting spatial variations. The Buffalo BD-129 (BuBD-129) was the most abundant and the longest CA-BD in the distal-MRT segments and was predicted to be heavily O-glycosylated.

CONCLUSIONS

The maintenance of the structural core, despite the sequence divergence, indicated the conservation of the molecular functions of the CA-BDs. The expression of the buffalo CA-BDs in both the distal-MRT segments and non-reproductive tissues indicate the retention the primordial microbicidal activity, which was also predicted by in silico sequence analyses. However, the observed spatial variations in their expression across the MRT hint at their region-specific roles. Their comparison across mammalian species revealed a pattern in which the various CA-BDs appeared to follow dissimilar evolutionary paths. This pattern appears to maintain only the highly efficacious CA-BD alleles and diversify their functional repertoire in the ruminants. Our preliminary results and analyses indicated that BuBD-129 could be the functional ortholog of the primate DEFB-126. Further studies are warranted to assess its molecular functions to elucidate its role in immunity, reproduction and fertility.

Muskox myoglobin: purification, characterization and kinetics studies compared with cattle and water buffalo myoglobins

BACKGROUND

The Arctic muskox has economic potential as an alternative meat species and is becoming increasingly popular. The present study aimed to determine the primary structure and pseudoperoxidase activity of muskox myoglobin (Mb) compared to cattle and water buffalo myoglobins.

RESULTS

The primary structure of muskox Mb was determined via a matrix-assisted laser desorption ionization-time of flight mass spectrometry-based mapping approach using the sheep Mb as a reference sequence. The muskox Mb consists of 153 amino acid residues and shows 100% identity with sheep Mb, whereas 98.69% and 97.38% identity is found with cattle and water buffalo Mbs, respectively. Muskox Mb has an autoxidation rate (MetMb formation) higher than both cattle and water buffalo Mbs at pH 7.2 (37 °C). Moreover, its pseudoperoxidase activity is higher than both cattle and water buffalo Mbs at pH 7.4 (physiological pH), whereas it is slightly lower than cattle Mb and higher than water buffalo at a lower pH (5.8), corresponding to the conditions in meat.

CONCLUSION

For the first time, the present study reports the purification of myoglobin from muskoxen and, furthermore, a comparative study is conducted on autoxidation and pseudoperoxidase activity with respect to cattle and water buffalo Mbs at both physiological and acid pH. Overall, the results of the current research provide novel information for future studies useful to the meat industry when considering the importance of myoglobin as a principal pigment in meat colour stability. © 2019 Society of Chemical Industry.

Genetic profiling of Hsp70 gene in Murrah buffalo (Bubalus bubalis) under sub-tropical climate of India

This study was aimed to genetic profiling of heat shock protein 70 (Hsp70) gene in Murrah buffalo investigating 50 unrelated adult animals at ICAR-Research Complex for Eastern Region, Patna (India) in winter, spring, and summer. PCR ready genomic DNA samples and season-wise total RNA samples were prepared. The PCR products of Hsp70 eluted from agarose gel were sequenced and analyzed. The first-strand cDNA was synthesized and concentration was equalized to 25 ng/μl. Expression kinetics of mRNA transcripts in different seasons was studied using Brilliant SYBR Green QPCR technique and the data retrieved was analyzed by least-squares ANOVA. DNA sequencing by primer walking revealed four allelic variants of Hsp70 gene. Alignment study revealed one substitution in 5'UTR, six substitutions in coding region, and one addition in 3'UTR. The highest percent identity and negligible phylogenetic distance were found among the alleles and reference bovine sequences. The relative mRNA expression was significantly higher in summer when THI ≥ 84 than the spring and winter; fold change increased by 4.5 times in summer than the spring whereas found nearly half in winter. These findings can be useful for heat stress management in buffaloes and help in understanding the mechanism of thermo-regulation well.

Complete mitogenome sequencing of Andaman buffalo: an endangered germplasm of Andaman and Nicobar Islands, India

Andaman buffalo is an indigenous buffalo of Andaman and Nicobar Islands, India. Over the last decade, it has witnessed a rapid decline in population, necessitating its immediate characterization and conservation. The present study reports the complete mitogenome profile of Andaman buffalo which is 16,359 bp in length and comprised of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs and two ribosomal RNAs. In addition, one A + T rich region (D-loop) was also present. A biasness towards A and T base was observed in all the genes. All the PCGs except ND6 were present on heavy strand. Start codons for all the 13 PCGs were ATN codon and abbreviated/truncated stop codons were observed in ND1, ND2, COX3, ND3 and ND4. The phylogenetic analysis revealed that the Andaman buffalo is closely related to buffalo from India and China. The results from this study will help in sketching the conservation plan of the threatened breed.

High-Throughput RNA Sequencing Reveals NDUFC2-AS lncRNA Promotes Adipogenic Differentiation in Chinese Buffalo (Bubalus bubalis L)

The buffalo (Bubalus bubalis L.) is prevalent in China and the increasing demand for meat production has changed its role from being a beast of burden to a meat source. The low fat deposition level has become one of the main barriers for its use in meat production. It is urgent to reveal factors involved in fat deposition in buffalo. This study performed RNA sequencing to investigate both long noncoding RNAs (lncRNAs) and mRNAs of adipose tissues in young and adult buffalos. A total of 124 lncRNAs and 2008 mRNAs showed differential expression patterns between young and adult samples. Coexpression analysis and functional enrichment revealed 585 mRNA–lncRNA pairs with potential function in fat deposition. After validation by qRT-PCR, we focused on a lncRNA transcribed from the ubiquinone oxidoreductase subunit C2 (NDUFC2) antisense (AS) strand which showed high correlation with thyroid hormone responsive protein (THRSP). NDUFC2-AS lncRNA is highly expressed in adipose tissue and maturation adipocytes and mainly exists in the nucleus. Functional assays demonstrated that NDUFC2-AS lncRNA promotes adipogenic differentiation by upregulating the expression levels of THRSP and CCAAT enhancer binding protein alpha (C/EBPα) in buffalo. These results indicate that NDUFC2-AS lncRNA promotes fat deposition in buffalo.

Polymorphism association and expression analysis of alpha-lactalbumin (LALBA) gene during lactation in Nili Ravi buffalo

Alpha-lactalbumin has been reported as a highly polymorphic gene that potentially alters the gene expression and is associated with milk composition in dairy breeds. Current study was conducted in two phases. In the first phase, polymorphisms identification in alpha-lactalbumin (LALBA) gene and its association with milk composition was performed. To identify the genetic polymorphism, Nili Ravi buffaloes at their second lactation were selected from Government livestock farm (Buffalo Research Institute, Pattoki). Genomic DNA was extracted from blood samples. After PCR amplification, products were sequenced, and data was analyzed. Results showed that the identified polymorphisms at chromosomal position 34310940 were found associated with major whey protein. In the second phase of study, milk samples were collected from five healthy mastitis-free Nili Ravi buffaloes in their second lactation for expression analysis of alpha-lactalbumin gene at their transition (day 15), mid (day 90), and late (day 250) lactation. Gene expression was observed highest in transition phase with a gradual decrease of expression in mid and late phase of lactation. Further studies are needed to explore the regulation of milk production genes and their translational efficiency during the course of lactation.

Heptaplex Polymerase Chain Reaction Assay for the Simultaneous Detection of Beef, Buffalo, Chicken, Cat, Dog, Pork, and Fish in Raw and Heat-Treated Food Products

Species authentication of meat and fish products is crucial to safeguard public health, economic investment, and religious sanctity. We developed a heptaplex polymerase chain reaction assay targeting short amplicon length (73-198 bp) for the simultaneous detection and differentiation of cow, buffalo, chicken, cat, dog, pig, and fish species in raw and processed food using species-specific primers targeting mitochondrial cytb, ND5, and 16s rRNA genes. Assay validation of adulterated and various heat-treated meatball matrices showed excellent stability and sensitivity under all processing conditions. The detection limit was 0.01-0.001 ng of DNA under pure states and 0.5% meat in meatball products. Buffalo was detected in 86.7% (13 out of 15) of tested commercial beef products, while chicken, pork, and fish products were found to be pure. The developed assay was efficient enough to detect target species simultaneously, even in highly degraded and processed food products at reduced time.

Reducing the cytoplasmic volume during hand-made cloning adversely affects the developmental competence and quality, and alters relative abundance of mRNA transcripts and epigenetic status of buffalo (Bubalus bubalis) embryos

Hand-made cloning (HMC) is a method of choice for somatic cell nuclear transfer (SCNT). There is 20% to 50% of cytoplasm lost during manual enucleation of oocytes with HMC. To compensate, two enucleated demicytoplasts, instead of one, are fused with each donor cell, which leads to cytoplasm pooling from two different demicytoplasts. In this study, effects of using one, instead of two demicytoplasts (controls) was examined, for production of embryos using HMC. Use of one demicytoplast decreased blastocyst development (12.7 ± 1.98% compared with 47.6 ± 3.49%, P < 0.001), total cell number (TCN, 167.6 ± 14.66 compared with 335.9 ± 58.96, P < 0.01), apoptotic index (2.11 ± 0.38 compared with 3.43±0.38, P < 0.05) but did not significantly alter inner cell mass:trophectoderm cell number ratio (0.17 ± 0.01 compared with 0.19 ± 0.02) and the global content of H3K9ac and H3K27me3 of blastocysts, compared to controls. There were gene expression alterations in pluripotency- (SOX2 and NANOG but not OCT4), epigenetic- (DNMT1 but not DNMT3a and HDAC1), apoptosis- (CASPASE3 but not BCL-2 and BAX), trophectoderm- (CDX2), development- (G6PD but not GLUT1) and cell cycle check point control-related related genes (P53) compared with controls. Transfer of cloned blastocysts from one demicytoplast (n = 8) to recipients resulted in a live calf birth that after 12 days died whereas, with transfer of control blastocysts (n = 14) there was birth of a healthy calf. In conclusion, use of one, instead of two demicytoplasts for HMC, compromises in vitro developmental competence, and alters expression of several important genes affecting embryo development.

Identification of novel polymorphism in buffalo stanniocalcin-1 gene and its expression analysis in mammary gland under different stages of lactation

Stanniocalcin-1 (STC1) is secreted by the variety of tissues having a major role in the regulation of calcium ions in the involuting mammary gland. The present work aims to sequence and structural characterization as well as expression profiling of STC1 gene in buffalo. Polymorphism identified in the 3-untranslated region (UTR) was analysed by polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) genotyping in riverine and swamp buffaloes. Expression profiling of STC1 was performed in different lactation stages of mammary gland and peripheral blood mononuclear cells to study the impact of 3'-UTR polymorphism on its expression. Different polymorphic sites were detected in the entire coding and noncoding regions of riverine and swamp buffaloes, including two INDELs. An identified polymorphic nucleotide locus A324G, having target sites for two miRNAs, namely bta-miR-2382 and bta-miR-1343, reported in cattle, was genotyped by PCR-RFLP to reveal variable allelic distribution among swamp and riverine buffaloes. Gene expression profiling across buffalo mammary tissues representing different lactation stages showed maximum expression of the STC1 gene in the involuting mammary gland. Ruminants' specific genetic variation has been observed in STC1 and its implication in buffalo mammary gland involution as well as coregulation of gene expression through miRNA binding in the 3'-UTR is suggested.

Study of whole genome linkage disequilibrium patterns of Iranian water buffalo breeds using the Axiom Buffalo Genotyping 90K Array

Accuracy of genome-wide association studies, and the successful implementation of genomic selection depends on the level of linkage disequilibrium (LD) across the genome and also the persistence of LD phase between populations. In the present study LD between adjacent SNPs and LD decay between SNPs was calculated in three Iranian water buffalo populations. Persistence of LD phase was evaluated across these populations and effective population size (Ne) was estimated from corrected r² information. A set of 404 individuals from three Iranian buffalo populations were genotyped with the Axiom Buffalo Genotyping 90K Array. Average r² and |D'| between adjacent SNP pairs across all chromosomes was 0.27 and 0.66 for AZI, 0.29 and 0.68 for KHU, and 0.32 and 0.72 for MAZ. The LD between the SNPs decreased with increasing physical distance from 100Kb to 1Mb between markers, from 0.234 to 0.018 for AZI, 0.254 to 0.034 for KHU, and 0.297 to 0.119 for MAZ, respectively. These results indicate that a density of 90K SNP is sufficient for genomic analyses relying on long range LD (e.g. GWAS and genomic selection). The persistence of LD phase decreased with increasing marker distances across all the populations, but remained above 0.8 for AZI and KHU for marker distances up to 100Kb. For multi-breed genomic evaluation, the 90K SNP panel is suitable for AZI and KHU buffalo breeds. Estimated effective population sizes for AZI, KHU and MAZ were 477, 212 and 32, respectively, for recent generations. The estimated effective population sizes indicate that the MAZ is at risk and requires careful management.