Diversification and Demography of the Oriental Garden Lizard (Calotes versicolor) on Hainan Island and the Adjacent Mainland

The Oriental garden lizard (Calotes versicolor) is one of the few non-gekkonid lizards that are geographically widespread in the tropics. We investigated its population dynamics on Hainan Island and the adjacent mainland of China and Vietnam, focusing on the impact of cyclic upheaval and submergence of land bridges during the Pleistocene. Our Bayesian phylogenetic analysis reveals two mitochondrial lineages, A and B, which are estimated to have coalesced about 0.26 million years ago (95% credibility interval: 0.05–0.61 million years ago). Lineage A contains individuals mainly from central and southern Wuzhi Mountain on Hainan Island, whereas lineage B mainly comprises individuals from other sites on the island plus the adjacent mainland. The estimated coalescence times within lineages A (0.05 million years ago) and B (0.13 million years ago) fall within a period of cyclical land-bridge formation and disappearance in the Pleistocene. A spatial analysis of molecular variance identified two distinct population groupings: I, primarily containing lineage A, and II, mainly consisting of lineage B. However, haplotypes from lineages A and B occur sympatrically, suggesting that gene flow is ongoing. Neither Wuzhi Mountain nor Qiongzhou Strait and Gulf of Tonkin act as barriers to gene flow among C. versicolor populations. Analyses of the data using mismatch distributions and extended Bayesian skyline plots provide evidence of a relatively stable population size through time for Group I, and moderate population expansions and contractions during the end of the Pleistocene for Group II. We conclude that the phylogeographical patterns of C. versicolor are the combined product of Pleistocene sea-level oscillations and nonphysical barriers to gene flow.     

Phylogeography and introgressive hybridization of the ground beetle Carabus yamato in Japan based on mitochondrial gene sequences

To study the phylogeography of the ground beetle Carabus yamato in Japan, we compared 1,020-bp sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene from 373 specimens from 37 localities with those of three parapatric species (C. albrechti, C. kimurai, and C. japonicus) that might share mitochondrial lineages with C. yamato through introgressive hybridization. We found 81 haplotypes from C. yamato. Of these, 17 haplotypes were considered to be of an introgressed lineage from C. albrechti, based on the phylogeny and geographic distribution. In addition, one haplotype of C. kimurai was likely an introgressant from C. yamato. Putative introgression events among the four species were restricted to these two directional cases. We analyzed the phylogeography of C. yamato using nested clade phylogeographical analysis and population genetic parameters. The mitochondrial lineages of C. yamato were estimated to have diverged no more than approximately 1.12 million years ago, implying that the estimated historical events occurred after the Early Pleistocene. Carabus yamato was inferred to have experienced a contraction of its distribution range, followed by recent range expansion. Populations in the western and eastern regions, segregated by Ise Bay and the Nobi and Okazaki Plains, diverged in the mitochondrial clades. The northern and most western populations possessed one clade only (except an introgressed lineage), whereas eastern and some southwestern populations possessed several diverged clades, which were considered to be ancestral; these populations may have been associated with refugia during glacial periods.     

Two mitochondrial lineages occur in the Asian corn borer, Ostrinia furnacalis (Lepidoptera: Crambidae), in Japan

The genealogy and diversity of the mitochondrial cytochrome oxidase subunit II (COII) gene were investigated for Ostrinia furnacalis in Japan. A preliminary examination of mitochondrial lineages in China and the Philippines was also made. Two lineages (A and B) were found in the COII gene. Lineage A was frequent throughout the Japanese main islands (Hokkaido, Honshu, Shikoku and Kyushu), while the frequency of lineage B varied among these islands. No clear patterns of geographical population structure were found. Population genetic features suggested that the O. furnacalis population harboring the lineage A mitochondria expanded in the recent past, while lineage B showed weak signals of a population expansion. It is not clear whether the two lineages of mtDNA evolved in separate or identical geographical regions. We discuss two hypotheses regarding the two lineages of mtDNA: a cryptic race/species hypothesis and a selective sweep hypothesis.     

Independent maternal origin of Chinese swamp buffalo (Bubalus bubalis)

To obtain more knowledge on the origin and genetic diversity of the swamp buffalo (Bubalus bubalis) in China, the complete mitochondrial D-loop sequences of 119 samples representing seven native types were compared. Two mitochondrial DNA (mtDNA) lineages (lineages A and B) were determined for the Chinese swamp buffalo. Examination of the diversity patterns suggest that lineage A has undergone a population expansion event. Divergence of lineages A and B was estimated at 18,000 years ago. Combined analyses of mtDNA sequences from Chinese, Indian, Brazilian/Italian and Southeast Asian/Australian buffalo samples showed independent domestication events in the swamp buffalo from China and the river buffalo from the India subcontinent. The spread of swamp and river buffalo from China and India respectively to mainland Southeast Asia suggests that Southeast Asia is a hybrid zone for buffalo. Our data support the hypothesis of the evolution of domesticated swamp and river buffalo from ancestral swamp-like animals. These ancestral animals were extensively distributed across mainland Asia and most likely are represented today by the wild Asian buffalo (Bubalus arnee).     

The impact of Quaternary climate oscillations on divergence times and historical population sizes in Thylamys opossums from the Andes

Climate oscillations during the Quaternary altered the distributions of terrestrial animals at a global scale. In mountainous regions, temperature fluctuations may have led to shifts in range size and population size as species tracked their shifting habitats upslope or downslope. This creates the potential for both allopatric speciation and population size fluctuations, as species are either constrained to smaller patches of habitat at higher elevations or able to expand into broader areas at higher latitudes. We considered the impact of climate oscillations on three pairs of marsupial species from the Andes (Thylamys opossums) by inferring divergence times and demographic changes. We compare four different divergence dating approaches, using anywhere from one to 26 loci. Each pair comprises a northern (tropical) lineage and a southern (subtropical to temperate) lineage. We predicted that divergences would have occurred during the last interglacial (LIG) period approximately 125 000 years ago and that population sizes for northern and southern lineages would either contract or expand, respectively. Our results suggest that all three north–south pairs diverged in the late Pleistocene during or slightly after the LIG. The three northern lineages showed no signs of population expansion, whereas two southern lineages exhibited dramatic, recent expansions. We attribute the difference in responses between tropical and subtropical lineages to the availability of ‘montane‐like’ habitats at lower elevations in regions at higher latitudes. We conclude that climate oscillations of the late Quaternary had a powerful impact on the evolutionary history of some of these species, both promoting speciation and leading to significant population size shifts.     

Mitochondrial genetic variation in European, African and Indian cattle populations

Mitochondrial DNA from representative animals of 13 different cattle breeds was assayed for restriction fragment length polymorphisms (RFLP) to determine phylogenetic relationships and levels of variation among breeds; 16 different mitotypes were found, described by 20 polymorphisms. Within these 16 mitotypes two major lineages were apparent: an Afro-European and an Asian type. These were found to differ at over 2.3% of sites surveyed. None of the mitotypes found in the Asian lineage was detectable in the Afro-European lineage and vice versa. Within each of the major mitotypes there were no further significant differences within or among breeds. Using rates of mitochondrial evolution estimated from other species, the two lineages were estimated to have diverged between 575000 and 1150000 years ago; well outside the 10000 years bp timeframe postulated by a single domestication hypothesis. The results presented are concordant with those generated in other studies and provide strong evidence for an independent domestication of Asian Bos indicus. Furthermore, the grouping of all African indicine populations within the clade containing all Bos taurus lineages points to the hybrid origins of the humped cattle of that continent.     

Genomic analysis of the emergent aquatic plant Sparganium stoloniferum provides insights into its clonality,local adaptation and demographic history

Clonal propagation and extensive dispersal of seeds and asexual propagules are two important features of aquatic plants that help them adapt to aquatic environments. Accurate measurements of clonality and effective clonal dispersal are essential for understanding the evolution of aquatic plants. Here, we first assembled a high-quality chromosome-level genome of a widespread emergent aquatic plant Sparganium stoloniferum to provide a reference for its population genomic study. We then performed high-depth resequencing of 173 individuals from 20 populations covering different basins across its range in China. Population genomic analyses revealed three genetic lineages reflecting the northeast (NE), southwest (SW) and northwest (NW) of its geographical distribution. The NE lineage diverged in the middle Pleistocene while the SW and NW lineages diverged until about 2400 years ago. Clonal relationship analyses identified nine populations as monoclonal population. Dispersal of vegetative propagules was identified between five populations covering three basins in the NE lineage, and dispersal distance was up to 1041 km, indicating high dispersibility in emergent aquatic plant species. We also identified lineage-specific positively selected genes that are likely to be involved in adaptations to saline wetlands and high-altitude environments. Our findings accurately measure the clonality, determine the dispersal range and frequency of vegetative propagules, and detect genetic signatures of local adaptation in a widespread emergent aquatic plant species, providing new perspectives on the evolution of aquatic plants.     

Pan-African Genetic Structure in the African Buffalo (Syncerus caffer): Investigating Intraspecific Divergence

The African buffalo (Syncerus caffer) exhibits extreme morphological variability, which has led to controversies about the validity and taxonomic status of the various recognized subspecies. The present study aims to clarify these by inferring the pan-African spatial distribution of genetic diversity, using a comprehensive set of mitochondrial D-loop sequences from across the entire range of the species. All analyses converged on the existence of two distinct lineages, corresponding to a group encompassing West and Central African populations and a group encompassing East and Southern African populations. The former is currently assigned to two to three subspecies (S. c. nanus, S. c. brachyceros, S. c. aequinoctialis) and the latter to a separate subspecies (S. c. caffer). Forty-two per cent of the total amount of genetic diversity is explained by the between-lineage component, with one to seventeen female migrants per generation inferred as consistent with the isolation-with-migration model. The two lineages diverged between 145 000 to 449 000 years ago, with strong indications for a population expansion in both lineages, as revealed by coalescent-based analyses, summary statistics and a star-like topology of the haplotype network for the S. c. caffer lineage. A Bayesian analysis identified the most probable historical migration routes, with the Cape buffalo undertaking successive colonization events from Eastern toward Southern Africa. Furthermore, our analyses indicate that, in the West-Central African lineage, the forest ecophenotype may be a derived form of the savanna ecophenotype and not vice versa, as has previously been proposed. The African buffalo most likely expanded and diverged in the late to middle Pleistocene from an ancestral population located around the current-day Central African Republic, adapting morphologically to colonize new habitats, hence developing the variety of ecophenotypes observed today.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

Genetic differentiation between Príncipe Island and mainland populations of the grey parrot (Psittacus erithacus), and implications for conservation

The range of the grey parrot (Psittacus erithacus), one of the most heavily harvested bird species for the international pet trade, spans the forest belt of Central and West Africa and includes the oceanic island of Príncipe (Gulf of Guinea). Morphological variation led to the recognition of two mainland subspecies (P. e. erithacus and P. e. timneh). The population from Príncipe was originally described as a separate species (P. princeps) but is currently included in the nominate race. We used 1932 bp of the mitochondrial genome to clarify the genealogical affinities between the two currently-recognized mainland subspecies and the Príncipe population. Sampling included 20 individuals from Príncipe, 17 from P. e. erithacus, and 13 from P. e. timneh. We found that the two mainland subspecies form two independent lineages, having diverged up to 2.4 million years ago (Ma), and that the Príncipe population is composed of two lineages that diverged in allopatry. The most common lineage is descended from the first colonizers and evolved in isolation for up to 1.4 Ma. Contrary to current understanding, this 'Príncipe lineage' is more closely related to the timneh than to the erithacus subspecies. The second lineage consists of P. e. erithacus birds from the nearby mainland that colonized the island in recent times. The evolutionary dynamics of the grey parrot population of Príncipe are primarily characterized by isolation, with new genetic variation being added through rare immigration events. The heavily harvested Príncipe grey parrot population should therefore be treated as an independent conservation unit.     

Population genetics of the high elevation black fly Simulium (Nevermannia) feuerborni Edwards in Thailand

In this study, we examined the genetic structure and population history of the high elevation black fly Simulium feuerborni in Thailand at both cytogenetic and molecular genetic levels. Cytological examination revealed two cytoforms differentiated by fixed chromosome inversions. The distributions of the cytoforms were associated with geographic origins. Cytoform A was found in the lower north and northeast, and cytoform B was found in the upper northern region of Thailand. Molecular data based on the mitochondrial cytochrome oxidase subunit I (COI) barcoding sequence supports the separation of the cytoforms. The average sequence divergence between the two cytoforms was 3.75%, which is higher than the threshold value for the species level based on a COI barcoding sequence. Median joining network clearly differentiated the haplotypes of the cytoforms into different lineages. Population pairwise F_ST and amova analyses reveal significant genetic differentiation between cytoforms. This indicates that the low land areas separating these populations act as a gene flow barrier. No genetic differentiation was detected within cytoforms. This could be due to a recent sharing of population history. Mismatch distribution analysis revealed population expansion in the northern lineage of the cytoform B approximately 220 000 years ago. More recent expansion (32 000 years ago) was found in the lower north and northeast (cytoform A) lineage. The demographic history of S. feuerborni mirrored previous findings in black flies and other insect species in Thailand. This indicates the important role of Pleistocene climatic change on genetic structure and diversity of Southeast Asian mainland species.     

Cryptic biodiversity and phylogeographical patterns in a snapping shrimp species complex

Recent investigations suggest that marine biodiversity may be much higher than earlier estimates, and an important hidden source of diversity in marine systems is the phenomenon of cryptic species complexes. Such complexes are informative models for research into the evolutionary processes that govern species compositions of marine fauna. The snapping shrimp genera Alpheus and Synalpheus are known to harbour large numbers of cryptic species; here, I characterize the genetic structure of the Alpheus armillatus species complex in the northern Caribbean, west Atlantic, and Gulf of Mexico using mitochondrial and nuclear sequence data. Over this geographical region, the complex harbours at least three lineages that are probable reproductively isolated species; all major lineages diverged subsequent to the close of the Isthmus of Panama. Only one lineage was present in the Gulf of Mexico, whereas outside the Gulf of Mexico there was no clear tendency for lineage dominance by geographical region, as most sites were populated by shrimp from at least two lineages. However, within each lineage, there was strong evidence of population genetic differentiation between geographical regions. All lineages showed strong signals of demographic expansion, and one lineage showed sharply reduced genetic diversity, suggestive of past population bottlenecks or recently founded populations with low gene flow from other sites. These results show that evolutionary processes leading to divergence and speciation have been common and recent in the snapping shrimp, and suggest that connectivity among shrimp populations may be limited.     

长爪沙鼠线粒体细胞色素b 基因的遗传变异及地理分化

本文基于6 个种群52 个个体的线粒体细胞色素b (Cyt b)基因全序列的遗传变异分析,探讨了我国长爪沙鼠种群的遗传结构和种群历史。在1140 bp Cyt b基因的碱基序列中,共发现了74个变异位点,其中转换和颠换位点分别为65和9个,共定义了37个单倍型。种群差异分析表明包头种群与其他种群之间都存在着显著的遗传差异,而其他种群间遗传差异则不显著。52 个个体在系统发生树中明显聚为两支(Clade A 和Clade B),其中Clade B 仅含有贺兰山东部种群中的3 种单倍型,而其余的单倍型分布于Clade A 中。来自Clade A 的49个个体的歧点分布分析结果呈单峰状,提示长爪沙鼠种群可能在历史上经历了种群增长和扩张事件,这一结果同时也得到了Tajima 检验结果(D=-1.86,P<0.05)和Fu 检验结果(Fs=-21.89,P <0.05) 的支持。通过种群扩张系数(τ)和分子钟的推算提示该种群扩张事件大约发生于11 万年前。我们的研究结果表明,长爪沙鼠在第四纪冰期中形成了明显的隔离分化,而贺兰山东部可能作为冰期避难所,包含于所有的谱系分支中,在末次间冰期种群进行了强烈的扩张,从而形成当今的分布格局。     

Contrasting genetic and morphological differentiation among geographical lineages of a stenotopic miniature rasborine,Boraras maculatus,in Peninsular Malaysia

The variability in the stenotopic miniature rasborine Boraras maculatus (Cypriniformes: Danionidae: Rasborinae) across acidic-water habitats of Peninsular Malaysia (PM) was investigated using two molecular markers (the mitochondrial cytochrome c oxidase subunit I [COI] gene and the nuclear rhodopsin gene), as well as morphological evidence. Molecular phylogenetic analyses revealed differentiation among populations of B. maculatus in PM with the distinction of four allopatric lineages. Each of them was recognized as a putative species by automatic species delimitation methods. These lineages diverged from each other between 7.4 and 1.9 million years ago. A principal component analysis (PCA) was conducted to examine the multivariate variation in 11 morphometric measurements among three of these lineages. PCA results showed a significant overlap in morphological characteristics among these lineages. Additionally, a photograph-based machine learning approach failed to fully differentiate these lineages, suggesting limited morphological differentiation. B. maculatus represents a case of morphological stasis in a stenotopic miniature species. Strong habitat preference, coupled with long-term habitat fragmentation, may explain why each lineage of B. maculatus has a restricted distribution and did not disperse to other regions within and outside of PM, despite ample possibilities when the Sunda shelf was emerged and drained by large paleodrainages for most of the past 7 million years. The conservation status of B. maculatus and its peat swamp habitats are discussed, and it is concluded that peat swamps comprise several evolutionary units. Each of these units is considered a conservation unit and deserves appropriate protection.     

Evidence of three maternal lineages in Near Eastern sheep supporting multiple domestication events

The variability of mtDNA was analysed in local sheep breeds reared throughout Turkey, for which a fragment of the D-loop region and the complete cytochrome b were sequenced. Phylogenetic analyses performed independently for the D-loop and the Cyt b gene revealed three clearly separated clusters indicating three major maternal lineages, two of which had been previously described as types B and A. The new type, C, was present in all the breeds analysed and showed considerable mtDNA variability. Divergence time was obtained on the basis of Cyt b gene and was estimated to be around 160,000-170,000 years ago for lineages B and A, whereas the divergence of lineage C proved to have occurred earlier (between 450,000 and 750,000 years ago). These times greatly predate domestication and suggest that the origin of modern sheep breeds was more complex than previously thought and that at least three independent sheep domestication events occurred. Our results, together with archaeological information and the current wild sheep populations in the Near East region support the high importance of this area in the sheep domestication process. Finally, the evidence of a third maternal lineage has important implications regarding the history of modern sheep.     

Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.     

Genetic variation and phylogeography of free-living mouse species (genus Mus) in the Balkans and the Middle East

This work presents a study of the distribution and pattern of variation throughout the ranges of three free-living mouse species of the genus Mus-M. macedonicus, M. spicilegus, and a M. cypriacus - based on sequencing of two segments of the mitochondrial DNA (mtDNA) control region. The study shows a similar level of variability in the three species and suggests their recent population expansion. The highest proportion of variation is found within populations indicating low genetic structuring. Phylogenetic analysis confirms the significant divergence of a mitochondrial lineage of M. macedonicus from Israel, recently described as a new subspecies, M. macedonicus spretoides. Conversely, no genetic hiatus is revealed between European and Asian populations of M. macedonicus macedonicus. Although phylogenetic relationships among M. spicilegus populations could not be unravelled precisely, the results suggest a recent westward expansion of the species. The mtDNA divergence between M. macedonicus and M. spicilegus is 7.3%, suggesting their split between c. 700,000 and 1 million years ago. These dates correspond with a coalescent estimate about 720,000 years ago. On the other hand, M. cypriacus appeared almost twice as divergent from the former species (4.5%) as from the latter (8.8%) suggesting a divergence of c. 430,000-610,000 years ago (coalescent approximately 490,000 years ago) and 830,000-1.2 million years ago (coalescent approximately 780,000 years ago), respectively. Approximate times of population expansion have also been estimated for all taxa and groups of populations. Existence of several glacial refuges and various colonization scenarios are discussed; since all estimated divergence times fall within interglacial periods it seems that climatic oscillations did not play a crucial role in the evolution of the three species.     

Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

R1 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of arthropods. One aspect of R1 evolution that has been difficult to explain is the presence of divergent lineages of R1 in the rDNA loci of the same species. Multiple lineages should compete for a limited number of insertion sites, in addition to being subject to the concerted evolution processes homogenizing the rRNA genes. The presence of multiple lineages suggests either the ability of the elements to overcome these factors and diverge within rDNA loci, or the introduction of new lineages by horizontal transmission. To address this issue, we attempted to characterize the complete set of R1 elements in the rDNA locus from five Drosophila species groups (melanogaster, obscura, testacea, quinaria, and repleta). Two major R1 lineages, A and B, that diverged about 100 MYA were found to exist in Drosophila. Elements of the A lineage were found in all 35 Drosophila species tested, while elements of the B lineage were found in only 11 species from three species groups. Phylogenetic analysis of the R1 elements, supported by comparison of their rates of nucleotide sequence substitution, revealed that both the A and the B lineages have been maintained by vertical descent. The B lineage was less stable and has undergone numerous, independent elimination events, while the A lineage has diverged into three sublineages, which were, in turn, differentially stable. We conclude that while the differential retention of multiple lineages greatly complicates its phylogenetic history, the available R1 data continue to be consistent with the strict vertical descent of these elements.     

Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.