Impact of mountain chains, sea straits and peripheral populations on genetic and taxonomic structure of a freshwater turtle, Mauremys leprosa (Reptilia, Testudines, Geoemydidae)

Mauremys leprosa, distributed in Iberia and North‐west Africa, contains two major clades of mtDNA haplotypes. Clade A occurs in Portugal, Spain and Morocco north of the Atlas Mountains. Clade B occurs south of the Atlas Mountains in Morocco and north of the Atlas Mountains in eastern Algeria and Tunisia. However, we recorded a single individual containing a clade B haplotype in Morocco from north of the Atlas Mountains. This could indicate gene flow between both clades. The phylogenetically most distinct clade A haplotypes are confined to Morocco, suggesting both clades originated in North Africa. Extensive diversity within clade A in south‐western Iberia argues for a glacial refuge located there. Other regions of the Iberian Peninsula, displaying distinctly lower haplotype diversities, were recolonized from within south‐western Iberia. Most populations in Portugal, Spain and northern Morocco contain the most common clade A haplotype, indicating dispersal from the south‐western Iberian refuge, gene flow across the Strait of Gibraltar, and reinvasion of Morocco by terrapins originating in south‐western Iberia. This hypothesis is consistent with demographic analyses, suggesting rapid clade A population increase while clade B is represented by stationary, fragmented populations. We recommend the eight, morphologically weakly diagnosable, subspecies of M. leprosa be reduced to two, reflecting major mtDNA clades: Mauremys l. leprosa (Iberian Peninsula and northern Morocco) and M. l. saharica (southern Morocco, eastern Algeria and Tunisia). Peripheral populations could play an important role in evolution of M. leprosa because we found endemic haplotypes in populations along the northern and southern range borders. Previous investigations in another western Palearctic freshwater turtle (Emys orbicularis) discovered similar differentiation of peripheral populations, and phylogeographies of Emys orbicularis and Mauremys rivulata underline the barrier status of mountain chains, in contrast to sea straits, suggesting common patterns for western Palearctic freshwater turtles.     

Phenotypic variation among local populations of phlebotomine sand flies (Diptera: Psychodidae) in southern Turkey.

The wing-shape morphology of local populations of the medically important phlebotomine sand flies, Phlebotomus sergenti, P. papatasi, P. tobbi, and P. similis, were examined in both sexes by using geometric morphometrics. There are three major mountain ranges that may serve as geographical barriers for species distribution in the study area and four main gaps were recognized among these barriers. We found no statistically important differences in wing morphology in all examined species in both sexes for all local populations. These results show that the barriers are not sufficient to stop gene flow among local populations of sand flies. The graphical depiction of PCA, CVA, and F-test confirmed our morphometric study suggesting that the difference in wing morphology between P. similis and P. sergenti indicates that these are clearly different species. These two show sympatric distribution in the Konya Plain of Anatolia.     

Mitochondrial differentiation and biogeography of Hyla meridionalis (Anura: Hylidae): an unusual phylogeographical pattern

Aim To study the patterns of genetic variation and the historical events and processes that influenced the distribution and intraspecific diversity in Hyla meridionalis Boettger, 1874. Location Hyla meridionalis is restricted to the western part of the Mediterranean region. In northern Africa it is present in Tunisia, Algeria and Morocco. In south‐western Europe it is found in the south of France, north‐western Italy and north‐eastern and south‐western Iberian Peninsula. There are also insular populations, as in the Canaries and Menorca. Methods Sampling included 112 individuals from 36 populations covering the range of the species. We used sequences of mitochondrial DNA Cytochrome Oxidase I (COI) for the phylogeographical analysis (841 bp) and COI plus a fragment including part of tRNA lysine, ATP synthase subunits 6 and 8 and part of Cytochrome Oxidase III for phylogenetic analyses (2441 bp). Phylogenetic analyses were performed with paup *4.0b10 (maximum likelihood, maximum parsimony) and Mr Bayes 3.0 (Bayesian analysis). Nested clade analysis was performed using tcs 1.18 and Geo Dis 2.2. A dispersal‐vicariant analysis was performed with diva 1.0 to generate hypotheses about the geographical distribution of ancestors. Results We found little genetic diversity within samples from Morocco, south‐western Europe and the Canary Islands, with three well‐differentiated clades. One is distributed in south‐western Iberia and the High Atlas, Anti‐Atlas and Massa River in Morocco. The second is restricted to the Medium Atlas Mountains. The third one is present in northern Morocco, north‐eastern Iberia, southern France and the Canaries. These three groups are also represented in the nested clade analysis. Sequences from Tunisian specimens are highly divergent from sequences of all other populations, suggesting that the split between the two lineages is ancient. diva analysis suggests that the ancestral distribution of the different lineages was restricted to Africa, and that an explanation of current distribution of the species requires three different dispersal events. Main conclusions Our results support the idea of a very recent colonization of south‐western Europe and the Canary Islands from Morocco. South‐western Europe has been colonized at least twice: once from northern Morocco probably to the Mediterranean coast of France and once from the western coast of Morocco to southern Iberia. Human transport is a likely explanation for at least one of these events. Within Morocco, the pattern of diversity is consistent with a model of mountain refugia during hyperarid periods within the Pleistocene. Evaluation of the phylogenetic relationships of Tunisian haplotypes will require an approach involving the other related hylid taxa in the area.     

Genetic signals of ancient decline in Aleppo pine populations at the species' southwestern margins in the Mediterranean Basin

Microsatellite markers were used to characterize the structure of genetic diversity in natural Moroccan Aleppo pine (Pinus halepensis Mill.) populations, the most southwesterly marginal populations of the species in the Mediterranean Basin. Twenty-two natural populations and one artificial population, located in four regions covering most of the natural range of P. halepensis in the country, were sampled. Across this range, towards the south and west (and towards high altitudes) the populations become increasingly discrete and discontinuous. The nuclear microsatellite marker analysis suggests that a large proportion of the Aleppo pines in Morocco have derived from a single genetic lineage, represented by a central group of 11 of the examined populations located in the High and Middle Atlas Mountains. In addition, two smaller groups, represented by the marginal southwestern High Atlas populations, and three still smaller north / northeastern groups of populations located in the Rif and northeast Middle Atlas Mountains, could be genetically distinguished. Further, coalescence analysis of historical demographic population patterns suggests that ancient bottlenecks occurred in all of the natural populations. However, the population differentiation and genetic diversity levels we found were good (F(st) =15.47), presumably because of the species' good potential for long-distance dispersal of seeds and high invasive capacity, which appear to have maintained a state of stable near-equilibrium, meta-population dynamics since ancient times.     

Phylogeographic patterns in Hypochaeris section Hypochaeris (Asteraceae, Lactuceae) of the western Mediterranean

Aim  To analyse phylogeographic patterns in the four species of Hypochaeris sect. Hypochaeris , evaluating possible areas of origin and the microevolutionary processes that have shaped their morphology, genetics and distribution.
Location  Western Mediterranean area.
Methods  We applied amplified fragment length polymorphism (AFLP) markers to a total of 494 individuals belonging to 82 populations of Hypochaeris arachnoidea , H. glabra , H. radicata and H. salzmanniana to determine population structure.
Results  Populations with the largest proportion of private and rare AFLP fragments were found in Morocco. This region was consequently inferred to be the ancestral area for H. arachnoidea , H. glabra , H. radicata and H. salzmanniana . The Guadalquivir River (southern Spain) was inferred to be an effective dispersal barrier for H. glabra and H. radicata. The Strait of Gibraltar was inferred to be a somewhat weaker barrier than the Guadalquivir River for H. radicata and a much weaker barrier for H. glabra . The main barrier for H. salzmanniana coincides with the extension of the Rif Mountains to the Atlantic coast in Morocco, and the Strait of Gibraltar is a much weaker barrier for this species. Hypochaeris arachnoidea appears to have originated in the Atlas Mountains.
Main conclusions  The highest levels of genetic variation in La Mamora forest ( H. glabra and H. salzmanniana ) or the adjacent central Middle Atlas ( H. arachnoidea and H. radicata ) in Morocco suggest that these areas were a centre of origin of Hypochaeris sect. Hypochaeris . All three potential barriers – the Guadalquivir River, the Strait of Gibraltar, and the Rif Mountains – have been important in shaping genetic diversity in species of section Hypochaeris .     

Otolith morphology varies between populations,sexes and male alternative reproductive tactics in a vocal toadfish Porichthys notatus

In this study, the morphology of sagittal otoliths of the plainfin midshipman fish Porichthys notatus was compared between populations, sexes and male alternative reproductive phenotypes (known as ‘type I males or guarders’ and ‘type II males or sneakers’). Sagitta size increased with P. notatus size and changes in shape were also detected with increasing body size. Porichthys notatus sagittae begin as simple rounded structures, but then elongate as they grow and take on a more triangular and complex shape with several prominent notches and indentations along the dorsal and caudal edges. Moreover, the sagittae of the two geographically and genetically distinct populations of P. notatus (northern and southern) differed in shape. Porichthys notatus from the north possessed taller sagittae with deeper caudal indentations compared to P. notatus from the south. Sagitta shape also differed between females and males of the conventional guarder tactic. Furthermore, guarder males had smaller sagittae for their body size than did sneaker males or females. These differences in sagittal otolith morphology are discussed in relation to ecological and life history differences between the sexes and male tactics of this species. This is the first study to investigate teleost otolith morphology from the perspective of alternative reproductive tactics.     

Intraspecific Sexual Size and Shape Dimorphism in an Australian Freshwater Fish Differs with Respect to a Biogeographic Barrier and Latitude

Geographically structured variation in morphology is a common phenomenon in animals with environmental factors covarying with both latitude and biogeographic barriers having profound impacts on body size and shape. The Pacific blue-eye (Pseudomugil signifer) is a freshwater fish that lives along Australia’s east coast and occurs on either side of a terrestrial barrier, the Burdekin Gap. By quantifying the size and shape of males and females from 10 populations we found that Pacific blue-eyes are not sexually size dimorphism north of the Burdekin Gap whereas the degree of dimorphism was dependent upon latitude south of the barrier. Rensch’s rule was not supported as the degree of male-biased size dimorphism did not increase with increasing population mean body size. Body shape was related to body size and was sexually dimorphic south of the Burdekin Gap but not north of it. Our study represents a rare case of identifying how both body size and shape differ with respect to latitude and a major terrestrial biogeographic barrier and lends further support to the notion that P. signifer may comprise two species, or incipient species, that are separated by the Burdekin Gap.     

Description of the Nymphs and Eggs of Acentrella almohades sp. n. from Morocco and Southern Spain (Ephemeroptera: Baetidae)

The nymph and egg of Acentrella almohades sp. n. are described and illustrated on the basis of material collected in southern Spain (in the upper reaches of the Guadalquivir River Basin) and Morocco (Rif Mountains and Atlas Mountains). Features distinguishing the new species from other European species are discussed and keyed.     

Cytisus purgans auct(Leguminosae-Papilionoideae comprises four distinct species)

Based on a study of the morphology and stem anatomy of a broad sample of specimens attributed to Cytisus purgans auct. non Cytisus purgans (L.) Boiss. four species are recognized. Two species are non-spiny: Cytisus oromediterraneus Rivas Martínez et at. is a widely distributed taxon ranging from central France to north-central Spain and Portugal, and C. galianoi Talavera & Gibbs, sp. nov., which is restricted to schistose sierras in south-east Spain. C. oromediterraneus has 8-ribbed stems and 2-ribbed stipular organs, whilst C. galianoi has 10-ribbed stems and 3-ribbed stipular organs. Of two spiny species, Cytisus balansae (Boiss.) Ball wkh 10-ribbed stems, occurs in calcareous mountain areas, 1800–2500 m in Algeria and southern Morocco, whilst C. valdesii Talavera & Gibbs, nom. nov. with 13-ribbed stems, occurs in die High and Middle Atlas of Morocco. Both north African species have 3-ribbed stipular organs.     

Randomly Amplified Polymorphic DNA Variation in Populations of Eastern Australian Koalas, Phascolarctos cinereus

Randomly amplified polymorphic DNA (RAPD) variation in populations of the koala, Phascolarctos cinereus, was investigated, revealing significant differences in the level of diversity between southern and northern regions of eastern Australia. Of the 20 polymorphic RAPD markers identified in koalas, 4-7 were polymorphic in southern populations, while 12-17 were polymorphic in northern populations. Analysis of molecular variance revealed a significant difference in the estimated variance between koalas from northern and those from southern regions (P < 0.001), where populations from the north were greater than twice as variable as their southern cousins. The total genetic diversity observed was attributed to regional differences (30.91%), population differences within a region (11.77%), and differences among individuals within a population (57.32%). For the within-region analyses, a large proportion of the genetic diversity was attributable to individual differences within a population, 80.34% for the north and 91.23% for the south. These results demonstrate that RAPD markers are useful for determining population structure among koalas.     

Cryptic differentiation and geographic variation in genetic diversity of Hall's Babbler Pomatostomus halli

Sequence variation was examined in domain I of the mitochondrial control region in three Queensland populations of Hall's Babbler Pomatostomus halli , a geographically restricted, monotypic songbird in eastern Australia. Surprisingly, we found that domain I sequences were strongly differentiated into two major clades differing by 3.29%. These two clades exhibited nearly complete geographic concordance with northern and southern populations, except for two haplotypes which were sampled in the north of the range but were phylogenetically allied to the southern clade. We also found a seven-fold higher level of genetic diversity in the northern than in the southern populations. Neutrality and molecular clock tests suggested that selection or differences in substitution rates were not responsible for this difference in diversity. However, a maximum likelihood analysis of gene flow between the north and south suggested that the difference in diversity could be due to both greater population size in the north and asymmetric gene flow dominated by south to north dispersal events. A likelihood ratio test rejected a model in which population sizes were equal and rates of gene flow symmetric, and came close to rejecting a model in which only population sizes were constrained to be equal. These results suggest that different population sizes and asymmetric gene flow could be a major source of differences in genetic variation between populations of Hall's Babbler, although ecological and biogeographic causes for these differences are obscure.     

Spatial distribution and breeding territories of Bonelli’s Eagle Aquila fasciata in the western Anti-Atlas Mountains,Morocco

The North African population of Bonelli’s Eagle Aquila fasciata (Vieillot, 1822) is limited to the south by the northern fringe of the Sahara Desert. This study provides the first data on the spatial distribution and density of breeding Bonelli’s Eagles in south-west Morocco, at the southern limit of their breeding distribution in the Western Palearctic. We used broad-scale sampling to investigate spatial patterns in occupancy of territorial pairs of Bonelli’s Eagles in an area of ～29 715 km2 in the western Anti-Atlas Mountains, southern Morocco, during 2016–2018. We found 28 nesting territories containing 40 used nests, heterogeneously distributed in areas of high topographic variation, from 60 to 1 890 m asl. The average nearest-neighbour nest distance was 14.12 ± 9.90 km and varied from 2.65 km in the north-west to 37.80 km in the pre-Saharan lands in the south-west portion of the study area. This work shows the importance of the western Anti-Atlas Mountains as one of the most significant strongholds of the species in Morocco. However, there is a need for systematic analyses of the different factors affecting the distribution of the species to implement conservation actions of this peripheral population.     

Historical biogeography of the Western Rattlesnake (Serpentes: viperidae: Crotalus viridis), inferred from mitochondrial DNA sequence information

We infer the phylogeography of the Western Rattlesnake (Crotalus viridis) using phylogenetic analysis of mitochondrial DNA sequences from 1345 bp of the genes for cytochrome b and NADH dehydrogenase subunit 4. Two main clades are revealed: one includes populations from east and south of the Rocky Mountains (conventionally referred to as Crotalus viridis viridis and C. v. nuntius), and the other consists of populations west of the Rocky Mountains. Within the western clade, a population from southern Arizona (C. v. cerberus) represents the sister taxon to the remaining western populations. The conventional subspecies recognized in this species do not fully correspond to the phylogenetic pattern, and a review of the systematic status of several populations is needed. Our data allow the inferences that small body size evolved twice and that the ability of one population (C. v. concolor) to secrete highly lethal toxins related to Mojave toxin arose within the complex. Our phylogeny should represent the basis for further studies on the causes of geographical variation in this complex.     

Phylogeography of North African Atlas cedar (Cedrus atlantica, Pinaceae): Combined molecular and fossil data reveal a complex Quaternary history

Northwest Africa is a major hotspot of plant biodiversity, but very little is known about the Quaternary range dynamics of plant species in this region. Here we investigate the range-wide population structure and phylogeography of Atlas cedar (Cedrus atlantica), an emblematic forest tree endemic to Morocco and Algeria. We genotyped 261 individuals from 11 populations using AFLP markers. Data were analyzed using both conventional F(ST)-based techniques and Bayesian clustering. Overall population differentiation was high (F(ST) = 0.25). Two major groups of populations were identified, one distributed through the Rif and Middle Atlas mountains in Morocco and the other through the Algerian Tell Atlas and Aurès mountains as well as the Middle Atlas. Combined molecular and fossil data indicate that C. atlantica survived the Last Glacial Maximum in at least three disjunct refugia along the coast of the Mediterranean Sea, whereas the Middle Atlas, today the core of the species range, has been colonized relatively recently (<10000 yr BP). The colonization history of individual populations has left clear imprints in their present-day diversity, which may vary greatly even between nearby stands. Our study illustrates how integrating different data sources and analytical approaches can help elucidate complex range dynamics that would otherwise remain undeciphered.     

The influence of landscape structure on occurrence, abundance and genetic diversity of the common frog, Rana temporaria

Human‐caused habitat destruction and modification constitute one of the largest threats to population persistence and biodiversity, and are also suspected to be the major cause behind the global decline of amphibian populations. We assessed the potential role of agriculture‐related habitat fragmentation on population size and genetic variability in the common frog (Rana temporaria) by recording the occurrence, population density and genetic diversity in three geographically disparate regions in Sweden – each containing landscapes of high and low agricultural activity – and related these to landscape variables extracted from digital maps. We found a highly significant region‐by‐landscape interaction in occurrence, population density and genetic diversity revealing a reversed response to agriculture from south to north: while the effects of agriculture on R. temporaria populations were negative in the south, there were no effects in the central region, whereas positive effects were observed in the north. Spatial autocorrelation analyses of genetic data revealed that populations in high agricultural activity areas were more isolated than populations in low activity areas both in the southern and central regions of Sweden. Landscape diversity showed a strong positive correlation with both density and occurrence of frogs in Sweden as a whole, as well as in the southern region. Also, negative effects of roads and positive effects of ditches on genetic diversity were found in the south. Overall, these results suggest clear but regionally opposite effects of habitat structure on the population size and genetic diversity of amphibian populations. This means that the management strategy aiming to maximize the size and genetic diversity of local common frog populations, and perhaps also those of other amphibian populations, should account for regional differences in existing land‐use patterns.     

Mixed support for an alignment between phenotypic plasticity and genetic differentiation in damselfly wing shape

The relationship between genetic differentiation and phenotypic plasticity can provide information on whether plasticity generally facilitates or hinders adaptation to environmental change. Here, we studied wing shape variation in a damselfly (Lestes sponsa) across a latitudinal gradient in Europe that differed in time constraints mediated by photoperiod and temperature. We reared damselflies from northern and southern populations in the laboratory using a reciprocal transplant experiment that simulated time-constrained (i.e. northern) and unconstrained (southern) photoperiods and temperatures. After emergence, adult wing shape was analysed using geometric morphometrics. Wings from individuals in the northern and southern populations differed significantly in shape when animals were reared in their respective native environment. Comparing wing shape across environments, we found evidence for phenotypic plasticity in wing shape, and this response differed across populations (i.e. G × E interactions). This interaction was driven by a stronger plastic response by individuals from the northern population and differences in the direction of plastic wing shape changes among populations. The alignment between genetic and plastic responses depended on the specific combination of population and rearing environment. For example, there was an alignment between plasticity and genetic differentiation under time-constrained, but not under non-time-constrained conditions for forewings. We thus find mixed support for the hypothesis that environmental plasticity and genetic population differentiation are aligned. Furthermore, although our laboratory treatments mimicked the natural climatic conditions at northern and southern latitudes, the effects of population differences on wing shape were two to four times stronger than plastic effects. We discuss our results in terms of time constraints and the possibility that natural and sexual selection is acting differently on fore- and hindwings.     

Range-wide molecular analysis of the western pond turtle (Emys marmorata): cryptic variation, isolation by distance, and their conservation implications

We analysed phylogeography and population genetic variation across the range of the western pond turtle (Emys marmorata) using rapidly evolving mitochondrial and nuclear DNA sequence data. Nuclear DNA sequences from two unlinked introns displayed extremely low levels of variation, but phylogenetic analyses based on mtDNA recovered four well-supported and geographically coherent clades. These included a large Northern clade composed of populations from Washington south to San Luis Obispo County, California, west of the Coast Ranges; a San Joaquin Valley clade from the southern Great Central Valley; a geographically restricted Santa Barbara clade from a limited region in Santa Barbara and Ventura counties; and a Southern clade that occurs south of the Tehachapi Mountains and west of the Transverse Range south to Baja California, Mexico. An analysis of molecular variance (amova) based on regional hydrographic units revealed that populations from the Sacramento Valley north to Washington were virtually invariant, with no evidence of population substructure among northern river drainage basins. In other areas, E. marmorata contains considerable unrecognized variation, particularly in central and southern California and in northern Baja California, Mexico. Our northern clade is congruent with the distribution of the subspecies Emys marmorata marmorata (Washington-central California). However, no clade is congruent with the distribution of the southern subspecies Emys marmorata pallida from central California-Baja. Thus, recognition of the current subspecies split is not warranted, based on the available genetic evidence. Our amova and phylogenetic results, in conjunction with a growing comparative database for other codistributed aquatic taxa, confirm the occurrence of genetic breaks across the Tehachapi Mountains and Transverse Range bounding the southern end of the Great Central Valley, and point to southern California as a rich source of cryptic genetic variation.     

Population structure of Cicada barbara Stål (Hemiptera, Cicadoidea) from the Iberian Peninsula and Morocco based on mitochondrial DNA analysis

We assess the genetic history and population structure of Cicada barbara in Morocco and the Iberian Peninsula, based on analysis of the mitochondrial cytochrome b gene. The divergence between Morocco and the Iberian Peninsula populations was strongly corroborated by the molecular data, suggesting genetically isolated populations with a low level of gene flow. The Ceuta population from Spanish North Africa was more similar to the Iberian populations than the surrounding Moroccan populations, suggesting that the Strait of Gibraltar has not been acting as a strict barrier to dispersal while the Rif Mountains have. The Iberian Peninsula specimens showed a signature of demographic expansion before that which occurred in Morocco, but some of the assumptions related to the demographic parameters should be considered with caution due to the small genetic variation found. The high haplotype diversity found in Morocco implies higher demographic stability than in the Iberian Peninsula populations. These results do not, however, suggest a Moroccan origin for Iberian cicadas; but the most northwest region in Africa, such as Ceuta, might have acted as a southern refuge for Iberian cicadas during the most severe climatic conditions, from where they could expand north when climate improved. The separation of two subspecies within C. barbara (C. barbara lusitanica and C. barbara barbara) finds support with these results.     

拟环纹豹蛛不同地理种群遗传多样性的AFLP标记研究

为查明拟环纹豹蛛Pardosa pseudoannulata不同地理种群的遗传多样性机制,应用AFLP技术对6个拟环纹豹蛛地理种群的遗传多样性进行了研究分析。8对引物组合扩增出1 038个AFLP条带,其中多态性条带占86.622%,全部个体显示了各自独特的AFLP图谱。AFLP标记的遗传多样性分析结果表明： 拟环纹豹蛛无论在物种水平(P=86.62%,H=0.2622,I=0.3101),还是在种群水平（P=73.0%,H=0.2155,I=0.2554）都表现出较高的遗传多样性。其中湖南长沙雷锋镇种群内遗传变异最大,云南高黎贡山福贡种群内遗传变异最小,华南北部（湖南、湖北、江西）地区拟环纹豹蛛遗传多样性明显高于华南南部（云南、海南）种群。据种群变异来源分析,有35.77%的遗传变异来自种群间,64.23%的变异来源于种群内（N_m=0.898）,不同地理种群显示出一定的遗传分化。分析认为海拔是影响拟环纹豹蛛遗传分化的重要因素,这为进一步明确我国稻田狼蛛优势种群在农药胁迫下的遗传适应性机制提供了实验依据。     

油松遗传结构与地理阻隔因素的相关性

油松(Pinus tabulaeformis)是我国特有广布树种。秦岭是我国亚热带和暖温带气候的重要分界线,贺兰山是我国西北干旱区与半干旱区的分界线。旨在以分布于秦岭、贺兰山和晋陕两省的10个天然油松种群为研究对象,通过对其线粒体DNA(nad1和matR内含子)的序列,探讨地理阻隔对油松遗传结构的影响。结果表明,10个天然油松种群的100个个体中共检测到27个单倍型,多态位点172个,简约信息位点35个,单一多态性位点137个。其中1个单倍型为全部种群共有,4个单倍型为2-5个种群共有,其余22个单倍型为各个种群所独有。在这10个种群中,有8个种群分别具有1-4个独有单倍型。尽管秦岭南侧种群和北侧种群分别具有两个和一个独有单倍型; 贺兰山东侧种群和西侧种群均具有3个独有单倍型;但是与晋陕种群(即油松分布区中心种群)相比,其种群独有单倍型平均数目则明显较少。与此同时,不仅秦岭南北两侧或贺兰山东西两侧种群的独有单倍型的进化关系往往较近,而且秦岭北侧或贺兰山东侧种群的独有单倍型与油松分布区中心种群的某些独有单倍型的关系亦较近。从而导致秦岭南北两侧或贺兰山东西两侧种群与油松分布区中心种群之间表现出比较复杂的聚类关系。由此可见,油松的遗传结构与山脉屏障的存在没有明显关系。