Sauria SINEs: Novel short interspersed retroposable elements that are widespread in reptile genomes

The Satellite DNA PcH-Sat, Isolated and Characterized in the Limpet Patella caerulea (Mollusca, Gastropoda), Suggests the Origin from a Nin-SINE Transposable Element

Satellite DNA (sat-DNA) was previously described as junk and selfish DNA in the cellular economy, without a clear functional role. However, during the last two decades, evidence has been accumulated about the roles of sat-DNA in different cellular functions and its probable involvement in tumorigenesis and adaptation to environmental changes. In molluscs, studies on sat-DNAs have been performed mainly on bivalve species, especially those of economic interest. Conversely, in Gastropoda (which includes about 80% of the currently described molluscs species), studies on sat-DNA have been largely neglected. In this study, we isolated and characterized a sat-DNA, here named PcH-sat, in the limpet Patella caerulea using the restriction enzyme method, particularly HaeIII. Monomeric units of PcH-sat are 179 bp long, AT-rich (58.7%), and with an identity among monomers ranging from 91.6 to 99.8%. Southern blot showed that PcH-sat is conserved in P. depressa and P. ulyssiponensis, while a smeared signal of hybridization was present in the other three investigated limpets (P. ferruginea, P. rustica and P. vulgata). Dot blot showed that PcH-sat represents about 10% of the genome of P. caerulea, 5% of that of P. depressa, and 0.3% of that of P. ulyssiponensis. FISH showed that PcH-sat was mainly localized on pericentromeric regions of chromosome pairs 2 and 4-7 of P. caerulea (2n = 18). A database search showed that PcH-sat contains a large segment (of 118 bp) showing high identity with a homologous trait of the Nin-SINE transposable element (TE) of the patellogastropod Lottia gigantea, supporting the hypothesis that TEs are involved in the rising and tandemization processes of sat-DNAs.

Interspecific Comparison of Orthologous Short Interspersed Elements Loci Using Whole-Genome Data

The polymorphism of SINE-containing loci reflects the evolutionary processes that occurred both during the period before the divergence of the taxa and after it. Orthologous loci containing SINE in two or more genomes indicate the relatedness of the taxa, while different copies may have a specific set of mutations and degree of difference. Polymorphic insertion can be interpreted with a high degree of confidence as a shared derived character in the phylogenetic reconstruction of the history of the taxon. The computational comparison of the entire set of SINE-containing loci between genomes is a challenging task, and we propose to consider it in detail using the genomes of representatives of squamate reptiles (lizards) as an example. Our approach allows us to extract copies of SINE from the genomes, find pairwise orthologous loci by using flanking genomic sequences, and analyze the resulting sets of loci for the presence or absence of SINE, the degree of similarity of the flanks, and the similarity of the SINE themselves. The workflow we propose allows us to efficiently extract and analyze orthologous SINE loci for the downstream analysis, as shown in our comparison of species- and genus-level taxa in lacertid lizards.

Characterization of Two Transposable Elements and an Ultra-Conserved Element Isolated in the Genome of Zootoca vivipara (Squamata, Lacertidae)

Transposable elements (TEs) constitute a considerable fraction of eukaryote genomes representing a major source of genetic variability. We describe two DNA sequences isolated in the lizard Zootoca vivipara, here named Zv516 and Zv817. Both sequences are single-copy nuclear sequences, including a truncation of two transposable elements (TEs), SINE Squam1 in Zv516 and a Tc1/Mariner-like DNA transposon in Zv817. FISH analyses with Zv516 showed the occurrence of interspersed signals of the SINE Squam1 sequence on all chromosomes of Z. vivipara and quantitative dot blot indicated that this TE is present with about 4700 copies in the Z. vivipara genome. FISH and dot blot with Zv817 did not produce clear hybridization signals. Bioinformatic analysis showed the presence of active SINE Squam 1 copies in the genome of different lacertids, in different mRNAs, and intronic and coding regions of various genes. The Tc1/Mariner-like DNA transposon occurs in all reptiles, excluding Sphenodon and Archosauria. Zv817 includes a trait of 284 bp, representing an amniote ultra-conserved element (UCE). Using amniote UCE homologous sequences from available whole genome sequences of major amniote taxonomic groups, we performed a phylogenetic analysis which retrieved Prototheria as the sister group of Metatheria and Eutheria. Within diapsids, Testudines are the sister group to Aves + Crocodylia (Archosauria), and Sphenodon is the sister group to Squamata. Furthermore, large trait regions flanking the UCE are conserved at family level.

Differences in Homomorphic Sex Chromosomes Are Associated with Population Divergence in Sex Determination in Carinascincus ocellatus (Scincidae: Lygosominae)

Sex determination directs development as male or female in sexually reproducing organisms. Evolutionary transitions in sex determination have occurred frequently, suggesting simple mechanisms behind the transitions, yet their detail remains elusive. Here we explore the links between mechanisms of transitions in sex determination and sex chromosome evolution at both recent and deeper temporal scales (<1 Myr; ~79 Myr). We studied a rare example of a species with intraspecific variation in sex determination, Carinascincus ocellatus, and a relative, Liopholis whitii, using c-banding and mapping of repeat motifs and a custom Y chromosome probe set to identify the sex chromosomes. We identified both unique and conserved regions of the Y chromosome among C. ocellatus populations differing in sex determination. There was no evidence for homology of sex chromosomes between C. ocellatus and L. whitii, suggesting independent evolutionary origins. We discuss sex chromosome homology between members of the subfamily Lygosominae and propose links between sex chromosome evolution, sex determination transitions, and karyotype evolution.

Novel Approaches for Phylogenetic Inference from Morphological Data and Total-Evidence Dating in Squamate Reptiles (Lizards, Snakes, and Amphisbaenians)

Here, I combine previously underutilized models and priors to perform more biologically realistic phylogenetic inference from morphological data, with an example from squamate reptiles. When coding morphological characters, it is often possible to denote ordered states with explicit reference to observed or hypothetical ancestral conditions. Using this logic, we can integrate across character-state labels and estimate meaningful rates of forward and backward transitions from plesiomorphy to apomorphy. I refer to this approach as MkA, for “asymmetric.” The MkA model incorporates the biological reality of limited reversal for many phylogenetically informative characters, and significantly increases likelihoods in the empirical data sets. Despite this, the phylogeny of Squamata remains contentious. Total-evidence analyses using combined morphological and molecular data and the MkA approach tend toward recent consensus estimates supporting a nested Iguania. However, support for this topology is not unambiguous across data sets or analyses, and no mechanism has been proposed to explain the widespread incongruence between partitions, or the hidden support for various topologies in those partitions. Furthermore, different morphological data sets produced by different authors contain both different characters and different states for the same or similar characters, resulting in drastically different placements for many important fossil lineages. Effort is needed to standardize ontology for morphology, resolve incongruence, and estimate a robust phylogeny. The MkA approach provides a preliminary avenue for investigating morphological evolution while accounting for temporal evidence and asymmetry in character-state changes.

[STRUCTURE AND DISTRIBUTION OF THE RETROTRANSPOSON BOV-B LINE]

The classification of mobile elements was discussed. Special attention was devoted to the retroelement of the LINE group: retrotransposon Bov-B LINE. The history of its origin and distribution in the nature was considered. The results of the phenomenon of horizontal transition of the retrotransposon Bov-B LINE between evolutionally distant classes were discussed.

The Organization of Repetitive DNA in the Genomes of Amazonian Lizard Species in the Family Teiidae

Repetitive DNA is the largest fraction of the eukaryote genome and comprises tandem and dispersed sequences. It presents variations in relation to its composition, number of copies, distribution, dynamics, and genome organization, and participates in the evolutionary diversification of different vertebrate species. Repetitive sequences are usually located in the heterochromatin of centromeric and telomeric regions of chromosomes, contributing to chromosomal structures. Therefore, the aim of this study was to physically map repetitive DNA sequences (5S rDNA, telomeric sequences, tropomyosin gene 1, and retroelements Rex1 and SINE) of mitotic chromosomes of Amazonian species of teiids (Ameiva ameiva, Cnemidophorus sp. 1, Kentropyx calcarata, Kentropyx pelviceps, and Tupinambis teguixin) to understand their genome organization and karyotype evolution. The mapping of repetitive sequences revealed a distinct pattern in Cnemidophorus sp. 1, whereas the other species showed all sequences interspersed in the heterochromatic region. Physical mapping of the tropomyosin 1 gene was performed for the first time in lizards and showed that in addition to being functional, this gene has a structural function similar to the mapped repetitive elements as it is located preferentially in centromeric regions and termini of chromosomes.

A Deluge of Complex Repeats: The Solanum Genome

Repetitive elements have lately emerged as key components of genome, performing varieties of roles. It has now become necessary to have an account of repeats for every genome to understand its dynamics and state. Recently, genomes of two major Solanaceae species, Solanum tuberosum and Solanum lycopersicum, were sequenced. These species are important crops having high commercial significance as well as value as model species. However, there is a reasonable gap in information about repetitive elements and their possible roles in genome regulation for these species. The present study was aimed at detailed identification and characterization of complex repetitive elements in these genomes, along with study of their possible functional associations as well as to assess possible transcriptionally active repetitive elements. In this study, it was found that ~50-60% of genomes of S. tuberosum and S. lycopersicum were composed of repetitive elements. It was also found that complex repetitive elements were associated with >95% of genes in both species. These two genomes are mostly composed of LTR retrotransposons. Two novel repeat families very similar to LTR/ERV1 and LINE/RTE-BovB have been reported for the first time. Active existence of complex repeats was estimated by measuring their transcriptional abundance using Next Generation Sequencing read data and Microarray platforms. A reasonable amount of regulatory components like transcription factor binding sites and miRNAs appear to be under the influence of these complex repetitive elements in these species, while several genes appeared to possess exonized repeats.

Genome-wide analysis of short interspersed nuclear elements SINES revealed high sequence conservation, gene association and retrotranspositional activity in wheat

Short interspersed nuclear elements (SINEs) are non-autonomous non-LTR retroelements that are present in most eukaryotic species. While SINEs have been intensively investigated in humans and other animal systems, they are poorly studied in plants, especially in wheat (Triticum aestivum). We used quantitative PCR of various wheat species to determine the copy number of a wheat SINE family, termed Au SINE, combined with computer-assisted analyses of the publicly available 454 pyrosequencing database of T. aestivum. In addition, we utilized site-specific PCR on 57 Au SINE insertions, transposon methylation display and transposon display on newly formed wheat polyploids to assess retrotranspositional activity, epigenetic status and genetic rearrangements in Au SINE, respectively. We retrieved 3706 different insertions of Au SINE from the 454 pyrosequencing database of T. aestivum, and found that most of the elements are inserted in A/T-rich regions, while approximately 38% of the insertions are associated with transcribed regions, including known wheat genes. We observed typical retrotransposition of Au SINE in the second generation of a newly formed wheat allohexaploid, and massive hypermethylation in CCGG sites surrounding Au SINE in the third generation. Finally, we observed huge differences in the copy numbers in diploid Triticum and Aegilops species, and a significant increase in the copy numbers in natural wheat polyploids, but no significant increase in the copy number of Au SINE in the first four generations for two of three newly formed allopolyploid species used in this study. Our data indicate that SINEs may play a prominent role in the genomic evolution of wheat through stress-induced activation.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, L-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

Integrated "omics" profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus

BACKGROUND

Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach.

RESULTS

Proteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a β-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9-18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulthood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus.

CONCLUSIONS

Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112-125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.

SINEBase: a database and tool for SINE analysis

SINEBase (http://sines.eimb.ru) integrates the revisited body of knowledge about short interspersed elements (SINEs). A set of formal definitions concerning SINEs was introduced. All available sequence data were screened through these definitions and the genetic elements misidentified as SINEs were discarded. As a result, 175 SINE families have been recognized in animals, flowering plants and green algae. These families were classified by the modular structure of their nucleotide sequences and the frequencies of different patterns were evaluated. These data formed the basis for the database of SINEs. The SINEBase website can be used in two ways: first, to explore the database of SINE families, and second, to analyse candidate SINE sequences using specifically developed tools. This article presents an overview of the database and the process of SINE identification and analysis.

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.

First draft of the genomic organization of a PIII-SVMP gene

The evolutionary pathway of highly toxic proteins expressed in snake venom glands from proteins without toxic function and expressed in non-parotid tissues remains poorly understood. Here we examine gene structure of a representative of a venom protein with an ADAMs metalloproteinase evolutionary origin. The structure of the 15,652 bp Echis ocellatus pre-pro EOC00089-like PIII-SVMP gene was assembled from PCR-amplified sequences of overlapping genomic fragments. The gene comprises 12 exons interrupted by 11 introns. In a homology model of the EOC00089-like protein, the insertion of introns interrupting coding regions lie just after or between secondary structure elements. Long interspersed nuclear retroelements (LINE) L2/CR1 and RTE/Bov-B, short interspersed nuclear retroelements SINE/Sauria, and a hobo-activator DNA (Charlie, hAT) transposon were identified within introns 1, 3, 7 and 8. Pairwise amino acid sequence comparisons between EOC00089-like PIII-SVMP and its closest orthologs, ADAM28, from a mammal, Homo sapiens, and the lizard, Anolis carolinensis, showed that the ORFs of these three proteins share 42%/59%, 49%/69%, and 48%/65% (identity/similarity), respectively. The protein-coding positions interrupted by each of the 11 introns of the Echis PIII-SVMP gene are entirely conserved in the A. carolinensis and human ADAM28 genes. However, the lizard and the human ADAM28 genes contain 5 introns not present in the E. ocellatus gene. Furthermore, Echis and Anolis introns exhibit quantitatively and qualitatively distinctions in their inserted retroelements. These findings identify introns as possible key elements in the recruitment and amplification process of SVMPs into the venom gland of extant snakes. Ongoing reptile genome sequencing projects may shed light on this intriguing aspect of the emergence and evolution of venom toxin genes. Furthermore, the organization of the PIII-SVMP reported here provides a genomic explanation for the emergence of dimeric disintegrin subunits encoded by short messengers.

The structural and functional diversification of the Toxicofera reptile venom system

The evolutionary origin and diversification of the reptilian venom system is described. The resolution of higher-order molecular phylogenetics has clearly established that a venom system is ancestral to snakes. The diversification of the venom system within lizards is discussed, as is the role of venom delivery in the behavioural ecology of these taxa (particularly Varanus komodoensis). The more extensive diversification of the venom system in snakes is summarised, including its loss in some clades. Finally, we discuss the contentious issue of a definition for "venom", supporting an evolutionary definition that recognises the homology of both the venom delivery systems and the toxins themselves.

Identification and characterization of a new member of the SINE Au retroposon family (GmAu1) in the soybean, Glycine max (L.) Merr., genome and its potential application

A plant short interspersed element (SINE) was identified in Glycine max after re-sequencing of the soybean sequence characterized amplified region (SCAR) markers. Detailed analysis revealed that this newly recognized SINE element consisted of a tRNA-related region, a tRNA non-related region, direct flanking repeat sequences, and a short stretch of Ts at the 3'-terminal region. These features are similar to previously characterized SINEs. To investigate the evolution of the SINE retroposon, BLASTN was used to search against genome sequences of other plants. Since it is homologous with the retroposon Au in Aegilops umbellulata (wheat) and its homology in soybean, the SINE is named as GmAu1. Genome analysis of the Glycine max var. Willimas 82 uncovered more than 847 copies of GmAu1 per haploid genome of soybean. Examination of the regions flanking the inserted GmAu1 sequences indicated a preference for introns over exons or other noncoding regions. Considering the flanking insertion sequences, 146 primers were designed in order to detect insertion mutations by a PCR-based method. Seventy-seven primers displayed polymorphism and were used to develop corresponding GmAu1-based SCAR markers. The retroposon GmAu1 and its related SCAR markers identified in this study will prove valuable to future investigations into the genetic mapping, phylogeny, and evolution of the Glycine genus.

Novel SINEs families in Medicago truncatula and Lotus japonicus: bioinformatic analysis

Although short interspersed elements (SINEs) were discovered nearly 30 years ago, the studies of these genomic repeats were mostly limited to animal genomes. Very little is known about SINEs in legumes--one of the most important plant families. Here we report identification, genomic distribution and molecular features of six novel SINE elements in Lotus japonicus (named LJ_SINE-1, -2, -3) and Medicago truncatula (MT_SINE-1, -2, -3), model species of legume. They possess all the structural features commonly found in short interspersed elements including RNA polymerase III promoter, polyA tail and flanking repeats. SINEs described here are present in low to moderate copy numbers from 150 to 3000. Bioinformatic analyses were used to searched public databases, we have shown that three of new SINE elements from M. truncatula seem to be characteristic of Medicago and Trifolium genera. Two SINE families have been found in L. japonicus and one is present in both M. truncatula and L. japonicus. In addition, we are discussing potential activities of the described elements.

Short interspersed elements (SINEs) of squamate reptiles (Squam1 and Squam2): structure and phylogenetic significance

Short interspersed elements (SINEs) are important nuclear molecular markers of the evolution of many eukaryotes. However, the SINEs of squamate reptile genomes have been little studied. We first identified two families of SINEs, termed Squam1 and Squam2, in the DNA of meadow lizard Darevskia praticola (Lacertidae) by performing DNA hybridization and PCR. Later, the same families of retrotransposons were found using the same methods in members of another 25 lizard families (from Iguania, Scincomorpha, Gekkota, Varanoidea, and Diploglossa infraorders) and two snake families, but their abundances in these taxa varied greatly. Both SINEs were Squamata-specific and were absent from mammals, birds, crocodiles, turtles, amphibians, and fish. Squam1 possessed some characteristics common to tRNA-related SINEs from fish and mammals, while Squam2 belonged to the tRNA(Ala) group of SINEs and had a more unusual and divergent structure. Squam2-related sequences were found in several unannotated GenBank sequences of squamate reptiles. Squam1 abundance in the Polychrotidae, Agamidae, Leiolepididae, Chamaeleonidae, Scincidae, Lacertidae, Gekkonidae, Varanidae, Helodermatidae, and two snake families were 10(2) -10(4) times higher than those in other taxa (Corytophanidae, Iguanidae, Anguidae, Cordylidae, Gerrhosauridae, Pygopodidae, and Eublepharidae). A less dramatic degree of copy number variation was observed for Squam2 in different taxa. Several Squam1 copies from Lacertidae, Chamaeleonidae, Gekkonidae, Varanidae, and Colubridae were sequenced and found to have evident orthologous features, as well as taxa-specific autapomorphies. Squam1 from Lacertidae and Chamaeleonidae could be divided into several subgroups based on sequence differences. Possible applications of these SINEs as Squamata phylogeny markers are discussed.

The evolution and diversity of DNA transposons in the genome of the Lizard Anolis carolinensis

DNA transposons have considerably affected the size and structure of eukaryotic genomes and have been an important source of evolutionary novelties. In vertebrates, DNA transposons are discontinuously distributed due to the frequent extinction and recolonization of these genomes by active elements. We performed a detailed analysis of the DNA transposons in the genome of the lizard Anolis carolinensis, the first non-avian reptile to have its genome sequenced. Elements belonging to six of the previously recognized superfamilies of elements (hAT, Tc1/Mariner, Helitron, PIF/Harbinger, Polinton/Maverick, and Chapaev) were identified. However, only four (hAT, Tc1/Mariner, Helitron, and Chapaev) of these superfamilies have successfully amplified in the anole genome, producing 67 distinct families. The majority (57/67) are nonautonomous and demonstrate an extraordinary diversity of structure, resulting from frequent interelement recombination and incorporation of extraneous DNA sequences. The age distribution of transposon families differs among superfamilies and reveals different dynamics of amplification. Chapaev is the only superfamily to be extinct and is represented only by old copies. The hAT, Tc1/Mariner, and Helitron superfamilies show different pattern of amplification, yet they are predominantly represented by young families, whereas divergent families are exceedingly rare. Although it is likely that some elements, in particular long ones, are subjected to purifying selection and do not reach fixation, the majority of families are neutral and accumulate in the anole genome in large numbers. We propose that the scarcity of old copies in the anole genome results from the rapid decay of elements, caused by a high rate of DNA loss.

BmSE, a SINE family with 3' ends of (ATTT) repeats in domesticated silkworm (Bombyx mori)

Short interspersed elements (SINEs), which are mainly composed of Bm1, are abundant in the domesticated silkworm. A 294 bp novel SINE family, designated as BmSE, was identified by mining the database of the complete Bombyx mori genome. A representational BmSE element is flanked by an 11 bp target site duplication sequence posterior poly (A) at the 3' end and has the sequence motifs of an internal promoter of RNA polymerase III, which are similar to that of Bm1. The repetitive elements of BmSE are widely distributed in all 28 chromosomes of the genome and share the common (ATTT) repeats at the ends. GC-content distribution shows that BmSE tends to accumulate preferably in the region of higher AT content than that of Bm1. A high proportion of the BmSEs are mapped to the coding sequence introns, whereas several elements are also present in the UTR of some transcripts, indicating that BmSEs are indeed exonized with UTRs. Of the 615 identified structural variants (SVs) of BmSE among the 40 domesticated and wild silkworms, only 230 SVs were found in the domesticated silkworms, indicating that many recent SV events of BmSE occurred after domestication, which was probably due to its mobilization. Our analysis might assist in developing BmSE as a potential marker and in understanding the evolutionary roles of SINEs in the domesticated silkworm.

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

Hox genes control many aspects of embryonic development in metazoans. Previous analyses of this gene family revealed a surprising diversity in terms of gene number and organization between various animal species. In vertebrates, Hox genes are grouped into tightly organized clusters, claimed to be devoid of repetitive sequences. Here, we report the genomic organization of the four Hox loci present in the green anole lizard and show that they have massively accumulated retrotransposons, leading to gene clusters larger in size when compared to other vertebrates. In addition, similar repeats are present in many other development-related gene-containing regions, also thought to be refractory to such repetitive elements. Transposable elements are major sources of genetic variations, including alterations of gene expression, and hence this situation, so far unique among vertebrates, may have been associated with the evolution of the spectacular realm of morphological variations in the body plans of Squamata. Finally, sequence alignments highlight some divergent evolution in highly conserved DNA regions between vertebrate Hox clusters, which may coincide with the emergence of mammalian-specific features.

The evolution of two partner LINE/SINE families and a full-length chromodomain-containing Ty3/Gypsy LTR element in the first reptilian genome of Anolis carolinensis

Transposable elements have been characterized in a number of vertebrates, including whole genomes of mammals, birds, and fishes. The Anolis carolinensis draft assembly provides the first opportunity to study retroposons in a reptilian genome. Here, we identified and reconstructed a number of retroposons based on database searches: Five Sauria short interspersed element (SINE) subfamilies, 5S-Sauria SINE chimeras, Anolis Bov-B long interspersed element (LINE), Anolis SINE 2, Anolis LINE 2, Anolis LINE 1, Anolis CR 1, and a chromodomain-containing Ty3/Gypsy LTR element. We focused on two SINE families (Anolis Sauria SINE and Anolis SINE 2) and their partner LINE families (Anolis Bov-B LINE and Anolis LINE 2). We demonstrate that each SINE/LINE pair is distributed similarly and predict that the retrotransposition of evolutionarily younger Sauria SINE members is via younger Bov-B LINE members while a correlation also exists between their respective evolutionarily older SINE/LINE members. The evolutionarily youngest Sauria SINE sequences evolved as part of novel rolling-circle transposons. The evolutionary time frame when Bov-B LINEs and Sauria SINEs were less active in their retrotransposition is characterized by a high retrotransposition burst of Anolis SINE 2 and Anolis LINE 2 elements. We also characterized the first full-length chromoviral LTR element in amniotes (Amn-ichi). This newly identified chromovirus is widespread in the Anolis genome and has been very well preserved, indicating that it is still active. Transposable elements in the Anolis genome account for approximately 20% of the total DNA sequence, whereas the proportion is more than double that in many mammalian genomes in which such elements have important biological functions. Nevertheless, 20% transposable element coverage is sufficient to predict that Anolis retroposons and other mobile elements also may have biologically and evolutionarily relevant functions. The new SINEs and LINEs and other ubiquitous genomic elements characterized in the Anolis genome will prove very useful for studies in comparative genomics, phylogenetics, and functional genetics.

Retroposons of salmonoid fishes (Actinopterygii: Salmonoidei) and their evolution

Short and long retroposons, or non-LTR retrotransposons (SINEs and LINEs, respectively) are two groups of interspersed repetitive elements amplifying in the genome via RNA and cDNA-mediated reverse transcription. In this process, SINEs entirely depend on the enzymatic machinery of autonomous LINEs. The impact of retroposons on the host genome is difficult to overestimate: their sequences account for significant portion of the eukaryotic genome, while propagation of their active copies gradually reshapes it. In this way, the retropositional activity plays a role of important evolutionary factor. More than 100 LINE and nearly 100 SINE families have been described to date from the genomes of various eukaryotes, and it is salmonoid fishes (Actinopterygii: Salmonoidei) that are particularly noticeable for the diversity of transposons they host in their genomes, including two LINE and seven SINE families. Moreover, this group of ray-finned fish represents an excellent opportunity to study such a rare evolutionary phenomenon as lateral gene transfer, due to a great variety of transposons and other sequences salmons share with a blood fluke, Schistosoma japonicum (Trematoda: Strigeiformes)--a parasitic helminth infecting various vertebrates. The aim of the present review is to structure all knowledge accumulated about salmonoid retroposons by now, as well as to complement it with the new data pertaining to the distribution of some SINE families.

Bov-B-mobilized SINEs in vertebrate genomes

Two new short retroposon families (SINEs) have been found in the genome of springhare Pedetes capensis (Rodentia). One of them, Ped-1, originated from 5S rRNA, while the other one, Ped-2, originated from tRNA-derived SINE ID. In contrast to most currently active mammalian SINEs mobilized by L1 long retrotransposon (LINE), Ped-1 and Ped-2 are mobilized by Bov-B, a LINE family of the widely distributed RTE clade. The 3' part of these SINEs originates from two sequences in the 5' and 3' regions of Bov-B. Such bipartite structure of the LINE-derived part has been revealed in all Bov-B-mobilized SINEs known to date (AfroSINE, Bov-tA, Mar-1, and Ped-1/2), which distinguishes them from other SINEs with only a 3' LINE-derived part. Structural analysis and the distribution of Bov-B LINEs and partner SINEs supports the horizontal transfer of Bov-B, while the SINEs emerged independently in lineages with this LINE.

Poxviruses as possible vectors for horizontal transfer of retroposons from reptiles to mammals

Poxviruses (Poxviridae) are a family of double-stranded DNA viruses with no RNA stage. Members of the genus Orthopoxvirus (OPV) are highly invasive and virulent. It was recently shown that the taterapox virus (TATV) from a West African rodent is the sister of camelpox virus and therefore belongs to the clade closest to the variola virus (VARV), the etiological agent of smallpox. Although these OPVs are among the most dreaded pathogens on Earth, our current knowledge of their genomes, their origins, and their possible hosts is still very limited. Here, we report the horizontal transfer of a retroposon (known only from reptilian genomes) to the TATV genome. After isolating and analyzing different subfamilies of short interspersed elements (SINEs) from lizards and snakes, we identified a highly poisonous snake (Echis ocellatus) from West Africa as the closest species from which the SINE sequence discovered in the TATV genome (TATV-SINE) was transferred to the virus. We discovered direct repeats derived from the virus flanking the TATV-SINE, and the absence of any snake-derived DNA flanking the SINE. These data provide strong evidence that the TATV-SINE was actually transferred within the snake to the viral genome by retrotransposition and not by any horizontal transfer at the DNA level. We propose that the snake is another host for TATV, suggesting that VARV-related epidemiologically relevant viruses may have derived from our cold-blooded ancestors and that poxviruses are possible vectors for horizontal transfer of retroposons from reptiles to mammals.

Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.

Extensive Morphological Convergence and Rapid Radiation in the Evolutionary History of the Family Geoemydidae (Old World Pond Turtles) Revealed by SINE Insertion Analysis

The family Geoemydidae is one of three in the superfamily Testudinoidea and is the most diversified family of extant turtle species. The phylogenetic relationships in this family and among related families have been vigorously investigated from both morphological and molecular viewpoints. The evolutionary history of Geoemydidae, however, remains controversial. Therefore, to elucidate the phylogenetic relationships of Geoemydidae and related species, we applied the SINE insertion method to investigate 49 informative SINE loci in 28 species. We detected four major evolutionary lineages (Testudinidae, Batagur group, Siebenrockiella group, and Geoemyda group) in the clade Testuguria (a clade of Geoemydidae + Testudinidae). All five specimens of Testudinidae form a monophyletic clade. The Batagur group comprises five batagurines. The Siebenrockiella group has one species, Siebenrockiella crassicollis. The Geoemyda group comprises 15 geoemydines (including three former batagurines, Mauremys reevesii, Mauremys sinensis, and Heosemys annandalii). Among these four groups, the SINE insertion patterns were inconsistent at four loci, suggesting that an ancestral species of Testuguria radiated and rapidly diverged into the four lineages during the initial stage of its evolution. Furthermore, within the Geoemyda group we identified three evolutionary lineages, namely Mauremys, Cuora, and Heosemys. The Heosemys lineage comprises Heosemys, Sacalia, Notochelys, and Melanochelys species, and its monophyly is a novel assemblage in Geoemydidae. Our SINE phylogenetic tree demonstrates extensive convergent morphological evolution between the Batagur group and the three species of the Geoemyda group, M. reevesii, M. sinensis, and H. annandalii.

Phylogenomic Investigation of CR1 LINE Diversity in Reptiles

It is unlikely that taxonomically diverse phylogenetic studies will be completed rapidly in the near future for nonmodel organisms on a whole-genome basis. However, one approach to advancing the field of "phylogenomics" is to estimate the structure of poorly known genomes by mining libraries of clones from suites of taxa, rather than from single species. The present analysis adopts this approach by taking advantage of megabase-scale end-sequence scanning of reptilian genomic clones to characterize diversity of CR1-like LINEs, the dominant family of transposable elements (TEs) in the sister group of mammals. As such, it helps close an important gap in the literature on the molecular systematics and evolution of retroelements in nonavian reptiles. Results from aligning more than 14 Mb of sequence from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta), Bahamian green anole (Anolis smaragdinus), Tuatara (Sphenodon punctatus), Emu (Dromaius novaehollandiae), and Zebra Finch (Taeniopygia guttata) against a comprehensive library approximately 3000 TE-encoding peptides reflect an increasing abundance of LINE and non-long-terminal-repeat (non-LTR) retrotransposon repeat types with the age of common ancestry among exemplar reptilian clades. The hypothesis that repeat diversity is correlated with basal metabolic rate was tested using comparative methods and a significant nonlinear relationship was indicated. This analysis suggests that the age of divergence between an exemplary clade and its sister group as well as metabolic correlates should be considered in addition to genome size in explaining patterns of retroelement diversity. The first phylogenetic analysis of the largely unexplored chicken repeat 1 (CR1) 3' reverse transcriptase (RT) conserved domains 8 and 9 in nonavian reptiles reveals a pattern of multiple lineages with variable branch lengths, suggesting presence of both old and young elements and the existence of several distinct well-supported clades not apparent from previous characterization of CR1 subfamily structure in birds and the turtle. This mode of CR1 evolution contrasts with historical patterns of LINE 1 diversification in mammals and hints toward the existence of a rich but still largely unexplored diversity of nonavian retroelements of importance to advancing both comparative vertebrate genomics and amniote systematics.

Toothed whale monophyly reassessed by SINE insertion analysis: the absence of lineage sorting effects suggests a small population of a common ancestral species

Morphological data have indicated that toothed whales form a monophyletic group. However, research published in the last several years has made the issue of the monophyly or paraphyly of toothed whales a subject of debate. Our group previously characterized three independent loci in which SINE insertions were shared among dolphins and sperm whales, thus supporting the traditional, morphologically based hypothesis of toothed whale monophyly. Although in recent years a few additional molecular works proposed this topology, there is still skepticism over this monophyly from the view point of molecular systematics. When the phylogeny of rapidly radiated taxa is examined using the SINE method, it is important to consider the ascertainment bias that arises when choosing a particular taxon for SINE loci screening. To overcome this methodological problem specific to the SINE method, we examined all possible topologies among sperm whales, dolphins and baleen whales by extensively screening SINE loci from species of all three lineages. We characterized nine independent SINE loci from the genomes of sperm whales and dolphins, all of which cluster sperm whales and dolphins but exclude baleen whales. Furthermore, we characterized ten independent loci from baleen whales, all of which were amplified in a common ancestor of these whales. From these observations, we conclude that toothed whales form a monophyletic group and that no ancestral SINE polymorphisms hinder their phylogenetic assignment despite the short divergence times of the major lineages of extant whales during evolution. These results suggest that a small population of common ancestors of all toothed whales ultimately diverged into the lineages of sperm whales and dolphins.