Genomic clues of association between clinical mastitis and SNPs identified by ddRAD sequencing in Murrah buffaloes

The total milk production of India is 209.96 MT out of which 45% is contributed by the indigenous buffalo and due to their high producing virtue, the prevalence of mastitis is 5-20%. Despite the increasing level of technological advancement, mastitis is still an issue of concern for dairy industry in India as well as across the world. Therefore, the present study aimed to identify the SNPs and associate them with the incidence of clinical mastitis in Murrah buffalo using the ddRAD sequencing approach taking mastitis incidence data of 96 Murrah buffaloes. A total of 246 million quality controlled reads were obtained with an average alignment rate of 99.01% and at a read depth of 10, quality controlled SNPs obtained were 18,056. The logistic regression model was used and a total of seven SNPs were found significantly associated (p < 0.001) with mastitis incidence and seven genes were identified viz., NCBP1, FOXN3, TPK1, XYLT2, CPXM2, HERC1, and OPCML. The majority of them were having tumor suppressing action, related to immunogenetics or glycolytic and energy production. Conclusively, the SNPs identified in this study may be useful for future studies on mastitis incidence in Murrah buffalo and the SNP associations can be further validated.

Salivary cell-free HSD17B1 and HSPA1A transcripts as potential biomarkers for estrus identification in buffaloes (Bubalus bubalis)

Estrus detection is a major problem in buffaloes because of the poor expression of estrus signs leading to low reproductive efficiency. Salivary transcripts analysis is a promising tool to identify biomarkers; therefore, the present study was carried out to evaluate their potential as estrus biomarkers. The levels of HSD17B1, INHBA, HSPA1A, TES transcripts were compared in saliva during estrous cycle stages [early proestrus (day -2, EP), late proestrus (day-1, LP), estrus (E), metestrus (ME) and diestrus (DE)] of cyclic heifers (n = 8) and pluriparous (n = 8) buffaloes by employing quantitative real-time polymerase chain reaction (qRT-PCR). The levels of HSD17B1 (EP/DE 1.46-2.43 fold, LP/DE 2.49-3.06 fold; E/DE 7.21-11.9-fold p < 0.01; ME/D 1.0-1.16 fold) and HSPA1A (EP/DE 0.93-2.39 fold, LP/DE 2.68-3.23 fold; E/DE 8.52-15.18 fold p < 0.01; ME/D 0.86-1.01 fold) were significantly altered during the estrus than other estrous cycle stages in both cyclic heifers and pluriparous buffaloes. Receiver operating characteristic curve analysis revealed the ability of salivary HSD17B1 (AUC 0.96; p < 0.001) and HSPA1A (AUC 0.99; p < 0.01) to differentiate E from other stages of the estrous cycle. Significantly higher levels of HSD17B1 and HSPA1A transcripts in saliva during the estrus phase suggest their biomarkers potential for estrus detection in buffaloes.

Comparative analysis of various methods for estimating expected genetic gain of production and reproduction traits in Murrah buffaloes

The present study was conducted to assess the expected genetic gain for first lactation production and reproduction traits in Murrah buffaloes, in addition to optimization of progenies/sire. Data for period 1971-2020 were used from National Dairy Research Institute. Performance traits considered were 305 days milk yield (305DMY), average daily milk yield (ADMY), peak yield (PY), lactation length (LL), calving to first insemination interval (CFI), days open (DO), and calving interval (CI). Expected ΔG was estimated and compared by three different methods; method I involved heritability and selection differential; method II involved selection intensity, phenotypic standard deviation, and heritability; method III involved estimation of ΔG through four paths of inheritance. Initially, eleven progenies/sire were utilized for assessing expected ΔG by method III, and expected ΔG was found as 34.33, 0.12, 0.12 kg, 2.63, 1.51, 2.74, and 2.80 days/year for 305DMY, ADMY, PY, LL, CFI, DO, and CI, respectively. Additionally, there was a significant increase in expected ΔG on increasing progenies/sire from 6 to 11 while subsequent increase upto 16 had little effect on expected ΔG. These findings will be helpful in formulating breeding strategies worldwide in small buffalo herds to obtain sustainable ΔG in production and reproduction traits.

Genetic parameters for reproductive traits in the Italian Mediterranean buffalo using milk yield as a correlated trait

Until now, the genetic evaluation of the Italian Mediterranean Buffalo has been mainly focused on production traits. However, female fertility affects the efficiency of the dairy industry as it is essential to maintain the profitability of dairy farms. Indeed, the estimation of its genetic component is crucial for its improvement. In this study, 3 measures of buffalo's fertility were analyzed: the age at first calving (AFC), the interval between first and second calving (CIV1), and the interval between second and successive calvings (CIV2_12). Milk yield at 270 d (MY270) was used as a correlated trait. First, genetic parameters were estimated using 7,915 buffalo cows with first calving from 1991 to 2018, then breeding values were calculated from 236,087 buffalo cows. Genetic parameters were estimated by Bayesian inference fitting a multiple-trait animal model using the GIBBS1F90 program, and BLUPF90 was used for estimation of breeding value. The heritability and repeatability estimates of fertility traits were low. The genetic correlations among fertility traits ranged from 0.10 (AFC-CIV1) to 0.92 (CIV1-CIV2_12). Genetic correlation between MY270 and fertility traits was unfavorable, ranging from 0.23 to 0.48. The results from this study can be used as a basis for the future genetic improvement of fertility traits in the Italian Mediterranean Buffaloes.

Genetic variants of CSN1S1, CSN2, CSN3, and BLG genes and their association with dairy production traits in Sahiwal cattle and Nili-Ravi buffaloes

Milk protein genes are associated with milk yield and composition in dairy animals. The present study aimed to identify milk protein genes (CSN1S1, CSN2, CSN3, and BLG) genetic variants and their association with milk yield in Sahiwal cattle and Nili-Ravi buffaloes. One hundred animals from each species were selected to collect blood samples and milk production records. Primers were designed for these milk protein genes for PCR amplification. Sequencing of resultant PCR products revealed a higher number of SNPs (13 vs. 7, 5 vs. 1, and 6 vs. 2) in Sahiwal as compared to Nili-Ravi animals in CSN1S1, CSN2, and CSN3 genes, respectively. However, a single SNP was observed in BLG gene of both species. Association analysis revealed that one SNP in BLG gene of Nili-Ravi was associated (p < 0.05) with 305-day milk yield. Two SNPs at CSN1S1 gene in Sahiwal were associated with dry-period. Similarly, one SNP at CSN1S1 and two SNPs at CSN3 gene showed significant association (p < 0.05) with average calving-interval in Sahiwal while two SNPs in CSN1S1 gene were associated (p < 0.05) with this trait in Nili-Ravi. These SNPs could be helpful as candidate variants for marker-assisted selection in cattle and buffaloes for improvement of lactation performance.

Electroporation-based CRISPR gene editing in adult buffalo fibroblast cells

Electroporation is a widely used method for delivering CRISPR components into cells; however, it presents challenges when applied to difficult-to-transfect cells like adult buffalo fibroblasts. In this study, the ITGB2 gene (encoding the CD18 protein), plays vital for cellular adhesion and immune responses, was selected for editing experiments. To optimize electroporation conditions, we investigated parameters such as electric field strength, pulse duration, plasmid DNA amount, cuvette type, and cell type. The best transfection rates were obtained in a 4 mm gap cuvette with a single 20-millisecond pulse of 300 V using a 10 μg of all-in-one CRISPR plasmid for 106 cells in 100 μL of electroporation buffer. Increasing DNA quantity enhanced transfection rates but compromised cell viability. The 4 mm cuvette gap had high transfection rates than the 2 mm gap, and newborn cells exhibited higher transfection rates than adult cells. We achieved transfection rates of 10-12% with a cell viability of 25-30% for adult fibroblast cells. Subsequently, successfully edited the ITGB2 gene with a 30% editing efficiency, confirmed through various analysis methods, including T7E1 assay, TIDE and ICE analysis, and TA cloning. In conclusion, electroporation conditions reported here can edit buffalo gene(s) for various biotechnological research applications.

Evolutionary analysis of buffalo sterol regulatory element-binding factor (SREBF) family genes and their affection on milk traits

The sterol regulatory element-binding factor (SREBF) genes are a vital group of proteins binding to the sterol regulatory element 1 (SRE-1) regulating the synthesis of fatty acid. Two potential candidate genes (SREBF1 and SREBF2) have been identified as affecting milk traits. This study aims to identify the SREBF family of genes and find candidate markers or SREBF genes influencing lactation production in buffalo. A genome-wide study was performed and identified seven SREBF genes randomly distributed on 7 chromosomes and 24 protein isoforms in buffalos. The SREBF family of genes were also characterized in cattle, goat, sheep and horse, and using these all-protein sequences, a phylogenetic tree was built. The SREBF family genes were homologous between each other in the five livestock. Eight single nucleotide polymorphisms (SNPs) within or near the SREBF genes in the buffalo genome were identified and at least one milk production trait was associated with three of the SNP. The expression of SREBF genes at different lactation stages in buffalo and cattle from published data were compared and the SREBF genes retained a high expression throughout lactation with the trend being the same for buffalo and cattle. These results provide valuable information for clarifying the evolutionary relationship of the SREBF family genes and determining the role of SREBF genes in the regulation of milk production in buffalo.

A long non-coding RNA lnc210 promotes adipogenic differentiation of buffalo intramuscular adipocytes

Intramuscular fat (IMF) content is one of the most significant factors influencing beef quality in terms of tenderness, flavor, and juiciness. Thus, internal factors affecting IMF deposition have received considerable attention for decades. In this study, we demonstrated a long non-coding RNA, lnc210, promoted adipogenic differentiation of buffalo intramuscular adipocytes. lnc210 was rich in adipose tissue and showed increased expression with the adipogenic differentiation of buffalo intramuscular adipocytes. lnc210 was mainly expressed in the nucleus of adipocytes. Full-length lnc210 was obtained by rapid amplification of cDNA ends technology. lnc210 overexpression promoted lipid accumulation by upregulating the mRNA expression of peroxisome proliferator-activated receptor gamma (PPARG) and CCAAT enhancer binding protein alpha (C/EBPα) in buffalo intramuscular adipocytes. These results provide a basis for an in-depth analysis of the role of lnc210 in accelerating IMF deposition in buffaloes.

Comparative RNA isolation methods from fresh ejaculated spermatozoa in Sahiwal cattle (Bos indicus) and Murrah buffalo (Bubalus bubalis) bulls for high quality and enhanced RNA yield

Sperm mRNA transcriptional profiling can be used to evaluate the fertility of breeding bulls. The aim of the study was to compare the modified RNA isolation methods for higher RNA yield and quality from freshly ejaculated sperm of cattle and buffalo bulls. Ten fresh ejaculates from each Sahiwal (n = 10 bulls × 10 ejaculates) and Murrah bulls (n = 10 bulls x 10 ejaculates) were used for RNA isolation. From the recovered live sperm, total sperm RNA was isolated by conventional methods (TRIzol, Double TRIzol), membrane-based methods combined with TRIzol (RNeasy + TRIzol) with the addition of β-mercaptoethanol (BME) and Kit (RNeasy mini) methods in fresh semen. Among different isolation methods; the membrane-based modified methods combined with TRIzol (RNeasy + TRIzol) with the addition of β-mercaptoethanol (BME) resulted significantly (p < .05) higher total RNA quantity (300-340 ng/µL) and better purity in different concentrations of spermatozoa viz., 30-40 million, 70-80 million and 300-400 million sperm. The study concluded that the inclusion of BME to the combined membrane-based methods with somatic cell lysis buffer solution was best for constant increased yield and purity of RNA isolation from Sahiwal cattle and Murrah buffalo bull sperm.

Sustainable multi-trait selection index based on production, reproduction, and health traits for genetic improvement of Murrah buffaloes

Selection of Murrah buffaloes based on production, reproduction, and health traits play an important role in sustainable genetic improvement. With the aim to construct an index comprising of production, reproduction and health traits of Murrah buffaloes, the production and reproduction data during a 50-year span (1971-2020) and for health traits data of 27 years (1994-2020) were analyzed. Various three, four, five, and six trait combinations of production, reproduction, and health traits were explored for construction of performance indices. The income and expenditure method was used to calculate relative economic values (REV) for traits under consideration. A total of 37 performance indices were constructed using a multi-trait combination of production, reproduction, and health traits. Out of 37 indices, the (ADMY-PY-DO-AFC-FF) based index was identified as the best, with 78%, 2%, and 20% weightage given to production, reproduction, and female fertility, respectively. Furthermore, when REVs of traits were increased by 25% and 50%, there was no change observed in the relative importance of traits in the best index, indicating that the index was robust to substantial increases in economic values. Therefore, it is recommended that the above constructed composite selection index can be employed in selection of female buffaloes to achieve sustainable genetic gain.

Identification of novel polymorphism in mammary-derived growth inhibitor gene of water buffalo and its expression analysis in the mammary gland

Mammary-derived growth inhibitor (MDGI), a member of the lipophilic family of fatty acid-binding proteins, plays an important role in the development, regulation, and differentiation of the mammary gland. The aim of the study was to identify polymorphism in the MDGI gene and its expression analysis in the mammary gland at various stages of lactation, in Indian buffalo. Nucleotide sequence analysis of MDGI gene in different breeds of riverine and swamp buffaloes revealed a total of 16 polymorphic sites and one Indel. Different transcription factor binding sites were predicted for buffalo MDGI gene promoter sequence, using online tools and in-silico analysis indicating that the SNPs in this region can impact the gene expression regulation. Phylogenetic analysis exhibited the MDGI of buffalo being closer to other ruminants like cattle, yak, sheep, and goats. Further, the expression analysis revealed that buffalo MDGI being highly expressed in well-developed mammary glands of lactating buffalo as compared to involution/non-lactating and before functional development to start the milk production stage in heifers. Stage-specific variation in expression levels signifies the important functional role of the MDGI gene in mammary gland development and milk production in buffalo, an important dairy species in Southeast Asia.

SIRT6 Knockdown in Buffalo Fetal Fibroblasts Exacerbates Premature Senescence Caused by DNA and Telomere Damage

As a gene with antiaging functions, sirtuin6 (SIRT6) belonging to the sirtuin family plays a vital role in DNA repair, telomerase function, and cellular senescence, as well as maintains epigenomic stability and promotes longevity. However, its role in cell senescence in large animals, such as buffaloes, remains unknown. Fibroblasts are commonly used for somatic reprogramming, and their physiological characteristics affect the efficiency of this process. We aimed to elucidate the role of SIRT6 in cellular senescence and proliferation and analyze its effect on the biological function of buffalo fibroblasts to help improve the efficiency of buffalo somatic cell reprogramming. The expression of SIRT6 and related DNA damage was measured in buffalo fibroblasts obtained at different developmental stages (in the fetus and at 3 and 10 years of age), and the effect of SIRT6 knockdown on the senescence of buffalo fetal fibroblast was investigated. An inverse relationship was observed between SIRT6 expression and senescence in buffalo fibroblasts obtained from animals of various ages. This was accompanied by decreased cell growth, viability, and increased DNA damage. Short hairpin RNA-mediated SIRT6 knockdown accelerated the senescence of buffalo fetal fibroblasts. It blocked the cell cycle during in vitro cell culture, which further enhanced DNA damage, particularly with respect to the telomeres. Collectively, our findings suggest that SIRT6 expression was closely associated with buffalo senescence in fibroblasts. These findings serve as a foundation to better understand the cellular functions of SIRT6 and also aid in selecting donor cells for buffalo somatic cell reprogramming.

Identification of core genes affecting IMF deposition in bovine

Intramuscular fat (IMF) content is an important economic factor in beef production. However, knowledge on the key factors controlling bovine IMF is limited. In this study, using weighted gene co-expression network analysis (WGCNA), nine modules were identified and the number of transcripts in these modules ranged from 36 to 3191. Two modules were found to be significantly associated with fat deposition and three genes (TCAP, MYH7, and TNNC1) were further identified by Protein-protein interaction (PPI), which may be the hub genes regulating bovine IMF deposition. In addition, considering the genetic variation, the PCK1 gene was found by functional enrichment analysis of overlapping genes, which was previously reported to be involved in IMF deposition. We noted that the core promoter region of buffalo PCK1 binds to transcription factors involved in lipid metabolism while cattle PCK1 binds transcription factors involved in muscle development. The results suggest that PCK1 participated in IMF deposition of buffalo and cattle in different ways. In summary, gene expression networks and new candidate genes associated with IMF deposition identified in this study. This would lay the foundation for further research into the molecular regulatory mechanisms underlying bovine IMF deposition.

Genetic analysis of milk and milk composition traits in Murrah buffaloes using Bayesian inference

Accurate and unbiased assessment of genetic parameters of milk and milk composition traits play an important role in formulating breeding program for genetic improvement of Murrah buffaloes. In this study, data spread over 28 years were utilized to estimate genetic parameters of traits viz., 305 d milk yield (305MY), 305 d fat yield (305FY), 305 d solid not fat yield (305SNFY), milk fat percentage (fat%) and solid not fat percentage (SNF) percentage (SNF%) in Murrah buffaloes kept at ICAR-National Dairy Research Institute, Karnal. Bayesian multiple-trait analysis was done using animal model and Gibbs sampling to estimate (co)variance components. Posterior means of heritability and posterior standard deviation for 305MY, 305FY, 305SNFY, fat% and SNF% were 0.18 ± 0.05, 0.17 ± 0.05, 0.18 ± 0.05, 0.07 ± 0.03 and 0.15 ± 0.06 and posterior means of repeatability estimates along with posterior standard deviation for corresponding traits were 0.33 ± 0.04, 0.32 ± 0.04, 0.33 ± 0.04, 0.14 ± 0.02 and 0.30 ± 0.04, respectively. Estimates of genetic correlation varied from -0.080 (305MY and fat %) to 0.999 (305MY and 305SNFY). Permanent environmental correlations varied from -0.060 (305MY and SNF%) to 0.999 (305FY and 305SNFY). This study indicated that all considered traits except fat% have ample genetic variability which can be exploited for selection and genetic improvement of Murrah buffaloes.

Sequence characterization and comparative expression profile of buffalo WNT10B gene in adult and fetal tissues

In this study, Wingless-type MMTV (mouse mammary tumor virus) integration site family member (WNT10B) gene was sequence characterized in the Indian water buffalo. Sequence analysis revealed an open reading frame of 1176 nucleotides in buffalo, encoding 391 amino acids long protein. Nineteen nucleotide variations were observed between cattle and buffalo resulting in six amino acid changes. Phylogenetic analysis showed the clustering of ruminant species together. Real-time expression analysis of WNT10B in tissues collected from different organs of fetal and adult buffalo, revealed, the gene being abundantly expressed in the rumen and liver of the fetus. The fetal ovary, heart, kidney, lung, testis and mammary gland showed moderate expression, while in adult tissues, expression was high in the ovary, testis, brain, kidney, small intestine and liver, whereas lower expression was observed in the adult rumen. Significant differences in WNT10B expression levels were found for the brain, small intestine, testes, kidney, heart, rumen, and ovary when adult and fetal tissues were compared. A moderate level of genetic variation was found between cattle and buffalo WNT10B and expression patterns in a variety of tissues in adult buffalo implies that in addition to possible roles in adipogenesis and hematopoiesis, the WNT10B gene might be playing a significant role in other regulatory pathways as well.

Genetic characterization and diversity assessment in 'Bhangor' indigenous swamp buffalo population using heterologous microsatellite markers

'Bhangor' newly identified swamp buffalo population from North East Indian, was characterized using microsatellite markers. Genomic DNA was isolated from blood samples of 76 unrelated animals, 15 microsatellite markers (CSSM33, BM1818, CSRM60, HEL13, ILSTS019, ILSTS025, ILSTS028, ILSTS029, ILSTS033, ILSTS036, ILSTS056, ILSTS058, ILSTS061, ILSTS089 and ETH003) were found to be highly polymorphic in the population of the selected markers. A total of 114 alleles were observed, which ranged from 3 in CSRM60 and ILSTS025 locus to 12 in ILSTS056 and ILSTS061. The mean effective number of alleles across all polymorphic loci was found to be 3.76. The overall mean expected heterozygosity and unbiased expected heterozygosity values were 0.67 and 0.68, ranging from 0.067 (ILSTS025) to 0.85 (ILSTS058) and 0.068 (ILSTS025) to 0.86 (ILSTS058), respectively. Within the population, the inbreeding estimates (FIS) ranged between -0.4352 and 0.804, with an average FIS of 0.114 ± 0.033. The outcome for infinite allele model (IAM), two-phase model (TPM) and test for mode shift revealed the absence of any recent bottleneck in the investigated buffalo population. The population was found to be in optimum diversity based on polymorphic microsatellite markers. With fast changing agro-climatic conditions; there is an urgent need to characterize the nondescript livestock populations.

Signatures of selection in riverine buffalo populations revealed by genome-wide SNP data

The detection of selection signatures assists in understanding domestication, evolution, and the identification of genomic regions related to adaptation and production traits in buffaloes. The emergence of high-throughput technologies like Next Generation Sequencing and SNP genotyping had expanded our ability to detect these signatures of selection. In this study, we sought to identify signatures of selection in five buffalo populations (Brazilian Murrah, Bulgarian Murrah, Indian Murrah, Nili-Ravi, and Kundi) using Axiom Buffalo 90 K Genotyping Array data. Using seven different methodologies (Tajima's D, CLR, ROH, iHS, FST, FLK and hapFLK), we identified selection signatures in 374 genomic regions, spanning a total of 381 genes and 350 quantitative trait loci (QTLs). Among these, several candidate genes were associated with QTLs for milk production, reproduction, growth and carcass traits. The genes and QTLs reported in this study provide insight into selection signals shaping the genome of buffalo breeds. Our findings can aid in further genomic association studies, genomic prediction, and the implementation of breeding programmes in Indian buffaloes.

Functional roles for AGPAT6 in milk fat synthesis of buffalo mammary epithelial cells

AGPAT6 plays a crucial role in the triglyceride (TG) synthesis pathway in mammals. However, its roles in buffalo lactation remain unknown. Therefore, we investigated the functional roles of AGPAT6 in milk fat synthesis by transfecting overexpression and lentivirus interference vectors in buffalo mammary epithelial cells (BuMECs) in vitro. AGPAT6 overexpression in BuMECs significantly enhanced the mRNA expression of FABP4, SLC27A6, ACSL1, DGAT1, DGAT2, LPIN1, INSIG1, CEBPA and SREBF1 genes, and significantly reduced that of XDH, CPT1A, LIPE, INSIG2 and PPARGC1A, but has no significant influence to the mRNA abundance of FABP3, GPAM, PPARG and SREBF2. However, the interference with AGPAT6, the mRNA expression of FABP4, SLC27A6, ACSL1, DGAT1, DGAT2, INSIG1, CEBPA, SREBF1, XDH, CPT1A, LIPE, INSIG2 and PPARGC1A genes in BuMECs changed contrary to the overexpression experiment, and that of GPAM, PPARG and SREBF2 also did not change significantly, but the expression of FABP3 was significantly decreased. In addition, the overexpression/interference of AGPAT6 gene significantly increased/decreased TG content in BuMECs. The results here indicate that AGPAT6 gene is involved in TG synthesis in BuMECs, and affects the expression of major genes associated with FA transport and activation, TG synthesis and transcription regulation, FA oxidation and TG degradation during the lipogenesis of milk.

Role of Different Members of the AGPAT Gene Family in Milk Fat Synthesis in Bubalus bubalis

During triacylglycerol synthesis, the acylglycerol-3-phosphate acyltransferase (AGPAT) family catalyzes the conversion of lysophosphatidic acid to phosphatidic acid and the acylation of sn-2 fatty acids. However, the catalytic activity of different AGPAT members is different. Therefore, this study aimed to investigate the mechanism through which different AGPATs affect the efficiency of TAG synthesis and fatty acid composition. The conservation of amino acid sequences and protein domains of the AGPAT family was analyzed, and the functions of AGPAT1, AGPAT3, and AGPAT4 genes in buffalo mammary epithelial cells (BMECs) were studied using RNA interference and gene overexpression. Prediction of the protein tertiary structure of the AGPAT family demonstrated that four conservative motifs (motif1, motif2, motif3, and motif6) formed a hydrophobic pocket in AGPAT proteins, except AGPAT6. According to cytological studies, AGPAT1, AGPAT3, and AGPAT4 were found to promote the synthesis and fatty acid compositions of triacylglycerol, especially UFA compositions of triacylglycerol, by regulating ACSL1, FASN, GPAM, DGAT2, and PPARG gene expression. This study provides new insights into the role of different AGPAT gene family members involved in TAG synthesis, and a reference for improving the fatty acid composition of milk.

Clonal spread of non-O157 Shiga toxigenic Escherichia coli O21:H25 in raw water buffalo milks

AIMS

This study was conducted to investigate the presence of Shiga toxin-producing O157 and non-O157 E. coli in raw water buffalo milk, as well as to determine the virulence gene profiles, phylogroups, sequence types, and serotypes of the isolated strains.

METHODS AND RESULTS

A total of 200 hand-milked raw water buffalo milk samples were collected from 200 different water buffaloes over a period of three months from 20 different farms. Isolation of STEC was performed using CHROMagar STEC. Presence of stx1, stx2, and eaeA genes were investigated by mPCR. Phylogroups and sequence types of E. coli strains were determined by Clermont phylotyping and MLST. Serotyping was performed using PCR or WGS. According to the results, two milk samples obtained from two different farms were found as STEC-positive. All Stx-positive E. coli isolates belonged to phylogenetic group A and were assigned to ST10. WGS results indicated that serotype of two isolates was O21:H25 and average nucleotide identity was detected at 99.99%. Thirteen additional registered E. coli O21:H25 assembled WGS data were obtained from EnteroBase and a phylogenetic tree was constructed.

CONCLUSIONS

With this study, the presence of stx2 harboring E. coli O21:H25 in milk was identified for the first time. Although the identified serotype is considered a non-pathogen seropathotype, we conclude it could play an important role in the environmental circulation of Stx-phages and consequently contribute to the emergence of new STEC-related outbreaks.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds

BACKGROUND

Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated.

RESULTS

Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds.

CONCLUSION

The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Circulatory extracellular vesicle derived miR-195-5p promotes cellular apoptosis and suppresses cell proliferation in the buffalo endometrial primary cell culture

In pregnant animals, communication between the mother and conceptus occurs via extracellular vesicles (EVs) that carry several biomolecules such as nucleic acids (miRNAs, mRNAs), proteins, and lipids. At the time of implantation, the endometrium undergoes several morphological and physiological changes, such as angiogenesis, apoptosis, and cell proliferation regulation at the implantation site, to attain a receptive state. This study was conducted to detect pregnancy-specific miRNAs derived from extracellular vesicles in the systemic circulation of Bubalus bubalis (water buffalo) and to assess their functional significance in the modulation of endometrial primary cells. The extracellular vesicles were isolated from the blood plasma using a precipitation-based method and further characterized by various methods such as Differential light scattering, Nanoparticle tracking assay, Western blot, and transmission electron microscopy. The relative expression of the selected extracellular vesicles associated miRNAs (EV-miRNA) at different intervals (days 15, 19, 25, and 30) post artificial insemination (AI) was analyzed using RT-qPCR, and expression of miR-195-5p was found to be significantly higher (P < 0.01) in pregnant animals on day 19 post AI (implantation window) as compared to day 15 post AI. The elevated expression might indicate the involvement of this miRNA in the maternal-conceptus cross-talk occurring during the implantation period. The KEGG pathway enrichment and Gene Ontology analyses of the miR-195-5p target genes revealed that these were mostly involved in the PI3-Akt, MAPK, cell cycle, ubiquitin-mediated proteolysis, and mTOR signaling pathways, which are related to the regulation of cell proliferation. Transfecting the in vitro cultured cells with miR-195-5p mimic significantly suppressed (P < 0.05) the expression of its target genes such as YWHAQ, CDC27, AKT-3, FGF-7, MAPK8, SGK1, VEGFA, CACAND1, CUL2, MKNK1, and CACAN2D1. Furthermore, the downregulation of the miR-195-5p target genes was positively correlated with a significant increase in the apoptotic rate and a decrease in the proliferation. In conclusion, the current findings provide vital information on the presence of EV miR-195-5p in maternal circulation during the implantation window indicating its important role in the modulation of buffalo endometrium epithelial cells via promoting cell death. Altogether, the milieu of miR-195-5p may serve as a novel and potential molecular factor facilitating the implantation of the early embryo during the establishment of pregnancy in buffaloes. Thus, miR-195-5p may be identified as a unique circulatory EV biomarker related to establishing pregnancy in buffaloes as early as day 19 post-AI.

Gene expression of pregnancy-associated glycoproteins-1 (PAG-1), interferon-tau (IFNt) and interferon stimulated genes (ISGs) as diagnostic and prognostic markers of maternal-fetal cellular interaction in buffalo cows

The aim of this study was to determine the presence of Pregnancy-associated glycoprotein -1 (PAG-1) mRNA expression in the maternal circulation of pregnant buffaloes during the early stage of pregnancy. Contemporaneously, the mRNA expression levels of Interferon-tau (IFNt) and some Interferon stimulated genes (ISGs) (interferon stimulated gene 15 ubiquitin-like modifier interferon, ISG15; Mixoviruses resistance 1 and 2, MX1 and MX2; 2',5'-oligoadenylate synthase 1,OAS1) were evaluated in order to expand our knowledge of the molecular processes involved in the early stages of pregnancy and to identify potential biomarkers of maternal-fetal cellular interaction in buffalo. The study was conducted on 38 synchronized and artificially inseminated buffalo cows (d 0), divided ex post into 3 groups: Pregnant (n = 17), Non-pregnant (n = 15) and Embryo mortality (n = 6). Blood samples were collected on d 14, 19, 28 and 40 after artificial insemination (AI) for peripheral blood mononuclear cells (PBMCs) isolation. Expression levels of mRNA of PAG-1, IFNt, ISG15. MX1, MX2 and OAS1 were measured using RT-qPCR. No significant changes were observed in IFNt and PAG gene expressions between groups, while significant differences (p < 0.001) were found for ISG15, MX1, MX2, and OAS1. Pairwise comparisons revealed that the differences between groups occurred on days 19 and 28 post-AI. ISG15 proved to have the best diagnostic performance for distinguishing between pregnant animals and animals that experienced embryo mortality with the ROC analysis. According to the results of the univariate analyses, day 19 was identified as the most indicative to discriminate between groups while the most reliable genes for this differentiation were ISG15, MX1 and MX2. MX2 proved to be the best gene for discriminating pregnant buffaloes using the discriminant analysis, while MX1 was the gene that best predicted embryo mortality. Our results showed that among PAG-1, IFNt and ISGs expression as diagnostic and prognostic markers of maternal-fetal cellular interaction in buffalo cows, ISGs proved to be the best peripheral biomarkers for predicting pregnancy and embryonic mortality during the peri-implantation period. These insights into the mechanisms behind maternal-fetal interaction and the development of a method for the early detection of embryo distress may enable us to implement effective strategies to support embryo survival.

Purification and characterization of cysteine protease of Sarcocystis fusiformis from infected Egyptian water buffaloes

Sarcocystis spp. infects water buffaloes (Bubalus bubalis) causing sarcocystosis. In the present study, Sarcocystis fusiformis was recognized in Egyptian water buffaloes based on histological observation and molecular analysis of internal transcribed spacer 1 (ITS1), 18S ribosomal RNA (18S rRNA) and cytochrome c oxidase subunit I (COX-1) gene fragments. Chemotherapy and vaccines against Sarcocystis spp. could potentially target proteases because they may play a crucial role in the infection. Cysteine proteases are multifunctional enzymes involved in vital metabolic processes. However, the involvement of proteases in S. fusiform infection has not yet been characterized. Here, the purification and study on some biochemical properties of protease isolated from cysts of S. fusiform were carried out. Protease with a molecular weight of 100 kDa was purified. LC-MS/MS analyzed the protein sequence of purified protease and the data suggested that the enzyme might be related to the cysteine protease. The purified protease exhibited maximum activity at pH 6 and a temperature of 50 °C. The Michaelis-Menten constant (Km), the maximum velocity (Vmax), and the turnover number (Kcat) were determined. The complete inhibition effect of cysteine inhibitors indicated that the purified enzyme is a cysteine protease. The results suggested that S. fusiform proteolytic enzyme may be necessary for parasite survival in water buffaloes by digesting host tissues. Therefore, cysteine protease could be a suitable target for vaccinations.

Genetic analysis for semen quality traits in buffalo bulls

This study was attempted to estimate the genetic parameters of semen quality traits in buffalo bulls. The study data consisted of 10975 ejaculates from 45 Murrah buffalo bulls (aged 24-72 months) used for breeding program during year 2010 to 2020. Semen quality traits (ejaculate volume, concentration of sperm, mass activity, initial and post-thaw motility, number of sperms per ejaculate, motile sperm number and discard rates) were studied. It was observed that average ejaculate volume was 2.82 ± 1.45 mL with mean concentration of 1040.12 ± 523.26 million/mL. Higher heritability was observed for number of sperms per ejaculate, number of motile sperm and sperm concentration. Significant phenotypic correlation was obtained between volume and number of sperms per ejaculate as well as volume and number of motile sperms. Likewise, significant phenotypic correlation was evident between sperm concentration with sperm number per ejaculate. Highest phenotypic correlation was obtained between sperm count per ejaculate and motile sperm count. Estimated genetic trends showed significant change in volume and motile sperm number. In conclusion, this study ascertains that genetic parameters of semen traits can be considered during the selection of buffalo bulls in breeding program.