Geography is more important than life history in the recent diversification of the tiger salamander complex.

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research.

The influence of locus number and information content on species delimitation: an empirical test case in an endangered Mexican salamander.

Perhaps the most important recent advance in species delimitation has been the development of model-based approaches to objectively diagnose species diversity from genetic data. Additionally, the growing accessibility of next-generation sequence data sets provides powerful insights into genome-wide patterns of divergence during speciation. However, applying complex models to large data sets is time-consuming and computationally costly, requiring careful consideration of the influence of both individual and population sampling, as well as the number and informativeness of loci on species delimitation conclusions. Here, we investigated how locus number and information content affect species delimitation results for an endangered Mexican salamander species, Ambystoma ordinarium. We compared results for an eight-locus, 137-individual data set and an 89-locus, seven-individual data set. For both data sets, we used species discovery methods to define delimitation models and species validation methods to rigorously test these hypotheses. We also used integrated demographic model selection tools to choose among delimitation models, while accounting for gene flow. Our results indicate that while cryptic lineages may be delimited with relatively few loci, sampling larger numbers of loci may be required to ensure that enough informative loci are available to accurately identify and validate shallow-scale divergences. These analyses highlight the importance of striking a balance between dense sampling of loci and individuals, particularly in shallowly diverged lineages. They also suggest the presence of a currently unrecognized, endangered species in the western part of A. ordinarium's range.

Parallel tagged amplicon sequencing reveals major lineages and phylogenetic structure in the North American tiger salamander (Ambystoma tigrinum) species complex.

Modern analytical methods for population genetics and phylogenetics are expected to provide more accurate results when data from multiple genome-wide loci are analysed. We present the results of an initial application of parallel tagged sequencing (PTS) on a next-generation platform to sequence thousands of barcoded PCR amplicons generated from 95 nuclear loci and 93 individuals sampled across the range of the tiger salamander (Ambystoma tigrinum) species complex. To manage the bioinformatic processing of this large data set (344 330 reads), we developed a pipeline that sorts PTS data by barcode and locus, identifies high-quality variable nucleotides and yields phased haplotype sequences for each individual at each locus. Our sequencing and bioinformatic strategy resulted in a genome-wide data set with relatively low levels of missing data and a wide range of nucleotide variation. structure analyses of these data in a genotypic format resulted in strongly supported assignments for the majority of individuals into nine geographically defined genetic clusters. Species tree analyses of the most variable loci using a multi-species coalescent model resulted in strong support for most branches in the species tree; however, analyses including more than 50 loci produced parameter sampling trends that indicated a lack of convergence on the posterior distribution. Overall, these results demonstrate the potential for amplicon-based PTS to rapidly generate large-scale data for population genetic and phylogenetic-based research.

A support vector machine based test for incongruence between sets of trees in tree space.

BACKGROUND: The increased use of multi-locus data sets for phylogenetic reconstruction has increased the need to determine whether a set of gene trees significantly deviate from the phylogenetic patterns of other genes. Such unusual gene trees may have been influenced by other evolutionary processes such as selection, gene duplication, or horizontal gene transfer. RESULTS: Motivated by this problem we propose a nonparametric goodness-of-fit test for two empirical distributions of gene trees, and we developed the software GeneOut to estimate a p-value for the test. Our approach maps trees into a multi-dimensional vector space and then applies support vector machines (SVMs) to measure the separation between two sets of pre-defined trees. We use a permutation test to assess the significance of the SVM separation. To demonstrate the performance of GeneOut, we applied it to the comparison of gene trees simulated within different species trees across a range of species tree depths. Applied directly to sets of simulated gene trees with large sample sizes, GeneOut was able to detect very small differences between two set of gene trees generated under different species trees. Our statistical test can also include tree reconstruction into its test framework through a variety of phylogenetic optimality criteria. When applied to DNA sequence data simulated from different sets of gene trees, results in the form of receiver operating characteristic (ROC) curves indicated that GeneOut performed well in the detection of differences between sets of trees with different distributions in a multi-dimensional space. Furthermore, it controlled false positive and false negative rates very well, indicating a high degree of accuracy. CONCLUSIONS: The non-parametric nature of our statistical test provides fast and efficient analyses, and makes it an applicable test for any scenario where evolutionary or other factors can lead to trees with different multi-dimensional distributions. The software GeneOut is freely available under the GNU public license.

A population genetic signature of human releases in an invasive ladybeetle.

Biological invasions have been accelerated by a variety of human activities. Propagule pressure, the number of introduced individuals and independent introductions, is probably to be influenced by these human activities and may be an important factor for successful range expansion in new environments. We tested whether the current distribution of the predatory ladybeetle Coccinella septempunctata in the introduced range (USA) is the result of multiple historical human introductions or natural range expansion from the first established populations in the USA. To test this hypothesis, we compared historical records of propagule size, propagule number, specific introduction locations and the date of each introduction, with estimates of genetic variation in mitochondrial DNA (cytochrome oxidase I). Our results indicated that genetic diversity in the introduced range was positively correlated with historical records of propagule size and number and negatively correlated with distance to nearest introduction point, suggesting that multiple human releases were successful. Higher genetic diversity in populations found near introduction points suggest that initial founder effects were limited, but lower genetic diversity found farther from introduction points is probably the result of serial founder effects during secondary range expansion. These results suggest that the current distribution of C. septempunctata in the introduced range is the result of a combination of human releases and short-range expansion from multiple established populations in the introduced range.

Cast adrift on an island: introduced populations experience an altered balance between selection and drift.

A long-standing question in evolutionary biology is what becomes of adaptive traits when a species expands its range into novel environments. Here, we report the results of a study on an adaptive colour pattern polymorphism (stripes) of the coqui frog, Eleutherodactylus coqui, following its introduction to Hawaii from Puerto Rico. We compared population differentiation (Φ ST and F ST ) for the stripes locus--which underlies this colour pattern polymorphism--with neutral microsatellite loci to test for a signature of selection among native and introduced populations. Among native populations, Φ ST and F ST for stripes were lower than expected under the neutral model, suggesting uniform balancing selection. Alternatively, among introduced populations, Φ ST and F ST for stripes did not differ from the neutral model. These results suggest that the evolutionary dynamics of this previously adaptive trait have become dominated by random genetic drift following the range expansion.

Genetic basis of a color pattern polymorphism in the Coqui Frog Eleutherodactylus coqui.

Many species of frog exhibit striking color and pattern polymorphisms, but the genetic bases of these traits are not known for most species. The coqui frog, Eleutherodactylus coqui, a species endemic to the island of Puerto Rico, exhibits a wide variety of color and pattern polymorphisms including 4 discrete stripe patterns on its dorsal surface and an unstriped morph. We conducted breeding experiments to determine the mode of inheritance for these 5 dorsal color patterns in E. coqui. We analyzed results from 14 different cross types, which included 1519 offspring from 71 clutches. We found that color patterns segregate at ratios consistent with a single autosomal locus, 5-allele model, in which all alleles coding for stripes are codominant and the allele coding for the unstriped morph is recessive. We propose that this locus be named "stripes" with alleles B (interocular bar), L (dorsolateral stripes), N (narrow middorsal stripe), W (wide middorsal stripe), and u (unstriped). The results of this experiment suggest the genetic basis of stripe patterns in this well-studied species and provide a model for studying the evolution and maintenance of this phenotypic polymorphism.

Strong founder effects and low genetic diversity in introduced populations of Coqui frogs.

The success of non-native species may depend on the genetic resources maintained through the invasion process. The Coqui (Eleutherodactylus coqui), a frog endemic to Puerto Rico, was introduced to Hawaii in the late 1980s via the horticulture trade, and has become an aggressive invader. To explore whether genetic diversity and population structure changed with the introduction, we assessed individuals from 15 populations across the Hawaiian Islands and 13 populations across Puerto Rico using six to nine polymorphic microsatellite loci and five dorsolateral colour patterns. Allelic richness (R(T)) and gene diversity were significantly higher in Puerto Rico than in Hawaii populations. Hawaii also had fewer colour patterns (two versus three to five per population) than Puerto Rico. We found no isolation by distance in the introduced range, even though it exists in the native range. Results suggest extensive mixing among frog populations across Hawaii, and that their spread has been facilitated by humans. Like previous research, our results suggest that Hawaiian Coquis were founded by individuals from sites around San Juan, but unlike previous research the colour pattern and molecular genetic data (nuclear and mtDNA) support two separate introductions, one on the island of Hawaii and one on Maui. Coquis are successful invaders in Hawaii despite the loss of genetic variation. Future introductions may increase genetic variation and potentially its range.

Isolation of microsatellite loci from the coqui frog, Eleutherodactylus coqui.

Thirteen microsatellite loci were isolated from the coqui frog (Eleutherodactylus coqui) and optimized for future research. The loci were screened across 37 individuals from two Puerto Rican populations. Loci were variable with the number of alleles per locus ranging from three to 38. Polymorphic information content ranged from 0.453 to 0.963 and observed heterozygosity for each population ranged from 0.320 to 0.920.