Ongoing declines for the world's amphibians in the face of emerging threats.

Systematic assessments of species extinction risk at regular intervals are necessary for informing conservation action1,2. Ongoing developments in taxonomy, threatening processes and research further underscore the need for reassessment3,4. Here we report the findings of the second Global Amphibian Assessment, evaluating 8,011 species for the International Union for Conservation of Nature Red List of Threatened Species. We find that amphibians are the most threatened vertebrate class (40.7% of species are globally threatened). The updated Red List Index shows that the status of amphibians is deteriorating globally, particularly for salamanders and in the Neotropics. Disease and habitat loss drove 91% of status deteriorations between 1980 and 2004. Ongoing and projected climate change effects are now of increasing concern, driving 39% of status deteriorations since 2004, followed by habitat loss (37%). Although signs of species recoveries incentivize immediate conservation action, scaled-up investment is urgently needed to reverse the current trends.

Middle Jurassic fossils document an early stage in salamander evolution.

Salamanders are an important group of living amphibians and model organisms for understanding locomotion, development, regeneration, feeding, and toxicity in tetrapods. However, their origin and early radiation remain poorly understood, with early fossil stem-salamanders so far represented by larval or incompletely known taxa. This poor record also limits understanding of the origin of Lissamphibia (i.e., frogs, salamanders, and caecilians). We report fossils from the Middle Jurassic of Scotland representing almost the entire skeleton of the enigmatic stem-salamander Marmorerpeton. We use computed tomography to visualize high-resolution three-dimensional anatomy, describing morphologies that were poorly characterized in early salamanders, including the braincase, scapulocoracoid, and lower jaw. We use these data in the context of a phylogenetic analysis intended to resolve the relationships of early and stem-salamanders, including representation of important outgroups alongside data from high-resolution imaging of extant species. Marmorerpeton is united with Karaurus, Kokartus, and others from the Middle Jurassic-Lower Cretaceous of Asia, providing evidence for an early radiation of robustly built neotenous stem-salamanders. These taxa display morphological specializations similar to the extant cryptobranchid "giant" salamanders. Our analysis also demonstrates stem-group affinities for a larger sample of Jurassic species than previously recognized, highlighting an unappreciated diversity of stem-salamanders and cautioning against the use of single species (e.g., Karaurus) as exemplars for stem-salamander anatomy. These phylogenetic findings, combined with knowledge of the near-complete skeletal anatomy of Mamorerpeton, advance our understanding of evolutionary changes on the salamander stem-lineage and provide important data on early salamanders and the origins of Batrachia and Lissamphibia.

A Triassic stem-salamander from Kyrgyzstan and the origin of salamanders.

The origin of extant amphibians remains largely obscure, with only a few early Mesozoic stem taxa known, as opposed to a much better fossil record from the mid-Jurassic on. In recent time, anurans have been traced back to Early Triassic forms and caecilians have been traced back to the Late Jurassic Eocaecilia, both of which exemplify the stepwise acquisition of apomorphies. Yet the most ancient stem-salamanders, known from mid-Jurassic rocks, shed little light on the origin of the clade. The gap between salamanders and other lissamphibians, as well as Paleozoic tetrapods, remains considerable. Here we report a new specimen of Triassurus sixtelae, a hitherto enigmatic tetrapod from the Middle/Late Triassic of Kyrgyzstan, which we identify as the geologically oldest stem-group salamander. This sheds light not only on the early evolution of the salamander body plan, but also on the origin of the group as a whole. The new, second specimen is derived from the same beds as the holotype, the Madygen Formation of southwestern Kyrgyzstan. It reveals a range of salamander characters in this taxon, pushing back the rock record of urodeles by at least 60 to 74 Ma (Carnian-Bathonian). In addition, this stem-salamander shares plesiomorphic characters with temnospondyls, especially branchiosaurids and amphibamiforms.

Species delimitation in endangered groundwater salamanders: Implications for aquifer management and biodiversity conservation.

Groundwater-dependent species are among the least-known components of global biodiversity, as well as some of the most vulnerable because of rapid groundwater depletion at regional and global scales. The karstic Edwards-Trinity aquifer system of west-central Texas is one of the most species-rich groundwater systems in the world, represented by dozens of endemic groundwater-obligate species with narrow, naturally fragmented distributions. Here, we examine how geomorphological and hydrogeological processes have driven population divergence and speciation in a radiation of salamanders (Eurycea) endemic to the Edwards-Trinity system using phylogenetic and population genetic analysis of genome-wide DNA sequence data. Results revealed complex patterns of isolation and reconnection driven by surface and subsurface hydrology, resulting in both adaptive and nonadaptive population divergence and speciation. Our results uncover cryptic species diversity and refine the borders of several threatened and endangered species. The US Endangered Species Act has been used to bring state regulation to unrestricted groundwater withdrawals in the Edwards (Balcones Fault Zone) Aquifer, where listed species are found. However, the Trinity and Edwards-Trinity (Plateau) aquifers harbor additional species with similarly small ranges that currently receive no protection from regulatory programs designed to prevent groundwater depletion. Based on regional climate models that predict increased air temperature, together with hydrologic models that project decreased springflow, we conclude that Edwards-Trinity salamanders and other codistributed groundwater-dependent organisms are highly vulnerable to extinction within the next century.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra.

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders.

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Comparative osteology of the hynobiid complex Liua-Protohynobius-Pseudohynobius (Amphibia, Urodela): â . Cranial anatomy of Pseudohynobius.

Hynobiidae are a clade of salamanders that diverged early within the crown radiation and that retain a considerable number of features plesiomorphic for the group. Their evolutionary history is informed by a fossil record that extends to the Middle Jurassic Bathonian time. Our understanding of the evolution within the total group of Hynobiidae has benefited considerably from recent discoveries of stem hynobiids but is constrained by inadequate anatomical knowledge of some extant forms. Pseudohynobius is a derived hynobiid clade consisting of five to seven extant species living endemic to southwestern China. Although this clade has been recognized for over 37 years, osteological details of these extant hynobiids remain elusive, which undoubtedly has contributed to taxonomic controversies over the hynobiid complex Liua-Protohynobius-Pseudohynobius. Here we provide a bone-by-bone study of the cranium in the five extant species of Pseudohynobius (Ps. flavomaculatus, Ps. guizhouensis, Ps. jinfo, Ps. kuankuoshuiensis and Ps. shuichengensis) based on x-ray computer tomography data for 18 specimens. Our results indicate that the cranium in each of these species has a combination of differences in morphology, proportions and articulation patterns in both dermal and endochondral bones. Our study establishes a range of intraspecific differences that will serve as organizing hypotheses for future studies as more extensive collections of these species become available. Morphological features in the cranium for terrestrial ecological adaptation in Hynobiidae are summarized. Based on the results, we also discuss the evolution and development of several potential synapomorphies of Hynobiidae, including features of the orbitosphenoid and articular.

Detecting Mosaic Patterns in Macroevolutionary Disparity.

Evolutionary biologists have long sought to understand the full complexity in pattern and process that shapes organismal diversity. Although phylogenetic comparative methods are often used to reconstruct complex evolutionary dynamics, they are typically limited to a single phenotypic trait. Extensions that accommodate multiple traits lack the ability to partition multidimensional data sets into a set of mosaic suites of evolutionarily linked characters. I introduce a comparative framework that identifies heterogeneity in evolutionary patterns across large data sets of continuous traits. Using a model of continuous trait evolution based on the differential accumulation of disparity across lineages in a phylogeny, the approach algorithmically partitions traits into a set of character suites that best explains the data, where each suite displays a distinct pattern in phylogenetic morphological disparity. When applied to empirical data, the approach revealed a mosaic pattern predicted by developmental biology. The evolutionary distinctiveness of individual suites can be investigated in more detail either by fitting conventional comparative models or by directly studying the phylogenetic patterns in disparity recovered during the analysis. This framework can supplement existing comparative approaches by inferring the complex, integrated patterns that shape evolution across the body plan from disparate developmental, morphometric, and environmental sources of phenotypic data.

Miniaturization, Genome Size, and Biological Size in a Diverse Clade of Salamanders.

AbstractGenome size (C-value) can affect organismal traits across levels of biological organization from tissue complexity to metabolism. Neotropical salamanders show wide variation in genome and body sizes, including several clades with miniature species. Because miniaturization imposes strong constraints on morphology and development and because genome size is strongly correlated with cell size, we hypothesize that body size has played an important role in the evolution of genome size in bolitoglossine salamanders. If this hypothesis is correct, then genome size and body size should be correlated in this group. Using Feulgen image analysis densitometry, we estimated genome sizes for 60 species of Neotropical salamanders. We also estimated the "biological size" of species by comparing genome size and physical body sizes in a phylogenetic context. We found a significant correlation between C-value and physical body size using optimal regression with an Ornstein-Uhlenbeck model and report the smallest salamander genome found to date. Our index of biological size showed that some salamanders with large physical body size have smaller biological body size than some miniature species and that several clades demonstrate patterns of increased or decreased biological size compared with their physical size. Our results suggest a causal relationship between physical body size and genome size and show the importance of considering the impact of both on the biological size of organisms. Indeed, biological size may be a more appropriate measure than physical size when considering phenotypic consequences of genome size evolution in many groups.

Vastly underestimated species richness of Amazonian salamanders (Plethodontidae: Bolitoglossa) and implications about plethodontid diversification.

We present data showing that the number of salamander species in Amazonia is vastly underestimated. We used DNA sequences of up to five genes (3 mitochondrial and 2 nuclear) of 366 specimens, 189 corresponding to 89 non-Amazonian nominal species and 177 Amazonian specimens, including types or topotypes, of eight of the nine recognized species in the region. By including representatives of all known species of Amazonian Bolitoglossa, except for one, and 73% of the currently 132 recognized species of the genus, our dataset represents the broadest sample of Bolitoglossa species, specimens, and geographic localities studied to date. We performed phylogenetic analyses using parsimony with tree-alignment and maximum likelihood (ML) with similarity alignment, with indels as binary characters. Our optimal topologies were used to delimit lineages that we assigned to nominal species and candidate new species following criteria that maximize the consilience of the current species taxonomy, monophyly, gaps in branch lengths, genetic distances, and geographic distribution. We contrasted the results of our species-delimitation protocol with those of Automated Barcode Gap Discovery (ABGD) and multi-rate Poisson Tree Processes (mPTP). Finally, we inferred the historical biogeography of South American salamanders by dating the trees and using dispersal-vicariance analysis (DIVA). Our results revealed a clade including almost all Amazonian salamanders, with a topology incompatible with just the currently recognized nine species. Following our species-delimitation criteria, we identified 44 putative species in Amazonia. Both ABGD and mPTP inferred more species than currently recognized, but their numbers (23-49) and limits vary. Our biogeographic analysis suggested a stepping-stone colonization of the Amazonian lowlands from Central America through the Chocó and the Andes, with several late dispersals from Amazonia back into the Andes. These biogeographic events are temporally concordant with an early land bridge between Central and South America (~10-15 MYA) and major landscape changes in Amazonia during the late Miocene and Pliocene, such as the drainage of the Pebas system, the establishment of the Amazon River, and the major orogeny of the northern Andes.

Phylogenomic data reveal reticulation and incongruence among mitochondrial candidate species in Dusky Salamanders (Desmognathus).

Gene flow between evolutionarily distinct lineages is increasingly recognized as a common occurrence. Such processes distort our ability to diagnose and delimit species, as well as confound attempts to estimate phylogenetic relationships. A conspicuous example is Dusky Salamanders (Desmognathus), a common model-system for ecology, evolution, and behavior. Only 22 species are described, 7 in the last 40 years. However, mitochondrial datasets indicate the presence of up to 45 "candidate species" and multiple paraphyletic taxa presenting a complex history of reticulation. Some authors have even suggested that the search for species boundaries in the group may be in vain. Here, we analyze nuclear and mitochondrial data containing 161 individuals from at least 49 distinct evolutionary lineages that we treat as candidate species. Concatenated and species-tree methods do not estimate fully resolved relationships among these taxa. Comparing topologies and applying methods for estimating phylogenetic networks, we find strong support for numerous instances of hybridization throughout the history of the group. We suggest that these processes may be more common than previously thought across the phylogeography-phylogenetics continuum, and that while the search for species boundaries in Desmognathus may not be in vain, it will be complicated by factors such as crypsis, parallelism, and gene-flow.

Wide Distribution but Low Differentiation: Phylogenetic Relationships and Phylogeography of Hynobius nigrescens (Amphibia: Caudata).

We surveyed the genetic structure of Hynobius nigrescens Stejneger, 1907, a lentic breeding salamander widespread throughout montane and lowland regions of northeastern Japan. We performed a mitochondrial DNA analysis to explore intraspecific genetic variation and infer the evolutionary population history of H. nigrescens. Complete 1141 bp sequences of the mitochondrial cytochrome b gene were studied for 134 adult and larval individuals collected from 62 localities, encompassing the known range of the species. Hynobius nigrescens proved to be monophyletic, including two major clades (Clade II from southwestern Hokuriku and Clade I from all other localities). The latter clade comprises four well-supported and geographically structured subclades, which show genetic distances smaller than those seen in the widely sympatric species Hynobius lichenatus. Results of population statistical analyses indicated that Clade II of the westernmost range of H. nigrescens seems to have maintained a constant population size, while Clade I from most of the northeastern species range shows a tendency of recent population expansion, which is evident in Subclades I-A from the northernmost range and I-B from southern Tohoku to northern Kanto and eastern Chubu. In contrast, Subclades I-C from northeastern Chubu and Sado Is. and I-D from northwestern Chubu to Hokuriku seem to have been relatively stable in population size. Hynobius nigrescens differs greatly from other salamander species from northeastern Japan in its much more recent periods of genetic differentiation and its pattern thereof, and is suggested to be a young faunal element in this region.

Long-term sky islands generate highly divergent lineages of a narrowly distributed stream salamander (Pachyhynobius shangchengensis) in mid-latitude mountains of East Asia.

BACKGROUND: Climate oscillation may have a profound effect on species distributions, gene flow patterns and population demography. In response to environmental change, those species restricted to montane habitats experienced expansions and contractions along elevation gradients, which can drive differentiation among sky islands. RESULTS: The Shangcheng stout salamander (Pachyhynobius shangchengensis) is a cool stream amphibian restricted to high-elevation areas in the Dabie Mountains, East China. In the present study, we used mtDNA genes (Cyt b and ND2) of 193 individuals and 12 nuclear microsatellite loci genotyped on 370 individuals, representing 6 populations (JTX, KHJ, MW, TTZ, BYM and KJY) across the taxon's distribution area, to investigate their genetic variation and evolutionary history of P. shangchengensis. Most populations showed unusually high levels of genetic diversity. Phylogenetic analyses revealed five monophyletic clades with divergence times ranging from 3.96 to 1.4 Mya. Accordingly, significant genetic differentiation was present between these populations. Bayesian skyline plot analyses provided that all populations underwent long-term population expansions since the last inter-glacial (0.13 Mya ~ 0.12 Mya). Msvar analyses found recent signals of population decline for two northern populations (JTX and KHJ) reflecting a strong bottleneck (approximately 15-fold decrease) during the mid-Holocene (about 6000 years ago). Ecological niche modelling has shown a discontinuity in suitable habitats for P. shangchengensis under different historical climatic conditions. CONCLUSIONS: Our results suggest that the niche conservatism of P. shangchengensis and sky island effects may have led to long-term isolation between populations. In sky island refuges, the mid-latitude Dabie Mountains have provided a long-term stable environment for P. shangchengensis, which has led to the accumulation of genetic diversity and has promoted genetic divergence.

Phylogenomic insights into the diversification of salamanders in the Isthmura bellii group across the Mexican highlands.

Mountain formation in Mexico has played an important role in the diversification of many Mexican taxa. The Trans-Mexican Volcanic Belt in particular has served as both a cradle of diversification and conduit for dispersal. We investigated the evolutionary history of the Isthmura bellii group of salamanders, a widespread amphibian across the Mexican highlands, using sequence capture of ultraconserved elements. Results suggest that the I. bellii group probably originated in southeastern Mexico in the late Miocene and later dispersed across the Trans-Mexican Volcanic Belt and into the Sierra Madre Occidental. Pre-Pleistocene uplift of the Trans-Volcanic Belt likely promoted early diversification by serving as a mesic land-bridge across central Mexico. These findings highlight the importance of the Trans-Volcanic Belt in generating Mexico's rich biodiversity.

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders.

BACKGROUND: Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication. RESULTS: We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years. CONCLUSIONS: The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species.

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census-based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark-recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent-based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark-recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome-based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.

Comparative microRNAome analysis of the testis and ovary of the Chinese giant salamander.

MicroRNAs (miRNAs) are 18-24 nucleotides non-coding RNAs that regulate gene expression by post-transcriptional suppression of mRNA. The Chinese giant salamander (CGS, Andrias davidianus), which is an endangered species, has become one of the important models of animal evolution; however, no miRNA studies on this species have been conducted. In this study, two small RNA libraries of CGS ovary and testis were constructed using deep sequencing technology. A bioinformatics pipeline was developed to distinguish miRNA sequences from other classes of small RNAs represented in the sequencing data. We found that many miRNAs and other small RNAs such as piRNA and tsRNA were abundant in CGS tissue. A total of 757 and 756 unique miRNAs were annotated as miRNA candidates in the ovary and testis respectively. We identified 145 miRNAs in CGS ovary and 155 miRNAs in CGS testis that were homologous to those in Xenopus laevis ovary and testis respectively. Forty-five miRNAs were more highly expressed in ovary than in testis and 21 miRNAs were more highly expressed in testis than in ovary. The expression profiles of the selected miRNAs (miR-451, miR-10c, miR-101, miR-202, miR-7a and miR-499) had their own different roles in other eight tissues and different development stages of testis and ovary, suggesting that these miRNAs play vital regulatory roles in sexual differentiation, gametogenesis and development in CGS. To our knowledge, this is the first study to reveal miRNA profiles that are related to male and female CGS gonads and provide insights into sex differences in miRNA expression in CGS.

Testing the Relationships between Diversification, Species Richness, and Trait Evolution.

Understanding which traits drive species diversification is essential for macroevolutionary studies and to understand patterns of species richness among clades. An important tool for testing if traits influence diversification is to estimate rates of net diversification for each clade, and then test for a relationship between traits and diversification rates among clades. However, this general approach has become very controversial. Numerous papers have now stated that it is inappropriate to analyze net diversification rates in groups in which clade richness is not positively correlated with clade age. Similarly, some have stated that variation in net diversification rates does not explain variation in species richness patterns among clades across the Tree of Life. Some authors have also suggested that strong correlations between richness and diversification rates are a statistical artifact and effectively inevitable. If this latter point is true, then correlations between richness and diversification rates would be uninformative (or even misleading) for identifying how much variation in species richness among clades is explained by variation in net diversification rates. Here, we use simulations (based on empirical data for plethodontid salamanders) to address three main questions. First, how is variation in net diversification rates among clades related to the relationship between clade age and species richness? Second, how accurate are these net diversification rate estimators, and does the age-richness relationship have any relevance to their accuracy? Third, is a relationship between species richness and diversification rates an inevitable, statistical artifact? Our simulations show that strong, positive age-richness relationships arise when diversification rates are invariant among clades, whereas realistic variation in diversification rates among clades frequently disrupts this relationship. Thus, a significant age-richness relationship should not be a requirement for utilizing net diversification rates in macroevolutionary studies. Moreover, we find no difference in the accuracy of net diversification rate estimators between conditions in which there are strong, positive relationships between clade age and richness and conditions in which these strong relationships are absent. We find that net diversification rate estimators are reasonably accurate under many conditions (true and estimated rates are strongly corrrelated, and typically differ by ∼10-20%), but become more accurate when clades are older and less accurate when using incorrect assumptions about extinction. We also find that significant relationships between richness and diversification rates fail to arise under many conditions, especially when there are faster rates in younger clades. Therefore, a significant relationship between richness and diversification rates is not inevitable. Given this latter result, we suggest that relationships between richness and diversification should be tested for when attempting to explain the causes of richness patterns, to avoid potential misinterpretations (e.g., high diversification rates associated with low-richness clades). Similarly, our results also provide some support for previous studies suggesting that variation in diversification rates might explain much of the variation in species richness among major clades, based on strong relationships between clade richness and diversification rates.

Enlarged Multilocus Data set Provides Surprisingly Younger Time of Origin for the Plethodontidae, the Largest Family of Salamanders.

Deep phylogenetic relationships of the largest salamander family Plethodontidae have been difficult to resolve, probably reflecting a rapid diversification early in their evolutionary history. Here, data from 50 independent nuclear markers (total 48,582 bp) are used to reconstruct the phylogeny and divergence times for plethodontid salamanders, using both concatenation and coalescence-based species tree analyses. Our results robustly resolve the position of the enigmatic eastern North American four-toed salamander (Hemidactylium) as the sister taxon of Batrachoseps + Tribe Bolitoglossini, thus settling a long-standing question. Furthermore, we statistically reject sister taxon status of Karsenia and Hydromantes, the only plethodontids to occur outside the Americas, leading us to new biogeographic hypotheses. Contrary to previous long-standing arguments that plethodontid salamanders are an old lineage originating in the Cretaceous (more than 90 Ma), our analyses lead to the hypothesis that these salamanders are much younger, arising close to the K-T boundary (~66 Ma). These time estimates are highly stable using alternative calibration schemes and dating methods. Our data simulation highlights the potential risk of making strong arguments about phylogenetic timing based on inferences from a handful of nuclear genes, a common practice. Based on the newly obtained timetree and ancestral area reconstruction results, we argue that (i) the classic "Out of Appalachia" hypothesis of plethodontid origins is problematic; (ii) the common ancestor of extant plethodontids may have originated in northwestern North America in the early Paleocene; (iii) origins of Eurasian plethodontids likely result from two separate dispersal events from western North America via Beringia in the late Eocene (~42 Ma) and the early Miocene (~23 Ma), respectively.

Testing models of refugial isolation, colonization and population connectivity in two species of montane salamanders.

Pleistocene glaciations have profoundly affected patterns of genetic diversity within many species. Temperate alpine organisms likely experienced dramatic range shifts, given that much of their habitat was glaciated during this time. While the effects of glaciations are relatively well understood, the spatial locations of refugia and processes that gave rise to current patterns of diversity are less well known. We use a microsatellite data set to test hypotheses of population connectivity and refugial isolation in the web-toed salamanders (Hydromantes) of the Sierra Nevada. We reject models of refugia with subsequent expansion into either the high southern Sierra or low-elevation Owens Valley, in favor of a simple isolation model with no migration between current populations. We find no evidence of migration at even moderate spatial scales using a variety of analyses in the southern Sierra, and limited migration in the northern Sierra. These results suggest that divergence in isolation following fragmentation is the dominant process structuring genetic variation in these salamander species. In the context of anthropogenic climate change and habitat degradation, these results imply that salamanders and other low-vagility alpine organisms are at risk of decline as they are unlikely to migrate across large distances.