Domesticated buffaloes have been integral to rice-paddy agro-ecosystems for millennia, yet relatively little is known about the buffalo genomics. Here, we sequenced and assembled reference genomes for both swamp and river buffaloes and we re-sequenced 230 individuals (132 swamp buffaloes and 98 river buffaloes) sampled from across Asia and Europe. Beyond the many actionable insights that our study revealed about the domestication, basic physiology and breeding of buffalo, we made the striking discovery that the divergent domestication traits between swamp and river buffaloes can be explained with recent selections of genes on social behavior, digestion metabolism, strengths and milk production.

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.

Beef is considered to be a good quality meat product because it contains linoleic acid and specific proteins, which can bring significant benefits to health. Circular RNAs (circRNAs) have been reported to regulate skeletal myogenesis. RNA-seq was used to investigate the circRNA molecular regulatory mechanisms with respect to differences in muscle quality between buffalo and cattle. A total of 10,449 circRNA candidates were detected, and 1128 of these were found to be differentially expressed between cattle and buffalo muscle tissue libraries. Differentially expressed 23 circRNAs were verified by qPCR. CircEch1, one of the most up-regulated circRNAs during muscle development, was subsequently characterized. CCK-8 (65.05 ± 2.33%, P < 0.0001), EdU (72.99 ± 0.04%, P < 0.001), and Western blotting assays showed that overexpression of circEch1 inhibited the proliferation of bovine myoblasts but promoted differentiation. In vivo studies suggested that circEch1 stimulates skeletal muscle regeneration. These results demonstrate that the novel regulator circEch1 induces myoblast differentiation and skeletal muscle regeneration. They also provide new insights into the mechanisms of circRNA regulation of beef quality.

The buffalo was domesticated around 3000–6000 years ago and has substantial economic significance as a meat, dairy, and draught animal. The buffalo has remained underutilized in terms of the development of a well-annotated and assembled reference genome de novo. It is mandatory to explore the genetic architecture of a species to understand the biology that helps to manage its genetic variability, which is ultimately used for selective breeding and genomic selection. Morphological and molecular data have revealed that the swamp buffalo population has strong geographical genomic diversity with low gene flow but strong phenotypic consistency, while the river buffalo population has higher phenotypic diversity with a weak phylogeographic structure. The availability of recent high-quality reference genome and genotyping marker panels has invigorated many genome-based studies on evolutionary history, genetic diversity, functional elements, and performance traits. The increasing molecular knowledge syndicate with selective breeding should pave the way for genetic improvement in the climatic resilience, disease resistance, and production performance of water buffalo populations globally.

Fasciola gigantica and Fasciola hepatica are causative pathogens of fascioliasis, with the widest latitudinal, longitudinal, and altitudinal distribution; however, among parasites, they have the largest sequenced genomes, hindering genomic research. In the present study, we used various sequencing and assembly technologies to generate a new high-quality Fasciola gigantica reference genome. We improved the integration of gene structure prediction, and identified two independent transposable element expansion events contributing to (1) the speciation between Fasciola and Fasciolopsis during the Cretaceous-Paleogene boundary mass extinction, and (2) the habitat switch to the liver during the Paleocene-Eocene Thermal Maximum, accompanied by gene length increment. Long interspersed element (LINE) duplication contributed to the second transposon-mediated alteration, showing an obvious trend of insertion into gene regions, regardless of strong purifying effect. Gene ontology analysis of genes with long LINE insertions identified membrane-associated and vesicle secretion process proteins, further implicating the functional alteration of the gene network. We identified 852 predicted excretory/secretory proteins and 3300 protein-protein interactions between Fasciola gigantica and its host. Among them, copper/zinc superoxide dismutase genes, with specific gene copy number variations, might play a central role in the phase I detoxification process. Analysis of 559 single-copy orthologs suggested that Fasciola gigantica and Fasciola hepatica diverged at 11.8 Ma near the Middle and Late Miocene Epoch boundary. We identified 98 rapidly evolving gene families, including actin and aquaporin, which might explain the large body size and the parasitic adaptive character resulting in these liver flukes becoming epidemic in tropical and subtropical regions.

Buffalo is a luxurious genetic resource with multiple utilities (as a dairy, draft, and meat animal) and economic significance in the tropical and subtropical regions of the globe. The excellent potential to survive and perform on marginal resources makes buffalo an important source for nutritious products, particularly milk and meat. This study was aimed to investigate the evolutionary relationship, physiochemical properties, and comparative genomic analysis of the casein gene family (CSN1S1, CSN2, CSN1S2, and CSN3) in river and swamp buffalo. Phylogenetic, gene structure, motif, and conserved domain analysis revealed the evolutionarily conserved nature of the casein genes in buffalo and other closely related species. Results indicated that casein proteins were unstable, hydrophilic, and thermostable, although αs1-CN, β-CN, and κ-CN exhibited acidic properties except for αs2-CN, which behaved slightly basic. Comparative analysis of amino acid sequences revealed greater variation in the river buffalo breeds than the swamp buffalo indicating the possible role of these variations in the regulation of milk traits in buffalo. Furthermore, we identified lower transcription activators STATs and higher repressor site YY1 distribution in swamp buffalo, revealing its association with lower expression of casein genes that might subsequently affect milk production. The role of the main motifs in controlling the expression of casein genes necessitates the need for functional studies to evaluate the effect of these elements on the regulation of casein gene function in buffalo.

Pangolins are one of nature’s most fascinating species being scales covered and myrmecophagous diet, yet relatively little is known about the molecular basis. Here, we combine the multi-omics, evolution, and fundamental proteins feature analysis of both Chinese and Malayan pangolins, highlighting the molecular mechanism of both myrmecophagous diet and scale formation, representing a fascinating evolutionary strategy to occupy the unique ecological niches. In contrast to conserved organization of epidermal differentiation complex (EDC), pangolin has undergone large scale variation and gene loss events causing expression pattern and function conversion that contributes to cornified epithelium structures on stomach to adapt myrmecophagous diet. Our assemblies also enable us to discover large copies number of high glycine-tyrosine keratin associated proteins. In addition, highly homogenized tandem array, amino content, and the specific expression pattern further validate the strong connection between the molecular mechanism of scale hardness and HGT-KRTAPs.

Fasciola gigantica and Fasciola hepatica are causative pathogens of fascioliasis , with the widest latitudinal, longitudinal, and altitudinal distribution; however, among parasites, they have the largest sequenced genomes, hindering genomic research. In the present study, we used various sequencing and assembly technologies to generate a new high-quality Fasciola gigantica reference genome. We improved the integration of gene structure prediction, and identified two independent transposable element expansion events contributing to (1) the speciation between Fasciola and Fasciolopsis during the Cretaceous-Paleogene boundary mass extinction, and (2) the habitat switch to the liver during the Paleocene-Eocene Thermal Maximum , accompanied by gene length increment. Long interspersed element (LINE) duplication contributed to the second transposon-mediated alteration, showing an obvious trend of insertion into gene regions, regardless of strong purifying selection. Gene ontology analysis of genes with long LINE insertions identified membrane-associated and vesicle secretion process proteins, further implicating the functional alteration of the gene network . We identified 852 excretory/secretory proteins and 3300 protein-protein interactions between Fasciola gigantica and its host. Among them, copper/zinc superoxide dismutase genes, with specific gene copy number variations, might play a central role in the phase I detoxification process. Analysis of 559 single-copy orthologs suggested that Fasciola gigantica and Fasciola hepatica diverged at 11.8 Ma near the Middle and Late Miocene Epoch boundary. We identified 98 rapidly evolving gene families, including actin and aquaporin, which might explain the large body size and the parasitic adaptive character resulting in these liver flukes becoming epidemic in tropical and subtropical regions.