Contrasting Phylogeography of Sandy vs. Rocky Supralittoral Isopods in the Megadiverse and Geologically Dynamic Gulf of California and Adjacent Areas

Phylogeographic studies of animals with low vagility and restricted to patchy habitats of the supralittoral zone, can uncover unknown diversity and shed light on processes that shaped evolution along a continent’s edge. The Pacific coast between southern California and central Mexico, including the megadiverse Gulf of California, offers a remarkable setting to study biological diversification in the supralittoral. A complex geological history coupled with cyclical fluctuations in temperature and sea level provided ample opportunities for diversification of supralittoral organisms. Indeed, a previous phylogeographic study of Ligia, a supralittoral isopod that has limited dispersal abilities and is restricted to rocky patches, revealed high levels of morphologically cryptic diversity. Herein, we examined phylogeographic patterns of Tylos, another supralittoral isopod with limited dispersal potential, but whose habitat (i.e., sandy shores) appears to be more extensive and connected than that of Ligia. We conducted Maximum Likelihood and Bayesian phylogenetic analyses on mitochondrial and nuclear DNA sequences. These analyses revealed multiple highly divergent lineages with discrete regional distributions, despite the recognition of a single valid species for this region. A traditional species-diagnostic morphological trait distinguished several of these lineages. The phylogeographic patterns of Tylos inside the Gulf of California show a deep and complex history. In contrast, patterns along the Pacific region between southern California and the Baja Peninsula indicate a recent range expansion, probably postglacial and related to changes in sea surface temperature (SST). In general, the phylogeographic patterns of Tylos differed from those of Ligia. Differences in the extension and connectivity of the habitats occupied by Tylos and Ligia may account for the different degrees of population isolation experienced by these two isopods and their contrasting phylogeographic patterns. Identification of divergent lineages of Tylos in the study area is important for conservation, as some populations are threatened by human activities.     

Multiple transisthmian divergences,extensive cryptic diversity,occasional long‐distance dispersal,and biogeographic patterns in a marine coastal isopod with an amphi‐American distribution

Excirolana braziliensis is a coastal intertidal isopod with a broad distribution spanning the Atlantic and Pacific tropical and temperate coasts of the American continent. Two separate regional studies (one in Panama and one in Chile) revealed the presence of highly genetically divergent lineages, implying that this taxon constitutes a cryptic species complex. The relationships among the lineages found in these two different regions and in the rest of the distribution, however, remain unknown. To better understand the phylogeographic patterns of E. braziliensis, we conducted phylogenetic analyses of specimens from much of its entire range. We obtained DNA sequences for fragments of four mitochondrial genes (16S rDNA, 12S rDNA, COI, and Cytb) and also used publicly available sequences. We conducted maximum likelihood and Bayesian phylogenetic reconstruction methods. Phylogeographic patterns revealed the following: (1) new highly divergent lineages of E. braziliensis; (2) three instances of Atlantic–Pacific divergences, some of which appear to predate the closure of the Isthmus of Panama; (3) the distributional limit of highly divergent lineages found in Brazil coincides with the boundary between two major marine coastal provinces; (4) evidence of recent long‐distance dispersal in the Caribbean; and (5) populations in the Gulf of California have closer affinities with lineages further south in the Pacific, which contrasts with the closer affinity with the Caribbean reported for other intertidal organisms. The high levels of cryptic diversity detected also bring about challenges for the conservation of this isopod and its fragile environment, the sandy shores. Our findings underscore the importance of comprehensive geographic sampling for phylogeographic and taxonomical studies of broadly distributed putative species harboring extensive cryptic diversity.     

Patterns of morphological variation among populations of the widespread annual killifish Nothobranchius orthonotus are independent of genetic divergence and biogeography

Populations of annual killifish of the genus Nothobranchius occur in patchily distributed temporary pools in the East African savannah. Their fragmented distribution and low dispersal ability result in highly structured genetic clustering of their populations. In this study, we examined body shape variation in a widely distributed species, Nothobranchius orthonotus with known phylogeographic structure. We tested whether genetic divergence of major mitochondrial lineages forming two candidate species is congruent with phenotypic diversification, using linear and geometric morphometry analyses of body shape in 23 wild populations. We also conducted a common‐garden experiment with two wild‐derived populations to control for the effect of local environmental conditions on body shape. We identified different allometric trajectories for different mitochondrial lineages and candidate species in both sexes. However, in a principal components analysis of population‐level body shape, the separation among mitochondrial lineages was incomplete. Higher similarity of mitochondrial lineages belonging to different candidate species than that of same candidate species prevented distinction of the two candidate species on the basis of body shape. Analysis at the individual level demonstrated that N. orthonotus express high intrapopulation variability, with major overlap among individuals from all populations. In conclusion, we suggest that N. orthonotus be considered as a single species with an extensive geographic range, strong population genetic structure and high morphological variability.     

Variation in Body Shape across Species and Populations in a Radiation of Diaptomid Copepods

Inter and intra-population variation in morphological traits, such as body size and shape, provides important insights into the ecological importance of individual natural populations. The radiation of Diaptomid species (~400 species) has apparently produced little morphological differentiation other than those in secondary sexual characteristics, suggesting sexual, rather than ecological, selection has driven speciation. This evolutionary history suggests that species, and conspecific populations, would be ecologically redundant but recent work found contrasting ecosystem effects among both species and populations. This study provides the first quantification of shape variation among species, populations, and/or sexes (beyond taxonomic illustrations and body size measurements) to gain insight into the ecological differentiation of Diaptomids. Here we quantify the shape of five Diaptomid species (family Diaptomidae) from four populations each, using morphometric landmarks on the prosome, urosome, and antennae. We partition morphological variation among species, populations, and sexes, and test for phenotype-by-environment correlations to reveal possible functional consequences of shape variation. We found that intraspecific variation was 18-35% as large as interspecific variation across all measured traits. Interspecific variation in body size and relative antennae length, the two traits showing significant sexual dimorphism, were correlated with lake size and geographic location suggesting some niche differentiation between species. Observed relationships between intraspecific morphological variation and the environment suggest that divergent selection in contrasting lakes might contribute to shape differences among local populations, but confirming this requires further analyses. Our results show that although Diaptomid species differ in their reproductive traits, they also differ in other morphological traits that might indicate ecological differences among species and populations.     

Phylogeographic patterns of a lower intertidal isopod in the Gulf of California and the Caribbean and comparison with other intertidal isopods

A growing body of knowledge on the diversity and evolution of intertidal isopods across different regions worldwide has enhanced our understanding on biological diversification at the poorly studied, yet vast, sea–land interface. High genetic divergences among numerous allopatric lineages have been identified within presumed single broadly distributed species. Excirolana mayana is an intertidal isopod that is commonly found in sandy beaches throughout the Gulf of California. Its distribution in the Pacific extends from this basin to Colombia and in the Atlantic from Florida to Venezuela. Despite its broad distribution and ecological importance, its evolutionary history has been largely neglected. Herein, we examined phylogeographic patterns of E. mayana in the Gulf of California and the Caribbean, based on maximum‐likelihood and Bayesian phylogenetic analyses of DNA sequences from four mitochondrial genes (16S rDNA, 12S rDNA, cytochrome oxidase I gene, and cytochrome b gene). We compared the phylogeographic patterns of E. mayana with those of the coastal isopods Ligia and Excirolana braziliensis (Gulf of California and Caribbean) and Tylos (Gulf of California). We found highly divergent lineages in both, the Gulf of California and Caribbean, suggesting the presence of multiple species. We identified two instances of Atlantic–Pacific divergences. Some geographical structuring among the major clades found in the Caribbean is observed. Haplotypes from the Gulf of California form a monophyletic group sister to a lineage found in Venezuela. Phylogeographic patterns of E. mayana in the Gulf of California differ from those observed in Ligia and Tylos in this region. Nonetheless, several clades of E. mayana have similar distributions to clades of these two other isopod taxa. The high levels of cryptic diversity detected in E. mayana also pose challenges for the conservation of this isopod and its fragile environment, the sandy shores.     

A new and highly divergent mitochondrial lineage in the Small Five-toed Jerboa,Allactaga elater,from Iran (Mammalia: Rodentia)

The Small Five-toed Jerboa, Allactaga elater, is a small rodent adapted to desert and semi-arid habitats with a widespread distribution around the Caucasus. Previous studies have suggested the occurrence of subspecific variation within the species but, except for a recent phylogeny of the genus Allactaga, most of the work done on the taxonomy of the group relies on morphological data only. To contribute to the current understanding of patterns of genetic diversity of A. elater we analysed one mitochondrial locus, cytochrome b, from 13 Iranian specimens. Comparing to a recent phylogeny, our results suggest the existence of two additional mitochondrial lineages, one that clusters within previously described lineages and a new and highly divergent one. The two novel mitochondrial lineages occur in the north and form two highly divergent monophyletic groups (Dxy = 14%), which likely separated during the Pleistocene.     

A Complex Evolutionary History in a Remote Archipelago: Phylogeography and Morphometrics of the Hawaiian Endemic Ligia Isopods

Compared to the striking diversification and levels of endemism observed in many terrestrial groups within the Hawaiian Archipelago, marine invertebrates exhibit remarkably lower rates of endemism and diversification. Supralittoral invertebrates restricted to specific coastal patchy habitats, however, have the potential for high levels of allopatric diversification. This is the case of Ligia isopods endemic to the Hawaiian Archipelago, which most likely arose from a rocky supralittoral ancestor that colonized the archipelago via rafting, and diversified into rocky supralittoral and inland lineages. A previous study on populations of this isopod from Oʻahu and Kauaʻi revealed high levels of allopatric differentiation, and suggested inter-island historical dispersal events have been rare. To gain a better understanding on the diversity and evolution of this group, we expanded prior phylogeographic work by incorporating populations from unsampled main Hawaiian Islands (Maui, Molokaʻi, Lanaʻi, and Hawaiʻi), increasing the number of gene markers (four mitochondrial and two nuclear genes), and conducting Maximum likelihood and Bayesian phylogenetic analyses. Our study revealed new lineages and expanded the distribution range of several lineages. The phylogeographic patterns of Ligia in the study area are complex, with Hawaiʻi, Oʻahu, and the Maui-Nui islands sharing major lineages, implying multiple inter-island historical dispersal events. In contrast, the oldest and most geographically distant of the major islands (Kauaʻi) shares no lineages with the other islands. Our results did not support the monophyly of all the supralittoral lineages (currently grouped into L. hawaiensis), or the monophyly of the terrestrial lineages (currently grouped into L. perkinsi), implying more than one evolutionary transition between coastal and inland forms. Geometric-morphometric analyses of three supralittoral clades revealed significant body shape differences among them. A taxonomic revision of Hawaiian Ligia is warranted. Our results are relevant for the protection of biodiversity found in an environment subject to high pressure from disturbances.     

Geometric morphometric and phylogenetic analyses of Arizona Sky Island populations of Scaphinotus petersi Roeschke (Coleoptera: Carabidae)

Scaphinotus petersi Roeschke, 1907 (Carabidae) is a ground beetle endemic to Sky Islands in south‐eastern Arizona. Previous taxonomic studies described several subspecies with morphological differences inhabiting geographically isolated mountain ranges. We combined molecular sequence data and morphometric data, especially head and pronotum shape analyses, to examine the variation and divergence in subspecies and isolated montane populations. In this study, we employ a combination of distance morphometrics as well as geometric morphometrics to quantify the level of morphological variation, and to test the hypothesis that geographically distinct populations of S. petersi are phenotypically distinct. Results suggest that these isolated populations have diverged morphologically and genetically. Phylogenetic analyses identified two monophyletic lineages within the species that correspond generally to pronotum shape. We observed significant morphological variation among most montane populations in of S. petersi, with the pronotum shape as the clearest delimiting trait. © 2015 The Linnean Society of London     

Predictors of body shape among populations of a stream fish (Cyprinella venusta,Cypriniformes: Cyprinidae)

Performance‐related variation in fitness can manifest as morphological responses to ecological and evolutionary pressures. Eco‐morphological studies often utilize stark binary comparisons, such as lentic to lotic populations of freshwater fishes, to characterize relationships between form and function despite possible complications from confounding factors. In the present study, we compared body shape variation among lotic populations of a stream fish (Cyprinella venusta Girard) to disentangle the influence of ecological and evolutionary drivers of phenotypic change. We assessed the extent to which body shape corresponded to three key environmental factors (mean channel velocity, mean discharge, and mean annual run‐off), phylogeny (mitochondrial DNA divergence), and body size (centroid size). We also examined relationships between these parameters and a fineness index, which is a measure of streamlining and morphological optimization for steady swimming performance. All three environmental variables had some explanatory power, although morphological characteristics were predominantly associated with variation in mean annual run‐off. Phylogeny was also a strong predictor of morphological variation, whereas body size had little predictive power. Populations experiencing higher mean annual run‐off exhibited a shorter base of the dorsal fin, a more slender body and caudal peduncle, a smaller head in both horizontal and vertical dimensions, and a more anterior placement of the eye. With some exceptions, such as variation in jaw length, differences in body shape associated with phylogenetic history were similar to those associated with run‐off. Notably, all clades exhibited parallel responses to variation in run‐off. Populations experiencing high mean annual run‐off approached a hydrodynamic optimum, suggesting a morphology optimized for steady swimming performance. In contrast to previous studies that emphasize the importance of average water velocity, the findings of the present study indicate that morphological variation among populations of stream fishes is tightly linked to more complex aspects of hydrology and evolutionary history. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, ●● , ●●–●●.     

DNA barcoding uncovers extensive cryptic diversity in the African long-fin tetra Bryconalestes longipinnis (Alestidae: Characiformes)

To investigate the presence of cryptic diversity in the African longfin-tetra Bryconalestes longipinnis, we employed DNA barcoding in a phylogeographic context, as well as geometric morphometrics, documenting for the first time genetic and body shape variation in the species. Analysis of cytochrome oxidase I gene (coI) sequence variation exposed extremely high levels of genetic differentiation among samples from across the geographic range of the species (up to 18%), certainly much greater than the traditionally employed c. 3% sequence divergence heuristic threshold for conspecifics. Phylogeographic analyses of coI data revealed eight clusters/clades that diverge by >4% and up to 18% (p-distance), potentially representing cryptic members of a species complex. A clear biogeographic pattern was also uncovered, in which the two main coI lineages corresponded geographically with the upper Guinea (UG) and lower Guinea (LG) ichthyofaunal provinces of continental Africa, respectively. Within each of these main lineages, however, no apparent phylogeographic structuring was found. Despite strong genetic differentiation, there is considerable overlap in body shape variation between UG and LG populations. For the most part, morphological variation does not match the strength of the molecular phylogeographic signal. Therefore, the ability to reliably utilise external body shape for regional delimitation remains elusive. Further anatomical investigation appears necessary to establish whether compelling diagnostic morphological features do exist between the divergent lineages of the B. longipinnis complex uncovered in this study.     

Cryptic lineage diversity,body size divergence,and sympatry in a species complex of Australian lizards (Gehyra)

Understanding the joint evolutionary and ecological underpinnings of sympatry among close relatives remains a key challenge in biology. This problem can be addressed through joint phylogenomic and phenotypic analysis of complexes of closely related lineages within, and across, species and hence representing the speciation continuum. For a complex of tropical geckos from northern Australia—Gehyra nana and close relatives—we combine mtDNA phylogeography, exon‐capture sequencing, and morphological data to resolve independently evolving lineages and infer their divergence history and patterns of morphological evolution. Gehyra nana is found to include nine divergent lineages and is paraphyletic with four other species from the Kimberley region of north‐west Australia. Across these 13 taxa, 12 of which are restricted to rocky habitats, several lineages overlap geographically, including on the diverse Kimberley islands. Morphological evolution is dominated by body size shifts, and both body size and shape have evolved gradually across the group. However, larger body size shifts are observed among overlapping taxa than among closely related parapatric lineages of G. nana, and sympatric lineages are more divergent than expected at random. Whether elevated body size differences among sympatric lineages are due to ecological sorting or character displacement remains to be determined.     

Phenotypic variation among populations of Atherinops affinis(Atherinopsidae) with insights from a geometric morphometric analysis

Morphological character variation was examined in Atherinops affinis , a temperate marine silverside with a broad geographic range and presumed limited powers of dispersal. Populations of this species were sampled from three California mainland sites, one Channel Island site and one site in the upper Gulf of California. A geometric morphometric analysis yielded higher resolution in the assessment of phenotypic divergence among the four Pacific coast populations than either body measurement or meristic analysis, and it showed that most of the shape variation among these populations occurs in the head region and body depth of the fish. All three analyses supported the hypothesis that populations of A. affinis from central and southern California coastal waters and from Santa Catalina Island are morphologically distinct from each other; the Santa Catalina Island population was found to be the most divergent. On the basis of meristic characters alone, the population of A. affinis from the upper Gulf of California was different from A. affinis populations along the Pacific coast of California. The analyses revealed variation in several morphological characters, e.g . body depth and meristics, known to vary in association with environmental conditions. Given that A. affinis appears to have low among‐population genetic variation, this species may be phenotypically plastic in response to the environmental conditions of the habitat of each population.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Morphological divergence of lake and stream Phoxinus of Northern Italy and the Danube basin based on geometric morphometric analysis

Minnows of the genus Phoxinus are promising candidates to investigate adaptive divergence, as they inhabit both still and running waters of a variety of altitudes and climatic zones in Europe. We used landmark‐based geometric morphometric methods to quantify the level of morphological variability in Phoxinus populations from streams and lakes of Northern Italy and the Danube basin. We analyzed body shape differences of populations in the dorsal, lateral, and ventral planes, using a large array of landmarks and semilandmarks. As the species identification of Phoxinus on morphological characters is ambiguous, we used two mitochondrial genes to determine the genetic background of the samples and to ensure we are comparing homogenous groups. We have found significant body shape differences between habitats: Minnow populations inhabiting streams had a deeper body and caudal peduncle and more laterally inserted pectoral fins than minnows inhabiting lakes. We have also found significant body shape differences between genetic groups: Italian minnows had deeper bodies, deeper and shorter caudal peduncles, and a shorter and wider gape than both groups from the Danube. Our results show that the morphology of Phoxinus is highly influenced by habitat and that body shape variation between habitats was within the same range as between genetic groups. These morphological differences are possibly linked to different modes of swimming and foraging in the respective habitats and are likely results of phenotypic plasticity. However, differences in shape and interlandmark distances between the groups suggest that some (though few) morphometric characters might be useful for separating Phoxinus species.     

Cryptic diversity in disjunct populations of Middle American Montane Pitvipers: a systematic reassessment of Cerrophidion godmani

Jadin, R.C., Townsend, J.H., Castoe, T.A. & Campbell, J.A. (2012). Cryptic diversity in disjunct populations of Middle American Montane Pitvipers: a systematic reassessment of Cerrophidion godmani. —Zoologica Scripta, 41, 455–470. The discovery and taxonomic recognition of cryptic species has become increasingly frequent with the application of molecular phylogenetic analyses, particularly for species with broad geographic distributions. In this study we focus on the venomous pitviper species Cerrophidion godmani that is widely distributed throughout the highlands of Central America. We provide evidence based on both molecular phylogenetic analyses and morphological data that C. godmani represents three deeply divergent lineages and is possibly non‐monophyletic. These three lineages are relatively conserved in their morphology and tend to be highly variable among individuals, but we do find sufficient morphological characters to diagnose them as evolutionarily distinct. We apply these data, together with known geographic distributions of populations, to infer boundaries of these three divergent evolutionary lineages. Based on the body of evidence, we formally name and describe two new species of Cerrophidion and redescribe C. godmani sensu stricto.     

Evolution of thermal physiology and growth rate between populations of the western fence lizard (Sceloporus occidentalis)

Summary Hatchling Sceloporus occidentalis from northern populations (central Oregon) grow more slowly than hatchlings from southern populations (southern California) in nature. In this study, I determine whether this difference in growth rate results from differences in thermal environment and/or in thermoregulatory behavior. To determine the degree to which the thermal environment affects growth rate among populations, I reared hatchings from the northern and southern populations in a cycling thermal regime in one of three experimental treatments differing in access to radiant heat (6, 9, or 12 h radiant heat; remainder of 24 h at 15°C). I also measured the body temperature that each individual voluntarily selected over the course of the daily activity cycle. Growth rate varied positively with duration of access to radiant heat. Within the three treatments, individual growth rate was positively correlated with body temperature. Moreover, the difference in growth rate between the northern and southern populations was due in part to differences in behavior — individuals from northern populations selected lower body temperatures. I found that significant variation in body temperature was associated with family membership, suggesting that thermal physiology has a genetic basis. Moreover, growth rate was correlated with body temperature among families in each population suggesting a genetic correlation underlies the phenotypic correlations. Thus, genetically based variation in thermal physiology contributes to differences in growth rate among individuals within a population as well as to differences among populations.     

Adaptive responses and invasion: the role of plasticity and evolution in snail shell morphology

Invasive species often exhibit either evolved or plastic adaptations in response to spatially varying environmental conditions. We investigated whether evolved or plastic adaptation was driving variation in shell morphology among invasive populations of the New Zealand mud snail (Potamopyrgus antipodarum) in the western United States. We found that invasive populations exhibit considerable shell shape variation and inhabit a variety of flow velocity habitats. We investigated the importance of evolution and plasticity by examining variation in shell morphological traits 1) between the parental and F₁ generations for each population and 2) among populations of the first lab generation (F₁) in a common garden, full‐sib design using Canonical Variate Analyses (CVA). We compared the F₁ generation to the parental lineages and found significant differences in overall shell shape indicating a plastic response. However, when examining differences among the F₁ populations, we found that they maintained among‐population shell shape differences, indicating a genetic response. The F₁ generation exhibited a smaller shell morph more suited to the low‐flow common garden environment within a single generation. Our results suggest that phenotypic plasticity in conjunction with evolution may be driving variation in shell morphology of this widespread invasive snail.     

A non‐invasive geometric morphometrics method for exploring variation in dorsal head shape in urodeles: sexual dimorphism and geographic variation in Salamandra salamandra

The study of morphological variation among and within taxa can shed light on the evolution of phenotypic diversification. In the case of urodeles, the dorso‐ventral view of the head captures most of the ontogenetic and evolutionary variation of the entire head, which is a structure with a high potential for being a target of selection due to its relevance in ecological and social functions. Here, we describe a non‐invasive procedure of geometric morphometrics for exploring morphological variation in the external dorso‐ventral view of urodeles' head. To explore the accuracy of the method and its potential for describing morphological patterns we applied it to two populations of Salamandra salamandra gallaica from NW Iberia. Using landmark‐based geometric morphometrics, we detected differences in head shape between populations and sexes, and an allometric relationship between shape and size. We also determined that not all differences in head shape are due to size variation, suggesting intrinsic shape differences across sexes and populations. These morphological patterns had not been previously explored in S. salamandra , despite the high levels of intraspecific diversity within this species. The methodological procedure presented here allows to detect shape variation at a very fine scale, and solves the drawbacks of using cranial samples, thus increasing the possibilities of using collection specimens and alive animals for exploring dorsal head shape variation and its evolutionary and ecological implications in urodeles. J. Morphol. 278:475–485, 2017. © 2017 Wiley Periodicals, Inc.     

GEOGRAPHIC VARIATION AND PHENOTYPIC PLASTICITY OF NUMBER OF TRUNK VERTEBRAE IN SLENDER SALAMANDERS,BATRACHOSEPS (CAUDATA: PLETHODONTIDAE)

To understand the evolutionary significance of geographic variation, one must identify the factors that generate phenotypic differences among populations. I examined the causes of geographic variation in and evolutionary history of number of trunk vertebrae in slender salamanders, Batrachoseps (Caudata: Plethodontidae). Number of trunk vertebrae varies at many taxonomic levels within Batrachoseps. Parallel clines in number occur along an environmental gradient in three lineages in the Coast Ranges of California. These parallel clines may signal either adaptation or a shared phenotypically plastic response to the environmental gradient. By raising eggs from 10 populations representing four species of Batrachoseps, I demonstrated that number of trunk vertebrae can be altered by the developmental temperature; however, the degree of plasticity is insufficient to account for geographic variation. Thus, the geographic variation results largely from genetic variation. Number of trunk vertebrae covaries with body size and shape in diverse vertebrate taxa, including Batrachoseps. I hypothesize that selection for different degrees of elongation, possibly related to fossoriality, has led to the extensive evolution of number of trunk vertebrae in Batrachoseps. Analysis of intrapopulational variation revealed sexual dimorphism in both body shape and number of trunk vertebrae, but no correlation between these variables in either sex. Females are more elongate than males, a pattern that has been attributed to fecundity selection in other taxa. Patterns of covariation among different classes of vertebrae suggest that some intrapopulational variation in number results from changes in vertebral identity rather than changes in segmentation.