Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution.

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

The cost of travel: How dispersal ability limits local adaptation in host-parasite interactions.

Classical theory suggests that parasites will exhibit higher fitness in sympatric relative to allopatric host populations (local adaptation). However, evidence for local adaptation in natural host-parasite systems is often equivocal, emphasizing the need for infection experiments conducted over realistic geographic scales and comparisons among species with varied life history traits. Here, we used infection experiments to test how two trematode (flatworm) species (Paralechriorchis syntomentera and Ribeiroia ondatrae) with differing dispersal abilities varied in the strength of local adaptation to their amphibian hosts. Both parasites have complex life cycles involving sequential transmission among aquatic snails, larval amphibians and vertebrate definitive hosts that control dispersal across the landscape. By experimentally pairing 26 host-by-parasite population infection combinations from across the western USA with analyses of host and parasite spatial genetic structure, we found that increasing geographic distance-and corresponding increases in host population genetic distance-reduced infection success for P. syntomentera, which is dispersed by snake definitive hosts. For the avian-dispersed R. ondatrae, in contrast, the geographic distance between the parasite and host populations had no influence on infection success. Differences in local adaptation corresponded to parasite genetic structure; although populations of P. syntomentera exhibited ~10% mtDNA sequence divergence, those of R. ondatrae were nearly identical (<0.5%), even across a 900 km range. Taken together, these results offer empirical evidence that high levels of dispersal can limit opportunities for parasites to adapt to local host populations.

The geographic mosaic in parallel: Matching patterns of newt tetrodotoxin levels and snake resistance in multiple predator-prey pairs.

The Geographic Mosaic Theory of Coevolution predicts that coevolutionary arms races will vary over time and space because of the diverse ecological settings and population histories of interacting species across the landscape. Thus, understanding coevolution may require investigating broad sets of populations sampled across the range of the interaction. In addition, comparing coevolutionary dynamics between similar systems may reveal the importance of specific factors that structure coevolution. Here, we examine geographic patterns of prey traits and predator traits in the relatively unstudied interaction between the Sierra garter snake (Thamnophis couchii) and sympatric prey, the rough-skinned newt (Taricha granulosa), Sierra newt (Ta. sierrae) and California newt (Ta. torosa). This system parallels, in space and phenotypes, a classic example of coevolution between predatory common garter snakes (Th. sirtalis) and their toxic newt prey exhibiting hotspots of newt tetrodotoxin (TTX) levels and matching snake TTX resistance. We quantified prey and predator traits from hundreds of individuals across their distributions, and functional trait matching at sympatric sites. We show strong regional patterns of trait covariation across the shared ranges of Th. couchii and newt prey. Traits differ significantly among localities, with lower newt TTX levels and snake TTX resistance at the northern latitudes, and higher TTX levels and snake resistance at southern latitudes. Newts and snakes in northern populations show the highest degree of functional trait matching despite possessing the least extreme traits. Conversely, newts and snakes in southern populations show the greatest mismatch despite possessing exaggerated traits, with some snakes so resistant to TTX they would be unaffected by any sympatric newt. Nevertheless, individual variation was substantial, and appears to offer the opportunity for continued reciprocal selection in most populations. Overall, the three species of newts appear to be engaged in a TTX-mediated arms race with Th. couchii. These patterns are congruent with those seen between newts and Th. sirtalis, including the same latitudinal gradient in trait covariation, and the potential 'escape' from the arms race by snake predators. Such concordance in broad scale patterns across two distinct systems suggests common phenomena might structure geographic mosaics in similar ways.

Dating and biogeographical patterns in the sea slug genus Acanthodoris Gray, 1850 (Mollusca, Gastropoda, Nudibranchia).

Recent studies investigating vicariance and dispersal have been focused on correlating major geological events with instances of taxonomic expansion by incorporating the fossil record with molecular clock analyses. However, this approach becomes problematic for soft-bodied organisms that are poorly represented in the fossil record. Here, we estimate the phylogenetic relationships of the nudibranch genus Acanthodoris Gray, 1850 using three molecular markers (16S, COI, H3), and then test two alternative geologically calibrated molecular clock scenarios in BEAST and their effect on ancestral area reconstruction (AAR) estimates employed in LAGRANGE. The global temperate distribution of Acanthodoris spans multiple geological barriers, including the Bering Strait (∼5.32 Mya) and the Baja Peninsula (∼5.5 Mya), both of which are used in our dating estimates. The expansion of the Atlantic Ocean (∼95-105 Mya) is also used to calibrate the relationship between A. falklandica Eliot, 1905 and A. planca Fahey and Valdés, 2005, which are distributed in southern Chile and South Africa respectively. Phylogenetic analyses recovered strong biogeographical signal and recovered two major clades representing northern and southern hemispheric distributions of Acanthodoris. When all three geological events are applied to the calibration analyses, the age for Acanthodoris is estimated to be mid-Cretaceous. When the expansion of the Atlantic Ocean is excluded from our analyses, however, Acanthodoris is estimated to be much younger, with a divergence time estimate during the Miocene. Regardless of divergence estimates, our AAR suggests that Acanthodoris may have origins in the Atlantic Ocean with the Atlantic acting as a dispersal point to the northeastern Pacific. These results suggest that Acanthodoris exhibits a rare instance of western trans-arctic expansion. This study also shows that northeast Pacific specimens of A. pilosa should be regarded as A. atrogriseata and that A. serpentinotus should be regarded as a synonym of A. pina.

Family matters: The first molecular phylogeny of the Onchidorididae Gray, 1827 (Mollusca, Gastropoda, Nudibranchia).

Recent investigations into the evolution of the Onchidorididae using morphological based methods have resulted in low support for relationships among genera. This study aims to determine if molecular data corroborates recent morphological interpretations of the evolution of Onchidorididae. Five genetic markers: 16S, 18S, 28S, cytochrome c oxidase 1 (COI) and histone 3 (H3), were sequenced from 32 species comprising Onchidorididae and five other families, three from Phanerobranchia and two from Cryptobranchia. Phylogenies were estimated using maximum likelihood, and Bayesian inference analyses; with both yielding similar topologies. Molecular analyses resulted in high support for the monophyly of the suctorian clade and the placement of the genera within Onchidorididae. However, the Onchidorididae forms a paraphyletic grouping due to the recovery of the Goniodorididae and the Akiodorididae nested within family. In addition, the placement of Corambe as the most derived member of Onchidorididae is contradicted by the present study. Rather it is sister to a large clade that includes Acanthodoris and the species traditionally placed in Onchidoris and Adalaria, now defined as Onchidorididae. We have chosen to maintain Corambidae as a distinct taxon (including Corambe and Loy), sister to Onchidorididae. We also maintain Goniodorididae, Akiodorididae and Calycidoridae (including Calycidoris and Diaphorodoris), which along with the Onchidorididae and Corambidae comprise the suctorian superfamily Onchidoridoidea. Ancestral character reconstruction also suggests that the formation of a gill pocket, a character that currently defines the Cryptobranchia, may have evolved multiple times from an ancestor that lacked the ability to retract its gills into a fully formed gill pocket. The diversity of gill morphology displayed by the Onchidoridoidea will help give new insight into the evolution of this complex character within the Nudibranchia.

Aligning evidence: concerns regarding multiple sequence alignments in estimating the phylogeny of the Nudibranchia suborder Doridina.

Molecular estimates of phylogenetic relationships rely heavily on multiple sequence alignment construction. There has been little consensus, however, on how to properly address issues pertaining to the alignment of variable regions. Here, we construct alignments from four commonly sequenced molecular markers (16S, 18S, 28S and cytochrome c oxidase subunit I) for the Nudibranchia using three different methodologies: (i) strict mathematical algorithm; (ii) exclusion of variable or divergent regions and (iii) manually curated, and examine how different alignment construction methods can affect phylogenetic signal and phylogenetic estimates for the suborder Doridina. Phylogenetic informativeness (PI) profiles suggest that the molecular markers tested lack the power to resolve relationships at the base of the Doridina, while being more robust at family-level classifications. This supports the lack of consistent resolution between the 19 families within the Doridina across all three alignments. Most of the 19 families were recovered as monophyletic, and instances of non-monophyletic families were consistently recovered between analyses. We conclude that the alignment of variable regions has some effect on phylogenetic estimates of the Doridina, but these effects can vary depending on the size and scope of the phylogenetic query and PI of molecular markers.