River drainage reorganization and reticulate evolution in the Two-Lined Salamander (Eurycea bislineata) species complex.

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge". Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for >100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical.

Biofluorescent sexual dimorphism revealed in a southern Appalachian endemic salamander, Plethodon metcalfi.

Biofluorescence occurs when a living organism absorbs high energy light and reemits it at longer wavelengths. Many species within clades of vertebrates are known to fluoresce including mammals, reptiles, birds, and fish. Most, if not all, amphibians exhibit biofluorescence when exposed to either blue (440-460 nm) or ultra-violet (360-380 nm) wavelengths of light. Salamanders (Lissamphibia: Caudata) appear to consistently fluoresce in green wavelengths (520-560 nm) when excited by blue light. Biofluorescence is theorized to have many ecological functions including mate signaling, camouflage, and mimicry. Despite the discovery of their biofluorescence, its role in salamander ecology and behavior remains unresolved. In this study we present the first case of biofluorescent sexual dimorphism within Amphibia and the first documentation of the biofluorescent pattern of a salamander within the Plethodon jordani species complex. This sexually dimorphic trait was discovered in the southern Appalachian endemic species, Southern Gray-Cheeked Salamander (Plethodon metcalfi, Brimley in Proc Biol Soc Wash 25:135-140, 1912), and may extend into other species within the Plethodon jordani and Plethodon glutinosus species complexes. We propose that this sexually dimorphic trait could be related to fluorescence of ventral modified granular glands used in plethodontid chemosensory communication.

Breeding behaviour predicts patterns of natural hybridization in North American minnows (Cyprinidae).

Premating barriers such as variation in reproductive behaviour can evolve quickly, but because gametic and postzygotic incompatibilities often evolve more slowly, circumstances that bring gametes into contact can breach the boundaries of premating isolation. In aquatic environments, the gametes of organisms with external fertilization are released into a constantly moving environment and may come into contact with heterospecific gametes. In fishes, nest association (spawning in another species' nest) is a behaviour that brings gametes from different species into close spatiotemporal proximity. These interactions might increase chances of hybridization, especially when multiple species associate with a single nest builder. This study addresses these interactions in the largest clade of North American freshwater fishes, the minnows (Cyprinidae). We compiled a list of over 17,000 hybrid specimens in conjunction with species distribution data, breeding behaviours, and an inferred phylogeny to test if breeding behaviour, in addition to evolutionary history, is an important predictor of hybridization. We find that breeding behaviour is a significant predictor of hybridization, even when phylogenetic relatedness and divergence time are accounted for. Specifically, nest associates are more likely to hybridize with other nest associates whereas non-nesting species had relatively low rates of hybridization.

Extensive Cryptic Diversity Within the Physalaemus cuvieri-Physalaemus ephippifer Species Complex (Amphibia, Anura) Revealed by Cytogenetic, Mitochondrial, and Genomic Markers.

Previous cytogenetic and phylogenetic analyses showed a high variability in the frog taxa Physalaemus cuvieri and Physalaemus ephippifer and suggested the presence of undescribed diversity in this species complex. Here, by 1) adding specimens from the Brazilian Amazon region, 2) employing sequence-based species delimitation approaches, and 3) including RADseq-style markers, we demonstrate that the diversity in the P. cuvieri-P. ephippifer species complex is even greater than previously suspected. Specimens from Viruá and Western Pará, located at the Guiana Amazonian area of endemism, were recovered as distinct from all previously identified lineages by the phylogenetic analyses based on mitochondrial DNA and RAD markers, a PCA from RAD data, and cytogenetic analysis. The sequence-based species delimitation analyses supported the recognition of one or two undescribed species among these Amazonian specimens and also supported the recognition of at least three other species in the P. cuvieri-P. ephippifer species complex. These new results reinforce the need for a comprehensive taxonomic revision.

Morphological Polymorphism Associated with Alternative Reproductive Tactics in a Plethodontid Salamander.

Understanding polymorphism is a central problem in evolution and ecology, and alternative reproductive tactics (ARTs) provide compelling examples for studying the origin and maintenance of behavioral and morphological variation. Much attention has been given to examples where "parasitic" individuals exploit the reproductive investment of "bourgeois" individuals, but some ARTs are instead maintained by environmental heterogeneity, with alternative tactics exhibiting differential fitness in discontinuous reproductive niches. We use genomic, behavioral, karyological, and field observational data to demonstrate one such example in plethodontid salamanders. These ARTs ("searching" and "guarding" males) are associated with different reproductive niches and, unlike most other examples in amphibians, demonstrate substantial morphological differences and inflexibility within a reproductive season. Evidence suggests the existence of these ARTs within three putative species in the two-lined salamander (Eurycea bislineata) species complex, with other members of this clade fixed for one of the two tactics. We highlight directions for future research in this system, including the relationship between these ARTs and parental care.

A hierarchical Bayesian model to incorporate uncertainty into methods for diversity partitioning.

Recently there have been major theoretical advances in the quantification and partitioning of diversity within and among communities, regions, and ecosystems. However, applying those advances to real data remains a challenge. Ecologists often end up describing their samples rather than estimating the diversity components of an underlying study system, and existing approaches do not easily provide statistical frameworks for testing ecological questions. Here we offer one avenue to do all of the above using a hierarchical Bayesian approach. We estimate posterior distributions of the underlying "true" relative abundances of each species within each unit sampled. These posterior estimates of relative abundance can then be used with existing formulae to estimate and partition diversity. The result is a posterior distribution of diversity metrics describing our knowledge (or beliefs) about the study system. This approach intuitively leads to statistical inferences addressing biologically motivated hypotheses via Bayesian model comparison. Using simulations, we demonstrate that our approach does as well or better at approximating the "true" diversity of a community relative to naïve or ad-hoc bias-corrected estimates. Moreover, model comparison correctly distinguishes between alternative hypotheses about the distribution of diversity within and among samples. Finally, we use an empirical ecological dataset to illustrate how the approach can be used to address questions about the makeup and diversities of assemblages at local and regional scales.

Pairwise beta diversity resolves an underappreciated source of confusion in calculating species turnover.

Beta diversity is an important metric in ecology quantifying differentiation or disparity in composition among communities, ecosystems, or phenotypes. To compare systems with different sizes (N, number of units within a system), beta diversity is often converted to related indices such as turnover or local/regional differentiation. Here we use simulations to demonstrate that these naive measures of dissimilarity depend on sample size and design. We show that when N is the number of sampled units (e.g., quadrats) rather than the "true" number of communities in the system (if such exists), these differentiation measures are biased estimators. We propose using average pairwise dissimilarity as an intuitive solution. That is, instead of attempting to estimate an N-community measure, we advocate estimating the expected dissimilarity between any random pair of communities (or sampling units)-especially when the "true" N is unknown or undefined. Fortunately, measures of pairwise dissimilarity or overlap have been used in ecology for decades, and their properties are well known. Using the same simulations, we show that average pairwise metrics give consistent and unbiased estimates regardless of the number of survey units sampled. We advocate pairwise dissimilarity as a general standardization to ensure commensurability of different study systems.

Extending the Concept of Diversity Partitioning to Characterize Phenotypic Complexity.

Most components of an organism's phenotype can be viewed as the expression of multiple traits. Many of these traits operate as complexes, where multiple subsidiary parts function and evolve together. As trait complexity increases, so does the challenge of describing complexity in intuitive, biologically meaningful ways. Traditional multivariate analyses ignore the phenomenon of individual complexity and provide relatively abstract representations of variation among individuals. We suggest adopting well-known diversity indices from community ecology to describe phenotypic complexity as the diversity of distinct subsidiary components of a trait. Using a hierarchical framework, we illustrate how total trait diversity can be partitioned into within-individual complexity (α diversity) and between-individual components (ß diversity). This approach complements traditional multivariate analyses. The key innovations are (i) addition of individual complexity within the same framework as between-individual variation and (ii) a group-wise partitioning approach that complements traditional level-wise partitioning of diversity. The complexity-as-diversity approach has potential application in many fields, including physiological ecology, ecological and community genomics, and transcriptomics. We demonstrate the utility of this complexity-as-diversity approach with examples from chemical and microbial ecology. The examples illustrate biologically significant differences in complexity and diversity that standard analyses would not reveal.

iteRates: An R Package for Implementing a Parametric Rate Comparison on Phylogenetic Trees.

Patterns of diversification rate variation detected in phylogenetic hypotheses are frequently used to infer historical, ecological, and evolutionary processes. The parametric rate comparison (PRC) is a method for detecting rate variation in trees that models branch lengths as random variables drawn from familiar statistical distributions. iteRates is a library of functions for the R statistical computing environment for implementing PRC on phylogenetic trees. Here, we describe some of the functions in iteRates for subtree identification, tree manipulation, applying the PRC and K-clades PRC analyses, and conducting a whole-tree randomization test.

Symbiote transmission and maintenance of extra-genomic associations.

Symbiotes can be transmitted from parents to offspring or horizontally from unrelated hosts or the environment. A key question is whether symbiote transmission is similar enough to Mendelian gene transmission to generate and maintain coevolutionary associations between host and symbiote genes. Recent papers come to opposite conclusions, with some suggesting that any horizontal transmission eliminates genetic association. These studies are hard to compare owing to arbitrary differences in modeling approach, parameter values, and assumptions about selection. I show that associations between host and symbiote genes (extra-genomic associations) can be described by the same dynamic model as conventional linkage disequilibria between genes in the same genome. Thus, covariance between host and symbiote genomes depends on population history, geographic structure, selection, and co-transmission rate, just as covariance between genes within a genome. The conclusion that horizontal transmission rapidly erodes extra-genomic associations is equivalent to the conclusion that recombination rapidly erodes associations between genes within a genome. The conclusion is correct in the absence of population structure or selection. However, population structure can maintain spatial associations between host and symbiote traits, and non-additive selection (interspecific epistasis) can generate covariances between host and symbiote genotypes. These results can also be applied to cultural or other non-genetic traits. This work contributes to a growing consensus that genomic, symbiotic, and gene-culture evolution can be analyzed under a common theoretical framework. In terms of coevolutionary potential, symbiotes can be viewed as lying on a continuum between the intimacy of genes and the indifference of casually co-occurring species.

Similarity and differentiation between bacteria associated with skin of salamanders (Plethodon jordani) and free-living assemblages.

All animals and plants have intimate associations with microbes. Opinion has shifted from viewing microbes primarily as pathogens to the idea that healthy animals and plants carry specialized communities of coevolving microorganisms. However, the generality of this proposition is unknown because surveys rarely compare host-associated microbes with samples from relevant microhabitats. Symbiotic communities might be assembled from local environments with little evolutionary specialization. We evaluated the specificity of bacteria associated with salamander skin in comparison with surfaces in their immediate environments using 16S rRNA sequences. Host-associated and free-living samples were significantly different. However, relative abundances were strongly correlated; the most abundant taxa on salamander skin were also most abundant on moist debris on the forest floor. Thus, although bacterial assemblages on salamander skin are statistically differentiated from those on inanimate surfaces, they are not entirely 'distinct'. Candidate salamander specialists were few in number and occurred at low relative abundances. Within some OTUs, differences in allele frequency suggested genetic specialization at finer levels. Although host-associated and free-living assemblages were similar, a range of more or less specialized symbiotes was evident and bacteria on salamander skin were often specific genotypes of OTUs commonly found on other moist surfaces in the environment.

Alternative forms for genomic clines.

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory-based logit-logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one-to-one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package "HIest" for the R statistical software (CRAN; http://cran.r-project.org/).

Transgressive Hybrids as Hopeful Monsters.

The origin of novelty is a critical subject for evolutionary biologists. Early geneticists speculated about the sudden appearance of new species via special macromutations, epitomized by Goldschmidt's infamous "hopeful monster". Although these ideas were easily dismissed by the insights of the Modern Synthesis, a lingering fascination with the possibility of sudden, dramatic change has persisted. Recent work on hybridization and gene exchange suggests an underappreciated mechanism for the sudden appearance of evolutionary novelty that is entirely consistent with the principles of modern population genetics. Genetic recombination in hybrids can produce transgressive phenotypes, "monstrous" phenotypes beyond the range of parental populations. Transgressive phenotypes can be products of epistatic interactions or additive effects of multiple recombined loci. We compare several epistatic and additive models of transgressive segregation in hybrids and find that they are special cases of a general, classic quantitative genetic model. The Dobzhansky-Muller model predicts "hopeless" monsters, sterile and inviable transgressive phenotypes. The Bateson model predicts "hopeful" monsters with fitness greater than either parental population. The complementation model predicts both. Transgressive segregation after hybridization can rapidly produce novel phenotypes by recombining multiple loci simultaneously. Admixed populations will also produce many similar recombinant phenotypes at the same time, increasing the probability that recombinant "hopeful monsters" will establish true-breeding evolutionary lineages. Recombination is not the only (or even most common) process generating evolutionary novelty, but might be the most credible mechanism for sudden appearance of new forms.

Effects of climatic and geological processes during the pleistocene on the evolutionary history of the northern cavefish, Amblyopsis spelaea (teleostei: amblyopsidae).

Climatic and geological processes associated with glaciation cycles during the Pleistocene have been implicated in influencing patterns of genetic variation and promoting speciation of temperate flora and fauna. However, determining the factors promoting divergence and speciation is often difficult in many groups because of our limited understanding of potential vicariant barriers and connectivity between populations. Pleistocene glacial cycles are thought to have significantly influenced the distribution and diversity of subterranean invertebrates; however, impacts on subterranean aquatic vertebrates are less clear. We employed several hypothesis-driven approaches to assess the impacts of Pleistocene climatic and geological changes on the Northern Cavefish, Amblyopsis spelaea, whose current distribution occurs near the southern extent of glacial advances in North America. Our results show that the modern Ohio River has been a significant barrier to dispersal and is correlated with patterns of genetic divergence. We infer that populations were isolated in two refugia located north and south of the Ohio River during the most recent two glacial cycles with evidence of demographic expansion in the northern isolate. Finally, we conclude that climatic and geological processes have resulted in the formation of cryptic forms and advocate recognition of two distinct phylogenetic lineages currently recognized as A. spelaea.

Evidence for repeated loss of selective constraint in rhodopsin of amblyopsid cavefishes (Teleostei: Amblyopsidae).

The genetic mechanisms underlying regressive evolution-the degeneration or loss of a derived trait--are largely unknown, particularly for complex structures such as eyes in cave organisms. In several eyeless animals, the visual photoreceptor rhodopsin appears to have retained functional amino acid sequences. Hypotheses to explain apparent maintenance of function include weak selection for retention of light-sensing abilities and its pleiotropic roles in circadian rhythms and thermotaxis. In contrast, we show that there has been repeated loss of functional constraint of rhodopsin in amblyopsid cavefishes, as at least three cave lineages have independently accumulated unique loss-of-function mutations over the last 10.3 Mya. Although several cave lineages still possess functional rhodopsin, they exhibit increased rates of nonsynonymous mutations that have greater effect on the structure and function of rhodopsin compared to those in surface lineages. These results indicate that functionality of rhodopsin has been repeatedly lost in amblyopsid cavefishes. The presence of a functional copy of rhodopsin in some cave lineages is likely explained by stochastic accumulation of mutations following recent subterranean colonization.

Tests of two methods for identifying founder effects in metapopulations reveal substantial type II error.

Genetic analysis has been promoted as a way to reconstruct recent historical dynamics ("historical demography") by screening for signatures of events, such as bottlenecks, that disrupt equilibrium patterns of variation. Such analyses might also identify "metapopulation" processes like extinction and recolonization or source-sink dynamics, but this potential remains largely unrealized. Here we use simulations to test the ability of two currently used strategies to distinguish between a set of interconnected subpopulations (demes) that have undergone bottlenecks or extinction and recolonization events (metapopulation dynamics) from a set of static demes. The first strategy, decomposed pairwise regression, provides a holistic test for heterogeneity among demes in their patterns of isolation-by-distance. This method suffered from a type II error rate of 59-100 %, depending on parameter conditions. The second strategy tests for deviations from mutation-drift equilibrium on a deme-by-deme basis to identify sites likely to have experienced recent bottlenecks or founder effects. Although bottleneck tests have good statistical power for single populations with recent population declines, their validity in structured populations has been called into question, and they have not been tested in a metapopulation context with immigration (or colonization) and population recovery. Our simulations of hypothetical metapopulations show that population recovery can rapidly eliminate the statistical signature of a bottleneck, and that moderate levels of gene flow can generate a false signal of recent population growth for demes in equilibrium. Although we did not cover all possible metapopulation scenarios, the performance of the tests was disappointing. Our results indicate that these methods might often fail to identify population bottlenecks and founder effects if population recovery and/or gene flow are influential demographic features of the study system.

A parametric method for assessing diversification-rate variation in phylogenetic trees.

Phylogenetic hypotheses are frequently used to examine variation in rates of diversification across the history of a group. Patterns of diversification-rate variation can be used to infer underlying ecological and evolutionary processes responsible for patterns of cladogenesis. Most existing methods examine rate variation through time. Methods for examining differences in diversification among groups are more limited. Here, we present a new method, parametric rate comparison (PRC), that explicitly compares diversification rates among lineages in a tree using a variety of standard statistical distributions. PRC can identify subclades of the tree where diversification rates are at variance with the remainder of the tree. A randomization test can be used to evaluate how often such variance would appear by chance alone. The method also allows for comparison of diversification rate among a priori defined groups. Further, the application of the PRC method is not restricted to monophyletic groups. We examined the performance of PRC using simulated data, which showed that PRC has acceptable false-positive rates and statistical power to detect rate variation. We apply the PRC method to the well-studied radiation of North American Plethodon salamanders, and support the inference that the large-bodied Plethodon glutinosus clade has a higher historical rate of diversification compared to other Plethodon salamanders.

Lethal effects of water quality on threatened California salamanders but not on co-occurring hybrid salamanders.

Biological invasions and habitat alteration are often detrimental to native species, but their interactions are difficult to predict. Interbreeding between native and introduced species generates novel genotypes and phenotypes, and human land use alters habitat structure and chemistry. Both invasions and habitat alteration create new biological challenges and opportunities. In the intensively farmed Salinas Valley, California (U.S.A.), threatened California tiger salamanders (Ambystoma californiense) have been replaced by hybrids between California tiger salamander and introduced barred tiger salamanders (Ambystoma tigrinum mavortium). We conducted an enclosure experiment to examine the effects habitat modification and relative frequency of hybrid and native California tiger salamanders have on recruitment of salamanders and their prey, Pacific chorus frogs (Pseudacris regilla). We tested whether recruitment differed among genetic classes of tiger salamanders (hybrid or native) and pond hydroperiod (seasonal or perennial). Roughly 6 weeks into the experiment, 70% (of 378 total) of salamander larvae died in 4 out of 6 ponds. Native salamanders survived (n = 12) in these ponds only if they had metamorphosed prior to the die-offs. During die-offs, all larvae of native salamanders died, whereas 56% of hybrid larvae died. We necropsied native and hybrid salamanders, tested water quality, and queried the California Department of Pesticide Regulation database to investigate possible causes of the die-offs. Salamander die-offs, changes in the abundance of other community members (invertebrates, algae, and cyanobacteria), shifts in salamander sex ratio, and patterns of pesticide application in adjacent fields suggest that pesticide use may have contributed to die-offs. That all survivors were hybrids suggests that environmental stress may promote rapid displacement of native genotypes.

Estimating ancestry and heterozygosity of hybrids using molecular markers.

BACKGROUND: Hybridization, genetic mixture of distinct populations, gives rise to myriad recombinant genotypes. Characterizing the genomic composition of hybrids is critical for studies of hybrid zone dynamics, inheritance of traits, and consequences of hybridization for evolution and conservation. Hybrid genomes are often summarized either by an estimate of the proportion of alleles coming from each ancestral population or classification into discrete categories like F1, F2, backcross, or merely "hybrid" vs. "pure". In most cases, it is not realistic to classify individuals into the restricted set of classes produced in the first two generations of admixture. However, the continuous ancestry index misses an important dimension of the genotype. Joint consideration of ancestry together with interclass heterozygosity (proportion of loci with alleles from both ancestral populations) captures all of the information in the discrete classification without the unrealistic assumption that only two generations of admixture have transpired. METHODS: I describe a maximum likelihood method for joint estimation of ancestry and interclass heterozygosity. I present two worked examples illustrating the value of the approach for describing variation among hybrid populations and evaluating the validity of the assumption underlying discrete classification. RESULTS: Naively classifying natural hybrids into the standard six line cross categories can be misleading, and false classification can be a serious problem for datasets with few molecular markers. My analysis underscores previous work showing that many (50 or more) ancestry informative markers are needed to avoid erroneous classification. CONCLUSION: Although classification of hybrids might often be misleading, valuable inferences can be obtained by focusing directly on distributions of ancestry and heterozygosity. Estimating and visualizing the joint distribution of ancestry and interclass heterozygosity is an effective way to compare the genetic structure of hybrid populations and these estimates can be used in classic quantitative genetic methods for assessing additive, dominant, and epistatic genetic effects on hybrid phenotypes and fitness. The methods are implemented in a freely available package "HIest" for the R statistical software ( http://cran.r-project.org/web/packages/HIest/index.html).

Delimiting species using multilocus data: diagnosing cryptic diversity in the southern cavefish, Typhlichthys subterraneus (Teleostei: Amblyopsidae).

A major challenge facing biodiversity conservation and management is that a significant portion of species diversity remains undiscovered or undescribed. This is particularly evident in subterranean animals in which species delimitation based on morphology is difficult because differentiation is often obscured by phenotypic convergence. Multilocus genetic data constitute a valuable source of information for species delimitation in such organisms, but until recently, few methods were available to objectively test species delimitation hypotheses using genetic data. Here, we use recently developed methods for discovering and testing species boundaries and relationships using a multilocus dataset in a widely distributed subterranean teleost fish, Typhlichthys subterraneus, endemic to Eastern North America. We provide evidence that species diversity in T. subterraneus is currently underestimated and that the picture of a single, widely distributed species is not supported. Rather, several morphologically cryptic lineages comprise the diversity in this clade, including support for the recognition of T. eigenmanni. The high number of cryptic species in Typhlichthys highlights the utility of multilocus genetic data in delimiting species, particularly in lineages that exhibit slight morphological disparity, such as subterranean organisms. However, results depend on sampling of individuals and loci; this issue needs further study.