Evolutionary genomics revealed interkingdom distribution of Tcn1-like chromodomain-containing Gypsy LTR retrotransposons among fungi and plants
© Novikova et al; licensee BioMed Central Ltd. 2010
Received: 15 October 2009
Accepted: 8 April 2010
Published: 8 April 2010
Chromodomain-containing Gypsy LTR retrotransposons or chromoviruses are widely distributed among eukaryotes and have been found in plants, fungi and vertebrates. The previous comprehensive survey of chromoviruses from mosses (Bryophyta) suggested that genomes of non-seed plants contain the clade which is closely related to the retrotransposons from fungi. The origin, distribution and evolutionary history of this clade remained unclear mainly due to the absence of information concerning the diversity and distribution of LTR retrotransposons in other groups of non-seed plants as well as in fungal genomes.
In present study we preformed in silico analysis of chromodomain-containing LTR retrotransposons in 25 diverse fungi and a number of plant species including spikemoss Selaginella moellendorffii (Lycopodiophyta) coupled with an experimental survey of chromodomain-containing Gypsy LTR retrotransposons from diverse non-seed vascular plants (lycophytes, ferns, and horsetails). Our mining of Gypsy LTR retrotransposons in genomic sequences allowed identification of numerous families which have not been described previously in fungi. Two new well-supported clades, Galahad and Mordred, as well as several other previously unknown lineages of chromodomain-containing Gypsy LTR retrotransposons were described based on the results of PCR-mediated survey of LTR retrotransposon fragments from ferns, horsetails and lycophytes. It appeared that one of the clades, namely Tcn1 clade, was present in basidiomycetes and non-seed plants including mosses (Bryophyta) and lycophytes (genus Selaginella).
The interkingdom distribution is not typical for chromodomain-containing LTR retrotransposons clades which are usually very specific for a particular taxonomic group. Tcn1-like LTR retrotransposons from fungi and non-seed plants demonstrated high similarity to each other which can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission.
Retrotransposons are a class of mobile genetic elements, which use reverse transcription in their transposition. Five orders of retrotransposons are recognized: those having long terminal repeats (LTRs) (LTR retrotransposons); those lacking LTRs (non-LTR retrotransposons); DIRS retrotransposons; Penelope-like retrotransposable elements; and short interspersed nuclear elements (SINEs). According to the modern classification, LTR retrotransposons are divided into several superfamilies: Copia (Pseudoviridae), Gypsy (Metaviridae), Bel-Pao, Retrovirus (Retroviridae), and ERV .
Chromodomain-containing LTR retrotransposons or chromoviruses are the most widespread lineage of Gypsy LTR retrotransposons and are present in genomes of fungi as well as in plants and vertebrates [2, 3]. The characteristic feature of chromoviruses is the presence of an additional domain - the chromodomain (CHD). CHDs are present in various eukaryotic proteins involved in chromatin remodeling and regulation of gene expression during development [4–6]. CHDs perform a wide range of diverse functions including chromatin targeting and proteinDNA/RNA interactions . Recently, it has been shown that the CHDs target integration of new LTR retrotransposon copies into heterochromatin by recognizing histone modifications .
Our previous comprehensive survey of chromoviruses from mosses (Bryophyta) suggested that the diversity of CHD-containing Gypsy LTR retrotransposons in plant genomes is underestimated [8, 9]. There are four wellknown CHD-containing Gypsy LTR retrotransposon clades widely distributed among gymnosperms and angiosperms: Tekay, CRM, Galadriel and Reina [2, 3]. Four novel clades were found to be present in mosses. Moreover, we showed that representatives from one of the moss-specific clades are more closely related to retrotransposons from fungi than to retrotransposons from plants. Although we proposed that the retrotransposons from this clade could have been 'retained' from the last common ancestor of Fungi/Metazoa lineage and plants, the origin of this clade remains unclear.
The questions addressed in current investigation are as follows: (1) what kind of Gypsy LTR retrotransposons from fungi are closely related to the LTR retrotransposons detected in mosses; (2) how widely those clades which were previously identified in mosses are distributed among other non-seed plants including lycophytes, ferns and horsetails; and (3) what is the origin of the clade which is common for non-seed plants and fungi. The in silico analysis of Gypsy LTR retrotransposons in 25 species of fungi, genomes of which available in public databases, along with survey of related LTR retrotransposons from whole genome sequences (WGS) and expressed sequence tags (ESTs) of diverse plants including spikemoss Selaginella moellendorffii (Lycopodiophyta) and PCR-based screening of ferns, horsetails and lycophytes showed that a common clade of CHD-containing Gypsy LTR retrotransposons can be found in mosses, lycophytes from genus Selaginella, and basidiomycetes. According to the classification of CHD-containing Gypsy LTR retrotransposons proposed by Gorinsek et al. (2004)  this clade has name Tcn1. It seems that Tcn1 is a unique clade of chromoviruses which has a wide inter-kingdom distribution. Tcn1-like LTR retrotransposons from fungi and non-seed plants demonstrated higher similarity to each other in comparison with LTR retrotransposons from other clades. This can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission.
Gypsy LTR retrotransposons survey from fungal genomes
List of fungal species, genomes of which were analyzed in silico in present study
Species and strain
Chaetomium globosum CBS 148.51
Fusarium oxysporum 4286 FGSC
Fusarium verticillioides 7600
Nectria haematococca MPVI
Podospora anserina S mat+
Trichoderma reesei QM6a
Trichoderma virens Gv29-8
Aspergillus clavatus NRRL 1
Aspergillus niger ATCC1015
Aspergillus terreus NIH2624
Coccidioides immitis RS
Histoplasma capsulatum NAm1
Uncinocarpus reesii 1704
Sclerotinia sclerotiorum 1980
Botrytis cinerea B05.10
Alternaria brassicicola ATCC 96866
Pyrenophora tritici-repentis Pt-1C-BFP
Stagonospora nodorum SN15
Coprinus cinereus Okayama7#130
Laccaria bicolor S238N
Postia placenta MAD-698
Puccinia graminis f. sp. tritici
Batrachochytrium dendrobatidis JEL423
In addition to the abovelisted enzymatic domains, a deoxyuridine triphosphatase domain (dUTPase) has been found in several LTR retrotransposons from the basidiomycete Postia placenta MAD-698. The location of this domain varied among diverse families of PosPlaTy3 elements. It can be found either at amino-terminus (PosPlaTy3-3) and carboxyl-terminus of Pol (PosPlaTy3-4) or between PR and RT domains (PosPlaTy3-5) (Figure 3). The presence of dUTPase in LTR retrotransposon sequences has been described earlier for the elements from a basidiomycete Phanerochaete chrysosporium and an ascomycete Tuber melanosporum [14, 15]. The role and origin of this domain remained unclear. Moreover, it seems that the described LTR retrotransposons acquire this domain independently from different sources (Additional file 2) . It was assumed that the presence of dUTPase allows viruses that contain this domain to replicate in non-dividing cells, in which cellular dUTPase activity is absent because replication of DNA does not occur .
Tnc1 clade is found in non-seed plants
The whole sequence of SM-Tcn1 LTR retrotransposon was obtained from WGS. SM-Tcn1 is 5704 bp in length and carries two putative ORFs. ORF1 or gag (969 bp in length) encodes a 323 amino acid (aa) protein with strong similarity to retroviral Gag proteins (pfam03732). ORF2 or pol (3714 bp) encodes a 1238 aa polyprotein, with characteristic retroviral aspartyl protease (PR), RT, Int, and CHD domains (Figure 3; Additional file 3). SM-Tcn1 possesses 440 bp LTRs with conserved features, including the dinucleotide end sequences (TG...CA). Target site duplication (AACAC...AACAC) was also detected for the described copy of SM-Tcn1. The LTRs contained TATA and CAAT boxes. No putative primer-binding site (PBS) was found. The PBS is necessary for initiation of reverse transcription and synthesis of the first strand complementary 5'LTR sequence . Typical PBS is located near the 5'LTR and complementary to the 3' terminal nucleotides of the primer tRNA used for initiation. Other known mechanism for initiation of the first strand synthesis is self-priming, in this case a sequence derived from LTR is located just downstream of the 5'LTR [18, 19]. However, evidences were found neither for tRNA priming nor for self-priming of SM-Tcn1 LTR retrotransposon. The sequence presented between 5' LTR and gag is conservative and thymine-rich (Figure 3). The possible mechanism for initiation of reverse transcription of SM-Tcn1 remained unclear. A polypurine tract (PPT) was detected immediately upstream of the 3'LTR. The PPT sequence is involved in second-strand DNA synthesis. The BLAST search (blastn) of full-length SM-Tcn1 retrotransposon indicated the presence of more than 200 hits in the S. moellendorffii genome. The close examination of identified copies showed that they have in average 91.7% nucleotide identity with the original SM-Tcn1 sequence. A BLAST search (tblastn) using SM-Tcn1 putative Pol protein as a query yielded more than 500 hits. The size of S. moellendorffii genome is only ~100 Mbp , thus SM-Tcn1 retrotransposons comprise ca. 1.5% of genomic sequence.
Gypsy LTR retrotransposons from non-seed vascular plants
List of monilophytes and lycophytes, which were analyzed experimentally in present study
Ophioglossum vulgatum L.
Botrychium multifidum (Gmelin) Rupr.
Pteridium aquilinum (L.) Kuhn
Adiantum pedatum L.
Asplenium viride Huds.
Asplenium ruta-muraria L.
Athyrium sinense Rupr.
Athyrium distentifolium Tausch ex Opiz
Athyrium monomachii (Kom.) Kom.
Cystopteris fragilis (L.) Bern.
Woodsia polystichoides DC Eaton
Woodsia subcordata Turcz.
Phegopteris connectilis (Michx.) Watt.
Matteuccia struthiopteris (L.) Tod.
Onoclea sensibilis L.
Dryopteris expansa (Presl) Fraser-Jenk. & Jermy
Dryopteris filix-mas (L.) Schott
Dryopteris crassirhizoma Nak.
Dryopteris carthusiana (Vill.) Fuchs
Polystichum braunii (Spenner) Fée
Polystichum tripteron (Kunze) Presl
Polypodium vulgare L.
Pyrrosia lingua (Thunb.) Farw.
Salvinia natans (L.) All.
Equisetum hiemale L.
Equisetum fluviatile L.
Isoetes histrix Bory
Selaginella kraussiana (Kunze) A. Braun
Selaginella pulvinata (Hook. & Grev.) Maxim.
Diphasiastrum complanatum (L.) Holub
Huperzia squarrosa (Forst.) Trevis.
Lycopodium alpinum L.
Lycopodium clavatum L.
Lycopodium magellanicum Hert. ex Nessel
Lycopodium japonicum Thunb. ex Murr.
Lycopodium annotinum L.
The estimated diversity of monilophytes (= Infradivision Moniliformopses) is about 9000 species and includes horsetails, whisk ferns, and all eusporangiate and leptosporangiate ferns . Most of the species examined in the present study were leptosporangiate ferns from the order Polypodiales, class Polypodiopsida. This order covers more than 80% of current known diversity of ferns . Additionally, one representative of heterosporous ferns, Salvinia natans (Polypodiopsida, Salviniales), two ophioglossoid ferns (Psilotopsida, Ophioglossales) and three horsetails (Equisetopsida, Equisetales) were included . Lycophytes are much less diverse in comparison with monilophytes and comprise less than 1% of extant land plants (around 1200 living species). Three major lineages are distinguished among lycophytes: clubmosses and firmosses (Lycopodiaceae), spikemosses (Selaginellaceae), and quillworts (Isoetaceae) . Among lycophytes included in the present study are Isoetes and Huperzia species (Isoetaceae) as well as two Selaginella species (Selaginellaceae), which belong to the class Isoetopsida, and seven diverse species from Lycopodiaceae (Lycopodiales, Lycopodiopsida).
The presence of CHD-containing Gypsy LTR retroelements among the listed plants was tested by amplifying genomic DNA with previously developed degenerate oligonucleotide primers [8, 24]. Consistent with the spacing of reverse transcriptase (RT) domains, the amplified PCR products were approximately 320 bp in length. In total, 98 clones with sequence similarity to known RT sequences were isolated, of which 76 were from monilophytes and 22 from lycophytes. The preliminary blastp search revealed that 10 clones were not from CHDcontaining LTR retrotransposons but were from Athila-like Gypsy elements (Additional file 4). Many representatives of this clade possess not only classical gag and pol sequences, but also an additional open reading frame that might encode an env-like protein [25, 26].
The phylogenetic relationships among obtained clones and known CHD-containing Gypsy LTR retrotransposons, extracted from databases, were reconstructed using neighbor-joining (NJ) analysis based on the multiple alignment of nucleotide sequences of RT fragment (Figure 4). The Gypsy LTR retrotransposons from Drosophila melanogaster were used as an outgroup. The newlyidentified LTR retrotransposon grouped into four large clusters on the phylogenetic tree. The group of clones from diverse monilophytes and one LTR retrotransposon from lycophytes Lycopodium alpinum (LycAlpTy3-1 clone) form a common group with previously described retroelements from mosses Tetraphis pellucida and Vesicularia dubyana . Although the bootstrap support is below 50% (data not shown), this new clade (named "Galahad") seems to be a sister group to the Galadriel clade. Galahad appears to be one of the oldest widely distributed clades of CHD-containing LTR retrotransposons from plants. Since Galahad clade was found in all non-seed plants including mosses, the probable age of this clade would be in the range of 400-700 Myr, which is estimated time divergence of liverworts and mosses from vascular plants .
The second group is formed by LTR retrotransposons from both monilophytes and lycophytes. The phylogenetic analysis did not provide support for a monophyletic origin of this cluster. Moreover, the relationships inside the cluster remained unclear, with the exception of several lineages. One of the lineages ('a' in Figure 4), contained members from lycophytes in the family Lycopodiaceae: Lycopodium clavatum (LycClavGty3 clones), L. japonicum (LycJapGty3-1 clone), and Diphasiastrum complanatum (DiphComGty3 clones). Another lineage seems to have had a long-term association with fern genomes, since it appears to be widely distributed among leptosporangiate ferns and can be found in Dennstaedtiaceae (Pteridium aquilinum), Pteridaceae (Adiantum pedatum), Aspleniaceae (Asplenium viride), Woodsiaceae (Athyrium distentifolium and Cystopteris fragilis), and Dryopteridaceae (Dryopteris crassirhizoma and Polystichum tripteron) (lineage 'b' on Figure 4). Additionally, five satellite lineages, represented mostly by single clones can be found on the phylogenetic tree.
The largest group is represented by 37 LTR retrotransposons. Three clearly separated clusters can be found inside this group (marked as 'd', 'f', and Mordred on Figure 4). One of these clusters is formed by LTR retrotransposons from Ophioglossaceae (Botrychium multifidum), Aspleniaceae (A. viride) and Dryopteridaceae (Dryopteris expansa, D. carthusiana, and D. crassirhizoma). The second cluster has a bootstrap support of 77% and is represented by clones isolated from Lycopodium annotinum (Lycopodiaceae) and Woodsia polystichoides (Woodsiaceae). The last monophyletic cluster, clade Mordred, is the largest, well-supported clade which is widely distributed among representatives of all investigated classes except Isoetopsida. It was found to be present in fern genomes from families Ophioglossaceae (Ophioglossum vulgatum and Botrychium multifidum), Pteridaceae (Adiantum pedatum), Woodsiaceae (Athyrium sinense, A. distentifolium, A. monomachii, C. fragilis, and W. polystichoides), Onocleaceae (Onoclea sensibilis), and Dryopteridaceae (Dryopteris filix-mas and D. crassirhizoma); in horsetails Equisetum hiemale and E. fluviatile (Equisetaceae); and lycophyte Lycopodium magellanicum (Lycopodiaceae).
The Tcn1-like LTR retrotransposons were detected only in Selaginella species, S. kraussiana (SelKraGty3 clones) and S. pulvinata (SelPulGty3 clones) in addition to the previously described LTR retrotransposons from mosses and SM-Tcn1 LTR retrotransposon from Selaginella moellendorffii (Figure 4) . The absence of Tcn1-like clones isolated from other lycophytes, ferns, and horsetails can be explained by failed PCR amplification due to the high divergence of Tcn1-like elements in genomes of these species or, more likely, by lack of these elements from their genomes.
Tcn1-like LTR retrotransposons: 'retained' or horizontally transmitted?
As a rule chromoviruses clades are specific for a particular group of eukaryotic organisms such as Ascomycota fungi (Nessie, Pyret, Maggy, Pyggy, MGLR3, Yeti, Coccy1, Coccy2, Polly, Afut1, Tf1, Ty3, and Afut4), Basidiomycota fungi (MarY1, Laccy1, Laccy2, Tcn2, Puccy1, and Puccy2), or plants (Reina, CRM, Tekay, Galadriel, and Chlamyvir as well as additional less investigated clades from mosses) [3, 8]. In the light of such specificity, it was unexpected to find a clade containing elements from basidiomycetes and non-seed plants. Nevertheless, it seems that Tcn1 clade has an interkingdom distribution and can be found in a number of fungi, diverse mosses (Bryophyta) as well as in lycophytes (genus Selaginella). Such a wide distribution makes the Tcn1 clade unique among the CHD-containing Gypsy LTR retrotransposons. The interkingdom distribution of Tcn1 clade could be the result of horizontal transmission (HT) of LTR retrotransposons among fungi and plants; otherwise, Tcn1-like LTR retrotransposons could have been 'retained' by mosses and lycophytes from the most recent common ancestor of plants and Fungi/Metazoa lineage of eukaryotes [8, 28].
The hypothesis of 'retained' genes is based on the observation that EST data of Physcomitrella contained a fraction of transcripts derived from putative genes ('retained' genes), which are not present in seed plants but can be found in other kingdoms including fungi. It was proposed that such retained genes along with Physcomitrella-specific (or moss-specific) genes encode functions that make mosses unique in terms of physiology and metabolism . We used these data and compared the levels of similarity for RT-Int fragments from Tcn1-like LTR retrotransposons and two putatively retained genes from Physcomitrella, which showed a high similarity with functional genes from fungi: uric acid-xanthine permease (uapA, TIGR00801) and inorganic phosphate transporter (Pho88, pfam10032).
Amino acid divergences of proteins and RT-Int fragments of CHD-containing Gypsy LTR retrotransposons from Tcn1, Pyggy and Pyret clades
Genes or LTR retrotransposons
Amino acid identity (%)
Evolutionary rate (10-9)b
inorganic phosphate transporter (Pho88)a
Physcomitrella patens (XM_001783642)/Schizosaccharomyces pombe (NM_001019478)
P. patens (XM_001783642)/Coprinus cinereus (XM_001836748)
P. patens/Phanerochaete chrysosporium (JGI: scaffold_7 [325337..326089])
P. patens/Cryptococcus neoformans (XM_569621)
C. cinereus/Ph. chrysosporium
C. cinereus/C. neoformans
Ph. chrysosporium/C. neoformans
uric acid-xanthine permease (uapA)a
P. patens (XM_001784081)/Ashbya gossypii (NM_212305)
P. patens/C. cinereus (XM_001839036)
P. patens/Ph. chrysosporium (JGI: scaffold_22 [347275..348820])
P. patens/C. neoformans (AF542528)
C. cinereus/Ph. chrysosporium
C. cinereus/C. neoformans
Ph. chrysosporium/C. neoformans
Tcn1 C. neoformans/Ccchromovir-1 C. cinereus
Tcn1 C. neoformans/PcMetavir6 Ph. chrysosporium
Tcn1 C. neoformans/BatDenTy3-1 B. dendrobatidis
PcMetavir6 Ph. chrysosporium/Ccchromovir-1 C. cinereus
PcMetavir6 Ph. chrysosporium/BatDenTy3-1 Batrachochytrium dendrobatidis
BatDenTy3-1 B. dendrobatidis/Ccchromovir-1 C. cinereus
Tcn1 C. neoformans/PpatensLTR1 P. patens
Tcn1 C. neoformans/SM-Tcn1 Selaginella moellendorffii
PcMetavir6 Ph. chrysosporium/PpatensLTR1 P. patens
PcMetavir6 Ph. chrysosporium/SM-Tcn1 S. moellendorffii
PpatensLTR1 P. patens/SM-Tcn1 S. moellendorffii
PyrTriTy3-2 Pyrenophora tritici-repentis/NecHaemTy3-4Nectria haematococca
PyrTriTy3-2 P. tritici-repentis/ChaGloTy3-8 Chaetomium globosum
PyrTriTy3-2P. tritici-repentis/grh Magnaporthe grisea
PyrTriTy3-2 P. tritici-repentis/Dane4Aspergillus nidulans
AltBraTy3-2Alternaria brassicicola/NecHaemTy3-4N. haematococca
AltBraTy3-2A. brassicicola/Dane4 A. nidulans
grh M. grisea/Dane4 A. nidulans
skippy Fusarium oxysporum/PyrTriTy3-1 P. tritici-repentis
skippy F. oxysporum/AFLAV Aaspergillus flavus
AFLAV A. flavus/PyrTriTy3-1 P. tritici-repentis
It seems to be that the hypothesis of 'retained' genes cannot be implemented as explanation for Tcn1 clade distribution since investigated RT-Int fragments of Tcn1-like LTR retrotransposons from fungi and plants have higher similarity to each other than functional proteins which were proposed to be 'retained' . RT-Int fragments from Tcn1-like LTR retrotransposons have average similarity 49%. Moreover, evolutionary rates estimated for Tcn1 LTR retrotransposons appeared to be less than evolutionary rates for 'retained' genes or other LTR retrotransposons (Table 3).
Despite a number of whole genome sequence studies, the distribution and diversity of CHD-containing Gypsy LTR retrotransposons is still poorly understood. Current knowledge of distribution and evolution of this group of mobile elements has been mainly obtained from diverse model organisms . The quickly generated massive data sets (for example, WGS and EST databases) provide a great opportunity to perform detailed analysis for non-model organisms. Experimental data accumulation also should not be neglected. In the absence of information about genomic sequences from monilophytes, PCR screening seems to be very useful for isolation and characterization of new LTR retrotransposons from this group of non-seed vascular plants. Two new well-supported clades, Galahad and Mordred, as well as several other previously unknown lineages of CHD-containing Gypsy LTR retrotransposons were described based on the results of PCR-mediated survey of RT fragments from ferns, horsetails and lycophytes.
Several criteria can be used for HT event recognition. The first criterion is inconsistencies between the phylogenies of transposable elements (TEs) and host species [29, 34]. There are potential problems with application of this criterion for HT detection. Multiple transposable element lineages can be present within genomes. Moreover, transposable elements are multicopy components of genomes. Comparisons of paralogous copies instead of orthologs along with varying rates of their sequence evolution are the main sources for incongruence in phylogenetic analysis, this could be misidentified as HT. The second criterion, which seems to offer the strongest evidence, is a higher degree of observed sequence similarity for transposable elements than for functional genes, so called 'slowdown effect on evolutionary rates'. Once inserted, a new copy of transposable element is presumed to evolve without functional constrains. Thus, all types of mutations should have an equal chance to be fixed . The lower than expected sequence divergence of TEs in comparison with non-mobile nuclear genes of the host species can be explained either by strong selective constraints in TE sequence coupled with a strict vertical transmission, or by horizontal transfer [31, 36, 37]. The third criterion of inferring HT is the discontinuous distribution of TEs among closely related taxa, i.e., presence of a TE in one lineage and its absence in a sister lineage. Such discontinuous distribution could be due to random loss of TEs, ancestral polymorphism, or independent sorting of copies into descendant species. By itself, this kind of evidence provides only weak support for HT since TE can be lost through population dynamics or ecological forces that are difficult to reconstruct [38, 39].
All three criteria are satisfied in case of Tcn1-like LTR retrotransposons. They demonstrated patchy distribution among fungi and plants (Figure 5). They were found in all investigated mosses, but only in a few lycophytes and they absent in basal lineages of green plants such as green and red algae as well as in all seed plants investigated so far. The hypothesis of 'retained' genes in moss Physcomitrella represented an attractive alternative to horizontal transmission as an explanation of the phylogenetic inconsistencies as well as the existence of a number of functional genes in Physcomitrella genome, which seem to have non-plant origin and can be found in bacteria, fungi and protozoa but not in higher plants . Nevertheless, the most important feature of Tcn1-like LTR retrotransposons in the context of HT is their lower evolutionary rates in comparison with other groups of CHD-containing LTR retrotransposons. The close examination and comparison of evolutionary rates for LTR retrotransposons including representatives of Tcn1, Pyggy and Pyret clades, and evolutionary rates estimated for putatively 'retained' genes suggests that a horizontal transmission of Tcn1-like LTR retrotransposons took place among fungi and the last common ancestor (LCA) of mosses and lycophytes (Table 3 and Figure 5) . Alternatively, it is possible, but highly unlikely, that two independent acts of HT occurred. First HT event could happen among fungi and LCA of mosses since all investigated mosses contain Tcn1-like LTR retrotransposons . The second HT could occur among fungi and LCA of Selaginella since only representatives of this genus carry this group of retrotransposons among all investigated lycophytes (Figure 4 and Figure 5). It is necessary to note that despite HT seeming to be a preferable explanation for the observed distribution. The evidence is not strong enough to discard other explanations; such as selective pressure coupled with vertical transmission of retrotransposons in genomes of non-seed plants and loss of these elements by other plants.
Another putative case of HT based on the results of present survey of LTR retrotransposons from fungal species was found for PyrTriTy3-2 LTR retrotransposon from Pyrenophora tritici-repentis Pt-1C-BFP (Dothideomycetes). PyrTriTy3-2 belongs to Pyggy clade and appeared to be more closely related to LTR retrotransposons from Sordariomycetes (NecHaemTy3-4 from Nectria haematococca MPVI and ChaGloTy3-8 from Chaetomium globosum CBS 148.51) than to the elements from other Dothideomycetes such as AltBraTy3-2 from Alternaria brassicicola ATCC 96866, REAL from Alternaria alternata (AB025309) , and PYGGY from Pyrenophora graminea (AF533704)  (Figure 2). The pairwise comparisons of RT-Int fragments and investigation of evolutionary rates for retrotransposons from Pyggy and Pyret clades revealed the unexpectedly high similarity between PyrTriTy3-2 and NecHaemTy3-4 (77% identical amino acids), much higher than between any other retrotransposons from Pyggy or Pyret clades, and at least two times lower evolutionary rate in the couple PyrTriTy3-2/NecHaemTy3-4 than in comparisons of other LTR retrotransposons (Table 3).
The high similarity, phylogenetic inconsistencies, as well as lower evolutionary rates could be explained by very strict evolutionary constraints or a HT event. However, taking into consideration that the high selective pressure could be implemented only in the case of functional importance of the PyrTriTy3-2 or NecHaemTy3-4, HT looks more preferable for the explanation of the described case. It is known that transposable elements can alter gene expression since they carry their own regulatory sequences and insertions can be selectively advantageous. However, only those transposable elements, which were involved in regulation, evolve under strict selective pressure [44, 45].
While extremely rare, horizontal transfer seems to be quite common and recurrent in eukaryotes. An incomplete list of putative HT events includes: HT as a key event in the evolution of several fungal genes [46–48]; HT from fungi to rice weevil Sitophilus oryzae proposed for pectinase gene ; numerous HT events described for eukaryotic transposable elements [30–38]; as well as HTs of mitochondrial genes, for example, multiple angiosperm-angiosperm HTs of homing group I intron in the mitochondrial cox1 gene (for a review, see ); and a HT of the intron II and two adjacent exons of the mitochondrial nad1 gene from the flowering plants (angiosperms) to Gnetum (gymnosperms) .
The actual mechanisms of horizontal transfer for eukaryotic genes and transposable elements are still unknown since it is not possible to show experimentally how HT can occur. Parasites, symbionts, bacteria, or viruses all could be suggested as potential vectors for horizontal transfer. Moreover, based on an example of massive HT from a land plant donor to the basal angiosperm Amborella trichopoda, it has been demonstrated that direct plant-to-plant transfer can take a place . The associations between biotrophic fungi and their plant hosts are ubiquitous in nature and range from mutually beneficial to potentially fatal pathogenic interactions. Mycorrhiza refers to an association or symbiosis between plants and fungi that colonize the cortical tissue of plant roots. Ectomycorrhizal fungi are mostly basidiomycetes that grow between root cortical cells of many tree species whereas arbuscular mycorrhizal (AM) fungi belong to the order Glomales (Glomeromycota) and form highly branched structures called arbuscules, within root cortical cells of wide range of land plant species [53–55]. Both types of mycorrhiza represent intimate association and could provide suitable conditions for HT of transposable elements. AM-like mycorrhiza is widely distributed among mosses, ferns and lycophytes (for review, ).
Tcn1-like LTR retrotransposons were found in basidiomycota fungi and non-seed plants, including all investigated mosses and lycophytes from genus Selaginella. Such interkingdom distribution is not typical for chromodomain-containing LTR retrotransposons clades which are usually very specific for a particular taxonomic group and can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission. The close examination and comparison of evolutionary rates for LTR retrotransposons including representatives of Tcn1 and two other clades of LTR retrotransposons, and evolutionary rates estimated for putatively 'retained' genes from mosses and fungi suggests that a horizontal transmission of Tcn1-like LTR retrotransposons took place among fungi and mosses/lycophytes. However the evidence is not strong enough to discard other explanations; such as selective pressure coupled with vertical transmission of retrotransposons in genomes of non-seed plants and loss of these elements by other plants.
Genomic sequences screening, sequence and phylogenetic analysis
Fungal genomic sequences are available at: Fungal Genome Initiative ; The DOE Joint Genome Institute ; and The Sanger Institute . The source of individual genomes can be found in table represented in Additional file 5.
We used UniPro uGENE software  for LTR retrotransposons identification. The designed pipeline for Gypsy LTR retrotransposons identification and classification included: loading genomic sequence, translation of genomic sequence over six possible reading frames to amino acids, and subsequent search for homologous regions performed using "HMMER search" options of UniPro uGENE. The algorithm of HMMER search is based on profile hidden Markov models, which can perform amino acid sequence searches by use of an appropriate profile . For the analyses, we used a multiple alignment consensus sequence, which contains Gypsy LTR retrotransposon reverse transcriptase (RT) and partial integrase (Int) domains. The profile HMM, based on this consensus sequence, was built using UniPro uGENE software. An additional test for the presence of RT and partial Int domains was performed using BLAST (blastp) which also was incorporated in the designed pipeline. All BLAST analysis was essentially performed using sequence databases accessible from the National Center for Biotechnology Information . The classification of the newly identified elements was performed by a comparative analysis of their sequences. Newly identified elements and their accession numbers in public databases are listed in Additional file 6.
The whole nucleotide sequences of the transposable elements, if possible, were also extracted with the assistance of UniPro uGENE software. After localization of amino acid sequences obtained during HMMER search in the initial genomes in its nucleotide representation, the sequences were expanded up to 15 Kb and used for long terminal repeats (LTRs) search. The algorithm for repeats search, 'Repeat Find', is included to the UniPro uGENE as well as the visualization feature and 'ORF Find' option which were used to identify the putatively intact copies of LTR retrotransposons. Structural features of newly identified LTR retrotransposons can be found in Additional file 1.
Tcn1-like LTR retrotransposon search was carried out using BLAST (blastp and blastx). BLAST analysis was performed using sequence databases accessible from the National Center for Biotechnology Information (NCBI) server , The U.S. Department of Energy Joint Genome Institute , and Broad Institute of MIT and Harvard  as well as Phytozome, a tool for green plant comparative genomics . The described copy of SM-Tcn1 from spikemoss Selaginella moellendorffii (Lycopodiophyta) is located in scaffold_0 (1426925-1421008) of genomic sequence version 1.0 which is available at The U.S. Department of Energy Joint Genome Institute web-site . The whole sequence of SM-Tcn1 with annotations can be found in Additional file 3. Other websites used in the present study were: Repbase , NCBI conserved domain database and search service , ESTs from Porphyra yezoensis at Kazusa DNA Research Institute , Cyanidioschyzon merolae Genome Project , The Plant Genomics Consortium , The Institute for Genomic Research , Cassava and Leafy Spurge EST Project .
All multiple DNA alignments were performed by ClustalW  and edited manually in UniPro uGENE. Phylogenetic analyses were performed using the Neighbor-Joining (NJ) method in MEGA 4.0 program . Statistical support for the NJ tree was evaluated by bootstrapping (number of replications, 1000) . Evolutionary rates were estimated by standard methods . Poisson correction distances (d) were estimated from the equation d = -ln(1 - p), where p represents the proportion of different amino acids. The rate of amino acid substitution (r) was estimated by the standard equation r = d/2T, where T is the divergence time of the last common ancestor of the compared species. The estimated divergence times used were: Plants/Fungi, 1500 Myr and Basidiomycetes/Ascomycetes, 1200 Myr according to Hedges (2002) ; Homobasidiomycetes/Chytridiomycetes, 900 Myr, Sordaryomycetes/Eurotiomycetes, 540 Myr, and Sordaryomycetes/Dothideomycetes, 490 Myr according to Padovan et al. (2005) ; and Heterobasidiomycetes (or Tremellomycetes)/Homobasidiomycetes (Agaricomycetes and Dacrymycetes), 700 Myr according to Hibbett et al. (2007)  and Taylor et al. (2004) .
Species collection and total DNA isolation
Table 2 lists plant species, and Table 3 lists fungal species used in present study. The taxonomy of vascular non-seed plants (monilophytes and lycophytes) is given after Pryer et al. (2004) , Smith et al. (2006) , and Korall et al. (2007) . Plant species (monilophytes) were collected in nature. The detailed label data are available from the authors. The genomic DNA of lycophytes was provided by the Royal Botanic Gardens, Kew, London, UK . Genomic DNA was isolated from the leaves. Extraction was performed using the QIAGEN DNeasy Plant Mini Kit (QIAGEN). Isolated DNA was used directly in PCR amplifications.
Gypsy LTR retrotransposons PCR amplification and sequencing
Previously designed degenerate PCR primers for chromodomain-containing Gypsy LTR retrotransposons were used in present study: GyRT1 = 5'-MRNATGTGYGTNGAYTAYMG-3'  and ty3-A = 5'-AATTCGCTGCCGCTAAGATNARNADRTCRTC-3' , where M = A + C, Y = C + T, R = A + G, D = A + G + T and N = A + G + C + T. These primers were designed to amplify the most conserved part of the reverse transcriptase (RT) domain of LTR retrotransposons and were proved to be efficient [8, 24]. The expected length of PCR products was about 320 bp. PCR amplification with degenerate primers was performed using 0.1 μg of genomic DNA in 10-μl volume of 10 mM Tris-HCl (pH 8.9), 1 mM (NH4)2SO4, 4 mM MgCl2, 200 μM each of four dNTPs, 0.5 μM primers, and 2.5 units of Taq polymerase. After an initial denaturation step for 3 min at 94°C, the PCR reactions were subjected to 30 cycles of amplification consisting of 30 sec denaturation at 94°C, 42 sec annealing at 50°C, and 1 min extension at 72°C. PCR products were separated by agarose gel electrophoresis. The resulting PCR products were directly ligated into a pGEM vector using a pGEM-T-Easy cloning kit (Promega) for sequence determination.
Clones were amplified by PCR with M13 primers, and 40 ng of the product was used in a 10 μl cycle sequencing reaction with the ABI BigDye Terminator Kit on an ABI 310 Genetic Analyser (Applied Biosystems); or sequencing reactions were performed with Dye Terminator Cycle Sequencing Kit (Beckman Coulter) and analyzed on CEQ 8000 Genetic Analysis System. Sequences were deposited to GenBank under Acc. Numbers GQ443314-GQ443445 and AY959294-AY959313.
List of abbreviations used
long terminal repeat
short interspersed nuclear element
open reading frame
deoxyuridine triphosphatase domain
horizontal transfer or horizontal transmission
last common ancestor
hidden Markov models
million years ago.
Author thanks Dr. Mark L. Farman (Department of Plant Pathology, University of Kentucky, USA) for the helpful comments and Dr. David Thornbury (Department of Plant Pathology, University of Kentucky, USA) for his stylistic suggestions. This work was supported by the Russian Foundation for Basic Research (grant number RFBR 09-04-00360-a) and by state contract 10002-251/П-25/155-270/200404-082 and Siberian Branch of the Russian Academy of Sciences (project No. 10.4).
The sequence data for Laccaria bicolor, Trichoderma reesei, Trichoderma virens, Aspergillus niger, Stagonospora nodorum, Postia placenta and Sporobolomyces roseus were produced by the US Department of Energy Joint Genome Institute http://www.jgi.doe.gov/. Preliminary sequence data for Alternaria brassicicola were obtained from Genome Sequencing Center at Washington University Medical School http://genome.wustl.edu/.
The genomic DNA of lycophytes was provided by the Royal Botanic Gardens, Kew, London, UK http://www.kew.org. Monilophytes material was kindly provided by Dr. Alexander Shmakov (Altai State University, Barnaul, Russia) and Dr. Elena Korolyuk (Central Siberian Botanical Garden, Novosibirsk, Russia).
- Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O, Paux E, SanMiguel P, Schulman AH: A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007, 8: 973-982. 10.1038/nrg2165.PubMedView ArticleGoogle Scholar
- Marin I, Llorens C: Ty3/Gypsy retrotransposons: description of new Arabidopsis thaliana elements and evolutionary perspectives derived from comparative genomic data. Mol Biol Evol. 2000, 17: 1040-1049.PubMedView ArticleGoogle Scholar
- Gorinsek B, Gubensek F, Kordis D: Evolutionary genomics of chromoviruses in eukaryotes. Mol Biol Evol. 2004, 21: 781-798. 10.1093/molbev/msh057.PubMedView ArticleGoogle Scholar
- Paro R, Hogness DS: The Polycomb protein shares a homologous domain with a heterochromatin-associated protein of Drosophila. Proc Natl Acad Sci USA. 1991, 88: 263-267. 10.1073/pnas.88.1.263.PubMed CentralPubMedView ArticleGoogle Scholar
- Koonin EV, Zhou S, Lucchesi JC: The chromo superfamily: new members, duplication of the chromo domain and possible role in delivering transcription regulators to chromatin. Nucleic Acids Res. 1995, 23: 4229-4233. 10.1093/nar/23.21.4229.PubMed CentralPubMedView ArticleGoogle Scholar
- Brehm A, Tufteland KR, Aasland R, Becker PB: The many colours of chromodomains. Bioessays. 2004, 26: 133-140. 10.1002/bies.10392.PubMedView ArticleGoogle Scholar
- Gao X, Hou Y, Ebina H, Levin HL, Voytas DF: Chromodomains direct integration of retrotransposons to heterochromatin. Genome Res. 2008, 18: 359-369. 10.1101/gr.7146408.PubMed CentralPubMedView ArticleGoogle Scholar
- Novikova O, Mayorov V, Smyshlyaev G, Fursov M, Adkison L, Pisarenko O, Blinov A: Novel clades of chromodomain-containing Gypsy LTR retrotransposons from mosses (Bryophyta). Plant J. 2008, 56: 562-574. 10.1111/j.1365-313X.2008.03621.x.PubMedView ArticleGoogle Scholar
- Novikova O: Chromodomains and LTR retrotransposons in plants. Comm & Integr Biol. 2009, 2: 158-162.View ArticleGoogle Scholar
- Neuvéglise C, Feldmann H, Bon E, Gaillardin C, Casaregola S: Genomic evolution of the long terminal repeat retrotransposons in hemiascomycetous yeasts. Genome Res. 2002, 12: 930-943. 10.1101/gr.219202.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmid-Berger N, Schmid B, Barth G: Ylt1, a highly repetitive retrotransposon in the genome of the dimorphic fungus Yarrowia lipolytica. J Bacteriol. 1994, 176: 2477-2482.PubMed CentralPubMedGoogle Scholar
- Goodwin TJ, Poulter RT: Multiple LTR-retrotransposon families in the asexual yeast Candida albicans. Genome Res. 2000, 10: 174-191. 10.1101/gr.10.2.174.PubMedView ArticleGoogle Scholar
- Goodwin TJ, Poulter RT: The diversity of retrotransposons in the yeast Cryptococcus neoformans. Yeast. 2001, 18: 865-880. 10.1002/yea.733.PubMedView ArticleGoogle Scholar
- Novikova OS, Blinov AG: dUTPase-containing Metaviridae LTR retrotransposons from the genome of Phanerochaete chrysosporium (Fungi: Basidiomycota). Dokl Bioch Bioph. 2008, 420: 146-149. 10.1134/S1607672908030137.View ArticleGoogle Scholar
- Riccioni C, Rubini A, Belfiori B, Passeri V, Paolocci F, Arcioni S: Tmt1: the first LTR-retrotransposon from a Tuber spp. Curr Genet. 2008, 53: 23-34. 10.1007/s00294-007-0155-9.PubMedView ArticleGoogle Scholar
- Payne SL, Elder JH: The role of retroviral dUTPases in replication and virulence. Curr Protein Pept Sci. 2001, 2: 381-388. 10.2174/1389203013381008.PubMedView ArticleGoogle Scholar
- Varmus H, Brown P: Retroviruses. Mobile DNA. Edited by: Berg DE, Howe MM. 1989, Washington: American Society for Microbiology, 53-108.Google Scholar
- Wilhelm M, Wilhelm FX: Reverse transcription of retroviruses and LTR retrotransposons. Cell Mol Life Sci. 2001, 58: 1246-1262. 10.1007/PL00000937.PubMedView ArticleGoogle Scholar
- Levin HL: A novel mechanism of self-primed reverse transcription defines a new family of retroelements. Mol Cell Biol. 1995, 15: 3310-3317.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang W, Tanurdzic M, Luo M, Sisneros N, Kim HR, Weng JK, Kudrna D, Mueller C, Arumuganathan K, Carlson J, Chapple C, de Pamphilis C, Mandoli D, Tomkins J, Wing RA, Banks JA: Construction of a bacterial artificial chromosome library from the spikemoss Selaginella moellendorffii: a new resource for plant comparative genomics. BMC Plant Biol. 2005, 5: 10-10.1186/1471-2229-5-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith AR, Pryer KM, Schuettpelz E, Korall P, Schneider H, Wolf PG: A classification for extant ferns. Taxon. 2006, 55: 705-731.View ArticleGoogle Scholar
- Pryer KM, Schuettpelz E, Wolf PG, Schneider H, Smith AR, Cranfill R: Phylogeny and evolution of ferns (monilophytes) with a focus on the early leptosporangiate divergences. Am J Bot. 2004, 91: 1582-1598. 10.3732/ajb.91.10.1582.PubMedView ArticleGoogle Scholar
- Wikstrom N, Kenrick P: Evolution of Lycopodiaceae (Lycopsida): estimating divergence times from rbcL gene sequences by use of nonparametric rate smoothing. Mol Phyl Evol. 2001, 19: 177-186. 10.1006/mpev.2001.0936.View ArticleGoogle Scholar
- Friesen N, Brandes A, Heslop-Harrison JS: Diversity, origin, and distribution of retrotransposons (gypsy and copia) in conifers. Mol Biol Evol. 2001, 18: 1176-1188.PubMedView ArticleGoogle Scholar
- Wright DA, Voytas DF: Potential retroviruses in plants: Tat1 is related to a group of Arabidopsis thaliana Ty3/gypsy retrotransposons that encode envelope-like proteins. Genetics. 1998, 149: 703-715.PubMed CentralPubMedGoogle Scholar
- Vitte C, Panaud O: LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model. Cytogenet Genome Res. 2005, 110: 91-107. 10.1159/000084941.PubMedView ArticleGoogle Scholar
- Hedges SB: The origin and evolution of model organisms. Nat Rev Genet. 2002, 3: 838-849. 10.1038/nrg929.PubMedView ArticleGoogle Scholar
- Rensing SA, Fritzowsky D, Lang D, Reski R: Protein encoding genes in an ancient plant: analysis of codon usage, retained genes and splice sites in a moss, Physcomitrella patens. BMC Genomics. 2005, 6: 43-10.1186/1471-2164-6-43.PubMed CentralPubMedView ArticleGoogle Scholar
- Hartl DL, Lohe AR, Lozovskaya ER: Modern thoughts on an ancyent marinere: function, evolution, regulation. Annu Rev Genet. 1997, 31: 337-358. 10.1146/annurev.genet.31.1.337.PubMedView ArticleGoogle Scholar
- Jordan IK, Matyunina LV, McDonald JF: Evidence for the recent horizontal transfer of long terminal repeat retrotransposon. Proc Natl Acad Sci USA. 1999, 96: 12621-12625. 10.1073/pnas.96.22.12621.PubMed CentralPubMedView ArticleGoogle Scholar
- Novikova O, Sliwiñska E, Fet V, Settele J, Blinov A, Woyciechowski M: CR1 clade of non-LTR retrotransposons from Maculinea butterflies (Lepidoptera: Lycaenidae): evidence for recent horizontal transmission. BMC Evol Biol. 2007, 7: 93-10.1186/1471-2148-7-93.PubMed CentralPubMedView ArticleGoogle Scholar
- Roulin A, Piegu B, Wing RA, Panaud O: Evidence of multiple horizontal transfers of the long terminal repeat retrotransposon RIRE1 within the genus Oryza. Plant J. 2008, 53: 950-959. 10.1111/j.1365-313X.2007.03388.x.PubMedView ArticleGoogle Scholar
- Roulin A, Piegu B, Fortune PM, Sabot F, D'Hont A, Manicacci D, Panaud O: Whole genome surveys of rice, maize and sorghum reveal multiple horizontal transfers of the LTR-retrotransposon Route66 in Poaceae. BMC Evol Biol. 2009, 9: 58-10.1186/1471-2148-9-58.PubMed CentralPubMedView ArticleGoogle Scholar
- Kidwell MG: Horizontal transfer. Curr Opin Genet Dev. 1992, 2: 868-873. 10.1016/S0959-437X(05)80109-1.PubMedView ArticleGoogle Scholar
- Volff JN, Körting C, Schartl M: Multiple lineages of the non-LTR retrotransposon Rex1 with varying success in invading fish genomes. Mol Biol Evol. 2000, 17: 1673-1684.PubMedView ArticleGoogle Scholar
- Kordis D, Gubensek F: Unusual horizontal transfer of a long interspersed nuclear element between distant vertebrate classes. Proc Natl Acad Sci USA. 1998, 95: 10704-10709. 10.1073/pnas.95.18.10704.PubMed CentralPubMedView ArticleGoogle Scholar
- Novikova O, Fet V, Blinov A: Non-LTR retrotransposons in fungi. Funct Integr Genomics. 2009, 9: 27-42. 10.1007/s10142-008-0093-8.PubMedView ArticleGoogle Scholar
- Kaplan N, Darden T, Langley CH: Evolution and extinction of transposable elements in Mendelian populations. Genetics. 1985, 109: 459-480.PubMed CentralPubMedGoogle Scholar
- Lohe AR, Moriyama EN, Lidholm DA, Hartl DL: Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements. Mol Biol Evol. 1995, 12: 62-72.PubMedView ArticleGoogle Scholar
- Bowman JL, Floyd SK, Sakakibara K: Green genes-comparative genomics of the green branch of life. Cell. 2007, 129: 229-234. 10.1016/j.cell.2007.04.004.PubMedView ArticleGoogle Scholar
- Berbee ML, Taylor JW: Fungal Molecular Evolution: Gene Trees and Geologic Time. The Mycota: a comprehensive treatise on fungi as experimental systems for basic and applied research. Systematics and Evolution, Part B. Edited by: McLaughlin DJ, McLaughlin EG, Lemke PA. 2001, New York: Springer-Verlag, VII: 229-246.Google Scholar
- Kaneko I, Tanaka A, Tsuge T: REAL, an LTR retrotransposon from the plant pathogenic fungus Alternaria alternata. Mol Gen Genet. 2000, 263: 625-634. 10.1007/s004380051210.PubMedView ArticleGoogle Scholar
- Taylor EJ, Konstantinova P, Leigh F, Bates JA, Lee D: Gypsy-like retrotransposons in Pyrenophora: an abundant and informative class of molecular markers. Genome. 2004, 47: 519-525. 10.1139/g04-008.PubMedView ArticleGoogle Scholar
- Medstrand P, Landry JR, Mager DL: Long terminal repeats are used as alternative promoters for the endothelin B receptor and apolipoprotein C-I genes in humans. J Biol Chem. 2001, 276: 1896-1903. 10.1074/jbc.M006557200.PubMedView ArticleGoogle Scholar
- Ono R, Kobayashi S, Wagatsuma H, Aisaka K, Kohda T, Kaneko-Ishino T, Ishino F: A retrotransposon-derived gene, PEG10, is a novel imprinted gene located on human chromosome 7q21. Genomics. 2001, 73: 232-237. 10.1006/geno.2001.6494.PubMedView ArticleGoogle Scholar
- Wenzl P, Wong L, Kwang-won K, Jefferson RA: A functional screen identifies lateral transfer of beta-glucuronidase (gus) from bacteria to fungi. Mol Biol Evol. 2005, 22: 308-316. 10.1093/molbev/msi018.PubMedView ArticleGoogle Scholar
- Slot JC, Hibbett DS: Horizontal transfer of a nitrate assimilation gene cluster and ecological transitions in fungi: a phylogenetic study. PLoS ONE. 2007, 2: e1097-10.1371/journal.pone.0001097.PubMed CentralPubMedView ArticleGoogle Scholar
- Khaldi N, Collemare J, Lebrun MH, Wolfe KH: Evidence for horizontal transfer of a secondary metabolite gene cluster between fungi. Genome Biol. 2008, 9: R18-10.1186/gb-2008-9-1-r18.PubMed CentralPubMedView ArticleGoogle Scholar
- Shen Z, Denton M, Mutti N, Pappan K, Kanost MR, Reese JC, Reeck GR: Polygalacturonase from Sitophilus oryzae: possible horizontal transfer of a pectinase gene from fungi to weevils. J Insect Sci. 2003, 3: 24-PubMed CentralPubMedView ArticleGoogle Scholar
- Richardson AO, Palmer JD: Horizontal gene transfer in plants. J Exp Bot. 2007, 58: 1-9. 10.1093/jxb/erl148.PubMedView ArticleGoogle Scholar
- Won H, Renner SS: Horizontal gene transfer from flowering plants to Gnetum. Proc Natl Acad Sci USA. 2003, 100: 10824-10829. 10.1073/pnas.1833775100.PubMed CentralPubMedView ArticleGoogle Scholar
- Bergthorsson U, Richardson AO, Young GJ, Goertzen LR, Palmer JD: Massive horizontal transfer of mitochondrial genes from diverse land plant donors to the basal angiosperm Amborella. Proc Natl Acad Sci USA. 2004, 101: 17747-17752. 10.1073/pnas.0408336102.PubMed CentralPubMedView ArticleGoogle Scholar
- Albrecht C, Geurts R, Bisseling T: Legume nodulation and mycorrhizae formation; two extremes in host specificity meet. EMBO J. 1999, 18: 281-288. 10.1093/emboj/18.2.281.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang B, Qiu YL: Phylogenetic distribution and evolution of mycorrhizas in land plants. Mycorrhiza. 2006, 16: 299-363. 10.1007/s00572-005-0033-6.PubMedView ArticleGoogle Scholar
- Hibbett DS, Matheny PB: The relative ages of ectomycorrhizal mushrooms and their plant hosts estimated using Bayesian relaxed molecular clock analyses. BMC Biol. 2009, 7: 13-10.1186/1741-7007-7-13.PubMed CentralPubMedView ArticleGoogle Scholar
- Broad Institute. [http://www.broad.mit.edu]
- The DOE Joint Genome Institute. [http://www.jcvi.org/]
- The Sanger Institute. [http://www.sanger.ac.uk]
- UniPro uGENE software. [http://genome.unipro.ru]
- Eddy SR: A probabilistic model of local sequence alignment that simplifies statistical significance estimation. PLoS Comput Biol. 2008, 4: e1000069-10.1371/journal.pcbi.1000069.PubMed CentralPubMedView ArticleGoogle Scholar
- National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov]
- Phytozome, a tool for green plant comparative genomics. [http://www.phytozome.net]
- Repbase. [http://www.girinst.org]
- National Center for Biotechnology Information conserved domain database and search service. [http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml]
- ESTs from Porphyra yezoensis at Kazusa DNA Research Institute. [http://est.kazusa.or.jp]
- Cyanidioschyzon merolae Genome Project. [http://merolae.biol.s.u-tokyo.ac.jp]
- The Plant Genomics Consortium. [http://nypg.bio.nyu.edu/]
- The Institute for Genomic Research. [http://www.tigr.org]
- Cassava and Leafy Spurge EST Project. [http://titan.biotec.uiuc.edu/cassava]
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Nei M, Kumar S: Molecular evolution and phylogenetics. 2000, New York: Oxford University PressGoogle Scholar
- Padovan AC, Sanson GF, Brunstein A, Briones MR: Fungi evolution revisited: application of the penalized likelihood method to a Bayesian fungal phylogeny provides a new perspective on phylogenetic relationships and divergence dates of Ascomycota groups. J Mol Evol. 2005, 60: 726-735. 10.1007/s00239-004-0164-y.PubMedView ArticleGoogle Scholar
- Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, Huhndorf S, James T, Kirk PM, Lücking R, Thorsten Lumbsch H, Lutzoni F, Matheny PB, McLaughlin DJ, Powell MJ, Redhead S, Schoch CL, Spatafora JW, Stalpers JA, Vilgalys R, Aime MC, Aptroot A, Bauer R, Begerow D, Benny GL, Castlebury LA, Crous PW, Dai YC, Gams W, Geiser DM, Griffith GW, Gueidan C, Hawksworth DL, Hestmark G, Hosaka K, Humber RA, Hyde KD, Ironside JE, Kõljalg U, Kurtzman CP, Larsson KH, Lichtwardt R, Longcore J, Miadlikowska J, Miller A, Moncalvo JM, Mozley-Standridge S, Oberwinkler F, Parmasto E, Reeb V, Rogers JD, Roux C, Ryvarden L, Sampaio JP, Schüssler A, Sugiyama J, Thorn RG, Tibell L, Untereiner WA, Walker C, Wang Z, Weir A, Weiss M, White MM, Winka K, Yao YJ, Zhang N: A higher-level phylogenetic classification of the Fungi. Mycol Res. 2007, 111: 509-547. 10.1016/j.mycres.2007.03.004.PubMedView ArticleGoogle Scholar
- Taylor JW, Spatafora J, O'Donnell K, Lutzoni F, James T, Hibbett DS, Geiser D, Bruns TD, Blackwell M: The relationships of fungi. Assembling the tree of life. Edited by: Cracraft J, Donoghue MJ. 2004, New York: Oxford University Press, 171-196.Google Scholar
- Korall P, Conant DS, Metzgar JS, Schneider H, Pryer KM: A molecular phylogeny of scaly tree ferns (Cyatheaceae). Am J Bot. 2007, 94: 873-886. 10.3732/ajb.94.5.873.PubMedView ArticleGoogle Scholar
- The Royal Botanic Gardens, Kew, London, UK. [http://www.kew.org]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.