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

First vertebrate BRICHOS antimicrobial peptides: β-hairpin host defense peptides in limbless amphibia lung resemble those of marine worms

Innate immunity of invertebrates offers potent antimicrobial peptides (AMPs) against drug-resistant infections. To identify new worm β-hairpin AMPs, we explored the sequence diversity of proteins with a BRICHOS domain, which comprises worm AMP precursors. Strikingly, we discovered new BRICHOS AMPs not in worms, but in caecilians, the least studied clade of vertebrates. Two precursor proteins from Microcaecilia unicolor and Rhinatrema bivittatum resemble SP-C lung surfactants and bear worm AMP-like peptides at C-termini. The analysis of M. unicolor tissue transcriptomes shows that the AMP precursor is highly expressed in the lung along with regular SP-C, suggesting a different, protective function. The peptides form right-twisted β-hairpins, change conformation upon lipid binding, and rapidly disrupt bacterial membranes. Both peptides exhibit broad-spectrum activity against multidrug-resistant ESKAPE pathogens with 1-4 μM MICs and remarkably low toxicity, giving 40-70-fold selectivity towards bacteria. These BRICHOS AMPs, previously unseen in vertebrates, reveal a novel lung innate immunity mechanism and offer a promising antibiotics template.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

Transcriptome mining extends the host range of the Flaviviridae to non-bilaterians

The flavivirids (family Flaviviridae) are a group of positive-sense RNA viruses that include well-documented agents of human disease. Despite their importance and ubiquity, the timescale of flavivirid evolution is uncertain. An ancient origin, spanning millions of years, is supported by their presence in both vertebrates and invertebrates and by the identification of a flavivirus-derived endogenous viral element in the peach blossom jellyfish genome (Craspedacusta sowerbii, phylum Cnidaria), implying that the flaviviruses arose early in the evolution of the Metazoa. To date, however, no exogenous flavivirid sequences have been identified in these hosts. To help resolve the antiquity of the Flaviviridae, we mined publicly available transcriptome data across the Metazoa. From this, we expanded the diversity within the family through the identification of 32 novel viral sequences and extended the host range of the pestiviruses to include amphibians, reptiles, and ray-finned fish. Through co-phylogenetic analysis we found cross-species transmission to be the predominate macroevolutionary event across the non-vectored flavivirid genera (median, 68 per cent), including a cross-species transmission event between bats and rodents, although long-term virus-host co-divergence was still a regular occurrence (median, 23 per cent). Notably, we discovered flavivirus-like sequences in basal metazoan species, including the first associated with Cnidaria. This sequence formed a basal lineage to the genus Flavivirus and was closer to arthropod and crustacean flaviviruses than those in the tamanavirus group, which includes a variety of invertebrate and vertebrate viruses. Combined, these data attest to an ancient origin of the flaviviruses, likely close to the emergence of the metazoans 750-800 million years ago.

Evolutionary diversification of epidermal barrier genes in amphibians

The epidermal differentiation complex (EDC) is a cluster of genes encoding components of the skin barrier in terrestrial vertebrates. EDC genes can be categorized as S100 fused-type protein (SFTP) genes such as filaggrin, which contain two coding exons, and single-coding-exon EDC (SEDC) genes such as loricrin. SFTPs are known to be present in amniotes (mammals, reptiles and birds) and amphibians, whereas SEDCs have not yet been reported in amphibians. Here, we show that caecilians (Amphibia: Gymnophiona) have both SFTP and SEDC genes. Two to four SEDC genes were identified in the genomes of Rhinatrema bivittatum, Microcaecilia unicolor and Geotrypetes seraphini. Comparative analysis of tissue transcriptomes indicated predominant expression of SEDC genes in the skin of caecilians. The proteins encoded by caecilian SEDC genes resemble human SEDC proteins, such as involucrin and small proline-rich proteins, with regard to low sequence complexity and high contents of proline, glutamine and lysine. Our data reveal diversification of EDC genes in amphibians and suggest that SEDC-type skin barrier genes have originated either in a common ancestor of tetrapods followed by loss in Batrachia (frogs and salamanders) or, by convergent evolution, in caecilians and amniotes.

Novel Expansion of Matrix Metalloproteases in the Laboratory Axolotl (Ambystoma mexicanum) and Other Salamander Species

Matrix metalloprotease (MMP) genes encode endopeptidases that cleave protein components of the extracellular matrix (ECM) as well as non-ECM proteins. Here we report the results of a comprehensive survey of MMPs in the laboratory axolotl and other representative salamanders. Surprisingly, 28 MMPs were identified in salamanders and 9 MMP paralogs were identified as unique to the axolotl and other salamander taxa, with several of these presenting atypical amino acid insertions not observed in other tetrapod vertebrates. Furthermore, as assessed by sequence information, all of the novel salamander MMPs are of the secreted type, rather than cell membrane anchored. This suggests that secreted type MMPs expanded uniquely within salamanders to presumably execute catalytic activities in the extracellular milieu. To facilitate future studies of salamander-specific MMPs, we annotated transcriptional information from published studies of limb and tail regeneration. Our analysis sets the stage for comparative studies to understand why MMPs expanded uniquely within salamanders.

The gastrin-releasing peptide/bombesin system revisited by a reverse-evolutionary study considering Xenopus

Bombesin is a putative antibacterial peptide isolated from the skin of the frog, Bombina bombina. Two related (bombesin-like) peptides, gastrin-releasing peptide (GRP) and neuromedin B (NMB) have been found in mammals. The history of GRP/bombesin discovery has caused little attention to be paid to the evolutionary relationship of GRP/bombesin and their receptors in vertebrates. We have classified the peptides and their receptors from the phylogenetic viewpoint using a newly established genetic database and bioinformatics. Here we show, by using a clawed frog (Xenopus tropicalis), that GRP is not a mammalian counterpart of bombesin and also that, whereas the GRP system is widely conserved among vertebrates, the NMB/bombesin system has diversified in certain lineages, in particular in frog species. To understand the derivation of GRP system in the ancestor of mammals, we have focused on the GRP system in Xenopus. Gene expression analyses combined with immunohistochemistry and Western blotting experiments demonstrated that GRP peptides and their receptors are distributed in the brain and stomach of Xenopus. We conclude that GRP peptides and their receptors have evolved from ancestral (GRP-like peptide) homologues to play multiple roles in both the gut and the brain as one of the ‘gut-brain peptide’ systems.

Gigantic Genomes Provide Empirical Tests of Transposable Element Dynamics Models

Transposable elements (TEs) are a major determinant of eukaryotic genome size. The collective properties of a genomic TE community reveal the history of TE/host evolutionary dynamics and impact present-day host structure and function, from genome to organism levels. In rare cases, TE community/genome size has greatly expanded in animals, associated with increased cell size and changes to anatomy and physiology. Here, we characterize the TE landscape of the genome and transcriptome in an amphibian with a giant genome — the caecilianIchthyophis bannanicus, which we show has a genome size of 12.2 Gb. Amphibians are an important model system because the clade includes independent cases of genomic gigantism. The I. bannanicus genome differs compositionally from other giant amphibian genomes, but shares a low rate of ectopic recombination-mediated deletion. We examine TE activity using expression and divergence plots; TEs account for 15% of somatic transcription, and most superfamilies appear active. We quantify TE diversity in the caecilian, as well as other vertebrates with a range of genome sizes, using diversity indices commonly applied in community ecology. We synthesize previous models that integrate TE abundance, diversity, and activity, and test whether the caecilian meets model predictions for genomes with high TE abundance. We propose thorough, consistent characterization of TEs to strengthen future comparative analyses. Such analyses will ultimately be required to reveal whether the divergent TE assemblages found across convergent gigantic genomes reflect fundamental shared features of TE/host genome evolutionary dynamics.

Conserved keratin gene clusters in ancient fish: An evolutionary seed for terrestrial adaptation

Type I and type II keratins are subgroups of intermediate filament proteins that provide toughness to the epidermis and protect it from water loss. In terrestrial vertebrates, the keratin genes form two major clusters, clusters 1 and 2, each of which is dominated by type I and II keratin genes. By contrast, such clusters are not observed in teleost fish. Although the diversification of keratins is believed to have made a substantial contribution to terrestrial adaptation, its evolutionary process has not been clarified. Here, we performed a comprehensive genomic survey of the keratin genes of a broad range of vertebrates. As a result, we found that ancient fish lineages such as elephant shark, reedfish, spotted gar, and coelacanth share both keratin gene clusters. We also discovered an expansion of keratin genes that form a novel subcluster in reedfish. Syntenic and phylogenetic analyses revealed that two pairs of krt18/krt8 keratin genes were shared among all vertebrates, thus implying that they encode ancestral type I and II keratin protein sets. We further revealed that distinct keratin gene subclusters, which show specific expressions in the epidermis of adult amphibians, stemmed from canonical keratin genes in non-terrestrial ancestors. Molecular evolutionary analyses suggested that the selective constraints were relaxed in the adult epidermal subclusters of amphibians as well as the novel subcluster of reedfish. The results of the present study represent the process of diversification of keratins through a series of gene duplications that could have facilitated the terrestrial adaptation of vertebrates.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Insights into the skin of caecilian amphibians from gene expression profiles

Background

Gene expression profiles can provide insights into the molecular machinery behind tissue functions and, in turn, can further our understanding of environmental responses, and developmental and evolutionary processes. During vertebrate evolution, the skin has played a crucial role, displaying a wide diversity of essential functions. To unravel the molecular basis of skin specialisations and adaptations, we compared gene expression in the skin with eight other tissues in a phylogenetically and ecologically diverse species sample of one of the most neglected vertebrate groups, the caecilian amphibians (order Gymnophiona).

Results

The skin of the five studied caecilian species showed a distinct gene expression profile reflecting its developmental origin and showing similarities to other epithelial tissues. We identified 59 sequences with conserved enhanced expression in the skin that might be associated with caecilian dermal specialisations. Some of the up-regulated genes shared expression patterns with human skin and potentially are involved in skin functions across vertebrates. Variation trends in gene expression were detected between mid and posterior body skin suggesting different functions between body regions. Several candidate biologically active peptides were also annotated.

Conclusions

Our study provides the first atlas of differentially expressed sequences in caecilian tissues and a baseline to explore the molecular basis of the skin functions in caecilian amphibians, and more broadly in vertebrates.

Morphological Evidence for an Oral Venom System in Caecilian Amphibians

Amphibians are known for their skin rich in glands containing toxins employed in passive chemical defense against predators, different from, for example, snakes that have active chemical defense, injecting their venom into the prey. Caecilians (Amphibia, Gymnophiona) are snake-shaped animals with fossorial habits, considered one of the least known vertebrate groups. We show here that amphibian caecilians, including species from the basal groups, besides having cutaneous poisonous glands as other amphibians do, possess specific glands at the base of the teeth that produce enzymes commonly found in venoms. Our analysis of the origin of these glands shows that they originate from the same tissue that gives rise to teeth, similar to the venom glands in reptiles. We speculate that caecilians might have independently developed mechanisms of production and injection of toxins early in their evolutionary history.

•

Amphibian caecilians have tooth-related glands in both upper and lower jaws

•

The glands have the same origin of reptile venom glands

•

The secretion contains proteins with enzymatic activities commonly found in venoms

•

Caecilians might have developed the ability to inject oral toxins early in evolution

POSTMORTEM FINDINGS IN EIGHT SPECIES OF CAPTIVE CAECILIAN (AMPHIBIA: GYMNOPHIONA) OVER A TEN-YEAR PERIOD

Between July 2007 and June 2017 there were 86 deaths in the populations of eight caecilian species at the Zoological Society of London (ZSL) London Zoo. The mortality rate (deaths per animal-year at risk) ranged from 0.03 in the Congo caecilian (Herpele squalostoma) to 0.85 in Kaup's caecilian (Potomotyphlus kaupii). Among the 73 individuals examined post mortem, no cause of death or primary diagnosis could be established in 35 cases, but of the others the most common cause of death was dermatitis (22 cases). When all significant pathological findings were considered, skin lesions of varying types were again the commonest (56 cases), particularly among the aquatic species: Typhlonectes compressicauda (18 out of 21 cases), T. natans (8/10) and P. kaupii (12/14). Other common findings were poor gut-fill (35 cases), kidney and gastrointestinal lesions (10 cases each), generalized congestion (8 cases) and poor body condition (6 cases). This review adds to the growing body of knowledge regarding the presentations and causes of disease in captive caecilians.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

What lies beneath? Molecular evolution during the radiation of caecilian amphibians

Background

Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians.

Results

A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata.

Conclusions

The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5694-1) contains supplementary material, which is available to authorized users.