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

Diversity of PBI-DdeI satellite DNA in snakes correlates with rapid independent evolution and different functional roles

To better understand PBI-DdeI satellite DNA located in the centromeric region of python, molecular evolution analysis was conducted on 40 snake species. A ladder-like pattern of DNA bands with repetition of the 194–210 bp monomer was observed in 15 species using PCR. Molecular cloning was performed to obtain 97 AT-rich monomer sequences. Phylogenetic and network analyses showed three PBI-DdeI subfamilies with sequences grouped in species-specific clusters, suggesting rapid evolution. Slow evolution was found in eight species with shared PBI-DdeI sequences, suggesting recent species diversification, allowing PBI-DdeI no time to diverge, with limited homogenization and fixation processes. Quantitative real-time PCR showed large differences in copy number between Python bivittatus and other snakes, consistent with repeat scanning of whole genome sequences. Copy numbers were significantly higher in female Naja kaouthia than in males, concurring with chromosomal distribution of PBI-DdeI specifically localized to female W chromosomes. PBI-DdeI might act as an evolutionary driver with several repeats to promote W chromosome differentiation and heterochromatinization in N. kaouthia. Analysis revealed PBI-DdeI with a reduced copy number, compared to P. bivittatus, in most snakes studied, and it is possible that it subsequently dispersed and amplified on W chromosomes with different functional roles in N. kaouthia.

Shifts in Selective Pressures on Snake Phototransduction Genes Associated with Photoreceptor Transmutation and Dim-Light Ancestry

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole-genome and targeted capture sequencing data. We used these data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures.

Rapid identification of common medicinal snakes and their adulterants using the Bar-HRM analysis method

Effective identification methods for snake species are lacking, exacerbating the extermination of medicinal and commercially valuable snake species. Hence, it is imperative to find fast and reliable methods to distinguish snake samples available on the market. Seventy-three samples from four families belonging to 13 genera were collected in China and found to contain common medicinal snakes and their adulterants. Cytochrome oxidase I (COI) was utilized as a DNA barcode to analyse these common snakes, and a DNA mini-barcode was employed for fast detection. Then, the DNA mini-barcode assays were coupled with a high-resolution melting (HRM) analysis (Bar-HRM) to realize the rapid discrimination of these snake species. The results showed the power of DNA barcoding with COI, which was capable of distinguishing all collected snake samples, and the combined Bar-HRM method can successfully identify the adulterants and different snake species. In particular, Bar-HRM revealed Bungarus fasciatus adulterants in B. multicinctus at concentrations as low as 1.6%. Moreover, the results of the study confirmed the effectiveness of the technique in terms of the rapid identification of snakes, which has great potential for ensuring the safety of commercially valuable snake species.

The genome of the tegu lizard Salvator merianae: combining Illumina, PacBio, and optical mapping data to generate a highly contiguous assembly

Background

Reptiles are a species-rich group with great phenotypic and life history diversity but are highly underrepresented among the vertebrate species with sequenced genomes.

Results

Here, we report a high-quality genome assembly of the tegu lizard, Salvator merianae, the first lacertoid with a sequenced genome. We combined 74X Illumina short-read, 29.8X Pacific Biosciences long-read, and optical mapping data to generate a high-quality assembly with a scaffold N50 value of 55.4 Mb. The contig N50 value of this assembly is 521 Kb, making it the most contiguous reptile assembly so far. We show that the tegu assembly has the highest completeness of coding genes and conserved non-exonic elements (CNEs) compared to other reptiles. Furthermore, the tegu assembly has the highest number of evolutionarily conserved CNE pairs, corroborating a high assembly contiguity in intergenic regions. As in other reptiles, long interspersed nuclear elements comprise the most abundant transposon class. We used transcriptomic data, homology- and de novo gene predictions to annotate 22,413 coding genes, of which 16,995 (76%) likely have human orthologs as inferred by CESAR-derived gene mappings. Finally, we generated a multiple genome alignment comprising 10 squamates and 7 other amniote species and identified conserved regions that are under evolutionary constraint. CNEs cover 38 Mb (1.8%) of the tegu genome, with 3.3 Mb in these elements being squamate specific. In contrast to placental mammal-specific CNEs, very few of these squamate-specific CNEs (<20 Kb) overlap transposons, highlighting a difference in how lineage-specific CNEs originated in these two clades.

Conclusions

The tegu lizard genome together with the multiple genome alignment and comprehensive conserved element datasets provide a valuable resource for comparative genomic studies of reptiles and other amniotes.

Comparative genomic investigation of high-elevation adaptation in ectothermic snakes

Several previous genomic studies have focused on adaptation to high elevations, but these investigations have been largely limited to endotherms. Snakes of the genus Thermophis are endemic to the Tibetan plateau and therefore present an opportunity to study high-elevation adaptations in ectotherms. Here, we report the de novo assembly of the genome of a Tibetan hot-spring snake (Thermophis baileyi) and then compare its genome to the genomes of the other two species of Thermophis, as well as to the genomes of two related species of snakes that occur at lower elevations. We identify 308 putative genes that appear to be under positive selection in Thermophis We also identified genes with shared amino acid replacements in the high-elevation hot-spring snakes compared with snakes and lizards that live at low elevations, including the genes for proteins involved in DNA damage repair (FEN1) and response to hypoxia (EPAS1). Functional assays of the FEN1 alleles reveal that the Thermophis allele is more stable under UV radiation than is the ancestral allele found in low-elevation lizards and snakes. Functional assays of EPAS1 alleles suggest that the Thermophis protein has lower transactivation activity than the low-elevation forms. Our analysis identifies some convergent genetic mechanisms in high-elevation adaptation between endotherms (based on studies of mammals) and ectotherms (based on our studies of Thermophis).

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.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

Fasting for 21days leads to changes in adipose tissue and liver physiology in juvenile checkered garter snakes (Thamnophis marcianus)

Snakes often undergo periods of prolonged fasting and, under certain conditions, can survive years without food. Despite this unique phenomenon, there are relatively few reports of the physiological adaptations to fasting in snakes. At post-prandial day 1 (fed) or 21 (fasting), brain, liver, and adipose tissues were collected from juvenile checkered garter snakes (Thamnophis marcianus). There was greater glycerol-3-phosphate dehydrogenase (G3PDH)-specific activity in the liver of fasted than fed snakes (P=0.01). The mRNA abundance of various fat metabolism-associated factors was measured in brain, liver, and adipose tissue. Lipoprotein lipase (LPL) mRNA was greater in fasted than fed snakes in the brain (P=0.04). Adipose triglyceride lipase (ATGL; P=0.006) mRNA was greater in the liver of fasted than fed snakes. In adipose tissue, expression of peroxisome proliferator-activated receptor (PPAR)γ (P=0.01), and fatty acid binding protein 4 (P=0.03) was greater in fed than fasted snakes. Analysis of adipocyte morphology revealed that cross-sectional area (P=0.095) and diameter (P=0.27) were not significantly different between fed and fasted snakes. Results suggest that mean adipocyte area can be preserved during fasting by dampening lipid biosynthesis while not changing rates of lipid hydrolysis. In the liver, however, extensive lipid remodeling may provide energy and lipoproteins to maintain lipid structural integrity during energy restriction. Because the duration of fasting was not sufficient to change adipocyte size, results suggest that the liver is important as a short-term provider of energy in the snake.

Chinese alligator genome illustrates molecular adaptations

In the September issue of Cell Research, Wan et al. analyze the Chinese Alligator genome and identify multiple instances of positively-selected changes in protein-coding genes and expansions of multi-gene families that appear to explain suites of phenotypes central to the unique habits and physiology of crocodilian reptiles. They demonstrate the effectiveness of comparative systems genomics in integrating information from comparative genomics, molecular evolution, and systems biology to understand the evolutionary dynamics of complex systems.

Inference of the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes from comparative gene mapping

Comparative genome analysis of non-avian reptiles and amphibians provides important clues about the process of genome evolution in tetrapods. However, there is still only limited information available on the genome structures of these organisms. Consequently, the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes in tetrapods remain poorly understood. We constructed chromosome maps of functional genes for the Chinese soft-shelled turtle (Pelodiscus sinensis), the Siamese crocodile (Crocodylus siamensis), and the Western clawed frog (Xenopus tropicalis) and compared them with genome and/or chromosome maps of other tetrapod species (salamander, lizard, snake, chicken, and human). This is the first report on the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes inferred from comparative genomic analysis of vertebrates, which cover all major non-avian reptilian taxa (Squamata, Crocodilia, Testudines). The eight largest macrochromosomes of the turtle and chicken were equivalent, and 11 linkage groups had also remained intact in the crocodile. Linkage groups of the chicken macrochromosomes were also highly conserved in X. tropicalis, two squamates, and the salamander, but not in human. Chicken microchromosomal linkages were conserved in the squamates, which have fewer microchromosomes than chicken, and also in Xenopus and the salamander, which both lack microchromosomes; in the latter, the chicken microchromosomal segments have been integrated into macrochromosomes. Our present findings open up the possibility that the ancestral amniotes and tetrapods had at least 10 large genetic linkage groups and many microchromosomes, which corresponded to the chicken macro- and microchromosomes, respectively. The turtle and chicken might retain the microchromosomes of the amniote protokaryotype almost intact. The decrease in number and/or disappearance of microchromosomes by repeated chromosomal fusions probably occurred independently in the amphibian, squamate, crocodilian, and mammalian lineages.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus)

Background

The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara.

Results

We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers.

Conclusions

This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct “Seq-to-SSR” approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

A multi-organ transcriptome resource for the Burmese Python (Python molurus bivittatus)

Background

Snakes provide a unique vertebrate system for studying a diversity of extreme adaptations, including those related to development, metabolism, physiology, and venom. Despite their importance as research models, genomic resources for snakes are few. Among snakes, the Burmese python is the premier model for studying extremes of metabolic fluctuation and physiological remodelling. In this species, the consumption of large infrequent meals can induce a 40-fold increase in metabolic rate and more than a doubling in size of some organs. To provide a foundation for research utilizing the python, our aim was to assemble and annotate a transcriptome reference from the heart and liver. To accomplish this aim, we used the 454-FLX sequencing platform to collect sequence data from multiple cDNA libraries.

Results

We collected nearly 1 million 454 sequence reads, and assembled these into 37,245 contigs with a combined length of 13,409,006 bp. To identify known genes, these contigs were compared to chicken and lizard gene sets, and to all Genbank sequences. A total of 13,286 of these contigs were annotated based on similarity to known genes or Genbank sequences. We used gene ontology (GO) assignments to characterize the types of genes in this transcriptome resource. The raw data, transcript contig assembly, and transcript annotations are made available online for use by the broader research community.

Conclusion

These data should facilitate future studies using pythons and snakes in general, helping to further contribute to the utilization of snakes as a model evolutionary and physiological system. This sequence collection represents a major genomic resource for the Burmese python, and the large number of transcript sequences characterized should contribute to future research in this and other snake species.

Sequencing the genome of the Burmese python (Python molurus bivittatus) as a model for studying extreme adaptations in snakes

The Consortium for Snake Genomics is in the process of sequencing the genome and creating transcriptomic resources for the Burmese python. Here, we describe how this will be done, what analyses this work will include, and provide a timeline.