Limited evidence for genetic differentiation or adaptation in two amphibian species across replicated rural-urban gradients.

Urbanization leads to complex environmental changes and poses multiple challenges to organisms. Amphibians are highly susceptible to the effects of urbanization, with land use conversion, habitat destruction, and degradation ranked as the most significant threats. Consequently, amphibians are declining in urban areas, in both population numbers and abundance, however, the effect of urbanization on population genetic parameters remains unclear. Here, we studied the genomic response to urbanization in two widespread European species, the common toad Bufo bufo (26 localities, 480 individuals), and the smooth newt Lissotriton vulgaris (30 localities, 516 individuals) in three geographic regions: southern and northern Poland and southern Norway. We assessed genome-wide SNP variation using RADseq (ca. 42 and 552 thousand SNPs in toads and newts, respectively) and adaptively relevant major histocompatibility complex (MHC) class I and II genes. The results linked most of the genetic differentiation in both marker types to regional (latitudinal) effects, which also correspond to historical biogeography. Further, we did not find any association between genetic differentiation and level of urbanization at local scales for either species. However, urban smooth newts, but not toads, have lower levels of within-population genome-wide diversity, suggesting higher susceptibility to the negative effects of urbanization. A decreasing level of genetic diversity linked to increasing urbanization was also found for MHC II in smooth newts, while the relationship between MHC class I diversity and urbanization differed between geographic regions. We did not find any effects of urbanization on MHC diversity in the toad populations. Although two genetic environment association analyses of genome-wide data, LFMM and BayPass, revealed numerous (219 in B. bufo and 7040 in L. vulgaris) SNPs statistically associated with urbanization, we found a marked lack of repeatability between geographic regions, suggesting a complex and multifaceted response to natural selection elicited by life in the city.

The genomic response to urbanization in the damselfly Ischnura elegans.

The complex and rapid environmental changes brought about by urbanization pose significant challenges to organisms. The multifaceted effects of urbanization often make it difficult to define and pinpoint the very nature of adaptive urban phenotypes. In such situations, scanning genomes for regions differentiated between urban and non-urban populations may be an attractive approach. Here, we investigated the genomic signatures of adaptation to urbanization in the damselfly Ischnura elegans sampled from 31 rural and urban localities in three geographic regions: southern and northern Poland, and southern Sweden. Genome-wide variation was assessed using more than 370,000 single nucleotide polymorphisms (SNPs) genotyped by ddRADseq. Associations between SNPs and the level of urbanization were tested using two genetic environment association methods: Latent Factors Mixed Models and BayPass. While we found numerous candidate SNPs and a highly significant overlap between candidates identified by the two methods within the geographic regions, there was a distinctive lack of repeatability between the geographic regions both at the level of individual SNPs and of genomic regions. However, we found "synapse organization" at the top of the functional categories enriched among the genes located in the proximity of the candidate urbanization SNPs. Interestingly, the overall significance of "synapse organization" was built up by the accretion of different genes associated with candidate SNPs in different geographic regions. This finding is consistent with the highly polygenic nature of adaptation, where the response may be achieved through a subtle adjustment of allele frequencies in different genes that contribute to adaptive phenotypes. Taken together, our results point to a polygenic adaptive response in the nervous system, specifically implicating genes involved in synapse organization, which mirrors the findings from several genomic and behavioral studies of adaptation to urbanization in other taxa.

Divergence history of the Carpathian and smooth newts modelled in space and time.

Information about demographic history is essential for the understanding of the processes of divergence and speciation. Patterns of genetic variation within and between closely related species provide insights into the history of their interactions. Here, we investigated historical demography and genetic exchange between the Carpathian (Lissotriton montandoni, Lm) and smooth (L. vulgaris, Lv) newts. We combine an extensive geographical sampling and multilocus nuclear sequence data with the approximate Bayesian computation framework to test alternative scenarios of divergence and reconstruct the temporal and spatial pattern of gene flow between species. A model of recent (last glacial period) interspecific gene flow was favoured over alternative models. Thus, despite the relatively old divergence (4-6 mya) and presumably long periods of isolation, the species have retained the ability to exchange genes. Nevertheless, the low migration rates (ca. 10(-6) per gene copy per generation) are consistent with strong reproductive isolation between the species. Models allowing demographic changes were favoured, suggesting that the effective population sizes of both species at least doubled as divergence reaching the current ca. 0.2 million in Lm and 1 million in Lv. We found asymmetry in rates of interspecific gene flow between Lm and one evolutionary lineage of Lv. We suggest that intraspecific polymorphism for hybrid incompatibilities segregating within Lv could explain this pattern and propose further tests to distinguish between alternative explanations. Our study highlights the importance of incorporating intraspecific genetic structure into the models investigating the history of divergence.

Molecular Inversion Probes for targeted resequencing in non-model organisms.

Applications that require resequencing of hundreds or thousands of predefined genomic regions in numerous samples are common in studies of non-model organisms. However few approaches at the scale intermediate between multiplex PCR and sequence capture methods are available. Here we explored the utility of Molecular Inversion Probes (MIPs) for the medium-scale targeted resequencing in a non-model system. Markers targeting 112 bp of exonic sequence were designed from transcriptome of Lissotriton newts. We assessed performance of 248 MIP markers in a sample of 85 individuals. Among the 234 (94.4%) successfully amplified markers 80% had median coverage within one order of magnitude, indicating relatively uniform performance; coverage uniformity across individuals was also high. In the analysis of polymorphism and segregation within family, 77% of 248 tested MIPs were confirmed as single copy Mendelian markers. Genotyping concordance assessed using replicate samples exceeded 99%. MIP markers for targeted resequencing have a number of advantages: high specificity, high multiplexing level, low sample requirement, straightforward laboratory protocol, no need for preparation of genomic libraries and no ascertainment bias. We conclude that MIP markers provide an effective solution for resequencing targets of tens or hundreds of kb in any organism and in a large number of samples.

Plant-herbivorous beetle networks: molecular characterization of trophic ecology within a threatened steppic environment.

DNA barcoding facilitates many evolutionary and ecological studies, including the examination of the dietary diversity of herbivores. In this study, we present a survey of ecological associations between herbivorous beetles and host plants from seriously threatened European steppic grasslands. We determined host plants for the majority (65%) of steppic leaf beetles (55 species) and weevils (59) known from central Europe using two barcodes (trnL and rbcL) and two sequencing strategies (Sanger for mono/oligophagous species and Illumina for polyphagous taxa). To better understand the ecological associations between steppic beetles and their host plants, we tested the hypothesis that leaf beetles and weevils differ in food selection as a result of their phylogenetic relations (within genera and between families) and interactions with host plants. We found 224 links between the beetles and the plants. Beetles belonging to seven genera feed on the same or related plants. Their preferences were probably inherited from common ancestors and/or resulted from the host plant's chemistry. Beetles from four genera feed on different plants, possibly reducing intrageneric competition and possibly due to an adaptation to different plant chemical defences. We found significant correlations between the numbers of leaf beetle and weevil species feeding on particular plants for polyphagous taxa, but not for nonpolyphagous beetles. Finally, we found that the previous identifications of host plants based on direct observations are generally concordant with host plant barcoding from insect gut. Our results expand basic knowledge about the trophic relations of steppic beetles and plants and are immediately useful for conservation purposes.

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure.

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre-eminent system for the study of selective pressures that arise from host-pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population-genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.

Accuracy of allele frequency estimation using pooled RNA-Seq.

For nonmodel organisms, genome-wide information that describes functionally relevant variation may be obtained by RNA-Seq following de novo transcriptome assembly. While sequencing has become relatively inexpensive, the preparation of a large number of sequencing libraries remains prohibitively expensive for population genetic analyses of nonmodel species. Pooling samples may be then an attractive alternative. To test whether pooled RNA-Seq accurately predicts true allele frequencies, we analysed the liver transcriptomes of 10 bank voles. Each sample was sequenced both as an individually barcoded library and as a part of a pool. Equal amounts of total RNA from each vole were pooled prior to mRNA selection and library construction. Reads were mapped onto the de novo assembled reference transcriptome. High-quality genotypes for individual voles, determined for 23,682 SNPs, provided information on 'true' allele frequencies; allele frequencies estimated from the pool were then compared with these values. 'True' frequencies and those estimated from the pool were highly correlated. Mean relative estimation error was 21% and did not depend on expression level. However, we also observed a minor effect of interindividual variation in gene expression and allele-specific gene expression influencing allele frequency estimation accuracy. Moreover, we observed strong negative relationship between minor allele frequency and relative estimation error. Our results indicate that pooled RNA-Seq exhibits accuracy comparable with pooled genome resequencing, but variation in expression level between individuals should be assessed and accounted for. This should help in taking account the difference in accuracy between conservatively expressed transcripts and these which are variable in expression level.

Development, validation and high-throughput analysis of sequence markers in nonmodel species.

DNA sequences derived from multiple regions of the nuclear genome are essential for historical inferences in the fields of phylogeography and phylogenetics. The appropriate markers should be single-copy, variable, easy to amplify from multiple samples and easy to sequence using high-throughput technologies. This may be difficult to achieve for species lacking sequenced genomes and particularly challenging for species possessing large genomes, which consist mostly of repetitive sequences. Here, we present a cost-effective, broadly applicable framework for designing, validating and high-throughput sequencing of multiple markers in nonmodel species without sequenced genomes. We demonstrate its utility in two closely related species of newts, representatives of urodeles, a vertebrate group characterized by large genomes. We show that over 80 markers, c. 600 bp each, developed mainly from 3' untranslated transcript regions (3'UTR) may be effectively multiplexed and sequenced. Data are further processed using standard, freely available bioinformatic tools, producing phase-resolved sequences. The approach does not require barcoded PCR primers, and the cost of library preparation is independent of the number of markers investigated. We hope that this approach will be of broad interest for researchers working at the interface of population genetics and phylogenetics, exploring deep intraspecific genetic structure, species boundaries and phylogeographies of closely related species.

No evidence for nuclear introgression despite complete mtDNA replacement in the Carpathian newt (Lissotriton montandoni).

Patterns of interspecific introgression may vary geographically, and the distribution of introgressed variants can yield insight into the historical dynamics of genetic interactions between hybridizing species. Urodele amphibians, often characterized by limited mobility, deep intraspecific genetic structuring and vulnerability to climatic changes, constitute suitable models for such historical inferences. Here, we combine an extensive survey of the mitochondrial (mtDNA) and nuclear (15 microsatellites) genomes in the Carpathian newt, Lissotriton montandoni (Lm) with species distribution modelling (SDM). Populations of the smooth newt, L. vulgaris (Lv) from the areas surrounding the Lm range were also sampled to test whether gene flow between these hybridizing species extends beyond the area of strict syntopy. The extent of introgression differs dramatically between the mitochondrial genome and the nuclear genome. While multiple, spatially and temporally distinct introgression events from Lv resulted in complete mtDNA replacement in Lm, there was little evidence of recent interspecific nuclear gene flow in the assayed markers. Microsatellite differentiation within Lm defines three units, probably derived from separate glacial refugia, located in the northern, eastern and southern part of the Carpathians. In situ survival and range fragmentation of Lm are supported by SDM, corroborating the role of the Carpathians as a major refugial area. Our results, in combination with previous reports of extensive introgression of the major histocompatibility complex (MHC) genes, emphasize the complexity of historical gene exchange between Lm and Lv.

Methods for MHC genotyping in non-model vertebrates.

Genes of the major histocompatibility complex (MHC) are considered a paradigm of adaptive evolution at the molecular level and as such are frequently investigated by evolutionary biologists and ecologists. Accurate genotyping is essential for understanding of the role that MHC variation plays in natural populations, but may be extremely challenging. Here, I discuss the DNA-based methods currently used for genotyping MHC in non-model vertebrates, as well as techniques likely to find widespread use in the future. I also highlight the aspects of MHC structure that are relevant for genotyping, and detail the challenges posed by the complex genomic organization and high sequence variation of MHC loci. Special emphasis is placed on designing appropriate PCR primers, accounting for artefacts and the problem of genotyping alleles from multiple, co-amplifying loci, a strategy which is frequently necessary due to the structure of the MHC. The suitability of typing techniques is compared in various research situations, strategies for efficient genotyping are discussed and areas of likely progress in future are identified. This review addresses the well established typing methods such as the Single Strand Conformation Polymorphism (SSCP), Denaturing Gradient Gel Electrophoresis (DGGE), Reference Strand Conformational Analysis (RSCA) and cloning of PCR products. In addition, it includes the intriguing possibility of direct amplicon sequencing followed by the computational inference of alleles and also next generation sequencing (NGS) technologies; the latter technique may, in the future, find widespread use in typing complex multilocus MHC systems.

Long-term survival of a urodele amphibian despite depleted major histocompatibility complex variation.

Depletion of polymorphism at major histocompatibility complex (MHC) genes has been hypothesized to limit the ability of populations to respond to emerging pathogens, thus putting their survival at risk. As pathogens contribute substantially to the global amphibian decline, assessing patterns of MHC variation is important in devising conservation strategies. Here, we directly compare levels of MHC class II and neutral variation between multiple populations of the great crested newt (Triturus cristatus) from refugial (REF: Romania) and postglacial expansion (PGE: Germany, Poland and UK) areas. REF populations harboured high levels of adaptive variation (24 expressed alleles), exhibiting clear signatures of historical positive selection, which points to the overall importance of MHC class II variation in this species. On the other hand, PGE populations were extremely depauperate (two alleles) but nevertheless have survived for c. 10,000 years, since the postglacial expansion. Variation in putative MHC class II pseudogenes, microsatellites and allozymes also showed a significant southern richness-northern purity pattern. The populations in the postglacial expansion area thus provide the clearest example to date of the long-term survival of populations in which MHC variation, historically under positive selection, has been depleted.

A combination of techniques proves useful in the development of nuclear markers in the newt genus Triturus.

To increase the number of markers available for study of phylogeny and phylogeography in the newt genus Triturus, we developed and tested 59 primer pairs using three different techniques. Primers were obtained from published sources, by designing exon-primed intron-crossing primers and from randomly cloned anonymous nuclear DNA fragments. Successful polymerase chain reaction products were cloned and sequenced. Five fragments were successfully amplified and sequenced for six species of Triturus: intron 7 of the ß-fibrinogen gene (ßfibint7), third intron of the calreticulin gene (CalintC), the 11th intron of the α-subunit of the platelet derived growth factor receptor (PDGFRα) and two anonymous markers (Cri1 and Cri4). The average percentage species divergence across all the markers is low (c. 3%), compared to what has been found in mitochondrial DNA (25-30%).

Contrasting patterns of variation in MHC loci in the Alpine newt.

Major histocompatibility complex (MHC) genes are essential in pathogen recognition and triggering an adaptive immune response. Although they are the most polymorphic genes in vertebrates, very little information on MHC variation and patterns of evolution are available for amphibians, a group known to be declining rapidly worldwide. As infectious diseases are invoked in the declines, information on MHC variation should contribute to devising appropriate conservation strategies. In this study, we examined MHC variation in 149 Alpine newts (Mesotriton alpestris) from three allopatric population groups in Poland at the northeastern margin of the distribution of this species. The genetic distinctiveness of the population groups has previously been shown by studies of skin graft rejection, allozymes and microsatellites. Two putative expressed MHC II loci with contrasting levels of variation and clear evidence of gene conversion/recombination between them were detected. The Meal-DAB locus is highly polymorphic (37 alleles), and shows evidence of historical positive selection for amino acid replacements and substantial geographical differentiation in allelic richness. On the contrary, the Meal-DBB locus exhibits low polymorphism (three alleles differing by up to two synonymous substitutions) and a uniform distribution of three alleles among geographical regions. The uniform frequencies of the presumptively neutral Meal-DBB alleles may be explained by linkage to Meal-DAB. We found differences in allelic richness in Meal-DAB between regions, consistent with the hypothesis that genetic drift prevails with increasing distance from glacial refugia. Pseudogene loci appear to have evolved neutrally. The level of DAB variation correlated with variation in microsatellite loci, implying that selection and drift interplayed to produce the pattern of MHC variation observed in marginal populations of the Alpine newt.

MHC-DRB3 variation in a free-living population of the European bison, Bison bonasus.

MHC genes play a crucial role in pathogen recognition and are the most polymorphic genes in vertebrates. Loss of variation in these genes in bottlenecked species is thought to put their survival at risk. We examined variation at the MHC II DRB3 locus in the European bison, Bison bonasus, a species that has undergone an extreme bottleneck: the current population originated from only 12 founders. We also tested for the association of DRB3 genes with the incidence of posthitis, a disease affecting the reproductive organs of bulls and posing a new threat to the survival of the species. We found very limited MHC diversity, with only four alleles segregating in a sample of 172 individuals from a free-ranging Bialowieza population. The alleles were highly divergent and revealed the hallmark of positive selection acting on them in the past, that is, a significant excess of nonsynonymous substitutions. This excess was concentrated in putative antigen-binding sites, suggesting that selection was driven by pathogens. However, we did not observe departures from Hardy-Weinberg equilibrium, an indicator of strong ongoing selection. Neither have we found a significant association between DRB3 alleles or genotypes and susceptibility to posthitis. Alleles conferring resistance to males may have been lost during the extreme bottleneck the species had undergone.

Genetic structure in northeastern populations of the Alpine newt (Triturus alpestris): evidence for post-Pleistocene differentiation.

Genetic variation in 13 populations of the Alpine newt, Triturus alpestris, was assessed at the northeastern margin of its range (southern Poland). Variation at six microsatellite loci was scored in 354 newts, and two mitochondrial DNA fragments (c. 2000 bp) were sequenced in a subset of 27 individuals. Significant differences in allele frequencies and the presence of private alleles determined genetic units corresponding to three separate mountain ranges, i.e. the Carpathian, Sudetes and Holy Cross Mountains. F(ST)'s were three times greater in among than in within mountain range pairwise comparisons. An assignment test and pairwise F(ST)'s suggested relatively high levels of gene flow at the local level, although the Sudetes populations revealed some subtle structuring. Genetic variation was lower in the Carpathians and Holy Cross Mountains. The geographic pattern of mitochondrial DNA variation indicated that these newt populations originated from a single glacial refugium/founder population, and that the colonization of southern Poland took place in an easterly direction. The data show that substantial neutral variation and between group divergence has accumulated relatively quickly in these low-vagility organisms. The Alpine newt case exemplifies species history as a factor determining patterns of genetic diversity in marginal populations.

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber).

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.

Phylogeography of two European newt species--discordance between mtDNA and morphology.

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5-1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long-term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.

Mitochondrial phylogeography of the moor frog, Rana arvalis.

The moor frog Rana arvalis is a lowland species with a broad Eurasiatic distribution, from arctic tundra through forest to the steppe zone. Its present-day range suggests that glacial refugia of this frog were located outside southern European peninsulas. We studied the species-wide phylogeographical pattern using sequence variation in a 682 base pairs fragment of mtDNA cytochrome b gene; 223 individuals from 73 localities were analysed. Two main clades, A and B, differing by c. 3.6% sequence divergence were detected. The A clade is further subdivided into two subclades, AI and AII differing by 1.0%. All three lineages are present in the Carpathian Basin (CB), whereas the rest of the species range, including huge expanses of Eurasian lowlands, are inhabited solely by the AI lineage. We infer that AII and B lineages survived several glacial cycles in the CB but did not expand, at least in the present interglacial, to the north. The geographical distribution and genealogical relationships between haplotypes from the AI lineage indicate that this group had two glacial refugia, one located in the eastern part of the CB and the other probably in southern Russia. Populations from both refugia contributed to the colonization of the western part of the range, whereas the eastern part was colonized from the eastern refugium only. The effective population size as evidenced by theta(ML) is an order of magnitude higher in the AI lineage than in the AII and B lineages. Demographic expansion was detected in all three lineages.

Nuclear markers, mitochondrial DNA and male secondary sexual traits variation in a newt hybrid zone (Triturus vulgaris x T. montandoni).

The smooth and the Montandon's newts (Triturus vulgaris and T. montandoni) are genetically similar sister species with highly divergent male secondary sexual traits involved in complex courtship behaviour. Their parapatric ranges overlap at moderate elevations in the Carpathian Mountains where they hybridize readily. Here we present a detailed study of genetic and morphological variation in populations from the area of sympatry. Analysis of variation at seven nuclear markers, mtDNA and male sexual secondary traits was complemented with an ecological survey of breeding sites characteristics. Extensive hybridization was revealed with back-cross individuals similar to either parental species predominating among hybrids. The hybrid zone exhibited a mosaic pattern: the genetic composition of the populations was correlated only weakly with their geographical position. No association with habitat type was found. Departures from Hardy-Weinberg proportions, significant linkage disequilibria and bimodal distribution of genotypes suggest strongly that assortative mating is an important factor shaping the genetic composition of hybrid populations. The pattern of cytonuclear disequilibria did not indicate much asymmetry in interspecific matings. Changes in the frequency of nuclear markers were highly concordant, whereas mtDNA showed much wider bidirectional introgression with 14% excess of T. montandoni haplotype. We argue that the mosaic structure of the newt hybrid zone results mainly from stochastic processes related to extinction and recolonization. Microgeographical differences in mtDNA introgression are explained by historical range shifts. Since morphologically intermediate males were underrepresented when compared to hybrid males identified by genetic markers, sexual selection acting against the morphological intermediates is implied. We discuss the implications of these findings in the context of reinforcement of prezygotic isolation in newts.

Genetic differentiation of the two subspecies of the smooth newt inhabiting Romania, Triturus vulgaris vulgaris and T.v. ampelensis (Urodela, Salamandridae) as revealed by enzyme electrophoresis.

Starch gel electrophoresis was applied to quantify genetic variation and divergence in samples from Romanian populations of the nominal form of the smooth newt Triturus vulgaris and those of the endemic Romanian subspecies T.v. ampelensis, a population from a parapatric area was additionally included. All the samples had similar levels of genetic variation measured by the mean heterozygosity, proportion of polymorphic loci, and mean number of alleles per locus. T.v. ampelensis samples were genetically clearly different from the nominal form samples, the mean genetic distance between the two subspecies was being estimated as DN = 0.114. No fixed differences in allele composition between the two subspecies were found, although some of the alleles were found either exclusively in the nominal form (Aat-1) or in T.v. ampelensis (Mpi a). Other alleles at these loci together with Mdh-1 differed markedly in frequency. The population from the parapatric area was intermediate in allelic composition, but grouped together with the T.v. ampelensis samples in a maximum likelihood tree (99.7% bootstrap support for this grouping). The data indicate that the two subspecies interbreed in a parapatric zone. The molecular clock applied to electrophoretic data indicates that these two forms split during the Pleistocene.