Evolution of GHF5 endoglucanase gene structure in plant-parasitic nematodes: no evidence for an early domain shuffling event
© Kyndt et al; licensee BioMed Central Ltd. 2008
Received: 17 October 2007
Accepted: 03 November 2008
Published: 03 November 2008
Endo-1,4-beta-glucanases or cellulases from the glycosyl hydrolase family 5 (GHF5) have been found in numerous bacteria and fungi, and recently also in higher eukaryotes, particularly in plant-parasitic nematodes (PPN). The origin of these genes has been attributed to horizontal gene transfer from bacteria, although there still is a lot of uncertainty about the origin and structure of the ancestral GHF5 PPN endoglucanase. It is not clear whether this ancestral endoglucanase consisted of the whole gene cassette, containing a catalytic domain and a carbohydrate-binding module (CBM, type 2 in PPN and bacteria) or only of the catalytic domain while the CBM2 was retrieved by domain shuffling later in evolution. Previous studies on the evolution of these genes have focused primarily on data of sedentary nematodes, while in this study, extra data from migratory nematodes were included.
Two new endoglucanases from the migratory nematodes Pratylenchus coffeae and Ditylenchus africanus were included in this study. The latter one is the first gene isolated from a PPN of a different superfamily (Sphaerularioidea); all previously known nematode endoglucanases belong to the superfamily Tylenchoidea (order Rhabditida). Phylogenetic analyses were conducted with the PPN GHF5 endoglucanases and homologous endoglucanases from bacterial and other eukaryotic lineages such as beetles, fungi and plants. No statistical incongruence between the phylogenetic trees deduced from the catalytic domain and the CBM2 was found, which could suggest that both domains have evolved together. Furthermore, based on gene structure data, we inferred a model for the evolution of the GHF5 endoglucanase gene structure in plant-parasitic nematodes. Our data confirm a close relationship between Pratylenchus spp. and the root knot nematodes, while some Radopholus similis endoglucanases are more similar to cyst nematode genes.
We conclude that the ancestral PPN GHF5 endoglucanase gene most probably consisted of the whole gene cassette, i.e. the GHF5 catalytic domain and the CBM2, rather than that it evolved by domain shuffling. Our evolutionary model for the gene structure in PPN GHF5 endoglucanases implies the occurrence of an early duplication event, and more recent gene duplications at genus or species level.
Endo-1,4-beta-glucanases (EC 184.108.40.206), or endoglucanases, can degrade the beta-1,4-linkages of cellulose, the most abundant component of plant cell walls. The endoglucanases are classified in different glycosyl hydrolase families (GHF) on the basis of sequence similarity and hydrophobic cluster analysis . Animal endoglucanases belong to three structurally and presumably phylogenetically unrelated families: GHF5, GHF9 and GHF45 . GHF5 genes are found in various plant-parasitic bacteria and nematodes [3–6] from the order Rhabditida (infra-order Tylenchomorpha), which is part of one of the three evolutionary independent plant-parasitic nematode clades . In the other two clades (Dorylaimida and Triplonchida), no GHF5 endoglucanases have been found so far. Next to plant-parasitic nematodes (PPN), some protists, plants and fungi, as well as two beetle species [8, 9] have GHF5 endoglucanases. The domain structure of GHF5 PPN endoglucanases shows remarkable similarity to particular bacterial genes that consist of a GHF5 catalytic domain, a linker and a carbohydrate-binding module of family 2 (CBM2) [10, 11]. Interestingly, some GHF5 PPN endoglucanases lack the CBM. Other bacterial endoglucanases as well as those from plants and fungi contain different types of carbohydrate-binding modules, while currently known GHF5 endoglucanases of protists and beetles comprise no CBM.
Horizontal gene transfer (HGT) from bacteria has been proposed to explain the origin of GHF5 endoglucanases in nematodes . However, although HGT was also suggested to lie at the origin of GHF9 cellulases, this hypothesis was rejected based on recent evidence showing that there was a GHF9 cellulase gene in the last common ancestor of bilaterian animals . At present, no hard evidence for the HGT hypothesis of GHF5 endoglucanases has been found. An evolutionary scheme proposed for GHF5 genes of sedentary nematodes by Ledger et al.  hypothesized that after the possible HGT event, the endoglucanase genes must have been duplicated several times, and in some cases a sequential loss of the linker and CBM occurred. This model assumes that the domain structure of the ancestral PPN GHF5 endoglucanase consisted of a catalytic domain, a linker and a CBM, although the structure of this ancestral gene has never been firmly established. A possible HGT event from bacteria could have involved the whole cassette as well as only the CBM2. The fact that the GHF5 domain has also been found in other eukaryotic lineages can be seen as evidence for its presence in the ancestral eukaryote, and a subsequent combination with a horizontally-acquired CBM2 through 'early' (i.e. in one of the ancestral nematodes) domain shuffling could have resulted in the presence of proteins with both domains in PPN. For bacterial genes, domain shuffling between GHF5 and CBM2 domains has already been proven .
Domain shuffling  is a special form of exon shuffling, in which new combinations of exons (encompassing a complete protein domain in the case of 'domain shuffling') from unrelated genes can be created by recombination within the intervening intron sequences, yielding rearranged genes with altered function . Since the discovery of introns there has been considerable debate concerning their origin and possible involvement in exon/domain shuffling events.
The original introns-early concept states that (nearly) all introns have been inherited by eukaryotic genes from prokaryotic ancestors [16–19]. The difference in exon/intron-structure among homologous genes is then considered to be mostly due to differential intron loss. Observed relationships between intron locations and protein structures (module boundaries) are used as evidence for the ancient origin of introns and for their role in exon/domain shuffling .
Cavalier-Smith strongly disagrees with the introns-early view . On the basis of phylogenetic considerations and the unique requirements of the eukaryotic nucleus, he suggests that the original genes were uninterrupted. He and other proponents of the introns-late theory argue that the introns observed in extant genes occur at far too many positions to have been present all together in a common ancestral gene . In addition, they argue that the introns-early idea is hard to reconcile with the absence of spliceosomal introns and the spliceosome itself in prokaryotes (explained by 'introns-early' debaters as the result of genome streamlining in prokaryotes). The introns-late concept states that introns are an invention of eukaryotes and new introns have been continuously emerging during eukaryotic evolution . This theory is supported by the discovery of a conserved sequence at exon-intron boundaries, called the protosplice-site (MAG/R) , which is believed to be a specific target site for newly gained introns, or a preferential site for intron fixation.
A synthetic theory of intron evolution has been put forward by de Souza . This intermediate position proposes that a subset of present-day introns is ancient, but in addition, allows the emergence of some new introns. The so-called 'many-introns-early-in-evolution' theory places the origin of numerous spliceosomal introns during (and maybe playing a role in) eukaryotic radiation . Recent studies provide high support for this theory [27–29]. The lack of conserved intron positions in ancient paralogues was suggested to demonstrate that the invasion of group II self-splicing elements (which were to become spliceosomal introns in the future) occurred during gene duplication, most probably in the earliest phase of eukaryotic evolution .
In order to determine the domain structure of the ancestral PPN GHF5 endoglucanase gene present in the PPN of the phylum Tylenchomorpha, the evolution of the GHF5 catalytic domain and the CBM2 cassette was investigated through phylogenetic analyses including sequences from fungi, plants, protista and bacteria.
To study the subsequent evolution of the GHF5 endoglucanase gene family within nematodes, intron positions were used as phylogenetic characters and a model was developed for their evolution.
Until now, evolutionary models for PPN GHF5 endoglucanases focused on genes from cyst and root knot nematodes and one migratory nematode, all belonging to the Tylenchoidea superfamily. In order to have a broader overview of the endoglucanase evolution in the infra-order Tylenchomorpha, the gene structure of six additional genes was incorporated in our study. These include the four Radopholus similis genes that were recently described  and two novel GHF5 endoglucanases from the migratory endoparasitic nematodes Pratylenchus coffeae and Ditylenchus africanus. The latter gene is the first GHF5 endoglucanase isolated from the superfamily Sphaerularioidea .
Results and discussion
Glycosyl hydrolase family 5 (GHF5) endo-1,4-beta-glucanase genes have been identified in various plant-parasitic nematodes (PPN). These genes have different structures: all have a signal peptide and a catalytic domain, some have an additional linker and carbohydrate-binding module (CBM, type 2 or CBM2) and others only have a linker but no CBM2. This modular domain structure could have arisen by exon or domain shuffling, a process by which a chimeric protein is created by domains or segments of two different genes through recombination in intron regions. In bacteria, all endo-1,4-beta-glucanases are supposed to originate from a few progenitor sequences by mutation and domain shuffling  and Quillet et al. has proven the validity of this hypothesis for an endo-1,4-beta-glucanase of the bacterium Myxococcus xanthus . The authors suggest that the catalytic domain and the CBM of this endoglucanase were acquired by independent horizontal transfers from different soil bacteria, and that these domains were assembled by exon shuffling to form a new gene.
While the ancient origin of the PPN GHF5 endoglucanases (horizontal gene transfer or not ) still is a matter of debate, in this study no support for an early exon/domain shuffling event after the independent gain of a catalytic domain (e.g. through evolution from an eukaryotic ancestor) and a CBM2 (e.g. through HGT from bacteria) in an ancestral PPN species could be found in the currently available sequence data.
Endoglucanase genes of Ditylenchus africanus and Pratylenchus coffeae
Endo-1,4-beta-glucanase fragments were amplified by a degenerate PCR on genomic DNA of D. africanus and P. coffeae. This resulted in a fragment of 605 bp for D. africanus and a fragment of 443 bp for P. coffeae. Sequence analysis confirmed them to be endoglucanase fragments and the genes were named Da-eng1 and Pc-eng1 respectively. Full length sequences of the genes were obtained by genome walking. The length of these sequences from start till stop codon is 2007 bp for Da-eng1 and 1527 bp for Pc-eng1. The corresponding coding sequences were amplified from a cDNA pool and have respective lengths of 1395 bp and 1377 bp. Da-eng1 contains 5 introns whereas Pc-eng1 only has 3.
The putative endoglucanase proteins Da-ENG1 and Pc-ENG1 have an estimated molecular weight of 49 and 47 kDa respectively. Both have a predicted N-terminal signal peptide for secretion. A conserved domain search confirmed that the proteins are endoglucanases belonging to the glycosyl hydrolase family 5 (GHF5, pfam00150). In addition to the catalytic domain, both proteins have a linker and a carbohydrate-binding module (CBM) of family 2 (pfam00553).
In addition to the already known endoglucanase genes from PPN and the newly identified sequences in this study, we searched for the most homologous endoglucanase genes from other major lineages. Based on their sequence similarity to the PPN endoglucanases, one or more representative sequences were selected per lineage for further analyses. These included GHF5 endoglucanase sequences from the kingdom of Plantae, Animalia, Fungi, Protista and Bacteria (See additional file 1). Some of these genes consist of both a catalytic domain and a CBM, while others lack the CBM.
Separate alignments were made for the catalytic domain and the carbohydrate-binding module (CBM) of the genes. Since not all GHF5 endoglucanase genes include a CBM, the alignment and hence also the phylogenetic tree for the CBM contained less taxa. For both domains, alignment of the sequences was performed at both DNA and protein sequence level.
Alignments showed that plant sequences are not sufficiently similar to PPN GHF5 endoglucanases to be useful for phylogenetic analyses. Preliminary unrooted phylogenetic tree construction with the other endoglucanases revealed fungal sequences to be most divergent from PPN GHF5 endoglucanases. Therefore they were chosen as outgroup to construct rooted trees. All alignments were analyzed using Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian methods.
No data from linker sequences have been included in this study (except for the presence of intron 22 in the linker of Mi-eng1, see further) because linkers show very little large-scale sequence homology. At DNA level, only linkers from very closely related genes show some degree of homology. This is not surprising as the same observation has been made for linkers from bacterial endoglucanases . In general, nematode linkers, as well as bacterial linkers, are short and mainly composed of small or tiny amino acids like glycine, serine, threonine, alanine, proline, and lysine.
All evolutionary trees, whether based on DNA or protein data, showed similar clustering, independent of the method used.
Within the PPN, Mi-ENG2, Da-ENG1 and Hg-ENG6 are at the base of the protein family but do not show a well-supported clustering. This grouping is correlated with their aberrant exon/intron structure (see further). The divergence of Da-ENG1 at sequence and structural level is not surprising since Ditylenchus africanus is the only species in our study belonging to the Sphaerularioidea and not to the Tylenchoidea.
A well-confirmed cluster groups the GHF5 endoglucanases of Pratylenchus coffeae and P. penetrans with those of the root knot nematodes (genus Meloidogyne). Cyst nematode (Heterodera and Globodera) endoglucanases are grouped together with some endoglucanases from Radopholus similis (classified in the Pratylenchidae family), Rs-ENG1B and Rs-ENG3. Although the positions of Rs-ENG1A and Rs-ENG2 are unresolved, these observations partially confirm recent phylogenetic trees derived from 28SrRNA data  revealing a closer relationship of Radopholus similis to Heteroderidae than to Pratylenchidae. This distinction was also observed on a morphological level by Luc , who found that the genus Radopholus can be distinguished from other genera in Pratylenchidae by its strong secondary sexual dimorphism and the distinctive lip pattern.
Paralogy between groups of proteins is well-confirmed by high posterior probabilities. For instance, Rs-ENG3 and Hg-ENG5 are paralogous to each other and to two orthologues from Globodera rostochiensis (Gr-ENG3 and Gr-ENG4). Also, Rs-ENG1B is paralogous to some intra-species or intra-genus duplicates from Heterodera and Globodera spp.
Investigation of gene/domain structure and intron properties
The ORF (Open Reading Frame) of PPN GHF5 endoglucanase genes is interrupted by maximum 8 introns (Hg-eng1) (Figure 4) with an average of 5.28 introns. The average length of each intron is 132 bp, making a total of 697 bp of intron per gene, which is smaller than the (already small) average intron size observed in the nematode lineage (467 bp per intron ).
The GC percentage per intron varies from 5% (intron 10 of Mi-eng5) to 58% (intron 2 of Gr-eng1), with an average for all introns (n = 110) of 31%. Apparently there are differences in GC percentage according to the species. For instance, the average GC% of all introns of M. incognita endoglucanase genes (n = 23) is 22%, whereas for R. similis introns (n = 13) it is 42%. This is probably related to the overall GC% of the coding sequences of a certain species. M. incognita coding sequences generally have a relatively low GC content (37%) , whereas R. similis coding sequences have a high GC content (54%) .
67% of the analysed introns are in phase 0, while only 17% correspond to phase 1 or phase 2. Phase 0 is highly over-represented in PPN endoglucanases, as also observed by Ledger et al. , with twice the expected value for an addition of intron with equal probability for the three codon sites. This is also considerably higher than the observed prevalence in C. elegans (47% ). 49% (54/110) of the exons are symmetrical. As a consequence of this high over-representation of phase 0 introns, symmetrical exons of type 0-0 (44%, 48/110) are largely in excess compared to symmetrical exons of type 1-1 (4%, 4/110) or type 2-2 (2%, 2/110). The observed excess of phase 0 introns can be due to the preferential retention of these introns by natural selection  or as a consequence of the burst of domain shuffling that has been shown to coincide with the big bang of metazoan radiation in the Cambrian period [15, 38].
Exon-intron boundaries were compared for all introns observed in the available PPN GHF5 endoglucanase gene sequences. The great majority of the detected introns are of the GU-AG type. Only three intron positions (intron 13, 17a and 17b) show the alternative, but rare (0.6% in C. elegans) GC-AG splice site in a minority of the homologues. This non-classical splice site was already reported by Yan et al. .
Protosplice site tendency in PPN GHF5 endoglucanases is moderately high as 46% (11/24) of the introns show a consensus splice site that is completely consistent with the protosplice site (MAG/R), pointing to a probable recent origin of these introns.
Even though many clusters of closely spaced introns were detected in the aligned sequence, no evidence was found in our study, or in other studies, that proves their possible homology. Although intron sliding cannot be completely ruled out, the influence of this process is considered to be negligible (except for sliding by one base), and intron position diversity is indicated to be primarily arising by the gain or loss of introns during eukaryotic evolution . Therefore, a pair of introns was required to occur in the exact same position in the aligned sequence of the orthologous/paralogous genes to be considered homologous.
Noteworthy, there is a strong conservation of the location and phase of introns within the PPN GHF5 gene family . This conservation is evident among orthologous genes, but also among paralogues. Some intron positions are very well conserved (intron 2, 10, 13, 16, 17b, 21 and 23) while others only occur in one gene (intron 3, 6, 7, 8, 11, 12, 15, 19, 22).
Since spliceosomal introns diverge so rapidly that sequence similarity indicative of homology quickly vanishes, significant intron sequence conservation can only be expected in genes resulting from very recent duplications. Sequence comparisons between introns that are conserved in the majority of the dataset indeed showed similarities between paralogous genes. For instance, introns of Gr-eng1 and Gr-eng2 as well as Gr-eng3 and Gr-eng4 reveal remarkable similarity. This is also the case for introns in Hg-eng2, Hg-eng3 and Hg-eng4. The conserved introns of the Meloidogyne endoglucanase genes are also very similar, confirming the homology of these conserved introns and the paralogy of the genes.
As could be expected based on the low number of informative sites (14), Wagner Parsimony of the 1/0-data for intron presence/absence shows very low resolution (See additional file 2). Only some paralogous and orthologous genes with very similar exon/intron structure are grouped together with relatively high bootstrap support. For instance, Mi-eng2 and Hg-eng6 form a well-confirmed cluster. This is also the case for some groups of paralogues and orthologues from H. glycines, G. rostochiensis and R. similis. Generally, resolved branches show similar clustering as catalytic domain and CBM-based trees implying that intron presence data reflect the evolutionary tract of PPN GHF5 genes.
Objections against an early domain shuffling event
With the currently available endoglucanase sequences, no evidence for the domain shuffling hypothesis could be found. No statistical incongruence between the evolutionary history of the GHF5 catalytic domain and the CBM2 in the whole dataset as well as in PPN only (Table 1 and Fig. 3) was detected. Some additional trends in the gene structure data also decline the domain shuffling idea.
Gilbert  proposed that correlation between intron locations and structural elements implies a role for introns in the evolution of proteins. The intronic architecture of the GHF5 endoglucanase genes was mapped onto a multiple protein sequence alignment in the context of the overall secondary structure of an archetypal GHF5 endoglucanase (See additional file 3), but there appears to be no general correlation.
However, some conserved intron positions are specifically located at the border of the catalytic domain and could therefore be implicated in domain shuffling. Introns 20 and 21 are found in the transition zone between the catalytic domain and the linker region in many genes. Intron 1 is located after the signal peptide, and before the start of the catalytic domain, but is only present in two genes. Successful domain shuffling requires that the domains are bordered by introns that are of the same phase, that is, that the domain is symmetrical , because shuffling of asymmetrical exons/domains will result in a frame-shift. The only possible symmetric, and hence movable, domain would be a 1-1 domain bordered by intron 1 and intron 21. Although 1-1 symmetrical domains are suggested to be frequently associated with domain shuffling events , the low conservation of intron 1 among the investigated gene structures (it is only present in 2 genes) as well as its well-conserved protosplice site point to a quite recent origin, and therefore make its involvement in a possible domain shuffling event very unlikely.
Finally, exon or domain shuffling is more likely to involve large introns as the frequency of recombination is known to be proportional to DNA length. Although we cannot exclude the possibility that larger introns were present in an ancestral GHF5 PPN endoglucanase, introns in the present-day nematode genes are extremely small (also for introns 20 and 21) and are therefore less likely to be involved in recombination.
A model for the evolution of introns in the PPN GHF5 endoglucanases
The evolutionary tree based on the catalytic domains shows the genes with the same structure grouped in the same terminal branches, although the length and sequence of introns have diverged.
Our working hypothesis proposes that the divergence of the gene structure of the PPN GHF5 gene family is associated with the expansion of the number of members from an ancient or early eukaryotic ancestral gene. The possible congruent evolutionary tract of modules, as observed in this study (Table 1 and Fig. 3), suggests that the ancestral PPN GHF5 gene probably contained a CBM2. Bacterial endoglucanases show significant sequence homology in both domains and if the HGT hypothesis is true then the analyzed GHF5 endoglucanases likely evolved from a bacterial orthologue containing both domains. Present-day PPN GHF5 endoglucanase genes exhibit a large diversity of exon/intron structures, but protein sequences in the catalytic domain are relatively well-conserved. Generally, PPN GHF5 endoglucanases contain a rather high number of introns. It is known that paralogues undergo functional diversification, which tends to be accompanied by a weakened selection and the acceleration of sequence evolution . Babenko et al.  observed two clear trends in the evolution of paralogous gene families: the evolution of their gene structure is notably dynamic, and this evolution involves more intron gains than losses. In addition, Logsdon et al. noticed that nematodes have a particularly high rate of intron turnover compared to other animals . These tendencies are also recognizable in the PPN GHF5 gene family. Based on the recent indications that introns emerged during the earliest phases of eukaryotic evolution  our model starts from an intronless ancestral PPN GHF5 endoglucanase gene.
The exon/intron structure of Hg-eng6 and Mi-eng2 and their sequence divergence from the other genes indicate that they have probably originated from a different type of GHF5 gene (type A). This implies that the ancestral PPN GHF5 endoglucanase must have been duplicated already early in the Tylenchomorpha lineage. This early duplication event was also suggested by Ledger et al. . During or shortly after this duplication one of the copies (type A) must have lost its linker and CBM2. All other genes are proposed to have evolved from the second copy (type B) containing linker and CBM2. Selection might have played a role in the preferential conservation of this type, possibly due to a positive effect of the CBM2. Whether intron 23 has been gained before, during or shortly after the duplication cannot be evaluated based on the available sequence data.
In the ancestral nematode of the lineage leading to the cyst nematodes and Radopholus similis, a second gene duplication took place. These nematodes all contain two divergent types of GHF5 endoglucanases within their genomes.
At intra-species level, the conserved gene structure of certain paralogues and the observation of some sequence conservation in their introns indicate a history of recent duplication events in the present-day genomes of plant-parasitic nematodes.
Ledger et al.  also proposed an evolutionary scheme for the evolution of GHF5 endoglucanases in plant-parasitic nematodes. Their model proposed two ancestral copies of the GHF5 endoglucanase gene: one with and one without CBM, which is in agreement with our data. However, the model of Ledger et al.  does not take duplication events in ancestral taxa (e.g. as we have detected in the ancestor of Radopholus similis and the cyst nematodes) into consideration and is therefore confusing in regard of gene and species evolution.
In some cases, intra-species duplicates lost their CBM and occasionally also their linker region. Generally, the observed variation in exon/intron structures between homologous and paralogous genes implies that introns have been gained and lost along all major lineages of the tree and at different time points during evolution. Concordant with the fact that the conservation of protosplice sites strongly supports the recent origin of an intron , recently gained introns in our model (intron 5, 6, 12, 15, 19 and 22) show a stronger conservation of the protosplice site (MAG|R) than older introns.
Our model implies some special evolutionary gene structure processes for which little information is available and that are therefore difficult to prove.
The phylogenetic hallmark of intron sliding is a distribution in which one intron is nested within the distribution of another . We propose the occurrence of intron sliding for two introns in our model, namely sliding of intron 17a to intron 17b and intron 2 to intron 3, and this hypothesis is based on their close vicinity (1 nucleotide between 17a and 17b and 3 nucleotides between 2 and 3) and the apparent nested distribution. While the nested clustering theory of Stoltzfus et al.  pertains for these introns and they are located in very close vicinity, the fast evolution of intron sequences precludes the ability of finding significant sequence similarity at intron-level. In literature, the few established reports of intron sliding all involve very recent events [47–50]. Alternative mechanisms like separate gain and loss events might equally likely have been involved in the observed pattern of intron distribution.
The genes with no or only one intron probably have lost their introns by reverse transcription and subsequent homologous reinsertion of the synthesized cDNA, affecting either the whole gene or only a part of it. Hg-eng5 is an example of a gene where this process might have happened. We propose a similar mechanism for Rs-eng1A and Rs-eng2, involving a reverse-transcribed mRNA containing only a catalytic domain, but in this case the recombination must have been followed by a duplication event. The loss of introns through reverse transcription of mRNAs followed by the recombination of the synthesized cDNA in the genome [51–53] is hard to prove and an alternative gene duplication event cannot be ruled out because the genomes of plant-parasitic nematodes probably contain more GHF5 endoglucanase genes than currently available in databases. If a gene duplication would have given rise to these genes, then a closely related paralogue, generated in the same duplication event, should be found in the genome . The impossibility to establish a clear relation between the intronless gene Hg-eng5 and any other gene from the same species suggests that reverse transcription followed by a homologous reinsertion is most probable. For Rs-eng1A and Rs-eng2, the reverse transcription event is supposed to have taken place earlier in evolution, involving only the 5' end, and should have been followed by a duplication event, the loss of linker and CBM and a gain of intron 21 in one of the paralogues.
This study integrates previously known sequence information from mainly sedentary plant-parasitic nematodes and new data from different migratory nematode species. With the available data from PPN GHF5 endoglucanases and other major prokaryotic and eukaryotic lineages no evidence for an early domain shuffling event between the two domains (GHF5 catalytic domain and CBM2) present in these genes was found, suggesting parallel evolutionary tracts. Moreover, the phylogenetic data confirmed that some endoglucanases of the migratory species Radopholus similis are more related to genes from the Heteroderidae, while Pratylenchus spp. show a close relationship to the root knot nematodes.
Our evolutionary model for the PPN GHF5 gene family proposes an intronless ancestral GHF5 endoglucanase gene in plant-parasitic nematodes, consisting of a catalytic domain and a CBM2. Early in the evolution of plant-parasitic nematodes, the ancestral gene must have been duplicated, and numerous introns were gained. Later in evolution, some introns were lost, gained or shifted in certain lineages, and additional duplication events frequently took place. Our data reveal that loss of the CBM2 and in some cases also the linker region almost certainly occurred independently in different lineages.
While our data show that exon shuffling most probably did not lie at the origin of the emergence of the PPN GHF5 endoglucanase domain structure, more recent exon/domain shuffling events at lower taxonomic levels cannot be ruled out.
In this study the HGT hypothesis can neither be accepted nor rejected. Additional sequence data from other major eukaryotic lineages are awaited before a thorough and sound investigation of this problem can be realized. However, evidence was provided that it is very unlikely that an ancestral eukaryotic GHF5 catalytic domain was joined with a horizontally-acquired CBM2 in the early evolution of the PPN GHF5 domain structure. Therefore, HGT of the whole structural cassette (GHF5-linker-CBM2) remains one of the most likely explanations for its presence in both bacteria and PPN.
Pratylenchus coffeae and Ditylenchus africanus were maintained on carrot disks in small petri dishes at a constant temperature of 25°C . Nematodes were collected by rinsing the petri dishes with sterile demineralized water approximately 5 weeks (for D. africanus) or 8 weeks (for P. coffeae) after inoculation of the carrot disks. Collected nematodes of mixed stages were stored as a pellet at -20°C for DNA or RNA extraction.
Cloning of endoglucanase genes
Genomic DNA was isolated from the nematodes according to Bolla et al. . RNA extraction was done with TRIzol (Invitrogen) according to the manufacturer's instructions. First strand cDNA was synthesized using 1 μg RNA, 4 mM dNTP's, 0.5 μM oligodT primer (T25VN), 10 mM DTT, 50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2 and 200 U Superscript II Reverse Transcriptase (Invitrogen). This reaction mixture was incubated for 2 h at 42°C.
A polymerase chain reaction (PCR) using degenerate primers ENG1 (TAYGTIATHGTIGAYTGGCA) and ENG2 (GTICCRTAYTCIGTIACRAA)  was performed on genomic DNA of P. coffeae and D. africanus. The reaction mixture contained 150 ng DNA, 0.5 μM of each primer, 4 mM each of dNTPs, 1.5 mM MgCl2, 20 mM Tris-HCl (pH8.3), 50 mM KCl and 1 U of Taq DNA polymerase (Invitrogen). The PCR conditions were as follows: 2' at 94°C followed by 35 cycles of 1' at 94°C, 1' at 52°C and 1' at 72°C. Resulting fragments were separated on a 0.5× TAE 1.5% agarose gel, excised and purified by the QIAquick Gel Extraction kit (Qiagen). Purified fragments were ligated into pGEM-T (Promega) and transformed in E. coli DH5α cells (Invitrogen). A colony PCR under standard conditions revealed the presence or absence of inserts. Plasmids of positive colonies were isolated using the Nucleobond AX kit (Machery-Nagel) and the insert was sequenced at the VIB Genetic Service Facility (VIB-GSF, Antwerp, Belgium).
To obtain the 5' and 3' genomic sequences of the endoglucanase fragments, the Genome Walker Universal kit (Clontech) was used according to the manufacturer's instructions. Primers used for the P. coffeae fragment were Pc-eng1-up1 (TTGGAGTCAATGGCCCGGATGG), Pc-eng1-up2 (GATGTTCGGATTGGAGCCATATTGC) to walk upstream en Pc-eng1-down1 (TCAAATCTCTGCTACACTCTCCAC) and Pc-eng1-down2 (CAAACAATCCCTCAGGGATAAGGC) to walk downstream. For the D. africanus fragment, Da-eng1-up1 (CCAGCGTATTTTTGTGCCATTTGG), Da-eng1-up2 (CTCTCAACCAATATCTTTGACTTAACAGC), Da-eng1-down1 (AATAGTCGCGTTTAGGACGTGGACGTC) and Da-eng1-down2 (TCCACTTCTACGCCGGCACTCAC) were used to walk up- and downstream respectively. The longest obtained fragments were cloned and sequenced as described above. Corresponding cDNA sequences of the endoglucanase genes were amplified by PCR using gene-specific primers Da-eng1-start (ATGAAATTCTTCGCCAGCCTCG) and Da-eng1-stop (TCAGCAATATCCATAAGACACAACGC) on the cDNA pool of D. africanus and Pc-eng1-start (ATGGCATTCACTTTGCTTTCC) and Pc-eng1-stop (TGTGAGCGCCACTGCCTGCTAA) on P. coffeae cDNA. These cDNA sequences were cloned and sequenced as described above. Obtained sequences were submitted to the GenBank database (Da-eng1: EU180235; Pc-eng1: EU176871).
Putative protein sequences were obtained by translating the cDNA sequences using the EMBOSS program "Transeq" . The molecular weight of the proteins was estimated by the "pI/Mw" tool of the Expasy server http://us.expasy.org. Signal peptides were predicted using SignalP 3.0 .
Sequence searches and alignments
Using BlastX, a GenBank search for all available endoglucanase gene sequences from nematode species was performed on May 21, 2007. For comparison, the most homologous endoglucanases from other eukaryotic and prokaryotic lineages were included in the dataset (See additional file 1).
Protein, DNA and mRNA sequences were downloaded. Domains were identified using the SMART tool , which is based on the PFAM database. Sequences were aligned using ClustalW  and then manually adjusted in BioEdit .
Exon-intron structure analyses
The presence, location and phase of introns were evaluated for each open reading frame (ORF) by comparing DNA and mRNA sequences. Intron positions were mapped onto the sequence alignments. Intron phase was named 0 when the intron was located in-between two consecutive codons; 1 if an intron was located between the first and second codon nucleotides and 2 if an intron was found between the second and third codon nucleotides. The intron positions were mapped onto the secondary structure (See additional file 3). Positions of alpha helices and beta strands in the catalytic domain were deduced from an alignment of all nematode endoglucanases with an endoglucanase from Erwinia chrysanthemi (1EGZC), which has a known 3D structure. For the carbohydrate-binding module, the location of the beta strands was inferred by aligning the sequences to a CBM from an endoglucanase from Cellulomonas fimi (1EXG). GC percentages of the introns were calculated with the internet tool GC calculator http://www.genomicsplace.com/gc_calc.html.
Maximum Parsimony (MP) analyses were executed using PAUP* v4.0b10 , on the informative characters of the mRNA and protein alignments using the heuristic search option with random sequence addition (100 random replications) and TBR branch-swapping. Support for the different clusters was evaluated by bootstrap analysis (100 replicates).
Intron presence was coded by binary scoring (1/0) for presence/absence. Introns were only considered to be homologous when they appeared at exactly the same position and the same phase in the orthologous/paralogous genes. Although Dollo Parsimony is frequently used in intron evolution , recent evidence shows that this method tends to underestimate the number of introns in ancestral nodes and to overestimate the number of gains in branches leading to extant species  since it is based on the assumption that intron presence originated only once in the tree. Considering the fact that no consensus has been reached about the dominance of intron gain or loss during eukaryotic evolution, our 1/0-data were analysed using Wagner Parsimony .
Congruence between the obtained Maximum Parsimony trees based on catalytic domain and CBM was tested using the partition homogeneity test (HTF) as designated by Johnson & Soltis , which is analogous to the ILD (incongruence length difference) test of Farris et al.  as implemented in PAUP* v4.0b10 using the heuristic search option with random sequence addition (100 random replications) and TBR branch-swapping.
Hierarchical Likelihood Ratio Tests were performed on all datasets using Modeltest 3.5  to determine the best-fitting evolutionary model. Estimated model parameters were applied in a Maximum Likelihood (ML) analysis in PAUP* v4.0b10 using the heuristic search option with random sequence addition, for DNA sequences. Tree-puzzle 5.2  was used for ML analysis of protein sequences. Support for the different clusters was evaluated by bootstrap analysis (100 replicates).
In order to compare trees obtained from ML analyses, congruence between tree topologies derived from catalytic domain and CBM, was statistically evaluated using the likelihood-based test of Shimodaira & Hasewaga (SH-test)  as implemented in Tree-puzzle 5.2.
Bayesian analysis was run using MrBayes version 3.1.2 .
For the DNA datasets the best-fitting evolutionary model, as identified using Modeltest 3.5, was applied. For the protein datasets a mixed model of amino acid evolution was applied to allow model-jumping for fixed-rate models. Bayesian inference was run for 1,000,000 generations, and the first 25,000 generations were discarded as burn-in.
The authors would like to thank Prof. A. McDonald and Dr. H. Fourie (ARC-GCI, Potchefstroom, South Africa) for providing the initial Ditylenchus africanus culture and Dr. A. Elsen (K.U. Leuven, Belgium) for providing Pratylenchus coffeae. We also thank Dr. Uehara for the DNA sequence information of the Pratylenchus penetrans endoglucanases. TK is a postdoctoral researcher at Ghent University (BOF) and AH has a grant from the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen). This research was supported by a grant from Ghent University (GOA).
- Henrissat B, Bairoch A: New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochemical Journal. 1993, 293: 781-788.PubMed CentralView ArticlePubMedGoogle Scholar
- Lo N, Watanabe H, Sugimura M: Evidence for the presence of a cellulase gene in the last common ancestor of bilaterian animals. Proceedings of the Royal Society of London Series B-Biological Sciences. 2003, 270: S69-S72. 10.1098/rsbl.2003.0016.View ArticleGoogle Scholar
- Smant G, Stokkermans JPWG, Yan YT, de Boer JM, Baum TJ, Wang XH, et al: Endogenous cellulases in animals: Isolation of beta-1,4-endoglucanase genes from two species of plant-parasitic cyst nematodes. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95: 4906-4911. 10.1073/pnas.95.9.4906.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosso MN, Favery B, Piotte C, Arthaud L, de Boer JM, Hussey RS, et al: Isolation of a cDNA encoding a beta-1,4-endoglucanase in the root-knot nematode Meloidogyne incognita and expression analysis during plant parasitism. Molecular Plant-Microbe Interactions. 1999, 12: 585-591. 10.1094/MPMI.19220.127.116.115.View ArticlePubMedGoogle Scholar
- Uehara T, Kushida A, Momota Y: PCR-based cloning of two beta-1,4-endoglucanases from the root-lesion nematode Pratylenchus penetrans. Nematology. 2001, 3: 335-341. 10.1163/156854101317020259.View ArticleGoogle Scholar
- Haegeman A, Jacob J, Vanholme B, Kyndt T, Gheysen G: A family of GHF5 endo-1,4-beta-glucanases in the migratory plant-parasitic nematode Radopholus similis. Plant Pathology. 2008, 57: 581-590. 10.1111/j.1365-3059.2007.01814.x.View ArticleGoogle Scholar
- Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, et al: A molecular evolutionary framework for the phylum Nematoda. Nature. 1998, 392: 71-75. 10.1038/32160.View ArticlePubMedGoogle Scholar
- Sugimura M, Watanabe H, Lo N, Saito H: Purification, characterization, cDNA cloning and nucleotide sequencing of a cellulase from the yellow-spotted longicorn beetle, Psacothea hilaris. European Journal of Biochemistry. 2003, 270: 3455-3460. 10.1046/j.1432-1033.2003.03735.x.View ArticlePubMedGoogle Scholar
- Wei YD, Lee KS, Gui ZZ, Yoon HJ, Kim I, Zhang GZ, et al: Molecular cloning, expression, and enzymatic activity of a novel endogenous cellulase from the mulberry longicorn beetle, Apriona germari. Comparative Biochemistry and Physiology B-Biochemistry & Molecular Biology. 2006, 145: 220-229. 10.1016/j.cbpb.2006.07.007.View ArticleGoogle Scholar
- Gilkes NR, Henrissat B, Kilburn DG, Miller RC, Warren RAJ: Domains in microbial beta-1,4-glycanases – Sequence conservation, function and enzyme families. Microbiological Reviews. 1991, 55: 303-315.PubMed CentralPubMedGoogle Scholar
- Boraston AB, Bolam DN, Gilbert HJ, Davies GJ: Carbohydrate-binding modules: fine-tuning polysaccharide recognition. Biochemical Journal. 2004, 382: 769-781. 10.1042/BJ20040892.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones JT, Furlanetto C, Kikuchi T: Horizontal gene transfer from bacteria and fungi as a driving force in the evolution of plant parasitism in nematodes. Nematology. 2005, 7: 641-646. 10.1163/156854105775142919.View ArticleGoogle Scholar
- Ledger TN, Jaubert S, Bosselut N, Abad P, Rosso MN: Characterization of a new beta-1,4-endoglucanase gene from the root-knot nematode Meloidogyne incognita and evolutionary scheme for phytonematode family 5 glycosyl hydrolases. Gene. 2006, 382: 121-128. 10.1016/j.gene.2006.06.023.View ArticlePubMedGoogle Scholar
- Quillet L, Barray S, Labedan B, Petit F, Guespinmichel J: The gene encoding the beta-1,4-endoglucanase (Cela) from Myxococcus xanthus – Evidence for independent acquisition by horizontal transfer of binding and catalytic domains from Actinomycetes. Gene. 1995, 158: 23-29. 10.1016/0378-1119(95)00091-J.View ArticlePubMedGoogle Scholar
- Patthy L: Genome evolution and the evolution of exon-shuffling – a review. Gene. 1999, 238: 103-114. 10.1016/S0378-1119(99)00228-0.View ArticlePubMedGoogle Scholar
- Gilbert W, Marchionni M, Mcknight G: On the antiquity of introns. Cell. 1986, 46: 151-154. 10.1016/0092-8674(86)90730-0.View ArticlePubMedGoogle Scholar
- Doolittle WF: Genes in pieces – Were they ever together. Nature. 1978, 272:Google Scholar
- Gilbert W: The exon theory of genes. Cold Spring Harbor symposia on quantitative biology. 1987, 52: 901-905.View ArticlePubMedGoogle Scholar
- Roy SW: Recent evidence for the exon theory of genes. Genetica. 2003, 118: 251-266. 10.1023/A:1024190617462.View ArticlePubMedGoogle Scholar
- de Souza SJ, Long M, Schoenbach L, Roy SW, Gilbert W: Intron positions correlate with module boundaries in ancient proteins. Proceedings of the National Academy of Sciences of the United States of America. 1996, 93: 14632-14636. 10.1073/pnas.93.25.14632.PubMed CentralView ArticlePubMedGoogle Scholar
- Cavalier-Smith T: Selfish DNA and the origin of introns. Nature. 1985, 315: 283-284. 10.1038/315283b0.View ArticlePubMedGoogle Scholar
- Stoltzfus A, Logsdon JM, Palmer JD, Doolittle WF: Intron "sliding" and the diversity of intron positions. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94: 10739-10744. 10.1073/pnas.94.20.10739.PubMed CentralView ArticlePubMedGoogle Scholar
- Stoltzfus A: Origin of introns – Early or late. Nature. 1994, 369: 526-527. 10.1038/369526b0.View ArticlePubMedGoogle Scholar
- Dibb NJ, Newman AJ: Evidence that introns arose at proto-splice sites. Embo Journal. 1989, 8: 2015-2021.PubMed CentralPubMedGoogle Scholar
- de Souza SJ: The emergence of a synthetic theory of intron evolution. Genetica. 2003, 118: 117-121. 10.1023/A:1024193323397.View ArticlePubMedGoogle Scholar
- Rogozin IB, Wolf YI, Sorokin AV, Mirkin BG, Koonin EV: Remarkable interkingdom conservation of intron positions and massive, lineage-specific intron loss and gain in eukaryotic evolution. Current Biology. 2003, 13: 1512-1517. 10.1016/S0960-9822(03)00558-X.View ArticlePubMedGoogle Scholar
- Slamovits CH, Keeling PJ: A high density of ancient spliceosomal introns in oxymonad excavates. Bmc Evolutionary Biology. 2006, 6: 34-10.1186/1471-2148-6-34.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins L, Penny D: Complex spliceosomal organization ancestral to extant eukaryotes. Molecular Biology and Evolution. 2005, 22: 1053-1066. 10.1093/molbev/msi091.View ArticlePubMedGoogle Scholar
- Vanacova S, Yan WH, Carlton JM, Johnson PJ: Spliceosomal introns in the deep-branching eukaryote Trichomonas vaginalis. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102: 4430-4435. 10.1073/pnas.0407500102.PubMed CentralView ArticlePubMedGoogle Scholar
- Sverdlov AV, Csuros M, Rogozin IB, Koonin EV: A glimpse of a putative pre-intron phase of eukaryotic evolution. Trends in Genetics. 2007, 23: 105-108. 10.1016/j.tig.2007.01.001.View ArticlePubMedGoogle Scholar
- De Ley P, Blaxter M: Systematic Position and Phylogeny. The Biology of Nematodes. Edited by: Lee DL. 2002, London: Taylor & Francis, 1-30.View ArticleGoogle Scholar
- Subbotin SA, Sturhan D, Chizhov VN, Vovlas N, Baldwin JG: Phylogenetic analysis of Tylenchida Thorne, 1949 as inferred from D2 and D3 expansion fragments of the 28S rRNA gene sequences. Nematology. 2006, 8: 455-474. 10.1163/156854106778493420.View ArticleGoogle Scholar
- Luc M: A reappraisal of Tylenchina (Nemata). 7. The family Pratylenchidae Thorne, 1949. Revue Nématologie. 1987, 10: 203-218.Google Scholar
- Deutsch M, Long M: Intron-exon structures of eukaryotic model organisms. Nucleic Acids Research. 1999, 27: 3219-3228. 10.1093/nar/27.15.3219.PubMed CentralView ArticlePubMedGoogle Scholar
- Mitreva M, Wendl MC, Martin J, Wylie T, Yin Y, Larson A, et al: Codon usage patterns in Nematoda: analysis based on over 25 million codons in thirty-two species. Genome Biology. 2006, 7:Google Scholar
- Jacob J, Mitreva M, Vanholme B, Gheysen G: Exploring the transcriptome of the burrowing nematode Radopholus similis. Molecular Genetics and Genomics. 2008, 280: 1-17. 10.1007/s00438-008-0340-7.View ArticlePubMedGoogle Scholar
- Sverdlov AV, Rogozin IB, Babenko VN, Koonin EV: Evidence of splice signal migration from exon to intron during intron evolution. Current Biology. 2003, 13: 2170-2174. 10.1016/j.cub.2003.12.003.View ArticlePubMedGoogle Scholar
- Kaessmann H, Zollner S, Nekrutenko A, Li WH: Signatures of domain shuffling in the human genome. Genome Research. 2002, 12: 1642-1650. 10.1101/gr.520702.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan YT, Smant G, Stokkermans J, Qin L, Helder J, Baum T, et al: Genomic organization of four beta-1,4-endoglucanase genes in plant-parasitic cyst nematodes and its evolutionary implications. Gene. 1998, 220: 61-70. 10.1016/S0378-1119(98)00413-2.View ArticlePubMedGoogle Scholar
- Blumenthal T, Steward K: RNA processing and gene structure. C. elegans II. Edited by: Riddle DL, Blumenthal T, Meyer BJ, Priess JR. 1997, New York, NY: Cold Spring Harbor Laboratory PressGoogle Scholar
- Irimia M, Penny D, Roy SW: Coevolution of genomic intron number and splice sites. Trends in Genetics. 2007, 23: 321-325. 10.1016/j.tig.2007.04.001.View ArticlePubMedGoogle Scholar
- Gilbert W: Why Genes in Pieces. Nature. 1978, 271: 501-10.1038/271501a0.View ArticlePubMedGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-1155. 10.1126/science.290.5494.1151.View ArticlePubMedGoogle Scholar
- Babenko VN, Rogozin IB, Mekhedov SL, Koonin EV: Prevalence of intron gain over intron loss in the evolution of paralogous gene families. Nucleic Acids Research. 2004, 32: 3724-3733. 10.1093/nar/gkh686.PubMed CentralView ArticlePubMedGoogle Scholar
- Logsdon JM, Tyshenko MG, Dixon C, Jafari JD, Walker VK, Palmer JD: Seven newly discovered intron positions in the triose-phosphate isomerase gene – Evidence for the introns-late theory. Proceedings of the National Academy of Sciences of the United States of America. 1995, 92: 8507-8511. 10.1073/pnas.92.18.8507.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoshihama M, Nakao A, Nguyen HD, Kenmochi N: Analysis of ribosomal protein gene structures: Implications for intron evolution. Plos Genetics. 2006, 2: e25-10.1371/journal.pgen.0020025.PubMed CentralView ArticlePubMedGoogle Scholar
- Kloek AP, Sherman DR, Goldberg DE: Novel gene structure and evolutionary context of Caenorhabditis elegans globin. Gene. 1993, 129: 215-221. 10.1016/0378-1119(93)90271-4.View ArticlePubMedGoogle Scholar
- Sato Y, Niimura Y, Yura K, Go M: Module-intron correlation and intron sliding in family F/10 xylanase genes. Gene. 1999, 238: 93-101. 10.1016/S0378-1119(99)00321-2.View ArticlePubMedGoogle Scholar
- Schafer UA, Reed DW, Hunter DG, Yao K, Weninger AM, Tsang EWT, et al: An example of intron junctional sliding in the gene families encoding squalene monooxygenase homologues in Arabidopsis thaliana and Brassica napus. Plant Molecular Biology. 1999, 39: 721-728. 10.1023/A:1006172120929.View ArticlePubMedGoogle Scholar
- Rogozin IB, Lyons-Weiler J, Koonin EV: Intron sliding in conserved gene families. Trends in Genetics. 2000, 16: 430-432. 10.1016/S0168-9525(00)02096-5.View ArticlePubMedGoogle Scholar
- Charlesworth D, Liu FL, Zangh L: The evolution of the alcohol dehydrogenase gene family by loss of introns in plants of the genus Laevenworthia (Brassicaceae). Molecular Biology and Evolution. 1998, 15: 552-559.View ArticlePubMedGoogle Scholar
- Liaud MF, Brinkmann H, Cerff R: The beta-tubulin gene family of Pea – Primary structures, genomic organization and intron-dependent evolution of genes. Plant Molecular Biology. 1992, 18: 639-651. 10.1007/BF00020007.View ArticlePubMedGoogle Scholar
- Sverdlov AV, Babenko VN, Rogozin IB, Koonin EV: Preferential loss and gain of introns in 3' portions of genes suggests a reverse-transcription mechanism of intron insertion. Gene. 2004, 338: 85-91. 10.1016/j.gene.2004.05.027.View ArticlePubMedGoogle Scholar
- Boudet N, Aubourg S, Toffano-Nioche C, Kreis M, Lecharny A: Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila. Genome Research. 2001, 11: 2101-2114. 10.1101/gr.200801.PubMed CentralView ArticlePubMedGoogle Scholar
- Moody EH, Lownsbery BF, Ahmed JM: Culture of root-lesion nematode Pratylenchus vulnus on carrot disks. Journal of Nematology. 1973, 5: 225-226.PubMed CentralPubMedGoogle Scholar
- Bolla RI, Weaver C, Winter REK: Genomic differences among pathotypes of Bursaphelenchus xylophilus. Journal of Nematology. 1988, 20: 309-316.PubMed CentralPubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: The European molecular biology open software suite. Trends in Genetics. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. Journal of Molecular Biology. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.View ArticlePubMedGoogle Scholar
- Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: Identification of signaling domains. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 Research. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: A Biological Sequence Alignment Editor for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 41: 95-98.Google Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony. 2003, Sunderland, Massachusetts, USA: Sinauer AssociatesGoogle Scholar
- Carmel L, Rogozin IB, Wolf YI, Koonin EV: Evolutionarily conserved genes preferentially accumulate introns. Genome Research. 2007, 17: 1045-1050. 10.1101/gr.5978207.PubMed CentralView ArticlePubMedGoogle Scholar
- Farris JS: Methods for computing Wagner trees. Systematic Zoology. 1970, 34: 21-34.Google Scholar
- Johnson LA, Soltis DE: Assessing congruence: Empirical examples from molecular data. Molecular Systematics of Plants. Edited by: Soltis DE, Soltis P, Doyle JJ. 1998, Boston, MA: Kluwer Academic Publishers, 297-349.View ArticleGoogle Scholar
- Farris JS, Kallersjo M, Kluge AG, Bult C: Testing significance of incongruence. Cladistics-the International Journal of the Willi Hennig Society. 1994, 10: 315-319. 10.1111/j.1096-0031.1994.tb00181.x.View ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, Von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Molecular Biology and Evolution. 1999, 16: 1114-1116.View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: A sequence logo generator. Genome Research. 2004, 14: 1188-1190. 10.1101/gr.849004.PubMed CentralView ArticlePubMedGoogle Scholar
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.