We present a taxonomic revision of the black salamander (Aneides flavipunctatus) complex of northwestern California and extreme southeastern Oregon. The revision is based on a number of published works as well as new molecular and morphological data presented herein. The subspecies Aneides flavipunctatus niger Myers & Maslin 1948 is raised in rank to a full species. It is isolated far to the south of the main range on the San Francisco Peninsula, south and west of San Francisco Bay. Another geographically isolated set of populations occurs well inland in Shasta County, northern CA, mainly in the vicinity of Shasta Lake. It is raised from synonymy and recognized as Aneides iecanus (Cope 1883). The remaining taxa occur mainly along and inland from the coast from the vicinity of the Russian River and Lake Berryessa/Putah Creek, north to the vicinity of the Smith River near the Oregon border and more inland along the Klamath and Trinity Rivers and tributaries into Oregon. The northern segment of this nearly continuous range is named Aneides klamathensis Reilly and Wake 2019. We use molecular data to provide a detailed examination of a narrow contact zone between the northern A. klamathensis and the more southern A. flavipunctatus in southern Humboldt County in the vicinity of the Van Duzen and main fork of the Eel rivers. To the south is the remnant of the former species and it takes the name Aneides flavipunctatus (Strauch 1870). It is highly diversified morphologically and genetically and requires additional study.

Abstract The ring species Ensatina represents a classic example of locally adapted lineages. The Monterey ensatina (Ensatina eschscholtzii eschscholtzii) is a cryptic subspecies with brown coloration, although a recently discovered polymorphic population within a wind-blown sand region also contains leucistic (pink) and xanthistic (orange) morphs. In the present study, the frequency of leucism/xanthism was mapped across the subspecies’ range, revealing that these morphs are generally rare or absent except within regions containing light-coloured substrate. Attack rates were estimated using clay models of the three morphs, deployed only at the crepuscular period and during the night, on both light and dark substrates at a site within the dune sand region. Model selection found that the interaction between morph and substrate colour best predicted attack rates. Typical morphs had equal attack rates on both substrates while xanthistic and leucistic morphs incurred significantly fewer attacks on light vs. dark substrate, and there was no significant difference in attack rates among morphs on light substrates. These results support the idea that xanthistic and leucistic morphs are poorly adapted for dark substrates compared to typical morphs, but they are more or less equally adapted for light substrates. We suggest that this microgeographic island of relaxed selection on light-coloured morphs helps to explain the existence of this polymorphic population.

Aim: Phylogeographical structure in the black salamander (Aneides flavipunctatus) was inferred using two independent genetic datasets. Concordance between the datasets was sought in order to evaluate earlier suggestions of species-level breaks and evidence of vicariance and long-term isolation within the complex. We hypothesized that major phylogeographical breaks would either correspond to current tectonic plate boundaries or to historical geological processes. Location: North-western California and southern Oregon (USA). Methods: Three mitochondrial DNA (mtDNA) genes were sequenced for 240 black salamanders from 136 localities, and up to 13 nuclear DNA loci were sequenced for 145 black salamanders representing 93 localities. Phylogenetic analysis of our mitochondrial dataset was performed to recover major lineages, while spatial clustering analysis of our nuclear dataset was utilized to identify points of concordance with our mtDNA phylogeny. Levels of gene flow were estimated for all contact zones. Results: Our mitochondrial phylogeny recovered four major lineages. Cluster analysis of our nuclear dataset is consistent with a four-population scheme, with the boundaries matching those of the mtDNA lineages. Gene flow across a contact zone in southern Humboldt between two of the populational units is extremely limited (2Nm < 1). In what is identified as the Central Core population, two distinctive subpopulations were delineated based on nuclear data, but mitochondrial data are discordant. Main conclusions: The Aneides flavipunctatus complex comprises at least four species-level units. Two of the boundaries between these units are associated with current tectonic plate boundaries. The contact zone between our Northwest and Central Core populations lies adjacent to the Mendocino Triple Junction (MTJ), where the Humboldt and North American plates meet, while the area separating the Santa Cruz and Central Core populations corresponds to the boundary between the Pacific and Humboldt tectonic plates. The phylogeographical break within the Central Core population lies in a region in which uplift occurred that is associated with the historical position of the migrating MTJ.

Background

Phylogeography is an important tool that can be used to reveal cryptic biodiversity and to better understand the processes that promote lineage diversification. We studied the phylogeographic history of the Arboreal Salamander (Aneides lugubris), a wide-ranging species endemic to the California floristic province. We used multi-locus data to reconstruct the evolutionary history of A. lugubris and to discover the geographic location of major genetic breaks within the species. We also used species distribution modeling and comparative phylogeography to better understand the environmental factors that have shaped the genetic history of A. lugubris.

Results

We found six major mitochondrial clades in A. lugubris. Nuclear loci supported the existence of at least three genetically distinct groups, corresponding to populations north of the San Francisco Bay and in the Sierra Nevada, in the Santa Cruz Mountains, and in the central coast and southern California. All of the genetic breaks in mitochondrial and nuclear loci corresponded to regions where historical barriers to dispersal have been observed in other species. Geologic or water barriers likely were the most important factors restricting gene flow among clades. Climatic unsuitability during glacial maximum may have contributed to the isolation of the mitochondrial clades in the central coast and southern California. A projection of our species distribution model to a future scenario with a moderate amount of climate change suggests that most of the range of A. lugubris will remain climatically suitable, but climatic conditions in the Sierra Nevada and low elevation areas in Southern California are likely to deteriorate.

Conclusions

Aneides lugubris contains substantial cryptic genetic diversity as a result of historical isolation of populations. At least two (and perhaps three) evolutionarily significant units in A. lugubris merit protection; all six mitochondrial clades should be considered as management units within the species.

(a) A map of Wallacea showing islands invaded by Duttaphrynus melanostictus in red, islands inhabited by Varanus komodoensis in blue, and localities of genetic samples in yellow points. (b) A D. melanostictus from Lombok Island. (c) Environmental niche model for the Sunda Islands clade of D. melanostictus projected into Wallacea. Green color indicates very high suitability, yellow color indicates high suitability, and orange color indicates moderate suitability.

The Lesser Sundas Archipelago is comprised of two parallel chains of islands that extend between the Asian continental shelf (Sundaland) and Australo-Papuan continental shelf (Sahul). These islands have served as stepping stones for taxa dispersing between the Asian and Australo-Papuan biogeographical realms. While the oceanic barriers have prevented many species from colonizing the archipelago, a number of terrestrial vertebrate species have colonized the islands either by rafting/swimming or by human introduction. Here, we examine phylogeographic structure within the Lesser Sundas for three snake, two lizard and two frog species that each has a Sunda Shelf origin. These species are suspected to have recently colonized the archipelago, though all have inhabited the Lesser Sundas for over 100 years. We sequenced mtDNA from 231 samples to test whether there is sufficiently deep genetic structure within any of these taxa to reject human-mediated introduction. Additionally, we tested for genetic signatures of population expansion consistent with recent introduction and estimated the ages of Lesser Sundas clades, if any exist. Our results show little to no genetic structure between populations on different islands in five species and moderate structure in two species. Nucleotide diversity is low for all species, and the ages of the most recent common ancestor for species with monophyletic Lesser Sundas lineages date to the Holocene or late Pleistocene. These results support the hypothesis that these species entered the archipelago relatively recently and either naturally colonized or were introduced by humans to most of the larger islands in the archipelago within a short time span.

Abstract: The Lesser Sunda Archipelago consists of hundreds of oceanic islands located in southern Wallacea. The Sunda ratsnake, Coelognathus subradiatus, is endemic to the Lesser Sundas and is found on most of the major islands. Mitochondrial DNA was sequenced from snakes representing five of the major islands revealing that levels of sequence divergence between islands range from 2-7%. Phylogenetic analyses recover what can be interpreted as a three-lineage polytomy consisting of lineages from 1) Alor, 2) Sumbawa + Flores, and 3) Timor + Wetar. The archipelago was colonized from the Sunda Shelf ∼7 Myr with subsequent population divergence occurring ∼4.5 Myr, likely resulting in insular species formation.

Aim: The Lesser Sunda Islands are situated between the Sunda and Sahul Shelves, with a linear arrangement that has functioned as a two-way filter for taxa dispersing between the Asian and Australo-Papuan biogeographical realms. Distributional patterns of many terrestrial vertebrates suggest a stepping-stone model of island colonization. Here we investigate the timing and sequence of island colonization in Asian-origin fanged frogs from the volcanic Sunda Arc islands with the goal of testing the stepping-stone model of island colonization. Location: The Indonesian islands of Java, Lombok, Sumbawa, Flores and Lembata. Taxon: Limnonectes dammermani and L. kadarsani (Family: Dicroglossidae). Methods: Mitochondrial DNA was sequenced from 153 frogs to identify major lineages and to select samples for an exon-capture experiment. We designed probes to capture sequence data from 974 exonic loci (1,235,981 bp) from 48 frogs including the outgroup species, L. microdiscus. The resulting data were analysed using phylogenetic, population genetic and biogeographical model testing methods. Results: The mtDNA phylogeny finds L. kadarsani paraphyletic with respect to L. dammermani, with a pectinate topology consistent with the stepping-stone model. Phylogenomic analyses of 974 exons recovered the two species as monophyletic sister taxa that diverged ~7.6 Ma with no detectable contemporary gene flow, suggesting introgression of the L. dammermani mitochondrion into L. kadarsani on Lombok resulting from an isolated ancient hybridization event ~4 Ma. Within L. kadarsani, the Lombok lineage diverged first while the Sumbawa and Lembata lineages are nested within a Flores assemblage composed of two parapatrically distributed lineages meeting in central Flores. Biogeographical model comparison found strict stepping-stone dispersal to be less likely than models involving leap-frog dispersal events. Main conclusions: These results suggest that the currently accepted stepping-stone model of island colonization might not best explain the current patterns of diversity in the archipelago. The high degree of genetic structure, large divergence times, and absent or low levels of migration between lineages suggests that L. kadarsani represents five distinct species.

Bent-toed Geckos, genus Cyrtodactylus, are one of the most diverse terrestrial vertebrate groups, and their range extends from South Asia into Australo-Papua and adjacent Pacific islands. Given the generally high faunal endemism on Wallacean islands, it is rather paradoxical that the diversity in these geckos appears to be so low (21 species in Wallacea, 15 in the Philippines) compared with continental shelf assemblages (>300 species on Sunda + Sahul Shelves + adjacent islands). To determine whether this shortfall was real or an artifact of historical undersampling, we analyzed mitochondrial DNA sequences of hundreds of southern Wallacean samples (Lesser Sundas + southern Maluku). After screening to guide sample selection for target capture data collection, we obtained a 1150-locus genomic dataset (1,476,505 bp) for 119 samples of southern Wallacean and closely related lineages. The results suggest that species diversity of Cyrtodactylus in southern Wallacea is vastly underestimated, with phylogenomic and clustering analyses suggesting as many as 25 candidate species, in contrast to the 8 currently described. Gene exchange between adjacent candidate species is absent or minimal across the archipelago with only one case of > 0.5 migrants per generation. Biogeographical analysis suggests that the hitherto unrecognized diversity is the result of at least three independent dispersals from Sulawesi or its offshore islands into southern Wallacea between 6 and 14 Ma, with one invasion producing small-bodied geckos and the other two or three producing larger-bodied geckos. The smaller-bodied laevigatus group appears to be able to coexist with members of either larger-bodied clade, but we have yet to find members of the two larger-bodied clades occurring in sympatry, suggesting that ecological partitioning or competitive exclusion may be shaping individual island assemblages.