Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus

Background

The mariner family of transposable elements is one of the most widespread in the Metazoa. It is subdivided into several subfamilies that do not mirror the phylogeny of these species, suggesting an ancient diversification. Previous hybridization and PCR studies allowed a partial survey of mariner diversity in the Metazoa. In this work, we used a comparative genomics approach to access the genus-wide diversity and evolution of mariner transposable elements in twenty Drosophila sequenced genomes.

Results

We identified 36 different mariner lineages belonging to six distinct subfamilies, including a subfamily not described previously. Wide variation in lineage abundance and copy number were observed among species and among mariner lineages, suggesting continuous turn-over. Most mariner lineages are inactive and contain a high proportion of damaged copies. We showed that, in addition to substitutions that rapidly inactivate copies, internal deletion is a major mechanism contributing to element decay and the generation of non-autonomous sublineages. Hence, 23% of copies correspond to several Miniature Inverted-repeat Transposable Elements (MITE) sublineages, the first ever described in Drosophila for mariner. In the most successful MITEs, internal deletion is often associated with internal rearrangement, which sheds light on the process of MITE origin. The estimation of the transposition rates over time revealed that all lineages followed a similar progression consisting of a rapid amplification burst followed by a rapid decrease in transposition. We detected some instances of multiple or ongoing transposition bursts. Different amplification times were observed for mariner lineages shared by different species, a finding best explained by either horizontal transmission or a reactivation process. Different lineages within one species have also amplified at different times, corresponding to successive invasions. Finally, we detected a preference for insertion into short TA-rich regions, which appears to be specific to some subfamilies.

Conclusions

This analysis is the first comprehensive survey of this family of transposable elements at a genus scale. It provides precise measures of the different evolutionary processes that were hypothesized previously for this family based on PCR data analysis. mariner lineages were observed at almost all “life cycle” stages: recent amplification, subsequent decay and potential (re)-invasion or invasion of genomes.

Electronic supplementary material

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

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

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

Whole genome comparative analysis of transposable elements provides new insight into mechanisms of their inactivation in fungal genomes

Background

Transposable Elements (TEs) are key components that shape the organization and evolution of genomes. Fungi have developed defense mechanisms against TE invasion such as RIP (Repeat-Induced Point mutation), MIP (Methylation Induced Premeiotically) and Quelling (RNA interference). RIP inactivates repeated sequences by promoting Cytosine to Thymine mutations, whereas MIP only methylates TEs at C residues. Both mechanisms require specific cytosine DNA Methyltransferases (RID1/Masc1) of the Dnmt1 superfamily.

Results

We annotated TE sequences from 10 fungal genomes with different TE content (1-70%). We then used these TE sequences to carry out a genome-wide analysis of C to T mutations biases. Genomes from either Ascomycota or Basidiomycota that were massively invaded by TEs (Blumeria, Melampsora, Puccinia) were characterized by a low frequency of C to T mutation bias (10-20%), whereas other genomes displayed intermediate to high frequencies (25-75%). We identified several dinucleotide signatures at these C to T mutation sites (CpA, CpT, and CpG). Phylogenomic analysis of fungal Dnmt1 MTases revealed a previously unreported association between these dinucleotide signatures and the presence/absence of sub-classes of Dnmt1.

Conclusions

We identified fungal genomes containing large numbers of TEs with many C to T mutations associated with species-specific dinucleotide signatures. This bias suggests that a basic defense mechanism against TE invasion similar to RIP is widespread in fungi, although the efficiency and specificity of this mechanism differs between species. Our analysis revealed that dinucleotide signatures are associated with the presence/absence of specific Dnmt1 subfamilies. In particular, an RID1-dependent RIP mechanism was found only in Ascomycota.

Electronic supplementary material

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

RiTE database: a resource database for genus-wide rice genomics and evolutionary biology

Background

Comparative evolutionary analysis of whole genomes requires not only accurate annotation of gene space, but also proper annotation of the repetitive fraction which is often the largest component of most if not all genomes larger than 50 kb in size.

Results

Here we present the Rice TE database (RiTE-db) - a genus-wide collection of transposable elements and repeated sequences across 11 diploid species of the genus Oryza and the closely-related out-group Leersia perrieri. The database consists of more than 170,000 entries divided into three main types: (i) a classified and curated set of publicly-available repeated sequences, (ii) a set of consensus assemblies of highly-repetitive sequences obtained from genome sequencing surveys of 12 species; and (iii) a set of full-length TEs, identified and extracted from 12 whole genome assemblies.

Conclusions

This is the first report of a repeat dataset that spans the majority of repeat variability within an entire genus, and one that includes complete elements as well as unassembled repeats. The database allows sequence browsing, downloading, and similarity searches. Because of the strategy adopted, the RiTE-db opens a new path to unprecedented direct comparative studies that span the entire nuclear repeat content of 15 million years of Oryza diversity.

Electronic supplementary material

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

Genome ecosystem and transposable elements species

Transposable elements are known to be “selfish DNA” sequences able to spread and be maintained in all genomes analyzed so far. Their evolution depends on the interaction they have with the other components of the genome, including genes and other transposable elements. These relationships are complex and have often been compared to those of species living and competing in an ecosystem. The aim of this current work is a proposition to fill the conceptual gap existing between genome biology and ecology, assuming that genomic components, such as transposable elements families, can be compared to species interacting in an ecosystem. Using this framework, some of the main models defined in the population genetics of transposable elements can then been reformulated, and some new kinds of realistic relationships, such as symbiosis between different genomic components, can then be modelled and explored.     

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.     

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

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

Insertion sites of the transposable elementmariner are fixed in the genome ofDrosophila sechellia

Summary The abundance of the transposable elementmariner differs dramatically in the genomes of the closely related speciesDrosophila simulans, D. mauritiana, D. sechellia, andD. melanogaster. Natural populations ofD. simulans andD. mauritiana have 1–10 and 20–30 copies per diploid genome, respectively, and the insertion sites are polymorphic. The element has not been found inD. melanogaster. In this paper we show thatD. sechellia, a species endemic to the Seychelles Islands, contains only twomariner elements that are at fixed sites in the genome. One element, inserted in chromosome 2R at 51A1–2, contains three deletions and is missing much of the 3 end. The other element, inserted in chromosome 3L at 64A10–11, is the full length of 1286 bp. Although the sequence of the full-length element is polymorphic in populations ofD. sechellia, at least some of the sequences are closely related to elements fromD. simulans andD. mauritiana that are known to be active. However, judging from the progeny of crosses betweenD. sechellia andD. simulans, the biological activity of the full-lengthD. sechellia element appears to be low, either because of the nucleotide sequence of the element or because of its position in the genome, or both.     

Chromosome labeling with transposable elements in maize

Transposable elements randomly insert into a targeted locus at a frequency of 10^-6 to 10^-5. The En element has been shown in previous studies to transpose more frequently into closely linked sites. Thus, it is appropriate to place an En element onto each of the 20 chromosome arms in maize to maximize tagging efficiency. This is called chromosome labeling for tagging purposes with transposons. After a chromosome arm has been labeled with a transposon, genes residing in that arm will have a greater chance to be tagged by the transposon. To date, all of the maize chromosome arms have been labeled with at least one of five Encontaining alleles. The elements were linked to the arms using reciprocal translocations. The usage of these arm-labeled lines is discussed in the context of gene tagging.Journal Paper No. 15224 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa; Project No. 3176     

Mu transposable elements are structurally diverse and distributed throughout the genusZea

Summary The Robertson's Mutator stock of maize exhibits a high mutation rate due to the transposition of theMu family of transposable elements. All characterizedMu elements contain similar 200-bp terminal inverted repeats, yet the internal sequences of the elements may be completely unrelated. Non-Mutator stocks of maize have a 20–100-fold lower mutation rate relative to Mutator stocks, yet they contain multiple sequences that hybridize to theMu terminal inverted repeats. Most of these sequences do not cohybridize to internal regions of previously clonedMu elements. We have cloned two such sequences from the maize line B37, a non-Mutator inbred line. These sequences, termedMu4 andMu5, have an organization characteristic of transposable elements and possess 200-bpMu terminal inverted repeats that flank internal DNA, which is unrelated to other clonedMu elements.Mu4 andMu5 are both flanked by 9-bp direct repeats as has been observed for otherMu elements. However, we have no direct evidence that they have recently transposed because they have not been found in known genes. Although the internal regions ofMu4 andMu5 are not related by sequence similarity, both elements share an unusual structural feature: the terminal inverted repeats extend more than 100 bp internally fromMu-similar termini. The distribution of these elements in maize lines and related species suggests thatMu elements are an ancient component of the maize genome. Moreover, the structure of theMu termini and the fact thatMu termini are found flanking different internal sequences leads us to speculate thatMu termini once may have been capable of transposing as independent entities.     

Detection of protein and carbohydrate in the extracellular matrix released byCochliobolus heterostrophus during germination

The extracellular matrix (ECM) ofCochliobolus heterostrophus (anamorph:Bipolaris maydis) was made visible by gold/silver and FITC-lectin staining at different stages of germ tube development. A proteinaceous material was released from conidia as germ tubes began to emerge and continued to be released from the germ tube tip throughout elongation. A material that did not stain for protein was observed to surround germ tubes upon their elongation. At later stages of maturation, germ tubes were surrounded by a sheath of proteinaceous material. After 15 h of incubation, staining with the FITC-labeled Concanavalin A revealed that a carbohydrate material surrounded and extended between hyphae. The ECM extract was separated into two fractions which were shown by SDS-PAGE and HPLC analyses to consist of proteins and carbohydrates. The results demonstrate that the composition and physical structure of the ECM change over time. Thus, the ECM is not a static material. Rather, the components of the ECM appear to be laid down at different stages of fungal morphogenesis, possibly related to germ tube emergence, elongation, and maturation.     

The splicing of transposable elements and its role in intron evolution

Recent studies have demonstrated that transposable elements in maize and Drosophila are spliced from pre-mRNA. These transposable element introns represent the first examples of recent addition of introns into nuclear genes. The eight reported examples of transposable element splicing include members of the maize Ac/Ds and Spm/dSpm and the Drosophila P and 412 element families. The details of the splicing of these transposable elements and their relevance to models of intron origin are discussed.     

The adaptive role of transposable elements in the Drosophila genome

Transposable elements (TEs) are short DNA sequences with the capacity to move between different sites in the genome. This ability provides them with the capacity to mutate the genome in many different ways, from subtle regulatory mutations to gross genomic rearrangements. The potential adaptive significance of TEs was recognized by those involved in their initial discovery although it was hotly debated afterwards. For more than two decades, TEs were considered to be intragenomic parasites leading to almost exclusively detrimental effects to the host genome. The sequencing of the Drosophila melanogaster genome provided an unprecedented opportunity to study TEs and led to the identification of the first TE-induced adaptations in this species. These studies were followed by a systematic genome-wide search for adaptive insertions that allowed for the first time to infer that TEs contribute substantially to adaptive evolution. This study also revealed that there are at least twice as many TE-induced adaptations that remain to be identified. To gain a better understanding of the adaptive role of TEs in the genome we clearly need to (i) identify as many adaptive TEs as possible in a range of Drosophila species as well as (ii) carry out in-depth investigations of the effects of adaptive TEs on as many phenotypes as possible.     

An interchromosomal gene conversion of the Drosophila dunce locus identified with restriction site polymorphisms: a potential involvement of transposable elements in gene conversion

Summary Females heterozygous for the two alleles dnc ² and dnc ^M14 of the X-linked gene dunce (dnc), and carrying a copy of dnc ⁺ progeny X-chromosomes from recombination experiments. Restriction site polymorphisms have been used as genetic markers to follow the parentage of dnc locus segments in these chromosomes. All six chromosomes are identical with respect to the spectrum of restriction site markers they carry in the dnc ⁺ chromosomal region. In the progeny chromosomes, this region is comprised of sequences like the dnc ^M14 X-chromosome and the translocation copy of dnc ⁺. Sequences flanking the dnc gene in the progeny chromosomes are like the dnc ^M14 chromosome. Internal to the gene but near the 5 end, is a segment from the dnc ⁺ translocation which has apparently originated from an interchromosomal and premeiotic gene conversion event. In addition, two transposable elements have inserted into the progeny chromosomes, one towards the 5 end of dnc and the other near the 3 end. The insertion of these elements occurred premeiotically since all six chromosomes are structurally identical. The data are interpreted with respect to a potential role of transposable element transposition in the process of gene conversion.     

Cohabitation of KP and full-length P elements in the genome of MR strains inducing P-M-like hybrid dysgenesis in Drosophila melanogaster

Summary P strains of Drosophila melanogaster are characterized by the presence of both full-length and deletion derivatives of the transposable element P in their genome, and by their ability to induce the syndrome of hybrid dysgenesis among the progeny of certain intra-strain crosses, when introduced through the male parents. In contrast, strains belonging to the M' class, and which were also found to bear P element-homologous sequences, lack this ability and this has been attributed to the presence in the genome of most of these strains of a distinct class of deletion derivatives termed KP, which can suppress the action of functional P factors. Here we demonstrate that KP elements are present, next to full-length ones, in the genome of at least three strains which induce P-M-like dysgenic symptoms, including GD sterility. KP elements form the majority of the P-homologous sequences in the strains MR-h12, 23.5/CyL ⁴ and the latter's derivative 23.5 ^*/Cy. While the first one is a genuine P strain and the second one depicts a strong P cytotype, the third is a genuine M' strain. The hybrid dysgenesis induced by the two 23.5 MRF strains seems to be due, not primarily to the P elements, but to the action of hobo elements.     

Rates of movement of transposable elements inDrosophila melanogaster

Mobilization rates of nine families of transposable elements (P, hobo, FB, gypsy, 412, copia, blood, 297, andjockey) were estimated by using 182 lines. Lines were started from a completely isogenic population ofDrosophila melanogaster, carrying the markersepia as an indicator of possible contamination, and have been accumulating spontaneous mutations independently for 80 generations of brother-sister (or two double-first-cousin) matings. Transposable element movements have been analyzed in complete genomes by the Southern technique. Mobilization was a rare event, with an average rate of 10^–5 per site per generation. The most active element wasFB. In contrast, the retroelementsgypsy andblood did not move at all. Most changes in restriction patterns were consistent with rearrangements rather than with true transposition. The euchromatic or heterochromatic location of elements was tested by comparing insertion patterns from adults and salivary glands. Certain putative rearrangements involved heterochromatic copies of the retroelements412, copia or297. Clustering of movement across families was observed, suggesting that movement of different families may be non-independent. An association between modified insertion patterns and mutant effects on quantitative traits shows that spontaneous transposition events cause continuous variation.     

Role of transposable elements in the propagation of minisatellites in the rice genome

A survey of minisatellites (MSs) in 5.3 Mb of randomly selected rice DNA sequences from public databases was carried out to clarify the role of transposable elements (TEs) in the dispersal of MSs in the rice genome. The estimated frequency of MSs in this sample was one per 23.4 kb, and this frequency is approximately equivalent to that of Class I microsatellites in the rice genome. Of the MSs in the 5.3-Mb sequence sample, 82% were found to be present in multiple copies in the rice genome, and all of these were a part of TE sequences. In this study at least 61 TE groups were identified as MS carriers. It was also shown that the GC-rich MS pOs6.2H, which was previously reported to be one of the interspersed MSs in the rice genome, is a component of an En / Spm -like element. These results indicate that the majority of MSs in the rice genome are maintained in TEs, and amplified and dispersed as components of the TEs. The G+C content of the multi-locus MS sequences reflected that of the TE sequences containing those MSs, but no obvious bias towards the high G+C content of DNA was observed. Single locus MSs also did not show any obvious bias towards the high G+C content of DNA in the rice genome. In this respect, the MSs in the rice genome are quite different from those in the human genome: in the latter, the majority of MSs show an obvious bias towards the high G+C content of DNA.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by M.-A. Grandbastien     

Accumulation of anthraquinones in the reddish brown-colored polyoxin-resistant mutants of Cochliobolus heterostrophus

 Cochliobolus heterostrophus Pol2 and Pol5 mutants are pleiotropic, and each mutant gene is responsible for alterations of both unrelated phenotypes: reddish-brown pigmentation and polyoxin resistance. The three pigments accumulated in these mutants were isolated. Structural characterization by spectroscopic analyses indicated that these three pigments were polyhydroxyanthraquinones: emodin, chrysophanol, and citreorosein. Emodin is known to be an antidote against benzimidazole fungicide, although no antidoting activity against polyoxin was observed. Received: February 12, 2002 / Accepted: April 27, 2002     

Mitomycin C induces genomic rearrangements involving transposable elements in Drosophila melanogaster

Summary Mitomycin C was injected into the abdomen of male flies of the y ² sc¹ w^aG strain of Drosophila melanogaster. They were mated with females bearing attached-X chromosomes, and the male offspring (F₁) were analysed for the appearance of mutations in the X chromosome. We observed y ¹ and sc ⁺ reversions induced either by excision of mdg4 (gypsy) with retention of one long terminal repeat (LTR) or by insertion of a foreign sequence into mdg4, partial reversion of the w ^aG mutation, w ^aGw ^aGd, and unstable f mutations. The overall mutation frequency was considerably higher than in control flies of the y ² sc¹ w^aG strain. Possible mechanisms of genomic rearrangements induced by Mitomycin C, in particular the role of homologous recombination, are discussed.     

A branching process model for the evolution of transposable elements

A discrete-time multitype branching process model is presented for the evolution of transposable elements in haploid populations. An individual is classified as type i if it possesses i copies of the TE, i0. The general model incorporates copy-dependent selection and transposition, and recursion relations are derived for the distribution of the number of individuals of the various types. The asymptotic relative proportions of individuals of the different types is studied in the neutral case. The behavior of this equilibrium distribution is examined for various patterns of regulated transposition and deletion.