Analysis of porcine body size variation using re-sequencing data of miniature and large pigs.

BACKGROUND: Domestication has led to substantial phenotypic and genetic variation in domestic animals. In pigs, the size of so called minipigs differs by one order of magnitude compared to breeds of large body size. We used biallelic SNPs identified from re-sequencing data to compare various publicly available wild and domestic populations against two minipig breeds to gain better understanding of the genetic background of the extensive body size variation. We combined two complementary measures, expected heterozygosity and the composite likelihood ratio test implemented in "SweepFinder", to identify signatures of selection in Minipigs. We intersected these sweep regions with a measure of differentiation, namely FST, to remove regions of low variation across pigs. An extraordinary large sweep between 52 and 61 Mb on chromosome X was separately analyzed based on SNP-array data of F2 individuals from a cross of Goettingen Minipigs and large pigs. RESULTS: Selective sweep analysis identified putative sweep regions for growth and subsequent gene annotation provided a comprehensive set of putative candidate genes. A long swept haplotype on chromosome X, descending from the Goettingen Minipig founders was associated with a reduction of adult body length by 3% in F2 cross-breds. CONCLUSION: The resulting set of genes in putative sweep regions implies that the genetic background of body size variation in pigs is polygenic rather than mono- or oligogenic. Identified genes suggest alterations in metabolic functions and a possible insulin resistance to contribute to miniaturization. A size QTL located within the sweep on chromosome X, with an estimated effect of 3% on body length, is comparable to the largest known in pigs or other species. The androgen receptor AR, previously known to influence pig performance and carcass traits, is the most obvious potential candidate gene within this region.

Reducing the interval of a growth QTL on chromosome 4 in laying hens.

In our previous research, we identified a QTL with an interval of 3.4 Mb for growth on chicken chromosome (GGA) 4 in an advanced intercross population of an initial cross between the New Hampshire inbred line (NHI) and the White Leghorn inbred line (WL77). In the current study, an association analysis was performed in a population of purebred white layers (WLA) with White Leghorn origin. Genotypic data of 130 SNPs within the previously identified 3.4-Mb region were obtained using a 60K SNP chip. In total, 24 significant SNPs (LOD ≥ 4.44) on GGA4 were detected for daily weigh gain from 8 to 14 weeks and two SNPs (LOD ≥ 4.80) for body weight at 14 weeks. The QTL interval was reduced by 1.9 Mb to an interval of 1.5 Mb (74.6-76.1 Mb) that harbors 15 genes. Furthermore, to identify additional loci for chicken growth, a genome-wide association study (GWAS) was carried out in a WLA population. The GWAS identified an additional QTL on GGA6 for body weight at six weeks (19.8-21.2 Mb). Our findings showed that by using a WLA population we were able to further reduce the QTL confidence interval previously detected using a NHI × WL77 advanced intercross population.

Ghrelin plasma concentration does not covary with energy demand in adult laying hens.

The peptide hormone ghrelin is suggested to be involved in food intake regulation in young growing chicken. Whether ghrelin is involved in the regulation of energetic balance associated with laying performance in adult laying hens was studied by use of 4 chicken lines that differ in laying performance and phylogeny (4 lines; 16 hens per line). As housing conditions are also known to affect energy demand, half of the hens per line were housed in single cages and the other half of hens were maintained in a floor housing system. Plasma samples were collected at 17 to 19, 33 to 35, 49 to 51, and 72 wk of age and analyzed with a chicken ghrelin ELISA Kit. From caged hens, individual food consumption and laying performance additionally was recorded. Due to its function in growth and its relationship with ghrelin, also GH plasma concentrations were analyzed. Ghrelin concentrations did not differ between the 4 lines at any of the test periods (all P > 0.05). Ghrelin was negatively related to food consumption only in the growing period of the high-performing lines (both P < 0.0001). During this phase, floor-housed hens showed greater ghrelin concentrations compared with caged hens (P < 0.0001). Our results suggest that in adult layers ghrelin is not involved in regulating energy intake related to laying performance but rather seems to be related to body growth and housing condition before start of lay, the latter possibly due to differences in hens' behavioral activity.

Chicken line-dependent mortality after experimental infection with three type IIxIII recombinant Toxoplasma gondii clones.

Three genetically different clones of Toxoplasma gondii, also different in mouse virulence, were studied by experimental infection in chickens. For the experiments, four chicken lines were used, which differed in phylogenetic origin and performance level: two white egg layer lines, one with high laying performance (WLA), one with low (R11) and two brown layer lines, also displaying high (BLA) and low (L68) egg number. Chickens were intraperitoneally infected with three different T. gondii isolates representing type IIxIII recombinant clones, i.e. showing both, type II- and type III-specific alleles. These clones (K119/2 2C10, B136/1 B6H6, K119/2 A7) had exhibited virulence differences in a mouse model. In chickens, a significantly higher mortality was observed in white layer lines, but not in brown layer lines, suggesting that differences in the phylogenetic background may influence the susceptibility of chickens for toxoplasmosis. In addition, antibody (IgY) levels varied in surviving chickens at 31 days post infection. While low to intermediate antibody levels were observed in white layers, intermediate to high levels were measured in brown layers. Infection with a T. gondii clone showing low chicken virulence resulted in higher antibody levels in all chicken lines compared to infection with T. gondii clones of intermediate or high chicken virulence. This was in agreement with the parasite load as determined by real-time PCR. Overall, results show that progeny resulting from natural sexual recombination of T. gondii clonal lineages, may differ in their virulence for mice and chickens.

Haematological and febrile response to Escherichia coli lipopolysaccharide in 12-week-old cockerels of genetically diverse layer lines fed diets with increasing L-arginine levels.

Due to its decisive function in the avian metabolic, endocrine and immune system L-arginine (Arg) is dietary indispensable for chickens. In 12-week-old cockerels of two high- and two low-performing purebred layer lines, the effects of increasing dietary Arg on the haematological and febrile response were studied over 48 h after single lipopolysaccharide (LPS) injection. The offered diets contained Arg equivalent to 70%, 100% and 200% of recommended supply. Pathophysiological alterations in weight gain, feed intake, body temperature and differential blood count were examined in comparison with their physiological initial values. Within the first 24 h after LPS injection, cockerels reduced feed intake and lost body weight subsequently. Thereby, low-performing genotypes lost body weight to a lesser extent than high-performing ones. The loss of body weight was further intensified by deficient dietary Arg. Within the following 24 h, cockerels recovered by improving feed intake and weight gain. Furthermore, LPS induced genotype-specific fever response: both brown genotypes showed initial hypothermia followed by longer lasting moderate hyperthermia, whereas the white genotypes exhibited biphasic hyperthermia. Fever response was accompanied by significant changes in differential blood counts. Characterized by lymphopenia and heterophilia, a severe leucopenia was observed from 4 to 8 h after LPS injection and replaced by a marked leucocytosis with longer lasting monocytosis up to 48 h after LPS injection. Under given pathophysiological conditions, deficiently Arg-supplied cockerels showed higher total leucocyte counts than adequately and excessively Arg-supplied cockerels. However, deficient and surplus dietary Arg tended to cause higher ratios between heterophils and lymphocytes. To conclude, present results confirmed that LPS induced numerous immunological changes in 12-week-old cockerels and emphasized that chicken's genotype is a source of variation to be considered for immunological studies. Deficient dietary Arg intensified acute changes in differential blood counts and weight gain during LPS-induced inflammation.

Metabolic and clinical response to Escherichia coli lipopolysaccharide in layer pullets of different genetic backgrounds supplied with graded dietary L-arginine.

L-arginine (Arg) is an essential amino acid in birds that plays a decisive role in avian protein synthesis and immune response. Effects of graded dietary Arg supply on metabolic and clinical response to Escherichia coli lipopolysaccharide (LPS) were studied over 48 hours after a single intramuscular LPS injection in 18-week-old genetically diverse purebred pullets. LPS induced a genotype-specific fever response within 4 hours post injectionem. Whereas brown genotypes showed an initial hypothermia followed by longer-lasting moderate hyperthermia, white genotypes exhibited a biphasic hyperthermia without initial hypothermia. Furthermore, within 2 hours after LPS injection, sickness behavior characterized by lethargy, anorexia, intensified respiration, and ruffled feathers appeared, persisted for 3 to 5 hours and recovered 12 hours post injectionem. The varying grades of Arg did not alter the examined traits named above, whereas insufficient Arg reduced body growth and increased relative weights of liver and pancreas significantly. At 48 hours post injectionem, increased relative weights of liver and spleen were also found in LPS treated pullets, whereas LPS decreased those of pancreas, bursa, thymus, and cecal tonsils. Moreover, LPS lowered the sum of plasma amino acids and decreased plasma concentrations of Arg, citrulline, glutamate, methionine, ornithine, phenylalanine, proline, tryptophan, and tyrosine, and increased those of aspartate, glutamine, lysine, 1- and 3-methyl-histidine. Elevating concentrations of dietary Arg led to increasing plasma concentrations of Arg, citrulline, ornithine, and 3-methyl-histidine subsequently. As quantitative expression of LPS-induced anorexia, proteolysis, and the following changes in plasma amino acids, pullets showed a significant decrease of feed and nitrogen intake and catabolic metabolism characterized by negative nitrogen balance and body weight loss in the first 24 hours post injectionem. Pullets recovered from the challenge within the second 24 hours post injectionem and changed to anabolism with re-increased feed and nitrogen intake, positive nitrogen retention, and weight gain. To conclude, present results confirmed that LPS induced numerous metabolic and physiological changes in pullet's genotypes, whereas dietary Arg affected the examined traits only slightly.

Using the variability of linkage disequilibrium between subpopulations to infer sweeps and epistatic selection in a diverse panel of chickens.

A whole-genome scan for identifying selection acting on pairs of linked loci is proposed and implemented. The scan is based on , one of Ohta's 1982 measures of between-population linkage disequilibrium (LD). An approximate empirical null distribution for the statistic is suggested. Although the partitioning of LD into between-population components was originally used to investigate epistatic selection, we demonstrate that values of may also be influenced by single-locus selective sweeps with linkage but no epistasis. The proposed scan is implemented in a diverse panel of chickens including 72 distinct breeds genotyped at 538 298 single-nucleotide polymorphisms. In all, 1723 locus pairs are identified as putatively corresponding to a selective sweep or epistatic selection. These pairs of loci generally cluster to form overlapping or neighboring signals of selection. Known variants that were expected to have been under selection in the panel are identified, as well as an assortment of novel regions that have putatively been under selection in chickens. Notably, a promising pair of genes located 8 MB apart on chromosome 9 are identified based on as demonstrating strong evidence of dispersive epistatic selection between populations.

Maternal genealogical patterns of chicken breeds sampled in Europe.

The aim of this study was to investigate the maternal genealogical pattern of chicken breeds sampled in Europe. Sequence polymorphisms of 1256 chickens of the hypervariable region (D-loop) of mitochondrial DNA (mtDNA) were used. Median-joining networks were constructed to establish evolutionary relationships among mtDNA haplotypes of chickens, which included a wide range of breeds with different origin and history. Chicken breeds which have had their roots in Europe for more than 3000 years were categorized by their founding regions, encompassing Mediterranean type, East European type and Northwest European type. Breeds which were introduced to Europe from Asia since the mid-19th century were classified as Asian type, and breeds based on crossbreeding between Asian breeds and European breeds were classified as Intermediate type. The last group, Game birds, included fighting birds from Asia. The classification of mtDNA haplotypes was based on Liu et al.'s (2006) nomenclature. Haplogroup E was the predominant clade among the European chicken breeds. The results showed, on average, the highest number of haplotypes, highest haplotype diversity, and highest nucleotide diversity for Asian type breeds, followed by Intermediate type chickens. East European and Northwest European breeds had lower haplotype and nucleotide diversity compared to Mediterranean, Intermediate, Game and Asian type breeds. Results of our study support earlier findings that chicken breeds sampled in Europe have their roots in the Indian subcontinent and East Asia. This is consistent with historical and archaeological evidence of chicken migration routes to Europe.

Properties of different selection signature statistics and a new strategy for combining them.

Identifying signatures of recent or ongoing selection is of high relevance in livestock population genomics. From a statistical perspective, determining a proper testing procedure and combining various test statistics is challenging. On the basis of extensive simulations in this study, we discuss the statistical properties of eight different established selection signature statistics. In the considered scenario, we show that a reasonable power to detect selection signatures is achieved with high marker density (>1 SNP/kb) as obtained from sequencing, while rather small sample sizes (~15 diploid individuals) appear to be sufficient. Most selection signature statistics such as composite likelihood ratio and cross population extended haplotype homozogysity have the highest power when fixation of the selected allele is reached, while integrated haplotype score has the highest power when selection is ongoing. We suggest a novel strategy, called de-correlated composite of multiple signals (DCMS) to combine different statistics for detecting selection signatures while accounting for the correlation between the different selection signature statistics. When examined with simulated data, DCMS consistently has a higher power than most of the single statistics and shows a reliable positional resolution. We illustrate the new statistic to the established selective sweep around the lactase gene in human HapMap data providing further evidence of the reliability of this new statistic. Then, we apply it to scan selection signatures in two chicken samples with diverse skin color. Our analysis suggests that a set of well-known genes such as BCO2, MC1R, ASIP and TYR were involved in the divergent selection for this trait.

Global diversity and genetic contributions of chicken populations from African, Asian and European regions.

Genetic diversity and population structure of 113 chicken populations from Africa, Asia and Europe were studied using 29 microsatellite markers. Among these, three populations of wild chickens and nine commercial purebreds were used as reference populations for comparison. Compared to commercial lines and chickens sampled from the European region, high mean numbers of alleles and a high degree of heterozygosity were found in Asian and African chickens as well as in Red Junglefowl. Population differentiation (FST ) was higher among European breeds and commercial lines than among African, Asian and Red Junglefowl populations. Neighbour-Net genetic clustering and structure analysis revealed two main groups of Asian and north-west European breeds, whereas African populations overlap with other breeds from Eastern Europe and the Mediterranean region. Broilers and brown egg layers were situated between the Asian and north-west European clusters. structure analysis confirmed a lower degree of population stratification in African and Asian chickens than in European breeds. High genetic differentiation and low genetic contributions to global diversity have been observed for single European breeds. Populations with low genetic variability have also shown a low genetic contribution to a core set of diversity in attaining maximum genetic variation present from the total populations. This may indicate that conservation measures in Europe should pay special attention to preserving as many single chicken breeds as possible to maintain maximum genetic diversity given that higher genetic variations come from differentiation between breeds.

Phylogenetic resolution power of microsatellites and various single-nucleotide polymorphism types assessed in 10 divergent chicken populations.

There has been some debate over the question of which types of DNA variation are most appropriate to accurately reconstruct evolutionary events. We compared the capacity of microsatellites (STRs) and various types of single-nucleotide polymorphism (SNP) loci in the chicken genome. The SNP types differ in their location: in exons, introns and promoters. Genetic distances between all possible pairs of 10 populations were calculated for each marker type. STR loci, which are much more polymorphic than are SNPs, are considered to have occurred at recent time compared with old evolutionary events of SNPs. Using structure software, STR loci assigned individuals to their population much more correctly than did any other marker types, whereas SNPs at promoter regions gave the poorest ascription. Furthermore, 29 STR markers were even better than all 152 SNPs together. Ancient evolutionary events that produced genetic differences between the most distant populations such as Red Jungle Fowl and domestic chicken were detected better by exons and introns than by STR loci and promoters. The significant interactions found between marker types and populations suggest that marker types had different phylogenetic histories, possibly related to a different timescale.

Genetic diversity and population structure of five Ethiopian chicken ecotypes.

This study aimed to analyse the genetic diversity and population structure of five Ethiopian chicken ecotypes (N = 155), which were compared with six commercial purebreds (N = 180). For the analysis of genetic diversity, 26 AVIANDIV microsatellite markers were used. The number of alleles in Ethiopian ecotypes ranged from 2 to 19 per locus, with a mean of 6.1. The average observed heterozygosity within ecotype varied between 0.53 and 0.57. The overall heterozygote deficiency (F(IT)) in Ethiopian ecotypes was 0.124 ± 0.037. Over 68% of F(IT) was because of within-ecotype deficiency (F(IS)). In the phylogenetic tree, Ethiopian ecotypes clustered into two groups. The analysis of the relationship between populations using the structure program provided further evidence for the occurrence of at least two subgroups in the Ethiopian ecotypes. Findings of this study may provide the background for future studies to identify the origin of the two gene pools representing the Ethiopian chicken ecotypes and to characterize the gene variants influencing economically important traits.

Molecular tools and analytical approaches for the characterization of farm animal genetic diversity.

Genetic studies of livestock populations focus on questions of domestication, within- and among-breed diversity, breed history and adaptive variation. In this review, we describe the use of different molecular markers and methods for data analysis used to address these questions. There is a clear trend towards the use of single nucleotide polymorphisms and whole-genome sequence information, the application of Bayesian or Approximate Bayesian analysis and the use of adaptive next to neutral diversity to support decisions on conservation.

Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations.

Many studies in human genetics compare informativeness of single-nucleotide polymorphisms (SNPs) and microsatellites (single sequence repeats; SSR) in genome scans, but it is difficult to transfer the results directly to livestock because of different population structures. The aim of this study was to determine the number of SNPs needed to obtain the same differentiation power as with a given standard set of microsatellites. Eight chicken breeds were genotyped for 29 SSRs and 9216 SNPs. After filtering, only 2931 SNPs remained. The differentiation power was evaluated using two methods: partitioning of the Euclidean distance matrix based on a principal component analysis (PCA) and a Bayesian model-based clustering approach. Generally, with PCA-based partitioning, 70 SNPs provide a comparable resolution to 29 SSRs. In model-based clustering, the similarity coefficient showed significantly higher values between repeated runs for SNPs compared to SSRs. For the membership coefficients, reflecting the proportion to which a fraction segment of the genome belongs to the ith cluster, the highest values were obtained for 29 SSRs and 100 SNPs respectively. With a low number of loci (29 SSRs or ≤100 SNPs), neither marker types could detect the admixture in the Gödöllö Nhx population. Using more than 250 SNPs allowed a more detailed insight into the genetic architecture. Thus, the admixed population could be detected. It is concluded that breed differentiation studies will substantially gain power even with moderate numbers of SNPs.

Genetic diversity of ten Egyptian chicken strains using 29 microsatellite markers.

In this study, we assessed the genetic diversity of three Egyptian local chicken strains (Fayoumi, Dandarawi and Sinai) and six synthetic breeds derived from Fayoumi and Sinai by intercrossing with Barren Plymouth Rock, Rhode Island Red or White Cornish. Diversity measures were based on interrogation of 29 microsatellites. We identified three main clusters of chicken populations encompassing selected Fayoumi lines and Doki-4 (cluster-1), native Dandarawi (cluster-2) and Sinai, and all six synthetic breeds (cluster-3). Dandarawi and Fayoumi lines exhibited lower intra-population genetic diversity and allelic privacy than Sinai and synthetic breeds. The global inbreeding (F(IT) ) was 0.11, among-population differentiation (F(ST) ) was 0.07, and within-population differentiation (F(IS) ) was 0.04. The between-population marker-estimated kinship was lower than within-population estimates. The cluster analysis classified the Fayoumi lines, Dandarawi and Gimmizah as clearly separated populations. The other strains were configured in mosaic admixed groups.

Diversity and origin of South African chickens.

The objectives of this study were to analyze the genetic diversity and structure of South African conserved and field chicken populations and to investigate the maternal lineages of these chicken populations. Four South African conserved chicken populations (n = 89), namely, Venda (VD_C), Ovambo, Naked Neck, and Potchefstroom Koekoek from the Animal Production Institute of the Agricultural Research Council, and 2 field populations, the Venda and Ovambo (OV_F), from which the Ovambo and the Venda conservation flocks were assumed to have been sampled, were genotyped for 460 bp of the mitochondrial DNA (mtDNA) D-loop sequence. Haplotypes of these chickens were aligned to 7 Japanese and 9 Chinese and Eurasian chicken mtDNA D-loop sequences taken from GenBank and reflecting populations from presumed centers of domestication. Sequence analysis revealed 48 polymorphic sites that defined 13 haplotypes in the South African chicken populations. All 6 South African conserved and field chicken populations observed were found to be polymorphic, with the number of haplotypes ranging from 3 for VD_C to 8 for OV_F. The lowest haplotype diversity, 0.54 ± 0.08, was observed in VD_C chickens, whereas the highest value, 0.88 ± 0.05, was observed in OV_F chickens. Genetic diversity between the 4 South African conserved and 2 field chicken populations constituted 12.34% of the total genetic variation, whereas within-population diversity constituted 87.66% of the total variation. The median network analysis of the mtDNA D-loop haplotypes observed in the South African conserved and field populations and the reference set resulted in 5 main clades. All 6 South African chickens were equally represented in the major clade, E, which is presumed to be of Indian subcontinent maternal origin and may have its roots in Southeast Asia. The results showed multiple maternal lineages of South African chickens. Conservation flocks and field chicken populations shared the major haplotypes A, D and E, which were presumed to be of Chinese, Southeast Asian, and Indian subcontinental origin.

Conservation priorities and optimum allocation of conservation funds for Vietnamese local chicken breeds.

The objectives of this study were to estimate conservation potential of Vietnamese local breeds and to investigate optimal allocation of conservation funds to maximize genetic diversity conserved between these breeds. Twenty-nine microsatellites were genotyped in 353 individuals from nine Vietnamese local chicken breeds and two chicken breeds of Chinese origin. The Vietnamese chicken breeds were sampled from the northern and southern parts of Vietnam while the two Chinese breeds have been kept as conservation flocks at the National Institute of Animal Sciences, Hanoi. The Weitzman approach was used to assess alternative strategies for conserving genetic diversity between breeds. Three different models, which reflect the range of possible functions in typical conservation situations, were applied. An average extinction probability of 48.5% was estimated for all Vietnamese chicken breeds. The highest conservation potential was found in the Te, Dong Tao and Ac chicken breeds, whereas the lowest corresponding values were observed in the Ri and Mia chicken breeds. The conservation funds were mainly allocated to the same three breeds (Te, Dong Tao and Ac) under all three models. This study suggests that conservation potential of the Vietnamese chicken breeds varies considerably. Population priorities for allocation of conservation funds in this study do not depend on the cost model used. The three breeds (Te, Dong Tao and Ac) with the highest conservation potential should be the prime candidates to be allocated conservation funds if the conservation budgets are limited.

Genetic diversity and conservation of South African indigenous chicken populations.

In this study, we compare the level and distribution of genetic variation between South African conserved and village chicken populations using microsatellite markers. In addition, diversity in South African chickens was compared to that of a reference data set consisting of other African and purebred commercial lines. Three chicken populations Venda, Ovambo and Eastern Cape and four conserved flocks of the Venda, Ovambo, Naked Neck and Potchefstroom Koekoek from the Poultry Breeding Resource Unit of the Agricultural Research Council were genotyped at 29 autosomal microsatellite loci. All markers were polymorphic. Village chicken populations were more diverse than conservation flocks. structure software was used to cluster individuals to a predefined number of 2 ≤ K ≤ 6 clusters. The most probable clustering was found at K = 5 (95% identical runs). At this level of differentiation, the four conservation flocks separated as four independent clusters, while the three village chicken populations together formed another cluster. Thus, cluster analysis indicated a clear subdivision of each of the conservation flocks that were different from the three village chicken populations. The contribution of each South African chicken populations to the total diversity of the chickens studied was determined by calculating the optimal core set contributions based on Marker estimated kinship. Safe set analysis was carried out using bootstrapped kinship values calculated to relate the added genetic diversity of seven South African chicken populations to a set of reference populations consisting of other African and purebred commercial broiler and layer chickens. In both core set and the safe set analyses, village chicken populations scored slightly higher to the reference set compared to conservation flocks. Overall, the present study demonstrated that the conservation flocks of South African chickens displayed considerable genetic variability that is different from that of the assumed founder populations (village chickens).

Genetic diversity in farm animals--a review.

Domestication of livestock species and a long history of migrations, selection and adaptation have created an enormous variety of breeds. Conservation of these genetic resources relies on demographic characterization, recording of production environments and effective data management. In addition, molecular genetic studies allow a comparison of genetic diversity within and across breeds and a reconstruction of the history of breeds and ancestral populations. This has been summarized for cattle, yak, water buffalo, sheep, goats, camelids, pigs, horses, and chickens. Further progress is expected to benefit from advances in molecular technology.

Assessing genetic diversity of Vietnamese local chicken breeds using microsatellites.

This study aimed to assess genetic diversity within and between nine Vietnamese local chicken breeds and two Chinese breeds included for comparison. Genotyping 29 microsatellites revealed high diversity of both Vietnamese and Chinese breeds. Cluster analysis using the STRUCTURE software suggested six clusters as the most likely grouping of the 11 breeds studied. These groups encompassed four homogeneous clusters, one formed by the two Chinese breeds and the other three representing a single breed each: the Mekong Delta breed Ac, the South Central Coast breed Choi, and the Red River Delta breed Dong Tao. The six remaining breeds formed two additional admixed clusters.