Changes in Hox genes' structure and function during the evolution of the squamate body plan

Vertebral evolution and the diversification of squamate reptiles

Taxonomic, morphological, and functional diversity are often discordant and independent components of diversity. A fundamental and largely unanswered question in evolutionary biology is why some clades diversify primarily in some of these components and not others. Dramatic variation in trunk vertebral numbers (14 to >300) among squamate reptiles coincides with different body shapes, and snake-like body shapes have evolved numerous times. However, whether increased evolutionary rates or numbers of vertebrae underlie body shape and taxonomic diversification is unknown. Using a supertree of squamates including 1375 species, and corresponding vertebral and body shape data, we show that increased rates of evolution in vertebral numbers have coincided with increased rates and disparity in body shape evolution, but not changes in rates of taxonomic diversification. We also show that the evolution of many vertebrae has not spurred or inhibited body shape or taxonomic diversification, suggesting that increased vertebral number is not a key innovation. Our findings demonstrate that lineage attributes such as the relaxation of constraints on vertebral number can facilitate the evolution of novel body shapes, but that different factors are responsible for body shape and taxonomic diversification.

Reptilian-transcriptome v1.0, a glimpse in the brain transcriptome of five divergent Sauropsida lineages and the phylogenetic position of turtles

Background

Reptiles are largely under-represented in comparative genomics despite the fact that they are substantially more diverse in many respects than mammals. Given the high divergence of reptiles from classical model species, next-generation sequencing of their transcriptomes is an approach of choice for gene identification and annotation.

Results

Here, we use 454 technology to sequence the brain transcriptome of four divergent reptilian and one reference avian species: the Nile crocodile, the corn snake, the bearded dragon, the red-eared turtle, and the chicken. Using an in-house pipeline for recursive similarity searches of >3,000,000 reads against multiple databases from 7 reference vertebrates, we compile a reptilian comparative transcriptomics dataset, with homology assignment for 20,000 to 31,000 transcripts per species and a cumulated non-redundant sequence length of 248.6 Mbases. Our approach identifies the majority (87%) of chicken brain transcripts and about 50% of de novo assembled reptilian transcripts. In addition to 57,502 microsatellite loci, we identify thousands of SNP and indel polymorphisms for population genetic and linkage analyses. We also build very large multiple alignments for Sauropsida and mammals (two million residues per species) and perform extensive phylogenetic analyses suggesting that turtles are not basal living reptiles but are rather associated with Archosaurians, hence, potentially answering a long-standing question in the phylogeny of Amniotes.

Conclusions

The reptilian transcriptome (freely available at http://www.reptilian-transcriptomes.org) should prove a useful new resource as reptiles are becoming important new models for comparative genomics, ecology, and evolutionary developmental genetics.

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.

Complex selection associated with Hox genes in a natural population of lizards

Hox genes are recognized for their explanatory power of bilateral development. However, relatively little is known about natural variation in, and the evolutionary dynamics of, Hox genes within wild populations. Utilizing a natural population of sand lizards (Lacerta agilis), we screened HoxA13 for genetic variation and an association with incidence of offspring malformations. We found significant effects of parental genetic similarity and offspring sex, and their interaction, on risk of hatching malformed as an offspring. We also found within population genetic variation in HoxA13, and identified a significant effect of a three-way interaction among Hox genotype, parental genetic similarity, and offspring sex on the risk of hatching malformation. Since malformed offspring in this population do not survive to maturity, this study reveals complex and ongoing selection associated with Hox genes in a wild reptile population. Importantly, this demonstrates the utility of natural populations in unveiling microevolutionary processes shaping variation in highly conserved genes.

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.

Rampant horizontal transfer of SPIN transposons in squamate reptiles

Transposable elements (TEs) are highly abundant in the genome and capable of mobility, two properties that make them particularly prone to transfer horizontally between organisms. Although the impact of horizontal transfer (HT) of TEs is well recognized in prokaryotes, the frequency of this phenomenon and its contribution to genome evolution in eukaryotes remain poorly appreciated. Here, we provide evidence that a DNA transposon called SPIN has colonized the genome of 17 species of reptiles representing nearly every major lineage of squamates, including 14 families of lizards, snakes, and amphisbaenians. Slot blot analyses indicate that SPIN has amplified to high copy numbers in most of these species, ranging from 2,000-28,000 copies per haploid genome. In contrast, we could not detect the presence of SPIN in any of the turtles (seven species from seven families) and crocodiles (four species) examined. Genetic distances between SPIN sequences from species belonging to different squamate families are consistently very low (average = 0.1), considering the deep evolutionary divergence of the families investigated (most are >100 My diverged). Furthermore, these distances fall below interfamilial distances calculated for two genes known to have evolved under strong functional constraint in vertebrates (RAG1, average = 0.24 and C-mos, average = 0.27). These data, combined with phylogenetic analyses, indicate that the widespread distribution of SPIN among squamates is the result of at least 13 independent events of HTs. Molecular dating and paleobiogeographical data suggest that these transfers took place during the last 50 My on at least three different continents (North America, South America and, Africa). Together, these results triple the number of known SPIN transfer events among tetrapods, provide evidence for a previously hypothesized transoceanic movement of SPIN transposons during the Cenozoic, and further underscore the role of HT in the evolution of vertebrate genomes.

The transposable element profile of the anolis genome: How a lizard can provide insights into the evolution of vertebrate genome size and structure

The recent sequencing of the lizard genome provides a unique opportunity to examine the evolution of vertebrate genomes in a phylogenetic context. The lizard genome contains an extraordinary diversity of active transposable elements that far exceeds the diversity reported in extant mammals and birds. Retrotransposons and DNA transposons are represented by multiple active families, contributing to the very diverse repetitive landscape of the lizard. Surprisingly, ancient transposon copies are relatively rare suggesting that the transposon copy number is tightly controlled in lizard. This bias in favor of young copies results from the joint effect of purifying selection acting on novel insertions and a high rate of DNA loss. Recent analyses have revealed that the repetitive landscape of reptiles differ drastically from other extant amniotes by their diversity but also by the dynamics of amplification of their transposons. Thus, from the point of view of mobile elements, reptile genomes show more similarity to fish and amphibians than to other amniotes.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

Evolution of gene regulatory networks controlling body plan development

Evolutionary change in animal morphology results from alteration of the functional organization of the gene regulatory networks (GRNs) that control development of the body plan. A major mechanism of evolutionary change in GRN structure is alteration of cis-regulatory modules that determine regulatory gene expression. Here we consider the causes and consequences of GRN evolution. Although some GRN subcircuits are of great antiquity, other aspects are highly flexible and thus in any given genome more recent. This mosaic view of the evolution of GRN structure explains major aspects of evolutionary process, such as hierarchical phylogeny and discontinuities of paleontological change.

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.

Involvement of the conserved Hox gene Antennapedia in the development and evolution of a novel trait

BACKGROUND

Hox proteins specify segment identity during embryogenesis and have typical associated expression patterns. Changes in embryonic expression and activity of Hox genes were crucial in the evolution of animal body plans, but their role in the post-embryonic development of lineage-specific traits remains largely unexplored. Here, we focus on the insect Hox genes Ultrabithorax (Ubx) and Antennapedia (Antp), and implicate the latter in the formation and diversification of novel, butterfly-specific wing patterns.

RESULTS

First, we describe a conserved pattern of Ubx expression and a novel pattern of Antp expression in wing discs of Bicyclus anynana butterflies. The discrete, reiterated domains of Antp contrast with the typical expression of Hox genes in single continuous regions in arthropod embryos. Second, we show that this pattern is associated with the establishment of the organizing centres of eyespots. Antp upregulation is the earliest event in organizer development described to date, and in contrast to all genes implicated in eyespot formation, is exclusive to those centres. Third, our comparative analysis of gene expression across nymphalids reveals unexpected differences in organizer determination.

CONCLUSIONS

We show that the Antp's recruitment for the formation of novel traits in butterfly wing discs involved the evolution of new expression domains, and is restricted to a particular lineage. This study contributes novel insights into the evolution of Antp expression, as well as into the genetic mechanisms underlying morphological diversification. Our results also underscore how a wider representation of morphological and phylogenetic diversity is essential in evolutionary developmental biology.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.

Comparative analyses of vertebrate posterior HoxD clusters reveal atypical cluster architecture in the caecilian Typhlonectes natans

BACKGROUND

The posterior genes of the HoxD cluster play a crucial role in the patterning of the tetrapod limb. This region is under the control of a global, long-range enhancer that is present in all vertebrates. Variation in limb types, as is the case in amphibians, can probably not only be attributed to variation in Hox genes, but is likely to be the product of differences in gene regulation. With a collection of vertebrate genome sequences available today, we used a comparative genomics approach to study the posterior HoxD cluster of amphibians. A frog and a caecilian were included in the study to compare coding sequences as well as to determine the gain and loss of putative regulatory sequences.

RESULTS

We sequenced the posterior end of the HoxD cluster of a caecilian and performed comparative analyses of this region using HoxD clusters of other vertebrates. We determined the presence of conserved non-coding sequences and traced gains and losses of these footprints during vertebrate evolution, with particular focus on amphibians. We found that the caecilian HoxD cluster is almost three times larger than its mammalian counterpart. This enlargement is accompanied with the loss of one gene and the accumulation of repeats in that area. A similar phenomenon was observed in the coelacanth, where a different gene was lost and expansion of the area where the gene was lost has occurred. At least one phylogenetic footprint present in all vertebrates was lost in amphibians. This conserved region is a known regulatory element and functions as a boundary element in neural tissue to prevent expression of Hoxd genes.

CONCLUSION

The posterior part of the HoxD cluster of Typhlonectes natans is among the largest known today. The loss of Hoxd-12 and the expansion of the intergenic region may exert an influence on the limb enhancer, by having to bypass a distance seven times that of regular HoxD clusters. Whether or not there is a correlation with the loss of limbs remains to be investigated. These results, together with data on other vertebrates show that the tetrapod Hox clusters are more variable than previously thought.

High level of structural polymorphism driven by mobile elements in the Hox genomic region of the Chaetognath Spadella cephaloptera

Little is known about the relationships between genome polymorphism, mobile element dynamics, and population size among animal populations. The chaetognath species Spadella cephaloptera offers a unique perspective to examine this issue because they display a high level of genetic polymorphism at the population level. Here, we have investigated in detail the extent of nucleotide and structural polymorphism in a region harboring Hox1 and several coding genes and presumptive functional elements. Sequencing of several bacterial artificial chromosome inserts representative of this nuclear region uncovered a high level of structural heterogeneity, which is mainly caused by the polymorphic insertion of a diversity of genetic mobile elements. By anchoring this variation through individual genotyping, we demonstrated that sequence diversity could be attributed to the allelic pool of a single population, which was confirmed by detection of extensive recombination within the genomic region studied. The high average level of nucleotide heterozygosity provides clues of selection in both coding and noncoding domains. This pattern stresses how selective processes remarkably cope with intense sequence turnover due to substitutions, mobile element insertions, and recombination to preserve the integrity of functional landscape. These findings suggest that genome polymorphism could provide pivotal information for future functional annotation of genomes.

Hox genes and regional patterning of the vertebrate body plan

Several decades have passed since the discovery of Hox genes in the fruit fly Drosophila melanogaster. Their unique ability to regulate morphologies along the anteroposterior (AP) axis (Lewis, 1978) earned them well-deserved attention as important regulators of embryonic development. Phenotypes due to loss- and gain-of-function mutations in mouse Hox genes have revealed that the spatio-temporally controlled expression of these genes is critical for the correct morphogenesis of embryonic axial structures. Here, we review recent novel insight into the modalities of Hox protein function in imparting specific identity to anatomical regions of the vertebral column, and in controlling the emergence of these tissues concomitantly with providing them with axial identity. The control of these functions must have been intimately linked to the shaping of the body plan during evolution.