Multi-tissue transcriptomes of caecilian amphibians highlight incomplete knowledge of vertebrate gene families

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians.

Phototactic tails: Evolution and molecular basis of a novel sensory trait in sea snakes

Dermal phototaxis has been reported in a few aquatic vertebrate lineages spanning fish, amphibians and reptiles. These taxa respond to light on the skin of their elongate hind-bodies and tails by withdrawing under cover to avoid detection by predators. Here, we investigated tail phototaxis in sea snakes (Hydrophiinae), the only reptiles reported to exhibit this sensory behaviour. We conducted behavioural tests in 17 wild-caught sea snakes of eight species by illuminating the dorsal surface of the tail and midbody skin using cold white, violet, blue, green and red light. Our results confirmed phototactic tail withdrawal in the previously studied Aipysurus laevis, revealed this trait for the first time in A. duboisii and A. tenuis, and suggested that tail photoreceptors have peak spectral sensitivities between blue and green light (457-514 nm). Based on these results, and an absence of photoresponses in five Aipysurus and Hydrophis species, we tentatively infer that tail phototaxis evolved in the ancestor of a clade of six Aipysurus species (comprising 10% of all sea snakes). Quantifying tail damage, we found that the probability of sustaining tail injuries was not influenced by tail phototactic ability in snakes. Gene profiling showed that transcriptomes of both tail skin and body skin lacked visual opsins but contained melanopsin (opn4x) in addition to key genes of the retinal regeneration and phototransduction cascades. This work suggests that a nonvisual photoreceptor (e.g., Gq rhabdomeric) signalling pathway underlies tail phototaxis, and provides candidate gene targets for future studies of this unusual sensory innovation in reptiles.

Morphological evolution and modularity of the caecilian skull

BACKGROUND

Caecilians (Gymnophiona) are the least speciose extant lissamphibian order, yet living forms capture approximately 250 million years of evolution since their earliest divergences. This long history is reflected in the broad range of skull morphologies exhibited by this largely fossorial, but developmentally diverse, clade. However, this diversity of form makes quantification of caecilian cranial morphology challenging, with highly variable presence or absence of many structures. Consequently, few studies have examined morphological evolution across caecilians. This extensive variation also raises the question of degree of conservation of cranial modules (semi-autonomous subsets of highly-integrated traits) within this clade, allowing us to assess the importance of modular organisation in shaping morphological evolution. We used an intensive surface geometric morphometric approach to quantify cranial morphological variation across all 32 extant caecilian genera. We defined 16 cranial regions using 53 landmarks and 687 curve and 729 surface sliding semilandmarks. With these unprecedented high-dimensional data, we analysed cranial shape and modularity across caecilians assessing phylogenetic, allometric and ecological influences on cranial evolution, as well as investigating the relationships among integration, evolutionary rate, and morphological disparity.

RESULTS

We found highest support for a ten-module model, with greater integration of the posterior skull. Phylogenetic signal was significant (K_mult = 0.87, p < 0.01), but stronger in anterior modules, while allometric influences were also significant (R² = 0.16, p < 0.01), but stronger posteriorly. Reproductive strategy and degree of fossoriality were small but significant influences on cranial morphology (R² = 0.03-0.05), after phylogenetic (p < 0.03) and multiple-test (p < 0.05) corrections. The quadrate-squamosal 'cheek' module was the fastest evolving module, perhaps due to its pivotal role in the unique dual jaw-closing mechanism of caecilians. Highly integrated modules exhibited both high and low disparities, and no relationship was evident between integration and evolutionary rate.

CONCLUSIONS

Our high-dimensional approach robustly characterises caecilian cranial evolution and demonstrates that caecilian crania are highly modular and that cranial modules are shaped by differential phylogenetic, allometric, and ecological effects. More broadly, and in contrast to recent studies, this work suggests that there is no simple relationship between integration and evolutionary rate or disparity.

Macroevolutionary shift in the size of amphibian genomes and the role of life history and climate

The evolution and great diversity of genome size has been of long-standing interest to biologists, but has seldom been investigated on a broad phylogenetic scale. Here we present a comparative quantitative analysis of factors shaping genome size evolution in amphibians, the extant class of vertebrates with the largest variation in genome size. We find that amphibian genomes have undergone saltations in size, although these are rare and the evolutionary history of genome size in amphibians has otherwise been one of gradual, time-dependent variation (that is, Brownian motion). This macroevolutionary homogeneity is remarkable given the evolutionary and ecological diversity of most other aspects of the natural history of amphibians. Contrary to previous claims, we find no evidence for associations between life cycle complexity and genome size despite the high diversity of reproductive modes and the multiple events of independent evolution of divergent life cycles in the group. Climate (temperature and humidity) affects genome size indirectly, at least in frogs, as a consequence of its effect on premetamorphic developmental period, although directionality of the relationship between developmental period and genome size is not unequivocal.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Recent Asian origin of chytrid fungi causing global amphibian declines

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

A new species of Uropeltis Cuvier, 1829 (Serpentes: Uropeltidae) from the Anaikatty Hills of the Western Ghats of India

A new species of Uropeltis is described from a series of six type specimens from the Anaikatty Hills of the Western Ghats of Tamil Nadu, peninsular India. Uropeltis bhupathyi sp. nov. is distinguished from congeners by having more than 200 ventral scales, 17 dorsal scale rows at midbody and by the size and shape of the rostral and frontal shields. Although tens of specimens have been seen in the vicinity of the type locality (and previously reported as U. ellioti), the new species is known only from this locality and faces threats from road traffic, habitat loss and change, and possibly a condition that deforms heads and head shields which is at least superficially similar to snake fungal disease reported from wild snakes in North America and Europe.

Terrestrial reproduction as an adaptation to steep terrain in African toads

How evolutionary novelties evolve is a major question in evolutionary biology. It is widely accepted that changes in environmental conditions shift the position of selective optima, and advancements in phylogenetic comparative approaches allow the rigorous testing of such correlated transitions. A longstanding question in vertebrate biology has been the evolution of terrestrial life histories in amphibians and here, by investigating African bufonids, we test whether terrestrial modes of reproduction have evolved as adaptations to particular abiotic habitat parameters. We reconstruct and date the most complete species-level molecular phylogeny and estimate ancestral states for reproductive modes. By correlating continuous habitat measurements from remote sensing data and locality records with life-history transitions, we discover that terrestrial modes of reproduction, including viviparity evolved multiple times in this group, most often directly from fully aquatic modes. Terrestrial modes of reproduction are strongly correlated with steep terrain and low availability of accumulated water sources. Evolutionary transitions to terrestrial modes of reproduction occurred synchronously with or after transitions in habitat, and we, therefore, interpret terrestrial breeding as an adaptation to these abiotic conditions, rather than an exaptation that facilitated the colonization of montane habitats.

A new species of Rhinophis Hemprich, 1820 (Serpentes: Uropeltidae) from the central hills of Sri Lanka

A new species of uropeltid snake, Rhinophis roshanpererai sp. nov., is described based on three specimens from Badulla District, Sri Lanka. The new species is distinguished from its congeners by the number of dorsal scale rows and ventral scales, a colour pattern that lacks pale stripes, and by having a very small tail shield with four or three prominent spines. Based on lack of records from similar habitats and elevations elsewhere in Sri Lanka, the new species likely has a very small range in the central highlands. The description of the new species brings the total nominal species of Rhinophis to 20, with four of the 16 Sri Lankan species having been described since 2009.

Recharacterization of Rhinophis dorsimaculatus Deraniyagala, 1941 (Serpentes: Uropeltidae), including description of new material

The Sri Lankan uropeltid (shieldtail) snake Rhinophis dorsimaculatus Deraniyagala, 1941 was described originally from two specimens that were subsequently lost. The small amount of previously published data and lack of published colour photographs made this one of South Asia's most poorly known snake species, and this resulted in at least one instance of taxonomic misidentification. An additional 10 specimens from a historical collection from the vicinity of the type locality recently came to light. This material is reviewed and documented and the species recharacterized. An additional locality for the species is reported. The newly reported material helps to corroborate the taxonomic validity and distinctiveness of Rhinophis dorsimaculatus. The species is readily distinguished from congeners by having 227 or more ventral scales; a large, dorsally carinate rostral shield; posterior margins of paired anals that are largely separated by the posteriormost ventral scale; and a distinctive colour pattern with bilaterally asymmetrical dark blotches within a broad, pale middorsal stripe and regularly punctate flanks.

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.

A new large species of Bitis Gray, 1842 (Serpentes: Viperidae) from the Bale Mountains of Ethiopia

A new species of viperine viperid snake is described, Bitis harenna sp. nov. The new species is a member of the subgenus Macrocerastes based on it having three scales separating the nasal and rostral shields, and on the combination of 'divisions' of dorsal scale rows on the upper flanks and 'fusions' of rows on the lower flanks. Bitis harenna sp. nov. is distinguished from other members of the subgenus by its unique colour pattern, posterior parietal flange on the lateral wall of the braincase, and possibly by differences in scalation and head proportions. Only a single museum specimen is known, a female collected from 'Dodola' in Ethiopia probably in the late 1960s and previously identified as a possibly unusually coloured and patterned B. parviocula. A live, presumably male, specimen very closely resembling the holotype of Bitis harenna sp. nov. was photographed on the Harenna escarpment of the Bale Mountains National Park, Ethiopia in 2013, providing secure occurrence data and evidence that the holotype is not a uniquely aberrant specimen. A revised key to the species of Bitis in Ethiopia is presented. Aspects of body scalation are compared among species of the subgenus Macrocerastes and between species of Macrocerastes and Bitis, and several systematic characters are highlighted and clarified.

A reassessment of Melanophidium Günther, 1864 (Squamata: Serpentes: <br />Uropeltidae) from the Western Ghats of peninsular India, with the description of a new species

A new species of the uropeltid snake genus, Melanophidium Günther, 1864 is described based on a series of eight specimens. Melanophidium khairei sp. nov. is the fourth species described in the genus, and the first for 144 years. Superficially M. khairei sp. nov. resembles M. punctatum Beddome, 1871 in being piebald and punctate (and it was previously misidentified as M. punctatum), but in many scalation characters it more closely resembles M. wynaudense (Beddome, 1863). The new species occurs in southern Maharashtra, Goa, and northern Karnataka, in the Western Ghats region of peninsula India. It is the most northerly member of its genus. Lectotypes and paralectotypes are designated for M. wynaudense, M. bilineatum Beddome, 1870, and M. punctatum. A new key to the species of Melanophidium is presented. Aspects of the morphology, taxonomy and distribution of the three previously described species of Melanophidium are reviewed and revised.

Multiple rod-cone and cone-rod photoreceptor transmutations in snakes: evidence from visual opsin gene expression

In 1934, Gordon Walls forwarded his radical theory of retinal photoreceptor 'transmutation'. This proposed that rods and cones used for scotopic and photopic vision, respectively, were not fixed but could evolve into each other via a series of morphologically distinguishable intermediates. Walls' prime evidence came from series of diurnal and nocturnal geckos and snakes that appeared to have pure-cone or pure-rod retinas (in forms that Walls believed evolved from ancestors with the reverse complement) or which possessed intermediate photoreceptor cells. Walls was limited in testing his theory because the precise identity of visual pigments present in photoreceptors was then unknown. Subsequent molecular research has hitherto neglected this topic but presents new opportunities. We identify three visual opsin genes, rh1, sws1 and lws, in retinal mRNA of an ecologically and taxonomically diverse sample of snakes central to Walls' theory. We conclude that photoreceptors with superficially rod- or cone-like morphology are not limited to containing scotopic or photopic opsins, respectively. Walls' theory is essentially correct, and more research is needed to identify the patterns, processes and functional implications of transmutation. Future research will help to clarify the fundamental properties and physiology of photoreceptors adapted to function in different light levels.

Visual system evolution and the nature of the ancestral snake

The dominant hypothesis for the evolutionary origin of snakes from 'lizards' (non-snake squamates) is that stem snakes acquired many snake features while passing through a profound burrowing (fossorial) phase. To investigate this, we examined the visual pigments and their encoding opsin genes in a range of squamate reptiles, focusing on fossorial lizards and snakes. We sequenced opsin transcripts isolated from retinal cDNA and used microspectrophotometry to measure directly the spectral absorbance of the photoreceptor visual pigments in a subset of samples. In snakes, but not lizards, dedicated fossoriality (as in Scolecophidia and the alethinophidian Anilius scytale) corresponds with loss of all visual opsins other than RH1 (λmax 490-497 nm); all other snakes (including less dedicated burrowers) also have functional sws1 and lws opsin genes. In contrast, the retinas of all lizards sampled, even highly fossorial amphisbaenians with reduced eyes, express functional lws, sws1, sws2 and rh1 genes, and most also express rh2 (i.e. they express all five of the visual opsin genes present in the ancestral vertebrate). Our evidence of visual pigment complements suggests that the visual system of stem snakes was partly reduced, with two (RH2 and SWS2) of the ancestral vertebrate visual pigments being eliminated, but that this did not extend to the extreme additional loss of SWS1 and LWS that subsequently occurred (probably independently) in highly fossorial extant scolecophidians and A. scytale. We therefore consider it unlikely that the ancestral snake was as fossorial as extant scolecophidians, whether or not the latter are para- or monophyletic.

A new species of Indian caecilian highlights challenges for species delimitation within Gegeneophis Peters, 1879 (Amphibia: Gymnophiona: Indotyphlidae)

A new species of indotyphlid caecilian amphibian, Gegeneophis tejaswini sp. nov., is described based on eight specimens from lowlands of the most northerly district of the state of Kerala in the southern part of the Western Ghats region, India. This species is distinguished from all other Gegeneophis in annulation characters and genetics (&gt; 6% different from most similar nominal species for 883 base pairs of mitochondrial 12S and 16S rRNA gene sequence data). The high degree of morphological similarity of G. krishni, G. mhadeiensis and the new species underlines that, for some Gegeneophis, larger samples and/or new characters will be needed to further advance the taxonomy of this genus.