Mating system variation and morph fluctuations in a polymorphic lizard

In polymorphic male painted dragon lizards (Ctenophorus pictus), red males win staged contests for females over yellow males, and yellow males have greater success in staged sperm competition trials than red males. This predicts different reproductive strategies in the wild with red males being more coercive or better mate guarders than yellow males. Yellow males would be expected to sire more offspring per copulation and have a greater proportion of offspring from clutches with mixed paternity. However, here we show using microsatellites that the frequency of mixed paternity in the wild is low (< 20% on average across years), that all morphs on average have the same number of offspring sired per year, and that mating system variation (polyandry vs. monandry) is strongly correlated with perch density on male territories. Furthermore, a logistic regression on male successful vs. unsuccessful mate acquisition showed that red males were under negative selection when they dominated the population, which suggests ongoing frequency dependent selection on male colouration.     

Population genetic structure and demographic history of the black fly vector,Simulium nodosum in Thailand

An understanding of the genetic structure and diversity of vector species is crucial for effective control and management. In this study, mitochondrial DNA sequences were used to examine the genetic structure, diversity and demographic history of a black fly vector, Simulium nodosum Puri (Diptera: Simuliidae), in Thailand. A total of 145 sequences were obtained from 10 sampling locations collected across geographical ranges in the country. Low genetic diversity was found in populations of S. nodosum that could be explained by the recent population history of this species. Demographic history analysis revealed a signature of demographic expansion dating back to only 2600–5200 years ago. Recent population expansion in S. nodosum possibly followed an increase in agriculture that enabled its hosts', humans and domestic animals, densities to increase. Alternatively, the Thai populations could be a derivative of an older expansion event in the more northern populations. Mitochondrial DNA genealogy revealed no genetically divergent lineages, which agrees with the previous cytogenetic study. Genetic structure analyses found that only 27% of the pairwise comparisons were significantly different. The most likely explanation for the pattern of genetic structuring is the effect of genetic drift because of recent colonization.     

Numts help to reconstruct the demographic history of the ocellated lizard (Lacerta lepida) in a secondary contact zone

In northwestern Iberia, two largely allopatric Lacerta lepida mitochondrial lineages occur, L5 occurring to the south of Douro River and L3 to the north, with a zone of putative secondary contact in the region of the Douro River valley. Cytochrome b sequence chromatograms with polymorphisms at nucleotide sites diagnostic for the two lineages were detected in individuals in the region of the Douro River and further north within the range of L3. We show that these polymorphisms are caused by the presence of four different numts (I-IV) co-occurring with the L3 genome, together with low levels of heteroplasmy. Two of the numts (I and II) are similar to the mitochondrial genome of L5 but are quite divergent from the mitochondrial genome of L3 where they occur. We show that these numts are derived from the mitochondrial genome of L5 and were incorporated in L3 through hybridization at the time of secondary contact between the lineages. The additional incidence of these numts to the north of the putative contact zone is consistent with an earlier postglacial northward range expansion of L5, preceding that of L3. We show that genetic exchange between the lineages responsible for the origin of these numts in L3 after secondary contact occurred prior to, or coincident with, the northward expansion of L3. This study shows that, in the context of phylogeographic analysis, numts can provide evidence for past demographic events and can be useful tools for the reconstruction of complex evolutionary histories.     

The danger within: the role of genetic,behavioural and ecological factors in population persistence of colour polymorphic species

Polymorphic species have been the focus of important work in evolutionary biology. It has been suggested that colour polymorphic species have specific evolutionary and population dynamics that enable them to persist through environmental changes better than less variable species. We suggest that recent empirical and theoretical work indicates that polymorphic species may be more vulnerable to extinction than previously thought. This vulnerability arises because these species often have a number of correlated sexual, behavioural, life history and ecological traits, which can have a simple genetic underpinning. When exacerbated by environmental change, these alternate strategies can lead to conflict between morphs at the genomic and population levels, which can directly or indirectly affect population and evolutionary dynamics. In this perspective, we identify a number of ways in which the nature of the correlated traits, their underpinning genetic architecture, and the inevitable interactions between colour morphs can result in a reduction in population fitness. The principles illustrated here apply to all kinds of discrete polymorphism (e.g. behavioural syndromes), but we focus primarily on colour polymorphism because they are well studied. We urge further empirical investigation of the genetic architecture and interactions in polymorphic species to elucidate the impact on population fitness.     

Phylogeographic structure in the chromosomally polymorphic rodent Cricetulus barabensis sensu lato (Mammalia,Cricetidae)

Striped hamsters (Cricetulus barabensis sensu lato) represent a complex of chromosomally distinct allopatric lineages/taxa of either species or subspecies rank. They are widely distributed across the steppes of eastern and central Palearctic. Phylogenetic analysis of cytochrome b gene sequences based on 496 specimens from 112 localities revealed five well‐supported lineages divergent at 2%–4%. Two of them correspond to “griseus” (2n = 22) and “pseudogriseus” (2n = 24) karyomorphs and are placed as sister taxa. The “barabensis” (2n = 20) karyomorph is represented by three other branches and appears non‐monophyletic. All mtDNA lineages are distributed allopatrically or parapatrically; no indications of gene flow between populations of different chromosomal races were found. The results of the molecular clock analysis suggest that the main lineages diverged in the late Middle Pleistocene. The inferred evolutionary scenario implies that the common ancestor of the recent lineages belonged to the 2n = 20 karyomorph and originated in the eastern part of the contemporary range.     

Global population structure and demographic history of the grey seal

Although the grey seal Halichoerus grypus is one of the most familiar and intensively studied of all pinniped species, its global population structure remains to be elucidated. Little is also known about how the species as a whole may have historically responded to climate‐driven changes in habitat availability and anthropogenic exploitation. We therefore analysed samples from over 1500 individuals collected from 22 colonies spanning the Western and Eastern Atlantic and the Baltic Sea regions, represented by 350 bp of the mitochondrial hypervariable region and up to nine microsatellites. Strong population structure was observed at both types of marker, and highly asymmetrical patterns of gene flow were also inferred, with the Orkney Islands being identified as a source of emigrants to other areas in the Eastern Atlantic. The Baltic and Eastern Atlantic regions were estimated to have diverged a little over 10 000 years ago, consistent with the last proposed isolation of the Baltic Sea. Approximate Bayesian computation also identified genetic signals consistent with postglacial population expansion across much of the species range, suggesting that grey seals are highly responsive to changes in habitat availability.     

Genetic structure is correlated with phenotypic divergence rather than geographic isolation in the highly polymorphic strawberry poison‐dart frog

Phenotypic and genetic divergence can be influenced by a variety of factors, including sexual and natural selection, genetic drift and geographic isolation. Investigating the roles of these factors in natural systems can provide insight into the relative influences of allopatric and ecological modes of biological diversification in nature. The strawberry poison frog, Dendrobates pumilio, presents an excellent opportunity for this kind of research, displaying a diverse array of colour morphs and inhabiting a heterogeneous landscape that includes oceanic islands, fragmented rainforest patches and wide expanses of suitable habitat. In this study, we use 15 highly polymorphic microsatellite loci to estimate population structure and gene flow among populations from across the range of D. pumilio and a causal modelling framework to statistically test 12 hypotheses regarding the geographic and phenotypic variables that explain genetic differentiation within this system. Our results demonstrate that the genetic distance between populations is most strongly associated with differences in dorsal coloration. Previous experimental studies have shown that phenotypic differences can result in sexual and natural selection against non‐native phenotypes, and our results now show that these forces lead to genetic isolation between different colour morphs in the wild, presenting a potential case of incipient speciation through selection.     

Phylogenomics at the tips: inferring lineages and their demographic history in a tropical lizard,Carlia amax

High‐throughput sequencing approaches offer opportunities to better understand the evolutionary processes driving diversification, particularly in nonmodel organisms. In particular, the 100–1000's of loci that can now be sequenced are providing unprecedented power in population, speciation and phylogenetic studies. Here, we apply an exon capture approach to generate >99% complete sequence and SNP data across >2000 loci from a tropical skink, Carlia amax, and exploit these data to identify divergent lineages and infer their relationships and demographic histories. This is especially relevant to low‐dispersal tropical taxa that often have cryptic diversity and spatially dynamic histories. For C. amax, clustering of nuclear SNPs and coalescent‐based species delimitation analyses identify four divergent lineages, one fewer than predicted based on geographically coherent mtDNA clades (>9.4% sequence divergence). Three of these lineages are widespread and parapatric on the mainland, whereas the most divergent is restricted to islands off the northeast Northern Territory. Tests for population expansion reject an equilibrium isolation‐by‐distance model for two of the three widespread lineages and infer refugial expansion sources in the relatively mesic northeast Top End and northwest Kimberley. The latter is already recognized as a hotspot of endemism, but our results also suggest that a stronger focus on the northeast Top End, and adjacent islands is warranted. More generally, our results show how genome‐reduction methods such as exon capture can yield insights into the pattern and dynamics of biodiversity across complex landscapes with as yet poorly understood biogeographic history and how exon data can link between population and phylogenetic questions.     

Sexual and geographical variation in life history parameters of the shorthorn sculpin

A total of 293 shorthorn sculpins Myoxocephalus scorpius from Tromsø, northern Norway, were sampled between November 1998 and April 1999 to determine sex, total length, age, growth, maturity and mortality. Females grew to larger sizes ( L _∞=26·9 v. 18·5 cm), matured later (2 v. 1 year of age) at larger size (maturation length=16 v. 14 cm L _T), and had lower instantaneous mortality rates (0·93 v. 1·20 year⁻¹) than males. The life history parameters of shorthorn sculpins in northern Norway were more similar to the parameters of short-lived central European populations than to the parameters of the long-lived population of Newfoundland. This study confirms that northern Norwegian shorthorn sculpins exhibit sexual dimorphism as in other shorthorn sculpin populations. The relationships between growth pattern, age at maturity and mortality rates observed in the Tromsø population and in other shorthorn sculpin populations, correspond well with the predictions from a published life history model.     

Geographic and clonal variation in the milkweed-oleander aphid,Aphis nerii (Homoptera: Aphididae), for winged morph production,life history,and morphology in relation to host plant permanence

Summary Populations of the milkweed-oleander aphid,Aphis nerii, were sampled in California, Iowa and Puerto Rico. Among these localities the aphid's host plants differ greatly in permanence. I compared populations for migratory potential, measured as the proportion of winged offspring produced in response to being crowded, and for life history and morphometric traits of the subsequent adult winged aphids. I predicted a negative correlation between degree of host plant permanence and migratory potential. As predicted, aphids from Iowa, where migration on to temporary hosts must occur each year, produce a greater proportion of winged offspring (37.7%) than those from California (25.7%) or Puerto Rico (31.6%) where hosts are more permanent. However, hosts in Puerto Rico appear to be more permanent than those in California, yet the difference between populations for migratory potential was opposite to that predicted. Within California the prediction again held: aphids collected from the most impermanent sites produce the greatest proportion of winged offspring. There were no population differences for any life history or morphometric traits of winged aphids that are important contributors to fitness or migratory ability such as time to reproductive maturity, fecundity or wing length. Nor did any traits covary with migratory potential. Thus, there does not appear to be an association of life history and morphology with migratory potential that could enhance the colonizing ability of migrant aphids. I was unable to detect population differentiation for life history and morphology even though there is ample genetic variation within populations on which selection could act and an absence of constraints arising from genetic correlations that could prevent appropriate evolution of traits within populations. The exploitation of temporary host plants therefore occurs by an increase in the number of colonists produced and not by change in life history or morphology of those colonists.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Immigration history and gene dispersal: allozyme variation in Nordic populations of the red campion, Silene dioica (Caryophyllaceae)

Most of the Nordic region was ice-covered during the last (Weichselian) glaciation. During the postglacial period, plant and animal species recolonized the region from several directions and the geographic structuring of genetic variation within Nordic species may still contain a historic component that reflects patterns of postglacial immigration. The present investigation of 69 populations of Silene dioica represents the first large-scale allozyme study of a widespread herbaceous plant in the Nordic region. Although the frequencies of individual alleles showed a range of different geographic patterns, mapping of the axis scores from an ordination of variation at eight polymorphic loci revealed a division into two main geographic groups of populations. The broadly south-western and north-eastern distributions of these two groups of populations suggest that immigration into the region may have involved both eastern and southern geographic sources. However, the geographic boundaries between the two groups of populations are diffuse, and the relatively low between-population component of genetic diversity (G_ST = 16.4%) suggests a history of extensive gene dispersal by pollen.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 77 , 23–34.     

Population structure,connectivity, and demographic history of an apex marine predator,the bull shark Carcharhinus leucas

Knowledge of population structure, connectivity, and effective population size remains limited for many marine apex predators, including the bull shark Carcharhinus leucas. This large‐bodied coastal shark is distributed worldwide in warm temperate and tropical waters, and uses estuaries and rivers as nurseries. As an apex predator, the bull shark likely plays a vital ecological role within marine food webs, but is at risk due to inshore habitat degradation and various fishing pressures. We investigated the bull shark's global population structure and demographic history by analyzing the genetic diversity of 370 individuals from 11 different locations using 25 microsatellite loci and three mitochondrial genes (CR, nd4, and cytb). Both types of markers revealed clustering between sharks from the Western Atlantic and those from the Western Pacific and the Western Indian Ocean, with no contemporary gene flow. Microsatellite data suggested low differentiation between the Western Indian Ocean and the Western Pacific, but substantial differentiation was found using mitochondrial DNA. Integrating information from both types of markers and using Bayesian computation with a random forest procedure (ABC‐RF), this discordance was found to be due to a complete lack of contemporary gene flow. High genetic connectivity was found both within the Western Indian Ocean and within the Western Pacific. In conclusion, these results suggest important structuring of bull shark populations globally with important gene flow occurring along coastlines, highlighting the need for management and conservation plans on regional scales rather than oceanic basin scale.     

Association of Mc1r variants with ecologically relevant phenotypes in the European ocellated lizard, Lacerta lepida

A comprehensive knowledge on the genetic basis of coloration is crucial to understand how new colour phenotypes arise and how they contribute to the emergence of new species. Variation in melanocortin‐1 receptor (Mc1r), a gene that has been reported as a target for repeated evolution in a wide range of vertebrate taxa, was assessed in European ocellated lizards (Lacerta lepida) to search for associations with melanin‐based colour phenotypes. Lacerta lepida subspecies’ distribution is associated with the three major bio‐climatic regions in the Iberian Peninsula. A nonconserved and derived substitution (T162I) was associated with the L. l. nevadensis phenotype (prevalence of brown scales). Another substitution (S172C) was associated with the presence of black scales in both L. l. lepida and L. l. iberica, but no mutations were found to be associated with the higher proportion of black in L. l. iberica. Extensive genotyping of Mc1r along the contact zone between L. l. nevadensis and L. l. lepida revealed low gene flow (only two hybrids detected). The implications of these findings are discussed in the context of previous knowledge about the evolutionary history of ocellated lizards.     

Barcoding,population structure,and demographic history of Prodiplosis longifila associated with the Andes

Prodiplosis longifila Gagné (Diptera: Cecidomyiidae) is an insect pest that attacks various types of crops, including tomato, Solanum lycopersicum L. (Solanaceae), a vegetable with substantial economic significance worldwide. Prodiplosis longifila is a widely distributed pest in Colombia, Ecuador, and Peru, countries characterized by the presence of significant geographic barriers like the Andes Mountains. It has been reported that geographic barriers affect the dynamics and genetic differentiation of insect populations. Therefore, the aim of this study was to assess the diversity, genetic structure, and demographic history of P. longifila through the analysis of sequences within the mitochondrial region of cytochrome oxidase I (COI) and rDNA‐ITS2 in 27 populations located in Colombia and Ecuador. Analyses were performed on populations distributed in three geographic groups separated by the presence of the Andes Mountains. A total of 11 haplotypes were identified with the COI gene and only one haplotype in the rDNA‐ITS2 was found. Analyses of population structure and demographic history revealed that there is a structure associated with the Andes, which is reflected in an uneven distribution of the haplotype frequencies between regions, but even so, gene flow between populations was detected which produces low genetic differentiation. Because P. longifila has a short‐range dispersion that determines its territorial nature, it would be expected that other factors are producing the genetic exchange between populations. We suggest that the anthropogenic effect produced by farming practices, such as the use of seedlings as seed, which may carry P. longifila larvae, cause passive dispersal of pest throughout the Andes, particularly in Colombia.     

Geographic variation in sexual selection among populations of an iguanid lizard,Sauromalus obesus (=ater)

Abstract.— Geographic variation in selection pressures may result in population divergence and speciation, especially if sexual selection varies among populations. Yet spatial variation in targets and intensity of sexual selection is well studied in only a few species. Even more rare are simultaneous studies of multiple populations combining observations from natural settings with controlled behavioral experiments. We investigated how sexual selection varies among populations of the chuckwalla, Sauromalus obesus. Chuckwallas are sexually dimorphic in color, and males vary in coloration among populations. Using field observations and multiple regression techniques, we investigated how sexual selection acts on various male traits in three populations in which males differed in coloration. The influence of sexual selection on male coloration was then investigated in more detail using controlled experiments. Results from field observations indicate that phenotypic selection was acting on territory quality in all three populations. In two populations, selection was also acting either directly or indirectly on male coloration. Male color likely functions as an indicator of food resources to females because male color is based partly on carotenoid pigments. In controlled experiments, significantly more females from these two populations chose males with brighter colors over dull males, a result consistent with studies on carotenoid pigments in other taxa. In a third population, no evidence of sexual selection on male coloration was found in either the field study or controlled experiment. Lack of female preferences for male color in this population, in which chuckwalla densities are low and home ranges are large, may result from searching costs to females.     

间颅鼠兔的种群遗传结构和历史演化动态

本研究基于多个线粒体基因（包括COI、Cyt b、tRNA-Thr、tRNA-Pro、D-loop控制区部分序列共计2417bp）对间颅鼠兔（Ochotona cansus）的种群遗传结构和历史演化动态进行了研究。基于线粒体基因联合数据构建的贝叶斯树显示,间颅鼠兔分为两大支系：来自岷山山系东部的甘肃白水江和四川唐家河的种群组成了支系一;其余采样点的种群组成了支系二。基于以上基因所构建的单倍型网络图显示的遗传结构与贝叶斯树一致。贝叶斯聚类分析（BAPS）将间颅鼠兔分为4个地理种群。A种群包括甘肃白水江（BSJ）和四川唐家河（TJH）的样本;B种群仅包括四川卧龙（WL）的样本;C种群包括青海北部种群（GC、MY、QL）、青海南部-四川北部种群（JW、JZ、ZK、RE）、四川康定（KD）和陕西太白山（TB）的样本;D种群包括山西娄烦（LF）以及陕西吴起（WQ）的样本。基于线粒体基因构建的物种树进行支系分化时间的推测,结果显示4个地理种群的分化时间约为0.72 Ma、0.32 Ma、0.12 Ma, 这与造成青藏高原剧烈隆起的昆黄运动（1.10-0.60Ma）和共和运动（0.15Ma）发生时间基本吻合。采用中性检验、错配分布、扩展贝叶斯天际线（EBSPs）3种方法对间颅鼠兔的种群历史动态进行预测,分析结果显示间颅鼠兔在历史上并没有经历过显著的种群扩张,而且在末次盛冰期后开始衰减。生态位模型（ENM）预测结果表明：末次盛冰期时间颅鼠兔的潜在分布区比其现今分布区小,末次盛冰期比末次间冰期潜在分布区略大,表明间颅鼠兔受末次盛冰期影响较大。气候因子分析显示降水比温度对间颅鼠兔种群影响更大,推测降水对植被的演化影响显著,从而影响间颅鼠兔的食物来源和栖息环境。     

Genetic variation in the life‐history traits of Epiphyas postvittana: population structure and local adaptation

The light brown apple moth, Epiphyas postvittana (Walker) shows high intraspecific variability in morphological, physiological, demographic and behavioural characters. To gain insight into the extent of adaptation and evolutionary changes in response to environmental heterogeneity in this species, quantitative genetic analyses of life‐history variation were conducted for two natural populations under two thermal conditions (23°C and 28°C). Paternal half‐sib heritability and genetic correlation in six life‐history traits (i.e. development time, adult body weight, adult lifespan, age at first reproduction, the number of eggs laid during the first 5 days after emergence, and total fecundity) were compared. Significant heritabilities were shown consistently in development time; this is also true for adult body weight, except for the Canberra population at 23°C. However, neither population differences nor the effect of temperature were statistically detectable for any of these heritabilities, confirming the genetically determined flexibility. Positive genetic correlations between development time and adult body weight, and negative genetic correlations between the number of eggs laid during the first 5 days and adult lifespan were present for these populations at both temperatures, indicating the presence of genetic constraints. Pairwise comparisons of genetic correlations revealed the heterogeneity of the two populations and across temperatures. These results suggest that the structure of genetic covariance might have changed significantly during the divergence of natural populations and in response to the alteration of environmental conditions in E. postvittana.     

Population genetic structure and demographic history of the endemic Formosan lesser horseshoe bat (Rhinolophus monoceros)

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. Taiwan formed around 5 million years ago from tectonic uplift, and has been connected to mainland Asia several times since its emergence. A central mountain range runs north to south, bisecting the island, and potentially impedes gene flow along an east-west axis. The Formosan lesser horseshoe bat (Rhinolophus monoceros) is endemic to Taiwan, where it is found mainly at low altitude. To determine the population structure and the demographic and colonization history of this species, we examined variation in the mitochondrial DNA control region in 203 bats sampled at 26 sites. We found very high haplotype and nucleotide diversity, which decreased from the centre to the south and north. Population differentiation followed a pattern of isolation by distance, though most regional genetic variance was attributable to differences between the relatively isolated southern population and those from other regions. A haplotype network was consistent with these findings and also suggested a southward colonization, followed by subsequent secondary contact between the south and other regions. Mismatch distributions were used to infer a past population expansion predating the last glacial maximum, and a neighbour-joining tree showed that R. monoceros formed a monophyletic grouping with respect to its sister taxa. Taken together, our results suggest that this taxon arose from a single period of colonization, and that demographic growth followed in the late Pleistocene. Current genetic structure reflects limited gene flow, probably coupled with stepwise colonization in the past. We consider explanations for the persistence of the species through multiple glacial maxima.     

Geographic variation in the structure of oak hybrid zones provides insights into the dynamics of speciation

Studying geographic variation in the rate of hybridization between closely related species could provide a useful window on the evolution of reproductive isolation. Reinforcement theory predicts greater prezygotic isolation in areas of prolonged contact between recently diverged species than in areas of recent contact, which implies that old contact zones would be dominated by parental phenotypes with few hybrids (bimodal hybrid zones), whereas recent contact zones would be characterized by hybrid swarms (unimodal hybrid zones). Here, we investigate how the hybrid zones of two closely related Chinese oaks, Quercus mongolica and Q. liaotungensis, are structured geographically using both nuclear and chloroplast markers. We found that populations of Q. liaotungensis located around the Changbai Mountains in Northeast China, an inferred glacial refugium, were introgressed by genes from Q. mongolica, suggesting historical contact between the two species in this region. However, these introgressed populations form sharp bimodal hybrid zones with Q. mongolica. In contrast, populations of Q. liaotungensis located in North China, which show no sign of ancient introgression with Q. mongolica, form unimodal hybrid zones with Q. mongolica. These results are consistent with the hypothesis that selection against hybrids has had sufficient time to reinforce the reproductive barriers between Q. liaotungensis and Q. mongolica in Northeast China but not in North China.