A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds

Conjugated linoleic acid (CLA) modulates bovine peripheral blood mononuclear cells (PBMC) proteome in vitro

Conjugated linoleic acid (CLA) is a group of natural isomers of the n-6 polyunsaturated fatty acid (PUFA) linoleic acid, exerting biological effects on cow physiology. This study assessed the impact of the mixture 50:50 (vol:vol) of CLA isomers (cis-9, trans-11 and trans-10, cis-12) on bovine peripheral blood mononuclear cells (PBMC) proteome, identifying 1608 quantifiable proteins. A supervised multivariate statistical analysis, sparse variant partial least squares - discriminant analysis (sPLS-DA) for paired data identified 407 discriminant proteins (DP), allowing the clustering between the CLA and controls. The ProteINSIDE workflow found that DP with higher abundance in the CLA group included proteins related to innate immune defenses (PLIN2, CD36, C3, C4, and AGP), with antiapoptotic (SERPINF2 and ITIH4) and antioxidant effects (HMOX1). These results demonstrated that CLA modulates the bovine PBMC proteome, supports the antiapoptotic and immunomodulatory effects observed in previous in vitro studies on bovine PBMC, and suggests a cytoprotective role against oxidative stress. SIGNIFICANCE: In this study, we report for the first time that the mixture 50:50 (vol:vol) of cis-9, trans-11, and trans-10, cis-12-CLA isomers modulates the bovine PBMC proteome. Our results support the immunomodulatory and antiapoptotic effects observed in bovine PBMC in vitro. In addition, the present study proposes a cytoprotective role of CLA mixture against oxidative stress. We suggest a molecular signature of CLA treatment based on combining a multivariate sparse discriminant analysis and a clustering method. This demonstrates the great value of sPLS-DA as an alternative option to identify discriminant proteins with relevant biological significance.

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.

SWATH-MS Analysis of Blood Plasma and Circulating Small Extracellular Vesicles Enables Detection of Putative Protein Biomarkers of Fertility in Young and Aged Dairy Cows

The development of biomarkers of fertility could provide benefits for the genetic improvement of dairy cows. Circulating small extracellular vesicles (sEVs) show promise as diagnostic or prognostic markers since their cargo reflects the metabolic state of the cell of origin; thus, they mirror the physiological status of the host. Here, we employed data-independent acquisition mass spectrometry to survey the plasma and plasma sEV proteomes of two different cohorts of Young (Peripubertal; n = 30) and Aged (Primiparous; n = 20) dairy cows (Bos taurus) of high- and low-genetic merit of fertility and known pregnancy outcomes (ProteomeXchange data set identifier PXD042891). We established predictive models of fertility status with an area under the curve of 0.97 (sEV; p value = 3.302e-07) and 0.95 (plasma; p value = 6.405e-08). Biomarker candidates unique to high-fertility Young cattle had a sensitivity of 0.77 and specificity of 0.67 (*p = 0.0287). Low-fertility biomarker candidates uniquely identified in sEVs from Young and Aged cattle had a sensitivity and specificity of 0.69 and 1.0, respectively (***p = 0.0005). Our bioinformatics pipeline enabled quantification of plasma and circulating sEV proteins associated with fertility phenotype. Further investigations are warranted to validate this research in a larger population, which may lead to improved classification of fertility status in cattle.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

A Study of Heavy Metal Status and its Relationship with Hematologic and Biochemical Indices in River Buffaloes in Southwest Iran

Heavy metals are among the most important environmental pollutants which accumulate in various organs and are associated with several toxic effects. This study was performed to determine the status of heavy metals in river buffaloes in Khuzestan province, Iran, and its relationship with hematologic and serum biochemical parameters. A total of 103 apparently healthy buffaloes were sampled from the region. The concentration of heavy metals, including lead (Pb), mercury (Hg), and cadmium (Cd), was determined in serum samples by atomic spectroscopy. In addition, complete blood counts and serum biochemical profiles were assessed. The serum concentration of Cd, Pb, and Hg in the sampled buffaloes, as mean±standard error, were 0.55±0.01, 6.51±0.10, and 6.28±0.09 µg/l, respectively, which are within the permissible serum levels in the livestock. Serum Cd and Hg levels showed no significant relationship with hematologic or biochemical analytes. However, there were significant negative correlations between Pb levels and phosphorus, magnesium, sodium, as well as potassium concentrations, while serum iron was positively correlated with lead (P<0.05). In addition, there was a significant positive correlation between Hg level and serum aspartate aminotransferase activity (P<0.05). Despite the fact that river buffaloes in Khuzestan spend a long time daily in the Karun River with high industrial pollution, no serum evidence of heavy metal toxicity was found in these animals. It can be suggested that river buffaloes in Khuzestan seem to be resistant to the environmental pollution caused by heavy metals. However, further studies are required to confirm this issue and identify its possible explanations.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Linkage Disequilibrium and Effective Population Size of Buffalo Populations of Iran, Turkey, Pakistan, and Egypt Using a Medium Density SNP Array

Linkage disequilibrium (LD) across the genome provides information to identify the genes and variations related to quantitative traits in genome-wide association studies (GWAS) and for the implementation of genomic selection (GS). LD can also be used to evaluate genetic diversity and population structure and reveal genomic regions affected by selection. LD structure and Ne were assessed in a set of 83 water buffaloes, comprising Azeri (AZI), Khuzestani (KHU), and Mazandarani (MAZ) breeds from Iran, Kundi (KUN) and Nili-Ravi (NIL) from Pakistan, Anatolian (ANA) buffalo from Turkey, and buffalo from Egypt (EGY). The values of corrected r² (defined as the correlation between two loci) of adjacent SNPs for three pooled Iranian breeds (IRI), ANA, EGY, and two pooled Pakistani breeds (PAK) populations were 0.24, 0.28, 0.27, and 0.22, respectively. The corrected r² between SNPs decreased with increasing physical distance from 100 Kb to 1 Mb. The LD values for IRI, ANA, EGY, and PAK populations were 0.16, 0.23, 0.24, and 0.21 for less than 100Kb, respectively, which reduced rapidly to 0.018, 0.042, 0.059, and 0.024, for a distance of 1 Mb. In all the populations, the decay rate was low for distances greater than 2Mb, up to the longest studied distance (15 Mb). The r² values for adjacent SNPs in unrelated samples indicated that the Affymetrix Axiom 90 K SNP genomic array was suitable for GWAS and GS in these populations. The persistency of LD phase (PLDP) between populations was assessed, and results showed that PLPD values between the populations were more than 0.9 for distances of less than 100 Kb. The Ne in the recent generations has declined to the extent that breeding plans are urgently required to ensure that these buffalo populations are not at risk of being lost. We found that results are affected by sample size, which could be partially corrected for; however, additional data should be obtained to be confident of the results.

Signatures of selection reveal candidate genes involved in production traits in Chinese crossbred buffaloes

Identification of selection signature is important for a better understanding of genetic mechanisms that affect phenotypic differentiation in livestock. However, the genome-wide selection responses have not been investigated for the production traits of Chinese crossbred buffaloes. In this study, an SNP data set of 133 buffaloes (Chinese crossbred buffalo, n = 45; Chinese local swamp buffalo, n = 88) was collected from the Dryad Digital Repository database (https://datadryad.org/stash/). Population genetics analysis showed that these buffaloes were divided into the following 2 groups: crossbred buffalo and swamp buffalo. The crossbred group had higher genetic diversity than the swamp group. Using 3 complementary statistical methods (integrated haplotype score, cross population extended haplotype homozygosity, and composite likelihood ratio), a total of 31 candidate selection regions were identified in the Chinese crossbred population. Here, within these candidate regions, 25 genes were under the putative selection. Among them, several candidate genes were reported to be associated with production traits. In addition, we identified 13 selection regions that overlapped with bovine QTLs that were mainly involved in milk production and composition traits. These results can provide useful insights regarding the selection response for production traits of Chinese crossbred buffalo, as identified candidate genes influence production performance.

Peroxisome proliferator-activated receptor gamma regulates genes involved in milk fat synthesis in mammary epithelial cells of water buffalo

Peroxisome proliferator-activated receptor gamma (PPARγ) is a critical transcription factor regulating lipid and glucose metabolism. However, the regulatory effect of PPARγ on milk fat synthesis in buffalo mammary gland is not clear. In order to explore the role of buffalo PPARG gene in milk fat synthesis, lentivirus-mediated interference was used to knock it down and then the recombinant fusion expression vector was transfected into buffalo mammary epithelial cell (BMEC) to overexpress it. PPARG gene knockdown significantly decreased the expression of CD36, FABP3, FABP4, ACSS2, ELOVL6, DGAT2, BTN1A1, AGPAT6, LPIN1, ABCG2, PPARGC1A, INSIG1, FASN, and SREBF2 genes and significantly upregulated the expression of INSIG2 gene but had no significant effect on the expression of ACSL1, GPAM, and SREBF1 genes. PPARG overexpression significantly increased the relative mRNA abundance of CD36, FABP3, FABP4, ACSS2, ELOVL6, DGAT2, BTN1A1, AGPAT6, LPIN1, PPARGC1A, INSIG1, and SREBF2 genes and significantly downregulated the expression of INSIG2 gene but had no significant effect on the expression of ACSL1, GPAM, ABCG2, FASN, and SREBF1 genes. In addition, knockdown/overexpression of PPARG gene significantly decreased/increased triacylglycerol (TAG) content in BMECs. This study revealed that buffalo PPARG gene is a key gene regulating buffalo milk fat synthesis.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events

Background

Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types.

Results

The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction.

Conclusions

Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00616-3.

Functional and Flow Cytometric Analysis of Buffalo Cryopreserved Spermatozoa: Comparison of Different Breeds and Incubation Times

BACKGROUND

The purpose of this research was to compare the functional parameters of frozen-thawed Iranian Azari buffalo spermatozoa with imported semen samples of Italian Mediterranean buffalo (IMB) after the thawing process and 4 hours of incubation.

MATERIALS AND METHODS

In this experimental study, a total of twenty-four ejaculates from four Iranian Azari buffalo bulls were collected. Semen samples were diluted in AndroMed extender at a concentration of 50×10⁶ spermatozoa/ ml. The diluted samples were filled in 0.5 ml straws and were frozen in a programmable freezer. For imported semen samples, twenty-four samples of four IMB were used, which were diluted in AndroMed extender and frozen by the same procedure. Frozen-thawed sperm motion patterns, mitochondrial activity, membrane integrity, DNA integrity, reactive oxygen species (ROS), and apoptosis status were evaluated immediately after thawing and 4 hours of incubation.

RESULTS

Post-thawed sperm motility, progressive motility (PM), mitochondrial activity, membrane integrity were significantly higher in imported semen samples in compare with Iranian Azari buffalo. After 4 hours of incubation, sperm velocity patterns were superior in Iranian Azari semen samples. Moreover, the percentage of sperm cells with un-damaged DNA was higher in Iranian semen samples compared to imported samples at the time 0 of incubation. Following 4 hours of incubation, a significant increase in intracellular ROS level leads to reduced membrane integrity, mitochondrial activity, and DNA integrity in both buffalo breeds. At time 4, Iranian samples showed significantly lower apoptosis and higher dead spermatozoa compared to imported semen samples.

CONCLUSION

Our study showed that the post-thawed quality of Iranian Azari buffalo semen was comparable with imported samples after 4 hours of incubation. Further investigations are recommended to assess the in vitro and in vivo fertility rate of both buffalo breeds.

Bubalus bubalis: A Short Story

The domestic buffalo (Bubalus bubalis), also known as water buffalo or Asian buffalo to prevent confusion with the American bison (Bison bison), wrongly named buffalo in North America, comprises two subspecies: the river buffalo (B. bubalis bubalis) and the swamp buffalo (B. bubalis kerebau). The swamp buffalo has a consistent phenotype and is considered as one type, even if many breeds are recognized within it; conversely, the river buffalo subspecies has many breeds. We found limited information available regarding the worldwide distribution of buffaloes. The best estimate is that 208,098,759 buffalo head are distributed in 77 countries in five continents. In this review, we presented the basic aspects of the water buffalo and unraveled the buffalo path followed from the origin of the species to its current global distribution. We reviewed several data sources to provide a better estimate of the world buffalo count and distribution.

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.

Whole-Genome Resequencing Reveals Adaptation Prior to the Divergence of Buffalo Subspecies

The application of high-throughput genotyping or sequencing data helps us to understand the genomic response to natural and artificial selection. In this study, we scanned the genomes of five indigenous buffalo populations belong to three recognized breeds, adapted to different geographical and agro-ecological zones in Iran, to unravel the extent of genomic diversity and to localize genomic regions and genes underwent past selection. A total of 46 river buffalo whole genomes, from West and East Azerbaijan, Gilan, Mazandaran, and Khuzestan provinces, were resequenced. Our sequencing data reached to a coverage above 99% of the river buffalo reference genome and an average read depth around 9.2× per sample. We identified 20.55 million SNPs, including 63,097 missense, 707 stop-gain, and 159 stop-loss mutations that might have functional consequences. Genomic diversity analyses showed modest structuring among Iranian buffalo populations following frequent gene flow or admixture in the recent past. Evidence of positive selection was investigated using both differentiation (Fst) and fixation (Pi) metrics. Analysis of fixation revealed three genomic regions in all three breeds with aberrant polymorphism contents on BBU2, 20, and 21. Fixation signal on BBU2 overlapped with the OCA2-HERC2 genes, suggestive of adaptation to UV exposure through pigmentation mechanism. Further validation using resequencing data from other five bovine species as well as the Axiom Buffalo Genotyping Array 90K data of river and swamp buffaloes indicated that these fixation signals persisted across river and swamp buffaloes and extended to taurine cattle, implying an ancient evolutionary event occurred before the speciation of buffalo and taurine cattle. These results contributed to our understanding of major genetic switches that took place during the evolution of modern buffaloes.

Genome-Wide Detection of Selection Signatures in Duroc Revealed Candidate Genes Relating to Growth and Meat Quality

With the development of high-throughput genotyping techniques, selection signatures in the genome of domestic pigs have been extensively interrogated in the last decade. The Duroc, a major commercial pig breed famous for its fast growth rate and high lean ratio, has not been extensively studied focusing on footprints of intensively artificial selection in their genomes by a lot of re-sequencing data. The goal of this study was to investigate genomic regions under artificial selection and their contribution to the unique phenotypic traits of the Duroc using whole-genome resequencing data from 97 pigs. Three complementary methods (d_i, CLR, and iHH12) were implemented for selection signature detection. In Total, 464 significant candidate regions were identified, which covered 46.4 Mb of the pig genome. Within the identified regions, 709 genes were annotated, including 600 candidate protein-coding genes (486 functionally annotated genes) and 109 lncRNA genes. Genes undergoing selective pressure were significantly enriched in the insulin resistance signaling pathway, which may partly explain the difference between the Duroc and other breeds in terms of growth rate. The selection signatures identified in the Duroc population demonstrated positive pressures on a set of important genes with potential functions that are involved in many biological processes. The results provide new insights into the genetic mechanisms of fast growth rate and high lean mass, and further facilitate follow-up studies on functional genes that contribute to the Duroc's excellent phenotypic traits.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences

BACKGROUND

Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal.

RESULTS

In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits.

CONCLUSIONS

The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

Identification, molecular characteristics, and tissue differential expression of DGAT2 full-CDS cDNA sequence in Binglangjiang buffalo (Bubalus bubalis)

It has been found that diacylglycerol acyltransferase-2 (DGAT2) plays a crucial role in the synthesis of triglycerides (TGs) in some mammals, but its role in buffalo lactation is unclear. In the present study, the DGAT2 full-CDS cDNA sequence of Binglangjiang buffalo was isolated, and the physicochemical characteristics and structure of its encoding protein were characterized. Furthermore, the differential expressions of this gene in 10 tissues of lactating and non-lactating buffalo were analyzed by real-time quantitative PCR (RT-qPCR). The results showed that the coding region (CDS) of this gene was 1086 bp in length, encoding a peptide composed of 361 amino acid residues. The deduced amino acid sequence shared more than 98.6 % identity with that of cattle, zebu, yak, and bison in the Bovidae family. Buffalo DGAT2 protein is a slightly hydrophobic protein with a transmembrane region, which functions in membrane of endoplasmic reticulum. Besides, this protein belongs to the LPLAT_MGAT-like family and contains a conserved domain of DAGAT that has a function in the synthesis of TGs. The multi-tissue differential expression analysis demonstrated that DGAT2 was expressed in the heart, liver, mammary gland, and muscle in both non-lactating and lactating buffalo. And its expression level in the heart, liver, and mammary gland during lactation was significantly higher than that during non-lactation. The results indicate that buffalo DGAT2 may be involved in milk fat synthesis. This study can establish a foundation for further elucidating mechanisms of the buffalo DGAT2 gene in milk fat synthesis.

Genome-wide analysis reveals genetic diversity, linkage disequilibrium, and selection for milk production traits in Chinese buffalo breeds

The water buffalo is an important dual-purpose livestock that is widespread throughout central and southern China. However, there has been no characterization of the population genetics of Chinese buffalo. Using an Axiom buffalo genotyping array (Thermo Fisher Scientific, Wilmington, DE), we analyzed the genetic diversity, linkage disequilibrium pattern, and signature of selection in 176 Chinese buffaloes from 13 breeds. A total of 35,547 SNP passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between swamp and river types. Ten Chinese indigenous breeds were clustered into the swamp group, the Murrah and Nili-Ravi breeds were clustered into the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosity. Linkage disequilibrium decay distance was much shorter in the swamp group compared with the river group, with an average square of correlation coefficient value of 0.2 of approximately 50 kb. Analysis of runs of homozygosity indicated extensive remote and recent inbreeding within swamp and river groups, respectively. Moreover, one genomic region under selection was detected between the river and swamp groups. Our findings contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be used in buffalo breeding programs.

Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis

Background

Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (F_ROH), excess of homozygosity (F_HOM), correlation between uniting gametes (F_UNI), and diagonal elements of the genomic relationship matrix (F_GRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies.

Results

In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on F_HOM, F_UNI and F_GRM were higher in AZ than they were in KHZ, which was in contrast to the F_ROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX).

Conclusion

The calculated F_ROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that F_ROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.

Association of SNPs in GnRH gene with sperm quality traits of Chinese water buffalo

Hypothalamic gonadotropin-releasing hormone (GnRH) controls the activity of hypothalamic-pituitary-gonadal axis and plays a key role in the reproductive performance of animals. In this study, five single nucleotide polymorphisms (SNPs), namely g.991T > C, g.1041T > C g.3424T > C, g.3462C > A and g.3463Inde A, were detected in the GnRH gene of 162 water buffaloes by Sanger sequencing. Each SNP was associated with more than two sperm quality traits of ejaculate volume, sperm concentration, post-thaw sperm motility and sperm abnormality. g.3424T > C and g.3462C > A were related to these four traits and had a remarkable effect on ejaculate volume. The three other SNPs were related to sperm concentration, post-thaw sperm motility and sperm abnormality. Moreover, six haplotypes (H1: TCCAI, H2: CTTC-, H3: TCCCI, H4: CTTA-, H5: CCTA- and H6: CTCC-) composed of five SNPs comprising seven different combined genotypes were generated by linkage disequilibrium analysis. Statistics followed by one-way ANOVA indicated that water buffaloes with the haplotype combination H1H1 had the highest genotypic frequency, and those with the H4H4 haplotype combination had the highest ejaculate volume. The sperm concentration of those with haplotype combination H1H5 was higher than that of the other genotypes. In summary, our study showed a remarkable association between the SNPs of GnRH and sperm quality traits of Chinese water buffalo.

Genome-wide association study applied to type traits related to milk yield in water buffaloes (Bubalus bubalis)

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

Study of whole genome linkage disequilibrium patterns of Iranian water buffalo breeds using the Axiom Buffalo Genotyping 90K Array

Accuracy of genome-wide association studies, and the successful implementation of genomic selection depends on the level of linkage disequilibrium (LD) across the genome and also the persistence of LD phase between populations. In the present study LD between adjacent SNPs and LD decay between SNPs was calculated in three Iranian water buffalo populations. Persistence of LD phase was evaluated across these populations and effective population size (Ne) was estimated from corrected r² information. A set of 404 individuals from three Iranian buffalo populations were genotyped with the Axiom Buffalo Genotyping 90K Array. Average r² and |D'| between adjacent SNP pairs across all chromosomes was 0.27 and 0.66 for AZI, 0.29 and 0.68 for KHU, and 0.32 and 0.72 for MAZ. The LD between the SNPs decreased with increasing physical distance from 100Kb to 1Mb between markers, from 0.234 to 0.018 for AZI, 0.254 to 0.034 for KHU, and 0.297 to 0.119 for MAZ, respectively. These results indicate that a density of 90K SNP is sufficient for genomic analyses relying on long range LD (e.g. GWAS and genomic selection). The persistence of LD phase decreased with increasing marker distances across all the populations, but remained above 0.8 for AZI and KHU for marker distances up to 100Kb. For multi-breed genomic evaluation, the 90K SNP panel is suitable for AZI and KHU buffalo breeds. Estimated effective population sizes for AZI, KHU and MAZ were 477, 212 and 32, respectively, for recent generations. The estimated effective population sizes indicate that the MAZ is at risk and requires careful management.