Microbotryum heliospermae, a new anther smut fungus parasitic on Heliosperma pusillum in the mountains of the European Alpine System

The members of the smut genus Microbotryum are pathogens of a wide range of host plant species from nine dicotyledonous families. Within the genus, the species sporulating in anthers of Caryophyllaceae form a monophyletic group that in recent years attracted much interest in various biological studies. The phylogenetic framework developed for species delimitation within Microbotryum revealed that high level host-specificity is a major feature of most caryophyllaceous anther smuts. However, the great number of anther smut specimens on diverse host plant species reported worldwide has still not been included in phylogenetic analyses due to the inaccessibility of recently collected specimens, and thus many species remain still undiscovered. In this study, anther smut specimens on Heliosperma pusillum originating from all main mountain ranges of the European Alpine System were examined using partial rDNA sequence and/or morphological analyses. The investigation revealed that all specimens are morphologically uniform and phylogenetically represent a monophyletic lineage, sister to Microbotryum lagerheimii complex on Atocion rupestre/Silene lacera/Silene vulgaris/Viscaria vulgaris. This lineage cannot be attributed to any of the previously described species, and therefore the smut in anthers of H. pusillum is described and illustrated here as a new species, Microbotryum heliospermae. The species is known from subalpine zone of the Alps, the Carpathians, the Dinaric Alps, and the Pyrenees, inhabiting host plants growing in open spring communities or semihumid mountain meadows.     

Chaos of Rearrangements in the Mating-Type Chromosomes of the Anther-Smut Fungus Microbotryum lychnidis-dioicae

Sex chromosomes in plants and animals and fungal mating-type chromosomes often show exceptional genome features, with extensive suppression of homologous recombination and cytological differentiation between members of the diploid chromosome pair. Despite strong interest in the genetics of these chromosomes, their large regions of suppressed recombination often are enriched in transposable elements and therefore can be challenging to assemble. Here we show that the latest improvements of the PacBio sequencing yield assembly of the whole genome of the anther-smut fungus, Microbotryum lychnidis-dioicae (the pathogenic fungus causing anther-smut disease of Silene latifolia), into finished chromosomes or chromosome arms, even for the repeat-rich mating-type chromosomes and centromeres. Suppressed recombination of the mating-type chromosomes is revealed to span nearly 90% of their lengths, with extreme levels of rearrangements, transposable element accumulation, and differentiation between the two mating types. We observed no correlation between allelic divergence and physical position in the nonrecombining regions of the mating-type chromosomes. This may result from gene conversion or from rearrangements of ancient evolutionary strata, i.e., successive steps of suppressed recombination. Centromeres were found to be composed mainly of copia-like transposable elements and to possess specific minisatellite repeats identical between the different chromosomes. We also identified subtelomeric motifs. In addition, extensive signs of degeneration were detected in the nonrecombining regions in the form of transposable element accumulation and of hundreds of gene losses on each mating-type chromosome. Furthermore, our study highlights the potential of the latest breakthrough PacBio chemistry to resolve complex genome architectures.     

Ultrastructural analysis of the behavior of the dimorphic fungus Microbotryum violaceum in fungus-induced anthers of female Silene latifolia flowers

Summary. The development of male organs is induced in female flowers of the dioecious plant Silene latifolia by infection with the fungus Microbotryum violaceum. Stamens in a healthy female flower grow only to stage 6, whereas those in an infected female flower develop to the mature stage (stage 12), at which the stamens are filled with fungal teliospores instead of pollen grains. To investigate these host–parasite interactions, young floral buds and fungus-induced anthers of infected female flowers were examined by electron microscopy following fixation by a high-pressure freezing method. Using this approach, we found that parasitic hyphae of this fungus contain several extracellular vesicles and have a consistent appearance up to stage 8. At that stage, parasitic hyphae are observed adjacent to dying sporogenous cells in the infected female anther. At stage 9, an increased number of dead and dying sporogenous cells is observed, among which the sporogenous hyphae of the fungus develop and form initial teliospores. Several types of electron-dense material are present in proximity to some fungi at this stage. The initial teliospores contain two types of vacuoles, and the fungus cell wall contains abundant carbohydrate, as revealed by silver protein staining. The sporogenous cell is probably sensitive to infection by the fungus, resulting in disruption. In addition, the fungus accelerates cell death in the anther and utilizes constituents of the dead host cell to form the mature teliospore. Correspondence and reprints (present address): Molecular Membrane Biology Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama 351-0198, Japan.     

DNA polymorphism in recombining and non-recombining mating-type-specific loci of the smut fungus Microbotryum

The population-genetic processes leading to the genetic degeneration of non-recombining regions have mainly been studied in animal and plant sex chromosomes. Here, we report population genetic analysis of the processes in the non-recombining mating-type-specific regions of the smut fungus Microbotryum violaceum. M. violaceum has A1 and A2 mating types, determined by mating-type-specific ‘sex chromosomes'' that contain 1–2 Mb long non-recombining regions. If genetic degeneration were occurring, then one would expect reduced DNA polymorphism in the non-recombining regions of this fungus. The analysis of DNA diversity among 19 M. violaceum strains, collected across Europe from Silene latifolia flowers, revealed that (i) DNA polymorphism is relatively low in all 20 studied loci (π∼0.15%), (ii) it is not significantly different between the two mating-type-specific chromosomes nor between the non-recombining and recombining regions, (iii) there is substantial population structure in M. violaceum populations, which resembles that of its host species, S. latifolia, and (iv) there is significant linkage disequilibrium, suggesting that widespread selfing in this species results in a reduction of the effective recombination rate across the genome. We hypothesise that selfing-related reduction of recombination across the M. violaceum genome negates the difference in the level of DNA polymorphism between the recombining and non-recombining regions, and may possibly lead to similar levels of genetic degeneration in the mating-type-specific regions of the non-recombining ‘sex chromosomes'' and elsewhere in the genome.     

Host-related genetic differentiation in the anther smut fungus Microbotryum violaceum in sympatric, parapatric and allopatric populations of two host species Silene latifolia and S. dioica

We investigated genetic diversity in West European populations of the fungal pathogen Microbotryum violaceum in sympatric, parapatric and allopatric populations of the host species Silene latifolia and S. dioica, using four polymorphic microsatellite loci. In allopatric host populations, the fungus was highly differentiated by host species, exhibiting high values of F(ST) and R(ST), and revealed clear and distinct host races. In sympatric and parapatric populations we found significant population differentiation as well, except for one sympatric population in which the two host species grew truly intermingled. The mean number of alleles per locus for isolates from each of the host species was significantly higher in sympatric/parapatric than in allopatric populations. This suggests that either gene flow between host races in sympatry, or in case of less neutral loci, selection in a more heterogeneous host environment can increase the level of genetic variation in each of the demes. The observed pattern of host-related genetic differentiation among these geographically spread populations suggest a long-term divergence between these host races. In sympatric host populations, both host races presumably come in secondary contact, and host-specific alleles are exchanged depending on the amount of fungal gene flow.     

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

Background

In recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.

Results

Here, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.

Conclusion

The significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Electronic supplementary material

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

The genome sequence of the biocontrol fungus Metarhizium anisopliae and comparative genomics of Metarhizium species

Background

Metarhizium anisopliae is an important fungal biocontrol agent of insect pests of agricultural crops. Genomics can aid the successful commercialization of biopesticides by identification of key genes differentiating closely related species, selection of virulent microbial isolates which are amenable to industrial scale production and formulation and through the reduction of phenotypic variability. The genome of Metarhizium isolate ARSEF23 was recently published as a model for M. anisopliae, however phylogenetic analysis has since re-classified this isolate as M. robertsii. We present a new annotated genome sequence of M. anisopliae (isolate Ma69) and whole genome comparison to M. robertsii (ARSEF23) and M. acridum (CQMa 102).

Results

Whole genome analysis of M. anisopliae indicates significant macrosynteny with M. robertsii but with some large genomic inversions. In comparison to M. acridum, the genome of M. anisopliae shares lower sequence homology. While alignments overall are co-linear, the genome of M. acridum is not contiguous enough to conclusively observe macrosynteny. Mating type gene analysis revealed both MAT1-1 and MAT1-2 genes present in M. anisopliae suggesting putative homothallism, despite having no known teleomorph, in contrast with the putatively heterothallic M. acridum isolate CQMa 102 (MAT1-2) and M. robertsii isolate ARSEF23 (altered MAT1-1). Repetitive DNA and RIP analysis revealed M. acridum to have twice the repetitive content of the other two species and M. anisopliae to be five times more RIP affected than M. robertsii. We also present an initial bioinformatic survey of candidate pathogenicity genes in M. anisopliae.

Conclusions

The annotated genome of M. anisopliae is an important resource for the identification of virulence genes specific to M. anisopliae and development of species- and strain- specific assays. New insight into the possibility of homothallism and RIP affectedness has important implications for the development of M. anisopliae as a biopesticide as it may indicate the potential for greater inherent diversity in this species than the other species. This could present opportunities to select isolates with unique combinations of pathogenicity factors, or it may point to instability in the species, a negative attribute in a biopesticide.

Electronic supplementary material

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

Integrated transcriptomic and proteomic analysis of the global response of Wolbachia to doxycycline-induced stress

The bacterium Wolbachia (order Rickettsiales), representing perhaps the most abundant vertically transmitted microbe worldwide, infects arthropods and filarial nematodes. In arthropods, Wolbachia can induce reproductive alterations and interfere with the transmission of several arthropod-borne pathogens. In addition, Wolbachia is an obligate mutualist of the filarial parasites that cause lymphatic filariasis and onchocerciasis in the tropics. Targeting Wolbachia with tetracycline antibiotics leads to sterilisation and ultimately death of adult filariae. However, several weeks of treatment are required, restricting the implementation of this control strategy. To date, the response of Wolbachia to stress has not been investigated, and almost nothing is known about global regulation of gene expression in this organism. We exposed an arthropod Wolbachia strain to doxycycline in vitro, and analysed differential expression by directional RNA-seq and label-free, quantitative proteomics. We found that Wolbachia responded not only by modulating expression of the translation machinery, but also by upregulating nucleotide synthesis and energy metabolism, while downregulating outer membrane proteins. Moreover, Wolbachia increased the expression of a key component of the twin-arginine translocase (tatA) and a phosphate ABC transporter ATPase (PstB); the latter is associated with decreased susceptibility to antimicrobials in free-living bacteria. Finally, the downregulation of 6S RNA during translational inhibition suggests that this small RNA is involved in growth rate control. Despite its highly reduced genome, Wolbachia shows a surprising ability to regulate gene expression during exposure to a potent stressor. Our findings have general relevance for the chemotherapy of obligate intracellular bacteria and the mechanistic basis of persistence in the Rickettsiales.     

Host adaptation in the anther smut fungus Ustilago violacea (Microbotryum violaceum): infection success,spore production and alteration of floral traits on two host species and their F1-hybrid

It is often assumed that host specialization is promoted by trade-offs in the performance of parasites on different host species, but experimental evidence for such trade-offs is scant. We studied differences in performance among strains of the anther smut fungus Ustilago violacea from two closely related host plant species, Silene alba and S. dioica, on progeny of (1) the host species from which they originated, (2) the alternative host species, and (3) inter-specific hybrids. Significant intra-specific variation in the pathogen was found for both infection success on a range of host genotypes (virulence) and components of spore production per infected host (aggressiveness) (sensu Burdon 1987). Strains did not have overall higher virulence on conspecifics of their host of origin than on strains from the heterospecific host, but they did have a significantly (c. 3 times) higher spore production per infected male host. This finding suggests that host adaptation may have evolved with respect to aggressiveness rather than virulence. The higher aggressiveness of strains on conspecifics of their host of origin resulted both from higher spore production per infected flower (spores are produced in the anthers), and greater ability to stimulate flower production on infected hosts. The latter indicates the presence of adaptive intraspecific variation in the ability of host manipulation. As transmission of the fungus is mediated by insects that are both pollinators of the host and vectors of the disease, we also assessed the effect of strains on host floral traits. Infection resulted in a reduction of inflorescence height, flower size, and nectar production per flower. Strains did not differ in their effect on nectar production, but infection with strains from S. alba resulted in a stronger reduction of inflorescence height and petal size on both host species. Vectors may therefore in principle discriminate among hosts infected by different strains and affect their efficiency of transmission. Contrary to assumptions of recent hypotheses about the role of host hybrids in the evolution of parasites, hybrids were not generally more susceptible than parental hosts. It is therefore unlikely that the rate of evolution of the pathogen on the parental species is slowed down by selection for specialization on the hybrids.     

The effect of sterilizing diseases on host abundance and distribution along environmental gradients

This study analyses the effect of host-specific pathogens on range restriction of their hosts across environmental gradients at population margins. Sterilizing diseases can limit host range by causing large reductions in population size in what would otherwise be the central area of a species range. Diseases showing frequency-dependent transmission can also pull back a population from its disease-free margin. A wide range of disease prevalence versus abundance patterns emerge which often differ from the classical expectation of increasing prevalence with increasing abundance. Surprisingly, very few empirical studies have investigated the dynamics of disease across environmental gradients or at range limits.     

Comparative genomic analysis reveals distinct genotypic features of the emerging pathogen Haemophilus influenzae type f

Background

The incidence of invasive disease caused by encapsulated Haemophilus influenzae type f (Hif) has increased in the post-H. influenzae type b (Hib) vaccine era. We previously annotated the first complete Hif genome from a clinical isolate (KR494) that caused septic shock and necrotizing myositis. Here, the full genome of Hif KR494 was compared to sequenced reference strains Hib 10810, capsule type d (Hid) Rd Kw20, and finally nontypeable H. influenzae 3655. The goal was to identify possible genomic characteristics that may shed light upon the pathogenesis of Hif.

Results

The Hif KR494 genome exhibited large regions of synteny with other H. influenzae, but also distinct genome rearrangements. A predicted Hif core genome of 1390 genes was shared with the reference strains, and 6 unique genomic regions comprising half of the 191 unique coding sequences were revealed. The majority of these regions were inserted genetic fragments, most likely derived from the closely-related Haemophilus spp. including H. aegyptius, H. haemolyticus and H. parainfluenzae. Importantly, the KR494 genome possessed several putative virulence genes that were distinct from non-type f strains. These included the sap2 operon, aef3 fimbriae, and genes for kanamycin nucleotidyltranserase, iron-utilization proteins, and putative YadA-like trimeric autotransporters that may increase the bacterial virulence. Furthermore, Hif KR494 lacked a hisABCDEFGH operon for de novo histidine biosynthesis, hmg locus for lipooligosaccharide biosynthesis and biofilm formation, the Haemophilus antibiotic resistance island and a Haemophilus secondary molybdate transport system. We confirmed the histidine auxotrophy and kanamycin resistance in Hif by functional experiments. Moreover, the pattern of unique or missing genes of Hif KR494 was similar in 20 Hif clinical isolates obtained from different years and geographical areas. A cross-species comparison revealed that the Hif genome shared more characteristics with H. aegyptius than Hid and NTHi.

Conclusions

The genomic comparative analyses facilitated identification of genotypic characteristics that may be related to the specific virulence of Hif. In relation to non-type f H. influenzae strains, the Hif genome contains differences in components involved in metabolism and survival that may contribute to its invasiveness.

Electronic supplementary material

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

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.     

Dissecting the fungal biology of Bipolaris papendorfii: from phylogenetic to comparative genomic analysis

Bipolaris papendorfii has been reported as a fungal plant pathogen that rarely causes opportunistic infection in humans. Secondary metabolites isolated from this fungus possess medicinal and anticancer properties. However, its genetic fundamental and basic biology are largely unknown. In this study, we report the first draft genome sequence of B. papendorfii UM 226 isolated from the skin scraping of a patient. The assembled 33.4 Mb genome encodes 11,015 putative coding DNA sequences, of which, 2.49% are predicted transposable elements. Multilocus phylogenetic and phylogenomic analyses showed B. papendorfii UM 226 clustering with Curvularia species, apart from other plant pathogenic Bipolaris species. Its genomic features suggest that it is a heterothallic fungus with a putative unique gene encoding the LysM-containing protein which might be involved in fungal virulence on host plants, as well as a wide array of enzymes involved in carbohydrate metabolism, degradation of polysaccharides and lignin in the plant cell wall, secondary metabolite biosynthesis (including dimethylallyl tryptophan synthase, non-ribosomal peptide synthetase, polyketide synthase), the terpenoid pathway and the caffeine metabolism. This first genomic characterization of B. papendorfii provides the basis for further studies on its biology, pathogenicity and medicinal potential.     

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.     

Effects of heat and UV radiation on the mobilization of transposon mariner-Mos1

There are many complex interactions between transposable elements (TEs) and host genomes. Environmental changes that induce stressful conditions help to contribute for increasing complexity of these interactions. The transposon mariner-Mos1 increases its mobilization under mild heat stress. It has putative heat shock elements (HSEs), which are probably activated by heat shock factors (HSFs). Ultraviolet radiation (UVC) is a stressor that has been suggested as able to activate heat shock protein genes (Hsp). In this study, we test the hypothesis that if UVC induces Hsp expression, as heat does, it could also promote mariner-Mos1 transposition and mobilization. The Drosophila simulans white-peach is a mutant lineage that indicates the mariner-Mos1 transposition phenotypically through the formation of mosaic eyes. This lineage was exposed to UVC or mild heat stress (28 °C) in order to evaluate the induction of mariner-Mos1 expression by RT-qPCR, as well as the mariner-Mos1 mobilization activity based on the count number of red spots in the eyes. The effects of both treatments on the developmental time of flies and cell cycle progression were also investigated. Both the analysis of eyes and mariner-Mos1 gene expression indicate that UVC radiation has no effect in mariner-Mos1 transposition, although heat increases the expression and mobilization of this TE soon after the treatment. However, the expression of Hsp70 gene increased after 24 h of UVC exposure, suggesting different pathway of activation. These results showed that heat promotes mariner-Mos1 mobilization, although UVC does not induce the expression or mobilization of this TE.

Electronic supplementary material

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

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.