Snake Genomes Research Articles

Appendages and gene regulatory networks: Lessons from the limbless.

Among squamate reptiles, dozens of lineages have independently evolved complete or partial limb reduction. This remarkable convergence of limbless and limb-reduced phenotypes provides multiple natural replicates of different ages to explore the evolution and development of the vertebrate limb and the gene regulatory network that controls its formation. The most successful and best known of the limb-reduced squamates are snakes, which evolved a limb-reduced body form more than 100 million years ago. Recent studies have revealed the unexpected finding that many ancient limb enhancers are conserved in the genomes of snakes. Analyses in limbed animals show that many of these limb enhancers are also active during development of the phallus, suggesting that these enhancers may have been retained in snakes due their importance in regulating transcription in the external genitalia. This hypothesis is substantiated by functional tests of snake enhancers, which demonstrate that snake enhancer elements have lost limb function while retaining genital enhancer function. The large degree of overlap in the gene regulatory networks deployed during limb and phallus development may act to constrain the divergence of shared gene network components and the evolution of appendage morphology. Future studies will reveal whether limb regulatory elements have undergone similar functional changes in other lineages of limb-reduced squamates.

A workflow of massive identification and application of intron markers using snakes as a model.

Relative to the commonly used mitochondrial and nuclear protein‐coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single‐copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next‐generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p‐distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.

Genomic regression of claw keratin, taste receptor and light-associated genes provides insights into biology and evolutionary origins of snakes

Regressive evolution of anatomical traits often corresponds with the regression of genomic loci underlying such characters. As such, studying patterns of gene loss can be instrumental in addressing questions of gene function, resolving conflicting results from anatomical studies, and understanding the evolutionary history of clades. The evolutionary origins of snakes involved the regression of a number of anatomical traits, including limbs, taste buds and the visual system, and by analyzing serpent genomes, I was able to test three hypotheses associated with the regression of these features. The first concerns two keratins that are putatively specific to claws. Both genes that encode these keratins are pseudogenized/deleted in snake genomes, providing additional evidence of claw-specificity. The second hypothesis is that snakes lack taste buds, an issue complicated by conflicting results in the literature. I found evidence that different snakes have lost one or more taste receptors, but all snakes examined retained at least one gustatory channel. The final hypothesis addressed is that the earliest snakes were adapted to a dim light niche. I found evidence of deleted and pseudogenized genes with light-associated functions in snakes, demonstrating a pattern of gene loss similar to other dim light-adapted clades. Molecular dating estimates suggest that dim light adaptation preceded the loss of limbs, providing some bearing on interpretations of the ecological origins of snakes.

A Step-by-Step Guide to Assemble a Reptilian Genome.

Multiple technologies and software are now available facilitating the de novo sequencing and assembly of any vertebrate genome. Yet the quality of most available sequenced genomes is substantially poorer than that of the golden standard in the field: the human genome. Here, we present a step-by-step protocol for the successful sequencing and assembly of a high-quality snake genome that can be applied to any other reptilian or avian species. We combine the great sequencing depth and accuracy of short reads with the use of different insert size libraries for extended scaffolding followed by optical mapping. We show that this procedure improved the corn snake scaffold N50 from 3.7kbp to 1.4Mbp, currently making it one of the snake genomes with the longest scaffolds. Short guidelines are also given on the extraction of long DNA molecules from reptilian blood and the necessary modifications in DNA extraction protocols. This chapter is accompanied by a website ( www.reptilomics.org/stepbystep.html ), where we provide links to the suggested software, examples of input and output files, and running parameters.

Snake Genome Sequencing: Results and Future Prospects.

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.

Non‐random accumulation of LINE1‐like sequences on differentiated snake W chromosomes

AbstractDue to their particular phylogenetic position and biological characteristics, squamate reptiles and, in particular, snakes are becoming an increasingly important model for fields such as evolutionary biology, molecular ecology and adaptation. Recently, during a study to analyze the evolutionary history of European whip snakes, we found a LINE1 (L1)‐like sequence (GenBank accession no. LM644476), herein called TRL1L, and while there are data on the abundance of L1 in snakes, their genomic and chromosome localization is still largely unexplored. We therefore performed a study to obtain information on TRL1L abundance, distribution and conservation in snake species, belonging to the Colubridae, Lamprophiidae and Viperidae families, using quantitative dot‐blot and fluorescence in situ hybridization (FISH). TRL1L showed a high identity with homologous segments of L1s of lizards the Anolis carolinensis and Lacerta agilis and the zebrafish Danio rerio. The discovered sequences are truncated L1 elements which occur with a low copy number, about 0.1% of the genome of the species studied. This evidence suggests that L1 retroposons have a similar landscape in lizard and snake genomes, probably because similar processes limited L1 distribution in their genomes. TRL1L showed a non‐random chromosome distribution pattern. It was scarcely located on autosomes and on the euchromatic W chromosome of Cerastes vipera, while mostly found on the heterochromatic W chromosome of Hierophis carbonarius and Elaphe quatuorlineata. Our data support the hypothesis that a ‘purifying selection’ against the accumulation of L1 elements takes place in recombining regions and highlight the possible role of these elements in the differentiation processes of the snake heterochromatic W chromosome. Interestingly, the preferential distribution of TRL1L on the heterochromatic W chromosomes of the studied snakes appears to be similar to that observed in mammals for L1 accumulation on differentiated Y chromosomes. This finding suggests that a convergent process may have taken place in the differentiation of vertebrate heterochromatic sex chromosomes.

Shared Enhancer Activity in the Limbs and Phallus and Functional Divergence of a Limb-Genital cis-Regulatory Element in Snakes

The amniote phallus and limbs differ dramatically in their morphologies but share patterns of signaling and gene expression in early development. Thus far, the extent to which genital and limb transcriptional networks also share cis-regulatory elements has remained unexplored. We show that many limb enhancers are retained in snake genomes, suggesting that these elements may function in non-limb tissues. Consistent with this, our analysis of cis-regulatory activity in mice and Anolis lizards reveals that patterns of enhancer activity in embryonic limbs and genitalia overlap heavily. In mice, deletion of HLEB, an enhancer of Tbx4, produces defects in hindlimbs and genitalia, establishing the importance of this limb-genital enhancer for development of these different appendages. Further analyses demonstrate that the HLEB of snakes has lost hindlimb enhancer function while retaining genital activity. Our findings identify roles for Tbx4 in genital development and highlight deep similarities in cis-regulatory activity between limbs and genitalia.

Expression of venom gene homologs in diverse python tissues suggests a new model for the evolution of snake venom.

Snake venom gene evolution has been studied intensively over the past several decades, yet most previous studies have lacked the context of complete snake genomes and the full context of gene expression across diverse snake tissues. We took a novel approach to studying snake venom evolution by leveraging the complete genome of the Burmese python, including information from tissue-specific patterns of gene expression. We identified the orthologs of snake venom genes in the python genome, and conducted detailed analysis of gene expression of these venom homologs to identify patterns that differ between snake venom gene families and all other genes. We found that venom gene homologs in the python are expressed in many different tissues outside of oral glands, which illustrates the pitfalls of using transcriptomic data alone to define "venom toxins." We hypothesize that the python may represent an ancestral state prior to major venom development, which is supported by our finding that the expansion of venom gene families is largely restricted to highly venomous caenophidian snakes. Therefore, the python provides insight into biases in which genes were recruited for snake venom systems. Python venom homologs are generally expressed at lower levels, have higher variance among tissues, and are expressed in fewer organs compared with all other python genes. We propose a model for the evolution of snake venoms in which venom genes are recruited preferentially from genes with particular expression profile characteristics, which facilitate a nearly neutral transition toward specialized venom system expression.

Presence of a glycine-cysteine-rich beta-protein in the oberhautchen layer of snake epidermis marks the formation of the shedding layer.

The complex differentiation of snake epidermis largely depends on the variation in the production of glycine-cysteine-rich versus glycine-rich beta-proteins (beta-keratins) that are deposited on a framework of alpha-keratins. The knowledge of the amino acid sequences of beta-proteins in the snake Pantherophis guttatus has allowed the localization of a glycine-cysteine-rich beta-protein in the spinulated oberhautchen layer of the differentiating shedding complex before molting takes place. This protein decreases in the beta-layer and disappears in mesos and alpha-layers. Conversely, while the mRNA for a glycine-rich beta-protein is highly expressed in differentiating beta-cells, the immunolocalization for this protein is low in these cells. This discrepancy between expression and localization suggests that the epitope in glycine-rich beta-proteins is cleaved or modified by posttranslational processes that take place during the differentiation and maturation of the beta-layer. The present study suggests that among the numerous beta-proteins coded in the snake genome to produce epidermal layers with different textures, the glycine-cysteine-rich beta-protein marks the shedding complex formed between alpha- and beta-layers that allows for molting while its disappearance between the beta- and alpha-layers (mesos region for scale growth) is connected to the formation of the alpha-layers.

Reptile genomes open the frontier for comparative analysis of amniote development and regeneration.

Developmental genetic studies of vertebrates have focused primarily on zebrafish, frog and mouse models, which have clear application to medicine and well-developed genomic resources. In contrast, reptiles represent the most diverse amniote group, but have only recently begun to gather the attention of genome sequencing efforts. Extant reptilian groups last shared a common ancestor ?280 million years ago and include lepidosaurs, turtles and crocodilians. This phylogenetic diversity is reflected in great morphological and behavioral diversity capturing the attention of biologists interested in mechanisms regulating developmental processes such as somitogenesis and spinal patterning, regeneration, the evolution of "snake-like" morphology, the formation of the unique turtle shell, and the convergent evolution of the four-chambered heart shared by mammals and archosaurs. The complete genome of the first non-avian reptile, the green anole lizard, was published in 2011 and has provided insights into the origin and evolution of amniotes. Since then, the genomes of multiple snakes, turtles, and crocodilians have also been completed. Here we will review the current diversity of available reptile genomes, with an emphasis on their evolutionary relationships, and will highlight how these genomes have and will continue to facilitate research in developmental and regenerative biology.

Genes for Extremes

The first two snake genomes, published this week, reflect the amazing evolutionary tales of a prey-crushing python and a venomous cobra.

Accelerated Evolution of Fetuin Family Proteins inProtobothrops flavoviridis(Habu Snake) Serum and the Discovery of an L1-Like Genomic Element in the Intronic Sequence of a Fetuin-Encoding Gene

Habu serum factor (HSF) and HSF-like protein (HLP) are fetuin family proteins isolated from Protobothrops flavoviridis (habu snake) serum with different physiological activities. A comparison of their cDNAs and intronic sequences revealed that nucleotide substitutions were primarily in protein-coding regions, and the substitution patterns indicated accelerated evolution of these proteins. Genomic DNA fragment analysis, including intron 1, revealed a 6.6-kb insertion homologous to the full-length mammalian LINE1 (L1) retrotransposable element (PfL1) only in the HLP gene. This segment retains an open reading frame (ORF) that encodes a reverse transcriptase (RT)-like protein (PfRT). We further found that a large number of homologous segments have dispersed in the habu snake genome, although we could not determine the enzymatic activities of their products. Moreover, an analysis of habu snake liver RNA indicated active transcription of the PfRT genes, suggesting that high levels of RT activity in this snake have driven the evolution of unique phenotypes of venom enzymes and serum inhibitors of them.

Complete mitochondrial genome sequence from an endangered Indian snake, Python molurus molurus (Serpentes, Pythonidae)

This paper reports the complete mitochondrial genome sequence of an endangered Indian snake, Python molurus molurus (Indian Rock Python). A typical snake mitochondrial (mt) genome of 17258 bp length comprising of 37 genes including the 13 protein coding genes, 22 tRNA genes, and 2 ribosomal RNA genes along with duplicate control regions is described herein. The P. molurus molurus mt. genome is relatively similar to other snake mt. genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases. The nucleotide composition of the genome shows that there are more A-C % than T-G% on the positive strand as revealed by positive AT and CG skews. Comparison of individual protein coding genes, with other snake genomes suggests that ATP8 and NADH3 genes have high divergence rates. Codon usage analysis reveals a preference of NNC codons over NNG codons in the mt. genome of P. molurus. Also, the synonymous and non-synonymous substitution rates (ka/ks) suggest that most of the protein coding genes are under purifying selection pressure. The phylogenetic analyses involving the concatenated 13 protein coding genes of P. molurus molurus conformed to the previously established snake phylogeny.

A proposal to sequence the genome of a garter snake (Thamnophis sirtalis)

Here we develop an argument in support of sequencing a garter snake (Thamnophis sirtalis) genome, and outline a plan to accomplish this. This snake is a common, widespread, nonvenomous North American species that has served as a model for diverse studies in evolutionary biology, physiology, genomics, behavior and coevolution. The anole lizard is currently the only genome sequence available for a non-avian reptile. Thus, the garter snake at this time would be the first available snake genome sequence and as such would provide much needed comparative representation of non-avian reptilian genomes, and would also allow critical new insights for vertebrate comparative genomic studies. We outline the major areas of discovery that the availability of the garter snake genome would enable, and describe a plan for whole-genome sequencing.

Discovery of Highly Divergent Repeat Landscapes in Snake Genomes Using High-Throughput Sequencing

We conducted a comprehensive assessment of genomic repeat content in two snake genomes, the venomous copperhead (Agkistrodon contortrix) and the Burmese python (Python molurus bivittatus). These two genomes are both relatively small (∼1.4 Gb) but have surprisingly extensive differences in the abundance and expansion histories of their repeat elements. In the python, the readily identifiable repeat element content is low (21%), similar to bird genomes, whereas that of the copperhead is higher (45%), similar to mammalian genomes. The copperhead's greater repeat content arises from the recent expansion of many different microsatellites and transposable element (TE) families, and the copperhead had 23-fold greater levels of TE-related transcripts than the python. This suggests the possibility that greater TE activity in the copperhead is ongoing. Expansion of CR1 LINEs in the copperhead genome has resulted in TE-mediated microsatellite expansion (“microsatellite seeding”) at a scale several orders of magnitude greater than previously observed in vertebrates. Snakes also appear to be prone to horizontal transfer of TEs, particularly in the copperhead lineage. The reason that the copperhead has such a small genome in the face of so much recent expansion of repeat elements remains an open question, although selective pressure related to extreme metabolic performance is an obvious candidate. TE activity can affect gene regulation as well as rates of recombination and gene duplication, and it is therefore possible that TE activity played a role in the evolution of major adaptations in snakes; some evidence suggests this may include the evolution of venom repertoires.

Evolution of the mitochondrial genome in snakes: Gene rearrangements and phylogenetic relationships

BackgroundSnakes as a major reptile group display a variety of morphological characteristics pertaining to their diverse behaviours. Despite abundant analyses of morphological characters, molecular studies using mitochondrial and nuclear genes are limited. As a result, the phylogeny of snakes remains controversial. Previous studies on mitochondrial genomes of snakes have demonstrated duplication of the control region and translocation of trnL to be two notable features of the alethinophidian (all serpents except blindsnakes and threadsnakes) mtDNAs. Our purpose is to further investigate the gene organizations, evolution of the snake mitochondrial genome, and phylogenetic relationships among several major snake families.ResultsThe mitochondrial genomes were sequenced for four taxa representing four different families, and each had a different gene arrangement. Comparative analyses with other snake mitochondrial genomes allowed us to summarize six types of mitochondrial gene arrangement in snakes. Phylogenetic reconstruction with commonly used methods of phylogenetic inference (BI, ML, MP, NJ) arrived at a similar topology, which was used to reconstruct the evolution of mitochondrial gene arrangements in snakes.ConclusionThe phylogenetic relationships among the major families of snakes are in accordance with the mitochondrial genomes in terms of gene arrangements. The gene arrangement in Ramphotyphlops braminus mtDNA is inferred to be ancestral for snakes. After the divergence of the early Ramphotyphlops lineage, three types of rearrangements occurred. These changes involve translocations within the IQM tRNA gene cluster and the duplication of the CR. All phylogenetic methods support the placement of Enhydris plumbea outside of the (Colubridae + Elapidae) cluster, providing mitochondrial genomic evidence for the familial rank of Homalopsidae.

Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity

As more data are generated from proteome and transcriptome analyses of snake venoms, we are gaining an appreciation of the complexity of the venoms and, to some degree, the various sources of such complexity. However, our knowledge is still far from complete. The translation of genetic information from the snake genome to the transcriptome and ultimately the proteome is only beginning to be appreciated, and will require significantly more investigation of the snake venom genomic structure prior to a complete understanding of the genesis of venom composition. Venom complexity, however, is derived not only from the venom genomic structure but also from transcriptome generation and translation and, perhaps most importantly, post-translation modification of the nascent venom proteome. In this review, we examine the snake venom metalloproteinases, some of the predominant components in viperid venoms, with regard to possible synthesis and post-translational mechanisms that contribute to venom complexity. The aim of this review is to highlight the state of our knowledge on snake venom metalloproteinase post-translational processing and to suggest testable hypotheses regarding the cellular mechanisms associated with snake venom metalloproteinase complexity in venoms.

Reappraisal of Vipera aspis Venom Neurotoxicity

BackgroundThe variation of venom composition with geography is an important aspect of intraspecific variability in the Vipera genus, although causes of this variability remain unclear. The diversity of snake venom is important both for our understanding of venomous snake evolution and for the preparation of relevant antivenoms to treat envenomations. A geographic intraspecific variation in snake venom composition was recently reported for Vipera aspis aspis venom in France. Since 1992, cases of human envenomation after Vipera aspis aspis bites in south-east France involving unexpected neurological signs were regularly reported. The presence of genes encoding PLA2 neurotoxins in the Vaa snake genome led us to investigate any neurological symptom associated with snake bites in other regions of France and in neighboring countries. In parallel, we used several approaches to characterize the venom PLA2 composition of the snakes captured in the same areas.Methodology/Principal FindingsWe conducted an epidemiological survey of snake bites in various regions of France. In parallel, we carried out the analysis of the genes and the transcripts encoding venom PLA2s. We used SELDI technology to study the diversity of PLA2 in various venom samples. Neurological signs (mainly cranial nerve disturbances) were reported after snake bites in three regions of France: Languedoc-Roussillon, Midi-Pyrénées and Provence-Alpes-Côte d'Azur. Genomes of Vipera aspis snakes from south-east France were shown to contain ammodytoxin isoforms never described in the genome of Vipera aspis from other French regions. Surprisingly, transcripts encoding venom neurotoxic PLA2s were found in snakes of Massif Central region. Accordingly, SELDI analysis of PLA2 venom composition confirmed the existence of population of neurotoxic Vipera aspis snakes in the west part of the Massif Central mountains.Conclusions/SignificanceThe association of epidemiological studies to genetic, biochemical and immunochemical analyses of snake venoms allowed a good evaluation of the potential neurotoxicity of snake bites. A correlation was found between the expression of neurological symptoms in humans and the intensity of the cross-reaction of venoms with anti-ammodytoxin antibodies, which is correlated with the level of neurotoxin (vaspin and/or ammodytoxin) expression in the venom. The origin of the two recently identified neurotoxic snake populations is discussed according to venom PLA2 genome and transcriptome data.

Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region

BackgroundThe mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.ResultsWe present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.ConclusionFluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.

GENETIC DIVERSITY IN SNAKE VENOM THREE-FINGER PROTEINS AND PHOSPHOLIPASE A2 ENZYMES

Based on genetic organization and high degree of sequence identity, three-finger protein genes and phospholipase A2 genes appear to represent two multigene families in snake genomes. Unlike intron insertions or deletions, the evolution of protein coding regions of these genes is accelerated by adaptive selection. Thus, the protein coding regions of the two multigene families are more diversified than their intron regions. Alternatively, genetic elements in regulating protein expression were mostly retained during the evolution of the two multigene families. These observations indicate that the multiplicity of the two multigene families expanded on their ancestral gene by gene duplication followed by divergence, and suggest a common route for the evolution of snake toxins.