Analysis of genetic diversity in Bolivian llama populations using microsatellites

South American camelids (SACs) have a major role in the maintenance and potential future of rural Andean human populations. More than 60% of the 3.7 million llamas living worldwide are found in Bolivia. Due to the lack of studies focusing on genetic diversity in Bolivian llamas, this analysis investigates both the genetic diversity and structure of 12 regional groups of llamas that span the greater part of the range of distribution for this species in Bolivia. The analysis of 42 microsatellite markers in the considered regional groups showed that, in general, there were high levels of polymorphism (a total of 506 detected alleles; average PIC across per marker: 0.66), which are comparable with those reported for other populations of domestic SACs. The estimated diversity parameters indicated that there was high intrapopulational genetic variation (average number of alleles and average expected heterozygosity per marker: 12.04 and 0.68, respectively) and weak genetic differentiation among populations (F_ST range: 0.003–0.052). In agreement with these estimates, Bolivian llamas showed a weak genetic structure and an intense gene flow between all the studied regional groups, which is due to the exchange of reproductive males between the different flocks. Interestingly, the groups for which the largest pairwise F_ST estimates were observed, Sud Lípez and Nor Lípez, showed a certain level of genetic differentiation that is probably due to the pattern of geographic isolation and limited communication infrastructures of these southern localities. Overall, the population parameters reported here may serve as a reference when establishing conservation policies that address Bolivian llama populations.     

Past climate change and recent anthropogenic activities affect genetic structure and population demography of the greater long‐tailed hamster in northern China

The genetic diversity and the spatial structure of a species are likely consequences of both past and recent evolutionary processes, but relevant studies are still rare in East Asia where the Pleistocene climate has unique influences. In this study, we examined the impact of past climate change and recent anthropogenic activities on the genetic structure and population size of the greater long‐tailed hamster (Tscherskia triton), an agricultural rodent pest species in northern China. DNA sequence data of 2 mitochondrial genes and genotypic data of 11 microsatellite DNA loci from 41 populations (545 individuals) were gathered. Phylogenetic and population genetic analyses, as well as species distribution modeling and coalescent simulations, were conducted to infer its historical and demographic patterns and processes. Two deeply diverged mitochondrial clades were recovered. A small one was restricted to the Shandong Peninsula while the main clade was further divided into 3 geographic clusters by their microsatellite DNA genotypes: Northwest, North‐center and Northeast. Divergence dating indicated a Middle‐to‐Late Pleistocene divergence between the 2 clades. Demographic analysis indicated that all 3 and pooled populations showed consistent long‐period expansions during last glacial period; but not during the Holocene, probably due to the impact of climate warming and human disturbances. Conflicting patterns between mtDNA and microsatellite markers imply an anthropogenic impact on North‐center populations due to intensified agricultural cultivation in this region. Our study demonstrated that the impact of past glaciation on organisms in East Asia significantly differs from that of Europe and North America, and human activity is an important factor in determining the genetic diversity of a species, as well as its spatial structure.     

Genetic diversity and differentiation of Mongolian and Russian yak populations

In this study we examined the genetic diversity of yak populations in the northernmost part of their current global distribution. Five Mongolian and one Russian yak populations as well as one Chinese yak population from the Qinghai–Tibetan Plateau, the putative centre of yak domestication, were analysed with 15 microsatellite loci to determine the level of genetic variation within populations as well as the genetic differentiation and relationship between populations. A total of 116 microsatellite alleles were identified. The mean number of alleles per locus (MNA) across populations was 7.73 ± 1.98 and the mean expected heterozygosity (H_E) was 0.696 ± 0.026. The relative magnitude of gene differentiation (F_ST) among populations was 4.1%, and all genetic differentiations (F_ST) between populations were significant (p < 0.001). A significant inbreeding effect (F_IS) was detected in the Hovsgol yak (p < 0.01). There was no indication of a recent bottleneck in any of the populations studied. The results showed that yak populations in Mongolia and Russia have maintained high genetic diversity within populations and a low, although significant, genetic differentiation between populations. Both phylogenetic and principal component analyses support a close genetic relationship between the Gobi Altai, south Gobi and north Hangai populations, and between the Hovsgol and Buryatia populations respectively. Our results indicate that these yak populations should be considered as distinct genetic entities in respect of conservation and breeding programmes.     

High genetic diversity in the endangered and narrowly distributed amphibian species Leptobrachium leishanense

Threatened species typically have a small or declining population size, which make them highly susceptible to loss of genetic diversity through genetic drift and inbreeding. Genetic diversity determines the evolutionary potential of a species; therefore, maintaining the genetic diversity of threatened species is essential for their conservation. In this study, we assessed the genetic diversity of the adaptive major histocompatibility complex (MHC) genes in an endangered and narrowly distributed amphibian species, Leptobrachium leishanense in Southwest China. We compared the genetic variation of MHC class I genes with that observed in neutral markers (5 microsatellite loci and cytochrome b gene) to elucidate the relative roles of genetic drift and natural selection in shaping the current MHC polymorphism in this species. We found a high level of genetic diversity in this population at both MHC and neutral markers compared with other threatened amphibian species. Historical positive selection was evident in the MHC class I genes. The higher allelic richness in MHC markers compared with that of microsatellite loci suggests that selection rather than genetic drift plays a prominent role in shaping the MHC variation pattern, as drift can affect all the genome in a similar way but selection directly targets MHC genes. Although demographic analysis revealed no recent bottleneck events in L. leishanense, additional population decline will accelerate the dangerous status for this species. We suggest that the conservation management of L. leishanense should concentrate on maximizing the retention of genetic diversity through preventing their continuous population decline. Protecting their living habitats and forbidding illegal hunting are the most important measures for conservation of L. leishanense.     

Assessment of genetic diversity and population structure of Vietnamese indigenous cattle populations by microsatellites

Cattle play a very important role in agriculture and food security in Vietnam. A high level of cattle diversity exists and serves different needs of Vietnamese cattle keepers but has not yet been molecularly characterized. This study evaluates the genetic diversity and structure of Vietnamese indigenous cattle populations, using microsatellite markers. A total of 410 individuals from six indigenous cattle populations and an exotic breed was characterized using 27 microsatellite markers A total of 362 alleles was detected and the number of alleles per locus ranged from 8 (INRA005 and ILSTS005) to 17 (ETH185). The level of gene diversity was high indicated by a mean expected heterozygosity (He) across populations and loci of 0.73. Level of inbreeding (mean F_IS=0.05) and genetic differentiation (mean F_ST=0.04) was moderate. The phylogenetic tree based on Reynolds genetic distance reflected geographic distances. Structure analysis indicated five homogeneous clusters. The Brahman, Lang Son, Ha Giang and U Dau Riu cattle were assigned to independent clusters while Nghe An, Thanh Hoa and Phu Yen cattle were grouped in a single cluster. We conclude that Vietnamese indigenous cattle have high levels of genetic diversity and distinct genetic structures. Based on these results, we recommend that for conservation homogenous populations (Nghe An, Thanh Hoa and Phu Yen) can be grouped to reduce costs and U Dau Riu, Lang Son and Ha Giang populations should be conserved separately to avoid loss of genetic diversity.     

Fluctuating fine-scale spatial genetic structure in the common shrews (Eulipotyphla,Mammalia)

The spatial genetic structure is a topical issue in the studies of various aspects of ecology and evolution. Using the multilocus autocorrelation method with hypervariable microsatellite genetic markers, we investigated a fine-scale pattern of genetic structure in 5 local populations of the common shrew Sorex araneus separated by distances of 300–1000 m (the Moscow chromosomal race). Spatial genetic autocorrelation analysis based on 5 microsatellite loci (expected heterozygosity >0.79) with 90 alleles revealed a consistent pattern of significant positive genetic structure. By testing the autocorrelation at multiple scales from 25 to 500 m, we found that positive spatial genetic structure is detectable in distance classes of <500 m. The weaker spatial genetic structure positively correlated with a higher ratio of nonresident individuals to residents’ activity (number of captures). In contrast to nonresident animals, the residents demonstrated prominent genetic structure. Genetic difference (F_ST) between the populations was significant (0.016–0.051) and comparable with that between populations of different races analyzed previously (0.016–0.038). F_ST was not associated with geographic distance. These demographic patterns allowed us to propose a scheme of genetic-structure dynamics involving periodic appearance of more related local groups and renewal of allelic profiles from а common pool where the alleles are mixed. The scheme predicts fluctuating genetic structure and random similar differences among local populations.     

Genetic structure and recent population demographic history of Taihangshan macaque (Macaca mulatta tcheliensis), North China

Massive actions have been and are being taken into protecting the world's primates from extinction, while the study of the properties of genetic diversity, demographic history, and ecological relationships will benefit the understanding of the long-term survival of a species. The Taihangshan macaque (Macaca mulatta tcheliensis), a subspecies of rhesus macaque (Macaca mulatta), is endemic to China and currently restricted to southern Mt. Taihangshan area. Herein, we evaluated the genetic diversity, population structure, and demographic history of this subspecies using mitochondrial (Cytb and high variable region I: HVR I) and nuclear markers (microsatellite loci) of 131 individuals collected from 9 localities covering the distribution range of this subspecies. Both phylogenetic analyses and genetic assignment revealed that the wild populations of Taihangshan macaques could be divided into 2 major highly divergent clades, THS-east and THS-west. Low genetic diversity (π: 0.00266 ± 0.00016) but high haplotype diversity (H_d: 0.80352 ± 0.015) were detected in the Taihangshan macaques, particularly in THS-east. Analyses of demographic history suggested that the Taihangshan macaques experienced first a stable historical population size from Holocene to early 19th century but a subtle decline and then slight growth in the recent 200 years. We suggest that bridging the neighbor populations (i.e. setting corridors) would facilitate the male-mediated gene flow and subsequently increase the genetic diversity of the Taihangshan macaque populations.     

Using genetic tools to inform conservation of fragmented populations of Asian elephants (Elephas maximus) across their range in China

A herd of 15 Chinese elephants attracted international attention during their 2021 northward trek, motivating the government to propose establishment of an Asian elephant national park. However, planning is hampered by a lack of genetic information on the remaining populations in China. We collected DNA from 497 dung samples from all 5 populations encompassing the entire range of elephants in China and used mitochondrial and microsatellite markers to investigate their genetic and demographic structure. We identified 237 unique genotypes (153 females, 84 males), representing 81% of the known population. However, the effective population size was small (28, range 25–32). Historic demographic contraction appeared to account for low haplotype diversity (H_d = 0.235), but moderate nucleotide and nuclear diversity (π = 0.6%, H_e = 0.55) was attributable to post-bottleneck recovery involving recent population expansion plus historical gene exchange with elephants in Myanmar, Lao PDR, and Vietnam. The 5 populations fell into 3 clusters, with Nangunhe elephants differing consistently from the other 4 populations (F_ST = 0.23); elephants from Mengyang, Simao, and Jiangcheng belonged to a single population (henceforth, MSJ), and differed from the Shangyong population (F_ST = 0.11). Interpopulation genetic variation reflected isolation by distance and female-biased dispersal. Chinese elephants should be managed as 2 distinct units: Nangunhe and another combining Shangyong and MSJ; their long-term viability will require restoring gene flow between Shangyong and MSJ, and between elephants in China and neighboring countries. Our results have the potential to inform conservation planning for an iconic megafaunal species.     

Analysis of genetic diversity and phylogenetic relationship of red deer subspecies in XinJiang,China

Polymorphisms for seven microsatellite loci in three red deer subspecies (9 populations) found in XinJiang were detected by polymerase chain reaction (PCR), 12% nondenaturation polyacrylamide gel electrophoresis and the Sanguinetti silver staining method. Numbers of alleles, average effective numbers of alleles (E) and the average rate of homozygosity, allelic frequencies of seven microsatellite loci, polymorphism information content (PIC), mean heterozygosity (H) and genetic distances among the populations were calculated for each population. Dendrograms were constructed based on genetic distances by the neighbor‐joining method (NJ), utilizing molecular evolutionary genetics analysis software PHYLIP (3.6). The phylogenetic tree was constructed based on allelic frequencies using maximum likelihood (ML); the bootstrap value was estimated by bootstrap test in the tree. Lastly, phylogenesis was analyzed. The results showed that four of the seven microsatellite loci were highly polymorphic, but BMS2508 and Celjp0023 showed no polymorphism and BM5004 was a neutral polymorphism. It is our conclusion that the four microsatellite loci are effective DNA markers for the analysis of genetic diversity and phylogenetic relationships among the three red deer subspecies. The mean PIC, H and E‐values across the microsatellite loci were 0.5393, 0.5736 and 2.64, which showed that these microsatellite loci are effective DNA markers for the genetic analysis of red deer. C.e. songaricus populations from Regiment 104, 151 and Hami are clustered together. C.e. yarkandensis populations from Regiment 35, Xaya and Alaer are clustered together. These two clusters also cluster together. Lastly, C.e. sibiricus populations from Burqin, Regiment 188 and the first two clusters were clustered together. The phylogenetic relationship among different red deer populations is consistent with the known origin, history of breeding and geographic distributions of populations.     

Mitochondrial DNA variation and genetic relationships in Spanish donkey breeds (Equus asinus)

In this study, the mitochondrial DNA diversity of six Spanish donkey breeds and two African donkey populations (one from Morocco and the other from Zimbabwe) was analysed. A total of 79 animals were sequenced for 313 bp of the cytochrome b gene, and 91 individuals for 383 bp of the D‐loop region or control‐region. Sequence comparisons and phylogenetic analyses of both Spanish and African populations revealed low diversity. Only six and seven haplotypes respectively were found in cytochrome b and the D‐loop region. Relatively low nucleotide diversity (π) values were detected in these populations. The π values, from the D‐loop region, ranged from 0.0006 to 0.0169 for the Catalana and Andaluza breeds, respectively. The obtained results seem to confirm the existence of two divergent maternal lineages of African origin (Equus asinus africanus and E. a. somaliensis). In this paper the genetic relationships between these breeds are analysed and compared with those obtained in other European populations. Also, the data on the genetic relationships between the populations, obtained from nuclear DNA (microsatellites) and mitochondrial DNA (cytochrome b and D‐loop region) is argued and interpreted.     

Genetic monitoring of ex situ African Penguin (Spheniscus demersus) populations in South Africa

The African Penguin (Spheniscus demersus) has suffered population declines and is listed in the IUCN Red List as Endangered. The species is endemic to the coast of southern Africa, and breeding colonies are distributed on the south-western coast of Africa. Currently, African Penguins are being kept in zoo and aquarium facilities throughout South Africa. In this study, molecular genetic data based on 12 microsatellite markers from 1 119 African Penguin samples from four facilities were generated in order to determine the level of genetic variation, population structure and differentiation, and effective population size to assist in the development of an effective captive management plan. Expected heterozygosity ranged from 0.57 to 0.62, and allelic richness from 4.2 to 5.1. However, based on differences between first- and second-generation captive birds, we conclude that the ex situ population is at risk of losing genetic variability in the future and management programmes should include exchange of birds between captive facilities in order to induce gene flow and increase effective population size. Adding individuals from in situ populations should also be considered in the future in cases where these birds cannot be rehabilitated. Molecular genetic analyses of wild penguin populations should be carried out for comparison, and to ascertain to what degree ‘in situ genetic diversity’ is represented among ex situ populations. With regular resampling and analyses, the extent of the effect of processes such as genetic drift on diversity in the ex situ penguin populations will become evident.     

Taxonomic status of Chinese blue sheep (Pseudois nayaur): new evidence of a distinct subspecies

The blue sheep is an endemic species to the Qinghai‐Tibet Plateau and surrounding regions. It has been regarded as having 2 subspecies: Pseudois nayaur nayaur and P. n. szechuanensis. However, such a classification remains controversial. Herein, we analyze 10 microsatellite loci and part of the mitochondrial control region for clarification in such taxonomic debates. We use samples from 168 individuals from 6 geographic populations covering almost all the distribution areas of the species in China to carry out comparisons. Phylogenetic trees derived from both the microsatellite and mitochondrial markers combined with the discriminant analysis of principal components (DAPC) and the STRUCTURE analysis reveal that the individuals in the Helan Mountains are well grouped with a distinct evolutionary lineage and are significantly different from the other populations of P. n. szechuanensis according to Fst values, implying that this isolated population should be categorized as a valid subspecies; namely, Pseudois nayaur alashanicus. The isolation‐by‐distance (IBD) analysis shows a significant positive relationship between genetic and geographical distances among the populations.     

Genomic differentiation with isolation by distance along a latitudinal gradient in the spotted-leg treefrog Polypedates megacephalus

The patterns of isolation by distance (IBD) entailing increased genetic differentiation among populations have aroused extensive concerns for evolutionary biologists. Although the IBD may act on spatial processes contributing to the genetic differentiation among populations in anuran species, the factors shaping the IBD of frogs among populations in natural systems are largely unknown. Here, we studied the genetic differentiation among six populations with 24 individuals of the spotted-leg treefrog along a latitudinal gradient (1860.31 km) based on 1020 single nucleotide polymorphisms from restriction site-associated DNA sequencing. The results showed that the genetic diversity differed significantly among populations and that the insular populations had higher genetic diversity than the mainland populations. Furthermore, we also found a significant genetic differentiation among populations (F_ST = 0.277) and no sign of inbreeding (F_UNI = −0.145). The IBD was detected for all populations, and a higher degree of the IBD was indicated when controlling for the effects of the isolation between Hainan and mainland populations caused by the Qiongzhou Strait. Our findings suggest that the form of the Qiongzhou Strait plays a key role in shaping the genetic diversity and population differentiation in treefrogs.     

Genetic diversity and structure in Asian native goat analyzed by newly developed SNP markers

In the current study, a total of 65 single nucleotide polymorphisms (SNPs) within the intron region were developed in goat (Capra hircus) by utilizing genomic information of cattle and sheep due to poor available genomic information on goat. Using these markers, we carried out genetic diversity and structure analyses for 10 Asian goat populations. The phylogenetic tree and principal components analysis showed good correspondence between clustered populations and their geographic locations. The STRUCTURE software analysis illustrated six divergent genetic structures among 10 populations. Myanmar and Cambodia populations showed high admixture patterns with different ancestry, suggesting genetic introgression into native goat populations. We also investigated the correlation between genetic diversity and geographic distance from a domestication center. This result showed a decreasing trend of genetic diversity according to the distance (P = 0.014). This result supported common consensus that western Asia is one of the centers of origin for modern Asian domestic goat.     

Is neutral genetic diversity related to quantitative variation in semen traits in bulls?

Conservation decisions based on neutral genetic diversity have been observed to promote retention of useful quantitative variation in biological populations. An experiment was undertaken to determine the association between microsatellite marker polymorphisms and phenotypic variation in semen production and cryosurvival traits in bulls. Thirty-five ejaculates were collected from ten bulls of two breeds and evaluated before and after cryopreservation for several semen traits. The bulls were also genotyped using a set of sixteen bovine-specific microsatellite marker loci. Fixation indices (F_ST), heterozygosity and Nei's genetic distance measures were computed from allele frequency data for each of the bulls. Molecular and phenotypic data were used to compute tri-distance matrices for the ten bulls and correlated using Mantel's test in GenAIEx 6.5. The study revealed extensive heterogeneity in semen traits, heterozygosity and F_ST values among the bulls. Large pairwise phenotypic and genetic distances were also observed. Correlation between pairwise genetic distances and phenotypic distances was significant and highly positive for sperm viability (r = .61, p < .001) and moderately positive for sperm motility (r = .40–42, p < .05) variables. For sperm morphology, ejaculate volume and sperm concentration, correlation with genetic distances was positive, low and not significantly different from zero (p > .05). A tendency for a triangular-shaped relationship between genetic and phenotypic distances for post-thaw motility and viability traits was also observed. Accordingly, association with neutral genetic diversity was absent for semen production traits and moderate to highly positive for sperm cryosurvival traits. Given these findings, conservation decisions based on neutral genetic diversity may capture variation in some adaptive traits, but not others.     

Diversity analysis of sheep breeds from Southern peninsular and Eastern regions of India

Information is presented on the genetic diversity and relationship among six Indian sheep breeds/populations belonging to the Southern peninsular and Eastern agroecological zones, based on microsatellite markers. Parameters of genetic variation, viz., allele diversity, observed heterozygosity, gene diversity and population inbreeding estimates, were calculated for the six breeds. The allele diversity ranged from 6.40 to 7.92, whereas the gene diversity varied from 0.617 to 0.727. The highest allele and gene diversity was observed for Nellore sheep, while the lowest was exhibited by Garole breed. Within population inbreeding estimate (F _IS) revealed a significant deficit of heterozygotes in Deccani, Madgyal, Nellore and Garole, whereas Ganjam and Chhotanagpuri sheep showed an excess of heterozygotes. The contribution of each breed to the total diversity of the breeds was quantified by the Weitzman approach. The marginal loss of diversity incurred with removal of Nellore and Garole breeds was higher (>27%), whereas removal of Deccani breed resulted in lowest loss of diversity (3.84%) from the set. Estimation of the genetic differentiation (F _ST) and genetic distance (D _A) between the pairs of breeds revealed a close relationship between Deccani and Madgyal sheep (F _ST = 0.017; D _A = 0.080) and greatest demarcation between Madgyal and Garole breeds (F _ST = 0.110; D _A = 0.622). The information generated would help in shaping genetic management and conservation programs for the sheep breeds under consideration.     

The assessment of genetic diversity within and among the eight subpopulations of Japanese Black cattle using 52 microsatellite markers

Japanese Black cattle are at risk for genetic homogeneity due to intensive use of a few sires. Therefore, assessment of the actual genetic diversity of this breed is important for future breeding plans. In the present study, we investigated the genetic diversity within and among eight subpopulations of Japanese Black cattle using 52 microsatellite markers. The parameters for genetic diversity of Japanese Black cattle were comparable to those of other cattle breeds, suggesting that the relatively high genetic diversity of the breed. However, upon comparison among the eight subpopulations, the Hyogo subpopulation showed markedly low genetic diversity. The results of the pairwise F_ST values, phylogenetic network and structure analysis indicated that the Hyogo population has remarkably high level of genetic differentiation from other populations, while Yamagata, Niigata, Hiroshima and Kagawa populations have low levels of genetic differentiation. Furthermore, multidimensional scaling plots indicated that individuals in some subpopulations were separated from individuals in the other subpopulations. We conclude that while the overall genetic diversity of Japanese Black cattle is still maintained at a relatively high level, that of a particular subpopulation is significantly reduced, and therefore the effective population size of the breed needs to be controlled by correct mating strategies.     

Genetic diversity,structure and individual assignment of Casta Navarra cattle: a well‐differentiated fighting bull population

The Casta Navarra lineage was one of the populations used to establish the fighting bull (FB) breed, and it has also been reproductively isolated from the others FBs. A total of 1284 individuals from two generations of 16 Casta Navarra herds were sampled to analyse their diversity, their genetic structure and the ability of 28 microsatellite markers to assign individuals to closely related populations. These animals were compared with closely related phylogenetic (FB) or geographical (Pirenaica and Monchina) populations. Hardy–Weinberg equilibrium analysis showed that 82% of the loci had a significant heterozygote deficit as a consequence of the Wahlund effect. The average proportion of genetic variation explained by farm differences was 9% by Wright's F_ST index. A phylogenetic tree constructed with a neighbour‐joining method based on Reynolds genetic distances and a Bayesian Markov chain Monte Carlo clustering approach revealed clear differences between farm groups that generally corresponded to historical information and could unambiguously differentiate Casta Navarra cattle from the other populations. The percentage of animals correctly assigned to the Casta Navarra population was 91.78% for a q threshold of >0.9. Admixture was only detected in 4.45% (q < 0.8) of the cattle. These results are relevant for the maintenance and development of diversity and conservation in the Casta Navarra population.     

Molecular genetic diversity and genetic structure of Vietnamese indigenous pig populations

The study characterized genetic diversity and genetic structure of five indigenous pig populations (Ha Lang, Muong Te, Mong Cai, Lung and Lung Pu), two wild pig populations (Vietnamese and Thai wild pigs) and an exotic pig breed (Yorkshire) using FAO/ISAG recommended 16 microsatellite markers in 236 samples. All estimated loci were very polymorphic indicated by high values of polymorphism information content (from 0.76 in S0225 to 0.92 in Sw2410). Indigenous populations had very high level of genetic diversity (mean He = 0.75); of all indigenous breeds, Lung Pu showed highest mean number of alleles (MNA = 10.1), gene diversity (He = 0.82), allele richness (5.33) and number of private alleles (10). Thirteen percentage of the total genetic variation observed was due to differences among populations. The neighbour‐joining dendrogram obtained from Nei's standard genetic distance differentiated eight populations into four groups including Yorkshire, two wild populations, Mong Cai population and a group of four other indigenous populations. The Bayesian clustering with the admixture model implemented in Structure 2.1 indicated seven possible homogenous clusters among eight populations. From 79% (Ha Lang) to 98% (Mong Cai). individuals in indigenous pigs were assigned to their own populations. The results confirmed high level of genetic diversity and shed a new light on genetic structure of Vietnam indigenous pig populations.