共查询到20条相似文献，搜索用时 15 毫秒
1.
Timothy J. Thurman Andre Szejner‐Sigal W. Owen McMillan 《Journal of evolutionary biology》2019,32(9):974-983
Hybrid zones have long been of interest to biologists as natural laboratories where we can gain insight into the processes of adaptation and speciation. Repeated sampling of individual hybrid zones has been particularly useful in elucidating the dynamic balance between selection and dispersal that maintains most hybrid zones. Here, we revisit a hybrid zone between Heliconius erato butterflies in Panamá for a third time over more than 30 years. We combine a novel Bayesian extension of stepped‐cline hybrid zone models with environmental data to understand the genetic and environmental causes of cline dynamics in this species. The cline has continued to move west, likely due to dominance drive, but has slowed and broadened. Environmental analyses suggest that widespread deforestation in Panamá could be leading to decreased avian predation and relaxed selection, causing the observed changes in cline dynamics. 相似文献
2.
Silu Wang Sievert Rohwer Kira Delmore Darren E. Irwin 《Journal of evolutionary biology》2019,32(11):1242-1251
Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long‐term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S. occidentalis) in the Washington Cascades was previously inferred to be moving from northern S. townsendi southwards towards S. occidentalis, based on plumage and behavioural patterns as well as a 2000‐km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2–3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987–94 and 2015–16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a ‘tension zone’ model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2–3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species. 相似文献
3.
K. E. Delmore R. A. Brenneman R. Lei C. A. Bailey A. Brelsford E. E. Louis Jr S. E. Johnson 《Journal of evolutionary biology》2013,26(8):1677-1690
Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south‐eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90‐km transect through the hybrid zone and examined variation in 26 microsatellites, the D‐loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection‐gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation. 相似文献
4.
Yann X. C. Bourgeois Joris A. M. Bertrand Boris Delahaie Hlne Holota Christophe Thbaud Borja Mil 《Molecular ecology》2020,29(6):1137-1153
Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process. 相似文献
5.
Identifying the mechanisms involved in the formation and maintenance of species is a central question in evolutionary biology, and distinguishing the selective drivers of populations’ divergence from demographic processes is of particular interest to better understand the speciation process. Hybrid zones are recognized to provide ideal places to investigate the genetic architecture of speciation and to identify the mechanisms allowing diverging species to maintain their integrity in the face of gene flow. Here, we studied two alpine butterfly species, Coenonympha macromma and C. gardetta, which can be found flying together and hybridizing in narrow contact zones in the southern French Alps. We characterized the genomic composition of individuals, their morphology and their local habitat requirements, within and around a hybrid zone. Genetic diversity analysis at 794 SNPs revealed that all individuals within the hybrid zone were highly admixed, which was not the case outside the hybrid zone. Cline analysis showed that, despite ongoing hybridization, 56 out of 122 loci differentially fixed or nearly so between the two species were impermeable to introgression across the sharp hybrid zone (9 km wide). We also found concordance in cline position and width among genetic, morphological and environmental variation, suggesting a coupling of different reproductive barriers. Habitat characteristics such as the presence of trees and shrubs and the start of the growing season were strongly associated with the genetic variation, and we found evidence of divergence at genetic markers associated with morphology and physiology, putatively involved in visual or environmental reproductive isolation. We discuss the various behavioural and ecological factors that might interplay to maintain current levels of divergence and gene flow between this species pair. 相似文献
6.
Ornamental traits typically advertise individual condition and can be costly to maintain. Plumage maintenance behaviour can increase plumage quality and positively influence female mate preference. We investigated this prediction by performing female mate‐choice trials and measuring male plumage maintenance behaviour in the European Serin Serinus serinus. More colourful males spent more time in plumage maintenance than less colourful males, but this was not correlated with body condition or ectoparasite load. Females preferred more colourful males, although this was not influenced by male plumage maintenance. Our findings highlight the fact that investment in plumage maintenance depends on individual coloration and this may reinforce the honesty of ornamental plumage and convey positive information for mate choice. 相似文献
7.
Elena Frederika Kappers Christiaan de Vries Anneke Alberda Sylvia Kuhn Mihai Valcu Bart Kempenaers Christiaan Both 《Journal of evolutionary biology》2020,33(9):1306-1315
How genetic polymorphisms are maintained in a population is a key question in evolutionary ecology. Previous work on a plumage colour polymorphism in the common buzzard Buteo buteo suggested heterozygote advantage as the mechanism maintaining the co‐existence of three morphs (light, intermediate and dark). We took advantage of 20 years of life‐history data collected in a Dutch population to replicate earlier studies on the relationship between colour morph and fitness in this species. We examined differences between morphs in adult apparent survival, breeding success, annual number of fledglings produced and cumulative reproductive success. We found that cumulative reproductive success differed among morphs, with the intermediate morph having highest fitness. We also found assortative mating for colour morph, whereby assortative pairs were more likely to produce offspring and had longer‐lasting pair bonds than disassortative pairs. Over the 20‐year study period, the proportion of individuals with an intermediate morph increased. This apparent evolutionary change did not just arise from selection on individual phenotypes, but also from fitness benefits of assortative mating. The increased frequency of intermediates might also be due to immigration or drift. We hypothesize that genetic variation is maintained through spatial variation in selection pressures. Further studies should investigate morph‐dependent dispersal behaviour and habitat choice. 相似文献
8.
Isolde van Riemsdijk Roger K. Butlin Ben Wielstra Jan W. Arntzen 《Molecular ecology》2019,28(5):1070-1083
Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes. 相似文献
9.
Guillaume de Lafontaine Julien Prunier Sébastien Gérardi Jean Bousquet 《Molecular ecology》2015,24(20):5229-5247
The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes. 相似文献
10.
Jessica L. Ward Mike J. Blum David M. Walters Brady A. Porter Noel Burkhead Byron Freeman 《Evolutionary Applications》2012,5(4):380-392
The erosion of species boundaries can involve rapid evolutionary change. Consequently, many aspects of the process remain poorly understood, including the formation, expansion, and evolution of hybrid swarms. Biological invasions involving hybridization present exceptional opportunities to study the erosion of species boundaries because timelines of interactions and outcomes are frequently well known. Here, we examined clinal variation across codominant and maternally inherited genetic markers as well as phenotypic traits to characterize the expansion and evolution of a hybrid swarm between native Cyprinella venusta and invasive Cyprinella lutrensis minnows. Discordant introgression of phenotype, microsatellite multilocus genotype, and mtDNA haplotype indicates that the observable expansion of the C. venusta × C. lutrensis hybrid swarm is a false invasion front. Both parental and hybrid individuals closely resembling C. lutrensis are numerically dominant in the expansion wake, indicating that the non‐native parental phenotype may be selectively favored. These findings show that cryptic introgression can extend beyond the phenotypic boundaries of hybrid swarms and that hybrid swarms likely expand more rapidly than can be documented from phenotypic variation alone. Similarly, dominance of a single parental phenotype following an introduction event may lead to instances of species erosion being mistaken for species displacement without hybridization. 相似文献
11.
12.
Gregory J. Glotzbecker David M. Walters Michael J. Blum 《Evolutionary Applications》2016,9(6):741-755
Unstable hybrid swarms that arise following the introduction of non‐native species can overwhelm native congeners, yet the stability of invasive hybrid swarms has not been well documented over time. Here, we examine genetic variation and clinal stability across a recently formed hybrid swarm involving native blacktail shiner (Cyprinella venusta) and non‐native red shiner (C. lutrensis) in the Upper Coosa River basin, which is widely considered to be a global hot spot of aquatic biodiversity. Examination of phenotypic, multilocus genotypic, and mitochondrial haplotype variability between 2005 and 2011 revealed that the proportion of hybrids has increased over time, with more than a third of all sampled individuals exhibiting admixture in the final year of sampling. Comparisons of clines over time indicated that the hybrid swarm has been rapidly progressing upstream, but at a declining and slower pace than rates estimated from historical collection records. Clinal comparisons also showed that the hybrid swarm has been expanding and contracting over time. Additionally, we documented the presence of red shiner and hybrids farther downstream than prior studies have detected, which suggests that congeners in the Coosa River basin, including all remaining populations of the threatened blue shiner (Cyprinella caerulea), are at greater risk than previously thought. 相似文献
13.
Katie L. Smith Joshua M. Hale Michael R. Kearney Jeremy J. Austin Jane Melville 《Molecular ecology》2013,22(7):1869-1883
Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization. 相似文献
14.
Tahani El Ayari Najoua Trigui El Menif Carlos Saavedra David Cordero Frédérique Viard Nicolas Bierne 《Ecology and evolution》2017,7(21):9016-9026
Molecular approaches have proven efficient to identify cryptic lineages within single taxonomic entities. Sometimes these cryptic lineages maybe previously unreported or unknown invasive taxa. The genetic structure of the marine gastropod Stramonita haemastoma has been examined in the Western Mediterranean and North‐Eastern Atlantic populations with mtDNA COI sequences and three newly developed microsatellite markers. We identified two cryptic lineages, differentially fixed for alternative mtDNA COI haplogroups and significantly differentiated at microsatellite loci. The mosaic distribution of the two lineages is unusual for a warm‐temperate marine invertebrate with a teleplanic larval stage. The Atlantic lineage was unexpectedly observed as a patch enclosed in the north of the Western Mediterranean Sea between eastern Spain and the French Riviera, and the Mediterranean lineage was found in Macronesian Islands. Although cyto‐nuclear disequilibrium is globally maintained, asymmetric introgression occurs in the Spanish region where the two lineages co‐occur in a hybrid zone. A first interpretation of our results is mito‐nuclear discordance in a stable postglacial hybrid zone. Under this hypothesis, though, the location of genetic discontinuities would be unusual among planktonic dispersers. An alternative interpretation is that the Atlantic lineage, also found in Senegal and Venezuela, has been introduced by human activities in the Mediterranean area and is introgressing Mediterranean genes during its propagation, as theoretically expected. This second hypothesis would add an additional example to the growing list of cryptic marine invasions revealed by molecular studies. 相似文献
15.
We investigated a hybrid zone between two major lineages of fence lizards (Sceloporus cowlesi and Sceloporus tristichus) in the Sceloporus undulatus species complex in eastern Arizona. This zone occurs in an ecotone between Great Basin Grassland and Conifer Woodland habitats. We analysed spatial variation in mtDNA (N=401; 969 bp), chromosomes (N=217), and morphology (N=312; 11 characters) to characterize the hybrid zone and assess species limits. A fine-scale population level phylogenetic analysis refined the boundaries between these species and indicated that four nonsister mtDNA clades (three belonging to S. tristichus and one to S. cowlesi) are sympatric at the centre of the zone. Estimates of cytonuclear disequilibria in the population closest to the centre of the hybrid zone suggest that the S. tristichus clades are randomly mating, but that the S. cowlesi haplotype has a significant nonrandom association with nuclear alleles. Maximum-likelihood cline-fitting analyses suggest that the karyotype, morphology, and dorsal colour pattern clines are all coincident, but the mtDNA cline is skewed significantly to the south. A temporal comparison of cline centres utilizing karyotype data collected in the early 1970s and in 2002 suggests that the cline may have shifted by approximately 1.5 km to the north over a 30-year period. The recent northward expansion of juniper trees into the Little Colorado River Basin resulting from intense cattle overgrazing provides a plausible mechanism for a shifting hybrid zone and the introgression of the mtDNA haplotypes, which appear to be selectively neutral. It is clear that complex interactions are operating simultaneously in this contact zone, including the formation of hybrids between populations within S. tristichus having diagnostic mtDNA, morphology, karyotypes, and dorsal colour patterns, and secondary contact between these and a distantly related yet morphologically cryptic mtDNA lineage (S. cowlesi). 相似文献
16.
Siu F. Lee Ying Chen Eyre‐Walker Rahul V. Rane Caroline Reuter Giovanna Vinti Lea Rako Linda Partridge Ary A. Hoffmann 《Molecular ecology》2013,22(10):2716-2725
In many invertebrates, body size shows genetically based clines, with size increasing in colder climates. Large body size is typically associated with prolonged development times. We consider variation in the CNS‐specific gene neurofibromin 1 (Nf1) and its association with body size and development time. We identified two major Nf1 haplotypes in natural populations, Nf1‐insertion‐A and Nf1‐deletion‐G. These haplotypes are characterized by a 45‐base insertion/deletion (INDEL) in Nf1 intron 2 and an A/G synonymous substitution (locus L17277). Linkage disequilibrium (LD) between the INDEL and adjacent sites is high but appears to be restricted within the Nf1 gene interval. In Australia, the frequency of the Nf1‐insertion‐A haplotype increases with latitude where wing size is larger, independent of the chromosomal inversion In(3R)Payne. Unexpectedly, the Nf1‐insertion‐A haplotype is negatively associated with wing size. We found that the Nf1‐insertion‐A haplotype is enriched in females with shorter development time. This suggests that the Nf1 haplotype cline may be driven by selection for development time rather than size; females from southern (higher latitude) D. melanogaster populations maintain a rapid development time despite being relatively larger, and the higher incidence of Nf1‐insertion‐A in Southern Australia may contribute to this pattern, whereas the effects of the Nf1 haplotypes on size may be countered by other loci with antagonistic effects on size and development time. Our results point to the potential complexity involved in identifying selection on genetic variants exhibiting pleiotropic effects when studies are based on spatial patterns or association studies. 相似文献
17.
Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. FST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at FST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (FST = 0.03), but seven loci (2%) were strong FST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated. 相似文献
18.
Ana P. Moraes Mariana Chinaglia Clarisse Palma‐Silva Fábio Pinheiro 《Ecology and evolution》2013,3(11):3824-3837
Interspecific hybridization is a primary cause of extensive morphological and chromosomal variation and plays an important role in plant species diversification. However, the role of interploidal hybridization in the formation of hybrid swarms is less clear. Epidendrum encompasses wide variation in chromosome number and lacks strong premating barriers, making the genus a good model for clarifying the role of chromosomes in postzygotic barriers in interploidal hybrids. In this sense, hybrids from the interploidal sympatric zone between E. fulgens (2n = 2x = 24) and E. puniceoluteum (2n = 4x = 56) were analyzed using cytogenetic techniques to elucidate the formation and establishment of interploidal hybrids. Hybrids were not a uniform group: two chromosome numbers were observed, with the variation being a consequence of severe hybrid meiotic abnormalities and backcrossing with E. puniceoluteum. The hybrids were triploids (2n = 3x = 38 and 40) and despite the occurrence of enormous meiotic problems associated with triploidy, the hybrids were able to backcross, producing successful hybrid individuals with broad ecological distributions. In spite of the nonpolyploidization of the hybrid, its formation is a long‐term evolutionary process rather than a product of a recent disturbance, and considering other sympatric zones in Epidendrum, these events could be recurrent. 相似文献
19.
ALEXANDRE ROULIN RETO BURRI SYLVAIN ANTONIAZZA 《Biological journal of the Linnean Society. Linnean Society of London》2011,102(3):573-582
Climate change acts as a major new selective agent on many organisms, particularly at high latitudes where climate change is more pronounced than at lower latitudes. Studies are required to predict which species are at a high risk of extinction and whether certain phenotypes may be more affected by climate change than others. The identification of susceptible phenotypes is important for evaluating the potential negative effect of climate change on biodiversity at the inter‐ and intraspecific levels. Melanin‐based coloration is an interesting and easily accessible candidate trait because, within certain species, reddish pheomelanin‐based coloration is associated with adaptations to warm climates. However, it is unclear whether the same holds among species. We tested one prediction of this hypothesis in four owl genera (wood, scops, screech, and pygmy owls), namely that darker reddish species are more prevalent near the equator than polewards. Our comparative analysis is consistent with this prediction for the northern hemisphere, suggesting that pale reddish species may be adapted to cold climates and dark reddish species to warmer climates. Thus, climate change may have a larger negative impact on pale pheomelanic owls and favour dark pheomelanic species. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 573–582. 相似文献
20.
House mice (Mus musculus) are human commensals and have served as a primary model in biomedical, ecological and evolutionary research. Although there is detailed knowledge of the biogeography of house mice in Europe, little is known of the history of house mice in China, despite the fact that China encompasses an enormous portion of their range. In the present study, 535 house mice caught from 29 localities in China were studied by sequencing the mitochondrial D‐loop and genotyping 10 nuclear microsatellite markers distributed on 10 chromosomes. Phylogenetic analyses revealed two evolutionary lineages corresponding to Mus musculus castaneus and Mus musculus musculus in the south and north, respectively, with the Yangtze River approximately representing the boundary. More detailed analyses combining published sequence data from mice sampled in neighbouring countries revealed the migration routes of the two subspecies into China: M. m. castaneus appeared to have migrated through a southern route (Yunnan and Guangxi), whereas M. m. musculus entered China from Kazakhstan through the north‐west border (Xinjiang). Bayesian analysis of mitochondrial sequences indicated rapid population expansions in both subspecies, approximately 4650–9300 and 7150–14 300 years ago for M. m. castaneus and M. m. musculus, respectively. Interestingly, the migration routes of Chinese house mice coincide with the colonization routes of modern humans into China, and the expansion times of house mice are consistent with the development of agriculture in southern and northern China, respectively. Finally, our study confirmed the existence of a hybrid zone between M. m. castaneus and M. m. musculus in China. Further study of this hybrid zone will provide a useful counterpart to the well‐studied hybrid zone between M. m. musculus and Mus musculus domesticus in central Europe. 相似文献