Two transects reveal remarkable variation in gene flow on opposite ends of a European toad hybrid zone.

Speciation entails a reduction in gene flow between lineages. The rates at which genomic regions become isolated varies across space and time. Barrier markers are linked to putative genes involved in (processes of) reproductive isolation, and, when observed over two transects, indicate species-wide processes. In contrast, transect-specific putative barrier markers suggest local processes. We studied two widely separated transects along the 900 km hybrid zone between Bufo bufo and B. spinosus, in northern and southern France, for ~1200 RADseq markers. We used genomic and geographic cline analyses to identify barrier markers based on their restricted introgression, and found that some markers are transect-specific, while others are shared between transects. Twenty-six barrier markers were shared across both transects, of which some are clustered in the same chromosomal region, suggesting that their associated genes are involved in reduced gene flow across the entire hybrid zone. Transect-specific barrier markers were twice as numerous in the southern than in the northern transect, suggesting that the overall barrier effect is weaker in northern France. We hypothesize that this is consistent with a longer period of secondary contact in southern France. The smaller number of introgressed genes in the northern transect shows considerably more gene flow towards the southern (B. spinosus) than the northern species (B. bufo). We hypothesize that hybrid zone movement in northern France and hybrid zone stability in southern France explain this pattern. The Bufo hybrid zone provides an excellent opportunity to separate a general barrier effect from localized gene flow-reducing conditions.

Phylogenomics of the adaptive radiation of Triturus newts supports gradual ecological niche expansion towards an incrementally aquatic lifestyle.

Newts of the genus Triturus (marbled and crested newts) exhibit substantial variation in the number of trunk vertebrae (NTV) and a higher NTV corresponds to a longer annual aquatic period. Because the Triturus phylogeny has thwarted resolution to date, the evolutionary history of NTV, annual aquatic period, and their potential coevolution has remained unclear. To resolve the phylogeny of Triturus, we generated a c. 6000 transcriptome-derived marker data set using a custom target enrichment probe set, and conducted phylogenetic analyses using: (1) data concatenation with RAxML, (2) gene-tree summary with ASTRAL, and (3) species-tree estimation with SNAPP. All analyses produce the same, highly supported topology, despite cladogenesis having occurred over a short timeframe, resulting in short internal branch lengths. Our new phylogenetic hypothesis is consistent with the minimal number of inferred changes in NTV count necessary to explain the diversity in NTV observed today. Although a causal relationship between NTV, body form, and aquatic ecology has yet to be experimentally established, our phylogeny indicates that these features have evolved together, and suggest that they may underlie the adaptive radiation that characterizes Triturus.