Development and validation of most efficient RNA isolation method from buffalo bull spermatozoa

BACKGROUND

Being highly fragmented and low in concentration, isolation of good quality RNA from sperm cells is a big challenge. Attempts have been made to evaluate various sperm RNA isolation methods from purified buffalo bull sperm cells.

METHODS

Both, non-membrane and membrane-based methods have been evaluated for isolating RNA from Murrah buffalo sperms and compared for their respective efficacies. The traditional TRIzol, TRIzol-heat lysed (H-TRIzol) and cocktail of TCEP-RLT lysis buffer (Qiagen RNeasy mini kit)-TRIzol (C-TRIzol) based isopropanol isolation methods have been evaluated.

RESULTS

H-TRIzol yielded best results among conventional methods. The combined T-RLT RNA isolation protocol yielded best quality and quantity compared to other membrane-based methods, due to high lytic property of cocktail of lysis reagents, necessary for complete breakdown of sperm membrane and RNA binding membrane for RNA isolation. Combined lysis performed by treatment with RLT-T and T-RLT differing in order of reagents used were also evaluated. T-RLT combination giving better results compared to RLT-T due to high gDNA contamination and membrane clogging in later protocol steps.

CONCLUSION

Overall, in terms of total RNA quantity and quality per million spermatozoa, the heat-lysed TRIzol method (H-TRIzol) performs best among RNA separation techniques employed and is also quite easy to perform. This comparative evaluation of sperm RNA isolation protocols can be useful in deciding the best protocol for isolation of good quality and high concentration sperm RNA from buffalo semen, for transcriptome and other downstream studies.

Estimation of genetic parameters for semen traits in Egyptian buffalo bulls

This study was conducted to characterize semen traits (ejaculate volume (VOL), mass motility (MM), sperm livability (LS), percentage of abnormal sperms (AS), and sperm concentration (CONC)) of Egyptian buffalo bulls and evaluate the importance of some nongenetic factors (year (YC) and season (SC) of semen collection and age of bull genetically and environmentally at collection (ABC)) affecting the investigated traits. A total of 7761 normal semen ejaculates were collected from 26 bulls from 2009 to 2019. Single-trait and bivariate repeatability animal models using Bayesian methods were used to estimate variance components, heritability, repeatability, and genetic correlations among the investigated semen traits. YC and ABC exerted significant effects on most semen traits, whereas SC exerted no significant effect on all the investigated semen traits. Heritability estimates were 0.08, 0.52, 0.51, 0.04, and 0.49 for VOL, MM, LS, AS, and CONC, respectively. Repeatability estimates were 0.14, 0.82, 0.79, 0.06, and 0.78 for VOL, MM, LS, AS, and CONC, respectively. The genetic correlations between MM and each of LS and CONC were highly significant (0.99 ± 0.01 and 0.95 ± 0.14, respectively), and that between LS and CONC was also highly significant (0.92 ± 0.20). The high heritability estimates for MM, LS, and CONC combined with the favorable high significant genetic correlations between these traits indicated that direct selection for MM may be an effective method to enhance semen quality in Egyptian buffalo bulls and consequently improve fertility.

Review: Development, adoption, and impact of assisted reproduction in domestic buffaloes

The domestic buffalo (Bubalus bubalis), also known as water buffalo, comprises two sub-species the River buffalo (B. bubalis ssp. bubalis; 50 chromosomes) and the Swamp buffalo (ssp. carabanensis; 48 chromosomes). Domestic buffaloes are a globally significant livestock species. In South Asia, the River buffalo is a primary source of milk and meat and has a very important role in food security. The River buffalo also supports high-value, differentiated food production in Europe and the Americas. The Swamp buffalo is an important draft animal and a source of food in Southeast Asia and East Asia. The growing importance of buffaloes requires that they undergo an accelerated rate of genetic gain for efficiency of production, product quality, and sustainability. This will involve the increased use of assisted reproduction. The initial application of reproductive technology in buffaloes had variable success as it relied on the adoption of procedures developed for cattle. This included artificial insemination (AI), sperm cryopreservation, and embryo technologies such as cloning and in vitro embryo production (IVEP). Reproductive technology has been progressively refined in buffaloes, and today, the success of AI and IVEP is comparable to cattle. Ovarian follicular superstimulation (superovulation) combined with in vivo embryo production results in low embryo recovery in buffaloes and has limited practical application. The contribution of elite female buffaloes to future genetic improvement will therefore rely mainly on oocyte pickup and IVEP. This will include IVEP from females before puberty to reduce generation intervals. This review provides for the first time a clear chronology on the development, adoption, and impact, of assisted reproduction in domestic buffaloes.

Non-genetic effects and repeatability estimates of semen production traits in Murrah buffalo bulls

The objectives of this study were to determine the repeatability estimates and the effect of ejaculate, season, period, age, and bull on semen production traits in Murrah buffalo bulls maintained in an organized semen station. A data set comprising 54,268 ejaculates from 76 Murrah buffalo bulls collected across 12 years were analyzed. The analysis of variables included non-genetic factors as the fixed effects and bull as the random effect. The repeatability was estimated by the restricted maximum likelihood method (REML) using WOMBAT program. The overall least-squares means for semen volume, sperm concentration, total sperm/ejaculate, mass activity, individual motility, post-thaw motility, and frozen semen straws/ejaculate were 2.65 ml, 1222.04 million/ml, 3030.10 million, 2.64, 67.45%, 51.73%, and 128.80 doses. The repeatability estimates were 0.27, 0.25, 0.22, 0.28, 0.34, 0.27, and 0.23, respectively. The studied non-genetic factors influenced (P < 0.01) all seminal traits. Likewise, variation (P < 0.01) due to individual bull was also observed. The first ejaculate was superior in terms of ejaculate volume, sperm concentration, total sperm/ejaculate, and frozen semen straws than the second ejaculate. The south-west monsoon and summer seasons were the best seasons with better semen quality and frozen semen production. Compared to the remaining periods, during Period V (2014-2016), greater values for most seminal traits were found. Bulls above 8 years of age expressed higher semen production than the younger bulls below 4 years, whereas the best performance was observed in the age group of 12 years and above. The repeatability estimates were low to moderate ranging from 0.22 for total sperm/ejaculate to 0.34 for individual motility. A significant (P < 0.01) positive correlation was observed between the seminal traits except the negative correlation in semen volume and total sperm/ejaculate. It was concluded that all the non-genetic factors considered in this study exerted influence on studied seminal traits. Therefore, appropriate management practices can be applied to improve semen quality and quantity.

Novel Insight into the Role of Squalene Epoxidase (SQLE) Gene in Determining Milk Production Traits in Buffalo

Understanding the genetic mechanisms underlying milk production traits contribute to improving the production potential of dairy animals. Squalene epoxidase (SQLE) is one of the rate-limiting enzymes for cholesterol biosynthesis and was highly expressed in the buffalo mammary. The objectives of the present study were to detect the polymorphisms within SQLE in buffalo, the genetic effects of these mutations on milk production traits, and to understand the gene regulatory effects on buffalo mammary epithelial cells (BuMECs). A total of five SNPs were identified by sequencing, g.18858G > A loci were significantly associated with fat yield, and g.22834C > T loci were significantly associated with peak milk yield, milk yield, fat yield, and protein yield. Notably, linkage disequilibrium analysis indicated that 2 SNPs (g.18858G > A and g.22834C > T) formed one haplotype block, which was found to be significantly associated with milk fat yield, fat percentage, and protein yield. Furthermore, expression of SQLE was measured in different tissues of buffalo and was found to be higher in the mammary. Knockdown of SQLE gene expression significantly affected the growth of BuMECs, including proliferation, cell cycle, and apoptosis, and significantly downregulated the expression of related genes MYC, PCNA, and P21. In addition, knockdown of the SQLE gene significantly reduces triglyceride concentrations and the signal intensity of oil red O staining. In addition, silencing of SQLE was also found to regulate the synthesis and secretion of β-casein and κ-casein negatively. Furthermore, SQLE knockdown is accompanied by the downregulation of critical genes (RPS6KB1, JAK2, eIF4E, and SREBP1) related to milk fat and protein synthesis. The current study showed the potential of the SQLE gene as a candidate for buffalo milk production traits. It provides a new understanding of the physiological mechanisms underlying buffalo milk production regulation.

Chromosome-level genome and recombination map of the male buffalo

BACKGROUND

The swamp buffalo (Bubalus bubalis carabanesis) is an economically important livestock supplying milk, meat, leather, and draft power. Several female buffalo genomes have been available, but the lack of high-quality male genomes hinders studies on chromosome evolution, especially Y, as well as meiotic recombination.

RESULTS

Here, a chromosome-level genome with a contig N50 of 72.2 Mb and a fine-scale recombination map of male buffalo were reported. We found that transposable elements (TEs) and structural variants (SVs) may contribute to buffalo evolution by influencing adjacent gene expression. We further found that the pseudoautosomal region (PAR) of the Y chromosome is subject to stronger purification selection. The meiotic recombination map showed that there were 2 obvious recombination hotspots on chromosome 8, and the genes around them were mainly related to tooth development, which may have helped to enhance the adaption of buffalo to inferior feed. Among several genomic features, TE density has the strongest correlation with recombination rates. Moreover, the TE subfamily, SINE/tRNA, is likely to play a role in driving recombination into SVs.

CONCLUSIONS

The male genome and sperm sequencing will facilitate the understanding of the buffalo genomic evolution and functional research.

Identification of differential methylome signatures of white pigmented skin patches in Nili Ravi buffalo of India

The DNA methylation events mark a major epigenetic change in the genome, reflecting non-genetic disease developments and varied phenotypes. The water buffalo is a dairy production animal with wide agro-climatic distribution in India. Breed-wise the coat color of water buffalo varies from ash-gray to jet black. A typical pigmentation pattern is found in one of the breeds of North India, Nili Ravi, with variedly distributed white patches. The DNA methylation pattern could potentially reveal the epigenetic factors responsible for the pigmentation patterns. To address this question, the DNA isolated from the skin tissues of Nili Ravi with varied white pigmentation and black Murrah buffaloes was subjected to reduced representation bisulfite sequencing. DNA methylation analysis revealed, 68.44%, 63.39%, and 47.94% of the promoter regions were hypermethylated in Nili Ravi over-white versus Murrah, Nili Ravi under-white versus Murrah, and Nili Ravi under-white versus Nili Ravi over-white, respectively. Major genes identified to be differentially methylated among over-white and under-white skin tissues in Nili Ravi included TBX2, SNAI2, HERC2, and CITED1. Overall the results have indicated differential methylation patterns to be potentially involved in hyper or hypopigmentation in Nili Ravi and Murrah buffaloes.

CircCLTH promotes skeletal muscle development and regeneration

Buffalo holds an excellent potential for beef production, and circRNA plays an important role in regulating myogenesis. However, the regulatory mechanism of circRNAs during buffalo skeletal muscle development has not been fully explored. In this study, circRNA expression profiles during the proliferation and differentiation stages of buffalo myoblasts were analysed by RNA-seq. Here, a total of 3,142 circRNAs candidates were identified, and 110 of them were found to be differentially expressed in the proliferation and differentiation stages of buffalo myoblast libraries. We focused on a 347 nt circRNA subsequently named circCLTH. It consists of three exons and is expressed specifically in muscle tissues. It is a highly conserved non-coding RNA with about 95% homology to both the human and the mouse circRNAs. The results of cell experiments and RNA pull-down assays indicated that circCLTH may capture PLEC protein, promote the proliferation and differentiation of myoblasts as well as inhibit apoptosis. Overexpression of circCLTH in vivo suggests that circCLTH is involved in the stimulation of skeletal muscle regeneration. In conclusion, we identified a novel noncoding regulator, circCLTH, that promotes proliferation and differentiation of myoblasts and skeletal muscles.

Milk Exosomal miR-27b Worsen Endoplasmic Reticulum Stress Mediated Colorectal Cancer Cell Death

The relationship between dietary constituents and the onset and prevention of colorectal cancer (CRC) is constantly growing. Recently, the antineoplastic profiles of milk and whey from Mediterranean buffalo (Bubalus bubalis) have been brought to attention. However, to date, compared to cow milk, the potential health benefits of buffalo milk exosome-miRNA are still little explored. In the present study, we profiled the exosomal miRNA from buffalo milk and investigated the possible anticancer effects in CRC cells, HCT116, and HT-29. Results indicated that buffalo milk exosomes contained higher levels of miR-27b, miR-15b, and miR-148a compared to cow milk. Mimic miR-27b transfection in CRC cells induced higher cytotoxic effects (p < 0.01) compared to miR-15b and miR-148a. Moreover, miR-27b overexpression in HCT116 and HT-29 cells (miR-27b+) induced apoptosis, mitochondrial reactive oxygen species (ROS), and lysosome accumulation. Exposure of miR-27b+ cells to the bioactive 3kDa milk extract aggravated the apoptosis rate (p < 0.01), mitochondrial stress (p < 0.01), and advanced endoplasmic reticulum (ER) stress (p < 0.01), via PERK/IRE1/XBP1 and CHOP protein modulation (p < 0.01). Moreover, GSK2606414, the ER-inhibitor (ER-i), decreased the apoptosis phenomenon and XBP1 and CHOP modulation in miR-27b+ cells treated with milk (p < 0.01 vs. miR-27b++Milk), suggesting the ER stress as a cell-death-aggravating mechanism. These results support the in vitro anticancer activity of 3kDa milk extract and unveil the contribution of miR-27b in the promising beneficial effect of buffalo milk in CRC prevention.

Estimation of genetic parameters of growth traits for direct and maternal effects in Murrah buffalo

An understanding of genetic principles and environmental factors affecting the growth traits is essential to implement optimal breeding and selection programs. Early growth is an indicator of future success in production and reproduction status of dairy animals. In this study, a total of 18,989 records of body weight were used to estimate genetic parameters of body weight at birth (BW), 3 months (3BW), 6 months (6BW), 9 months (9BW),12 months (12BW), 18 months (18BW), 24 months (24 BW), 30 months (3BW), and 36 months (36BW) in Murrah buffalo at ICAR-NDRI Karnal, Haryana, for the period 1974-2019. The genetic parameters were estimated using the average information restricted maximum likelihood (AIREML) procedure by excluding or including maternal effects. Six analytical models were fitted in order to optimize the model for each trait. The most appropriate univariate model was selected based on the log likelihood ratio test (LRT). Influencing factors like calf sex, period of birth, season of birth, and dam's parity were investigated. The results showed that the maternal genetic effects, in addition to direct additive effects, were important for unbiased and accurate genetic parameter estimates of growth traits in Murrah buffaloes. Total heritability estimates h²_T1 for BW, 3BW, 6BW, 9BW, 12BW, 18BW, 24BW, 30BW, and 36BW were 0.25, 0.04, 0.14, 0.16, 0.10, 0.15, 0.21, 0.24, and 0.23, respectively. Maternal effect was significant for birth weight and accounted for 13% variation through maternal genetic and 5% variability through maternal permanent environmental effect. Maternal genetic effect was also important for other traits. However, it interfered with the estimates of variance ratios in live weight traits owing to large and negative covariance between direct and maternal genetic effects. Direct genetic correlations between body weight traits were positive and high ranging from 0.10 to 0.94. Results revealed that the Murrah herd has a sizable genetic variability for growth traits and hence, there is sufficient scope for selection for achieving better growth rate if selection in this direction is applied. Owing to higher positive genetic correlation of 6BW with later ages, the scope of indirect selection for optimum growth in later ages can be aimed at by selecting animals for their higher 6-month live weight.

Individual variation in buffalo somatic cell cloning efficiency is related to glycolytic metabolism

Mammalian individuals differ in their somatic cell cloning efficiency, but the mechanisms leading to this variation is poorly understood. Here we found that high cloning efficiency buffalo fetal fibroblasts (BFFs) displayed robust energy metabolism, looser chromatin structure, high H3K9 acetylation and low heterochromatin protein 1α (HP1α) expression. High cloning efficiency BFFs had more H3K9ac regions near to the upstream of glycolysis genes by ChIP-seq, and involved more openness loci related to glycolysis genes through ATAC-seq. The expression of these glycolysis genes was also found to be higher in high cloning efficiency BFFs by qRT-PCR. Two key enzymes of glycolysis, PDKs and LDH, were confirmed to be associated with histone acetylation and chromatin openness of BFFs. Treatment of low cloning efficiency BFFs with PS48 (activator of PDK1) resulted in an increase in the intracellular lactate production and H3K9 acetylation, decrease in histone deacetylase activity and HP1α expression, less condensed chromatin structure and more cloning embryos developing to blastocysts. These results indicate that the cloning efficiency of buffalo somatic cells is associated with their glycolytic metabolism and chromatin structure, and can be improved by increasing glycolytic metabolism.

Genetic parameters of production and reproduction traits of Egyptian buffaloes under subtropical conditions

This research was aimed to investigate the production and reproduction traits and their genetic parameters of Egyptian buffaloes under subtropical environment. Heritability, phenotypic, and genetic parameters were estimated for productive and reproduction traits of first parity This study revealed the averages TMY, 305-dMY, LL, BW, DP, SP, CI, and AFC were 2260, 2150 kg, 271 days, 38.63 kg, 210.50 days, 195.20 days, 586.86 days, and 36.61 months, respectively. The heritability estimates for productive traits of the first lactation were 0.61, 0.52, 0.47, 0.20, and 0.23 for TMY, 305-dMY, LL, DP, and BW, respectively. Heritability estimates for reproductive traits SP, CI, and AFC were 0.07, 0.09, and 0.35, respectively. Genetic correlations of most of the investigated first parity traits were positive, with an exception of DP. The genetic correlation was negative between SP and BW, and between AFC and LL. Phenotypic correlations of the first parity among all investigated traits were positive except between DP and TMY, 305-dMY, LL, BW, and AFC. High and moderate heritability estimates indicated that the opportunity of genetic upgrading of these parameters could be achieved owing to sire selection. Selection for high milk yield will be associated with genetic improvement in lactation period and birth weight.

Comparative Genomic Characterization of Insulin-Like Growth Factor Binding Proteins in Cattle and Buffalo

The somatotropic axis consists of genes associated with economic traits like muscle growth and carcass traits in livestock. Insulin-like growth factor binding proteins (IGFBPs) are the major proteins that play a vital role in the somatotropic axis. The present study performed a genome-wide characterization of IGFBP genes in cattle. Genomic sequences of the IGFBP gene family for different mammals (cattle, buffalo, goat, and sheep) were recovered from the NCBI database. Sequence analyses were performed to investigate cattle's genomic variations in the IGFBP gene family. Phylogenetic analysis, gene structure, motif patterns, and conserved domain analysis (CDA) of the IGFBP family revealed the evolutionarily conserved nature of the IGFBP genes in cattle and other studied species. Physicochemical properties of IGFBP proteins in cattle revealed that most of these proteins are unstable, hydrophilic, thermostable, and acidic. Comparative amino acid analysis revealed variations in all protein sequences with single indels in IGFBP3 and IGFBP6. Mutation analysis revealed only one nonsynonymous mutation D212 > E in the IGFBP6 protein of cattle. A total of 245 nuclear hormone receptor (NHRs) sites were detected, including 139 direct repeats (DR), 65 everted repeats (ER), and 41 inverted repeats (IR). Out of 133 transcription factors (TFs), 10 TFs (AHR, AHRARNT, AP4, CMYB, E47, EGR2, GATA, SP1, and SRF) had differential distribution (P value < 0.05) of putative transcriptional binding sites (TFBS) in cattle compared to buffalo. Synteny analysis revealed the conserved nature of genes between cattle and buffalo. Two gene pairs (IGFBP1/IGFBP3 and IGFBP2/IGFBP5) showed tandem duplication events in cattle and buffalo. This study highlights the functional importance of genomic variation in IGFBP genes and necessitates further investigations better to understand the role and mechanisms of IGFBPs in bovines.

Meprin A1 subunit beta gene polymorphism is associated with the length of post-partum anestrus interval in Murrah buffaloes

Postpartum anestrus interval (PPAI) is the interval between parturition and the first postpartum estrus exhibition in animals. Appearance of both normal and PPA buffaloes under the same farm conditions indicates the role of possible genetic predisposition to PPA. To identify the genetic and non-genetic factors associated with PPA in buffaloes, we collected data on PPAI and other 38 non-genetic variables from 575 Murrah buffaloes in the field conditions and identified the PPA associated non-genetic factors in our previous study. To explore the genetic factors associated with the unexplained variation in PPAI residuals, the present study identified 41 single nucleotide polymorphisms (SNPs) in 13 candidate genes using Sanger sequencing. Exploration of their association with the PPAI residuals of 50 extreme PPA and 50 normal buffaloes identified the significant (P < 0.01) association of the SNP (g.37219977A>G) in the 3'-UTR region of the Meprin A 1 subunit beta (Mep1b) gene with PPAI, which was further validated (P = 0.058) in a large population sample (n = 417). Bioinformatics analysis of the 3'-UTR region has identified three miRNA, bta-miR-2420, bta-miR-2325b and bta-miR-453 that could regulate Igf-1 in the plasma of animals with different genotypes (GG, AG and AA). The higher Igf-1 levels in the GG genotypes than that of AA and AG genotypes of this SNP (g.37219977A>G) further suggest the association of Mep1b gene with PPA condition in Murrah buffaloes. As a result of this study, we propose that buffaloes with protective alleles at this SNP be selected to improve the herd's reproductive efficiency.

Whole-Genome Sequencing Reveals Lignin-Degrading Capacity of a Ligninolytic Bacterium (Bacillus cereus) from Buffalo (Bubalus bubalis) Rumen

The buffalo is an amazing ruminant. Its ability to degrade lignin, which has been recently reported, is most likely due to unique rumen microorganisms with lignin-degradation potential. Our goal was to explore the lignin-degradation potential of ruminal microorganisms, in which ligninolytic enzyme encoding genes were involved to provide ideas for revealing the mechanism of lignin degradation by buffalo. In this study, a bacterium strain identified as Bacillus cereus AH7-7 was isolated from the buffalo (Bubalus bubalis) rumen. After whole-genome sequencing, the results demonstrated that B. cereus AH7-7 had laccase, cytochrome P450 and vanillin alcohol oxidase-encoding genes. Sixty-four genes of B. cereus AH7-7 were involved in multiple aromatic metabolic pathways, such as phenylalanine metabolism and aminobenzoate degradation. A positive reaction resulting in guaiacol medium indicated that laccase secretion from B. cereus AH7-7 increased with time. A biodegradation experiment revealed that a significant reduction in kraft lignin content (25.9%) by B. cereus AH7-7 occurred at the end of 6 days of incubation, which confirmed its lignin-degradation capacity. Overall, this is the first report showing that B. cereus AH7-7 from the buffalo rumen can degrade lignin, and revealing the encoding genes of lignin-degrading enzymes from genome level.

Epidemiology and molecular characterization of Staphylococcus aureus causing bovine mastitis in water buffaloes from the Hazara division of Khyber Pakhtunkhwa, Pakistan

Buffalo represent a major source of milk in Pakistan. However, production is impacted by the disease bovine mastitis. Mastitis causes significant economic losses, with Staphylococcus aureus (S. aureus) being one of its major causative agents. While much work has been done understanding the epidemiology of bovine mastitis in Pakistan, detailed molecular characterization of the associated S. aureus is unavailable. In the current study both the epidemiological and molecular characterization of S. aureus from bovine mastitis in the Hazara division of Pakistan are examined. S. aureus was isolated from 18.41% of the animals, and left quarters more prone to infection (69.6%) than right quarters (30.4%). Sub-clinical mastitis (75.31%) was more prevalent than clinical mastitis (24.69%), with infections evenly distributed amongst the eight districts. Molecular characterization revealed that only 19.6% of the isolates were methicillin-resistant, and four strains types identified, including ST9-t7867-MSSA, ST9-MSSA, ST101-t2078-MSSA, and ST22-t8934-MRSA-IVa. Antiseptic resistance genes were not detected in the isolates, and low levels of antibiotic resistance were also noted, however the methicillin-resistant strains had higher overall antibiotic resistance. This study represents the most complete molecular typing data for S. aureus causing bovine mastitis in the Hazara district of Pakistan, and the country as a whole.

Transcriptional profile of cumulus associated GJA1, PTX3, PRSS35, and SERPINE2 genes with oocytes and embryonic development in water buffalo

BACKGROUND

In the present study, the potential of different groups of cumulus-oocyte complexes (COC's) for in vitro maturation (IVM) and embryonic development was assessed in two groups of COC's of water buffalo. Further, the expression pattern of cumulus-associated GJA1, PTX3, PRSS35, and SERPINE2 genes and their effects on embryonic development was analyzed. Slaughterhouse-derived buffalo COC's were graded into groups A and B. An equal number of 410 COC's were taken in both groups. IVM was carried out using Slaughterhouse-derived buffalo epididymis. A remarkable degree of cumulus expansion was noticed in group A (92.68%) as compared to group B (81.25%) oocytes. On in vitro fertilization (IVF) and embryo culture, group A produced a significantly higher rate of cleavage and blastocyst (92.682 ± 0.7179% and 42.682 ± 0.9683%) as compared to group B (85.365 ± 0.7608% and 31.707 ± 0.9688%). Also, the transcriptional analysis of cumulus-associated genes revealed significantly higher expression in group A as compared to group B.

RESULTS

It was revealed that oocytes having good cumulus mass had a higher developmental potential. Based on differential gene expression of cumulus-associated genes, different quality of COC's, and the resultant embryos after IVF, it was concluded that these genes could be used as a marker for predicting the developmental competence of the oocytes.

CONCLUSION

We concluded that morphologically good quality of COC's had a higher developmental competence, and also the differential expressions of cumulus-associated genes in cumulus cells and embryos. So, we can conclude that these genes could be used as marker genes for predicting the developmental competence of buffalo's oocytes.

Evolutionary and Association Analysis of Buffalo FABP Family Genes Reveal Their Potential Role in Milk Performance

The fatty acid-binding protein (FABP) family gene encode a group of proteins that affect long-chain fatty acid (LCFAs) trafficking and play a crucial function in the regulation of milk fat synthesis. Nevertheless, little is known regarding the identification of members, theevolutionary background, and functional characteristics of FABP genes in buffalo. In this study, in silico analysis was performed to identify the members of FABPs in buffalo. The results revealed that a total of 17 FABP genes were identified. Based on their phylogenetic relationships, these sequences clustered into five groups with similar motif patterns and gene structures. According to positive selection analyses, all duplicated gene pairs containing FABPs in buffalo had Ka/Ks (nonsynonymous/synonymous) ratios that were less than 1, suggesting that they were under purifying selection. Association analysis showed that one SNP in LOC102401361 was found significantly associated with buffalo milk yield. The expression levels of several FABPs in buffalo mammary epithelial cells were regulated by palmitic and stearic acid treatment. The findings of this study provide valuable information for further research on the role of FABPs in regulating buffalo milk synthesis.

Polymorphism of the leptin gene in buffalo breed groups from eastern Amazon

The objective of this work was to characterize genetically some buffalo herds raised in Varzea (VA) and Terra-Firme (upland) (TF) ecosystems through polymorphism of the intron 2 of the leptin gene ( LEP-1620). Two hundred seventy-nine animals from four distinct populations were evaluated using the PCR-RFLP method for LEP-1620 polymorphism (SNP) of the leptin gene with restriction enzyme BsaAI. The animal samples were sorted into 4 groups, according to their breed and environmental origin: Mediterranean TF, Murrah TF, Mediterranean VA and Crossbreed VA. Two alleles (A and G) were detected and their frequencies were analyzed. Allele A frequency ranged from 0.395 (Mediterranean TF) to 0.850 (Murrah TF), with AA genotype ranged from 0.114 (Murrah TF) to 0.700 (Mediterranean TF). The observed and expected heterozygosities ranged from 0.268 (Mediterranean TF) to 0.562 (Murrah TF), and 0.255 (Mediterranean TF) to 0.478 (Murrah TF), respectively. The Hardy-Weinberg probabilities were greater than 0.05. The Crossbred herd in Varzea was the only population with significant inbreeding and the Shannon index ranged between 0.423 (Mediterranean TF) and 0.671 (Murrah TF).

Identification, Molecular Characterization, and Tissue Expression Profiles of Three Smad Genes from Water Buffalo (Bubalus bubalis)

Smads are involved in a variety of biological activities by mediating bone morphogenetic protein (BMP) signals. The full-length coding sequences (CDSs) of buffalo Smads 1, 4, and 5 were isolated and identified through RT-PCR in this study. Their lengths are 1398 bp, 1662 bp, and 1398 bp, respectively. In silico analysis showed that their transcriptional region structures, as well as their amino acid sequences, physicochemical characteristics, motifs, conserved domains, and three-dimensional structures of their encoded proteins are highly consistent with their counterparts in the species of Bovidae. The three Smad proteins are all hydrophilic without the signal peptides and transmembrane regions. Each of them has an MH1 domain and an MH2 domain. A nuclear localization sequence was found in the MH1 domain of buffalo Smads 1 and 5. Prediction showed that the function of the three Smads is mainly protein binding, and they can interact with BMPs and their receptors. The three genes were expressed in all 10 buffalo tissues assayed, and their expression in the mammary gland, gonad, and spleen was relatively high. The results here indicate that the three buffalo Smads may be involved in the transcriptional regulation of genes in a variety of tissues.

Integration of transcriptomics and non-targeted metabolomics reveals the underlying mechanism of follicular atresia in Chinese buffalo

Follicular atresia is a complex physiological process, which results in the waste of follicles and oocytes from the ovary. Elucidating the physiological mechanism of follicular atresia will hopefully reverse the fate of follicles, thereby improve the reproductive efficiency of female animals. However, there are still many gaps to be filled during the follicular atresia process. In this study, we first comprehensively summarized and compared a variety of methods to classify Chinese buffalo follicles with different extent of atresia. Then follicular fluid and granulosa cells from the corresponding follicles with different extent of atresia were collected for non-targeted metabolomics and transcriptomics analysis, respectively. After the detection and analysis of 129 follicles, a reasonable classification standard was formed: on the basis of morphological classification, the relative concentrations of estradiol (E2) and progesterone (PROG) in the follicular fluid were determined, follicles with an estradiol-to-progesterone (E2/PROG) ratio >5 were classified as healthy follicles (HF), 1≤ E2/PROG ≤5 as early atretic follicles (EF) and E2/PROG <1 as late atretic follicles (LF). Correspondingly, follicles with granulosa cells apoptosis rate less than 15 % were divided into HF, 15%-25% were classified as EF and more than 25 % were classified as LF. The integration analysis of non-targeted metabolomics and transcriptomics highlights the following three aspects: (1) Atresia seriously damaged the lipid metabolism homeostasis of follicle, in which PPARγ play important roles. (2) Energy metabolism and nucleotide metabolism of atretic follicles were inhibited. (3) Bilirubin is involved in follicular atresia, and it may be the main force to prevent lipid peroxidation in follicular cells. In summary, results of this study provide new understanding of the molecular mechanisms of Chinese buffalo follicular atresia.

Delineation of temperature-humidity index (THI) as indicator of heat stress in riverine buffaloes (Bubalus bubalis) of a sub-tropical Indian region

The erstwhile developed temperature-humidity index (THI) has been popularly used to indicate heat stress in dairy cattle and often in buffaloes. However, scientific literature suggests differences in thermotolerance and physiological responses to heat stress between cattle and buffalo. Therefore, THI range used to indicate degree of heat stress (mild, moderate, and severe) in cattle should be recalibrated for indicating heat stress in buffaloes. The present study was carried out to delineate THI range to indicate onset and severity of heat stress in buffaloes based on physiological, biochemical, and expression profiling of heat shock response (HSR) genes in animals at different THI. The result indicated early onset of heat stress in buffaloes as compared to cattle. Physiological and biochemical parameters indicated onset of mild signs of heat stress in buffaloes at THI 68-69. Significant deviation in these parameters was again observed at THI range 73-76. At THI 77-80, the physiological and biochemical responses of animals were further intensified indicating extreme alteration in homeostasis. The in vivo expression profiling of HSR genes indicated that members of Hsp70 gene family are expressed in a temporal pattern over different THIs, whereas expressions of Hsf genes were evident during intense heat stress. Overall, the study established that amplitude of heat shock response and THI range for indicating severity of thermal stress for buffaloes are not in unison to cattle. The study also suggests skin temperature of the poll region could be used as non-invasive tool for monitoring heat stress in dairy buffaloes.

Profiling and integrated analysis of whole-transcriptome changes in uterine caruncles of pregnant and non-pregnant buffaloes

Improved reproductive performance in buffaloes can be achieved by understanding the basic mechanism governing the embryonic attachment and feto-maternal communication. Considering this, trascriptomic profiling and integrative analysis of long intergenic non-coding RNAs were carried out in the uterine caruncles of pregnant and non-pregnant buffaloes. Transcriptome data of pregnant and non-pregnant uterine caruncles after quality control was used to perform the analysis. Total of 86 novel lincRNAs expressed in uterine caruncular tissues were identified and characterized. Differential expression analysis revealed that 447 mRNAs and 185 mRNAs were up- and down- regulated, respectively. The number of up- and down- regulated lincRNAs were 114 and 13, respectively. Of the identified 86 novel lincRNAs, six novel lincRNAs were up-regulated in the pregnant uterine caruncles. GO terms (biological process) and PANTHER pathways associated with reproduction and embryogenesis were over-represented in differentially expressed genes. Through miRNA interaction analysis, interactions of 16 differentially expressed lincRNAs with mi-RNAs involved in reproduction were identified. This study has provided a catalogue of differentially expressed genes and novel regions previously unknown to play a significant role in buffalo reproduction. The results from the current study extends the buffalo uterine lncRNAs database and provides candidate regulators for future molecular genetic studies on buffalo uterine physiology to improve the embryo implantation and successful completion of pregnancy.

A genome-wide scan of copy number variants in three Iranian indigenous river buffaloes

BACKGROUND

In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K.

RESULTS

CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database.

CONCLUSIONS

This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo's natural adaptive history together to a recent selection for milk used as primary food source from this species.

Genomic Analysis Revealed a Convergent Evolution of LINE-1 in Coat Color: A Case Study in Water Buffaloes (Bubalus bubalis)

Visible pigmentation phenotypes can be used to explore the regulation of gene expression and the evolution of coat color patterns in animals. Here, we performed whole-genome and RNA sequencing and applied genome-wide association study, comparative population genomics and biological experiments to show that the 2,809-bp-long LINE-1 insertion in the ASIP (agouti signaling protein) gene is the causative mutation for the white coat phenotype in swamp buffalo (Bubalus bubalis). This LINE-1 insertion (3' truncated and containing only 5' UTR) functions as a strong proximal promoter that leads to a 10-fold increase in the transcription of ASIP in white buffalo skin. The 165 bp of 5' UTR transcribed from the LINE-1 is spliced into the first coding exon of ASIP, resulting in a chimeric transcript. The increased expression of ASIP prevents melanocyte maturation, leading to the absence of pigment in white buffalo skin and hairs. Phylogenetic analyses indicate that the white buffalo-specific ASIP allele originated from a recent genetic transposition event in swamp buffalo. Interestingly, as a similar LINE-1 insertion has been identified in the cattle ASIP gene, we discuss the convergent mechanism of coat color evolution in the Bovini tribe.