We used double digestion restriction site associated DNA (ddRAD) sequencing to discover SNPs in samples across a transect including a hybrid zone between Podarcis carbonelli and Podarcis carbonelli. We used P. bocagei and P. carbonelli samples from the locations at the extremes of the transect as references. We obtained a SNP dataset including all SNPs after removing loci with depth coverage <8, missing data >20%, removing loci containing more than five SNPs, and with more than 70% heterozygosity (complete dataset; 6905 SNPs, 329 individuals). Additionally, we obtained from the complete dataset two other datasets, prior to apply a missing data filter. One dataset contained loci with allele frequencies higher than 0.8 in the reference population containing only parental individuals of one species and lower than 0.2 in the reference population of the other species ("80/20" dataset; 2300 SNPs, 329 individuals); the other dataset comprised diagnostic SNPs between reference populations (diagnostic dataset; 1241 SNPs, 236 individuals) but excluding private alleles from references, i.e. excluding alleles that are not present in the populations of contact. Individuals with missing data >35% were removed from all datasets (the number of individuals reported for each dataset is after applying this filter, but note that the 80/20 and the diagnostic datasets were obtained before applying this filter to the complete dataset). Across datasets, average depth of coverage by individuals was 28 (median = 26.8, min = 12.5, max = 85.8) and by loci was 29 (median = 28.8; min = 15.6; max = 48.6). The analysis of replicate samples (four samples were replicated twice, i.e. were amplified and sequenced in independent libraries and SNP calling was performed independently) showed high levels (99.87%) of multilocus genotype replicability.

AbstractDifferences in seasonal migratory behaviours are thought to be an important component of reproductive isolation in many organisms. Stable isotopes have been used with success in estimating the location and qualities of disjunct breeding and wintering areas. However, few studies have used isotopic data to estimate the movements of hybrid offspring in species that form hybrid zones. Here, we use stable hydrogen to estimate the wintering locations and migratory patterns of two common and widespread migratory birds, Audubon's (Setophaga auduboni) and myrtle (S. coronata) warblers, as well as their hybrids. These two species form a narrow hybrid zone with extensive interbreeding in the Rocky Mountains of British Columbia and Alberta, Canada, which has been studied for over four decades. Isotopes in feathers grown on the wintering grounds or early on migration reveal three important patterns: (1) Audubon's and myrtle warblers from allopatric breeding populations winter in isotopically different environments, consistent with band recovery data and suggesting that there is a narrow migratory transition between the two species, (2) most hybrids appear to overwinter in the south-eastern USA, similar to where myrtle warblers are known to winter, and (3) some hybrid individuals, particularly those along the western edge of the hybrid zone, show Audubon's-like isotopic patterns. These data suggest there is a migratory divide between these two species, but that it is not directly coincident with the centre of the hybrid zone in the breeding range. We interpret these findings and discuss them within the context of previous research on hybrid zones, speciation and migratory divides.

We used double digestion restriction site associated DNA (ddRAD) sequencing to discover SNPs in samples from four contact zones between Podarcis carbonelli and four other Podarcis species. We obtained a panel of SNPs for each for each contact zone and reference populations and a dataset of diagnostic SNPs between reference populations for each contact zone but excluding private alleles from references, i.e. excluding alleles that are not present in the populations of contact. The final datasets (complete and diagnostic) were obtained after removing loci with depth coverage <8, missing data >20% and removing individuals with more than 35% of missing data. Across complete and diagnostic datasets, mean coverage by individuals ranged from 28 to 47 and by loci from 28 to 44. The analysis of replicate samples (about 6% of samples were replicated, i.e. were amplified and sequenced in independent libraries and SNP calling was performed independently) showed high levels (>99%) of multilocus genotype replicability.

Parasites can play a role in speciation, by exerting different selection pressures on different host lineages, leading to reproductive barriers in regions of possible interbreeding. Hybrid zones therefore offer an ideal system to study the effect of parasites on speciation. Here we study a hybrid zone in the foothills of the Rocky Mountains where two yellow-rumped warbler subspecies, Setophaga coronata coronata and S. c. auduboni, interbreed. There is partial reproductive isolation between them, but no evidence of strong assortative mating within the hybrid zone, suggesting the existence of a postzygotic selection against hybrids. Here, we test whether haemosporidian parasites might play a role in selecting against hybrids between S. c. coronata and S. c. auduboni. We screened birds from 5 transects across the hybrid zone for three phylogenetic groupings of avian haemosporidians Plasmodium, Haemoproteus and Leucocytozoon parasites and quantified intensity of infection. Contrary to our prediction, hybrids did not have higher haemosporidian parasite prevalence. Variation in Haemoproteus prevalence was best explained by an interaction between a birds’ hybrid index and elevation, while the probability of infection with Leucocytozoon parasites was only influenced by elevation. We also found no significant difference in the diversity of haemosporidian lineages between the warbler subspecies and their hybrids. Finally, intensity of infection by Haemoproteus increased significantly with elevation, but was not significantly linked to birds’ hybrid index. In conclusion, our data suggest that haemosporidian parasites do not seem to play a major role in selecting against hybrids in this system.

Restriction site-associated DNA sequencing (RAD-seq) provides high-resolution population genomic data at low cost, and has become an important component in ecological and evolutionary studies. As with all high-throughput technologies, analytic strategies require critical validation to ensure accurate and unbiased interpretation. To test for the impact of bioinformatic data processing on downstream population genetic inferences, we analysed mammalian RAD-seq data (>100 individuals) with 312 combinations of methodology (de novo vs. mapping to references of increasing divergence) and filtering criteria (missing data, HWE, FIS, coverage, mapping, genotype quality). In an effort to identify commonalities and biases in all pipelines, we computed summary statistics (nr. loci, nr. SNP, π, Hetobs, FIS, FST, Ne, m) and compared the results to independent null expectations (isolation-by-distance correlation, expected transition-to-transversion ratio Ts/Tv, Mendelian mismatch rates of known parent-offspring trios). We observed large differences between reference-based and de novo approaches, the former generally calling more SNPs and reducing FIS and Ts/Tv. Data completion levels showed little impact on most summary statistics, and FST estimates were robust across all pipelines. The site-frequency spectrum (SFS) was highly sensitive to the chosen approach as reflected in large variance of parameter estimates across demographic scenarios (single-population bottlenecks and isolation-with-migration model). Null-expectations were best met by reference-based approaches, though contingent on the specific criteria. We recommend RAD-seq studies employ reference-based approaches to a closely related genome, and due to the high stochasticity associated with the pipeline advocate the use of multiple pipelines to ensure robust population genetic and demographic inferences.

Patterns of sex-chromosome differentiation and gonadal development have been shown to vary among populations of Rana temporaria along a latitudinal transect in Sweden. Frogs from the northern-boreal population of Ammarnäs displayed well-differentiated X and Y haplotypes, early gonadal differentiation, and a perfect match between phenotypic and genotypic sex. In contrast, no differentiated Y haplotypes could be detected in the southern population of Tvedöra, where juveniles furthermore showed delayed gonadal differentiation. Here, we show that Dmrt1, a gene that plays a key role in sex determination and sexual development across all metazoans, displays significant sex differentiation in Tvedöra, with a Y-specific haplotype distinct from Ammarnäs. The differential segment is not only much shorter in Tvedöra than in Ammarnäs, it is also less differentiated and associates with both delayed gonadal differentiation and imperfect match between phenotypic and genotypic sex. Whereas Tvedöra juveniles with a local Y haplotype tend to ultimately develop as males, those without it may nevertheless become functional XX males, but with strongly female-biased progeny. Our findings suggest that the variance in patterns of sex determination documented in common frogs might result from a genetic polymorphism within a small genomic region that contains Dmrt1. They also substantiate the view that recurrent convergences of sex determination toward a limited set of chromosome pairs may result from the co-option of small genomic regions that harbor key genes from the sex-determination pathway.

We demonstrate a genotyping-by-sequencing approach to identify homomorphic sex chromosomes and their homolog in a distantly related reference genome, based on noninvasive sampling of wild-caught individuals, in the moor frog Rana arvalis. Double-digest RADseq libraries were generated using buccal swabs from 30 males and 21 females from the same population. Search for sex-limited markers from the unfiltered data set (411 446 RAD tags) was more successful than searches from a filtered data set (33 073 RAD tags) for markers showing sex differences in heterozygosity or in allele frequencies. Altogether, we obtained 292 putatively sex-linked RAD loci, 98% of which point to male heterogamety. We could map 15 of them to the Xenopus tropicalis genome, all but one on chromosome pair 1, which seems regularly co-opted for sex determination among amphibians. The most efficient mapping strategy was a three-step hierarchical approach, where R. arvalis reads were first mapped to a low-coverage genome of Rana temporaria (17 My divergence), then the R. temporaria scaffolds to the Nanorana parkeri genome (90 My divergence), and finally the N. parkeri scaffolds to the X. tropicalis genome (210 My). We validated our conclusions with PCR primers amplifying part of Dmrt1, a candidate sex determination gene mapping to chromosome 1: a sex-diagnostic allele was present in all 30 males but in none of the 21 females. Our approach is likely to be productive in many situations where biological samples and/or genomic resources are limited.

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

Background: Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications. Findings: We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement. Conclusions: The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

Identifying homology between sex chromosomes of different species is essential to understanding the evolution of sex determination. Here, we show that the identity of a homomorphic sex chromosome pair can be established using a linkage map, without information on offspring sex. By comparing sex-specific maps of the European tree frog Hyla arborea, we find that the sex chromosome (linkage group 1) shows a threefold difference in marker number between the male and female maps. In contrast, the number of markers on each autosome is similar between the two maps. We also find strongly conserved synteny between H. arborea and Xenopus tropicalis across 200 million years of evolution, suggesting that the rate of chromosomal rearrangement in anurans is low. Finally, we show that recombination in males is greatly reduced at the centers of large chromosomes, consistent with previous cytogenetic findings. Our research shows the importance of high-density linkage maps for studies of recombination, chromosomal rearrangement and the genetic architecture of ecologically or economically important traits.