Phylogeography and domestication of Chinese swamp buffalo

Maternal genetic diversity and phylogenetic analysis of Indian riverine and swamp buffaloes: insights from complete mitochondrial genomes

This study explored the genetic diversity and evolutionary history of riverine and swamp buffaloes in India, utilizing complete mitochondrial genome sequences. Through comprehensive sampling across varied agro-climatic zones, including 91 riverine buffaloes from 12 breeds and 6 non-descript populations, along with 16 swamp buffaloes of the Luit breed, this study employed next-generation sequencing techniques to map the mitogenomic landscape of these subspecies. Sequence alignments were performed with the buffalo mitochondrial reference genome to identify mitochondrial DNA (mtDNA) variations and distinct maternal haplogroups among Indian buffaloes. The results uncovered the existence of 212 variable sites in riverine buffaloes, yielding 67 haplotypes with high haplotype diversity (0.991), and in swamp buffaloes, 194 variable sites resulting in 12 haplotypes, displaying haplotype diversity of 0.950. Phylogenetic analyses elucidated the genetic relationships between Indian buffaloes and the recognized global haplogroups, categorizing Indian swamp buffaloes predominantly into the SA haplogroup. Intriguingly, the haplogroup SB2b was observed for the first time in swamp buffaloes. Conversely, riverine buffaloes conformed to established sub-haplogroups RB1, RB2, and RB3, underscoring the notion of Northwestern India as a pivotal domestication site for riverine buffaloes. The study supports the hypothesis of independent domestication events for riverine and swamp buffaloes, highlighting the critical role of genetic analysis in unraveling the complex evolutionary pathways of domestic animals. This investigation contributes to the global understanding of buffalo mitogenome diversity, offering insights into this important livestock species' domestication and dispersal patterns.

Space allowance: a tool for improving behavior, milk and meat production, and reproduction performance of buffalo in different housing systems-a review

Buffalo population has dramatically increased during the last two decades, especially in tropical and subtropical regions. Although buffalo are important milk and meat-producing animal, still practices of buffalo farming and welfare aspects are not well established. Housing system and stocking density are significant factors that affect the welfare and production of animals; however, no space allowance standards have been demonstrated for buffalo at different ages. This review article presents the following: (1) an overview of buffalo subtypes and the geographical distribution of buffalo populations and their production; (2) the effect of housing systems and space allowance on the social behavior and welfare indices; (3) the effects of space allowance on milk production and growth performance of buffalo; and (4) the relationship between space allowance and reproductive performance. Although the limited data in this area of research, it can be driven that a larger space allowance with access to a pool, especially during the hot season, maintains buffalo production at optimal levels. Moreover, optimal floor space improves the welfare and social indices of buffalo; however, there are discrepancies in aggressive and agonistic behavior results. Surprisingly, the reproductive performance of buffalo was not affected by space allowance. Therefore, further research is needed to identify the impact of the housing aspects, including space allowance and enrichment tools, on the productive performance, and welfare indices of buffalo. This would assist in implementing welfare-economic standards for buffalo production and reveal the potentiality of this eco-friendly animal.

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.

Genetic diversity of mitochondrial cytochrome b gene in swamp buffalo

The swamp buffalo (Bubalus carabanensis) is mainly bred for meat, transport and rice cultivation in China and Southeast Asian countries. In the current study, we investigated the genetic diversity, maternal origin and phylogenetic relationship of swamp buffalo by analyzing 1,786 mitochondrial cytochrome b (cytb) sequences from China, Vietnam, Laos, Thailand, India and Bangladesh. Our results indicated that the swamp buffalo can be divided into two major lineages (SA and SB with the sublineages) and three rare lineages (SC, SD and SE), which showed strong geographic differentiation. The SA1 lineage represented a major domestication event, which involved population expansion. Regions III and V showed higher haplotype diversity than the other regions, indicating that the regions of Southwest China and IndoChina are potential domestication centers for the swamp buffalo. In addition, the swamp buffalo showed a closer phylogenetic relationship to tamaraw. In conclusion, our findings revealed a high level of genetic diversity and the phylogenetic pattern of the swamp buffalo.

Mitogenome Diversity and Maternal Origins of Guangxi Buffalo Breeds

Guangxi Province, in the southwest of China, is one of the putative migratory corridors or domestication centers for swamp buffalo. In this study, we investigated the evolutionary status of two Guangxi native buffalo breeds (Fuzhong buffalo, n = 15; Xilin buffalo, n = 25) based on the complete mitogenome sequencing. Our results revealed rich genetic diversity in the two buffalo breeds. We detected five haplogroups (SA1, SA2, SB1, SB2, SB3) in the two Guangxi buffalo breeds, and the haplogroup SB3 in the Fuzhong buffalo. Our results showed that the haplogroup SA1 was associated with the major domestication event that involved population expansion in Guangxi buffalo. In conclusion, our findings revealed a high level of maternal genetic diversity and the phylogenetic pattern of the two Guangxi buffalo breeds.

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.

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.

Productive, reproductive, and estrus characteristics of different breeds of buffalo cows in Bangladesh

Objective:

The objective of this research work is to know the productive and reproductive performances and problems of local, crossbred, Nilli, and Murrah buffalo cows in selected study areas in Bangladesh.

Methodology:

A total of 1,241 local, crossbred, Nilli, and Murrah buffalo cows were surveyed in the selected areas with a pre-set questionnaire. Among 1,241 buffalo cows, 112 buffalo cows were randomly selected at day 0 of the estrus cycle for studying ovarian features.

Results:

Results showed that the average age, body condition score, and body weight were significantly (p < 0.05) different among the studied breeds. Milk production in Murrah and lactation length in Nilli cows were significantly (p < 0.05) higher than indigenous, crossbred, Nilli, and indigenous, crossbred, Murrah buffalo cows, respectively. Results also illustrated that sexual maturity, estrus cycle length, insemination time after the onset of estrus, and gestation length insignificantly (p > 0.05) varied among the surveyed breed. But, the fallout of the study denoted that estrus duration, first calving age, parity number, number of service per conception, calving interval, and voluntary waiting period varied significantly (p < 0.05) in different breeds. Ovarian physiological characteristics such as vaginal electrical resistance, average number of follicles in two ovaries, and largest follicular diameter, estrogen, and progesterone at day 0 of the estrus cycle of local, crossbred, Nilli, and Murrah buffalo cows showed insignificantly (p > 0.05) differences.

Conclusion:

The study will help the veterinarian and researcher to identify the constraints for the reproductive efficiency of buffalo in Bangladesh.

Multiple domestication of swamp buffalo in China and South East Asia

The domestication of swamp buffalo (Bubalus bubalis carabanesis) has been discussed for years and still remained unclear. To obtain a better understanding of where, when and how the swamp buffaloes were domesticated, the complete mitochondrial DNA D-loop sequences of 1,788 individuals from China, Vietnam, Laos, Thailand, Burma, Bangladesh and India were investigated. Our results revealed swamp buffalo with abundant genetic diversity. The lineage SA of swamp buffalo may be first domesticated in Lower Yangtze and then migrated following two different routes: one migrated along the Upper Yangtze to Southwest China (L1); the other migrated to Southeast China to Southeast Asia through Guangxi province (L2). During the migration process several later domestication events may be inspired in the Upper Yangtze (SB3), Southwest China and North of Southeast Asia (SB1 and SB2). In this study, we proposed a non-independent multiple domestication pattern in swamp buffalo.

Karyotypic and mtDNA based characterization of Malaysian water buffalo

Background

In Malaysia, the domestic water buffaloes (Bubalus bubalis) are classified into the swamp and the murrah buffaloes. Identification of these buffaloes is usually made via their phenotypic appearances. This study characterizes the subspecies of water buffaloes using karyotype, molecular and phylogenetic analyses. Blood of 105 buffaloes, phenotypically identified as swamp, murrah and crossbred buffaloes were cultured, terminated and harvested using conventional karyotype protocol to determine the number of chromosomes. Then, the D-loop of mitochondrial DNA of 10 swamp, 6 crossbred and 4 murrah buffaloes which were identified earlier by karyotyping were used to construct a phylogenetic tree was constructed.

Results

Karyotypic analysis confirmed that all 93 animals phenotypically identified as swamp buffaloes with 48 chromosomes, all 7 as crossbreds with 49 chromosomes, and all 5 as murrah buffaloes with 50 chromosomes. The D-loop of mitochondrial DNA analysis showed that 10 haplotypes were observed with haplotype diversity of 0.8000 ± 0.089. Sequence characterization revealed 72 variables sites in which 67 were parsimony informative sites with sequence diversity of 0.01906. The swamp and murrah buffaloes clearly formed 2 different clades in the phylogenetic tree, indicating clear maternal divergence from each other. The crossbreds were grouped within the swamp buffalo clade, indicating the dominant maternal swamp buffalo gene in the crossbreds.

Conclusion

Thus, the karyotyping could be used to differentiate the water buffaloes while genotypic analysis could be used to characterize the water buffaloes and their crossbreds.

Genome-Wide SNP Data Revealed the Extent of Linkage Disequilibrium, Persistence of Phase and Effective Population Size in Purebred and Crossbred Buffalo Populations

Linkage disequilibrium (LD) is a useful parameter for guiding the accuracy and power of both genome-wide association studies (GWAS) and genomic selection (GS) among different livestock species. The present study evaluated the extent of LD, persistence of phase and effective population size (Ne) for the purebred (Mediterranean buffalo; n = 411) and crossbred [Mediterranean × Jianghan × Nili-Ravi buffalo, n = 9; Murrah × Nili-Ravi × local (Xilin or Fuzhong) buffalo, n = 36] buffalo populations using the 90K Buffalo SNP genotyping array. The results showed that the average square of correlation coefficient (r ²) between adjacent SNP was 0.13 ± 0.19 across all autosomes for purebred and 0.09 ± 0.13 for crossbred, and the most rapid decline in LD was observed over the first 200 kb. Estimated r ² ≥ 0.2 extended up to ~50 kb in crossbred and 170 kb in purebred populations, while average r ² values ≥0.3 were respectively observed in the ~10 and 60 kb in the crossbred and purebred populations. The largest phase correlation (R _{P, C} = 0.47) was observed at the distance of 100 kb, suggesting that this phase was not actively preserved between the two populations. Estimated Ne for the purebred and crossbred population at the current generation was 387 and 113 individuals, respectively. These findings may provide useful information to guide the GS and GWAS in buffaloes.

Bayesian Statistics in Archaeology

Null hypothesis significance testing (NHST) is the most common statistical framework used by scientists, including archaeologists. Owing to increasing dissatisfaction, however, Bayesian inference has become an alternative to these methods. In this article, we review the application of Bayesian statistics to archaeology. We begin with a simple example to demonstrate the differences in applying NHST and Bayesian inference to an archaeological problem. Next, we formally define NHST and Bayesian inference, provide a brief historical overview of their development, and discuss the advantages and limitations of each method. A review of Bayesian inference and archaeology follows, highlighting the applications of Bayesian methods to chronological, bioarchaeological, zooarchaeological, ceramic, lithic, and spatial analyses. We close by considering the future applications of Bayesian statistics to archaeological research.

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.

Whole Mitogenomes Reveal the History of Swamp Buffalo: Initially Shaped by Glacial Periods and Eventually Modelled by Domestication

The newly sequenced mitochondrial genomes of 107 Asian swamp buffalo (Bubalus bubalis carabensis) allowed the reconstruction of the matrilineal divergence since ~900 Kya. Phylogenetic trees and Bayesian skyline plots suggest a role of the glacial periods in the demographic history of swamp buffalo. The ancestral swamp-buffalo mitogenome is dated ~232 ± 35 Kya. Two major macro-lineages diverged during the 2^nd Pleistocene Glacial Period (~200-130 Kya), but most (~99%) of the current matrilines derive from only two ancestors (SA1'2 and SB) that lived around the Last Glacial Maximum (~26-19 Kya). During the late Holocene optimum (11-6 Kya) lineages differentiated further, and at least eight matrilines (SA1, SA2, SB1a, SB1b, SB2a, SB2b, SB3 and SB4) were domesticated around 7-3 Kya. Haplotype distributions support an initial domestication process in Southeast Asia, while subsequent captures of wild females probably introduced some additional rare lineages (SA3, SC, SD and SE). Dispersal of domestic buffaloes created local population bottlenecks and founder events that further differentiated haplogroup distributions. A lack of maternal gene flow between neighboring populations apparently maintained the strong phylogeography of the swamp buffalo matrilines, which is the more remarkable because of an almost complete absence of phenotypic differentiation.

De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.

Origins and domestication of cultivated banana inferred from chloroplast and nuclear genes

BACKGROUND

Cultivated bananas are large, vegetatively-propagated members of the genus Musa. More than 1,000 cultivars are grown worldwide and they are major economic and food resources in numerous developing countries. It has been suggested that cultivated bananas originated from the islands of Southeast Asia (ISEA) and have been developed through complex geodomestication pathways. However, the maternal and parental donors of most cultivars are unknown, and the pattern of nucleotide diversity in domesticated banana has not been fully resolved.

METHODOLOGY/PRINCIPAL FINDINGS

We studied the genetics of 16 cultivated and 18 wild Musa accessions using two single-copy nuclear (granule-bound starch synthase I, GBSS I, also known as Waxy, and alcohol dehydrogenase 1, Adh1) and two chloroplast (maturase K, matK, and the trnL-F gene cluster) genes. The results of phylogenetic analyses showed that all A-genome haplotypes of cultivated bananas were grouped together with those of ISEA subspecies of M. acuminata (A-genome). Similarly, the B- and S-genome haplotypes of cultivated bananas clustered with the wild species M. balbisiana (B-genome) and M. schizocarpa (S-genome), respectively. Notably, it has been shown that distinct haplotypes of each cultivar (A-genome group) were nested together to different ISEA subspecies M. acuminata. Analyses of nucleotide polymorphism in the Waxy and Adh1 genes revealed that, in comparison to the wild relatives, cultivated banana exhibited slightly lower nucleotide diversity both across all sites and specifically at silent sites. However, dramatically reduced nucleotide diversity was found at nonsynonymous sites for cultivated bananas.

CONCLUSIONS/SIGNIFICANCE

Our study not only confirmed the origin of cultivated banana as arising from multiple intra- and inter-specific hybridization events, but also showed that cultivated banana may have not suffered a severe genetic bottleneck during the domestication process. Importantly, our findings suggested that multiple maternal origins and a reduction in nucleotide diversity at nonsynonymous sites are general attributes of cultivated bananas.