Sequence organization of barley centromeres

By sequencing, fingerprinting and in situ hybridization of a centromere-specific large insert clone (BAC 7), the sequence organization of centromeric DNA of barley could be elucidated. Within 23 kb, three copies of the Ty3/gypsy-like retroelement cereba were present. Two elements of ~7 kb, arranged in tandem, include long terminal repeats (LTRs) (~1 kb) similar to the rice centromeric retrotransposon RIRE 7 and to the cereal centromeric sequence family, the primer binding site, the complete polygene flanked by untranslated regions, as well as a polypurine tract 5′ of the downstream LTR. The high density (~200 elements/centromere) and completeness of cereba elements and the absence of internally deleted elements and solo LTRs from the BAC 7 insert represent unique features of the barley centromeres as compared to those of other cereals. Obviously, the conserved cereba elements together with barley-specific G+C-rich satellite sequences constitute the major components of centromeric DNA in this species.     

The Cassandra retrotransposon landscape in sugar beet (Beta vulgaris) and related Amaranthaceae: recombination and re-shuffling lead to a high structural variability

Background and AimsPlant genomes contain many retrotransposons and their derivatives, which are subject to rapid sequence turnover. As non-autonomous retrotransposons do not encode any proteins, they experience reduced selective constraints leading to their diversification into multiple families, usually limited to a few closely related species. In contrast, the non-coding Cassandra terminal repeat retrotransposons in miniature (TRIMs) are widespread in many plants. Their hallmark is a conserved 5S rDNA-derived promoter in their long terminal repeats (LTRs). As sugar beet (Beta vulgaris) has a well-described LTR retrotransposon landscape, we aim to characterize TRIMs in beet and related genomes.MethodsWe identified Cassandra retrotransposons in the sugar beet reference genome and characterized their structural relationships. Genomic organization, chromosomal localization, and distribution of Cassandra-TRIMs across the Amaranthaceae were verified by Southern and fluorescent in situ hybridization.Key resultsAll 638 Cassandra sequences in the sugar beet genome contain conserved LTRs and thus constitute a single family. Nevertheless, variable internal regions required a subdivision into two Cassandra subfamilies within B. vulgaris. The related Chenopodium quinoa harbours a third subfamily. These subfamilies vary in their distribution within Amaranthaceae genomes, their insertion times and the degree of silencing by small RNAs. Cassandra retrotransposons gave rise to many structural variants, such as solo LTRs or tandemly arranged Cassandra retrotransposons. These Cassandra derivatives point to an interplay of template switch and recombination processes – mechanisms that likely caused Cassandra’s subfamily formation and diversification.ConclusionsWe traced the evolution of Cassandra in the Amaranthaceae and detected a considerable variability within the short internal regions, whereas the LTRs are strongly conserved in sequence and length. Presumably these hallmarks make Cassandra a prime target for unequal recombination, resulting in the observed structural diversity, an example of the impact of LTR-mediated evolutionary mechanisms on the host genome.     

Sobo, a recently amplified satellite repeat of potato, and its implications for the origin of tandemly repeated sequences

Highly repetitive satellite DNA sequences are main components of heterochromatin in higher eukaryotic genomes. It is well known that satellite repeats can expand and contract dramatically, which may result in significant genome size variation among genetically related species. The origin of satellite repeats, however, is elusive. Here we report a satellite repeat, Sobo, from a diploid potato species, Solanum bulbocastanum. The Sobo repeat is mapped to a single location in the pericentromeric region of chromosome 7. This single Sobo locus spans approximately 360 kb of a 4.7-kb monomer. Sequence analysis revealed that the major part of the Sobo monomer shares significant sequence similarity with the long terminal repeats (LTRs) of a retrotransposon. The Sobo repeat was not detected in other Solanum species and is absent in some S. bulbocastanum accessions. Sobo monomers are highly homogenized and share >99% sequence identity. These results suggest that the Sobo repeat is a recently emerged satellite and possibly originated by a sudden amplification of a genomic region including the LTR of a retrotransposon and its flanking genomic sequences.     

A lineage‐specific centromere retrotransposon in Oryza brachyantha

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3‐gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo‐LTRs to full‐length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3‐associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.     

skippy,a retrotransposon from the fungal plant pathogen Fusarium oxysporum

A retrotransposon from the fungal plant pathogen Fusarium oxysporum f. sp. lycopersici has been isolated and characterized. The element, designated skippy (skp) is 7846 by in length, flanked by identical long terminal repeats (LTR) of 429 by showing structural features characteristic of retroviral and retrotransposon LTRs. Target-site duplications of 5 bp were found. Two long overlapping open reading frames (ORF) were identified. The first ORF, 2562 by in length, shows homology to retroviral gag genes. The second ORF, 3888 bp in length, has homology to the protease, reverse transciptase. RNase H and integrase domains of retroelement pol genes in that order. Sequence comparisons and the order of the predicted proteins from skippy indicate that the element is closely related to the gypsy family of LTR-retrotransposons. The element is present in similar copy numbers in the two races investigated, although RFLP analysis showed differences in banding patterns. The number of LTR sequences present in the genome is higher than the number of copies of complete elements, indicating excision by homologous recombination between LTR sequences.     

The Tvv1 retrotransposon family is conserved between plant genomes separated by over 100 million years

Key message

Combining several different approaches, we have examined the structure, variability, and distribution of Tvv1 retrotransposons. Tvv1 is an unusual example of a low-copy retrotransposon metapopulation dispersed unevenly among very distant species and is promising for the development of molecular markers.

Abstract

Retrotransposons are ubiquitous throughout the genomes of the vascular plants, but individual retrotransposon families tend to be confined to the level of plant genus or at most family. This restricts the general applicability of a family as molecular markers. Here, we characterize a new plant retrotransposon named Tvv1_Sdem, a member of the Copia superfamily of LTR retrotransposons, from the genome of the wild potato Solanum demissum. Comparative analyses based on structure and sequence showed a high level of similarity of Tvv1_Sdem with Tvv1-VB, a retrotransposon previously described in the grapevine genome Vitis vinifera. Extending the analysis to other species by in silico and in vitro approaches revealed the presence of Tvv1 family members in potato, tomato, and poplar genomes, and led to the identification of full-length copies of Tvv1 in these species. We were also able to identify polymorphism in UTL sequences between Tvv1_Sdem copies from wild and cultivated potatoes that are useful as molecular markers. Combining different approaches, our results suggest that the Tvv1 family of retrotransposons has a monophyletic origin and has been maintained in both the rosids and the asterids, the major clades of dicotyledonous plants, since their divergence about 100 MYA. To our knowledge, Tvv1 represents an unusual plant retrotransposon metapopulation comprising highly similar members disjointedly dispersed among very distant species. The twin features of Tvv1 presence in evolutionarily distant genomes and the diversity of its UTL region in each species make it useful as a source of robust molecular markers for diversity studies and breeding.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

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

Characterization of a Novel Class of Interspersed LTR Elements in Primate Genomes: Structure, Genomic Distribution, and Evolution

Retrovirus-like sequences and their solitary (solo) long terminal repeats (LTRs) are common repetitive elements in eukaryotic genomes. We reported previously that the tandemly arrayed genes encoding U2 snRNA (the RNU2 locus) in humans and apes contain a solo LTR (U2-LTR) which was presumably generated by homologous recombination between the two LTRs of an ancestral provirus that is retained in the orthologous baboon RNU2 locus. We have now sequenced the orthologous U2-LTRs in human, chimpanzee, gorilla, orangutan, and baboon and examined numerous homologs of the U2-LTR that are dispersed throughout the human genome. Although these U2-LTR homologs have been collectively referred to as LTR13 in the literature, they do not display sequence similarity to any known retroviral LTRs; however, the structure of LTR13 closely resembles that of other retroviral LTRs with a putative promoter, polyadenylation signal, and a tandemly repeated 53-bp enhancer-like element. Genomic blotting indicates that LTR13 is primate-specific; based on sequence analysis, we estimate there are about 2,500 LTR13 elements in the human genome. Comparison of the primate U2-LTR sequences suggests that the homologous recombination event that gave rise to the solo U2-LTR occurred soon after insertion of the ancestral provirus into the ancestral U2 tandem array. Phylogenetic analysis of the LTR13 family confirms that it is diverse, but the orthologous U2-LTRs form a coherent group in which chimpanzee is closest to the humans; orangutan is a clear outgroup of human, chimpanzee, and gorilla; and baboon is a distant relative of human, chimpanzee, gorilla, and orangutan. We compare the LTR13 family with other known LTRs and consider whether these LTRs might play a role in concerted evolution of the primate RNU2 locus. Received: 29 September 1997 / Accepted: 16 January 1998     

Cloning and molecular characterization of a novel rolling-circle replicating plasmid, pK1S-1, from Bacillus thuringiensis subsp. kurstaki K1

Bacillus thuringiensis, an entomopathogenic bacterium belonging to the B. cereus group, harbors numerous extra-chromosomal DNA molecules whose sizes range from 2 to 250 kb. In this study, we used a plasmid capture system (PCS) to clone three small plasmids from B. thuringiensis subsp. kurstaki Kl which were not found in B. thuringiensis subsp. kurstaki HD-1, and determined the complete nucleotide sequence of plasmid pKlS-1 (5.5 kb). Of the six putative open reading frames (ORF2-ORF7) in pKlS-1, ORF2 (MobKl) showed approximately 90% aa identity with the Mob-proteins of pGI2 and pTX14-2, which are rolling circle replicating group VII (RCR group VII) plasmids from B. thuringiensis. In addition, a putative origin of transfer (oriT) showed 95.8% identity with those of pGI2 and pTX14-2. ORF3 (RepK1) showed relatively low aa identity (17.8～25.2%) with the Rep protein coded by RCR plasmids, however. The putative double-strand origin of replication (dso) and single-strand origin of replication (sso) of pKlS-1 exhibited approximately 70% and 64% identities with those of pGI2 and pTX14-2. ORF6 and 7 showed greater than 50% similarities with alkaline serine protease, which belongs to the subtilase family. The other 2 ORFs were identified as hypothetical proteins. To determine the replicon of pKlS-1, seven subclones were contructed in the B. thuringiensis ori-negative pHTIK vector and were electroporated into a plasmid cured B. thuringiensis strain. The 1.6 kb region that included the putative ORF3 (ReplK), dso and ORF4, exhibited replication ability. These findings identified pKlS-1 as a new RCR group VII plasmid, and determined its replication region.     

Retrotransposons and Tandem Repeat Sequences in the Nuclear Genomes of Cryptomonad Algae

The cryptomonads are an enigmatic group of unicellular eukaryotic algae that possess two nuclear genomes, having acquired photosynthesis by the uptake and retention of a eukaryotic algal endosymbiont. The endosymbiont nuclear genome, or nucleomorph, of the cryptomonad Guillardia theta has been completely sequenced: at only 551 kilobases (kb) and with a gene density of ∼1 gene/kb, it is a model of compaction. In contrast, very little is known about the structure and composition of the cryptomonad host nuclear genome. Here we present the results of two small-scale sequencing surveys of fosmid clone libraries from two distantly related cryptomonads, Rhodomonas salina CCMP1319 and Cryptomonas paramecium CCAP977/2A, corresponding to ∼150 and ∼235 kb of sequence, respectively. Very few of the random end sequences determined in this study show similarity to known genes in other eukaryotes, underscoring the considerable evolutionary distance between the cryptomonads and other eukaryotes whose nuclear genomes have been completely sequenced. Using a combination of fosmid clone end-sequencing, Southern hybridizations, and PCR, we demonstrate that Ty3-gypsy long-terminal repeat (LTR) retrotransposons and tandem repeat sequences are a prominent feature of the nuclear genomes of both organisms. The complete sequence of a 30.9-kb genomic fragment from R. salina was found to contain a full-length Ty3-gypsy element with near-identical LTRs and a chromodomain, a protein module suggested to mediate the site-specific integration of the retrotransposon. The discovery of chromodomain-containing retroelements in cryptomonads further expands the known distribution of the so-called chromoviruses across the tree of eukaryotes. [Reviewing Editor: Dr. Debashish Bhattacharya]     

Aurora,a non-mobile retrotransposon in Drosophila melanogaster heterochromatin

A novel retrotransposon, aurora, containing 324 by long terminal repeats (LTRs) was detected in Drosophila melanogaster as a 5 kb insertion in the heterochromatic Stellate gene. This insertion causes a 5 bp duplication of the integration site. Southern analysis and in situ hybridization data show that all detectable copies of aurora are immobilized in the D. melanogaster heterochromatin. However, mobile copies of aurora were revealed in the cuchromatin of D. simulans. The element was also found in various species of the melanogaster subgroup and in the D. virilis genome.     

Molecular markers based on LTR retrotransposons BARE-1 and Jeli uncover different strata of evolutionary relationships in diploid wheats

Molecular markers based on retrotransposon insertions are widely used for various applications including phylogenetic analysis. Multiple cases were described where retrotransposon-based markers, namely sequence-specific amplification polymorphism (SSAP), were superior to other marker types in resolving the phylogenetic relationships due to their higher variability and informativeness. However, the patterns of evolutionary relationships revealed by SSAP may be dependent on the underlying retrotransposon activity in different periods of time. Hence, the proper choice of retrotransposon family is essential for obtaining significant results. We compared the phylogenetic trees for a diverse set of diploid A-genome wheat species (Triticum boeoticum, T. urartu and T. monococcum) based on two unrelated retrotransposon families, BARE-1 and Jeli. BARE-1 belongs to Copia class and has a uniform distribution between common wheat (T. aestivum) genomes of different origin (A, B and D), indicating similar activity in the respective diploid genome donors. Gypsy-class family Jeli was found by us to be an A-genome retrotransposon with >70% copies residing in A genome of hexaploid common wheat, suggesting a burst of transposition in the history of A-genome progenitors. The results indicate that a higher Jeli transpositional activity was associated with T. urartu versus T. boeoticum speciation, while BARE-1 produced more polymorphic insertions during subsequent intraspecific diversification; as an outcome, each retrotransposon provides more informative markers at the corresponding level of phylogenetic relationships. We conclude that multiple retroelement families should be analyzed for an image of evolutionary relationships to be solid and comprehensive.     

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.     

How Hot Are Drosophila Hotspots? Examining Recombination Rate Variation and Associations with Nucleotide Diversity,Divergence, and Maternal Age in Drosophila pseudoobscura

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5kb regions of highly elevated recombination rates (10–100x the background rate) called “hotspots.” Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40Kb region of Drosophila pseudoobscura chromosome 2, with one 20kb interval assayed every 5Kb and the adjacent 20kb interval bisected into 10kb pieces. We found that recombination events across the 40kb stretch were relatively evenly distributed across each of the 5kb and 10kb intervals, rather than concentrated in a single 5kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.     

IRAP and REMAP: two new retrotransposon-based DNA fingerprinting techniques

 The BARE-1 retrotransposon is an active, dispersed, and highly abundant component of the genome of barley (Hordeum vulgare) and other species in its genus. Like all retrotransposons of its kind, BARE-1 is bounded by long terminal repeats (LTRs). We have developed two amplification-based marker methods based on the position of given LTRs within the genome. The IRAP (Inter-Retrotransposon Amplified Polymorphism) markers are generated by the proximity of two LTRs using outward-facing primers annealing to LTR target sequences. In REMAP (REtrotransposon-Microsatellite Amplified Polymorphism), amplification between LTRs proximal to simple sequence repeats such as constitute microsatellites produces markers. The methods can distinguish between barley varieties and produce fingerprint patterns for species across the genus. The patterns indicate that although the BARE-1 family of retrotransposons is disperse, these elements are locally clustered or nested and often found near tandem arrays of a simple sequence repeat. Received: 30 June 1998 / Accepted: 21 August 1998     

A family of developmentally excised DNA elements in Tetrahymena is under selective pressure to maintain an open reading frame encoding an integrase-like protein

Tlr1 is a member of a family of ~20–30 DNA elements that undergo developmentally regulated excision during formation of the macronucleus in the ciliated protozoan Tetrahymena. Analysis of sequence internal to the right boundary of Tlr1 revealed the presence of a 2 kb open reading frame (ORF) encoding a deduced protein with similarity to retrotransposon integrases. The ORFs of five unique clones were sequenced. The ORFs have 98% sequence conservation and align without frameshifts, although one has an additional trinucleotide at codon 561. Nucleotide changes among the five clones are highly non-random with respect to the position in the codon and 93% of the nucleotide changes among the five clones encode identical or similar amino acids, suggesting that the ORF has evolved under selective pressure to preserve a functional protein. Nineteen T/C transitions in T/CAA and T/CAG codons suggest selection has occurred in the context of the Tetrahymena genome, where TAA and TAG encode Gln. Similarities between the ORF and those encoding retrotransposon integrases suggest that the Tlr family of elements may encode a polynucleotide transferase. Possible roles for the protein in transposition of the elements within the micronuclear genome and/or their developmentally regulated excision from the macronucleus are discussed.     

The Mitochondrial Genome of Raphanus sativus and Gene Evolution of Cruciferous Mitochondrial Types

To explore the mitochondrial genes of the Cruciferae family, the mitochondrial genome of Raphanus sativus (sat) was sequenced and annotated. The circular mitochondrial genome of sat is 239,723 bp and includes 33 protein-coding genes, three rRNA genes and 17 tRNA genes. The mitochondrial genome also contains a pair of large repeat sequences 5.9 kb in length, which may mediate genome reorga-nization into two sub-genomic circles, with predicted sizes of 124.8 kb and 115.0 kb, respectively. Furthermore, gene evolution of mitochondrial genomes within the Cruciferae family was analyzed using sat mitochondrial type (mitotype), together with six other re-ported mitotypes. The cruciferous mitochondrial genomes have maintained almost the same set of functional genes. Compared with Cycas taitungensis (a representative gymnosperm), the mitochondrial genomes of the Cruciferae have lost nine protein-coding genes and seven mitochondrial-like tRNA genes, but acquired six chloroplast-like tRNAs. Among the Cruciferae, to maintain the same set of genes that are necessary for mitochondrial function, the exons of the genes have changed at the lowest rates, as indicated by the numbers of single nucleotide polymorphisms. The open reading frames (ORFs) of unknown function in the cruciferous genomes are not conserved. Evolutionary events, such as mutations, genome reorganizations and sequence insertions or deletions (indels), have resulted in the non- conserved ORFs in the cruciferous mitochondrial genomes, which is becoming significantly different among mitotypes. This work represents the first phylogenic explanation of the evolution of genes of known function in the Cruciferae family. It revealed significant variation in ORFs and the causes of such variation.