Finding of Bov-B LINE retroelement in parthenogenetic and bisexual lizard species of the genus Darevskia (Lacertidae)

The Bov-B LINE retrotransposon was first discovered in Ruminantia and was long considered to be specific for this order. Later, this mobile element was described in snakes and some lizard species. Analysis of phylogenetic relationships of Bov-B LINE elements from different ruminants, snakes, and lizard species led to the suggestion on horizontal transfer of this retrotransposon from Squamata to Ruminantia. In the Squamata group, Bov-B LINE element was found in all snakes and some lizard species examined. The element was not detected in the genomes of some species of the genera Lacerta and Podarcis. In the present study, using PCR amplification and sequencing of PCR products, Bov-B LINE element was identified in the genomes of parthenogenetic and bisexual species of the genus Darevskia (Lacertidae), as well as in such species as Lacerta agilis and Zootoca vivipara, where this retrotransposon had not been not detected before.     

Molecular evolution of Bov-B LINEs in vertebrates

Since their discovery in family Bovidae (bovids), Bov-B LINEs, believed to be order-specific SINEs, have been found in all ruminants and recently also in Viperidae snakes. The distribution and the evolutionary relationships of Bov-B LINEs provide an indication of their origin and evolutionary dynamics in different species. The evolutionary origin of Bov-B LINE elements has been shown unequivocally to be in Squamata (squamates). The horizontal transfer of Bov-B LINE elements in vertebrates has been confirmed by their discontinuous phylogenetic distribution in Squamata (Serpentes and two lizard infra-orders) as well as in Ruminantia, by the high level of nucleotide identity, and by their phylogenetic relationships. The direction of horizontal transfer from Squamata to the ancestor of Ruminantia is evident from the genetic distances and discontinuous phylogenetic distribution of Bov-B LINE elements. The ancestor of Colubroidea snakes has been recognized as a possible donor of Bov-B LINE elements to Ruminantia. The timing of horizontal transfer has been estimated from the distribution of Bov-B LINE elements in Ruminantia and the fossil data of Ruminantia to be 40-50 My ago. The phylogenetic relationships of Bov-B LINE elements from the various Squamata species agrees with that of the species phylogeny, suggesting that Bov-B LINE elements have been stably maintained by vertical transmission since the origin of Squamata in the Mesozoic era.     

Horizontal transfer of non-LTR retrotransposons in vertebrates

Since their discovery in family Bovidae (bovids), Bov-B LINEs, believed to be order-specific SINEs, have been found in all ruminants and recently also in Viperidae snakes. The distribution and the evolutionary relationships of Bov-B LINEs provide an indication of their origin and evolutionary dynamics in different species. The evolutionary origin of Bov-B LINE elements has been shown unequivocally to be in Squamata (squamates). The horizontal transfer of Bov-B LINE elements in vertebrates has been confirmed by their discontinuous phylogenetic distribution in Squamata (Serpentes and two lizard infra-orders) as well as in Ruminantia, by the high level of nucleotide identity, and by their phylogenetic relationships. The direction of horizontal transfer from Squamata to the ancestor of Ruminantia is evident from the genetic distances and discontinuous phylogenetic distribution of Bov-B LINE elements. The ancestral snake lineage (Boidae) has been recognized as a possible donor of Bov-B LINE elements to Ruminantia. The timing of horizontal transfer has been estimated from the distribution of Bov-B LINE elements in Ruminantia and the fossil data of Ruminantia to be 40–50mya. The phylogenetic relationships of Bov-B LINE elements from the various Squamata species agrees with that of the species phylogeny, suggesting that Bov-B LINE elements have been stably maintained by vertical transmission since the origin of Squamata in the Mesozoic era. This revised version was published online in July 2006 with corrections to the Cover Date.     

Retrotransposon populations of <Emphasis Type="Italic">Vicia</Emphasis> species with varying genome size

The (non-LTR) LINE and Ty3-gypsy-type LTR retrotransposon populations of three Vicia species that differ in genome size (Vicia faba, Vicia melanops and Vicia sativa) have been characterised. In each species the LINE retrotransposons comprise a complex, very heterogeneous set of sequences, while the Ty3-gypsy elements are much more homogeneous. Copy numbers of all three retrotransposon groups (Ty1-copia, Ty3-gypsy and LINE) in these species have been estimated by random genomic sequencing and Southern hybridisation analysis. The Ty3-gypsy elements are extremely numerous in all species, accounting for 18–35% of their genomes. The Ty1-copia group elements are somewhat less abundant and LINE elements are present in still lower amounts. Collectively, 20–45% of the genomes of these three Vicia species are comprised of retrotransposons. These data show that the three retrotransposon groups have proliferated to different extents in members of the Vicia genus and high proliferation has been associated with homogenisation of the retrotransposon population.Electronic Supplementary Material Supplementary material is available for this article at .     

Specific organization of ribosomal DNA arrays in Squamata

A first report on structural organization of ribosomal DNA arrays in some members of the order Squamata is presented. The data obtained point to unusually small (for vertebrates) size of the rDNA repetitive unit (approximately, 10 to 15 kb) in the lizard species examined. Analysis of BAC library of Uta stansburiana (Iguania) showed that haploid genome of this lizard contained a single cluster, consisting of about ten rDNA repeats. Determination of the extent of rDNA unit repetition in some other representatives of the order Squamata, using the method of comparative real-time PCR, showed that the number of rDNA units varied from one or several dozens in Iguanina to several hundred repeats in Scincomorpha and Varonoidea. The results are discussed in terms of an ambiguous position of the infraorder Iguania on the evolutionary trees constructed based on morphological and molecular data.     

A mitogenomic study on the phylogenetic position of snakes

Phylogenetic relationships of squamates (lizards, amphisbaenians and snakes) have received considerable attention, although no consensus has been reached concerning some basal divergences. This paper focuses on the Serpentes (snakes), whose phylogenetic position within the Squamata remains uncertain despite a number of morphological and molecular studies. Some mitogenomic studies have suggested a sister-group relationship between snakes and varanid lizards, while other studies have identified snakes and lizards as sister groups. However, recent studies using nuclear data have presented a different scenario, with snakes being more closely related to anguimorph and iguanian lizards. In this mitogenomic study we have examined the above hypotheses with the inclusion of amphisbaenians, one gekkotan and one acrodont lizard, taxa not represented in previous mitogenomic studies. To this end we have also extended the representation of snakes by sequencing five additional snake genomes: two scolecophidians ( Ramphotyphlops australis and Typhlops mirus ) two henophidians ( Eunectes notaeus and Boa constrictor ) and one caenophidian ( Elaphe guttata ). The phylogenetic analysis recovered snakes and amphisbaenians as sister groups, thereby differing from previous hypotheses. In addition to a discussion on previous morphological and molecular studies in light of the results presented here, the current study also provides some details regarding features of the new snake mitochondrial genomes described.     

Identification and characterization of expressed retrotransposons in the genome of the Paracoccidioides species complex

Background

Species from the Paracoccidioides complex are thermally dimorphic fungi and the causative agents of paracoccidioidomycosis, a deep fungal infection that is the most prevalent systemic mycosis in Latin America and represents the most important cause of death in immunocompetent individuals with systemic mycosis in Brazil. We previously described the identification of eight new families of DNA transposons in Paracoccidioides genomes. In this work, we aimed to identify potentially active retrotransposons in Paracoccidioides genomes.

Results

We identified five different retrotransposon families (four LTR-like and one LINE-like element) in the genomes of three Paracoccidioides isolates. Retrotransposons were present in all of the genomes analyzed. P. brasiliensis and P. lutzii species harbored the same retrotransposon lineages but differed in their copy numbers. In the Pb01, Pb03 and Pb18 genomes, the number of LTR retrotransposons was higher than the number of LINE-like elements, and the LINE-like element RtPc5 was transcribed in Paracoccidioides lutzii (Pb01) but could not be detected in P. brasiliensis (Pb03 and Pb18) by semi-quantitative RT-PCR.

Conclusion

Five new potentially active retrotransposons have been identified in the genomic assemblies of the Paracoccidioides species complex using a combined computational and experimental approach. The distribution across the two known species, P. brasiliensis and P. lutzii, and phylogenetics analysis indicate that these elements could have been acquired before speciation occurred. The presence of active retrotransposons in the genome may have implications regarding the evolution and genetic diversification of the Paracoccidioides genus.

Electronic supplementary material

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

Functional and Structural Divergence of an Unusual LTR Retrotransposon Family in Plants

Retrotransposons with long terminal repeats (LTRs) more than 3 kb are not frequent in most eukaryotic genomes. Rice LTR retrotransposon, Retrosat2, has LTRs greater than 3.2 kb and two open reading frames (ORF): ORF1 encodes enzymes for retrotransposition whereas no function can be assigned to ORF0 as it is not found in any other organism. A variety of experimental and in silico approaches were used to determine the origin of Retrosat2 and putative function of ORF0. Our data show that not only is Retrosat2 highly abundant in the Oryza genus, it may yet be active in rice. Homologs of Retrosat2 were identified in maize, sorghum, Arabidopsis and other plant genomes suggesting that the Retrosat2 family is of ancient origin. Several putatively cis-acting elements, some multicopy, that regulate retrotransposon replication or responsiveness to environmental factors were found in the LTRs of Retrosat2. Unlike the ORF1, the ORF0 sequences from Retrosat2 and homologs are divergent at the sequence level, 3D-structures and predicted biological functions. In contrast to other retrotransposon families, Retrosat2 and its homologs are dispersed throughout genomes and not concentrated in the specific chromosomal regions, such as centromeres. The genomic distribution of Retrosat2 homologs varies across species which likely reflects the differing evolutionary trajectories of this retrotransposon family across diverse species.     

Retrotransposon insertion polymorphisms in six rice genes and their evolutionary history

Retrotransposons are abundant in higher plant genomes. Although retrotransposons associated with plant genes have been identified, little is known about their evolutionary conservation at the level of species and subspecies. In the present study, we investigated the phylogenetic distribution of long terminal repeat (LTR) retrotransposon, long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) insertions in six genes in 95 cultivated and wild rice genotypes. These six genes are likely to be functional based on nonsynonymous (Ka) to synonymous (Ks) substitution ratios which were found to be significantly <1. Different conservation patterns of these retrotransposons in genes were observed in cultivated and wild rice species. Four out of seven retrotransposon insertions appear to predate the ancestral Oryza AA genome. Two of these insertions in genes 4 and 5 occurred early in the evolutionary history of Oryza. Two retrotransposon insertions in gene 1 arose after the divergence of Asian cultivated rice from its wild ancestor. Furthermore, the retrotransposon insertion in gene 3 appears to have occurred in the ancestral lineage leading to temperate japonicas. Conservation of retrotransposon insertions in genes in specific groups, species, and lineages might be related to their specific function.     

Patterns of genetic diversity in Hepatozoon spp. infecting snakes from North Africa and the Mediterranean Basin

Species of Hepatozoon Miller, 1908 are blood parasites most commonly found in snakes but some have been described from all tetrapod groups and a wide variety of hematophagous invertebrates. Previous studies have suggested possible associations between Hepatozoon spp. found in predators and prey. Particularly, some saurophagous snakes from North Africa and the Mediterranean region have been found to be infected with Hepatozoon spp. similar to those of various sympatric lizard hosts. In this study, we have screened tissue samples of 111 North African and Mediterranean snakes, using specific primers for the 18S rRNA gene. In the phylogenetic analysis, the newly-generated Hepatozoon spp. sequences grouped separately into five main clusters. Three of these clusters were composed by Hepatozoon spp. also found in snakes and other reptiles from the Mediterranean Basin and North Africa. In the other two clusters, the new sequences were not closely related to geographically proximate known sequences. The phylogeny of Hepatozoon spp. inferred here was not associated with intermediate host taxonomy or geographical distribution. From the other factors that could explain these evolutionary patterns, the most likely seems series of intermediate hosts providing similar ribotypes of Hepatozoon and a high prevalence of host shifts for Hepatozoon spp. This is indicated by ribotypes of high similarity found in different reptile families, as well as by divergent ribotypes found in the same host species. This potentially low host specificity has profound implications for the systematics of Hepatozoon spp.     

A family of short retroposons (Squam1) from squamate reptiles (Reptilia: Squamata): Structure,evolution, and correlation with phylogeny

Sequences of the SINE family specific to squamate reptiles have been isolated from the genomes of lacertid lizards and sequenced. These retroposons, which we called Squam1, are 360–390 bp long and contain a region similar to the tRNA gene sequence at the 5’ end. This family has also been detected in representatives of other reptile families (varanids, iguanids (Anolis), gekkonids, and snakes), being absent from the genomes of crocodiles as well as amphibians, birds, and mammals. The primary structures of Squam1 copies have been comprehensively analyzed and compared with GenBank sequences. The genomes of most taxa contain two to three SINE subfamilies with specific diagnostic features in their primary structures. Individual similarity between the copies within each taxon is about 85%, with intrageneric similarity being only slightly higher. A comparison of consensus sequences between different lizard families has shown that Squam1 may be a convenient phylogenetic marker for this group of reptiles, having a number of both apomorphic and more or less pronounced synapomorphic features. By this criterion, snakes slightly differ from lizards but obviously belong to the same clade. However, they show no special affinity to varanids as the putative closest relatives of snakes, compared to other lizards.     

Precise excision of long terminal repeats of the gypsy (mdg4) retrotransposon of Drosophila melanogaster detected in Escherichia coli cells is explained by its integrase function

An Escherichia coli model system was developed to estimate the capacity of the integrase of the Drosophila melanogaster retrotransposon gypsy (mdg4) for precise excision of the long terminal repeat (LTR) and, hence, the entire gypsy. The gypsy retrotransposon was cloned in the form of a PCR fragment in the pBlue-Script II KS⁺ vector (pBSLTR), and the region of the second open reading frame (INT ORF2) of this element encoding integrase was cloned under the lacZ promoter in the pUC19 vector and then recloned in pACYC184 compatible with pBSLTR. The LTR was cloned in such a manner that its precise excision from the recombinant plasmid led to the restoration of the nucleotide sequence and the function of the lacZ gene; therefore, it was detected by the appearance of blue colonies on a medium containing X-gal upon IPTG induction. Upon IPTG induction of E. coli XL-1 Blue cells obtained by cotransformation with plasmids pACYCint and pBSLTR on an X-gal-containing medium, blue clones appeared with a frequency of 10^?4 to 10^?3, the frequency of spontaneously appearing blue colonies not exceeding 10^?9 to 10^?8. The presence of blue colonies indicated that that the integrase encoded by the INT ORF2 (pACYCint) fragment was active. After the expression of the integrase, it recognized and excised the gypsy LTR from pBSLTR, precisely restoring the nucleotide sequence and the function of the lacZ gene, which led to the expression of the β-galactosidase enzymatic activity. PCR analysis confirmed that the LTR was excised precisely. Thus, the resultant biplasmid model system allowed precise excisions of the gypsy LTR from the target site to be detected. Apparently, the gypsy integrase affected not only the LTR of this mobile element, but also the host genome nucleotide sequences. The system is likely to have detected only some of the events occurring in E. coli cells. Thus, the integrase of gypsy is actually capable of not only transposing this element by inserting DNA copies of the gypsy retrotransposon to chromosomes of Drosophila, but also excising them. gypsy is excised via a precise mechanism, with the original nucleotide sequence of the target site being completely restored. The obtained data demonstrate the existence of alternative ways of the transposition of retrotransposons and, possibly, retroviruses, including gypsy (mdg4).     

Genetic Diversity of Cryptosporidium spp. in Captive Reptiles

The genetic diversity of Cryptosporidium in reptiles was analyzed by PCR-restriction fragment length polymorphism and sequence analysis of the small subunit rRNA gene. A total of 123 samples were analyzed, of which 48 snake samples, 24 lizard samples, and 3 tortoise samples were positive for Cryptosporidium. Nine different types of Cryptosporidium were found, including Cryptosporidium serpentis, Cryptosporidium desert monitor genotype, Cryptosporidium muris, Cryptosporidium parvum bovine and mouse genotypes, one C. serpentis-like parasite in a lizard, two new Cryptosporidium spp. in snakes, and one new Cryptosporidium sp. in tortoises. C. serpentis and the desert monitor genotype were the most common parasites and were found in both snakes and lizards, whereas the C. muris and C. parvum parasites detected were probably the result of ingestion of infected rodents. Sequence and biologic characterizations indicated that the desert monitor genotype was Cryptosporidium saurophilum. Two host-adapted C. serpentis genotypes were found in snakes and lizards.     

Genomic organization of repetitive DNAs in the cichlid fish <Emphasis Type="Italic">Astronotus ocellatus</Emphasis>

To contribute to the knowledge of fish genomes, we identified and characterized by means of nucleotide sequencing and physical chromosome mapping, three classes of repetitive DNAs in the genome of the South American cichlid fish Astronotus ocellatus. The first class corresponds to a satellite DNA family (AoSat) that shares similarity with a centromeric satellite DNA of the pufferfish Tetraodon nigroviridis. The second repetitive DNA class (AoRex3) is related to the retrotransposon Rex3, which is widely distributed among teleost fishes. The last repetitive element (AoLINE) shows a high similarity to the CR1-like LINE element of other teleosts. The three isolated repetitive elements are clustered in the centromeric heterochromatin of all chromosomes of the complement. The repetitive sequences are not randomly distributed in the genome, suggesting a pattern of compartmentalization on chromosomes.     

Stress-induced rearrangement of Fusarium retrotransposon sequences

Rearrangement of Fusarium oxysporum retro- transposon skippy was induced by growth in the presence of potassium chlorate. Three fungal strains, one sensitive to chlorate (Co60) and two resistant to chlorate and deficient for nitrate reductase (Co65 and Co94), were studied by Southern analysis of their genomic DNA. Polymorphism was detected in their hybridization banding pattern, relative to the wild type grown in the absence of chlorate, using various enzymes with or without restriction sites within the retrotransposon. Results were consistent with the assumption that three different events had occurred in strain Co60: genomic amplification of skippy yielding tandem arrays of the element, generation of new skippy sequences, and deletion of skippy sequences. Amplification of Co60 genomic DNA using the polymerase chain reaction and divergent primers derived from the retrotransposon generated a new band, corresponding to one long terminal repeat plus flanking sequences, that was not present in the wild-type strain. Molecular analysis of nitrate reductase-deficient mutants showed that generation and deletion of skippy sequences, but not genomic amplification in tandem repeats, had occurred in their genomes.