Genetic differentiation of water buffalo (Bubalus bubalis) populations in China, Nepal and south-east Asia: inferences on the region of domestication of the swamp buffalo

Data from three published studies of genetic variation at 18 microsatellite loci in water buffalo populations in China (18 swamp type, two river type), Nepal (one wild, one domestic river, one hybrid) and south-east Asia (eight swamp, three river) were combined so as to gain a broader understanding of genetic relationships among the populations and their demographic history. Mean numbers of alleles and expected heterozygosities were significantly different among populations. Estimates of θ (a measure of population differentiation) were significant among the swamp populations for all loci and among the river populations for most loci. Differentiation among the Chinese swamp populations (which was due primarily to just one population) was much less than among the south-east Asian. The Nepal wild animals, phenotypically swamp type but genetically like river type, are significantly different from all the domestic river populations and presumably represent the ancestral Bubalus arnee (possibly with some river-type introgression). Relationships among the swamp populations (D(A) genetic distances, principal component analysis and structure analyses) show the south-east Asian populations separated into two groups by the Chinese populations. Given these relationships and the patterns of genetic variability, we postulate that the swamp buffalo was domesticated in the region of the far south of China, northern Thailand and Indochina. Following domestication, it spread south through peninsular Malaysia to Sumatra, Java and Sulawesi, and north through China, and then to Taiwan, the Philippines and Borneo.     

Genetic variation within and relationships among populations of Asian water buffalo (Bubalus bubalis)

Genetic variation at 53 protein-coding loci (25 polymorphic) was analysed for 17 water buffalo populations – 12 swamp, three Lankan and two of the Murrah breed (river type), to determine the magnitude of genetic differentiation and the genetic relationships among the populations. In accord with previous cytological studies, the Lankan buffalo clearly are river type. Significant deviations from Hardy–Weinberg equilibrium were shown for a number of locus–population combinations, with all populations but one showing significant heterogeneity in these deviations among loci. By contrast, heterogeneity among populations for each locus was much less, indicating locus-specific deviations, which suggest selection affecting allele frequencies at some loci. There was significant genetic differentiation among populations of both the swamp and river types. The differentiation among the swamp populations may reflect the geography of south-east Asia and the presumed spread of the swamp buffalo through this region. Phylogenies derived from pairwise genetic distance estimates show the clear separation of swamp and river types, but the topology of the swamp populations shows rather poor consistency with their geographic locations. For at least one population (Australia), it is clear that bottleneck effects have distorted the phylogenetic topology. Average genetic distances for both the swamp and river types, as compared with previous studies of livestock breeds, show that the genetic differentiation of each of these sets of populations is of the same order of magnitude as that among well-recognized and established breeds of other species.     

Genetic diversity and differentiation of Chinese domestic buffalo based on 30 microsatellite markers

To determine genetic diversity and evolutionary relationships among Chinese buffalo populations, 18 indigenous swamp buffalo populations and two introduced river buffalo breeds were genotyped for 30 microsatellite loci. The mean number of alleles across the 30 loci was 8.13, and the expected heterozygosity ranged from 0.517 (Yanjin) to 0.609 (Dehong). Although there was abundant genetic variation, genetic differentiation between Chinese buffalo populations was low, with only 2.8% of the total genetic variance among populations. The genetic differentiation pattern and genetic relationships among Chinese buffalo populations were consistent with their geographic distribution. The Dehong population was discerned as a distinct indigenous population, but suffered genetic admixture with river buffalo. The remaining populations were divided into four major clusters, i.e. the Upper and Middle Reaches of Yangtze Valley cluster (Guizhou, Guizhoubai, Yanjin, Fuling, Enshi and Jianghan), the Lower Reaches of Yangtze Valley cluster (Haizi, Shanqu and Dongliu), the South China cluster (Fuan and Xinfeng) and the Southwest China cluster (Xinglong, Xilin, Diandongnan and Dechang).     

Development and characterization of microsatellite markers to study population genetic diversity of Przewalski's naked carp Gymnocypris przewalskii in China

Two populations of Przewalski's naked carp Gymnocypris przewalskii, 30 individuals per population, were screened for 10 microsatellite loci. Moderate allele variation was found in these loci with two to eight alleles per locus. The expected and observed heterozygosity ranged from 0·019 to 0·805 and from 0·160 to 0·575, respectively.     

Genetic diversity of Asian water buffalo (Bubalus bubalis): mitochondrial DNA D-loop and cytochrome b sequence variation

Swamp and river buffalo mitochondrial DNA (mtDNA) was sequenced for 303 bp of the cytochrome b gene for 54 animals from 14 populations, and for 158 bp of the D-loop region for 80 animals from 11 populations. Only one cytochrome b haplotype was found in river buffalo. Of the four haplotypes identified in swamp buffalo, one found in all populations is apparently ancestral both to the other swamp haplotypes and to the river haplotype. The phylogenetic relationships among the 33 D-loop haplotypes, with a cluster of 11 found in swamp buffalo only, also support the evolution of domesticated swamp and river buffalo from an ancestral swamp-like animal, most likely represented today by the wild Asian buffalo ( Bubalus arnee ). The time of divergence of the swamp and river types, estimated from the D-loop data, is 28 000 to 87 000 years ago. We hypothesise that the species originated in mainland south-east Asia, and that it spread north to China and west to the Indian subcontinent, where the rive type evolved and was domesticated. Following domestication in China, the domesticated swamp buffalo spread through two separate routes, through Taiwan and the Philippines to the eastern islands of Borneo and Sulawesi, and south through mainland southeast Asia and then to the western islands of Indonesia.     

Development and characterization of 14 microsatellite loci in the Neotropical fish Geophagus brasiliensis (Perciformes,Cichlidae)

Fourteen polymorphic microsatellite loci were isolated and characterized for the Neotropical cichlid Geophagus brasiliensis and tested on 30 individuals belonging to a single population. Among the 14 loci described, four showed potential presence of null alleles, inferred from the excess of homozygous genotypes, and three of these loci showed significant deviations from Hardy–Weinberg equilibrium. Fifty‐nine different alleles were detected (ranging from two to eight alleles per locus), with estimates of observed and expected heterozygosity ranging from 0·167 to 0·700 and from 0·269 to 0·825. Cross‐amplification of primers was successful in five other cichlid species.     

Isolation and characterization of tri‐ and tetra‐repeat microsatellite loci in the white‐spotted charr Salvelinus leucomaenis (Salmonidae)

Tri‐ and tetra‐motif repeat microsatellite marker loci were developed for the white‐spotted charr Salvelinus leucomaenis. The 454 pyrosequencing was used to discover repeat arrays, and eight microsatellite‐primer sets, available for the estimation of polymorphisms, were identified. The number of alleles in a wild population ranged from two to four and the observed and expected heterozygosities were 0·180–0·600 and 0·188–0·599, respectively.     

Population genetic structure and comparative diversity of smallmouth bass Micropterus dolomieu: congruent patterns from two genomes

Genetic diversity and divergence patterns of smallmouth bass Micropterus dolomieu spawning groups are analysed across its northern native range with mtDNA cytochrome b gene sequences and eight unlinked nuclear DNA microsatellite loci. Results reveal high levels of genetic variability and significant differences in allelic representation among populations (mtDNA: mean ± s.e ., H_D = 0·50 ± 0·06, mean ± s.e ., θ_ST = 0·41 ± 0·02 and microsatellites: mean ± s.e . H_O = 0·46 ± 0·03, mean ± s.e . θ_ST = 0·25 ± 0·01). The distributions of 28 variant mtDNA haplotypes, which differ by an average of 3·94 nucleotides (range = 1–8), denote divergent representation among geographic areas. Microsatellite data support nine primary population groups, whose high self‐assignment probabilities likewise display marked divergence. Genetic patterns demonstrate: (1) high genetic diversity in both genomes, (2) significant divergence among populations, probably resulting from natal site homing and low lifetime migration, (3) support for three post‐glacial refugia that variously contributed to the current northern populations, which remain evident today despite waterway connectivity and (4) a weak yet significant genetic isolation by geographic distance pattern, indicating that other processes affect the differences among populations, such as territoriality and site fidelity.     

德宏水牛微卫星标记分析的群体遗传变异

德宏水牛是云南省地方水牛的优良品种之一,为了进一步阐明其群体遗传变异和遗传结构,筛选了分别位于水牛14条染色体上的15对微卫星引物,对德宏水牛81个个体进行了检测分析.共检测到62个等位基因,每个座位等位基因数目从2到6个不等,平均等位基因数为4.13,该水牛群体期望杂合度和多态信息含量分别为0.6520±0.1526和0.5863±0.1789,各座位的遗传分化系数在0～0.0919之间,平均值为0.0202.每个座位的基因流较大,平均12.1502.研究结果表明德宏水牛群体遗传多样性较丰富,亚群间的遗传分化程度低,基因流较大,且很少发生近交.     

Isolation and characterization of eight microsatellite loci in the brook lamprey Lampetra planeri (Petromyzontiformes) using 454 sequence data

Eight polymorphic microsatellite loci were developed for the brook lamprey Lampetra planeri through 454 sequencing and their usefulness was tested in 45 individuals of both L. planeri and the river lamprey Lampetra fluviatilis. The number of alleles per loci ranged between two and five; the Italian and Irish populations had a mean expected heterozygosity of 0·388 and 0·424 and a mean observed heterozygosity of 0·418 and 0·411, respectively.     

A set of polymorphic DNA microsatellites useful in swamp and river buffalo (Bubalus bubalis)

DNA microsatellites have found widespread application in gene mapping, pedigree determination and population genetics. In closely related species such as bovids, heterologous polymerase chain reaction (PCR) primers may in some cases be used, bypassing the need to isolate and characterize microsatellite-containing sequences and design PCR primers. We report on the ability of a set of eighty bovine derived DNA microsatellite primers to amplify sequences in the two types (swamp and river) of water buffalo ( Bubalus bubalis ). Number of alleles and per cent heterozygosities in a large number of animals were determined on a subset of microsatellite loci selected on the robustness of the primers. These loci will form the basis of a set of polymorphic DNA markers for use in water buffalo.     

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.     

Mitochondrial DNA analyses of Indian water buffalo support a distinct genetic origin of river and swamp buffalo

Water buffalo (Bubalus bubalis) is broadly classified into river and swamp categories, but it remains disputed whether these two types were independently domesticated, or if they are the result of a single domestication event. In this study, we sequenced the mitochondrial D-loop region and cytochrome b gene of 217 and 80 buffalo respectively from eight breeds/locations in northern, north-western, central and southern India and compared our results with published Mediterranean and swamp buffalo sequences. Using these data, river and swamp buffalo were distinguished into two distinct clades. Based upon the existing knowledge of cytogenetic, ecological and phenotypic parameters, molecular data and present-day distribution of the river and swamp buffalo, we suggest that these two types were domesticated independently, and that classification of the river and swamp buffalo as two related subspecies is more appropriate.     

Development and characterization of 11 novel microsatellite loci for the roundscale spearfish Tetrapturus georgii and their cross‐species amplification among other istiophorid species

Eleven novel polymorphic microsatellite loci were developed and characterized for the recently validated roundscale spearfish Tetrapturus georgii. Characterization of these markers, based on 35 roundscale spearfish from the western North Atlantic, revealed two to 21 alleles per locus with an average expected heterozygosity (H_E) of 0·09–0·94, and all loci conformed to Hardy–Weinberg expectations. Cross‐amplification of these 11 loci against all other eight known istiophorid species indicates promising prospects for the utility of these markers for istiophorids in general.     

Mitochondrial sequence-based evolutionary analysis of riverine–swamp hybrid buffaloes of India indicates novel maternal differentiation and domestication patterns

In this study, mitochondrial D-loop sequence data on riverine, swamp and hybrid buffaloes from India have been generated and compared with other reported Indian riverine, Chinese and Bangladeshi swamp buffalo populations. Sequence analysis revealed the presence of 132 haplotypes, with a haplotype diversity of 0.9611 ± 0.0045 and a nucleotide diversity of 0.04801 ± 0.00126. For the first time, the existence of riverine–swamp hybrids among the Indian Chilika buffalo population has been recorded, having 49 chromosomes, which was also confirmed by mitochondrial haplotype sharing between Chilika and Indian swamp as well as Chinese swamp buffalo populations in the network analysis. Phylogenetic analysis documents the sharing of reported pre-domestication haplogroups ‘SA1’, ‘SA2’, ‘SA3’ and ‘SB1’ between the Chilika and swamp buffalo populations of India, China and Bangladesh, an indication of the migration of swamp buffaloes towards Bangladesh and adjoining lower parts of India and north towards Chinese domestication sites. The results have also been supplemented by multidimension scaling, grouping Indian and Chinese swamp buffaloes more closely together with Bangladeshi buffaloes, but into a separate quadrant, whereas Chilika grouped away from other riverine as well as swamp buffaloes. These findings thus confirm the previous reports that the northeast region of India, close to the Indo-China border, is the point of evolution of swamp buffaloes with multiple sites of domestication.     

Genomic differentiation between swamp and river buffalo using a cattle high-density single nucleotide polymorphisms panel

Buffalo (Bubalus bubalis) is an important livestock species in many tropical and subtropical regions. In recent decades, the interest in buffalo’s milk have expanded and intensive buffalo farms start to emerge. However, breeding programs and population genetics information for this species is scarce or inexistent. The present study aims to test the suitability of the commercial high-density single nucleotide polymorphisms (SNP) genotyping panel, the Illumina BovineHD BeadChip, to estimate population genetics parameters, pedigree control and identification of common variants in major production candidate genes. From a total of 777 962 SNPs included in the panel, 20 479 were polymorphic in water buffalo at a call rate of 86% and an average expected heterozygosity (H_E) of 0.306. From these, 357 were mapped within or around the flanking regions of several major candidate genes. A principal components analysis identified three different clusters, each representing pure swamp buffalo type, pure river buffalo type and admixed river buffalo. The hybrids between swamp and river buffalo were clearly identified as an intermediary cluster. The suitability of these SNPs data set for parentage and identity testing demonstrated that the combination of just 30 to 50 SNPs were enough to attain high probabilities of parentage exclusion (0.9999) in both types and identity (2.3×10⁻⁵ and 2.0×10⁻⁷) for river and swamp buffalo, respectively. Our analysis confirms the suitability of the BovineHD BeadChip to assess population structure, hybridization and identity of the water buffalo populations.     

Genetic diversity of mitochondrial cytochrome b gene in Chinese native buffalo

The origins of the domestic water buffalo remain contentious. To better understand the origins of Chinese water buffalo, we sequenced the complete mitochondrial cytochrome b (MT-CYB) gene from 270 individuals representing 13 Chinese domestic swamp buffalo populations. We found genetic evidence of introgression of river buffalo into Chinese swamp buffalo herds. Swamp buffalo haplotypes can be divided into two highly divergent lineages (A and B), suggesting that Chinese native swamp buffalo have two maternal origins. We found that the A→G transition in the buffalo MT-CYB gene stop codon resulted in buffalo haplotypes being terminated by one of two stop codons: AGA or AGG. AGA is common to river buffalo and lineage A of swamp buffalo, while AGG is specific to lineage B of swamp buffalo. Lineage A appears to have been domesticated in China. Further genetic evidence is required to clarify the origins of lineage B.     

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.     

Strong and stable geographic differentiation of swamp buffalo maternal and paternal lineages indicates domestication in the China/Indochina border region

The swamp type of the Asian water buffalo is assumed to have been domesticated by about 4000 years BP, following the introduction of rice cultivation. Previous localizations of the domestication site were based on mitochondrial DNA (mtDNA) variation within China, accounting only for the maternal lineage. We carried out a comprehensive sampling of China, Taiwan, Vietnam, Laos, Thailand, Nepal and Bangladesh and sequenced the mtDNA Cytochrome b gene and control region and the Y‐chromosomal ZFY, SRY and DBY sequences. Swamp buffalo has a higher diversity of both maternal and paternal lineages than river buffalo, with also a remarkable contrast between a weak phylogeographic structure of river buffalo and a strong geographic differentiation of swamp buffalo. The highest diversity of the swamp buffalo maternal lineages was found in south China and north Indochina on both banks of the Mekong River, while the highest diversity in paternal lineages was in the China/Indochina border region. We propose that domestication in this region was later followed by introgressive capture of wild cows west of the Mekong. Migration to the north followed the Yangtze valley as well as a more eastern route, but also involved translocations of both cows and bulls over large distances with a minor influence of river buffaloes in recent decades. Bayesian analyses of various migration models also supported domestication in the China/Indochina border region. Coalescence analysis yielded consistent estimates for the expansion of the major swamp buffalo haplogroups with a credibility interval of 900 to 3900 years BP. The spatial differentiation of mtDNA and Y‐chromosomal haplotype distributions indicates a lack of gene flow between established populations that is unprecedented in livestock.     

Isolation and characterization of 13 polymorphic microsatellite loci in Hypostomus ancistroides (Teleostei,Loricariidae) and cross‐amplification in related species

This study details 13 novel polymorphic microsatellite loci in the armoured catfish Hypostomus ancistroides, and assesses their utility for population genetic studies. The analysis of 30 individuals revealed a total of 99 different alleles (ranging from two to 15 alleles per locus), with an average of 7·62 alleles per locus, with observed and expected heterozygosities ranging from 0·103 to 0·931 and from 0·102 to 0·906, respectively. One of the 13 loci showed significant deviations from Hardy–Weinberg equilibrium, probably due to the presence of null alleles, inferred from the excess of homozygotes.