Genomic Identification, Evolution, and Expression Analysis of Collagen Genes Family in Water Buffalo during Lactation

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 copy number variant analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis)

Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variants (CNV) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 CNV regions (CNVR), with 1,993 shared CNVR being found within the studied buffalo types. Analyzing CNVR highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVR that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci analysis revealed differentially expressed CNVR-derived genes (DECG) associated with milk production traits. Notably, known milk production-related genes were among these DECG, validating their relevance. Last, a GWAS identified 3 CNVR significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.

Functional transcriptome analysis revealed upregulation of MAPK-SMAD signalling pathways in chronic heat stress in crossbred cattle

Animal geneticists and breeders have the impending challenge of enhancing the resilience of Indian livestock to heat stress through better selection strategies. Climate change's impact on livestock is more intense in tropical countries like India where dairy cattle crossbreeds are more sensitive to heat stress. The main reason for this study was to find the missing relative changes in transcript levels in thermo-neutral and heat stress conditions in crossbred cattle through whole-transcriptome analysis of RNA-Seq data. Differentially expressed genes (DEGs) identified based on the minimum log twofold change value and false discovery rate 0.05 revealed 468 up-regulated genes and 2273 down-regulated significant genes. Functional annotation and pathway analysis of these significant DEGs were compared based on Gene Ontology (Biological process), Kyoto Encyclopedia of Genes and Genome (KEGG), and Reactome pathways using g: Profiler, ShinyGO v0.76, and iDEP.951 web tools. On finding network visualization, the most over-represented and correlated pathways were neuronal and sensory organ development, calcium signalling pathway, Mitogen-activated protein kinase (MAPK) and Smad signalling pathway, Ras-proximate-1, or Ras-related protein 1 (Rap 1) signalling pathway, apoptosis, and oxidative stress. Similarly, down-regulated genes were most expressed in mRNA processing, immune system, B-cell receptor signalling pathway, Nucleotide oligomerization domain (NOD)-like receptors (NLRs) signalling pathway and nonsense-mediated decay (NMD) pathway. The heat stress-responsive genes identified in this study will facilitate our understanding of the molecular basis for climate resilience and heat tolerance in Indian dairy crossbreeds.

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.

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.

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.

Comparative Genomic Analysis of the Thiolase Family and Functional Characterization of the Acetyl-Coenzyme A Acyltransferase-1 Gene for Milk Biosynthesis and Production of Buffalo and Cattle

Cattle and buffalo served as the first and second largest dairy animals, respectively, providing 96% milk products worldwide. Understanding the mechanisms underlying milk synthesis is critical to develop the technique to improve milk production. Thiolases, also known as acetyl-coenzyme A acetyltransferases (ACAT), are an enzyme family that plays vital roles in lipid metabolism, including ACAT1, ACAT2, ACAA1, ACAA2, and HADHB. Our present study showed that these five members were orthologous in six livestock species including buffalo and cattle. Transcriptomic data analyses derived from different lactations stages showed that ACAA1 displayed different expression patterns between buffalo and cattle. Immunohistochemistry staining revealed that ACAA1 were dominantly located in the mammary epithelial cells of these two dairy animals. Knockdown of ACAA1 inhibited mammary epithelial cell proliferation and triglyceride and β-casein secretion by regulating related gene expressions in cattle and buffalo. In contrast, ACAA1 overexpression promoted cell proliferation and triglyceride secretion. Finally, three novel SNPs (g.-681A>T, g.-23117C>T, and g.-24348G>T) were detected and showed significant association with milk production traits of Mediterranean buffaloes. In addition, g.-681A>T mutation located in the promoter region changed transcriptional activity significantly. Our findings suggested that ACAA1 play a key role in regulating buffalo and cattle milk synthesis and provided basic information to further understand the dairy animal lactation physiology.

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.

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.

Integrative Systems Biology Analysis Elucidates Mastitis Disease Underlying Functional Modules in Dairy Cattle

Background: Mastitis is the most prevalent disease in dairy cattle and one of the most significant bovine pathologies affecting milk production, animal health, and reproduction. In addition, mastitis is the most common, expensive, and contagious infection in the dairy industry. Methods: A meta-analysis of microarray and RNA-seq data was conducted to identify candidate genes and functional modules associated with mastitis disease. The results were then applied to systems biology analysis via weighted gene coexpression network analysis (WGCNA), Gene Ontology, enrichment analysis for the Kyoto Encyclopedia of Genes and Genomes (KEGG), and modeling using machine-learning algorithms. Results: Microarray and RNA-seq datasets were generated for 2,089 and 2,794 meta-genes, respectively. Between microarray and RNA-seq datasets, a total of 360 meta-genes were found that were significantly enriched as "peroxisome," "NOD-like receptor signaling pathway," "IL-17 signaling pathway," and "TNF signaling pathway" KEGG pathways. The turquoise module (n = 214 genes) and the brown module (n = 57 genes) were identified as critical functional modules associated with mastitis through WGCNA. PRDX5, RAB5C, ACTN4, SLC25A16, MAPK6, CD53, NCKAP1L, ARHGEF2, COL9A1, and PTPRC genes were detected as hub genes in identified functional modules. Finally, using attribute weighting and machine-learning methods, hub genes that are sufficiently informative in Escherichia coli mastitis were used to optimize predictive models. The constructed model proposed the optimal approach for the meta-genes and validated several high-ranked genes as biomarkers for E. coli mastitis using the decision tree (DT) method. Conclusion: The candidate genes and pathways proposed in this study may shed new light on the underlying molecular mechanisms of mastitis disease and suggest new approaches for diagnosing and treating E. coli mastitis in dairy cattle.

Genome-Wide Association Studies for Milk Somatic Cell Score in Romanian Dairy Cattle

Mastitis is one of the most frequently encountered diseases in dairy cattle, negatively affecting animal welfare and milk production. For this reason, contributions to understanding its genomic architecture are of great interest. Genome-wide association studies (GWAS) have identified multiple loci associated with somatic cell score (SCS) and mastitis in cattle. However, most of the studies have been conducted in different parts of the world on various breeds, and none of the investigations have studied the genetic architecture of mastitis in Romanian dairy cattle breeds up to this point in time. In this study, we report the first GWAS for SCS in dairy cattle breeds from Romania. For GWAS, we used an Axiom Bovine v3 SNP-chip (>63,000 Single Nucleotide Polymorphism -SNPs) and 33,330 records from 690 cows belonging to Romanian Spotted (RS) and Romanian Brown (RB) cattle. The results found one SNP significantly associated with SCS in the RS breed and 40 suggestive SNPs with -log₁₀ (p) from 4 to 4.9 for RS and from 4 to 5.4 in RB. From these, 14 markers were located near 12 known genes (AKAP8, CLHC1, MEGF10, SATB2, GATA6, SPATA6, COL12A1, EPS8, LUZP2, RAMAC, IL12A and ANKRD55) in RB cattle, 3 markers were close to ZDHHC19, DAPK1 and MMP7 genes, while one SNP overlapped the HERC3 gene in RS cattle. Four genes (HERC3, LUZP2, AKAP8 and MEGF10) associated with SCS in this study were previously reported in different studies. The most significant SNP (rs110749552) associated with SCS was located within the HERC3 gene. In both breeds, the SNPs and position of association signals were distinct among the three parities, denoting that mastitis is controlled by different genes that are dependent according to parity. The current results contribute to an expansion in the body of knowledge regarding the proportion of genetic variability explained by SNPs for SCS in dairy cattle.

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

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