Structural and sequence diversity of the transposon Galileo in the Drosophila willistoni genome

Background

Galileo is one of three members of the P superfamily of DNA transposons. It was originally discovered in Drosophila buzzatii, in which three segregating chromosomal inversions were shown to have been generated by ectopic recombination between Galileo copies. Subsequently, Galileo was identified in six of 12 sequenced Drosophila genomes, indicating its widespread distribution within this genus. Galileo is strikingly abundant in Drosophila willistoni, a neotropical species that is highly polymorphic for chromosomal inversions, suggesting a role for this transposon in the evolution of its genome.

Results

We carried out a detailed characterization of all Galileo copies present in the D. willistoni genome. A total of 191 copies, including 133 with two terminal inverted repeats (TIRs), were classified according to structure in six groups. The TIRs exhibited remarkable variation in their length and structure compared to the most complete copy. Three copies showed extended TIRs due to internal tandem repeats, the insertion of other transposable elements (TEs), or the incorporation of non-TIR sequences into the TIRs. Phylogenetic analyses of the transposase (TPase)-encoding and TIR segments yielded two divergent clades, which we termed Galileo subfamilies V and W. Target-site duplications (TSDs) in D. willistoni Galileo copies were 7- or 8-bp in length, with the consensus sequence GTATTAC. Analysis of the region around the TSDs revealed a target site motif (TSM) with a 15-bp palindrome that may give rise to a stem-loop secondary structure.

Conclusions

There is a remarkable abundance and diversity of Galileo copies in the D. willistoni genome, although no functional copies were found. The TIRs in particular have a dynamic structure and extend in different ways, but their ends (required for transposition) are more conserved than the rest of the element. The D. willistoni genome harbors two Galileo subfamilies (V and W) that diverged ~9 million years ago and may have descended from an ancestral element in the genome. Galileo shows a significant insertion preference for a 15-bp palindromic TSM.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-792) contains supplementary material, which is available to authorized users.     

Cloning and sequencing of the breakpoint regions of inversion <Emphasis Type="Italic">5g</Emphasis> fixed in <Emphasis Type="Italic">Drosophila buzzatii</Emphasis>

Chromosomal inversions are ubiquitous in Drosophila both as intraspecific polymorphisms and interspecific differences. Many gaps still remain in our understanding of the mechanisms that generate them. Previous work has shown that in Drosophila buzzatii, three polymorphic inversions were generated by ectopic recombination between copies of the transposon Galileo. In this study, we have characterized the breakpoint regions of inversion 5g, fixed in D. buzzatii and absent in Drosophila koepferae and other closely related species. A novel approach comprising four experimental steps was used. First, D. buzzatii BAC clones encompassing the breakpoints were identified and their ends sequenced. Then, breakpoint regions were mapped at high resolution in the Drosophila mojavensis genome sequence. Finally, breakpoint regions were isolated by polymerase chain reaction in D. buzzatii and D. koepferae and sequenced. Our aim was to shed light on the mechanism that generated inversion 5g and specifically to test for an implication of the transposon Galileo. No evidence implicates Galileo or other transposable elements in the origin of inversion 5g that was generated most likely by two independent breaks and non-homologous end-joining repair. Our results show that different inversion-generating mechanisms may coexist within the same lineage and suggest a hypothesis for the evolutionary time and mode of their operation. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Josefa González and Alfredo Ruiz contributed equally.     

Abundance and chromosomal distribution of six Drosophila buzzatii transposons: BuT1, BuT2, BuT3, BuT4, BuT5, and BuT6

The abundance and chromosomal distribution of six class-II transposable elements (TEs) of Drosophila buzzatii have been analyzed by Southern blotting and in situ hybridization. These six transposons had been previously found at the breakpoints of inversions 2j and 2q ⁷ of D. buzzatii. These two polymorphic inversions were generated by an ectopic recombination event between two copies of Galileo, a Foldback element. The four breakpoints became hotspots for TE insertions after the generation of the inversion and the transposons analyzed in this work are considered to be secondary invaders of these regions. Insertions of the six transposons are present in the euchromatin but show an increased density in the pericentromeric euchromatin–heterochromatin transition region and the dot chromosome. They are also more abundant in the inverted segments of chromosome 2 rearrangements. We further observed that the accumulation of TE insertions varies between elements and is correlated between dot, proximal regions, and inverted segments. These observations fully agree with previous data in Drosophila melanogaster and support recombination rate as the chief force explaining the chromosomal distribution of TEs.Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.Sequence data from this article have been deposited in the EMBL/GenBank Data Libraries under accession number DQ402469.     

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.     

INVERSION LENGTH AND BREAKPOINT DISTRIBUTION IN THE DROSOPHILA BUZZATII SPECIES COMPLEX: IS INVERSION LENGTH A SELECTED TRAIT?

Length and position of breakpoints are characteristics of inversions that can be precisely determined on the polytene chromosomes of Drosophila species, and they provide crucial information about the processes that govern the origin and evolution of inversions. Eighty-six paracentric inversions described in the Drosophila buzzatii species complex and 18 inversions induced by introgressive hybridization in D. buzzatii were analyzed. In contrast to previous studies, inversion length and breakpoint distribution have been considered simultaneously. We conclude that: (1) inversion length is a selected trait; rare inversions are predominantly small while evolutionarily successful inversions, polymorphic and fixed, are predominantly intermediate in length; a nearly continuous variation in length, from small to medium sized, is found between less and more successful inversions; (2) there exists a significant negative correlation between length and number of polymorphic inversions per species which explains 39% of the inversion length variance; (3) natural selection on inversion length seems the main factor determining the relative position of breakpoints along the chromosomes; (4) the distribution of breakpoints according to their band location is non-random, with chromosomal segments that accumulate up to eight breakpoints.     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Stage-Specific Effects of Candidate Heterochronic Genes on Variation in Developmental Time along an Altitudinal Cline of Drosophila melanogaster

Background

Previously, we have shown there is clinal variation for egg-to-adult developmental time along geographic gradients in Drosophila melanogaster. Further, we also have identified mutations in genes involved in metabolic and neurogenic pathways that affect development time (heterochronic genes). However, we do not know whether these loci affect variation in developmental time in natural populations.

Methodology/Principal Findings

Here, we constructed second chromosome substitution lines from natural populations of Drosophila melanogaster from an altitudinal cline, and measured egg-adult development time for each line. We found not only a large amount of genetic variation for developmental time, but also positive associations of the development time with thermal amplitude and altitude. We performed genetic complementation tests using substitution lines with the longest and shortest developmental times and heterochronic mutations. We identified segregating variation for neurogenic and metabolic genes that largely affected the duration of the larval stages but had no impact on the timing of metamorphosis.

Conclusions/Significance

Altitudinal clinal variation in developmental time for natural chromosome substitution lines provides a unique opportunity to dissect the response of heterochronic genes to environmental gradients. Ontogenetic stage-specific variation in invected, mastermind, cricklet and CG14591 may affect natural variation in development time and thermal evolution.     

Atmospheric Hypoxia Limits Selection for Large Body Size in Insects

Background

The correlations between Phanerozoic atmospheric oxygen fluctuations and insect body size suggest that higher oxygen levels facilitate the evolution of larger size in insects.

Methods and Principal Findings

Testing this hypothesis we selected Drosophila melanogaster for large size in three oxygen atmospheric partial pressures (aPO₂). Fly body sizes increased by 15% during 11 generations of size selection in 21 and 40 kPa aPO₂. However, in 10 kPa aPO₂, sizes were strongly reduced. Beginning at the 12^th generation, flies were returned to normoxia. All flies had similar, enlarged sizes relative to the starting populations, demonstrating that selection for large size had functionally equivalent genetic effects on size that were independent of aPO₂.

Significance

Hypoxia provided a physical constraint on body size even in a tiny insect strongly selected for larger mass, supporting the hypothesis that Triassic hypoxia may have contributed to a reduction in insect size.     

Intraspecific Inversions Pose a Challenge for the trnH-psbA Plant DNA Barcode

Background

The chloroplast trnH-psbA spacer region has been proposed as a prime candidate for use in DNA barcoding of plants because of its high substitution rate. However, frequent inversions associated with palindromic sequences within this region have been found in multiple lineages of Angiosperms and may complicate its use as a barcode, especially if they occur within species.

Methodology/Principal Findings

Here, we evaluate the implications of intraspecific inversions in the trnH-psbA region for DNA barcoding efforts. We report polymorphic inversions within six species of Gentianaceae, all narrowly circumscribed morphologically: Gentiana algida, Gentiana fremontii, Gentianopsis crinita, Gentianopsis thermalis, Gentianopsis macrantha and Frasera speciosa. We analyze these sequences together with those from 15 other species of Gentianaceae and show that typical simple methods of sequence alignment can lead to misassignment of conspecifics and incorrect assessment of relationships.

Conclusions/Significance

Frequent inversions in the trnH-psbA region, if not recognized and aligned appropriately, may lead to large overestimates of the number of substitution events separating closely related lineages and to uniting more distantly related taxa that share the same form of the inversion. Thus, alignment of the trnH-psbA spacer region will need careful attention if it is used as a marker for DNA barcoding.     

Comparative Analysis of the Recently Discovered hAT Transposon TcBuster in Human Cells

Background

Transposons are useful tools for creating transgenic organisms, insertional mutagenesis, and genome engineering. TcBuster, a novel hAT-family transposon system derived from the red flour beetle Tribolium castaneum, was shown to be highly active in previous studies in insect embryoes.

Methodology/Principal Findings

We tested TcBuster for its activity in human embryonic kidney 293 (HEK-293) cells. Excision footprints obtained from HEK-293 cells contained small insertions and deletions consistent with a hAT-type repair mechanism of hairpin formation and non-homologous end-joining. Genome-wide analysis of 23,417 piggyBac, 30,303 Sleeping Beauty, and 27,985 TcBuster integrations in HEK-293 cells revealed a uniquely different integration pattern when compared to other transposon systems with regards to genomic elements. TcBuster experimental conditions were optimized to assay TcBuster activity in HEK-293 cells by colony assay selection for a neomycin-containing transposon. Increasing transposon plasmid increased the number of colonies, whereas gene transfer activity dependent on codon-optimized transposase plasmid peaked at 100 ng with decreased colonies at the highest doses of transposase DNA. Expression of the related human proteins Buster1, Buster3, and SCAND3 in HEK-293 cells did not result in genomic integration of the TcBuster transposon. TcBuster, Tol2, and piggyBac were compared directly at different ratios of transposon to transposase and found to be approximately comparable while having their own ratio preferences.

Conclusions/Significance

TcBuster was found to be highly active in mammalian HEK-293 cells and represents a promising tool for mammalian genome engineering.     

Characterization of Transposable Elements in the Ectomycorrhizal Fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.     

PacBio-LITS: a large-insert targeted sequencing method for characterization of human disease-associated chromosomal structural variations

Background

Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high.

Results

We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki–Lupski syndrome (PTLS; MIM#610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants.

Conclusions

The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology.

Electronic supplementary material

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

Horizontal transfer of a non-autonomous Helitron among insect and viral genomes

Background

The movement of mobile elements among species by horizontal transposon transfer (HTT) influences the evolution of genomes through the modification of structure and function. Helitrons are a relatively new lineage of DNA-based (class II) transposable elements (TEs) that propagate by rolling-circle replication, and are capable of acquiring host DNA. The rapid spread of Helitrons among animal lineages by HTT is facilitated by shuttling in viral particles or by unknown mechanisms mediated by close organism associations (e.g. between hosts and parasites).

Results

A non-autonomous Helitron independently annotated as BmHel-2 from Bombyx mori and the MITE01 element from Ostrinia nubilalis was predicted in the genomes of 24 species in the insect Order Lepidoptera. Integrated Helitrons retained ≥ 65% sequence identity over a 250 bp consensus, and were predicted to retain secondary structures inclusive of a 3′-hairpin and a 5′-subterminal inverted repeat. Highly similar Hel-2 copies were predicted in the genomes of insects and associated viruses, which along with a previous documented case of real-time virus-insect cell line transposition suggests that this Helitron has likely propagated by HTT.

Conclusions

These findings provide evidence that insect virus may mediate the HTT of Helitron-like TEs. This movement may facilitate the shuttling of DNA elements among insect genomes. Further sampling is required to determine the putative role of HTT in insect genome evolution.

Electronic supplementary material

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

RNAi Experiments in D. melanogaster: Solutions to the Overlooked Problem of Off-Targets Shared by Independent dsRNAs

Background

RNAi technology is widely used to downregulate specific gene products. Investigating the phenotype induced by downregulation of gene products provides essential information about the function of the specific gene of interest. When RNAi is applied in Drosophila melanogaster or Caenorhabditis elegans, often large dsRNAs are used. One of the drawbacks of RNAi technology is that unwanted gene products with sequence similarity to the gene of interest can be down regulated too. To verify the outcome of an RNAi experiment and to avoid these unwanted off-target effects, an additional non-overlapping dsRNA can be used to down-regulate the same gene. However it has never been tested whether this approach is sufficient to reduce the risk of off-targets.

Methodology

We created a novel tool to analyse the occurance of off-target effects in Drosophila and we analyzed 99 randomly chosen genes.

Principal Findings

Here we show that nearly all genes contain non-overlapping internal sequences that do show overlap in a common off-target gene.

Conclusion

Based on our in silico findings, off-target effects should not be ignored and our presented on-line tool enables the identification of two RNA interference constructs, free of overlapping off-targets, from any gene of interest.     

Identification of multiple binding sites for the THAP domain of the Galileo transposase in the long terminal inverted-repeats

Galileo is a DNA transposon responsible for the generation of several chromosomal inversions in Drosophila. In contrast to other members of the P-element superfamily, it has unusually long terminal inverted-repeats (TIRs) that resemble those of Foldback elements. To investigate the function of the long TIRs we derived consensus and ancestral sequences for the Galileo transposase in three species of Drosophilids. Following gene synthesis, we expressed and purified their constituent THAP domains and tested their binding activity towards the respective Galileo TIRs. DNase I footprinting located the most proximal DNA binding site about 70 bp from the transposon end. Using this sequence we identified further binding sites in the tandem repeats that are found within the long TIRs. This suggests that the synaptic complex between Galileo ends may be a complicated structure containing higher-order multimers of the transposase. We also attempted to reconstitute Galileo transposition in Drosophila embryos but no events were detected. Thus, although the limited numbers of Galileo copies in each genome were sufficient to provide functional consensus sequences for the THAP domains, they do not specify a fully active transposase. Since the THAP recognition sequence is short, and will occur many times in a large genome, it seems likely that the multiple binding sites within the long, internally repetitive, TIRs of Galileo and other Foldback-like elements may provide the transposase with its binding specificity.     

Comparative polytene chromosome maps of <Emphasis Type="Italic">D. montana</Emphasis> and <Emphasis Type="Italic">D. virilis</Emphasis>

Chromosomal inversion polymorphism was characterized in Finnish Drosophila montana populations. A total of 14 polymorphic inversions were observed in Finnish D. montana of which nine had not been described before. The number of polymorphic inversions in each chromosome was not significantly different from that expected, assuming equal chance of occurrence in the euchromatic genome. There was, however, no correlation between the number of polymorphic inversions and that of fixed inversions in each chromosome. Therefore, a simple neutral model does not explain the evolutionary dynamics of inversions. Furthermore, in contrast to results obtained by others, no significant correlation was found between the two transposable elements (TEs) Penelope and Ulysses and inversion breakpoints in D. montana. This result suggests that these TEs were not involved in the creation of the polymorphic inversions seen in D. montana. A comparative analysis of D. montana and Drosophila virilis polytene chromosomes 4 and 5 was performed with D. virilis bacteriophage P1 clones, thus completing the comparative studies of the two species.     

Complete genome sequence and analysis of Lactobacillus hokkaidonensis LOOC260T,a psychrotrophic lactic acid bacterium isolated from silage

Background

Lactobacillus hokkaidonensis is an obligate heterofermentative lactic acid bacterium, which is isolated from Timothy grass silage in Hokkaido, a subarctic region of Japan. This bacterium is expected to be useful as a silage starter culture in cold regions because of its remarkable psychrotolerance; it can grow at temperatures as low as 4°C. To elucidate its genetic background, particularly in relation to the source of psychrotolerance, we constructed the complete genome sequence of L. hokkaidonensis LOOC260^T using PacBio single-molecule real-time sequencing technology.

Results

The genome of LOOC260^T comprises one circular chromosome (2.28 Mbp) and two circular plasmids: pLOOC260-1 (81.6 kbp) and pLOOC260-2 (41.0 kbp). We identified diverse mobile genetic elements, such as prophages, integrated and conjugative elements, and conjugative plasmids, which may reflect adaptation to plant-associated niches. Comparative genome analysis also detected unique genomic features, such as genes involved in pentose assimilation and NADPH generation.

Conclusions

This is the first complete genome in the L. vaccinostercus group, which is poorly characterized, so the genomic information obtained in this study provides insight into the genetics and evolution of this group. We also found several factors that may contribute to the ability of L. hokkaidonensis to grow at cold temperatures. The results of this study will facilitate further investigation for the cold-tolerance mechanism of L. hokkaidonensis.

Electronic supplementary material

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