Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

Significance Snake venoms are toxic protein cocktails used for prey capture. To investigate the evolution of these complex biological weapon systems, we sequenced the genome of a venomous snake, the king cobra, and assessed the composition of venom gland expressed genes, small RNAs, and secreted venom proteins. We show that regulatory components of the venom secretory system may have evolved from a pancreatic origin and that venom toxin genes were co-opted by distinct genomic mechanisms. After co-option, toxin genes important for prey capture have massively expanded by gene duplication and evolved under positive selection, resulting in protein neofunctionalization. This diverse and dramatic venom-related genomic response seemingly occurs in response to a coevolutionary arms race between venomous snakes and their prey.

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.

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.

No abstract