Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis )

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log₂ ratios between copy number estimates and the log₂ ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Comparative molecular characterization of Forkhead box protein 3 (FoxP3) gene of swamp-type (Bubalus carabanensis) and riverine-type (Bubalus bubalis) water buffaloes

FoxP3 is a forkhead family member that plays an important role in the development and function of a type of CD4 + T cell called T regulatory cells. Molecular characterization of FoxP3 gene in swamp- and riverine-type water buffaloes was conducted to determine its homology and compare it to the FoxP3 gene of other animal species (cattle, goat, sheep, horse, pig, cat, and dog), determine its unique characteristics in water buffaloes, and provide a reference for future studies to analyze its immunological function. FoxP3 nucleotide sequence of swamp- and riverine-type water buffaloes was 99% identical, whereas its protein translation revealed 97% homology. FoxP3 of swamp- and riverine-type water buffaloes were compared to FoxP3 of other animal species and revealed a high degree of homology which suggests that they may have the same biological properties. This study is the first report that describes the genetic characteristic of FoxP3 gene in water buffalo.

Genotype of MTNR1A gene regulates the conception rate following melatonin treatment in water buffalo

Buffaloes have tendency to show seasonal reproduction and remain in anestrus due to limited ovarian activity during summer. The seasonal reproductive behavior is ascribed the effect of melatonin related to photoperiod. Treating animals with melatonin could be a possible strategy to overcome the problem. The role of MTNR1A gene has not been fully explained in the buffalo. Therefore, we conducted a study on 114 buffalo heifers to detect the polymorphic site in MTNR1A gene and further treated them with melatonin implants to investigate the role of most frequent genotype following melatonin treatment on pregnancy. The present investigation is the first to investigate the association between melatonin treated different MTNR1A genotype buffalo and pregnancy. We confirmed SNP at position 72 in 812 bp fragment exon II of MTNR1A gene. RFLP of PCR products with Hpa I enzyme resulted in three genotypes: TT (812bp), CT (812, 743, 69bp) and CC (743, 69bp). Next, buffaloes of each genotype (TT, CC, CT; n = 28 for each) were treated with melatonin implants to compare the conception rate with their corresponding untreated control (n = 10 for each genotype). Melatonin concentrations were higher (P < 0.05) for the treatment groups of all genotypes compared to their respective untreated control from day 1-28. The pregnancy rate was significantly associated with the MTNR1A genotype. The conception rate was higher (P < 0.05) for TT genotype than for the other genotypes of buffaloes treated with melatonin. Furthermore, buffaloes of TT genotype treated with melatonin started exhibiting estrus activity soon from second week of melatonin treatment (14.1 ± 2.1; range: 10-17 days) and were found to be 7.8 times more likely to become pregnant compared to other genotypes following melatonin treatment. In conclusion, TT genotype of MTNR1A gene is more sensitive to melatonin treatment that favours pregnancy in buffaloes during summer.

Economic values and selection indices for production and reproduction traits of Iranian buffaloes (Bubalus bubalis)

Iranian buffalo plays a critical role in supplying a portion of the income and the necessities of the rural population. The first step to design a breeding program is difinition of breeding goal (BG), a linear combination of breeding values for various traits and their economic values (EV). The current study was aimed at determining EVs for important traits of Iranian buffaloes, namely milk yield (MY), milk fat (MF), age at the first calving (AFC), and calving interval (CI), as well as at estimating the genetic response of applying various types of selection indices. Economic and management data of 50 buffalo herds from various main regions of buffalo rearing in Iran were collected. The EVs were estimated using a simple profit function. Five selection indices were constructed by combining information on various traits. The EVs for BG traits of MY, MF, AFC, and CI were 0.18, 4.66, - 0.36, and - 1.87 US$, respectively. The highest predicted genetic gain in BG was 16.95 and came from applying the selection index that included all traits. The smallest genetic gain (4.93) was predicted for the index with only AFC included. Predicted genetic gain from an index that included production traits and AFC as a reproduction trait (16.9) was higher than that from the index with only production traits (16.15). Results showed that inclusion of reproductive traits in the selection index had a positive effect on genetic gain for breeding goal.

Metabolomics Reveals that Crossbred Dairy Buffaloes Are More Thermotolerant than Holstein Cows under Chronic Heat Stress

Heat stress (HS) threatens the worldwide dairy industry by decreasing animal production performance and health. Holstein cows and dairy buffaloes are the most important dairy animals, but their differences in the metabolic mechanism of thermotolerance remain elusive. In this study, we used serum metabolomics to evaluate the differences in thermotolerance between Holstein cows and crossbred dairy buffaloes under chronic heat stress (HS) and thermal-neutral conditions. In response to HS, the body temperatures and respiratory rates were increased more for Holstein cows than for dairy buffaloes (38.78 vs 38.24 °C, p < 0.001; 43.6 vs 32.5 breaths/min, p < 0.001). HS greatly affected serum metabolites associated with amino acids, fatty acids, and bile acids. The enriched metabolic pathways of these serum metabolites are closely related to HS. We demonstrated that buffaloes adapt to HS by adopting a metabolism of branched-chain amino acids and ketogenic amino acids and gluconeogenesis, but Holstein cows decrease the effect of HS with citrulline and proline metabolism. Both physiological parameters and serum metabolic profiles indicate that dairy buffaloes are more thermotolerant than Holstein cows, providing the feasibility to vigorously develop the buffalo dairy industry in tropical and subtropical regions.

Direct pentaplex PCR assay: An adjunct panel for meat species identification in Asian food products

A direct pentaplex PCR assay was developed for the identification of meat from sources other than those declared on the packaging. Species-specific primers were designed, based on the mitochondrial cytochrome oxidase I (COI) gene. The assay amplified specific DNA fragments from dog (230 bp), duck (283 bp), buffalo (363 bp), goat (396 bp), and sheep (477 bp). The proposed method is capable of identifying target species accurately and is reproducible, sensitive and robust for use with real-world foods and food products. In total, 26 of 117 meat and commercial food products tested were shown to contain DNA from species not declared on the label.

First molecular epidemiological study of Mycobacterium avium subsp. paratuberculosis in cattle and buffalo from different regions of Brazil

Paratuberculosis is an incurable disease in ruminants with great worldwide economic impact, caused by Mycobacterium avium subsp. paratuberculosis (MAP). The objective of this study was to carry out a study of the molecular epidemiology of the MAP using the restriction enzyme analysis (REA) technique of IS1311 MAP region in biological samples of feces, intestinal tissue, and mesenteric lymph nodes of cattle and buffaloes from six Brazilian states. In total, 109 samples of feces and tissues of cattle and buffaloes were collected from animal paratuberculosis suspected. Twenty-five samples were positive in the detection of the DNA of the IS900 region of MAP and it was possible to type 18 strains in the analysis of the region IS1311, being 100% of them identified as belonging to subtype Bison MAP strain. This is the first epidemiological molecular study of MAP in Brazil. The results indicate that paratuberculosis is widespread in cattle and in buffaloes in several regions of Brazil, and the subtype Bison MAP strain was the only one identified in the samples analyzed in this study, demonstrating the similarity between the strains from different states tested. These results provide the necessary support for the implementation of paratuberculosis control strategies in cattle and buffaloes in Brazil.

Comparative Analysis of V-Akt Murine Thymoma Viral Oncogene Homolog 3 (AKT3) Gene between Cow and Buffalo Reveals Substantial Differences for Mastitis

AKT3 gene is a constituent of the serine/threonine protein kinase family and plays a crucial role in synthesis of milk fats and cholesterol by regulating activity of the sterol regulatory element binding protein (SREBP). AKT3 is highly conserved in mammals and its expression levels during the lactation periods of cattle are markedly increased. AKT3 is highly expressed in the intestine followed by mammary gland and it is also expressed in immune cells. It is involved in the TLR pathways as effectively as proinflammatory cytokines. The aims of this study were to investigate the sequences differences between buffalo and cow. Our results showed that there were substantial differences between buffalo and cow in some exons and noteworthy differences of the gene size in different regions. We also identified the important consensus sequence motifs, variation in 2000 upstream of ATG, substantial difference in the "3'UTR" region, and miRNA association in the buffalo sequences compared with the cow. In addition, genetic analyses, such as gene structure, phylogenetic tree, position of different motifs, and functional domains, were performed to establish their correlation with other species. This may indicate that a buffalo breed has potential resistance to disease, environment changes, and airborne microorganisms and some good production and reproductive traits.

Selection and evaluation of housekeeping genes as endogenous controls for quantification of mRNA transcripts in Theileria parva using quantitative real-time polymerase chain reaction (qPCR)

The reliability of any quantitative real-time polymerase chain reaction (qPCR) experiment can be seriously compromised by variations between samples as well as between PCR runs. This usually result from errors in sample quantification, especially with samples that are obtained from different individuals and tissues and have been collected at various time intervals. Errors also arise from differences in qPCR efficiency between assays performed simultaneously to target multiple genes on the same plate. Consequently, the derived quantitative data for the target genes become distorted. To avoid this grievous error, an endogenous control, with relatively constant transcription levels in the target individual or tissue, is included in the qPCR assay to normalize target gene expression levels in the analysis. Several housekeeping genes (HKGs) have been used as endogenous controls in quantification studies of mRNA transcripts; however, there is no record in the literature of the evaluation of these genes for the tick-borne protozoan parasite, Theileria parva. Importantly, the expression of these genes should be invariable between different T. parva stocks, ideally under different experimental conditions, to gain extensive application in gene expression studies of this parasite. Thus, the expression of several widely used HKGs was evaluated in this study, including the genes encoding β-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 28S rRNA, cytochrome b and fructose-2.6-biphosphate aldolase (F6P) proteins. The qPCR analysis revealed that the expression of genes encoding cytochrome b, F6P and GAPDH varied considerably between the two T. parva stocks investigated, the cattle-derived T. parva Muguga and the buffalo-derived T. parva 7014. 28S rRNA and β-actin gene expression was the most stable; thus, these genes were considered suitable candidates to be used as endogenous control genes for mRNA quantification studies in T. parva.

Transcriptional and translational abundance of Bone morphogenetic protein (BMP) 2, 4, 6, 7 and their receptors BMPR1A, 1B and BMPR2 in buffalo ovarian follicle and the role of BMP4 and BMP7 on estrogen production and survival of cultured granulosa cells

BMPs and their receptors modulate the granulosa cell (GC) function in the follicle of domestic animals. Since little is known on BMPs in the buffalo, the present study was aimed to investigate the expression of BMP2, 4, 6, 7 and their receptors BMPR1A, BMPR1B, BMPR2 in the GC and theca cells (TC) of ovarian follicles and the role of BMP4 and BMP7 on buffalo GC. Follicles were classified into four groups based on size and E₂ level in the follicular fluid as follows: (i) Group1(4-6 mm; <0.5 ng/mL) (ii) Group 2 (7-9 mm; 0.5-5 ng/mL) (iii) Group 3 (10-13 mm; 5-40 ng/mL) and (iv) Group 4 (dominant follicle) (>13 mm; >180 ng/mL). The results revealed that except BMP6, BMP2, 4 7 and receptors BMPR1A, BMPR1B and BMPR2 showed a minimum of 1.5-2 fold increase in mRNA expression in the GC of dominant follicle as compared to other follicle classes. In the dominant follicle, a two-fold increase in BMP4 and BMP7 expression was observed in the TC. At 100 ng/mL, the BMP4 and BMP7 either alone or in combination maximally down-regulated CASPASE3 and stimulated the transcripts of PCNA, FSHR and CYP19A1 that was supported by E₂ secretion in the granulosa cell culture suggesting their role in cell survival and E2 production. In conclusion, GC and TC of dominant follicles express BMP 2, 4, 6, 7 and their receptors BMPR1A, BMPR1B and BMPR2. BMP4 and BMP7 stimulate E₂ production and promote GC survival.

Efficient genome editing in cultured cells and embryos of Debao pig and swamp buffalo using the CRISPR/Cas9 system

Myostatin (MSTN), a protein encoded by growth differentiation factor 8 (GDF8), is primarily expressed in skeletal muscle and negatively regulates the development and regeneration of muscle. Accordingly, myostatin-deficient animals exhibit a double-muscling phenotype. The CRISPR/Cas9 system has proven to be an efficient genome-editing tool and has been applied to gene modification in cells from many model organisms such as Drosophila melanogaster, zebrafish, mouse, rat, sheep, and human. Here, we edited the GDF8 gene in fibroblasts and embryos of Debao pig and swamp buffalo using the CRISPR/Cas9 system. The CRISPR/Cas9-mediated mutation efficiency in fibroblasts was as high as 87.5% in pig and 78.9% in buffalo. We then obtained single-cell clones with mutations at the specific sites of the GDF8 gene by screening with G418 in fibroblasts of pig and buffalo. In addition, the frequencies of Cas9/gRNA-mediated mutations were at 36 and 25% in the intracytoplasmic sperm injection embryos of pig and in vitro fertilization embryos of buffalo, respectively. Our work demonstrates that the Cas9/gRNA system is a highly efficient and fast tool for genome editing in cultured cells and embryos of Debao pig and swamp buffalo. These results can be helpful for the establishment of a new animal strain that can generate more meat.

Selection of suitable reference genes for normalization of quantitative RT-PCR (RT-qPCR) expression data across twelve tissues of riverine buffaloes (Bubalus bubalis)

Selection of reference genes has become an integral step in any real time quantitative PCR (RT-qPCR) based expression studies. The importance of this study stems from the fact that riverine buffaloes are major dairy species of Indian sub-continent and the information generated here will be of great interest to the investigators engaged in functional genomic studies of this important livestock species. In this study, an effort was made to evaluate a panel of 10 candidate reference genes (glyceraldehyde 3-phosphate dehydrogenase (GAPDH), beta- actin (ACTB), ubiquitously expressed transcript (UXT), ribosomal protein S15 (RPS15), ribosomal protein L-4 (RPL4), ribosomal protein S9 (RPS9), ribosomal protein S23 (RPS23), hydroxymethylbilane synthase (HMBS), β2 Microglobulin (β2M) and eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) across 12 tissues (mammary gland, kidney, spleen, liver, heart, intestine, ovary, lung, muscle, brain, subcutaneous fat and testis) of riverine buffaloes. In addition to overall analysis, tissue wise evaluation of expression stability of individual RG was also performed. Three different algorithms provided in geNorm, NormFinder and BestKeeper softwares were used to evaluate the stability of 10 potential reference genes from different functional classes. The M-value given by geNorm ranged from 0.9797 (RPS9 and UXT) to 1.7362 (RPS15). From the most stable to the least stable, genes were ranked as: UXT/RPS9> RPL4> RPS23> EEF1A1> ACTB> HMBS> GAPDH> B2M> RPS15. While NormFinder analysis ranked the genes as: UXT> RPS23> RPL4> RPS9> EEF1A1> HMBS> ACTB> β2M> GAPDH> RPS15. Based on the crossing point SD value and range of fold change expression, BestKeeper analysis ranked the genes as: RPS9> RPS23/UXT> RPL4> GAPDH> EEF1A1> ACTB> HMBS> β2M> RPS15. Overall the study has identified RPS23, RPS9, RPL4 and UXT genes to be the most stable and appropriate RGs that could be utilized for normalization of transcriptional data in various tissues of buffaloes. This manuscript thus provide useful information on panel of reference genes that could be helpful for researchers conducting functional genomic studies in riverine buffaloes.

Testosterone levels and the genetic variation of sex hormone-binding globulin gene of Bubalus bubalis, bulls in Egypt

Testosterone is themain plasma androgen produced and secreted by male testis. To be bioactive, it binds to sex hormonebinding globulin (SHBG), i.e. the major transporter protein of sex steroids in the blood of mammals and other vertebrate species.Firstly, this study aimed to determine the levels of testosterone in different-age groups of buffalo bulls (Bubalus bubalis) in Egypt, and secondly to screen the genetic polymorphisms in their SHBG gene and to investigatewhether these polymorphisms are associated with the level of the circulating serum testosterone. Blood samples were collected from 67 Egyptian buffalo bulls representing four different age-groups and testosterone concentration was measured by radioimmunoassay. Selected DNA samples from different age-groups were sequenced for the full length of theirSHBGgene.Asignificant increase in the concentration of serum testosterone in buffalo bulls was recorded with the advancement of age from 12/18 to 24/36 months. Based on the number of polymorphisms in SHBG gene, the current study reveals the presence of three bull genotypes, of which one is likely to be associated with low testosterone concentration, while other is likely to be associated with an increased testosterone concentration. Such genetic associations can provide a good tool for the selection of bull genotypes with higher testosterone concentrations, which are indispensible for breeding purposes.

Genetic responsiveness of African buffalo to environmental stressors: A role for epigenetics in balancing autosomal and sex chromosome interactions?

In the African buffalo (Syncerus caffer) population of the Kruger National Park (South Africa) a primary sex-ratio distorter and a primary sex-ratio suppressor have been shown to occur on the Y chromosome. A subsequent autosomal microsatellite study indicated that two types of deleterious alleles with a negative effect on male body condition, but a positive effect on relative fitness when averaged across sexes and generations, occur genome-wide and at high frequencies in the same population. One type negatively affects body condition of both sexes, while the other acts antagonistically: it negatively affects male but positively affects female body condition. Here we show that high frequencies of male-deleterious alleles are attributable to Y-chromosomal distorter-suppressor pair activity and that these alleles are suppressed in individuals born after three dry pre-birth years, likely through epigenetic modification. Epigenetic suppression was indicated by statistical interactions between pre-birth rainfall, a proxy for parental body condition, and the phenotypic effect of homozygosity/heterozygosity status of microsatellites linked to male-deleterious alleles, while a role for the Y-chromosomal distorter-suppressor pair was indicated by between-sex genetic differences among pre-dispersal calves. We argue that suppression of male-deleterious alleles results in negative frequency-dependent selection of the Y distorter and suppressor; a prerequisite for a stable polymorphism of the Y distorter-suppressor pair. The Y distorter seems to be responsible for positive selection of male-deleterious alleles during resource-rich periods and the Y suppressor for positive selection of these alleles during resource-poor periods. Male-deleterious alleles were also associated with susceptibility to bovine tuberculosis, indicating that Kruger buffalo are sensitive to stressors such as diseases and droughts. We anticipate that future genetic studies on African buffalo will provide important new insights into gene fitness and epigenetic modification in the context of sex-ratio distortion and infectious disease dynamics.

Design and validation of a 90K SNP genotyping assay for the water buffalo (Bubalus bubalis)

Background

The availability of the bovine genome sequence and SNP panels has improved various genomic analyses, from exploring genetic diversity to aiding genetic selection. However, few of the SNP on the bovine chips are polymorphic in buffalo, therefore a panel of single nucleotide DNA markers exclusive for buffalo was necessary for molecular genetic analyses and to develop genomic selection approaches for water buffalo. The creation of a 90K SNP panel for river buffalo and testing in a genome wide association study for milk production is described here.

Methods

The genomes of 73 buffaloes of 4 different breeds were sequenced and aligned against the bovine genome, which facilitated the identification of 22 million of sequence variants among the buffalo genomes. Based on frequencies of variants within and among buffalo breeds, and their distribution across the genome, inferred from the bovine genome sequence, 90,000 putative single nucleotide polymorphisms were selected to create an Axiom® Buffalo Genotyping Array 90K.

Results

This 90K “SNP-Chip” was tested in several river buffalo populations and found to have ∼70% high quality and polymorphic SNPs. Of the 90K SNPs about 24K were also found to be polymorphic in swamp buffalo. The SNP chip was used to investigate the structure of buffalo populations, and could distinguish buffalo from different farms. A Genome Wide Association Study identified genomic regions on 5 chromosomes putatively involved in milk production.

Conclusion

The 90K buffalo SNP chip described here is suitable for the analysis of the genomes of river buffalo breeds, and could be used for genetic diversity studies and potentially as a starting point for genome-assisted selection programmes. This SNP Chip could also be used to analyse swamp buffalo, but many loci are not informative and creation of a revised SNP set specific for swamp buffalo would be advised.

PCR-SSCP analysis of MDGI gene and its association with milk production traits in river buffalo (Bubalus bubalis)

In this study, we investigated the genetic variation within 3'UTR of Mammary-Derived Growth Inhibitor (MDGI) gene of buffalo using PCR-SSCP and sequencing; and also analyzed association of polymorphism with the milk production traits. The study revealed two conformational patterns, 'A' and 'B' among 234 Mehsana buffaloes maintained with their records in the field and at farm. The frequency of SSCP variant 'A' was found to be invariably high in the buffalo population under study. Further, association analysis of SSCP variants with various milk production and milk quality traits indicated no significant effect on any of the traits investigated. Sequencing of SSCP variant 'A' showed homozygous G/G and A/A and 'B' had heterozygous G/C and A/G at positions +124 and +140 respectively, in the 3'UTR of buffalo MDGI. The preliminary results showed the substantial variations in the distribution of SSCP variants' frequencies within Mehsana buffaloes, however these variants had non-significant association with milk yield, fat yield and fat percentage in Mehsana buffaloes.

Quantitative Tetraplex Real-Time Polymerase Chain Reaction Assay with TaqMan Probes Discriminates Cattle, Buffalo, and Porcine Materials in Food Chain

Cattle, buffalo, and porcine materials are widely adulterated, and their quantification might safeguard health, religious, economic, and social sanctity. Recently, conventional polymerase chain reaction (PCR) and PCR-restriction fragment length polymorphism (RFLP) assays have been documented but they are just suitable for identification, cannot quantify adulterations. We described here a quantitative tetraplex real-time PCR assay with TaqMan Probes to quantify contributions from cattle, buffalo, and porcine materials simultaneously. Amplicon-sizes were very short (106-, 90-, and 146-bp for cattle, buffalo, and porcine) because longer targets could be broken down, bringing serious ambiguity in molecular diagnostics. False negative detection was eliminated through an endogenous control (141-bp site of eukaryotic 18S rRNA). Analysis of 27 frankfurters and 27 meatballs reflected 84-115% target recovery at 0.1-10% adulterations. Finally, a test of 36 commercial products revealed 71% beef frankfurters, 100% meatballs, and 85% burgers contained buffalo adulteration, but no porcine was found in beef products.

Meat quality of buffaloes finished in traditional or silvopastoral system in the Brazilian Eastern Amazon

BACKGROUND

The present study aimed to assess the physical, chemical and sensory characteristics of meat from buffaloes finished in a Traditional System (TS) or Traditional + Silvopastoral System (TSPS) with dietary supplementation. Crossbred Murrah × Mediterranean buffaloes were raised from weaning to slaughter in the TS (n = 15) or were raised in the traditional system and finished in the TSPS (n = 15). After finishing, animals were slaughtered and their carcasses refrigerated for 24 h. The right side of each half-carcass was cut between the 12th and 13th ribs and the Longissimus thoracis muscle was removed. The cranial part underwent analyses of pH, color, weight loss as a result of cooking, water holding capacity, texture and sensory characteristics, whereas the rest of the muscle underwent microbiological analyses and determination of the chemical composition, fatty acid profile and mineral content.

RESULTS

No major difference between finishing systems was found (P > 0.05) in the physical analyses and chemical composition of meat. The percentage of myristic acid (C14:0) and the sum of polyunsaturated fatty acids differed between treatments. The TS meat had the best 'characteristic meat aroma'.

CONCLUSION

Considering the quality of meat produced in the TS or TSPS, it is concluded that finishing buffaloes in the pasture still represents the best alternative. © 2016 Society of Chemical Industry.

Impact of levels of total digestible nutrients on microbiome, enzyme profile and degradation of feeds in buffalo rumen

The present study was aimed at understanding a shift in rumen microbiome of buffaloes fed various levels of total digestible nutrients. To understand the process, the metagenomics of rumen microbes, in vivo and in vitro rumen fermentation studies were carried out. Three rumen fistulated adult male Murrah buffaloes were fed three isonitrogenous diets varying in total digestible nutrients (70, 85 and 100% of TDN requirement) in 3X3 switch over design. On dry matter basis, wheat straw/ roughage content were 81, 63 and 51% and that of maize grain was 8, 16 and 21% in three diets respectively. After 20 d of feeding, rumen liquor and rumen contents were sampled just before (0h) and 4h post feeding. Ruminococcus flavefaciens and R. albus (estimated with real time PCR) were higher in high roughage diets. The predominant phyla in all the three groups were Bacteroidetes, Firmicutes followed by Proteobacteria, Actinobacteria and Fibrobacteres. A core group of more than fifty rumen bacteria was present in all the animals with very little variations due to level of TDN. The most predominant bacterial genera reported in order of decreasing abundance were: Prevotella, Bacteroides, Clostridium, Ruminococcus, Eubacterium, Parabacteroides, Fibrobacter, Butyrivibrio etc. The higher diversity of the enyzmes families GH 23, GH 28, GH 39, GH 97, GH 106, and GH 127 (the enzymes active in fibre and starch degradation) were significantly higher on 100%TDN diet while CE 14 (required for the hydrolysis of bond between carbohydrate and lignin) was higher on low TDN (70%) diet, indicating ester bond cleavage was better in animals fed high roughage (wheat straw) diet.

Expression and functional role of Bone Morphogenetic Proteins (BMPs) in cyclical corpus luteum in buffalo (Bubalus bubalis)

The role of growth factors in the modulation of ovarian function is an interesting area of research in reproductive biology. Recently, we have shown the expression and role of IGF, EGF, VEGF and FGF in the follicle and CL. Here, we report the presence of Bone Morphogenetic Proteins (BMPs) and their functional receptors in the corpus luteum (CL) of buffalo. The bubaline CL was classified into four stages according to the morphology and progesterone (P₄) concentration. The qPCR, immunoblot and immunohistochemistry studies revealed that BMP2 and BMP Receptors (BMPR1A, BMPR1B and BMPR2) were significantly upregulated during the mid stage whereas BMP4 and BMP7 were upregulated during the early stage of CL (P<0.05). Studies on primary luteal cell culture (LCC) using mid CL showed a significant time and concentration dependent effect of BMP4 and BMP7 (P<0.05). At 100ngml^-1, the BMPs maximally stimulated the transcripts of StAR, CYP11A1 and 3βHSD that paralleled with P₄ accretion in the media (P<0.05). Further, the BMP4 as well as BMP7 upregulated the transcripts of PCNA and downregulated CASPASE3 in the LCC at the same concentration (P<0.05). Though the combined effect of BMP4 and 7 was significantly higher (P<0.05) than that of individual one, it was not additive. In conclusion, the expression of BMPs and their receptors were dependent on the stages of CL in the buffalo. Treatment of LCC with BMPs in vitro confirmed the presence of functional receptors that stimulated the P₄ production and luteal cell survival. Moreover, the results support the concept that the upregulation of P₄ and its biosynthetic pathway enzymes such as CYP11A1, StAR and 3βHSD in the CL is likely due to the autocrine and /or paracrine effects of BMP4 and BMP7 under physiological milieu.

Structural, Functional and Phylogenetic Analysis of Sperm Lysozyme-Like Proteins

Sperm lysozyme-like proteins belonging to c-type lysozyme family evolved in multiple forms. Lysozyme-like proteins, viz., LYZL2, LYZL3 or SLLP1, LYZL4, LYZL5 and LYZL6 are expressed in the testis of mammals. Not all members of LYZL family have been uniformly and unambiguously identified in the genome and proteome of mammals. Some studies suggested a role of SLLP1 and LYZL4 in fertilization; however, the function of other LYZL proteins is unknown. We identified all known forms of LYZL proteins in buffalo sperm by LC-MS/MS. Cloning and sequence analysis of the Lyzl cDNA showed 38-50% identity at amino acid level among the buffalo LYZL paralogs, complete conservation of eight cysteines and other signature sequences of c-type lysozyme family. Catalytic residues in SLLP1, LYZL4 and LYZL5 have undergone replacement. The substrate binding residues showed significant variation in LYZL proteins. Residues at sites 62, 101, 114 in LYZL4; 101 in SLLP1; 37, 62, and 101 in LYZL6 were more variable among diverse species. Sites 63 and 108 occupied by tryptophan were least tolerant to variation. Site 37 also showed lower tolerance to substitution in SLLP1, LYZL4 and LYZL5, but more variable in non-testicular lysozymes. Models of LYZL proteins were created by homology modeling and the substrate binding pockets were analyzed in term of binding energies and contacting residues of LYZL proteins with tri-N-acetylglucosamine (NAG)3 in the A-B-C and B-C-D binding mode. Except LYZL6, LYZL proteins did not show significant difference in binding energies in comparison to hen egg white lysozyme in the A-B-C mode. (NAG)3 binding energy in the B-C-D mode was higher by 1.3-2.2 kcal/mol than in A-B-C mode. Structural analysis indicated that (NAG)3 was involved in making more extensive interactions including hydrogen bonding with LYZL proteins in B-C-D mode than in A-B-C mode. Despite large sequence divergence among themselves and with respect to c-type lysozymes, substrate binding residues as well as hydrogen bonding network between (NAG)3 and proteins were mostly conserved. LYZL5 in buffalo and other mammalian species contained additional 10-12 amino acid sequence at c-terminal that matched with ankyrin repeat domain-containing protein 27. Phylogenetic analysis indicated LYZL2 to be most ancient among all the LYZL proteins and that the evolution of LYZL proteins occurred through several gene duplications preceding the speciation of mammals from other vertebrates as distant as reptiles and amphibians.

Double Gene Targeting Multiplex Polymerase Chain Reaction-Restriction Fragment Length Polymorphism Assay Discriminates Beef, Buffalo, and Pork Substitution in Frankfurter Products

Beef, buffalo, and pork adulteration in the food chain is an emerging and sensitive issue. Current molecular techniques to authenticate these species depend on polymerase chain reaction (PCR) assays involving long and single targets which break down under natural decomposition and/or processing treatments. This novel multiplex polymerase chain reaction-restriction fragment length polymorphism assay targeted two different gene sites for each of the bovine, buffalo, and porcine materials. This authentication ensured better security, first through a complementation approach because it is highly unlikely that both sites will be missing under compromised states, and second through molecular fingerprints. Mitochondrial cytochrome b and ND5 genes were targeted, and all targets (73, 90, 106, 120, 138, and 146 bp) were stable under extreme boiling and autoclaving treatments. Target specificity and authenticity were ensured through cross-amplification reaction and restriction digestion of PCR products with AluI, EciI, FatI, and CviKI-1 enzymes. A survey of Malaysian frankfurter products revealed rampant substitution of beef with buffalo but purity in porcine materials.

Corrigendum Polymorphism in the melatonin receptor gene in buffalo populations of the Brazilian Amazon - Genet. Mol. Res. 15 (2): gmr.15027960

The correction is only in the name of the first author and should be: E.M. Barbosa(1), B.B. Souza(2), R.C. Guimarães(2), J.S.N. Azevedo(3), E.C. Gonçalves(4), H.F.L. Ribeiro(2), S.T. Rolim Filho(2), E. Silva Filho(2).

De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.