Multiple contact zones and karyotypic evolution in a neotropical frog species complex.

Previous studies of DNA sequence and karyotypic data have revealed high genetic diversity in the Physalaemus cuvieri - Physalaemus ephippifer species complex-a group of small leptodactylid frogs in South America. To date, seven major genetic lineages have been recognized in this group, with species delimitation tests supporting four to seven of them as valid species. Among these, only P. ephippifer shows heteromorphic sex chromosomes, but the implications of cytogenetic divergence for the evolution of this group are unknown. We analyzed karyotypic, mitochondrial DNA, and 3RAD genomic data to characterize a putative contact zone between P. ephippifer and P. cuvieri Lineage 1, finding evidence for admixture and karyotypic evolution. We also describe preliminary evidence for admixture between two other members of this species complex-Lineage 1 and Lineage 3 of P. cuvieri. Our study sheds new light on evolutionary relationships in the P. cuvieri - P. ephippifer species complex, suggesting an important role of karyotypic divergence in its evolutionary history and underscoring the importance of hybridization as a mechanism of sex chromosome evolution in amphibians.

A Dataset of Amphibian Species in U.S. National Parks.

National parks and other protected areas are important for preserving landscapes and biodiversity worldwide. An essential component of the mission of the United States (U.S.) National Park Service (NPS) requires understanding and maintaining accurate inventories of species on protected lands. We describe a new, national-scale synthesis of amphibian species occurrence in the NPS system. Many park units have a list of amphibian species observed within their borders compiled from various sources and available publicly through the NPSpecies platform. However, many of the observations in NPSpecies remain unverified and the lists are often outdated. We updated the amphibian dataset for each park unit by collating old and new park-level records and had them verified by regional experts. The new dataset contains occurrence records for 292 of the 424 NPS units and includes updated taxonomy, international and state conservation rankings, hyperlinks to a supporting reference for each record, specific notes, and related fields which can be used to better understand and manage amphibian biodiversity within a single park or group of parks.

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.

Using salamanders as model taxa to understand vertebrate feeding constraints during the late Devonian water-to-land transition.

The vertebrate water-to-land transition and the rise of tetrapods brought about fundamental changes for the groups undergoing these evolutionary changes (i.e. stem and early tetrapods). These groups were forced to adapt to new conditions, including the distinct physical properties of water and air, requiring fundamental changes in anatomy. Nutrition (or feeding) was one of the prime physiological processes these vertebrates had to successfully adjust to change from aquatic to terrestrial life. The basal gnathostome feeding mode involves either jaw prehension or using water flows to aid in ingestion, transportation and food orientation. Meanwhile, processing was limited primarily to simple chewing bites. However, given their comparatively massive and relatively inflexible hyobranchial system (compared to the more muscular tongue of many tetrapods), it remains fraught with speculation how stem and early tetrapods managed to feed in both media. Here, we explore ontogenetic water-to-land transitions of salamanders as functional analogues to model potential changes in the feeding behaviour of stem and early tetrapods. Our data suggest two scenarios for terrestrial feeding in stem and early tetrapods as well as the presence of complex chewing behaviours, including excursions of the jaw in more than one dimension during early developmental stages. Our results demonstrate that terrestrial feeding may have been possible before flexible tongues evolved. This article is part of the theme issue 'Food processing and nutritional assimilation in animals'.

Adaptive divergence despite low effective population size in a peripherally isolated population of the pygmy rabbit, Brachylagus idahoensis.

Local adaptation can occur when spatially separated populations are subjected to contrasting environmental conditions. Historically, understanding the genetic basis of adaptation has been difficult, but increased availability of genome-wide markers facilitates studies of local adaptation in non-model organisms of conservation concern. The pygmy rabbit (Brachylagus idahoensis) is an imperiled lagomorph that relies on sagebrush for forage and cover. This reliance has led to widespread population declines following reductions in the distribution of sagebrush, leading to geographic separation between populations. In this study, we used >20,000 single nucleotide polymorphisms, genotype-environment association methods, and demographic modeling to examine neutral genetic variation and local adaptation in the pygmy rabbit in Nevada and California. We identified 308 loci as outliers, many of which had functional annotations related to metabolism of plant secondary compounds. Likewise, patterns of spatial variation in outlier loci were correlated with landscape and climatic variables including proximity to streams, sagebrush cover, and precipitation. We found that populations in the Mono Basin of California probably diverged from other Great Basin populations during late Pleistocene climate oscillations, and that this region is adaptively differentiated from other regions in the southern Great Basin despite limited gene flow and low effective population size. Our results demonstrate that peripherally isolated populations can maintain adaptive divergence.

Re-evaluation of the Wehrle's salamander (Plethodon wehrlei Fowler and Dunn) species group (Caudata: Plethodontidae) using genomic data, with the description of a new species.

Woodland salamanders of the genus Plethodon are characterized by strong ecological and morphological conservatism. One assemblage, the Wehrle's salamander (Plethodon wehrlei Fowler Dunn) species group, is distributed from New York to Tennessee, USA, and includes several morphological variants, four of which are sufficiently distinct to have been recognized as species in the past. For many years after two of these species were placed in synonymy, only P. wehrlei and P. punctatus Highton were recognized. A recent phylogeographic study using mitochondrial DNA and nuclear DNA uncovered considerable genetic diversity within the group and conservatively resurrected one of the previously synonymized forms (P. dixi Pope Fowler). However, their analysis could not resolve all relationships among remaining populations of P. wehrlei, leaving the taxon paraphyletic. We re-evaluated the evolutionary history of this group using genomic data, recovered strong support for at least five distinct clades, and corroborated previously reported relationships. We also collected morphological data and demonstrated morphological distinctiveness for four of the five clades that we herein recognize as species. We resurrect the synonymized name P. jacksoni Newman to represent the southern clades of P. wehrlei in southwestern Virginia and North Carolina exclusive of P. dixi. In addition, we describe a yellow-spotted form of P. wehrlei endemic to the Cumberland Plateau as a new species. Although our proposed changes rectify the paraphyly of P. wehrlei, our sampling was not sufficient to resolve potential taxonomic issues remaining within the species herein recognized as P. jacksoni.

Adapterama III: Quadruple-indexed, double/triple-enzyme RADseq libraries (2RAD/3RAD).

Molecular ecologists frequently use genome reduction strategies that rely upon restriction enzyme digestion of genomic DNA to sample consistent portions of the genome from many individuals (e.g., RADseq, GBS). However, researchers often find the existing methods expensive to initiate and/or difficult to implement consistently, especially because it is difficult to multiplex sufficient numbers of samples to fill entire sequencing lanes. Here, we introduce a low-cost and highly robust approach for the construction of dual-digest RADseq libraries that build on adapters and primers designed in Adapterama I. Major features of our method include: (1) minimizing the number of processing steps; (2) focusing on a single strand of sample DNA for library construction, allowing the use of a non-phosphorylated adapter on one end; (3) ligating adapters in the presence of active restriction enzymes, thereby reducing chimeras; (4) including an optional third restriction enzyme to cut apart adapter-dimers formed by the phosphorylated adapter, thus increasing the efficiency of adapter ligation to sample DNA, which is particularly effective when only low quantity/quality DNA samples are available; (5) interchangeable adapter designs; (6) incorporating variable-length internal indexes within the adapters to increase the scope of sample indexing, facilitate pooling, and increase sequence diversity; (7) maintaining compatibility with universal dual-indexed primers and thus, Illumina sequencing reagents and libraries; and, (8) easy modification for the identification of PCR duplicates. We present eight adapter designs that work with 72 restriction enzyme combinations. We demonstrate the efficiency of our approach by comparing it with existing methods, and we validate its utility through the discovery of many variable loci in a variety of non-model organisms. Our 2RAD/3RAD method is easy to perform, has low startup costs, has increased utility with low-concentration input DNA, and produces libraries that can be highly-multiplexed and pooled with other Illumina libraries.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix).

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext).

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.

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.