Tracing species replacement in Iberian marbled newts

Secondary contact between closely related species can lead to the formation of hybrid zones, allowing for interspecific gene flow. Hybrid zone movement can take place if one of the species possesses a competitive advantage over the other, ultimately resulting in species replacement. Such hybrid zone displacement is predicted to leave a genomic footprint across the landscape in the form of asymmetric gene flow (or introgression) of selectively neutral alleles from the displaced to the advancing species. Hybrid zone movement has been suggested for marbled newts in the Iberian Peninsula, supported by asymmetric gene flow and a distribution relict (i.e., an enclave) of Triturus marmoratus in the range of T. pygmaeus. We developed a panel of nuclear and mitochondrial SNP markers to test for the presence of a T. marmoratus genomic footprint in the Lisbon peninsula, south of the enclave. We found no additional populations of T. marmoratus. Analysis with the software Structure showed no genetic traces of T. marmoratus in T. pygmaeus. A principal component analysis showed some variation within the local T. pygmaeus, but it is unclear if this represents introgression from T. marmoratus. The results may be explained by (a) species replacement without introgressive hybridization and (b) displacement with hybridization followed by the near‐complete erosion of the footprint by purifying selection. We predict that testing for a genomic footprint north of the reported enclave would confirm that species replacement in these marbled newts occurred with hybridization.     

Is dispersal guided by the environment? A comparison of interspecific gene flow estimates among differentiated regions of a newt hybrid zone

How frequently genes pass through a hybrid zone may be influenced by the environment. Accordingly, in long hybrid zones that span more than one environmental setting, different patterns may emerge. The varied conditions allow testing of hypotheses on dispersal as a function of the environment. We reconstruct the amount and direction of gene flow across a heterogeneous hybrid zone of two species of marbled newts (Triturus marmoratus and Triturus pygmaeus), in four widely separated areas of the Iberian Peninsula from one mitochondrial and three nuclear genes. The main variables associated with the position of the contact zone are precipitation, rivers, altitude and relief. In some sections of the contact zone, however, its position is not correlated with any environmental factor and is instead determined by the shortest geographical distance between fixed positions at either side (mountains in the East and river in the West). In areas where the position of the zone is stable, gene flow was bidirectional. External data show that T. pygmaeus has superseded T. marmoratus over a large area and here gene flow was unidirectional. The prediction that a major river would reduce gene flow was not confirmed.     

Stabilization of a salamander moving hybrid zone

When related species meet upon postglacial range expansion, hybrid zones are frequently formed. Theory predicts that such zones may move over the landscape until equilibrium conditions are reached. One hybrid zone observed to be moving in historical times (1950–1979) is that of the pond‐breeding salamanders Triturus cristatus and Triturus marmoratus in western France. We identified the ecological correlates of the species hybrid zone as elevation, forestation, and hedgerows favoring the more terrestrial T. marmoratus and pond density favoring the more aquatic T. cristatus. The past movement of the zone of ca. 30 km over three decades has probably been driven by the drastic postwar reduction of the “bocage” hedgerow landscape, favoring T. cristatus over T. marmoratus. No further hybrid zone movement was observed from 1979 to the present. To explain the changing dynamics of the hybrid zone, we propose that it stalled, either because an equilibrium was found at an altitude of ca. 140 m a.s.l. or due to pond loss and decreased population densities. While we cannot rule out the former explanation, we found support for the latter. Under agricultural intensification, ponds in the study area are lost at an unprecedented rate of 5.5% per year, so that remaining Triturus populations are increasingly isolated, hampering dispersal and further hybrid zone movement.     

Testing the hybrid superiority hypothesis in crested and marbled newts

We tested the hybrid superiority hypothesis in the zone of overlap and hybridization of the newts Triturus cristatus and T. marmoratus. To do so, we compared size, age, and growth-related parameters in F₁ hybrids and both parental species in Mayenne, France. We found significant differences in snout-vent length (SVL), body mass and average lifespan between the parental species and hybrids, increasing from T. cristatus – T. marmoratus – hybrids. The relation between age and SVL fitted von Bertalanffy's growth model and showed that SVL_max was significantly larger in hybrids than in the parental species, while the growth coefficient was lower in hybrids and T. marmoratus than in T. cristatus. Triturus cristatus appears to be the better competitor, since it attains sexual maturity faster and thus achieves more annual breeding opportunities. At the evolutionary level, the observed heterosis appears not to have further consequences as the hybrids are largely infertile. Our results support the hypothesis raised for the genus Triturus, that infertile hybrids allocate resources to growth.     

Characterisation of a short,highly repeated and centromerically localised DNA sequence in crested and marbled newts of the genus Triturus

A 32–33 bp highly repeated DNA sequence, TkS1, has been isolated from genomic DNA of the newt Triturus karelini digested with the restriction endonucleases HaeIII or AluI. TkS1 is known to be localised in the centromeric heterochromatin of all the chromosomes in T. karelini and the related species T. cristatus TkS1 has been shown to be present in varying amounts in the genomic DNA of a range of species of Triturus, including representatives of the two main subgenera Triturus and Palaeotriton. A programme of sequencing of monomers, dimers and trimers of TkS1 was carried out in order to determine the level of conservation of the sequence within and between species of Triturus. Altogether 204 monomer (32/33 bp) clones were made of TkS1 from three individuals of T. karelini, and one individual each of T. cristatus, T. carnifex, T. dobrogicus and T. marmoratus, all members of the subgenus Triturus and the cristatus species group. A number of dimer (64 bp) and trimer (96 bp) clones were also made from DNA of a single specimen of T. karelini digested with HaeIII or AluI. Three distinct types of TkS1 were identified in all species examined, except for T. marmoratus where only two of the types were found. The types were distinguished on the basis of certain recurring divergent patterns in monomers sequenced from T. karelini. Type 1 is mainly characterised by the presence of an AluI site at positions 24–27 and type 3 mainly by the presence of an additional base (C) at position 14. Type 2 normally lacks the AluI site and the C at position 14, as well as having a number of other distinguishing features. TkS1 and its three types have remained remarkably constant in sequence since before the divergence of T. marmoratus from other species in the cristatus species group, about 10 million years ago. Examination of all 204 monomer clones and comparison with consensus sequences for the three types shows less than 5% divergence at any one position in the sequence. There is good evidence from examination of dimer and trimer clones of TkS1 that the different types are intermingled with each other, and all three types are likely to be present on all chromosomes. Dimeric (64 bp) TkS1 clones constructed from AluI fragments of T. karelini DNA show evidence of a trimeric (96 bp) supertype with the pattern type 1-type 3-type 1 that is much more common than would be expected on a random basis. The properties of TkS1 are discussed in relation to models for the involvement of highly repeated satellite DNA in chromosomal growth and evolution.W. Hennig     

Ecogeographical variation of body size in the newt Triturus carnifex: comparing the hypotheses using an information‐theoretic approach

Aim Ecogeographical variation of body size in vertebrates (e.g. Bergmann's rule) has long been recognized. However, the patterns and causes of intra‐specific ecogeographical variation of body size in ectotherms, and in amphibians in particular, are strongly debated. We identified the relationship between bioclimatic variables and body size predicted a priori by alternative hypotheses (heat balance, endurance, seasonality, starvation resistance, water availability, primary productivity, parental investment) proposed to explain ecogeographical patterns of body size in ectotherms, and we evaluated the relative support of these hypotheses in explaining variation in body size of the Italian crested newt, Triturus carnifex. Location Twenty‐three populations covering the whole range of T. carnifex (Austria, Croatia, Italy and Slovenia). Methods We obtained data on body size (snout–vent length, SVL) of 2639 adult newts from direct measurements and the literature; we obtained high‐resolution environmental data for the sampled localities. We used an information‐theoretic approach to evaluate the support of the data for the different hypotheses. We also integrated information on population genetics in our models. Results We observed strong geographical variation of body size. The best Akaike information criterion (AIC) models showed that populations with larger body size are associated with cold climates and secondarily with high primary productivity. Furthermore, sexual dimorphism increases in cold climates, as the increase in body size was stronger for females. When taking into account population genetics, we did not find support for relationships with the other variables. Main conclusion Our results are consistent with three hypotheses proposed to explain ecogeographical variation in amphibians: heat balance, increased parental investment of females and productivity. Information theory provides the framework for comparing hypotheses rather than looking for patterns. We suggest that evaluating the support for mechanisms can provide better insights than simply assessing whether ecogeographical variation is in agreement with some ‘rule’.     

Improving representativeness of genebank collections through species distribution models, gap analysis and ecogeographical maps

An efficient germplasm collecting method was evaluated using six Lupinus species and the Spanish Lupinus collection as a study case. This method includes the application of geographic information systems, ecogeographical land characterization maps, species distribution models and gap analysis to identify prioritized collecting sites. To evaluate the efficiency of this collecting method, field collecting expeditions were carried out focusing on prioritized sites and the results of these collections were analyzed. Prioritized sites were identified using spatial and ecogeographical gaps, and potential species richness maps. The spatial gaps corresponded to populations non-included in the collection but recorded by other information sources while ecogeographical gaps corresponded to spatial gaps that were located in ecogeographical categories (obtained from the ecogeographical map) that were scarcely represented in the collection. A potential Lupinus species richness map was obtained by adding the information of single maps of Lupinus species distribution models. Subsequently, prioritized sites were obtained in ecogeographical gaps with high potential species richness values. Collecting expeditions were made in Spain in 2006, 2007 and 2008. Results showed that using the efficient germplasm collecting methodology was highly positive not only from a quantitative viewpoint (between 7.8 and 11% increase) but also in qualitative terms, focusing collection efforts in ecogeographical categories with low or null representation in the Spanish Lupinus collection (41% of the new accessions). Phenotypic differences related to adaptation to environment were observed in the field between the populations that grow in low or null represented categories and those that grow in highly represented categories.     

Some features of the behaviour of the marbled newt Triturus marmoratus (Latreille, 1800) in captivity

The behaviour in groups of both sexes of marbled newts (Triturus marmoratus) from Sintra, Portugal, kept for a few weeks in a large aquarium, is described. Behaviour patterns described for this species and for other species in this genus were recognised. The following behaviours are described for the first time: Walk-over, Land-over, Stand-on, Show, Erect-body, Yawn, Circle, and Dig. Females seem to have a more active role in courtship than is usually considered. The behaviour frequencies for each sex were analysed using simulation statistics. Some elements usually classified as male sexual behaviours occur in both sexes and are likely to play other social roles. Received in revised form: 31 May 2001 Electronic Publication     

High aquatic niche overlap in the newts Triturus cristatus and T. marmoratus (Amphibia,Urodela)

We studied spatial niche metrics of large-bodied newts (Triturus cristatus and T. marmoratus) in three breeding ponds in western France. Adults and larvae were sampled with underwater funnel traps. Larvae were identified to the species with diagnostic microsatellite DNA markers. The distribution of adult T. cristatus and T. marmoratus across pond regions differed in one out of six cases, no differences were observed between larvae (two ponds studied). Niche overlap and niche breadth indices across resource states defined as pond regions or individual traps were high (Schoener's C: pond regions 0.60–0.98, traps 0.35–0.71; Levins' B: pond regions 0.71–0.98, traps 0.35–0.76). Adults of large-bodied newts significantly differed in resource use from small-bodied newts (T. helveticus). The results are discussed in view of the occurrence of interspecific breeding attempts, and the unpredictable ecological characteristics of newt breeding ponds.     

Taxonomy, biogeography and evolution of Euproctus (Amphibia: Salamandridae), with the resurrection of the genus Calotriton and the description of a new endemic species from the Iberian Peninsula

A phylogeny of West Palearctic Salamandridae based on 1208 bp of mtDNA sequences (300 bp of cytb, 346 bp of 12S rRNA and 562 bp of 16S rRNA) indicates the European brook newts (Euproctus) are polyphyletic. To reflect revised relationships, the Tyrrhenian species (E. montanus ( Savi, 1838 ) and E. platycephalus ( Gravenhorst, 1829 )) are retained in Euproctus Genè, 1839 , while the genus Calotriton Gray, 1858 is resurrected to include the Pyrenean brook newt ( Calotriton asper ( Dugès, 1852 ) comb. nov. ) and a new species from the massif of El Montseny, Catalonia, Spain, described herein as Calotriton arnoldi sp. nov. , which is both morphologically and genetically distinct. Although according to the principle of priority Megapterna Savi, 1838 should take precedence over Euproctus Genè, 1839 , for the sake of nomenclatural stability and in line with Art. 23.9.1 of the International Code of Zoological Nomenclature, Megapterna is considered a nomen oblitum and Euproctus a nomen protectum. The polyphyly of Euproctus (s.l.) contradicts previous, well‐accepted, biogeographical hypotheses and represents a clear case of convergence, involving several morphological traits and a unique reproductive behaviour that is advantageous in stream situations. Molecular dating suggests the Western brook newt lineage (C. asper+C. arnoldi) originated towards the end of the Miocene (8.3 ± 0.11 Mya) and is part of a well‐supported monophyletic assemblage, which also includes Neurergus kaiseri ( Schmidt, 1952 ) and a clade formed by Triturus karelinii ( Strauch, 1870 ), T. carnifex ( Laurenti, 1768 ), T. pygmaeus ( Wolterstoff, 1905 ) and T. marmoratus ( Latreille, 1800 ). Speciation separating E. montanus and E. platycephalus might have coincided with the onset of the Messinian salinity crisis. © 2005 The Linnean Society of London, Zoological Journal of the Linnean Society, 2005, 145 , 555–582.     

Chromosome 1 in crested and marbled newts (Triturus)

The heteromorphic chromosomes 1 of Triturus cristatus carnifex and T. marmoratus were studied in mitotic metaphase after staining with the Giemsa C-banding technique and with the fluorochromes, DAPI (AT-specific) and mithramycin (GC-specific). They were also examined in the lampbrush form under phase-contrast before fixation and after fixation and staining with Giemsa. Chromosomes 1 of T.c. carnifex are asynaptic and achiasmatic throughout most of their long arms. They are also heteromorphic in most of their long arms for the patterns of Giemsa and fluorochrome staining and the distribution of distinctive lampbrush loops. The heteromorphic regions correspond to the regions that are asynaptic and achiasmatic. They stain more strongly with mithramycin and more weakly with DAPI than the remainder of the chromosomes, signifying that their DNA is relatively rich in GC. The patterns of staining with Giemsa and fluorochromes and the distributions of distinctive lateral loops vary from one animal to another in the same species and even in the same population. The asynaptic and achiasmatic regions of chromosomes 1 in T. marmoratus extend throughout the whole of the long arms and well beyond the heterochromatic region. Chiasmata form only in the short arm and occasionally in the short euchromatic segment at the tip of the long arms. The staining patterns of chromosomes 1 in T. marmoratus differ from those in T.c. carnifex although, like carnifex, their DNA is relatively GC-rich. The chromosomes 1 of T. marmoratus are more submetacentric than those of T.c. carnifex. In T. marmoratus chromosome 1B is about 12% shorter than 1A. There is a short paracentric inversion heterozygosity in the long arm of chromosome 1B in T. marmoratus which probably accounts for the lack of chiasmata in the euchromatin that separates the centromere from the start of the heterochromatin. In both carnifex and marmoratus, embryos that are homomorphic for chromosome 1 arrest and die at the late tailbud stage of development. The same applies to F₁ hybrid embryos T.c. carnifex x T. marmoratus, and this has permitted identification of chromosomes 1A and 1B in both species. There is no correspondence between patterns of Giemsa or fluorochrome staining of the heteromorphic regions of chromosome 1 and any feature of the lampbrush chromosomes. However, the short euchromatic ends of the long arms of chromosomes 1 in both species are distinguished in the lampbrush form by a series of uniformly small loops of fine texture associated with very small chromomeres. The Giemsa C-staining patterns of both chromosomes 1A and 1B are different in each of the four subspecies of T. cristatus. T.c. karelinii stands out by having unusually large masses of Giemsa C-staining centromeric heterochromatin on all but 1 of its 12 chromosomes. A scheme is proposed for the evolution of chromosome 1 in T. cristatus and T. marmoratus, based on all available cytological and molecular data.     

RESTRICTED GENE FLOW IN A MOVING HYBRID ZONE OF THE NEWTS TRITURUS CRISTATUS AND T. MARMORATUS IN WESTERN FRANCE

Two hybridizing species of newts, Triturus cristatus and T. marmoratus, with overlapping distributions show a parapatric distribution when surveyed in detail. The factors that govern the distribution of cristatus vs. marmoratus in the département (province) of Mayenne in western France are identified as forestation and relief. The parapatric hybrid zone running through Mayenne is narrow but widens to approximately 20 km in an area with mixed habitat. In this area most breeding sites are shared and F₁ hybrids form about 4% of the total population. Analysis of survey data collected about 30 years previously also shows an essentially parapatric distribution. Comparison of past and present distribution maps reveals that cristatus has superseded marmoratus over large areas in the south of Mayenne. An area where marmoratus replaced cristatus also exists, but it is more limited in size. Gene flow between cristatus and marmoratus is analyzed using 10 diagnostic genetic markers [9 protein loci and mitochondrial (mt) DNA]. In syntopic populations nuclear gene flow is bidirectional with a mean frequency of introgressed alleles (f) of 0.3%. In allotopic populations of cristatus and marmoratus gene flow is present in areas of species replacement (f = 0.3%), while gene flow appears to be absent in those areas that have been continuously occupied by a single species. At the biogeographic level, the presence or absence of introgression is paralleled by the persistence or absence, respectively, of pockets of cristatus–marmoratus syntopy. All F₁ hybrids possess the cristatus type mtDNA. This may be due to asymmetric interspecific mate choice and would explain the observed absence of introgression of the maternally inherited mtDNA genome in areas where cristatus replaced marmoratus. The cristatus–marmoratus hybrid zone bears characteristics of both the clinal (parapatric) hybrid zone model and the mosaic hybrid zone model. Such a mixed model—for which we propose the term “reticulate hybrid zone”—can be appreciated only if studied over a two-dimensional geographic area and also through time.     

Phylogeographic analysis reveals northerly refugia for the riverine amphibian Triturus dobrogicus (Caudata: Salamandridae)

We investigated the recent evolutionary history of the Danube crested newt, Triturus dobrogicus through reconstructions of: (1) the number and position of refugia at the last glacial maximum, (2) the role of major central European rivers in pattern of post‐glacial dispersal, and (3) the present‐day distribution pattern. We analysed sequences of mitochondrial DNA (ND2, 1065 bp) and six microsatellite loci in 363 T. dobrogicus individuals from 58 populations covering the range of the species. Our analyses suggested that T. dobrogicus survived the last glacial maximum in two separate refugia positioned in northwestern Pannonia and in Southern Pannonia from where its range expanded along the Danube and Tisza Rivers. Our findings also confirmed that rivers played an important role in shaping the evolutionary history of amphibian species in Central Europe. We compared the T. dobrogicus range with another lowland amphibian, the fire‐bellied toad Bombina bombina, using species distribution modelling. In line with these models, the isolated mountains inside Pannonia are occupied not by B. bombina, but by its mountain‐dwelling sister‐species B. variegata. However, in contrast to the model, crested newts (including T. dobrogicus) are absent from these mountains. We attribute this biogeographical discrepancy to the positioning of the species’ refugia at the last glacial maximum.     

Integrating detection probabilities in species distribution models of amphibians breeding in Mediterranean temporary ponds

Aim The imperfect detection of species may lead to erroneous conclusions about species–environment relationships. Accuracy in species detection usually requires temporal replication at sampling sites, a time‐consuming and costly monitoring scheme. Here, we applied a lower‐cost alternative based on a double‐sampling approach to incorporate the reliability of species detection into regression‐based species distribution modelling. Location Doñana National Park (south‐western Spain). Methods Using species‐specific monthly detection probabilities, we estimated the detection reliability as the probability of having detected the species given the species‐specific survey time. Such reliability estimates were used to account explicitly for data uncertainty by weighting each absence. We illustrated how this novel framework can be used to evaluate four competing hypotheses as to what constitutes primary environmental control of amphibian distribution: breeding habitat, aestivating habitat, spatial distribution of surrounding habitats and/or major ecosystems zonation. The study was conducted on six pond‐breeding amphibian species during a 4‐year period. Results Non‐detections should not be considered equivalent to real absences, as their reliability varied considerably. The occurrence of Hyla meridionalis and Triturus pygmaeus was related to a particular major ecosystem of the study area, where suitable habitat for these species seemed to be widely available. Characteristics of the breeding habitat (area and hydroperiod) were of high importance for the occurrence of Pelobates cultripes and Pleurodeles waltl. Terrestrial characteristics were the most important predictors of the occurrence of Discoglossus galganoi and Lissotriton boscai, along with spatial distribution of breeding habitats for the last species. Main conclusions We did not find a single best supported hypothesis valid for all species, which stresses the importance of multiscale and multifactor approaches. More importantly, this study shows that estimating the reliability of non‐detection records, an exercise that had been previously seen as a naïve goal in species distribution modelling, is feasible and could be promoted in future studies, at least in comparable systems.     

Distribution and population genetic structure of the Mediterranean pine shoot beetle Tomicus destruens in the Iberian Peninsula and Southern France

1 The Mediterranean pine shoot beetle Tomicus destruens has long been indistinguishable from its congeneric Tomicus piniperda. Both species attack pines, and can be found in sympatry. The geographical distribution of T. destruens is still unclear in most of the Mediterranean Basin. 2 We aimed to describe the geographical distribution and zones of sympatry of both species in the Iberian Peninsula and France, and to study the molecular phylogeographical pattern of T. destruens. 3 Tomicus spp. adults were sampled in Portugal, Spain and France, and a portion of the mitochondrial genes COI and COII was sequenced for 84 individuals. Sequences were aligned to a data set previously obtained from French localities. 4 Tomicus destruens was found in all populations, except for one locality in Portugal and in the Landes (France). It was in sympatry with T. piniperda in two locations on Pinus pinaster and one location on Pinus radiata. 5 Within‐population genetic diversity was high, but we found a significant pattern of spatial distribution of genetic variation, as well as a significant effect of the host tree. 6 The data suggest the existence of two glacial refugia, from which T. destruens recolonized its current range. One refugium was located in Portugal where the beetle probably evolved on P. pinaster. The corresponding haplotypes show a West–East frequency gradient. The other refugium was probably in the eastern range, where the beetles evolved on Pinus halepensis and P. pinea. The corresponding haplotypes show an East–West frequency gradient.     

The crested newt Triturus cristatus recolonized temperate Eurasia from an extra‐Mediterranean glacial refugium

We assess the role of the Carpathians as an extra‐Mediterranean glacial refugium for the crested newt Triturus cristatus. We combine a multilocus phylogeography (one mitochondrial protein‐coding gene, three nuclear introns, and one major histocompatibility complex gene) with species distribution modelling (projected on current and Last Glacial Maximum climate layers). All genetic markers consistently show extensive genetic variation within and genetic depletion outside the Carpathians. The species distribution model suggests that most of the current range was unsuitable at the Last Glacial Maximum, but a small suitable area remained in the Carpathians. Triturus cristatus dramatically expanded its postglacial range, colonizing much of temperate Eurasia from a glacial refugium in the Carpathians. Within the Carpathians, T. cristatus persisted in multiple geographically discrete regions, providing further support for a Carpathian ‘refugia within refugia’ scenario. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 574–587.     

Banding patterns on lampbrush chromosomes of Triturus marmoratus (Amphibia Urodela) by the Giemsa stain

Lampbrush chromosome preparations from the newt species Triturus marmoratus have been submitted to a banding procedure by using a Giemsa stain technique (C-banding) as well as variants of the method. Centromeres, most of telomeres, the nucleolus organizing region and some segments along the chromosome axes appear to be differently stained. The centromere positions have been indicated on the maps of the lampbrush complement of the species. The possible relationships between banding and chromosome structure and organization are briefly discussed.     

Descriptive biogeography of Tomicus (Coleoptera: Scolytidae) species in Spain

Aim, location Tomicus (Coleoptera, Scolytidae) species are some of the principal pests of Eurasian forest and are represented by three coexisting species in Spain, Tomicus piniperda (Linnaeus, 1758), Tomicus destruens (Wollaston, 1865) and Tomicus minor (Harting, 1834). The distribution of two taxa are unknown as they have until recently been considered separate species. Therefore, we model the potential distribution centres and establish the potential distribution limits of Tomicus species in Iberia. We also assess the effectiveness of different models by comparing predicted results with observed data. These results will have application in forest pest management. Methods Molecular and morphological techniques were used to identify species from 254 specimens of 81 plots. For each plot, a Geographical Information System was used to extract a set of 14 environmental (one topographic, six climatic) and biotic variables (seven host tree distributions). General Additive Models and Ecological Niche Factor Analysis models are applied for modelling and predicting the potential distribution of the three especies of Tomicus. Results The results of both modelling methodologies are in agreement. Tomicus destruens is the predominant species in Spain, living in low and hot areas. Tomicus piniperda occurs in lower frequency and prefers wet and cold areas of north‐central Spain. We detected sympatric populations of T. destruens and T. piniperda in Northern coast of Spain, infesting mainly P. pinaster. Tomicus minor is the rarest species, and it occupies a fragmented distribution located in high and wet areas. The remarkable biotic variable is the distribution of P. sylvestris, incorporated into the models of T. destruens and T. piniperda. Main conclusions These results indicate that in wet areas of north‐central Spain where T. piniperda occurs (and possibly the high altitudes of the southern mountains), T. destruens has a climatic distribution limit. In the northern border of this area, both species overlap their distributions and some co‐occurrences were detected. Tomicus minor potentially occurs in high and wet fragmented areas.