Design and validation of a 90K SNP genotyping assay for the water buffalo (Bubalus bubalis)

Genomic analysis of genotype-environment interaction in age at first calving of Murrah buffaloes

Age at first calving (AFC) is a measure of sexual maturity associated with the start of productive life of dairy animals. Additionally, a lower AFC reduces the generation interval and early culling of females. However, AFC has low heritability, making it a trait highly influenced by environmental factors. In this scenario, one way to improve the reproductive performance of buffalo cows is to select robust animals according to estimated breeding value (EBV) using models that include genotype-environment interaction (GEI) with the application of reaction norm models (RNMs). This can be achieved by understanding the genomic basis related to GEI of AFC. Thus, in this study, we aimed to predict EBV considering GEI via the RNM and identify candidate genes related to this component in dairy buffaloes through genome-wide association studies (GWAS). We used 1795 AFC records from three Murrah buffalo herds and formed environmental gradients (EGs) from contemporary group solutions obtained from genetic analysis of 270-day cumulative milk yield. Heritability estimates ranged from 0.15 to 0.39 along the EG. GWAS of the RNM slope parameter identified important genomic regions. The genomic window that explained the highest percentage of genetic variance of the slope (0.67%) was located on BBU1. After functional analysis, five candidate genes were detected, involved in two biological processes. The results suggested the existence of a GEI for AFC in Murrah buffaloes, with reclassification of animals when different environmental conditions were considered. The inclusion of genomic information increased the accuracy of breeding values for the intercept and slope of the reaction norm. GWAS analysis suggested that important genes associated with the AFC reaction norm slope were possibly also involved in biological processes related to lipid metabolism and immunity.

CSN1S1, CSN3 and LPL: Three Validated Gene Polymorphisms Useful for More Sustainable Dairy Production in the Mediterranean River Buffalo

The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes and in DGAT1. Nowadays, similar results have not yet been reached in river buffaloes, despite the availability of advanced genomic technologies and accurate phenotype records. The aim of the present study was to investigate and validate the effect of four single nucleotide polymorphisms (SNP) in the CSN1S1, CSN3, SCD and LPL genes on seven milk traits in a larger buffalo population. These SNPs have previously been reported to be associated with, or affect, dairy traits in smaller populations often belonging to one farm. A total of 800 buffaloes were genotyped. The following traits were individually recorded, monthly, throughout each whole lactation period from 2010 to 2021: daily milk yield (dMY, kg), protein yield (dPY, kg) and fat yield (dFY, kg), fat and protein contents (dFP, % and dPP, %), somatic cell count (SCC, 103 cell/mL) and urea (mg/dL). A total of 15,742 individual milk test day records (2496 lactations) were available for 680 buffalo cows, with 3.6 ± 1.7 parities (from 1 to 13) and an average of 6.1 ± 1.2 test day records per lactation. Three out four SNPs in the CSN1S1, CSN3 and LPL genes were associated with at least one of analyzed traits. In particular, the CSN1S1 (AJ005430:c.578C>T) gave favorable associations with all yield traits (dMY, p = 0.022; dPY, p = 0.014; dFY, p = 0.029) and somatic cell score (SCS, p = 0.032). The CSN3 (HQ677596: c.536C>T) was positively associated with SCS (p = 0.005) and milk urea (p = 0.04). Favorable effects on daily milk yield (dMY, p = 0.028), fat (dFP, p = 0.027) and protein (dPP, p = 0.050) percentages were observed for the LPL. Conversely, the SCD did not show any association with milk traits. This is the first example of a confirmation study carried out in the Mediterranean river buffalo for genes of economic interest in the dairy field, and it represents a very important indication for the preselection of young bulls destined for breeding programs aimed at more sustainable dairy production.

A genotyping array for the globally invasive vector mosquito, Aedes albopictus

BACKGROUND

Although whole-genome sequencing (WGS) is the preferred genotyping method for most genomic analyses, limitations are often experienced when studying genomes characterized by a high percentage of repetitive elements, high linkage, and recombination deserts. The Asian tiger mosquito (Aedes albopictus), for example, has a genome comprising up to 72% repetitive elements, and therefore we set out to develop a single-nucleotide polymorphism (SNP) chip to be more cost-effective. Aedes albopictus is an invasive species originating from Southeast Asia that has recently spread around the world and is a vector for many human diseases. Developing an accessible genotyping platform is essential in advancing biological control methods and understanding the population dynamics of this pest species, with significant implications for public health.

METHODS

We designed a SNP chip for Ae. albopictus (Aealbo chip) based on approximately 2.7 million SNPs identified using WGS data from 819 worldwide samples. We validated the chip using laboratory single-pair crosses, comparing technical replicates, and comparing genotypes of samples genotyped by WGS and the SNP chip. We then used the chip for a population genomic analysis of 237 samples from 28 sites in the native range to evaluate its usefulness in describing patterns of genomic variation and tracing the origins of invasions.

RESULTS

Probes on the Aealbo chip targeted 175,396 SNPs in coding and non-coding regions across all three chromosomes, with a density of 102 SNPs per 1 Mb window, and at least one SNP in each of the 17,461 protein-coding genes. Overall, 70% of the probes captured the genetic variation. Segregation analysis found that 98% of the SNPs followed expectations of single-copy Mendelian genes. Comparisons with WGS indicated that sites with genotype disagreements were mostly heterozygotes at loci with WGS read depth < 20, while there was near complete agreement with WGS read depths > 20, indicating that the chip more accurately detects heterozygotes than low-coverage WGS. Sample sizes did not affect the accuracy of the SNP chip genotype calls. Ancestry analyses identified four to five genetic clusters in the native range with various levels of admixture.

CONCLUSIONS

The Aealbo chip is highly accurate, is concordant with genotypes from WGS with high sequence coverage, and may be more accurate than low-coverage WGS.

Transcriptome-Wide Association Study Reveals Potentially Candidate Genes Responsible for Milk Production Traits in Buffalo

Identifying key causal genes is critical for unraveling the genetic basis of complex economic traits, yet it remains a formidable challenge. The advent of large-scale sequencing data and computational algorithms, such as transcriptome-wide association studies (TWASs), offers a promising avenue for identifying potential causal genes. In this study, we harnessed the power of TWAS to identify genes potentially responsible for milk production traits, including daily milk yield (MY), fat percentage (FP), and protein percentage (PP), within a cohort of 100 buffaloes. Our approach began by generating the genotype and expression profiles for these 100 buffaloes through whole-genome resequencing and RNA sequencing, respectively. Through comprehensive genome-wide association studies (GWAS), we pinpointed a total of seven and four single nucleotide polymorphisms (SNPs) significantly associated with MY and FP traits, respectively. By using TWAS, we identified 55, 71, and 101 genes as significant signals for MY, FP, and PP traits, respectively. To delve deeper, we conducted protein-protein interaction (PPI) analysis, revealing the categorization of these genes into distinct PPI networks. Interestingly, several TWAS-identified genes within the PPI network played a vital role in milk performance. These findings open new avenues for identifying potentially causal genes underlying important traits, thereby offering invaluable insights for genomics and breeding in buffalo populations.

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

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

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.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

An Integrative Genomic Prediction Approach for Predicting Buffalo Milk Traits by Incorporating Related Cattle QTLs

Background: The 90K Axiom Buffalo SNP Array is expected to improve and speed up various genomic analyses for the buffalo (Bubalus bubalis). Genomic prediction is an effective approach in animal breeding to improve selection and reduce costs. As buffalo genome research is lagging behind that of the cow and production records are also limited, genomic prediction performance will be relatively poor. To improve the genomic prediction in buffalo, we introduced a new approach (pGBLUP) for genomic prediction of six buffalo milk traits by incorporating QTL information from the cattle milk traits in order to help improve the prediction performance for buffalo. Results: In simulations, the pGBLUP could outperform BayesR and the GBLUP if the prior biological information (i.e., the known causal loci) was appropriate; otherwise, it performed slightly worse than BayesR and equal to or better than the GBLUP. In real data, the heritability of the buffalo genomic region corresponding to the cattle milk trait QTLs was enriched (fold of enrichment > 1) in four buffalo milk traits (FY270, MY270, PY270, and PM) when the EBV was used as the response variable. The DEBV as the response variable yielded more reliable genomic predictions than the traditional EBV, as has been shown by previous research. The performance of the three approaches (GBLUP, BayesR, and pGBLUP) did not vary greatly in this study, probably due to the limited sample size, incomplete prior biological information, and less artificial selection in buffalo. Conclusions: To our knowledge, this study is the first to apply genomic prediction to buffalo by incorporating prior biological information. The genomic prediction of buffalo traits can be further improved with a larger sample size, higher-density SNP chips, and more precise prior biological information.

Meprin A1 subunit beta gene polymorphism is associated with the length of post-partum anestrus interval in Murrah buffaloes

Postpartum anestrus interval (PPAI) is the interval between parturition and the first postpartum estrus exhibition in animals. Appearance of both normal and PPA buffaloes under the same farm conditions indicates the role of possible genetic predisposition to PPA. To identify the genetic and non-genetic factors associated with PPA in buffaloes, we collected data on PPAI and other 38 non-genetic variables from 575 Murrah buffaloes in the field conditions and identified the PPA associated non-genetic factors in our previous study. To explore the genetic factors associated with the unexplained variation in PPAI residuals, the present study identified 41 single nucleotide polymorphisms (SNPs) in 13 candidate genes using Sanger sequencing. Exploration of their association with the PPAI residuals of 50 extreme PPA and 50 normal buffaloes identified the significant (P < 0.01) association of the SNP (g.37219977A>G) in the 3'-UTR region of the Meprin A 1 subunit beta (Mep1b) gene with PPAI, which was further validated (P = 0.058) in a large population sample (n = 417). Bioinformatics analysis of the 3'-UTR region has identified three miRNA, bta-miR-2420, bta-miR-2325b and bta-miR-453 that could regulate Igf-1 in the plasma of animals with different genotypes (GG, AG and AA). The higher Igf-1 levels in the GG genotypes than that of AA and AG genotypes of this SNP (g.37219977A>G) further suggest the association of Mep1b gene with PPA condition in Murrah buffaloes. As a result of this study, we propose that buffaloes with protective alleles at this SNP be selected to improve the herd's reproductive efficiency.

Whole-Genome-Based Web Genomic Resource for Water Buffalo (Bubalus bubalis)

Water buffalo (Bubalus bubalis), belonging to the Bovidae family, is an economically important animal as it is the major source of milk, meat, and drought in numerous countries. It is mainly distributed in tropical and subtropical regions with a global population of approximately 202 million. The advent of low cost and rapid sequencing technologies has opened a new vista for global buffalo researchers. In this study, we utilized the genomic data of five commercially important buffalo breeds, distributed globally, namely, Mediterranean, Egyptian, Bangladesh, Jaffrarabadi, and Murrah. Since there is no whole-genome sequence analysis of these five distinct buffalo breeds, which represent a highly diverse ecosystem, we made an attempt for the same. We report the first comprehensive, holistic, and user-friendly web genomic resource of buffalo (BuffGR) accessible at http://backlin.cabgrid.res.in/buffgr/, that catalogues 6028881 SNPs and 613403 InDels extracted from a set of 31 buffalo tissues. We found a total of 7727122 SNPs and 634124 InDels distributed in four breeds of buffalo (Murrah, Bangladesh, Jaffarabadi, and Egyptian) with reference to the Mediterranean breed. It also houses 4504691 SSR markers from all the breeds along with 1458 unique circRNAs, 37712 lncRNAs, and 938 miRNAs. This comprehensive web resource can be widely used by buffalo researchers across the globe for use of markers in marker trait association, genetic diversity among the different breeds of buffalo, use of ncRNAs as regulatory molecules, post-transcriptional regulations, and role in various diseases/stresses. These SNPs and InDelscan also be used as biomarkers to address adulteration and traceability. This resource can also be useful in buffalo improvement programs and disease/breed management.

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

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.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN)

BACKGROUND

German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement.

RESULTS

Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%.

CONCLUSION

The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.

Accuracy of Genomic Prediction for Milk Production Traits in Philippine Dairy Buffaloes

The objective of this study was to compare the accuracies of genomic prediction for milk yield, fat yield, and protein yield from Philippine dairy buffaloes using genomic best linear unbiased prediction (GBLUP) and single-step GBLUP (ssGBLUP) with the accuracies based on pedigree BLUP (pBLUP). To also assess the bias of the prediction, the regression coefficient (slope) of the adjusted phenotypes on the predicted breeding values (BVs) was also calculated. Two data sets were analyzed. The GENO data consisting of all female buffaloes that have both phenotypes and genotypes (n = 904 with 1,773,305-days lactation records) were analyzed using pBLUP and GBLUP. The ALL data, consisting of the GENO data plus females with phenotypes but not genotyped (n = 1,975 with 3,821,305-days lactation records), were analyzed using pBLUP and ssGBLUP. Animals were genotyped with the Affymetrix 90k buffalo genotyping array. After quality control, 60,827 single-nucleotide polymorphisms were used for downward analysis. A pedigree file containing 2,642 animals was used for pBLUP and ssGBLUP. Accuracy of prediction was calculated as the correlation between the predicted BVs of the test set and adjusted phenotypes, which were corrected for fixed effects, divided by the square root of the heritability of the trait, corrected for the number of lactations used in the test set. To assess the bias of the prediction, the regression coefficient (slope) of the adjusted phenotypes on the predicted BVs was also calculated. Results showed that genomic methods (GBLUP and ssGBLUP) provide more accurate predictions compared to pBLUP. Average GBLUP and ssGBLUP accuracies were 0.24 and 0.29, respectively, whereas average pBLUP accuracies (for GENO and ALL data) were 0.21 and 0.22, respectively. Slopes of the two genomic methods were also closer to one, indicating lesser bias, compared to pBLUP. Average GBLUP and ssGBLUP slopes were 0.89 and 0.84, respectively, whereas the average pBLUP (for GENO and ALL data) slopes were 0.80 and 0.54, respectively.

Genome-Wide Association Studies in Indian Buffalo Revealed Genomic Regions for Lactation and Fertility

Murrah breed of buffalo is an excellent dairy germplasm known for its superior milk quality in terms of milk fat and solids-not-fat (SNF); however, it is often reported that Indian buffaloes had lower lactation and fertility potential compared to the non-native cattle of the country. Recent techniques, particularly the genome-wide association studies (GWAS), to identify genomic variations associated with lactation and fertility traits offer prospects for systematic improvement of buffalo. DNA samples were sequenced using the double-digestion restriction-associated DNA (RAD) tag genotyping-by-sequencing. The bioinformatics pipeline was standardized to call the variants, and single-nucleotide polymorphisms (SNPs) qualifying the stringent quality check measures were retained for GWAS. Over 38,000 SNPs were used to perform GWAS on the first two principal components of test-day records of milk yields, fat percentages, and SNF percentages, separately. GWAS was also performed on 305 days’ milk yield; lactation persistency was estimated through the rate of decline after attaining the peak yield method, along with three other standard methods; and breeding efficiency, post-partum breeding interval, and age at sexual maturity were considered fertility traits. Significant association of SNPs was observed for the first principal component, explaining the maximum proportion of variation in milk yield. Furthermore, some potential genomic regions were identified to have a potential role in regulating milk yield and fertility in Murrah. Identification of such genomic regions shall help in carrying out an early selection of high-yielding persistent Murrah buffaloes and, in the long run, would be helpful in shaping their future genetic improvement programs.

An updated Axiom buffalo genotyping array map and mapping of cattle quantitative trait loci to the new water buffalo reference genome assembly

The objectives of this study were to provide the buffalo research community with an updated SNP map for the Axiom Buffalo Genotyping (ABG) array with genomic positions for SNP currently unmapped and to map all cattle QTL from the CattleQTLdb onto the buffalo reference assembly. To update the ABG array map, all SNP probe sequences from the ABG array were re-aligned against the UOA_WB_1 assembly. With the new map, the number of mapped markers increased by approximately 10% and went from 106 778 to 116 708, which reduced the average marker spacing by approximately 2 kb. A comparison of results between signatures of autozygosity study using the ABG and the new map showed that, when the additional markers were used there was an increase in the autozygosity peaks and additional peaks in BBU5 and BBU11 could be identified. After sequence alignment and quality control, 64 650 (UMD3.1) and 76 530 (ARS_UCD1.2) cattle QTL were mapped onto the buffalo genome. The mapping of the bovine QTL database onto the buffalo genome should be useful for genome-wide association studies in buffalo and, given the high homology between the two species, the positions of cattle QTL on the buffalo genome can serve as a stepping stone towards a water buffalo QTL database.

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.

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

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.

Opportunities and Challenges for Improving the Productivity of Swamp Buffaloes in Southeastern Asia

The swamp buffalo is a domesticated animal commonly found in Southeast Asia. It is a highly valued agricultural animal for smallholders, but the production of this species has unfortunately declined in recent decades due to rising farm mechanization. While swamp buffalo still plays a role in farmland cultivation, this species' purposes has shifted from draft power to meat, milk, and hide production. The current status of swamp buffaloes in Southeast Asia is still understudied compared to its counterparts such as the riverine buffaloes and cattle. This review discusses the background of swamp buffalo, with an emphasis on recent work on this species in Southeast Asia, and associated genetics and genomics work such as cytogenetic studies, phylogeny, domestication and migration, genetic sequences and resources. Recent challenges to realize the potential of this species in the agriculture industry are also discussed. Limited genetic resource for swamp buffalo has called for more genomics work to be done on this species including decoding its genome. As the economy progresses and farm mechanization increases, research and development for swamp buffaloes are focused on enhancing its productivity through understanding the genetics of agriculturally important traits. The use of genomic markers is a powerful tool to efficiently utilize the potential of this animal for food security and animal conservation. Understanding its genetics and retaining and maximizing its adaptability to harsher environments are a strategic move for food security in poorer nations in Southeast Asia in the face of climate change.

Genetic Association of PPARGC1A Gene Single Nucleotide Polymorphism with Milk Production Traits in Italian Mediterranean Buffalo

PPARGC1A gene plays an important role in the activation of various important hormone receptors and transcriptional factors involved in the regulation of adaptive thermogenesis, gluconeogenesis, fiber-type switching in skeletal muscle, mitochondrial biogenesis, and adipogenesis, regulating the reproduction and proposed as a candidate gene for milk-related traits in cattle. This study identified polymorphisms in the PPARGC1A gene in Italian Mediterranean buffaloes and their associations to milk production and quality traits (lactation length, peak milk yield, fat and protein yield, and percentage). As a result, a total of seven SNPs (g.-78A>G, g.224651G>C, g.286986G>A, g.304050G>A, g.325647G>A, g.325817T>C, and g.325997G>A) were identified by DNA pooled sequencing. Analysis of productivity traits within the genotyped animals revealed that the g.286986G>A located at intron 4 was associated with milk production traits, but the g.325817T>C had no association with milk production. Polymorphisms in g.-78A>G was associated with peak milk yield and milk yield, while g.304050G>A and g.325997 G>A were associated with both milk yield and protein percentage. Our findings suggest that polymorphisms in the buffalo PPARGC1A gene are associated with milk production traits and can be used as a candidate gene for milk traits and marker-assisted selection in the buffalo breeding program.

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.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Genetic Diversity of Bubalus bubalis in Germany and Global Relations of Its Genetic Background

This is the first study to explore the genetic diversity and population structure of domestic water buffalo (Bubalus bubalis) in Germany and their potential relations to herds in other parts of Europe or worldwide. To this end, animals from different herds in Germany, Bulgaria, Romania, and Hungary were genotyped and compared to genotypes from other populations with worldwide distribution and open to the public. The pilot study analyzed population structure, phylogenetic tree, and inbreeding events in our samples. In buffalos from Germany, a mixed genetic make-up with contributions from Bulgaria (Murrah breed), Romania, and Italy was found. All in all, a high degree of genetic diversity was identified in European buffalos, and a novel genotype was described in Hungarian buffalos by this study. We demonstrate that European buffalos stand out from other buffalo populations worldwide, supporting the idea that buffalos have not completely disappeared from the European continent during the late Pleistocene. The high genetic diversity in European buffalos seems to be an excellent prerequisite for the establishment of local breeds characterized by unique traits and features. This study may also be considered as an initial step on the way to genome characterization for the sustainable development of the buffalo economy in Germany and other parts of Europe in the future.

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.

Prospecting genomic regions associated with milk production traits in Egyptian buffalo

The objectives of the current study were to detect putative genomic loci and to identify candidate genes associated with milk production traits in Egyptian buffalo. A total number of 161 479 daily milk yield (DMY) records and 60 318 monthly measures for fat and protein percentages (FP and PP, respectively), along with fat and protein yields (FY and PY, respectively) from 1670 animals were used. Genotyping was performed using Axiom® Buffalo Genotyping 90 K array. Genome-wide association study (GWAS) for each trait was performed using PLINK. After Bonferroni correction, 47 SNPs were associated with one or more milk production traits. These SNPs were distributed over 36 quantitative trait loci (QTL) and located on 20 buffalo chromosomes (BBU). For the 47 SNPs, one was overlapped for three traits (DMY, FY, and PY), six were associated with two traits (one for PP and PY and five for FY and PY) while the rest were associated with only one trait. Out of 36 identified QTL, eleven were overlapped with previously reported loci in buffalo and/or cattle populations. Some of these SNPs are placed within or close to potential candidate genes, for example: TPD52, ZBTB10, RALYL and SNX16 on BBU15, ADGRD1 on BBU17, ESRRG on BBU5 and GRIP1 on BBU4. This is the first reported study between genome-wide markers and milk components in Egyptian buffalo. Our findings provide useful information to explore the genetic mechanisms and relevant genes contributing to the variation in milk production traits. Further confirmation studies with larger population size are necessary to validate the findings and detect the causal genetic variants.

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.

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.

Whole genome analysis of water buffalo and global cattle breeds highlights convergent signatures of domestication

More people globally depend on the water buffalo than any other domesticated species, and as the most closely related domesticated species to cattle they can provide important insights into the shared evolutionary basis of domestication. Here, we sequence the genomes of 79 water buffalo across seven breeds and compare patterns of between breed selective sweeps with those seen for 294 cattle genomes representing 13 global breeds. The genomic regions under selection between cattle breeds significantly overlap regions linked to stature in human genetic studies, with a disproportionate number of these loci also shown to be under selection between water buffalo breeds. Investigation of potential functional variants in the water buffalo genome identifies a rare example of convergent domestication down to the same mutation having independently occurred and been selected for across domesticated species. Cross-species comparisons of recent selective sweeps can consequently help identify and refine important loci linked to domestication.

The comparative genomics of domesticated lineages can yield insights into the signatures of artificial selection. This study sequences 79 water buffalo genomes from 7 breeds and reveals examples of convergent domestication at the genetic level between water buffalo and cattle.

Linkage Disequilibrium-Based Inference of Genome Homology and Chromosomal Rearrangements Between Species

The aim of this study was to analyze the genomic homology between cattle (Bos taurus) and buffaloes (Bubalus bubalis) and to propose a rearrangement of the buffalo genome through linkage disequilibrium analyses of buffalo SNP markers referenced in the cattle genome assembly and also compare it to the buffalo genome assembly. A panel of bovine SNPs (single nucleotide polymorphisms) was used for hierarchical, non-hierarchical and admixture cluster analyses. Thus, the linkage disequilibrium information between markers of a specific panel of buffalo was used to infer chromosomal rearrangement. Haplotype diversity and imputation accuracy of the submetacentric chromosomes were also analyzed. The genomic homology between the species enabled us to use the bovine genome assembly to recreate a buffalo genomic reference by rearranging the submetacentric chromosomes. The centromere of the submetacentric chromosomes exhibited high linkage disequilibrium and low haplotype diversity. It allowed hypothesizing about chromosome evolution. It indicated that buffalo submetacentric chromosomes are a centric fusion of ancestral acrocentric chromosomes. The chronology of fusions was also suggested. Moreover, a linear regression between buffalo and cattle rearranged assembly and the imputation accuracy indicated that the rearrangement of the chromosomes was adequate. When using the bovine reference genome assembly, the rearrangement of the buffalo submetacentric chromosomes could be done by SNP BTA (chromosome of Bos taurus) calculations: shorter BTA (shorter arm of buffalo chromosome) was given as [(shorter BTA length - SNP position in shorter BTA)] and larger BTA length as [shorter BTA length + (larger BTA length - SNP position in larger BTA)]. Finally, the proposed linkage disequilibrium-based method can be applied to elucidate other chromosomal rearrangement events in other species with the possibility of better understanding the evolutionary relationship between their genomes.

Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High-throughput high-density genotyping arrays continue to be a fast, accurate, and cost-effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high-density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole-genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives-honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.

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

Collagens, as extracellular matrix proteins, support cells for structural integrity and contribute to support mammary basic structure and development. This study aims to perform the genomic identification, evolution, and expression analyses of the collagen gene family in water buffalo (Bubalus bubalis) during lactation. A total of 128 buffalo collagen protein sequences were deduced from the 45 collagen genes identified in silico from buffalo genome, which classified into six groups based on their phylogenetic relationships, conserved motifs, and gene structure analyses. The identified collagen sequences were unequally distributed on 16 chromosomes. The tandem duplicated genes were found within three chromosomes, while only one segmental event occurred between Chr3 and Chr8. Collinearity analysis revealed that a total of 36 collagen gene pairs were orthologous between buffalo and cattle genomes despite having different chromosome numbers. Comparative transcription analyses revealed that a total of 23 orthologous collagen genes were detected in the milk samples at different lactation periods between the two species. Notably, the duplicated gene pair of COL4A1-COL4A2 during lactation had a higher mRNA expression level than that of cattle, while a higher expression level of COL6A1-COL6A2 pair was found in cattle compared with that of buffalo. The present study provides useful information for investigating the potential functions of the collagen family in buffalo during lactation and helps in the functional characterization of collagen genes in additional research.

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.

Asian water buffalo: domestication, history and genetics

The domestic Asian water buffalo (Bubalus bubalis) is found on all five continents, with a global population of some 202 million. The livelihoods of more people depend on this species than on any other domestic animal. The two distinct types (river and swamp) descended from different wild Asian water buffalo (Bubalus arnee) populations that diverged some 900 kyr BP and then evolved in separate geographical regions. After domestication in the western region of the Indian subcontinent (ca. 6300 years BP), the river buffalo spread west as far as Egypt, the Balkans and Italy. Conversely, after domestication in the China/Indochina border region ca. 3000-7000 years BP, swamp buffaloes dispersed through south-east Asia and China as far as the Yangtze River valley. Molecular and morphological evidence indicates that swamp buffalo populations have strong geographic genetic differentiation and a lack of gene flow, but strong phenotypic uniformity. In contrast, river buffalo populations show a weaker phylogeographic structure, but higher phenotypic diversity (i.e. many breeds). The recent availability of a high-quality reference genome and of a medium-density marker panel for genotyping has triggered a number of genome-wide investigations on diversity, evolutionary history, production traits and functional elements. The growing molecular knowledge combined with breeding programmes should pave the way to improvements in production, environmental adaptation and disease resistance in water buffalo populations worldwide.

Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it's uniqueness with riverine buffaloes

The present investigation was focused to study genomic diversity of Indian swamp buffalo populations through reduced representation approach (ddRAD). The heterozygosity (F_ST) among the swamp buffaloes was 0.11 between Assam and Manipuri; 0.20 between swamp (Manipuri) and riverine buffaloes; 0.30 between swamp (Manipuri) and cattle. The average observed and expected heterozygosity in swamp buffalo populations was 0.254 and 0.221 respectively. The Inbreeding coefficient (F_IS) value was 0.02 among the swamp buffaloes. PCA and structure analysis revealed Manipuri swamp buffalo was genetically distinct and closely related to Nagaland swamp buffalo and least to Assam swamp buffalo. Identification of selective sweeps revealed 1087 regions to have undergone selection related to immune response, adaptation and nervous system. A total of 3451 SSRs were identified in the genome of swamp buffaloes. The study evidenced the genomic diversity in the swamp buffalo populations and its uniqueness in comparison with riverine buffalo and cattle.

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.

A new chicken 55K SNP genotyping array

BACKGROUND

China has the richest local chicken breeding resources in the world and is the world's second largest producer of meat-type chickens. Development of a moderate-density SNP array for genetic analysis of chickens and breeding of meat-type chickens taking utility of those resources is urgently needed for conventional farms, breeding industry, and research areas.

RESULTS

Eight representative local breeds or commercial broiler lines with 3 pools of 48 individuals within each breed/line were sequenced and supplied the major SNPs resource. There were 7.09 million - 9.41 million SNPs detected in each breed/line. After filtering using multiple criteria such as preferred incorporation of trait-related SNPs and uniformity of distribution across the genome, 52.18 K SNPs were selected in the final array. It consists of: (i) 19.22 K SNPs from the genomes of yellow-feathered, cyan-shank partridge and white-feathered chickens; (ii) 5.98 K SNPs related to economic traits from the Illumina 60 K SNP Bead Chip, which were found as significant associated SNPs with 15 traits in a Beijing-You crossed Cobb F2 resource population by genome-wide association study analysis; (iii) 7.63 K SNPs from 861 candidate genes of economic traits; (iv) the 0.94 K SNPs related to residual feed intake; and (v) 18.41 K from chicken SNPdb. The polymorphisms of 9 extra local breeds and 3 commercial lines were examined with this array, and 40 K - 47 K SNPs were polymorphic (with minor allele frequency > 0.05) in those breeds. The MDS result showed that those breeds can be clearly distinguished by this newly developed genotyping array.

CONCLUSIONS

We successfully developed a 55K genotyping array by using SNPs segregated from typical local breeds and commercial lines. Compared to the existing Affy 600 K and Illumina 60 K arrays, there were 21,41 K new SNPs included on our Affy 55K array. The results of the 55K genotyping data can therefore be imputed to high-density SNPs genotyping data. The array offers a wide range of potential applications such as genomic selection breeding, GWAS of interested traits, and investigation of diversity of different chicken breeds.

Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis)

Water buffalo (Bubalus bubalis) is of great economic importance as a provider of milk and meat in many countries. However, the milk yield of buffalo is much lower than that of Holstein cows. Selection of candidate genes related to milk production traits can be applied to improve buffalo milk performance. A systematic review of studies of these candidate genes will be greatly beneficial for researchers to timely and efficiently understand the research development of molecular markers for buffalo milk production traits. Here, we identified and classified the candidate genes associated with buffalo milk production traits. A total of 517 candidate genes have been identified as being associated with milk performance in different buffalo breeds. Nineteen candidate genes containing 47 mutation sites have been identified using the candidate gene approach. In addition, 499 candidate genes have been identified in six genome-wide association studies (GWASes) including two studies performed with the bovine SNP chip and four studies with the buffalo SNP chip. Genes CTNND2 (catenin delta 2), APOB (apolipoprotein B), FHIT (fragile histidine triad) and ESRRG (estrogen related receptor gamma) were identified in at least two GWASes. These four genes, especially APOB, deserve further study to explore regulatory roles in buffalo milk production. With growth in the number of buffalo genomic studies, more candidate genes associated with buffalo milk production traits will be identified. Therefore, future studies, such as those investigating gene location and functional analyses, are necessary to facilitate the exploitation of genetic potential and the improvement of buffalo milk performance.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis )

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log₂ ratios between copy number estimates and the log₂ ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Integrative Analysis of Transcriptome and GWAS Data to Identify the Hub Genes Associated With Milk Yield Trait in Buffalo

The mammary gland is the production organ in mammals that is of great importance for milk production and quality. However, characterization of the buffalo mammary gland transcriptome and identification of the valuable candidate genes that affect milk production is limited. Here, we performed the differential expressed genes (DEGs) analysis of mammary gland tissue on day 7, 50, 140, and 280 after calving and conducted gene-based genome-wide association studies (GWAS) of milk yield in 935 Mediterranean buffaloes. We then employed weighted gene co-expression network analysis (WGCNA) to identify specific modules and hub genes related to milk yield based on gene expression profiles and GWAS data. The results of the DEGs analysis showed that a total of 1,420 DEGs were detected across different lactation points. In the gene-based analysis, 976 genes were found to have genome-wide association (P ≤ 0.05) that could be defined as the nominally significant GWAS geneset (NSGG), 9 of which were suggestively associated with milk yield (P < 10⁻⁴). Using the WGCNA analysis, 544 and 225 genes associated with milk yield in the turquoise module were identified from DEGs and NSGG datasets, respectively. Several genes (including BNIPL, TUBA1C, C2CD4B, DCP1B, MAP3K5, PDCD11, SRGAP1, GDPD5, BARX2, SCARA3, CTU2, and RPL27A) were identified and considered as the hub genes because they were involved in multiple pathways related to milk production. Our findings provide an insight into the dynamic characterization of the buffalo mammary gland transcriptome, and these potential candidate genes may be valuable for future functional characterization of the buffalo mammary gland.

Chromosome-level assembly of the water buffalo genome surpasses human and goat genomes in sequence contiguity

Rapid innovation in sequencing technologies and improvement in assembly algorithms have enabled the creation of highly contiguous mammalian genomes. Here we report a chromosome-level assembly of the water buffalo (Bubalus bubalis) genome using single-molecule sequencing and chromatin conformation capture data. PacBio Sequel reads, with a mean length of 11.5 kb, helped to resolve repetitive elements and generate sequence contiguity. All five B. bubalis sub-metacentric chromosomes were correctly scaffolded with centromeres spanned. Although the index animal was partly inbred, 58% of the genome was haplotype-phased by FALCON-Unzip. This new reference genome improves the contig N50 of the previous short-read based buffalo assembly more than a thousand-fold and contains only 383 gaps. It surpasses the human and goat references in sequence contiguity and facilitates the annotation of hard to assemble gene clusters such as the major histocompatibility complex (MHC).

Despite technological advances, chromosome-level assemblies of mammalian genomes are still rare. Here, the authors use PacBio, Chicago and Hi-C approaches to generate a highly contiguous and partially-phased genome assembly for the water buffalo, Bubalus bubalis

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

Background

Identification of genomic regions that have been targets of selection may shed light on the genetic history of livestock populations and help to identify variation controlling commercially important phenotypes. The Azeri and Kuzestani buffalos are the most common indigenous Iranian breeds which have been subjected to divergent selection and are well adapted to completely different regions. Examining the genetic structure of these populations may identify genomic regions associated with adaptation to the different environments and production goals.

Results

A set of 385 water buffalo samples from Azeri (N = 262) and Khuzestani (N = 123) breeds were genotyped using the Axiom® Buffalo Genotyping 90 K Array. The unbiased fixation index method (F_ST) was used to detect signatures of selection. In total, 13 regions with outlier F_ST values (0.1%) were identified. Annotation of these regions using the UMD3.1 Bos taurus Genome Assembly was performed to find putative candidate genes and QTLs within the selected regions. Putative candidate genes identified include FBXO9, NDFIP1, ACTR3, ARHGAP26, SERPINF2, BOLA-DRB3, BOLA-DQB, CLN8, and MYOM2.

Conclusions

Candidate genes identified in regions potentially under selection were associated with physiological pathways including milk production, cytoskeleton organization, growth, metabolic function, apoptosis and domestication-related changes include immune and nervous system development. The QTL identified are involved in economically important traits in buffalo related to milk composition, udder structure, somatic cell count, meat quality, and carcass and body weight.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4759-x) contains supplementary material, which is available to authorized users.

A high-density SNP chip for genotyping great tit (Parus major) populations and its application to studying the genetic architecture of exploration behaviour

High-density SNP microarrays ("SNP chips") are a rapid, accurate and efficient method for genotyping several hundred thousand polymorphisms in large numbers of individuals. While SNP chips are routinely used in human genetics and in animal and plant breeding, they are less widely used in evolutionary and ecological research. In this article, we describe the development and application of a high-density Affymetrix Axiom chip with around 500,000 SNPs, designed to perform genomics studies of great tit (Parus major) populations. We demonstrate that the per-SNP genotype error rate is well below 1% and that the chip can also be used to identify structural or copy number variation. The chip is used to explore the genetic architecture of exploration behaviour (EB), a personality trait that has been widely studied in great tits and other species. No SNPs reached genomewide significance, including at DRD4, a candidate gene. However, EB is heritable and appears to have a polygenic architecture. Researchers developing similar SNP chips may note: (i) SNPs previously typed on alternative platforms are more likely to be converted to working assays; (ii) detecting SNPs by more than one pipeline, and in independent data sets, ensures a high proportion of working assays; (iii) allele frequency ascertainment bias is minimized by performing SNP discovery in individuals from multiple populations; and (iv) samples with the lowest call rates tend to also have the greatest genotyping error rates.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.