A mitogenomic perspective on the phylogeny and biogeography of living caecilians (Amphibia: Gymnophiona)

Resistance Is Not Futile: Widespread Convergent Evolution of Resistance to Alpha-Neurotoxic Snake Venoms in Caecilians (Amphibia: Gymnophiona)

Predatory innovations impose reciprocal selection pressures upon prey. The evolution of snake venom alpha-neurotoxins has triggered the corresponding evolution of resistance in the post-synaptic nicotinic acetylcholine receptors of prey in a complex chemical arms race. All other things being equal, animals like caecilians (an Order of legless amphibians) are quite vulnerable to predation by fossorial elapid snakes and their powerful alpha-neurotoxic venoms; thus, they are under strong selective pressure. Here, we sequenced the nicotinic acetylcholine receptor alpha-1 subunit of 37 caecilian species, representing all currently known families of caecilians from across the Americas, Africa, and Asia, including species endemic to the Seychelles. Three types of resistance were identified: (1) steric hindrance from N-glycosylated asparagines; (2) secondary structural changes due to the replacement of proline by another amino acid; and (3) electrostatic charge repulsion of the positively charged neurotoxins, through the introduction of a positively charged amino acid into the toxin-binding site. We demonstrated that resistance to alpha-neurotoxins convergently evolved at least fifteen times across the caecilian tree (three times in Africa, seven times in the Americas, and five times in Asia). Additionally, as several species were shown to possess multiple resistance modifications acting synergistically, caecilians must have undergone at least 20 separate events involving the origin of toxin resistance. On the other hand, resistance in non-caecilian amphibians was found to be limited to five origins. Together, the mutations underlying resistance in caecilians constitute a robust signature of positive selection which strongly correlates with elapid presence through both space (sympatry with caecilian-eating elapids) and time (Cenozoic radiation of elapids). Our study demonstrates the extent of convergent evolution that can be expected when a single widespread predatory adaptation triggers parallel evolutionary arms races at a global scale.

Illustrated keys and a DNA barcode reference library of the amphibians and terrestrial reptiles (Amphibia, Reptilia) of São Tomé and Príncipe (Gulf of Guinea, West Africa)

The herpetofauna of São Tomé and Príncipe consists of nine species of amphibians, all endemic, and 21 species of terrestrial reptiles, of which 17 are endemic. Our current knowledge regarding its natural history, ecology, and distribution is limited. Here two important tools are provided to support researchers, conservationists, and local authorities in the identification of the country's herpetofauna: an illustrated key to the herpetofauna of the two islands and surroundings islets and a DNA barcode reference library. The keys allow a rapid and unambiguous morphological identification of all occurring species. The DNA barcodes for the entire herpetofauna of the country were produced from 79 specimens, all of which are deposited in museum collections. The barcodes generated are available in online repositories and can be used to provide unambiguous molecular identification of most of the species. Future applications and use of these tools are briefly discussed.

An aryl hydrocarbon receptor from the caecilian Gymnopis multiplicata suggests low dioxin affinity in the ancestor of all three amphibian orders

The aryl hydrocarbon receptor (AHR) plays pleiotropic roles in the development and physiology of vertebrates in conjunction with xenobiotic and endogenous ligands. It is best known for mediating the toxic effects of dioxin-like pollutants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While most vertebrates possess at least one AHR that binds TCDD tightly, amphibian AHRs bind TCDD with very low affinity. Previous analyses of AHRs from Xenopus laevis (a frog; order Anura) and Ambystoma mexicanum (a salamander; order Caudata) identified three amino acid residues in the ligand-binding domain (LBD) that underlie low-affinity binding. In X. laevis AHR1β, these are A354, A370, and N325. Here we extend the analysis of amphibian AHRs to the caecilian Gymnopis multiplicata, representing the remaining extant amphibian order, Gymnophiona. G. multiplicata AHR groups with the monophyletic vertebrate AHR/AHR1 clade. The LBD includes all three signature residues of low TCDD affinity, and a structural homology model suggests that its architecture closely resembles those of other amphibians. In transactivation assays, the EC50 for reporter gene induction by TCDD was 17.17 nM, comparable to X. laevis AhR1β (26.23 nM) and Ambystoma AHR (34.09 nM) and dramatically higher than mouse AhR (0.13 nM), a trend generally reflected in direct measures of TCDD binding. These shared properties distinguish amphibian AHRs from the high-affinity proteins typical of both vertebrate groups that diverged earlier (teleost fish) and those that appeared more recently (other tetrapods). These findings suggest the hypothesis that AHRs with low TCDD affinity represent a characteristic that evolved in a common ancestor of all three extant amphibian groups.

A review of the fossil record of caecilians (Lissamphibia: Gymnophionomorpha) with comments on its use to calibrate molecular timetrees

Gymnophiona, popularly known as caecilians, the most poorly known major taxon of extant amphibians, are elongate and limbless tetrapods, with compact ossified skulls and reduced eyes, mainly adapted to fossorial life as adults. Caecilians are poorly represented in the fossil record, but despite the scarcity of fossil specimens described (only four named taxa, in addition to indeterminate fragmentary material), their fossils play a key role in our knowledge of the origin and evolution of Lissamphibia, as well as contribute directly to a better understanding of the phylogeny, taxonomy and biogeography of extant gymnophionan taxa. These records are scattered throughout geological time (from the Jurassic to the sub-Recent) and space (North and South America and Africa). Here, we revisit the caecilian fossil record, providing a brief description of all known extinct taxa described so far, along with general remarks about their impact on systematics, time range, and geographical distribution of the clade, as well as prospects for future research. Possible calibration constraints based on the caecilian fossil record are provided.

Morphology and histology of the male reproductive tract of Caecilia thompsoni (amphibia: Gymnophiona)

We studied the male reproductive tract of individuals of different body sizes of Caecilia thompsoni to describe morphological characteristics in comparison to other Gymnophiona. The reproductive tract consists of paired testes segmented into chains of primary and secondary lobes, sperm ducts that empty to Wolffian ducts, the cloaca that receives the Wolffian ducts and possesses a phallodeum. Müllerian ducts are present and develop into paired glands that empty into the cloacal urodeum. Testicular secondary lobes contain lobules with cysts of the entire germinal cell line, whereas primary lobes, in the terminal ends of the chains, only have spermagonia, Sertoli cells, and connective tissue. The smallest individual examined (21 cm body length) was immature and only possessed a few testicular primary lobes. Once the individuals reach sexual maturity, the morphological characteristics are quite consistent at macroscopic and histological level among males of very different body sizes. The histological features of the Wolffian and Müllerian glands suggest a complementary secretory role between the two ducts. In the cloaca we found the propulsor muscle, venous sinuses, and blind sacs in the phallodeum, which differentiate C. thompsoni from other species of the genus. Despite these slight differences, the general morphological characteristics, both macroscopic and microscopic, of the reproductive tracts of adult males of C. thompsoni follow the pattern known for the reproductively active males of Gymnophiona.

A new species of Caecilia (Gymnophiona, Caeciliidae) from the Magdalena valley region of Colombia

A new species of the genus Caecilia (Caeciliidae) from the western foothills of the Serranía de los Yariguíes in Colombia is described. Caeciliapulchraserranasp. nov. is similar to C.degenerata and C.corpulenta but differs from these species in having fewer primary annular grooves and a shorter body length. With this new species, the currently recognized species in the genus are increased to 35. Mitochondrial DNA sequences, including newly sequenced terminals representing two additional, previously unanalyzed species, corroborate the phylogenetic position of the new species within Caecilia and the monophyly of the genus. This analysis also included newly sequenced terminals of Epicrionopsaff.parkeri (Rhinatrematidae) and trans-Andean Microcaecilianicefori (Siphonopidae). Evidence was found for the non-monophyly of the family Siphonopidae and the siphonopid genera Microcaecilia and Siphonops. The implications of these results for caecilian systematics are discussed and the status of the trans-Andean populations of Caeciliadegenerata is commented upon.

Stem caecilian from the Triassic of Colorado sheds light on the origins of Lissamphibia

The origin of the limbless caecilians remains a lasting question in vertebrate evolution. Molecular phylogenies and morphology support that caecilians are the sister taxon of batrachians (frogs and salamanders), from which they diverged no later than the early Permian. Although recent efforts have discovered new, early members of the batrachian lineage, the record of pre-Cretaceous caecilians is limited to a single species, Eocaecilia micropodia The position of Eocaecilia within tetrapod phylogeny is controversial, as it already acquired the specialized morphology that characterizes modern caecilians by the Jurassic. Here, we report on a small amphibian from the Upper Triassic of Colorado, United States, with a mélange of caecilian synapomorphies and general lissamphibian plesiomorphies. We evaluated its relationships by designing an inclusive phylogenetic analysis that broadly incorporates definitive members of the modern lissamphibian orders and a diversity of extinct temnospondyl amphibians, including stereospondyls. Our results place the taxon confidently within lissamphibians but demonstrate that the diversity of Permian and Triassic stereospondyls also falls within this group. This hypothesis of caecilian origins closes a substantial morphologic and temporal gap and explains the appeal of morphology-based polyphyly hypotheses for the origins of Lissamphibia while reconciling molecular support for the group's monophyly. Stem caecilian morphology reveals a previously unrecognized stepwise acquisition of typical caecilian cranial apomorphies during the Triassic. A major implication is that many Paleozoic total group lissamphibians (i.e., higher temnospondyls, including the stereospondyl subclade) fall within crown Lissamphibia, which must have originated before 315 million years ago.

Next-Generation Mitogenomics: A Comparison of Approaches Applied to Caecilian Amphibian Phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a ‘traditional’ Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina’s HiSeq and MiSeq, Roche’s 454 GS FLX, and Life Technologies’ Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Population genetics of the São Tomé caecilian (Gymnophiona: Dermophiidae: Schistometopum thomense) reveals strong geographic structuring

Islands provide exciting opportunities for exploring ecological and evolutionary mechanisms. The oceanic island of São Tomé in the Gulf of Guinea exhibits high diversity of fauna including the endemic caecilian amphibian, Schistometopum thomense. Variation in pigmentation, morphology and size of this taxon over its c. 45 km island range is extreme, motivating a number of taxonomic, ecological, and evolutionary hypotheses to explain the observed diversity. We conducted a population genetic study of S. thomense using partial sequences of two mitochondrial DNA genes (ND4 and 16S), together with morphological examination, to address competing hypotheses of taxonomic or clinal variation. Using Bayesian phylogenetic analysis and Spatial Analysis of Molecular Variance, we found evidence of four geographic clades, whose range and approximated age (c. 253 Kya – 27 Kya) are consistent with the spread and age of recent volcanic flows. These clades explained 90% of variation in ND4 (φ_CT = 0.892), and diverged by 4.3% minimum pairwise distance at the deepest node. Most notably, using Mismatch Distributions and Mantel Tests, we identified a zone of population admixture that dissected the island. In the northern clade, we found evidence of recent population expansion (Fu's F_s = −13.08 and Tajima's D = −1.80) and limited dispersal (Mantel correlation coefficient = 0.36, p = 0.01). Color assignment to clades was not absolute. Paired with multinomial regression of chromatic data, our analyses suggested that the genetic groups and a latitudinal gradient together describe variation in color of S. thomense. We propose that volcanism and limited dispersal ability are the likely proximal causes of the observed genetic structure. This is the first population genetic study of any caecilian and demonstrates that these animals have deep genetic divisions over very small areas in accordance with previous speculations of low dispersal abilities.

Influence of fossoriality on inner ear morphology: insights from caecilian amphibians

It is widely accepted that a relationship exists between inner ear morphology and functional aspects of an animal's biology, such as locomotor behaviour. Animals that engage in agile and spatially complex behaviours possess semicircular canals that morphologically maximise sensitivity to correspondingly complex physical stimuli. Stemming from the prediction that fossorial tetrapods require a well-developed sense of spatial awareness, we investigate the hypothesis that fossoriality leads to inner ear morphology that is convergent with other spatially adept tetrapods. We apply morphometrics to otic capsule endocasts of 26 caecilian species to quantify aspects of inner ear shape, and compare these with a sample of frog and salamander species. Our results reveal caecilians (and also frogs) possess strongly curved canals, a feature in common with spatially adept species. However, significantly shorter canals in caecilians suggest reduced sensitivity, possibly associated with reduced reliance on vestibulo-ocular reflexes in this group of visually degenerate tetrapods. An elaboration of the sacculus of caecilians is interpreted as a unique adaptation among amphibians to increase sensitivity to substrate-borne vibrations transmitted through the head. This study represents the first quantitative analyses of inner ear morphology of limbless fossorial tetrapods, and identifies features within a new behavioural context that will contribute to our understanding of the biological consequences of physical stimuli on sensory function and associated morphological evolution.

Life-history evolution and mitogenomic phylogeny of caecilian amphibians

We analyze mitochondrial genomes to reconstruct a robust phylogenetic framework for caecilian amphibians and use this to investigate life-history evolution within the group. Our study comprises 45 caecilian mitochondrial genomes (19 of them newly reported), representing all families and 27 of 32 currently recognized genera, including some for which molecular data had never been reported. Support for all relationships in the inferred phylogenetic tree is high to maximal, and topology tests reject all investigated alternatives, indicating an exceptionally robust molecular phylogenetic framework of caecilian evolution consistent with current morphology-based supraspecific classification. We used the mitogenomic phylogenetic framework to infer ancestral character states and to assess correlation among three life-history traits (free-living larvae, viviparity, specialized pre-adult or vernal teeth), each of which occurs only in some caecilian species. Our results provide evidence that an ancestor of the Seychelles caecilians abandoned direct development and re-evolved a free-living larval stage. This study yields insights into the concurrent evolution of direct development and of vernal teeth in an ancestor of Teresomata that likely gave rise to skin-feeding (maternal dermatophagy) behavior and subsequently enabled evolution of viviparity, with skin feeding possibly a homologous precursor of oviduct feeding in viviparous caecilians.

Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs

Anura (frogs and toads) constitute over 88% of living amphibian diversity but many important questions about their phylogeny and evolution remain unresolved. For this study, we developed an efficient method for sequencing anuran mitochondrial DNAs (mtDNAs) by amplifying the mitochondrial genome in 12 overlapping fragments using frog-specific universal primer sets. Based on this method, we generated 47 nearly complete, new anuran mitochondrial genomes and discovered nine novel gene arrangements. By combining the new data and published anuran mitochondrial genomes, we assembled a large mitogenomic data set (11,007 nt) including 90 frog species, representing 39 of 53 recognized anuran families, to investigate their phylogenetic relationships and evolutionary history. The resulting tree strongly supported a paraphyletic arrangement of archaeobatrachian (=nonneobatrachian) frogs, with Leiopelmatoidea branching first, followed by Discoglossoidea, Pipoidea, and Pelobatoidea. Within Neobatrachia, the South African Heleophrynidae is the sister-taxon to all other neobatrachian frogs and the Seychelles-endemic Sooglossidae is recovered as the sister-taxon to Ranoidea. These phylogenetic relationships agree with many nuclear gene studies. The chronogram derived from two Bayesian relaxed clock methods (MultiDivTime and BEAST) suggests that modern frogs (Anura) originated in the early Triassic about 244 Ma and the appearance of Neobatrachia took place in the late Jurassic about 163 Ma. The initial diversifications of two species-rich superfamilies Hyloidea and Ranoidea commenced 110 and 133 Ma, respectively. These times are older than some other estimates by approximately 30-40 My. Compared with nuclear data, mtDNA produces compatible time estimates for deep nodes (>150 Ma), but apparently older estimates for more shallow nodes. Our study shows that, although it evolves relatively rapidly and behaves much as a single locus, mtDNA performs well for both phylogenetic and divergence time inferences and will provide important reference hypotheses for the phylogeny and evolution of frogs.

Chromosomal homology of Uraeotyphlus oxyurus group of species (Amphibia, Gymnophiona, Ichthyophiidae)

Uraeotyphlus oxyurus (Dumeril et Bibron, 1841), Uraeotyphlus interruptus Pillai et Ravichandran, 1999, Uraeotyphlus narayani Seshachar, 1939 and Uraeotyphlus menoni Annandale, 1913 were cytogenetically analysed following conventional and differential staining techniques. These species show similar karyotypes with 2n=36 (FN=58). There were no traces of species-specific features in regard to C-banding and NOR staining. The comparative study of karyotypes shows chromosomal homologies among the four species. Chromosomal data seem to support the concept that two species groups exist in the genus Uraeotyphlus.

Is solid always best? Cranial performance in solid and fenestrated caecilian skulls

Caecilians (Lissamphibia: Gymnophiona) are characterized by a fossorial lifestyle that appears to play a role in the many anatomical specializations in the group. The skull, in particular, has been the focus of previous studies because it is driven into the substrate for burrowing. There are two different types of skulls in caecilians: (1) stegokrotaphic, where the squamosal completely covers the temporal region and the jaw closing muscles, and (2) zygokrotaphic, with incomplete coverage of the temporal region by the squamosal. We used 3-D imaging and modeling techniques to explore the functional consequences of these skull types in an evolutionary context. We digitally converted stegokrotaphic skulls into zygokrotaphic skulls and vice versa. We also generated a third, akinetic skull type that was presumably present in extinct caecilian ancestors. We explored the benefits and costs of the different skull types under frontal loading at different head angles with finite element analysis (FEA). Surprisingly, the differences in stress distributions and bending between the three tested skull types were minimal and not significant. This suggests that the open temporal region in zygokrotaphic skulls does not lead to poorer performance during burrowing. However, the results of the FEA suggest a strong relationship between the head angle and skull performance, implying there is an optimal head angle during burrowing.

Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia

Were molecular data available for extinct taxa, questions regarding the origins of many groups could be settled in short order. As this is not the case, various strategies have been proposed to combine paleontological and neontological data sets. The use of fossil dates as node age calibrations for divergence time estimation from molecular phylogenies is commonplace. In addition, simulations suggest that the addition of morphological data from extinct taxa may improve phylogenetic estimation when combined with molecular data for extant species, and some studies have merged morphological and molecular data to estimate combined evidence phylogenies containing both extinct and extant taxa. However, few, if any, studies have attempted to estimate divergence times using phylogenies containing both fossil and living taxa sampled for both molecular and morphological data. Here, I infer both the phylogeny and the time of origin for Lissamphibia and a number of stem tetrapods using Bayesian methods based on a data set containing morphological data for extinct taxa, molecular data for extant taxa, and molecular and morphological data for a subset of extant taxa. The results suggest that Lissamphibia is monophyletic, nested within Lepospondyli, and originated in the late Carboniferous at the earliest. This research illustrates potential pitfalls for the use of fossils as post hoc age constraints on internal nodes and highlights the importance of explicit phylogenetic analysis of extinct taxa. These results suggest that the application of fossils as minima or maxima on molecular phylogenies should be supplemented or supplanted by combined evidence analyses whenever possible.

Deciphering morphological variation in the braincase of caecilian amphibians (Gymnophiona)

High levels of morphological homoplasy have hindered progress in understanding morphological evolution within gymnophione lissamphibians. Stemming from the hypothesis that the braincase has the potential to yield phylogenetic information, the braincases of 27 species (23 genera) of gymnophione amphibians were examined using high-resolution micro-computed tomography and histologically prepared specimens. Morphology of the brain and its relationship to features of the braincase is described, and it is shown that eight different patterns exist in the distribution of foramina in the antotic region. The distribution of variants is congruent with molecule-based phylogeny. Additionally, all variants are shown to correspond directly to stages along developmental continua, suggesting that the evolutionary truncation of development in the antotic region at various stages has driven the evolution of morphology in this region. Attempts to correlate the observed morphology with proxies of putative heterochronic events (including those attributable to burrowing, life history, and size) fail to explain the distribution of morphology if each proxy is considered separately. Thus, it is concluded that either currently unrecognized causes of heterochrony or combinations thereof have influenced morphology in different lineages independently. These data identify clades whose morphology can now be reconsidered in light of previously unrecognized heterochronic events, thereby providing a foundation for future analyses of the evolution of morphology within Gymnophiona as a whole. Most significantly, these data confirm, for the first time in a lissamphibian group, that the braincase can preserve important phylogenetic information that is otherwise obscured in regions of the skull that experience strong influences from functional constraints.