Atypical relaxation of structural constraints in Hox gene clusters of the green anole lizard

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

First Genome of Rock Lizard Darevskia valentini Involved in Formation of Several Parthenogenetic Species

The extant reptiles are one of the most diverse clades among terrestrial vertebrates and one of a few groups with instances of parthenogenesis. Due to the hybrid origin of parthenogenetic species, reference genomes of the parental species as well as of the parthenogenetic progeny are indispensable to explore the genetic foundations of parthenogenetic reproduction. Here, we report on the first genome assembly of rock lizard Darevskia valentini, a paternal species for several parthenogenetic lineages. The novel genome was used in the reconstruction of the comprehensive phylogeny of Squamata inferred independently from 7369 trees of single-copy orthologs and a supermatrix of 378 conserved proteins. We also investigated Hox clusters, the loci that are often regarded as playing an important role in the speciation of animal groups with drastically diverse morphology. We demonstrated that Hox clusters of D. valentini are invaded with transposons and contain the HoxC1 gene that has been considered to be lost in the amniote ancestor. This study provides confirmation for previous works and releases new genomic data that will contribute to future discoveries on the mechanisms of parthenogenesis as well as support comparative studies among reptiles.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

Structure encoding in DNA

It is proposed that transposons and related long non-coding RNA define the fine structure of body parts. Although morphogens have long been known to direct the formation of many gross structures in early embryonic development, they do not have the necessary precision to define a structure down to the individual cellular level. Using the distinction between procedural and declarative knowledge in information processing as an analogy, it is hypothesized that DNA encodes fine structure in a manner that is different from the genetic code for proteins. The hypothesis states that repeated or near-repeated sequences that are in transposons and non-coding RNA define body part structures. As the cells in a body part go through the epigenetic process of differentiation, the action of methylation serves to inactivate all but the relevant structure definitions and some associated cell type genes. The transposons left active will then physically modify the DNA sequence in the heterochromatin to establish the local context in the three-dimensional body part structure. This brings the encoded definition of the cell type to the histone. The histone code for that cell type starts the regulatory cascade that turns on the genes associated with that particular type of cell, transforming it from a multipotent cell to a fully differentiated cell. This mechanism creates structures in the musculoskeletal system, the organs of the body, the major parts of the brain, and other systems.

Comparative analysis of squamate brains unveils multi-level variation in cerebellar architecture associated with locomotor specialization

Ecomorphological studies evaluating the impact of environmental and biological factors on the brain have so far focused on morphology or size measurements, and the ecological relevance of potential multi-level variations in brain architecture remains unclear in vertebrates. Here, we exploit the extraordinary ecomorphological diversity of squamates to assess brain phenotypic diversification with respect to locomotor specialization, by integrating single-cell distribution and transcriptomic data along with geometric morphometric, phylogenetic, and volumetric analysis of high-definition 3D models. We reveal significant changes in cerebellar shape and size as well as alternative spatial layouts of cortical neurons and dynamic gene expression that all correlate with locomotor behaviours. These findings show that locomotor mode is a strong predictor of cerebellar structure and pattern, suggesting that major behavioural transitions in squamates are evolutionarily correlated with mosaic brain changes. Furthermore, our study amplifies the concept of ‘cerebrotype’, initially proposed for vertebrate brain proportions, towards additional shape characters.

The cerebellum is critical in sensory-motor control and is structurally diverse across vertebrates. Here, the authors investigate the evolutionary relationship between locomotory mode and cerebellum architecture across squamates by integrating study of gene expression, cell distribution, and 3D morphology.

Evolutionary lability in Hox cluster structure and gene expression in Anolis lizards

Hox genes orchestrate development by patterning the embryonic axis. Vertebrate Hox genes are arranged in four compact clusters, and the spacing between genes is assumed to be crucial for their function. The genomes of squamate reptiles are unusually rich and variable in transposable elements (TEs), and it has been suggested that TE invasion is responsible for the Hox cluster expansion seen in snakes and lizards. Using de novo TE prediction on 17 genomes of lizards and snakes, I show that TE content of Hox clusters are generally 50% lower than genome‐wide TE levels. However, two distantly related lizards of the species‐rich genus Anolis have Hox clusters with a TE content that approaches genomic levels. The age distribution of TEs in Anolis lizards revealed that peaks of TE activity broadly coincide with speciation events. In accordance with theoretical models of Hox cluster regulation, I find that Anolis species with many TEs in their Hox clusters show aberrant Hox gene expression patterns, suggesting a functional link between TE accumulation and embryonic development. These results are consistent with the hypothesis that TEs play a role in developmental processes as well as in evolutionary diversifications.

Lizards possess the most complete tetrapod Hox gene repertoire despite pervasive structural changes in Hox clusters

Hox genes are a remarkable example of conservation in animal development and their nested expression along the head-to-tail axis orchestrates embryonic patterning. Early in vertebrate history, two duplications led to the emergence of four Hox clusters (A-D) and redundancy within paralog groups has been partially accommodated with gene losses. Here we conduct an inventory of squamate Hox genes using the genomes of 10 lizard and 7 snake species. Although the HoxC1 gene has been hypothesized to be lost in the amniote ancestor, we reveal that it is retained in lizards. In contrast, all snakes lack functional HoxC1 and -D12 genes. Varying levels of degradation suggest differences in the process of gene loss between the two genes. The vertebrate HoxC1 gene is prone to gene loss and its functional domains are more variable than those of other Hox1 genes. We describe for the first time the HoxC1 expression patterns in tetrapods. HoxC1 is broadly expressed during development in the diencephalon, the neural tube, dorsal root ganglia, and limb buds in two lizard species. Our study emphasizes the value of revisiting Hox gene repertoires by densely sampling taxonomic groups and its feasibility owing to growing sequence resources in evaluating gene repertoires across taxa.

The ecological origins of snakes as revealed by skull evolution

The ecological origin of snakes remains amongst the most controversial topics in evolution, with three competing hypotheses: fossorial; marine; or terrestrial. Here we use a geometric morphometric approach integrating ecological, phylogenetic, paleontological, and developmental data for building models of skull shape and size evolution and developmental rate changes in squamates. Our large-scale data reveal that whereas the most recent common ancestor of crown snakes had a small skull with a shape undeniably adapted for fossoriality, all snakes plus their sister group derive from a surface-terrestrial form with non-fossorial behavior, thus redirecting the debate toward an underexplored evolutionary scenario. Our comprehensive heterochrony analyses further indicate that snakes later evolved novel craniofacial specializations through global acceleration of skull development. These results highlight the importance of the interplay between natural selection and developmental processes in snake origin and diversification, leading first to invasion of a new habitat and then to subsequent ecological radiations.

Accumulation of transposable elements in Hox gene clusters during adaptive radiation of Anolis lizards

Transposable elements (TEs) are DNA sequences that can insert elsewhere in the genome and modify genome structure and gene regulation. The role of TEs in evolution is contentious. One hypothesis posits that TE activity generates genomic incompatibilities that can cause reproductive isolation between incipient species. This predicts that TEs will accumulate during speciation events. Here, I tested the prediction that extant lineages with a relatively high rate of speciation have a high number of TEs in their genomes. I sequenced and analysed the TE content of a marker genomic region (Hox clusters) in Anolis lizards, a classic case of an adaptive radiation. Unlike other vertebrates, including closely related lizards, Anolis lizards have high numbers of TEs in their Hox clusters, genomic regions that regulate development of the morphological adaptations that characterize habitat specialists in these lizards. Following a burst of TE activity in the lineage leading to extant Anolis, TEs have continued to accumulate during or after speciation events, resulting in a positive relationship between TE density and lineage speciation rate. These results are consistent with the prediction that TE activity contributes to adaptive radiation by promoting speciation. Although there was no evidence that TE density per se is associated with ecological morphology, the activity of TEs in Hox clusters could have been a rich source for phenotypic variation that may have facilitated the rapid parallel morphological adaptation to microhabitats seen in extant Anolis lizards.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

Tandem amino acid repeats in the green anole (Anolis carolinensis) and other squamates may have a role in increasing genetic variability

Background

Tandem amino acid repeats are characterised by the consecutive recurrence of a single amino acid. They exhibit high rates of length mutations in addition to point mutations and have been proposed to be involved in genetic plasticity. Squamate reptiles (lizards and snakes) diversify in both morphology and physiology. The underlying mechanism is yet to be understood. In a previous phylogenomic analysis of reptiles, the density of tandem repeats in an anole lizard diverged heavily from that of the other reptiles. To gain further insight into the tandem amino acid repeats in squamates, we analysed the repeat content in the green anole (Anolis carolinensis) proteome and compared the amino acid repeats in a large orthologous protein data set from six vertebrates (the Western clawed frog, the green anole, the Chinese softshell turtle, the zebra finch, mouse and human).

Results

Our results revealed that the number of amino acid repeats in the green anole exceeded those found in the other five species studied. Species-only repeats were found in high proportion in the green anole but not in the other five species, suggesting that the green anole had gained many amino acid repeats in either the Anolis or the squamate lineage. Since the amino acid repeat containing genes in the green anole were highly enriched in genes related to transcription and development, an important family of developmental genes, i.e., the Hox family, was further studied in a wide collection of squamates. Abundant amino acid repeats were also observed, implying the general high tolerance of amino acid repeats in squamates. A particular enrichment of amino acid repeats was observed in the central class Hox genes that are known to be responsible for defining cervical to lumbar regions.

Conclusions

Our study suggests that the abundant amino acid repeats in the green anole, and possibly in other squamates, may play a role in increasing the genetic variability, and contribute to the evolutionary diversity of this clade.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2430-y) contains supplementary material, which is available to authorized users.

Hox cluster characterization of Banna caecilian (Ichthyophis bannanicus) provides hints for slow evolution of its genome

BACKGROUND

Caecilians, with a discrete lifestyle, are the least explored group of amphibians. Though with distinct traits, many aspects of their biology are poorly investigated. Obtaining the caecilian genomic sequences will offer new perspectives and aid the fundamental studies in caecilian biology. The caecilian genomic sequences are also important and practical in the comparative genomics of amphibians. Currently, however, only sparse genomic sequences of caecilians are available. Hox genes, an old family of transcription factors playing central roles in the establishment of metazoan body plan. Understanding their structure and genomic organization may provide insights into the animal's genome, which is valuable for animals without a sequenced genome.

RESULTS

We sequenced and characterized the Hox clusters of Banna caecilian (Ichthyophis bannanicus) with a strategy combining long range PCR and genome walking. We obtained the majority of the four caecilian Hox clusters and identified 39 Hox genes, 5 microRNA genes and 1 pseudogene (ψHoxD12). There remained seven intergenic gaps we were unable to fill. From the obtained sequences, the caecilian Hox clusters contained less repetitive sequences and more conserved noncoding elements (CNEs) than the frog counterparts. We found that caecilian and coelacanth shared many more CNEs than frog and coelacanth did. Relative rate of sequence evolution showed that caecilian Hox genes evolved significantly more slowly than the other tetrapod species used in this study and were comparable to the slowly evolving coelacanth Hox genes. Phylogenetic tree of the four Hox clusters also revealed shorter branch length especially for the caecilian HoxA, HoxB and HoxD clusters. These features of the caecilian Hox clusters suggested a slowly evolving genome, which was supported by further analysis of a large orthologous protein dataset.

CONCLUSIONS

Our analyses greatly extended the knowledge about the caecilian Hox clusters from previous PCR surveys. From the obtained Hox sequences and the orthologous protein dataset, the caecilian Hox loci and its genome appear evolving comparatively slowly. As the basal lineage of amphibians and land vertebrate, this characteristic of the caecilian genome is valuable in the study concerning the genome biology and evolution of amphibians and early tetrapods.

Restriction and recruitment-gene duplication and the origin and evolution of snake venom toxins

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive "just-so story" in evolutionary biology.

Hox gene regulation and timing in embryogenesis

Hox genes are critical regulators of embryonic development in bilaterian animals. They exhibit a unique mode of transcriptional regulation where the position of the genes along the chromosome corresponds to the time and place of their expression during development. The sequential temporal activation of these genes in the primitive streak helps determining their subsequent pattern of expression along the anterior-posterior axis of the embryo, yet the precise correspondence between these two collinear processes is not fully understood. In addition, vertebrate Hox genes evolved similar modes of regulation along secondary body axes, such as the developing limbs. We review the current understanding of the mechanisms operating during activation, maintenance and silencing of Hox gene expression in these various contexts, and discuss the evolutionary significance of their genomic organization.

Lizards and LINEs: selection and demography affect the fate of L1 retrotransposons in the genome of the green anole (Anolis carolinensis)

Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila, applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard (Anolis carolinensis). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis. We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

Hox genes regulation in vertebrates

Hox genes encode transcription factors defining cellular identities along the major and secondary body axes. Their coordinated expression in both space and time is critical for embryonic patterning. Accordingly, Hox genes transcription is tightly controlled at multiple levels, and involves an intricate combination of local and long-range cis-regulatory elements. Recent studies revealed that in addition to transcription factors, dynamic patterns of histone marks and higher-order chromatin structure are important determinants of Hox gene regulation. Furthermore, the emerging picture suggests an involvement of various species of non-coding RNA in targeting activating and repressive complexes to Hox clusters. I review these recent developments and discuss their relevance to the control of Hox gene expression in vivo, as well as to our understanding of transcriptional regulatory mechanisms.

Evolution of Hox gene clusters in deuterostomes

Hox genes, with their similar roles in animals as evolutionarily distant as humans and flies, have fascinated biologists since their discovery nearly 30 years ago. During the last two decades, reports on Hox genes from a still growing number of eumetazoan species have increased our knowledge on the Hox gene contents of a wide range of animal groups. In this review, we summarize the current Hox inventory among deuterostomes, not only in the well-known teleosts and tetrapods, but also in the earlier vertebrate and invertebrate groups. We draw an updated picture of the ancestral repertoires of the different lineages, a sort of “genome Hox bar-code” for most clades. This scenario allows us to infer differential gene or cluster losses and gains that occurred during deuterostome evolution, which might be causally linked to the morphological changes that led to these widely diverse animal taxa. Finally, we focus on the challenging family of posterior Hox genes, which probably originated through independent tandem duplication events at the origin of each of the ambulacrarian, cephalochordate and vertebrate/urochordate lineages.

Salamander Hox clusters contain repetitive DNA and expanded non-coding regions: a typical Hox structure for non-mammalian tetrapod vertebrates?

Hox genes encode transcription factors that regulate embryonic and post-embryonic developmental processes. The expression of Hox genes is regulated in part by the tight, spatial arrangement of conserved coding and non-coding sequences. The potential for evolutionary changes in Hox cluster structure is thought to be low among vertebrates; however, recent studies of a few non-mammalian taxa suggest greater variation than originally thought. Using next generation sequencing of large genomic fragments (>100 kb) from the red spotted newt (Notophthalamus viridescens), we found that the arrangement of Hox cluster genes was conserved relative to orthologous regions from other vertebrates, but the length of introns and intergenic regions varied. In particular, the distance between hoxd13 and hoxd11 is longer in newt than orthologous regions from vertebrate species with expanded Hox clusters and is predicted to exceed the length of the entire HoxD clusters (hoxd13–hoxd4) of humans, mice, and frogs. Many repetitive DNA sequences were identified for newt Hox clusters, including an enrichment of DNA transposon-like sequences relative to non-coding genomic fragments. Our results suggest that Hox cluster expansion and transposon accumulation are common features of non-mammalian tetrapod vertebrates.

Evolutionary rate of human tissue-specific genes are related with transposable element insertions

The influence of transposable elements (TEs) on genome evolution has been widely studied. However, it remains unclear whether TE insertions also impact on evolutionary rate of human genes. In this study, we have compared the differences in TEs and evolutionary rates between human tissue-specific genes. Our results showed that various functional categories of human tissue-specific genes contained different TE numbers and divergent values of Ka/Ks, with human nucleic acid binding transcription factor activity genes having the fewest TE density and Ka/Ks value. Interestingly, we also found that human tissue-specific genes with TEs have also undergone faster evolution than those without TEs. Therefore, TEs have significant impact on the evolutionary rates of human tissue-specific genes. Furthermore, local genomic properties such as gene length, GC content and recombination rate may reflect a true transpositional bias for the particular TEs. Our results may provide important insights for further elucidating the evolution of human tissue-specific genes.

From lizard to snake; behind the evolution of an extreme body plan

The elongated, snake-like skeleton, as it has convergently evolved in numerous reptilian and amphibian lineages, is from a developmental biologist’s point of view amongst the most fascinating anatomical peculiarities in the animal kingdom. This type of body plan is characterized by a greatly increased number of vertebrae, a reduction of skeletal regionalization along the primary body axis and loss of the limbs. Recent studies conducted on both mouse and snakes now hint at how changes inside the gene regulatory circuitries of the Hox genes and the somitogenesis clock likely underlie these striking departures from standard tetrapod morphology, suggesting scenarios by which snakes and other elongated species may have evolved from more ordinarily bodied ancestors.

Transposable elements: from DNA parasites to architects of metazoan evolution

One of the most unexpected insights that followed from the completion of the human genome a decade ago was that more than half of our DNA is derived from transposable elements (TEs). Due to advances in high throughput sequencing technologies it is now clear that TEs comprise the largest molecular class within most metazoan genomes. TEs, once categorised as "junk DNA", are now known to influence genomic structure and function by increasing the coding and non-coding genetic repertoire of the host. In this way TEs are key elements that stimulate the evolution of metazoan genomes. This review highlights several lines of TE research including the horizontal transfer of TEs through host-parasite interactions, the vertical maintenance of TEs over long periods of evolutionary time, and the direct role that TEs have played in generating morphological novelty.

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.

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.

An evolving NGF-Hoxd1 signaling pathway mediates development of divergent neural circuits in vertebrates

Species are endowed with unique sensory capabilities encoded by divergent neural circuits. One potential explanation for how divergent circuits have evolved is that conserved extrinsic signals are differentially interpreted by developing neurons of different species to yield unique patterns of axonal connections. Although NGF controls survival, maturation and axonal projections of nociceptors of different vertebrates, whether the NGF signal is differentially transduced in different species to yield unique features of nociceptor circuits is unclear. We identified a species-specific signaling module induced by NGF and mediated by a rapidly evolving Hox transcription factor, Hoxd1. Mice lacking Hoxd1 display altered nociceptor circuitry which resembles that normally found in chicks. Conversely, ectopic expression of Hoxd1 in developing chick nociceptors promotes a pattern of axonal projections reminiscent of the mouse. We propose that conserved growth factors control divergent neuronal transcriptional events which mediate interspecies differences in neural circuits and the behaviors they control.

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.

Evolution of Hox gene clusters in gnathostomes: insights from a survey of a shark (Scyliorhinus canicula) transcriptome

It is now well established that there were four Hox gene clusters in the genome of the last common ancestor of extant gnathostomes. To better understand the evolution of the organization and expression of these genomic regions, we have studied the Hox gene clusters of a shark (Scyliorhinus canicula). We sequenced 225,580 expressed sequence tags from several embryonic cDNA libraries. Blast searches identified corresponding transcripts to almost all the HoxA, HoxB, and HoxD cluster genes. No HoxC transcript was identified, suggesting that this cluster is absent or highly degenerate. Using Hox gene sequences as probes, we selected and sequenced seven clones from a bacterial artificial chromosome library covering the complete region of the three gene clusters. Mapping of cDNAs to these genomic sequences showed extensive alternative splicing and untranslated exon sharing between neighboring Hox genes. Homologous noncoding exons could not be identified in transcripts from other species using sequence similarity. However, by comparing conserved noncoding sequences upstream of these exons in different species, we were able to identify homology between some exons. Some alternative splicing variants are probably very ancient and were already coded for by the ancestral Hox gene cluster. We also identified several transcripts that do not code for Hox proteins, are probably not translated, and all but one are in the reverse orientation to the Hox genes. This survey of the transcriptome of the Hox gene clusters of a shark shows that the high complexity observed in mammals is a gnathostome ancestral feature.

Paucity and preferential suppression of transgenes in late replication domains of the D. melanogaster genome

Background

Eukaryotic genomes are organized in extended domains with distinct features intimately linking genome structure, replication pattern and chromatin state. Recently we identified a set of long late replicating euchromatic regions that are underreplicated in salivary gland polytene chromosomes of D. melanogaster.

Results

Here we demonstrate that these underreplicated regions (URs) have a low density of P-element and piggyBac insertions compared to the genome average or neighboring regions. In contrast, Minos-based transposons show no paucity in URs but have a strong bias to testis-specific genes. We estimated the suppression level in 2,852 stocks carrying a single P-element by analysis of eye color determined by the mini-white marker gene and demonstrate that the proportion of suppressed transgenes in URs is more than three times higher than in the flanking regions or the genomic average. The suppressed transgenes reside in intergenic, genic or promoter regions of the annotated genes. We speculate that the low insertion frequency of P-elements and piggyBacs in URs partially results from suppression of transgenes that potentially could prevent identification of transgenes due to complete suppression of the marker gene. In a similar manner, the proportion of suppressed transgenes is higher in loci replicating late or very late in Kc cells and these loci have a lower density of P-elements and piggyBac insertions. In transgenes with two marker genes suppression of mini-white gene in eye coincides with suppression of yellow gene in bristles.

Conclusions

Our results suggest that the late replication domains have a high inactivation potential apparently linked to the silenced or closed chromatin state in these regions, and that such inactivation potential is largely maintained in different tissues.

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

Hox genes are central to the specification of structures along the anterior-posterior body axis, and modifications in their expression have paralleled the emergence of diversity in vertebrate body plans. Here we describe the genomic organization of Hox clusters in different reptiles and show that squamates have accumulated unusually large numbers of transposable elements at these loci, reflecting extensive genomic rearrangements of coding and non-coding regulatory regions. Comparative expression analyses between two species showing different axial skeletons, the corn snake and the whiptail lizard, revealed major alterations in Hox13 and Hox10 expression features during snake somitogenesis, in line with the expansion of both caudal and thoracic regions. Variations in both protein sequences and regulatory modalities of posterior Hox genes suggest how this genetic system has dealt with its intrinsic collinear constraint to accompany the substantial morphological radiation observed in this group.

Complete HOX cluster characterization of the coelacanth provides further evidence for slow evolution of its genome

The living coelacanth is a lobe-finned fish that represents an early evolutionary departure from the lineage that led to land vertebrates, and is of extreme interest scientifically. It has changed very little in appearance from fossilized coelacanths of the Cretaceous (150 to 65 million years ago), and is often referred to as a "living fossil." An important general question is whether long-term stasis in morphological evolution is associated with stasis in genome evolution. To this end we have used targeted genome sequencing for acquiring 1,612,752 bp of high quality finished sequence encompassing the four HOX clusters of the Indonesian coelacanth Latimeria menadoensis. Detailed analyses were carried out on genomic structure, gene and repeat contents, conserved noncoding regions, and relative rates of sequence evolution in both coding and noncoding tracts. Our results demonstrate conclusively that the coelacanth HOX clusters are evolving comparatively slowly and that this taxon should serve as a viable outgroup for interpretation of the genomes of tetrapod species.

Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes

We have sequenced and analyzed Hox gene clusters from elephant shark, a holocephalian cartilaginous fish. Elephant shark possesses 4 Hox clusters with 45 Hox genes that include orthologs for a higher number of ancient gnathostome Hox genes than the 4 clusters in tetrapods and the supernumerary clusters in teleost fishes. Phylogenetic analysis of elephant shark Hox genes from 7 paralogous groups that contain all of the 4 members indicated an ((AB)(CD)) topology for the order of Hox cluster duplication, providing support for the 2R hypothesis (i.e., 2 rounds of whole-genome duplication during the early evolution of vertebrates). Comparisons of noncoding sequences of the elephant shark and human Hox clusters have identified a large number of conserved noncoding elements (CNEs), which represent putative cis-regulatory elements that may be involved in the regulation of Hox genes. Interestingly, in fugu more than 50% of these ancient CNEs have diverged beyond recognition in the duplicated (HoxA, HoxB, and HoxD) as well as the singleton (HoxC) Hox clusters. Furthermore, the b-paralogs of the duplicated fugu Hox clusters are virtually devoid of unique ancient CNEs. In contrast to fugu Hox clusters, elephant shark and human Hox clusters have lost fewer ancient CNEs. If these ancient CNEs are indeed enhancers directing tissue-specific expression of Hox genes, divergence of their sequences in vertebrate lineages might have led to altered expression patterns and presumably the functions of their associated Hox genes.

Comparative genomics of chondrichthyan Hoxa clusters

BACKGROUND

The chondrichthyan or cartilaginous fish (chimeras, sharks, skates and rays) occupy an important phylogenetic position as the sister group to all other jawed vertebrates and as an early lineage to diverge from the vertebrate lineage following two whole genome duplication events in vertebrate evolution. There have been few comparative genomic analyses incorporating data from chondrichthyan fish and none comparing genomic information from within the group. We have sequenced the complete Hoxa cluster of the Little Skate (Leucoraja erinacea) and compared to the published Hoxa cluster of the Horn Shark (Heterodontus francisci) and to available data from the Elephant Shark (Callorhinchus milii) genome project.

RESULTS

A BAC clone containing the full Little Skate Hoxa cluster was fully sequenced and assembled. Analyses of coding sequences and conserved non-coding elements reveal a strikingly high level of conservation across the cartilaginous fish, with twenty ultraconserved elements (100%,100 bp) found between Skate and Horn Shark, compared to three between human and marsupials. We have also identified novel potential non-coding RNAs in the Skate BAC clone, some of which are conserved to other species.

CONCLUSION

We find that the Little Skate Hoxa cluster is remarkably similar to the previously published Horn Shark Hoxa cluster with respect to sequence identity, gene size and intergenic distance despite over 180 million years of separation between the two lineages. We suggest that the genomes of cartilaginous fish are more highly conserved than those of tetrapods or teleost fish and so are more likely to have retained ancestral non-coding elements. While useful for isolating homologous DNA, this complicates bioinformatic approaches to identify chondrichthyan-specific non-coding DNA elements.