Whole-Genome Sequencing and Characterization of Buffalo Genetic Resources: Recent Advances and Future Challenges

Omics approaches to understand impact of heat stress on semen quality and fertility in bovines

Recent surge in global climate change poses unprecedented challenges to traditional livestock breeding methods. In tropical countries, where large bovines are vital source of meat and milk, the escalating global warming significantly impact their productivity. Central to this challenge is the resilient ability of the animals, which directly influences the productivity. Artificial insemination (AI) programs, a cornerstone in modern livestock production, have also been severely hampered by climate change-induced heat stress, affecting the quality semen production. This stress not only affects the blood biochemical profiles of the animals but also their testicular physiology, leading to the issues such as low-quality semen with compromised freezability and fertility. Researchers have identified specific markers, including single nucleotide polymorphisms (SNPs), copy number variations (CNV), and epigenetic signatures like histone modifications, DNA methylation, and noncoding RNAs, influencing the semen quality in livestock species. Furthermore, the issue of heat stress has been addressed in a very precise way, and biomarkers have been identified, which can be integrated into the breeding programme to keep up the sire summary. Transcriptomic studies have further illuminated the temporal expression patterns of genes related to sperm quality during heat stress, pinpointing candidate genes for further exploration. This review comprehensively summarizes the progress made in understanding the intricacies of sperm biology in bovines, with a specific focus on cattle and buffalo delving into a spectrum of changes, from biochemical shifts to profound cellular alterations, including genomic, transcriptomic, and epigenetic modifications.

BuffExDb: web-based tissue-specific gene expression resource for breeding and conservation programmes in Bubalus bubalis

Amidst the global challenge of extreme poverty, the livestock sector can significantly contribute to global sustainable development goals by enhancing resilience, smallholder productivity, and market participation. The Indian livestock sector is one of the largest in the world with a total livestock population of 535.82 million, ∼10.7% of the world's livestock population. Buffalo (Bubalus bubalis) holds significant importance in India and other Asian countries, notably contributing to their economies by surpassing cattle in milk production and providing various valuable products. The limited availability of genomic and transcriptomic resources for buffaloes hinders the efforts to enhance their traits for increased milk and meat production. To address this gap, this study adopted the state-of-the-art bioinformatics tools to analyse 2429 transcriptomes representing 438 BioSamples from 23 BioProjects obtained from a public domain database, representing 76 different types of tissues and cell types from all major organ systems in buffalo species (river and swamp). The outcome of this exhaustive genomic data led to the development of a relational buffalo expression database based on a three-tier architecture named as BuffExDb (http://46.202.167.198/buffex/). The user-friendliness and flexibilities in retrieval of tissue-specific genes (TSGs) and their functional annotation are the major characteristics of BuffExDb. This is the first of its kind that offers an effortlessly navigable and filterable database, enabling users to examine and visualize the expression levels of each tissue across multiple samples, simultaneously. It also provides the Tau score parameter for the identification of TSGs along with their essential roles in tissue development, maintenance, and function as observed through the enrichment test for gene ontologies. The exhaustive outcome of this work would pave the way for the biological, functional, and evolutionary studies for easy access. This prior information based on tissue-specific mechanisms can be used for future genomic research, especially in association studies in endeavour of enhanced buffalo breeding and conservation programmes. Database URL: http://46.202.167.198/buffex/.

A multi-omics database of buffaloes from Yangtze valley reveals diversity of water buffalo (Bubalus bubalis)

Asian water buffalo (Bubalus bubalis) is the fundamental livestock resource for local rural populations and holds a promising prospect of their milk and meat. Xuyi mountain (XYM) and Haizi (HZ) buffaloes from Yangtze valley comprises of species diversity of Asian water buffaloes. Current multi-omics enables identification of causal genes and elucidation of genetic regulatory mechanisms underlying complex traits in buffaloes. Here, we conducted the integrated analysis of metabolome and metagenome of rumen fluid, transcriptome and metabolome of blood, and whole genome sequence data from XYM (n = 7) and HZ (n = 10) male buffaloes. Our results revealed the apparent diversity of multi-layer omics profiles between two buffalo species. The built-up multi-omics database supports the discoveries of diversity in Asian water buffalo and potentially serves valuable resources for studying causal regulatory variants and their mechanisms.

The genome landscape of the Xinglong buffalo

BACKGROUND

Xinglong buffalo, as an indigenous breed in Hainan province of China, possesses characteristics such as high humidity tolerance, disease resistance and high reproductive capacity. Combined with whole genome sequencing technology, comprehensive investigation can be undertaken to elucidate the genomic characteristics, functions and genetic variation of Xinglong buffalo population.

RESULTS

Xinglong buffalo has the highest genetic diversity, lowest runs of homozygosity average length, and fasted decay of linkage disequilibrium in our study population. Phylogenetic tree results revealed that Xinglong buffalo was gathered together with Fuzhong buffalo firstly. The population genetic structure analysis indicates that at K = 3, the Xinglong buffalo for the first time showed a distinct ancestral origin from other water buffalo. Furthermore, compared to different populations, candidate genes displaying significantly distinct patterns of single nucleotide polymorphisms (SNPs) (e.g., RYR2, COX15, PCDH9, DTWD2, FCRL5) distribution have been identified in the Xinglong buffalo.

CONCLUSIONS

Based on the whole genome sequencing data, this study identified a substantial number of SNPs and assessed the genetic diversity and selection signatures within the Xinglong buffalo population. These results contribute to understanding the genomic characteristics of Xinglong buffalo and their genetic evolutionary status. However, the practical significance of these signatures for genetic enhancement still requires confirmation through additional samples and further experimental validation.

In Silico Analysis: Molecular Characterization and Evolutionary Study of CLCN Gene Family in Buffalo

Chloride channels (ClCs) have received global interest due to their significant role in the regulation of ion homeostasis, fluid transport, and electrical excitability of tissues and organs in different mammals and contributing to various functions, such as neuronal signaling, muscle contraction, and regulating the electrolytes' balance in kidneys and other organs. In order to define the chloride voltage-gated channel (CLCN) gene family in buffalo, this study used in silico analyses to examine physicochemical properties, evolutionary patterns, and genome-wide identification. We identified eight CLCN genes in buffalo. The ProtParam tool analysis identified a number of important physicochemical properties of these proteins, including hydrophilicity, thermostability, in vitro instability, and basic nature. Based on their evolutionary relationships, a phylogenetic analysis divided the eight discovered genes into three subfamilies. Furthermore, a gene structure analysis, motif patterns, and conserved domains using TBtool demonstrated the significant conservation of this gene family among selected species over the course of evolution. A comparative amino acid analysis using ClustalW revealed similarities and differences between buffalo and cattle CLCN proteins. Three duplicated gene pairs were identified, all of which were segmental duplications except for CLCN4-CLCN5, which was a tandem duplication in buffalo. For each gene pair, the Ka/Ks test ratio findings showed that none of the ratios was more than one, indicating that these proteins were likely subject to positive selection. A synteny analysis confirmed a conserved pattern of genomic blocks between buffalo and cattle. Transcriptional control in cells relies on the binding of transcription factors to specific sites in the genome. The number of transcription factor binding sites (TFBSs) was higher in cattle compared to buffalo. Five main recombination breakpoints were identified at various places in the recombination analysis. The outcomes of our study provide new knowledge about the CLCN gene family in buffalo and open the door for further research on candidate genes in vertebrates through genome-wide studies.

Genome-wide comparative analyses highlight selection signatures underlying saline adaptation in Chilika buffalo

Chilika, a native buffalo breed of the Eastern coast of India, is mainly distributed around the Chilika brackish water lake connected with the Bay of Bengal Sea. This breed possesses a unique ability to delve deep into the salty water of the lake and stay there to feed on local vegetation of saline nature. Adaptation to salinity is a genetic phenomenon; however, the genetic basis underlying salinity tolerance is still limited in animals, specifically in livestock. The present study explores the genetic evolution that unveils the Chilika buffalo's adaptation to the harsh saline habitat, including both water and food systems. For this study, whole genome resequencing data on 18 Chilika buffalo and for comparison 10 Murrah buffalo of normal habitat were generated. For identification of selection sweeps, intrapopulation and interpopulation statistics were used. A total of 709, 309, 468, and 354 genes were detected to possess selection sweeps in Chilika buffalo using the nucleotide diversity (θπ), Tajima's D, nucleotide diversity ratio (θπ-ratio), and FST methods, respectively. Further analysis revealed a total of 23 genes including EXOC6B, VPS8, LYPD1, VPS35, CAMKMT, NCKAP5, COMMD1, myosin light chain kinase 3 (MYLK3), and B3GNT2 were found to be common by all the methods. Furthermore, functional annotation study of identified genes provided pathways such as MAPK signaling, renin secretion, endocytosis, oxytocin signaling pathway, etc. Gene network analysis enlists that hub genes provide insights into their interactions with each other. In conclusion, this study has highlighted the genetic basis underlying the local adaptive function of Chilika buffalo under saline environment.NEW & NOTEWORTHY Indian Chilika buffaloes are being maintained on extensive grazing system and have a unique ability to convert local salty vegetation into valuable human food. However, adaptability to saline habitat of Chilika buffalo has not been explored to date. Here, we identified genes and biological pathways involved, such as MAPK signaling, renin secretion, endocytosis, and oxytocin signaling pathway, underlying adaptability of Chilika buffalo to saline environment. This investigation shed light on the mechanisms underlying the buffalo's resilience in its native surroundings.

RNA sequencing and gene co-expression network of in vitro matured oocytes and blastocysts of buffalo

In reproductive technologies, uncovering the molecular aspects of oocyte and embryo competence under different conditions is crucial for refining protocols and enhancing efficiency. RNA-seq generates high-throughput data and provides transcriptomes that can undergo additional computational analyses. This study presented the transcriptomic profiles of in vitro matured oocytes and blastocysts produced in vitro from buffalo crossbred (Bubalus bubalis), coupled with gene co-expression and module preservation analysis. Cumulus Oophorus Complexes, obtained from slaughterhouse-derived ovaries, were subjected to in vitro maturation to yield metaphase II oocytes (616) or followed in vitro fertilization and culture to yield blastocysts for sequencing (526). Oocyte maturation (72%, ±3.34 sd) and embryo development (21.3%, ±4.18 sd) rates were obtained from three in vitro embryo production routines following standard protocols. Sequencing of 410 metaphase II oocytes and 70 hatched blastocysts (grade 1 and 2) identified a total of 13,976 genes, with 62% being ubiquitously expressed (8,649). Among them, the differentially expressed genes (4,153) and the strongly variable genes with the higher expression (fold-change above 11) were highlighted in oocytes (BMP15, UCHL1, WEE1, NLRPs, KPNA7, ZP2, and ZP4) and blastocysts (APOA1, KRT18, ANXA2, S100A14, SLC34A2, PRSS8 and ANXA2) as representative indicators of molecular quality. Additionally, genes exclusively found in oocytes (224) and blastocysts (2,200) with specific biological functions were identified. Gene co-expression network and module preservation analysis revealed strong preservation of functional modules related to exosome components, steroid metabolism, cell proliferation, and morphogenesis. However, cell cycle and amino acid transport modules exhibited weak preservation, which may reflect differences in embryo development kinetics and the activation of cell signaling pathways between buffalo and bovine. This comprehensive transcriptomic profile serves as a valuable resource for assessing the molecular quality of buffalo oocytes and embryos in future in vitro embryo production assays.

Susceptibility of Mediterranean Buffalo (Bubalus bubalis) following Experimental Infection with Lumpy Skin Disease Virus

Lumpy skin disease (LSD) is a viral disease of cattle and water buffalo characterized by cutaneous nodules, biphasic fever, and lymphadenitis. LSD is endemic in Africa and the Middle East but has spread to different Asian countries in recent years. The disease is well characterized in cattle while little is known about the disease in buffaloes in which no experimental studies have been conducted. Six buffaloes and two cattle were inoculated with an Albanian LSD virus (LSDV) field strain and clinically monitored for 42 days. Only two buffaloes showed fever, skin nodules, and lymphadenitis. All samples collected (blood, swabs, biopsies, and organs) were tested in real-time PCR and were negative. Between day 39 and day 42 after inoculation, anti-LSDV antibodies were detected in three buffaloes by ELISA, but all sera were negative by virus neutralization test (VNT). Cattle showed severe clinical signs, viremia, virus shedding proven by positive real-time PCR results, and seroconversion confirmed by both ELISA and VNT. Clinical findings suggest that susceptibility in buffaloes is limited compared to in cattle once experimentally infected with LSDV. Virological results support the hypothesis of buffalo resistance to LSD and its role as an accidental non-adapted host. This study highlights that the sensitivity of ELISA and VNT may differ between animal species and further studies are needed to investigate the epidemiological role of water buffalo.

A New Method to Detect Buffalo Mastitis Using Udder Ultrasonography Based on Deep Learning Network

Mastitis is one of the most predominant diseases with a negative impact on ranch products worldwide. It reduces milk production, damages milk quality, increases treatment costs, and even leads to the premature elimination of animals. In addition, failure to take effective measures in time will lead to widespread disease. The key to reducing the losses caused by mastitis lies in the early detection of the disease. The application of deep learning with powerful feature extraction capability in the medical field is receiving increasing attention. The main purpose of this study was to establish a deep learning network for buffalo quarter-level mastitis detection based on 3054 ultrasound images of udders from 271 buffaloes. Two data sets were generated with thresholds of somatic cell count (SCC) set as 2 × 105 cells/mL and 4 × 105 cells/mL, respectively. The udders with SCCs less than the threshold value were defined as healthy udders, and otherwise as mastitis-stricken udders. A total of 3054 udder ultrasound images were randomly divided into a training set (70%), a validation set (15%), and a test set (15%). We used the EfficientNet_b3 model with powerful learning capabilities in combination with the convolutional block attention module (CBAM) to train the mastitis detection model. To solve the problem of sample category imbalance, the PolyLoss module was used as the loss function. The training set and validation set were used to develop the mastitis detection model, and the test set was used to evaluate the network's performance. The results showed that, when the SCC threshold was 2 × 105 cells/mL, our established network exhibited an accuracy of 70.02%, a specificity of 77.93%, a sensitivity of 63.11%, and an area under the receiver operating characteristics curve (AUC) of 0.77 on the test set. The classification effect of the model was better when the SCC threshold was 4 × 105 cells/mL than when the SCC threshold was 2 × 105 cells/mL. Therefore, when SCC ≥ 4 × 105 cells/mL was defined as mastitis, our established deep neural network was determined as the most suitable model for farm on-site mastitis detection, and this network model exhibited an accuracy of 75.93%, a specificity of 80.23%, a sensitivity of 70.35%, and AUC 0.83 on the test set. This study established a 1/4 level mastitis detection model which provides a theoretical basis for mastitis detection in buffaloes mostly raised by small farmers lacking mastitis diagnostic conditions in developing countries.

Biodiversity of Indonesian indigenous buffalo: First review of the status, challenges, and development opportunities

In Indonesia, the buffalo is important for small and marginal farmers' livelihood and economic development as a source of food, working animal, and tourist attraction. Therefore, an in-depth study is needed to examine challenges and opportunities for buffalo development in Indonesia. In Indonesia, the buffalo is divided into two types: swamp buffalo and river buffalo. The buffalo population in Indonesia has declined significantly. A decrease of approximately 39.35% was recorded from 2022 to 2017. The decline occurred due to low reproduction rate and suboptimal rearing management systems. There are three buffalo-rearing systems: Intensive, semi-intensive, and extensive. The productivity of buffalo is diverse and closely related to the characteristics of the regional agroecosystem, consistent with existing natural resources and rearing management systems. The diversity of buffalo productivity provides a good opportunity to improve productivity. Improvement of buffalo genetics is urgently needed, by improving mating management, etc., especially to reduce potential inbreeding. In recent years, genetic and molecular research on Indonesian buffalo has made progress, including use of molecular markers, such as microsatellites and single-nucleotide polymorphisms, to evaluate genetic diversity within and among buffalo populations across Indonesia. In addition, studies are being conducted on the relationship of genotype mutations that contribute to appearance and phenotypic performance (heat stress, reproduction, behavior, coat color, and production attributes) in buffaloes. Identification of genetic diversity in local buffaloes can be improved using various genetic and genomic techniques. These findings will form a basis for the targeted conservation of local buffaloes in Indonesia. This study aimed to collect information on the genetic resources of the local buffalo, particularly its status and production system and provide recommendations for developing buffalo production in Indonesia.

Effects of different load weights on the work performance and physiological and hematobiochemical responses in working water buffalo

Background and Aim

Working animals are important in agriculture because they play a role in various agricultural activities, including milk and meat production. Thus, they contribute to the development of rural communities. In this study, we aimed to evaluate the effects of different load weights on the work performance and the physiological and hematological responses of working water buffalo (Bubalus bubalis).

Materials and Methods

The work performances of 12 buffaloes (average weight 782.16 ± 21.62 kg) transporting 200, 350, and 500 kg of African palm fruits in metal baskets placed on their backs were evaluated. Work performance variables evaluated immediately after work were as follows: total number of trips (n), total weight (kg), distance traveled (km/day), working period (h), stopped time (h), and average speed (km/h). In addition, we evaluated physiological, biochemical, and hematological variables at three different times: before the start of work, immediately after work, and on the rest day.

Results

Among the load weights, 500 kg (total load carried = 4,138.88 kg) improved work efficiency compared to 200 kg loads (total load = 3,322.22 kg) (p = 0.0281). However, 500 kg loads resulted in slower average speed (2.4 km, p = 0.0164), shorter working period (2.39 h, p < 0.0001) and distance traveled (7.29 km, p < 0.0001), and less total number of trips (8.27 trips, p < 0.0001) compared to 350 and 200 kg load weights (3.45 and 3.52 km/h, 2.55 and 2.79 h, 8.71 and 9.75 km, 10.94 and 16.61 trips, respectively); and the heaviest loads resulted in significantly higher (p < 0.005) respiratory rate, pulse, heart rate, rectal temperature, glucose, lactate dehydrogenase, creatine kinase, total protein, white blood cell count, neutrophils, lymphocytes, monocytes, eosinophils, and basophils. In contrast, was associated with lower levels of red blood cells, hemoglobin, and hematocrit compared to lower loads. All differences were more pronounced in 500 kg (p < 0.005) compared to 200 kg loads.

Conclusion

Working buffaloes responded to work related to the transport of African palm fruits through various physiological, biochemical, and hematological changes. However, some variables remained close to the reference values reported in the literature for water buffaloes, and in general, all variables were reestablished during the rest day, indicating that these animals have adapted to working conditions.

Genome-wide identification and evolutionary analysis of the FGF gene family in buffalo

Fibroblast growth factors (FGFs) are important polypeptide growth factors that play a critical role in many developmental processes, including differentiation, cell proliferation, and migration in mammals. This study employs in silico analyses to characterize the FGF gene family in buffalo, investigating their genome-wide identification, physicochemical properties, and evolutionary patterns. For this purpose, genomic and proteomic sequences of buffalo, cattle, goat, and sheep were retrieved from NCBI database. We identified a total of 22 FGF genes in buffalo. Physicochemical properties observed through ProtParam tool showed notable features of these proteins including in-vitro instability, thermostability, hydrophilicity, and basic nature. Phylogenetic analysis grouped 22 identified genes into nine sub-families based on evolutionary relationships. Additionally, analysis of gene structure, motif patterns, and conserved domains using TBtools revealed the remarkable conservation of this gene family across selected species throughout the course of evolution. Comparative amino acid analysis performed through ClustalW demonstrated significant conservation between buffalo and cattle FGF proteins. Mutational analysis showed three non-synonymous mutations at positions R103 > G, P7 > L, and E98 > Q in FGF4, FGF6, and FGF19, respectively in buffalo. Duplication events revealed only one segmental duplication (FGF10/FGF22) in buffalo and two in cattle (FGF10/FGF22 and FGF13/FGF13-like) with Ka/Ks values <1 indicating purifying selection pressure for these duplications. Comparison of protein structures of buffalo, goat, and sheep exhibited more similarities in respective structures. In conclusion, our study highlights the conservation of the FGF gene family in buffalo during evolution. Furthermore, the identified non-synonymous mutations may have implications for the selection of animals with better performance.Communicated by Ramaswamy H. Sarma.

Genome-Wide Identification, Evolutionary and Mutational Analysis of the Buffalo Sox Gene Family

The Sox gene family constitutes transcription factors with a conserved high mobility group box (HMG) that regulate a variety of developmental processes, including sex differentiation, neural, cartilage, and early embryonic development. In this study, we systematically analyzed and characterized the 20 Sox genes from the whole buffalo genome, using comparative genomic and evolutionary analyses. All the buffalo Sox genes were divided into nine sub-groups, and each gene had a specific number of exons and introns, which contributed to different gene structures. Molecular phylogeny revealed more sequence similarity of buffalo Sox genes with those of cattle. Furthermore, evolutionary analysis revealed that the HMG domain remained conserved in the all members of the Sox gene family. Similarly, all the genes are under strong purifying selection pressure; seven segmental duplications occurred from 9.65 to 21.41 million years ago (MYA), and four potential recombination breakpoints were also predicted. Mutational analysis revealed twenty non-synonymous mutations with potential effects on physiological functions, including embryonic development and cell differentiation in the buffalo. The present study provides insights into the genetic architecture of the Sox gene family in buffalo, highlights the significance of mutations, and provides their potential utility for marker-assisted selection for targeted genetic improvement in buffalo.

Comparative evolutionary and molecular genetics based study of Buffalo lysozyme gene family to elucidate their antibacterial function

Lysozyme is used as a food preservative, biological medicine, and infant food additive as a natural anti-infective chemical having bactericidal activity and abundantly secreted in mammals' milk, saliva, etc. We systematically analyzed the 16 coding LYZ genes (C and G-type) in buffalo and cattle to elucidate their evolutionary perspective thoroughly by evaluating an evolutionary relationship, motif patterning, physicochemical attributes, gene, and protein structure, as well as the functional role of the mammary gland-specific expressed buffalo and cattle LYZ genes precisely while considering expression levels difference and the interaction sites variation with bacteria envisaged the potential ability of buffalo LYZ protein with enhanced antibacterial effect. Thus, we speculated that the buffalo mammary glands expressed lysozyme has good antibacterial activity. This study on the buffalo lysozyme gene family not only provides comprehensive insights into the genetic architecture and their antibacterial effect but also offers a theoretical basis for the development of new veterinary drugs and animal health care for mastitis, as well as a new molecular genetic basis to study food or medical lysozyme.

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.

The freezability of Mediterranean buffalo sperm is associated with lysine succinylation and lipid metabolism

Semen cryopreservation is used for the propagation of variety among species and domestic breeding. Mitochondria are implicated in sperm freezability, and their proteins are prone to succinylation, but the relationship between sperm freezability and mitochondrial protein succinylation is unclear. In this study, six bulls were classified as having good or poor freezability ejaculates (GFE or PFE, each 3 bulls). The fresh sperm mitochondrial membrane potential (MMP) and pan succinylation level of the two groups were first detected. Then the lysine succinylome and fatty acid content of the two groups were analyzed using label-free LC-MS/MS and GC-MS/MS in multiple reaction monitoring (MRM) modes, respectively. The results indicated that the GFE sperm had significantly higher MMPs than the PFE group (p < 0.05). A total of 1393 succinylation sites corresponding to 426 proteins were assessed and 5 succinylated peptides of the GFE group were markedly upregulated, while 3 were significantly downregulated (FC > 2.0 - < 0.5 and p-value < 0.05) when compared to the PFE group. Forty-six succinylated proteins were identified to have consistent presence/absence expression. The upregulated succinylated proteins in the GFE sperm were enriched in lipid metabolic processes. A total of 31 fatty acids were further subjected to quantitative analysis of which 23 including arachidic (C20:0), linolenic (C18:3n3), and docosahexaenoic acids (C22:6n3) were decreased in GFE sperm when compared with PFE (p < 0.05). These results suggest that lysine succinylation can potentially influence the sperm freezability of Mediterranean buffaloes through mitochondrial lipid metabolism. This novel study provides our understanding of sperm succinylation and the molecular basis for the mechanism of sperm freezability.

De Novo Transcriptome Dataset Generation of the Swamp Buffalo Brain and Non-Brain Tissues

The sequenced data availability opened new horizons related to buffalo genetic control of economic traits and genomic diversity. The visceral organs (brain, liver, etc.) significantly involved in energy metabolism, docility, or social interactions. We performed swamp buffalo transcriptomic profiling of 24 different tissues (brain and non-brain) to identify novel transcripts and analyzed the differentially expressed genes (DEGs) of brain vs. non-brain tissues with their functional annotation. We obtained 178.57 Gb clean transcriptomic data with GC contents 52.77%, reference genome alignment 95.36%, exonic coverage 88.49%. Totally, 26363 mRNAs transcripts including 5574 novel genes were obtained. Further, 7194 transcripts were detected as DEGs by comparing brain vs. non-brain tissues group, of which 3,999 were upregulated and 3,195 downregulated. These DEGs were functionally associated with cellular metabolic activities, signal transduction, cytoprotection, and structural and binding activities. The related functional pathways included cancer pathway, PI3k-Akt signaling, axon guidance, JAK-STAT signaling, basic cellular metabolism, thermogenesis, and oxidative phosphorylation. Our study provides an in-depth understanding of swamp buffalo transcriptomic data including DEGs potentially involved in basic cellular activities and development that helped to maintain their working capacity and social interaction with humans, and also, helpful to disclose the genetic architecture of different phenotypic traits and their gene expression regulation.

Comparative Genomic Characterization of Relaxin Peptide Family in Cattle and Buffalo

Relaxin family peptides significantly regulate reproduction, nutrient metabolism, and immune response in mammals. The present study aimed to identify and characterize the relaxin family peptides in cattle and buffalo at the genome level. The genomic and proteomic sequences of cattle, buffalo, goat, sheep, horse, and camel were accessed through the NCBI database, and BLAST was performed. We identified four relaxin peptides genes (RLN3, INSL3, INSL5, and INSL6) in Bos taurus, whereas three relaxin genes (RLN3, INSL3, and INSL6) in Bubalus bubalis. Evolutionary analysis revealed the conserved nature of relaxin family peptides in buffalo and cattle. Physicochemical properties revealed that relaxin proteins were thermostable, hydrophilic, and basic peptides except for INSL5 which was an acidic peptide. Three nonsynonymous mutations (two in RLN3 at positions A16 > T and P29 > A, and one in INSL6 at position R32 > Q) in Bos taurus, whereas two nonsynonymous mutations (one in RLN3 at positions G105 > w and one in INSL3 at position G22 > R) in Bubalus bubalis, were identified. INSL3 had one indel (insertion) at position 55 in Bos taurus. Gene duplication analysis revealed predominantly segmental duplications (INSL5/RLN3 and INSL6/INSL3 gene pairs) that helped expand this gene family, whereas Bubalus bubalis showed primarily tandem duplication (INSL3/RLN3). Our study concluded that relaxin family peptides remained conserved during the evolution, and gene duplications might help to adapt and enrich specific functions like reproduction, nutrient metabolism, and immune response. Further, the nonsynonymous mutations identified potentially affect these functions in buffalo.

Genomic analysis of arginine vasopressin gene in riverine buffalo reveals its potential association with silent estrus behavior

BACKGROUND

Poor estrus expression behavior causes suboptimal reproductive efficiency through poor conception rate. Various signaling pathways are involved in estrus expression but arginine vasopressin (AVP) gene with oxytocin predominantly regulates estrus behavior. This study aimed to perform genomic characterization and evolutionary dynamics of AVP gene through association testing of the novel polymorphic loci and comparative genomic analysis to explore the potential effect of AVP gene on estrus behavior of Nili-Ravi buffaloes.

METHODS AND RESULTS

198 Nili-Ravi buffaloes were screened for the quest of novel polymorphism in the AVP gene. In exon-1, five polymorphic sites were detected including deletion of two (c.47delA and c.57delA) nucleotides that caused drastic variation in subsequent amino acid sequence due to frame shift including functional short peptide of nine residues. The 3-D structure revealed a loss of transmembrane loop between 16 and 31 residues in Nili-Ravi buffalo AVP protein sequence, suggesting that missing loop apparently reduced the gene functionality in Nili-Ravi buffalo by inhibiting cellular reactions and muting the animal estrus cyclicity. Three polymorphisms detected in AVP gene were significantly associated with silent estrus (P < 0.05). The comparative genomic analysis revealed that AVP gene is present on chromosome 14 having one conserved motif (Neurohypophysial) in buffalo.

CONCLUSIONS

This study suggested the potential use of polymorphic sites as promising genetic markers for selection of buffaloes with pronounced estrus expression.

An Age Effect of Rumen Microbiome in Dairy Buffaloes Revealed by Metagenomics

Age is an important factor in shaping the gut microbiome. However, the age effect on the rumen microbial community for dairy buffaloes remains less explored. Using metagenomics, we examined the microbial composition and functions of rumen microbiota in dairy Murrah buffaloes of different ages: Y (1 year old), M (3−5 years old), E (6−8 years old), and O (>9 years old). We found that Bacteroidetes and Firmicutes were the predominant phyla, with Prevotella accounting for the highest abundance at the genus level. The proportion of Bacteroides and Methanobrevibacter significantly increased with age, while the abundance of genus Lactobacillus significantly decreased with age (LDA > 3, p < 0.05). Most differed COG and KEGG pathways were enriched in Y with carbohydrate metabolism, while older buffaloes enriched more functions of protein metabolism and the processing of replication and repair (LDA > 2, p < 0.05). Additionally, the functional contribution analysis revealed that the genera Prevotella and Lactobacillus of Y with more functions of CAZymes encoded genes of glycoside hydrolases and carbohydrate esterases for their roles of capable of metabolizing starch and sucrose-associated oligosaccharide enzyme, hemicellulase, and cellulase activities than the other three groups (LDA > 2, p < 0.05), thus affecting the 1-year-old dairy buffalo rumen carbohydrate metabolism. This study provides comprehensive dairy buffalo rumen metagenome data and assists in manipulating the rumen microbiome for improved dairy buffalo production.

Molecular Characterization of TGF-Beta Gene Family in Buffalo to Identify Gene Duplication and Functional Mutations

The TGF-β superfamily is ubiquitously distributed from invertebrates to vertebrates with diverse cellular functioning such as cell adhesion, motility, proliferation, apoptosis, and differentiation. The present study aimed to characterize the TGF-β gene superfamily in buffalo through evolutionary, structural, and single nucleotide polymorphism (SNPs) analyses to find the functional effect of SNPs in selected genes. We detected 32 TGF-β genes in buffalo genome and all TGF-β proteins exhibited basic nature except INHA, INHBC, MSTN, BMP10, and GDF2, which showed acidic properties. According to aliphatic index, TGF-β proteins were thermostable but unstable in nature. Except for GDF1 and AMH, TGF-β proteins depicted hydrophilic nature. Moreover, all the detected buffalo TGF-β genes showed evolutionary conserved nature. We also identified eight segmental and one tandem duplication event TGF-β gene family in buffalo, and the ratio of Ka/Ks demonstrated that all the duplicated gene pairs were under selective pressure. Comparative amino acid analysis demonstrated higher variation in buffalo TGF-β gene family, as a total of 160 amino acid variations in all the buffalo TGF-β proteins were detected. Mutation analysis revealed that 13 mutations had an overall damaging effect that might have functional consequences on buffalo growth, folliculogenesis, or embryogenesis.

Generation of Heritable Prominent Double Muscle Buttock Rabbits via Novel Site Editing of Myostatin Gene Using CRISPR/Cas9 System

Rabbits have been domesticated for meat, wool, and fur production, and have also been cherished as a companion, artistic inspiration, and an experimental model to study many human diseases. In the present study, the muscle mass negative regulator gene myostatin (MSTN) was knocked out in rabbits at two novel sites in exon3, and the function of these mutations was determined in subsequent generations. The prominent double muscle phenotype with hyperplasia or hypertrophy of muscle fiber was observed in the MSTN-KO rabbits, and a similar phenotype was confirmed in the F1 generation. Moreover, the average weight of 80-day-old MSTN-KO rabbits (2,452 ± 63 g) was higher than that of wild-type rabbits (2,393.2 ± 106.88 g), and also the bodyweight of MSTN-KO rabbits (3,708 ± 43.06g) was significantly higher (P < 0.001) at the age of 180 days than wild-type (WT) rabbits (3,224 ± 48.64g). In MSTN-KO rabbits, fourteen rabbit pups from the F1 generation and thirteen from the F2 generation stably inherited the induced MSTN gene mutations. Totally, 194 pups were produced in the F1 generation of which 49 were MSTN-KO rabbits, while 47 pups were produced in the F2 generation of which 20 were edited rabbits, and the ratio of edited to wild-type rabbits in the F2 generation was approximately 1:1. Thus, we successfully generated a heritable double muscle buttocks rabbits via myostatin mutation with CRISPR/Cas9 system, which could be valuable in rabbit's meat production and also a useful animal model to study the development of muscles among livestock species and improve their important economic traits as well as the human muscle development-related diseases.

Genomic Identification, Evolution, and Expression Analysis of Bromodomain Genes Family in Buffalo

Bromodomain (BRD) is an evolutionarily conserved protein-protein interaction module that is critical in gene regulation, cellular homeostasis, and epigenetics. This study aimed to conduct an identification, evolution, and expression analysis of the BRD gene family in the swamp buffalo (Bubalus bubalis). A total of 101 BRD protein sequences deduced from 22 BRD genes were found in the buffalo genome. The BRD proteins were classified into six groups based on phylogenetic relationships, conserved motifs, and conserved domains. The BRD genes were irregularly distributed in 13 chromosomes. Collinearity analysis revealed 20 BRD gene pairs that had remarkable homologous relationships between the buffalo and cattle, although no tandem or segmental duplication event was found in the buffalo BRD genes. Comparative transcriptomics using a 10x sequencing platform analysis showed that 22 BRD genes were identified in the Sertoli cells (SCs) at different developmental stages of buffalo. Further, the mRNA expression levels of bromodomain and the extraterminal (BET) family in SCs at the pubertal stage were higher than that at the prepubertal stage of buffalo. However, the SMARCA2, PHIP, BRD9, and TAF1 genes exhibited the opposite trend. The maturation process of SCs may be regulated by the BRD family members expressed differentially in SCs at different developmental stages of buffalo. In summary, our findings provide an understanding of the evolutionary, structural, and functional properties of the buffalo BRD family members, and further characterize the function of the BRD family in the maturation of SCs. It also provides a theoretical basis for further understanding in the future of the mechanism of SCs regulating spermatogenesis.

Comparative Genomic Characterization of Buffalo Fibronectin Type III Domain Proteins: Exploring the Novel Role of FNDC5/Irisin as a Ligand of Gonadal Receptors

FN-III proteins are widely distributed in mammals and are usually involved in cellular growth, differentiation, and adhesion. The FNDC5/irisin regulates energy metabolism and is present in different tissues (liver, brain, etc.). The present study aimed to investigate the physiochemical characteristics and the evolution of FN-III proteins and FNDC5/irisin as a ligand targeting the gonadal receptors including androgen (AR), DDB1 and CUL4 associated factor 6 (DCAF6), estrogen-related receptor β (ERR-β), estrogen-related receptor γ (ERR-γ), Krüppel-like factor 15 (KLF15), and nuclear receptor subfamily 3 group C member 1 (NR3C1). Moreover, the putative role of irisin in folliculogenesis and spermatogenesis was also elucidated. We presented the molecular structure and function of 29 FN-III genes widely distributed in the buffalo genome. Phylogenetic analysis, motif, and conserved domain pattern demonstrated the evolutionary well-conserved nature of FN-III proteins with a variety of stable to unstable, hydrophobic to hydrophilic, and thermostable to thermo-unstable properties. The comparative structural configuration of FNDC5 revealed amino acid variations but still the FNDC5 structure of humans, buffalo, and cattle was quite similar to each other. For the first time, we predicted the binding scores and interface residues of FNDC5/irisin as a ligand for six representative receptors having a functional role in energy homeostasis, and a significant involvement in folliculogenesis and spermatogenesis in buffalo.