- Research article
- Open Access
The promastigote surface antigen gene family of the Leishmania parasite: differential evolution by positive selection and recombination
- Alain Devault1, 2Email author and
- Anne-Laure Bañuls1
https://doi.org/10.1186/1471-2148-8-292
© Devault and Bañuls; licensee BioMed Central Ltd. 2008
- Received: 21 April 2008
- Accepted: 24 October 2008
- Published: 24 October 2008
Abstract
Background
PSA (promastigote surface antigen) is one of the major classes of membrane proteins present at the surface of the parasitic protozoan Leishmania. While it harbours leucine rich repeats, which are suggestive of its involvement in parasite-to-host physical interactions, its exact role is largely unknown. Furthermore, the extent of diversity of this gene family, both in copy number and sequence has not been established.
Results
From the newly available complete genome sequences of L. major, L. infantum and L. braziliensis, we have established the complete list of PSA genes, based on the conservation of specific domain architecture. The latter includes an array of leucine rich repeats of unique signature flanked by conserved cysteine-rich domains. All PSA genes code either for secreted or membrane-anchored surface proteins. Besides the few previously identified PSA genes, which are shown here to be part of a relatively large subclass of PSA genes located on chromosome 12, this study identifies seven other PSA subtypes. The latter, whose genes lie on chromosomes 5, 9, 21 and 31 in all three species, form single gene (two genes in one instance) subfamilies, which phylogenetically cluster as highly related orthologs. On the other hand, genes found on chromosome 12 generally show high diversification, as reflected in greater sequence divergence between species, and in an extended set of divergent paralogs. Moreover, we show that the latter genes are submitted to strong positive selection. We also provide evidence that evolution of these genes is driven by intra- and intergenic recombination, thereby modulating the number of LRRs in protein and generating chimeric genes.
Conclusion
PSA is a Leishmania family of membrane-bound or secreted proteins, whose main signature consists in a specific LRR sequence. All PSA genes found in the genomes of three sequenced Leishmania species unambiguously distribute into eight subfamilies of orthologs. Seven of these are evolving relatively slowly and could correspond to basic functions related to parasite/host interactions. On the opposite, the other PSA gene class, which include all so far experimentally studied PSA genes, could be involved in more specialised adaptative functions.
Keywords
- Leucine Rich Repeat
- Leishmania Species
- Chr31A Cluster
- Leucine Rich Repeat Domain
- Last Common Ancestor
Background
Leishmania is a parasitic protozoan of the Trypanosomatidae family and is the pathogenic agent of leishmaniases in humans [1, 2]. These diseases are endemic in 88 countries on five continents, and cover a wide range of symptoms, from asymptomatic form and benign localised cutaneous lesions to mucocutaneous lesions and visceral outcome. While localised cutaneous lesions can be self-healing, the outcome of the other forms of the disease is more disabling, and the visceral forms are generally fatal if untreated by chemotherapy. The nature of these symptoms appears statistically associated with the parasite species. L. major and L. tropica mainly produce benign cutaneous forms, L. braziliensis, mucocutaneous forms and L. infantum and L.donovani generally cause visceral disorders. Leishmania undergoes two major developmental stages in its life cycle: one stage in an invertebrate host, the phlebotomine sand fly and one stage in a vertebrate host (wide range of hosts such as humans, canids, rodents, etc.). In the former, it proliferates in the sandfly gut as a flagellated form, promastigote, and later differentiates into the infectious metacyclic promastigote form while in vertebrate hosts it proliferates intracellularly as an a-flagellated form, amastigotes, after infection of macrophages. The parasite transits from one host to the other when the sandfly feeds by sucking the blood of the vertebrate host.
To understand host-parasite interactions at the molecular level, many studies have focused on the potential role of macromolecules present at the surface of the parasite. These include glycoinositolphospholipids (GIPL), lipophosphoglycan (LPG), as well as the membrane proteins proteophosphoglycan (PPG), MSP/GP63 endopeptidases, PSA (PSA-2/GP46) and amastins. Most of these membrane proteins are glycosylated and anchored via a covalent glycosylphosphatidylinositol (GPI) moiety [3]. Some isoforms of PPG and PSA are also exported outside the parasite as soluble proteins [4, 5]. LPG plays an important role in binding to the phlebotomine midgut. Furthermore, in conjunction with the MSP proteins, it is involved in the first steps of the vertebrate host infection. They promote resistance to complement, binding to and internalisation in macrophages through the CR3 ligand and fibronectin receptors, and finally inhibition of the oxidative burst [6, 7]. GIPLs and the cysteine proteases CPB are involved in later stages of infection like protection against the nitric oxide production and modulation of the immune response [8].
PSA proteins have been detected in both promastigote and amastigote stages [9, 10]. However, the levels drastically increase from exponentially growing to stationary phase parasites, which suggests that PSA are strongly overexpressed in metacyclic promastigotes. This overexpression seems linked with the virulence status of the parasites since serial passaged in vitro cultured promastigotes with low virulence do not show this burst of PSA expression in stationary phase [10]. This expression profile is consistent with the only known role for PSA: resistance to complement lysis [11]. But other roles are possible, either in the late parasite-phlebotomine or metacyclic/amastigote-macrophage interactions. The most suggestive functional determinant in the PSA2 primary structure is the presence of leucine rich repeats (LRR). LRRs are primarily known to be involved in protein-protein and protein-glycolipids interactions. Interestingly, many known LRR containing proteins are involved in host-pathogen interactions. Yersinia invasin, Listeria internalins and Streptococci LrrG protein are surface proteins that play key roles in binding and internalisation of these bacteria in their mammalian hosts [12, 13]. Likewise, the mammalian Toll-like receptor and Nod families [14], and the plant NBS-LRR resistance proteins [15] are primary actors of the innate immunological system that recognize and bind proteins, glycolipids and nucleic acids of pathogens.
As for the other membrane proteins above, PSAs are encoded by multicopy genes. Up to now, only one or two PSA genes have been described in L. amazonensis [16], L. major [9] and L. infantum [17], and before the complete sequencing of L. major genome, the number of copies in a single species was largely unknown. Here we have taken advantage of the newly available complete genome sequences of the three species L. major, L. infantum and L. braziliensis [18, 19] to draw the precise structure of the PSA gene family. After analysis of the gene phylogeny, which might serve in the identification of distinct PSA functions, we have investigated for the presence of recombination and positive selection as driving forces acting on the evolution of this multigene family.
Results
The PSAgene family
Analysis of the primary structure of PSA proteins. A) Consensus sequence of the PSA LRR repeats. Numbers indicate the percentage of occurrence in the LRR repeats of the corresponding residues above. Residues in position #3, 7, 10, 13, 16, 18, 23(or 24) are always hydrophobic. B) Schematic representation of the PSA domain architecture. S.P: signal peptide; T/S: threonine/serine TM: transmembrane domain.
Phylogenetic analysis of the PSA gene family. A phylogenetic tree was inferred by PhyML for a DNA alignment of the PSA ORFs of L. braziliensis, L. major and L. infantum, excluding the Thr/Ser and Cys-rich/transmembrane C-terminal domains as well as all but the C-ter most LRR repeat (see Additional file 1). Only bootstrap values (from 500 replicates) greater than 74% are indicated. Very similar trees (with same significant clustering of genes) were obtained with DNA or protein sequence alignments, using both maximum likelihood or distance methods (not shown). SP: signal peptide; #LRR: number of LRR repeats (number in brackets refers to additional LRRs that diverge from the consensus LRR sequence). T/S: number of Thr and Ser residues in the Thr-rich domain (or Ser-rich, indicated by (S)). CRTM: cysteine-rich and C-ter transmembrane domains. GPI: GPI anchor. Filled and open squares indicate the presence or absence of the above features, respectively.
The main signature of PSA proteins is their specific 24–25 amino-acids LRR motif, of consensus sequence G(T/S)LPxxWxx (M/I)xxLxxLxLxxxx(x)(V/L/I)(S/T) (Fig. 1). The most striking feature of this motif relative to other LRRs is the presence of quasi invariant glycine, proline and tryptophane residues at positions 1, 4 and 7, respectively. This motif is not present in the predicted proteins from the fully sequenced genomes of two other trypanosomatidae, T. brucei and T. cruzi, and no clear homologs of PSA could be identified in these organisms. The only PSA domain for which homology was found in the other trypanosomatidae is the Thr/Ser-rich domain. The T(n)KP2 and T(n)EAPT repeats found in most of these domains are present in T. cruzi mucins. However, while it makes up the bulk of the mucin protein architecture, this domain is not always present in PSA proteins (see Fig. 2), and mucins do not contain any LRR repeat. Therefore, PSA proteins and mucins are likely to play different roles. Actually, the LRRs mostly similar to those of Leishmania PSA proteins belong to the plant NBS-LRR resistance proteins.
PSA genes are dispersed on chromosomes (Chr) 5, 9, 12, 21 and 31 for all three Leishmania species. Chromosomes 9 and 21 contain a single gene, while Chr 5 contains 2 genes, and Chr 31, 3 genes (4 in L. major). PSA genes localised on different chromosomes (and also those of the different Chr31 loci) are sufficiently divergent to define PSA gene classes (see below). Inside each chromosome, all PSA genes are clustered in tandem, except those of Chr 5 and gene LbrM31_V2.3700, which are well separated. The PSA subfamily present on Chr 12 is the most complex (except in L. braziliensis where there is only one gene). In L. major, the 24 Chr 12 PSA genes belong to a single cluster of tandem repeats, with many of these genes separated by a single non-related intervening gene. This pattern is also observed for the 7 Chr12 PSA genes of L. infantum. All PSA genes arranged in clusters have the same orientation relative to transcription. The only PSA genes studied so far in vivo belong to the Chr12 array. Although the PSA genes we have identified on the other chromosomes show significant sequence divergence with the Chr12 genes (see below), the conservation of the precise domain organisation between all these genes justifies their being included in the same unique PSA family.
Phylogeny of the PSAgene family
An alignment of all PSA protein sequences from all three Leishmania species was generated. For phylogenic analysis, only the domains present in all PSA proteins were included in the alignment. The Thr/Ser and Cys-rich/transmembrane C-terminal domains were thus discarded, as well as the central LRR domain (except for the last LRR motif of each sequence, which was retained). LRR units probably evolve not only by nucleotide substitution but also by birth and loss of repeats following recombination events. Therefore, alignments of rather distant LRR domains (with low similarity between repeats and varying number of repeats) may not reflect properly their evolution. This truncated alignment was transformed into its DNA coding version (Additional file 1), which was used to generate a phylogenetic tree by PhyML (Fig. 2). This tree shows the existence of eight PSA subfamilies which can be defined by their chromosome localisation: Chr5 (A and B), 9, 12, 21 and three subfamilies on Chr31 (A, B, C). Subfamilies Chr5A and B are the most distantly related and pair apart from all other clusters. Among the latter, only the Chr21/Chr31A subfamilies can be phylogenetically linked with high confidence bootstrap values. Each of the three Leishmania species has at least one gene in each of these eight subfamilies. The L. infantum gene orthologs always group with those of L. major, which is expected from the known species tree. The relatively short branch length separating orthologs in comparison to those separating the inparalogs (ancestors of the seven gene subfamilies) might suggest that speciation of the three Leishmania is a late event in the evolution history of the PSA genes. In keeping with this interpretation, a tree with similar relative branch lengths was inferred when using the rate of synonymous substitutions (dS), a more neutral distance metric (not shown).
Chr5A, B, 9, 21, 31A, B, C subfamilies of PSAgenes
While the Chr12 subfamily is most diversified, having very different numbers (1 to 24) of rather distantly related paralogs in the three species, the PSA genes lying on the other chromosomes form very simple clusters. For each of the latter subfamilies, each species contains a single gene (except L. major, which has two highly similar paralogs for the Chr31A cluster) and the orthologs or pseudo-orthologs are quite similar (median = 75% amino acid identity for all three species and 93% for L. infantum/L. major orthologs). This degree of homology between orthologs, which was calculated on the same domains used above for phylogenetic tree determination, was also observed over the entire LRR domain. Indeed, confident alignments covering this domain can be obtained inside each of the Chr5A and B, 9, 21, 31A, B, and C clusters since the number of LRRs differ at most by one inside each subfamily and the similarity between orthologs is (median = 80% amino acid identity for all three species and 96% for L. infantum/L. major orthologs) (see Additional files 2, 3, 4, 5, 6, 7, 8 for sequence alignments). On the other hand, similarity between genes of different subfamilies is rather low over this domain (median = 32% amino acid identity, and much difference in number of repeats). This may suggest that the LRR domains characteristic of each cluster represent specific binding properties, which could confer specialized roles to their respective PSA proteins. Another link between these subfamilies and functional specialisation is given by the nature of the C-terminal domain. Indeed, all members of the same subfamily share the same membrane anchorage/secretion determinant. In particular, all genes of the phylogenetically linked Chr21 and 31A subgroups are predicted to code for secreted proteins (Fig. 2).
Chr12 subfamily of PSAgenes
Phylogenetic tree of the LRR domains of PSA genes located on chromosome 12. A phylogenetic tree was inferred by PhyML for a DNA alignment of the LRR domain (see Additional files 9 and 13 for protein and DNA alignments, respectively) of all Chr 12 PSA genes of L. braziliensis, L.major and L. infantum. Bootstrap values (from 500 replicates) greater than 74% are indicated. Very similar trees (with same significant clustering of genes) were obtained with DNA or protein sequence alignments and by maximum likelihood or distance methods (not shown). SP: signal peptide; #LRR: number of LRR repeats. T/S: number of Thr and Ser residues in the Thr-rich domain (or Ser-rich, indicated by (S)). CRTM: cysteine-rich and C-ter transmembrane domains. GPI: GPI anchor. Filled and open squares indicate the presence or absence of the above features, respectively.
Recombination in PSAgenes
Phylogenetic clustering of individual PSA domains. Trees were inferred by Neighbor-Joining (NJ) from nucleotide distance matrix (TN93 model) calculated individually for all the indicated domains of the L. major Chr12 PSA genes. Individual LRR domains that were present in only a few genes were not analysed. For each domain, those genes that cluster together in the tree (with bootstrap values greater than 74%) are boxed with the same colour. Grey boxes indicate genes that do not cluster with any other genes. Absence of a box indicates that the corresponding domain is absent in the gene. N-ter: PSA N-terminal domain up to the first LRR repeat. LRRs are numbered from N-ter to C-ter according to the LRR alignment (see Additional file 9). ω : refers to the last (C-ter most) LRR of each gene; therefore, LRRω + (ω-1) corresponds to the sequence comprising the last and penultimate LRRs in each gene. (part.): partial; for example, 4 (part.) means portion of LRR 4. CGC-C2: the CGCx(23–26)CxxxxxC domain at the C-ter end of the LRR domain.
Analysis of chimeric PSA genes by Simplot. A) LmjF12.0990. A DNA alignment of the indicated genes covering LRR #2, 9 and the last 4 C-ter most LRRs, together with the N-ter and CGC-C2 flanking domains (Additional file 14) was analysed with Simplot. Single or multiple genes from each phylogenetic PSA gene cluster of L. major Chr 12 (see Fig. 3) were first used for similarity analysis (Simplot) against the LmjF12.0990 putative chimeric gene. The six genes of cluster 12C were reduced to a consensus sequence. The best pair of candidate parental genes (LmjF12.0760 and consensus sequence for cluster 12C) and a control gene were then used for bootscanning and informative sites analysis (MaxChi2). B) LmjF12.0960. The same analysis was performed as in A) using a DNA alignment covering LRR #2, 9 and the last 3 C-ter most LRRs (see Additional file 15).
More profound rearrangements can be created by intergenic recombinations (both at coding and non-coding sites), which lead to whole gene duplication and deletions. One such rearrangement can be traced in the L. major Chr12 PSA gene cluster, where PSA genes are organised in tandem, tough often separated by one of two unrelated genes of type A or B. The distribution of these intervening genes together with the phylogenic signature of each PSA genes suggest that the two pairs of very similar paralogs LmjF12.0810/LmjF12.0940 and LmjF12.0830/LmjF12.0960 were created by a same segmental duplication involving a bloc of five genes (see Additional file 10).
Positive selection of PSAgenes
d N / d S vs d S plots for PSA genes pairwise comparisons. DNA alignments covering the LRR domain were generated for each of the Chr5A, B, 9, 12, 21, 31A, B, C PSA gene clusters (including genes of all three Leishmania species; see Additional file 13 for alignment of Chr12 genes). After reducing very homologous genes (with dS<0.02) to single entities, dN/dS and dS values were calculated for all pairwise comparions inside a given cluster. Orange triangles: Chr5A, B, 9, 21, 31A, B, C clusters; blue squares: Chr12 cluster. Arrows point to values for L. major-L. infantum ortholog comparisons.
Distribution of the predicted positively selected sites along the LRR amino acid consensus sequence. A) A DNA alignment covering the LRR domain plus the flanking N-ter and CGC-C2 domains of all Chr12 PSA genes (Additional file 16) was analysed by codeml. A) The predicted positively selected sites were plotted against their positioning in the LRR domain. Green bars represent the number of predicted positive sites with probabilities greater than 95% and blue bars, eventual increments when considering sites with probabilities greater than 90%. o, i and # refer to positions predicted to reside outside, inside or "in between" in the PSA LRR 3D structure. B) Means of dN/dS (ω) at each residue position in the LRR were calculated from the posterior dN/dS values of all sites sharing the same position in the LRR. Filled bars: variable residues; open bars: conserved residues.
3' Intergenic region
Expression of some PSA genes of Chr12 has been shown to be stage specific. Strong increases in the levels of both mRNA and protein were observed during metacyclogenesis in vitro [10]. This regulation was shown to depend on the 3' intergenic region (IR) of the PSA gene [17]. The question arises, then, as to whether the different classes of PSA might be regulated the same way. We therefore wished to compare the phylogeny of the 3' (IR) to that obtained above for the coding sequence. Around 2 kb of the 3' IR of PSA genes from all three Leishmania species (with few exceptions for L. infantum and L. braziliensis due to gaps in their sequences) were analysed by local BLAST and aligned by Dialign2. We found that the 3' IRs are conserved between all orthologs (and pseudo-orthologs for he Chr12 clade), but show no significant homology between genes of the 8 different PSA clades (05, 09, 12, 21, 31A, B, C). While the similarities found for pairwise comparisons involving L. braziliensis sequences are variable and relatively low, those of the other comparisons are high. The L. infantum/L. major orthologs of clades 5, 9, 21, 31A, B, C are equally well conserved (ranging from 87 to 92% identity; see Additional file 11 for sequence alignments), as are the set of pseudo-orthologs and paralogs of clade 12 (median = 95% identity; see Additional file 12 for sequence alignment). Yet, the alignment of the 3' IR of these Chr12 PSA genes is quite mosaic, showing large blocks of deletions/insertions. Phylogenetic analysis of an alignment retaining 1.2 kb of domains conserved in at least 70% of the genes revealed that Chr12 genes cluster according to the species (not shown). This is reminiscent of the extremities of the coding regions (Fig. 2). Functional studies will be necessary to determine if this high degree of conservation is the result of gene conversion or rather, of purifying selection driven by imperatives of gene expression. In any case, if indeed the 3' IR of the different PSA clades play a role in the control of gene expression, either the pattern of expression is variable between the clades or the cues governing this control are subtler than the primary sequence per se.
Discussion and Conclusion
Thanks to the whole genome sequencing of three Leishmania species, we could draw the complete picture of the PSA gene family. Up to now, only a few genes belonging to the Chr12 subfamily had been characterised. Also, L. braziliensis was thought not to contain any gene for that subfamily. We have proposed a nomenclature for all PSA genes (see Fig. 2 and 3). We have kept the PSA letters, even tough it is now known that these genes are expressed both in promastigote and amastigote stages and code not only for membrane proteins at the surface of the parasite but also for excreted soluble proteins. Follow the chromosome number and other symbols that refer to the phylogenetic PSA gene clusters.
PSAs are basically LRR repeats containing proteins, which are either membrane-bound at the surface of Leishmania or soluble as excreted proteins. These features suggest that this protein family could be involved in host-parasite interactions. The diversity found in PSAs, which concerns mainly the LRR domain, would be consistent with this role. This diversity is not accounted only by the sequence divergence of the LRR repeats but also by the number of these repeats present in each protein, which varies from 2 to 25 and could influence both the shape and the space occupied by PSAs, and thereby their binding properties. Two other factors of diversity in PSA reside outside the LRR domain. One concerns the Thr/Ser-rich domain, whose length varies considerably and can be totally absent in some PSAs. Although there is some indication that PSAs could be glycosylated [24, 25], both the sites of glycosylation and the nature of the amino acid-glycan linkage are unknown. Most of these Thr/Ser-rich domains contain a signature present in T. cruzi mucins which is known to be the site of O-glycosylation. Finally, membrane attachment determinants (which result mostly in anchoring through GPI) define two destinations for PSA proteins: excretion or attachment to the parasite's surface. Although the majority are predicted to be membrane bound, PSAs predicted to be secreted are found both as strongly conserved proteins between species (PSA21 in all three species; LiPSA12A and LmPSA12A) and as species specific proteins (LiPSA12I and most LmPSA12C). This suggests that the presence of these C-terminal truncated forms of PSAs is neither fortuitous nor the products of pseudogenes, but that secretion plays a fundamental role in Leishmania.
PSA genes distribute into eight subfamilies. In view of the found orthologies, the most parsimonious scenario would be that the last common ancestor (LCA) of L. major, L. infantum and L. braziliensis had nine genes: one in each of chromosomes 9, 12 and 21, two on chromosome 5 and three on chromosome 31. Among the PSA subfamilies, the Chr12 subfamily is the most diversified one, at least in L. major and L. infantum. This is brought about by nucleotide substitutions as well as unequal gene recombination, which modulates both the number of LRRs in PSA proteins and the number of genes in the family, and generates chimeric genes. The evolution of the Chr12 subfamily of PSA genes is probably best explained by a birth-and-loss process, by which the gene family is either expanded or contracted, while PSA genes are individually submitted to diversifying selection. We have detected clear evidence of positive selection for the Chr12 genes in both L. major and L. infantum, while the existence of orthologs, divergent paralogs, as well as probable pseudogenes (LmjF12.1005, LinJ12_v4.0662) all point to a birth-and-loss model. We note, however, that the extremities of the coding regions and the 3' IR of the Chr12 PSA genes are all species specific (no orthologies can be inferred) and remarkably similar (some stretches are even devoid of synonymous substitutions), suggesting that some gene homogenisation might have taken place by gene conversion. In contrast to the Chr12 subfamily, the seven other subfamilies seem to evolve in a more conservative way. Not only all genes of Chr5, 9, 21 and 31 from the three Leishmania species cluster as more similar orthologs but evidence of recombination is scarce. One finds only one pair of paralogs (the L. major Chr31A subfamily), and the number of LRR repeats is the same for each group of orthologs (except one difference in Chr21 subfamily). Furthermore, diversifying selection apparently does not act on these genes, as evidenced by dN/dS values well under those obtained for the Chr12 genes. We would therefore postulate the latter subfamilies to be involved in basic slowly evolving host-parasite interactions, whereas the Chr12 subfamily could be linked to more dynamic adaptative responses like optimisation of cellular infection (internalisation and maintenance in the macrophage) and escape from the immune system. The latter could also be related to the Leishmania species identities like the nature of the sand fly species infected and the clinical symptoms. Still, diversification of Chr12 genes apparently does not offer any advantage for L. braziliensis. At this stage, it is hard to find out what kind of adaptative role these genes could play in L. major and L. infantum, which would be dispensable for L. braziliensis.
Diversifying selection has been shown to drive evolution of genes involved in host-pathogen interactions. These include virulence genes, as found in Streptococci [26] and Listeria monocytogenes [27], as well as resistance genes, like the numerous plant NBS-LRR genes [28], the primate Trim5α [29] genes involved in restriction of retroviruses replication and many gene classes involved in the Drosophila innate immune system [30]. Interestingly, the Listeria virulence gene inlA and the NBS-LRR genes possess LRR domains, which were shown to evolve under positive selection. In the present work, we found that in the Chr12 PSA gene subfamily of two Leishmania species, LRR domains evolve under strong positive selection. Most notably, the several residues submitted to diversifying selection are predicted to reside at the protein's surface. The importance of this evolving LRR protein interface in the parasite's life cycle awaits the determination of the precise PSA biological function and the identification of eventual host PSA interactors.
In conclusion, we have identified eight PSA subfamilies. The phylogenetic distance separating them, which is paralleled by the divergence of their LRR domains, appears large enough to reflect functional specificities. The strong conservation of the orthologs of the subfamilies 5A, 5B, 9, 21, 31A, B and C suggest that they play essential roles. Generation of knock outs and gene expression analysis of all individual PSA genes will be needed to explore the importance of these genes in the different stages of the life cycle of Leishmania, and eventually in adaptation of the parasites to their specific insect or mammalian host. In comparison to the seven subfamilies above, the PSA12 subfamily is evolving much more rapidly, (at least in L. major and L. infantum) as revealed by the high number of paralogs and the presence of strong diversifying selection. Sequencing of these genes in several strains of L. major and L. infantum is expected to show a high incidence of polymorphism.
Methods
Alignment of PSA amino acid and nucleotide sequences and identification of membrane anchorage determinants
PSA protein and DNA sequences were extracted from the L. major, L. infantum and L. braziliensis database http://www.genedb.org/ after BLAST analysis. The latest versions of the genome sequences were used, i.e. 5.2 for L. major, 2.0 for L. braziliensis and 3.0 (and 4.0 for chromosome 12) for L. infantum. Protein sequences were aligned with ClustalX 2.0 [31] followed by manual correction. These alignments were transposed to the corresponding genomic DNA alignments by protal2dna http://mobyle.pasteur.fr/cgi-bin/MobylePortal/portal.py?form=tranalign. The presence of signal peptide and transmembrane domains was reported as found in annotations of the Leishmania genomes http://www.genedb.org/. The prediction of GPI-anchorage was performed by GPI-SOM http://gpi.unibe.ch/[32]. Genomic sequences flanking the PSA ORFs in 3' were aligned with Dialign 2.2.1 [33].
Phylogenetic analysis
Phylogenetic trees of Fig. 2 and 3 were obtained from nucleotide sequence alignments by maximum likelihood (PhyML v2.4.4) [34]. The selection of the nucleotide substitution model (and eventually of amino acids substitution) was determined by ModelGenerator http://distributed.cs.nuim.ie/multiphyl.php[35]. The HKY+I+G (Fig. 2) and TN93 +G+I (Fig. 3) nucleotide models were used, with the base frequency, Ts/Tv ratio, and gamma distribution parameter estimated by PhyML. These trees were compared to those obtained by PhyML from the corresponding amino acid sequences (WAG+G model of amino acids substitution) or with a distance method implemented in MEGA v3.1 [36] using both DNA (TN93) or amino acids (JTT model) alignments. Occasionally (see text in Results section), a codon model (dS, number of synonymous substitutions per synonymous sites; Nei-Gojobori method, JC model) as implemented in MEGA was used to calculate distance matrix. Trees were annotated and visualised by Treedyn (v194.3) [37].
Positive selection
One level of positive selection analysis involved calculation of dN/dS ratios between pairs of sequences (MEGA, Nei-Gojobori method, JC model). ω = dN/dS < 1 corresponds to purifying (negative) selection, ω = 1 to neutral evolution (absence of selection), and ω > 1 indicating adaptive evolution (positive selection). The maximum likelihood based method implemented in the codeml program of the PAML package (v3.15) [38] was used to detect positive selection as well as identifying positively selected sites, when the number of related sequences was sufficiently high. To detect positive selection, we performed likelihood ratio test (LRT) analysis between three pairs of models to determine whether particular models provided a significantly better fit to the data: M1a ("Nearly-Neutral" model) vs. M2a ("PositiveSelection" model), M7 ("beta" model) vs. M8 ("beta&ω" model) and M8a vs. M8. Only models M2a and M8 allow for site categories evolving under positive selection (see Yang et al. 2005 [39] for further details on the models). Positive selection is inferred when the (2Δλ) statistic is greater than critical values 1) of the Chi square distribution for a degree of freedom of 2 (M1a/M2a. and M7/M8 comparisons) or 2) of a 50:50 mixture of point mass 0 and Chi square (critical values are 2.71 at 5% and 5.41 at 1%; M8a/M8 comparison). Only predicted positively selected sites with posterior BEB (Bayes empirical Bayes) probabilities greater than 0.95 or 0.90 were retained.
Declarations
Acknowledgements
We thank F.Chevenet for advice in running the Treedyn software, A. Kajava for structure modelling of PSA LRRs and E. Desmarais for critical reading of the manuscript.
Authors’ Affiliations
References
- Schlein Y: Leishmania and Sandflies: interactions in the life cycle and transmission. Parasitol Today. 1993, 9 (7): 255-258.View ArticlePubMedGoogle Scholar
- Banuls AL, Hide M, Prugnolle F: Leishmania and the leishmaniases: a parasite genetic update and advances in taxonomy, epidemiology and pathogenicity in humans. Adv Parasitol. 2007, 64: 1-109.View ArticlePubMedGoogle Scholar
- Frommel TO, Button LL, Fujikura Y, McMaster WR: The major surface glycoprotein (GP63) is present in both life stages of Leishmania. Mol Biochem Parasitol. 1990, 38 (1): 25-32.View ArticlePubMedGoogle Scholar
- Jimenez-Ruiz A, Boceta C, Bonay P, Requena JM, Alonso C: Cloning, sequencing, and expression of the PSA genes from Leishmania infantum. Eur J Biochem. 1998, 251 (1–2): 389-397.View ArticlePubMedGoogle Scholar
- Ilg T: Proteophosphoglycans of Leishmania. Parasitol Today. 2000, 16 (11): 489-497.View ArticlePubMedGoogle Scholar
- Rittig MG, Bogdan C: Leishmania-host-cell interaction: complexities and alternative views. Parasitol Today. 2000, 16 (7): 292-297.View ArticlePubMedGoogle Scholar
- Santarem N, Silvestre R, Tavares J, Silva M, Cabral S, Maciel J, Cordeiro-da-Silva A: Immune response regulation by leishmania secreted and nonsecreted antigens. J Biomed Biotechnol. 2007, 2007 (6): 85154-PubMed CentralPubMedGoogle Scholar
- Olivier M, Gregory DJ, Forget G: Subversion mechanisms by which Leishmania parasites can escape the host immune response: a signaling point of view. Clin Microbiol Rev. 2005, 18 (2): 293-305.PubMed CentralView ArticlePubMedGoogle Scholar
- Handman E, Osborn AH, Symons F, van Driel R, Cappai R: The Leishmania promastigote surface antigen 2 complex is differentially expressed during the parasite life cycle. Mol Biochem Parasitol. 1995, 74 (2): 189-200.View ArticlePubMedGoogle Scholar
- Beetham JK, Donelson JE, Dahlin RR: Surface glycoprotein PSA (GP46) expression during short- and long-term culture of Leishmania chagasi. Mol Biochem Parasitol. 2003, 131 (2): 109-117.View ArticlePubMedGoogle Scholar
- Lincoln LM, Ozaki M, Donelson JE, Beetham JK: Genetic complementation of Leishmania deficient in PSA (GP46) restores their resistance to lysis by complement. Mol Biochem Parasitol. 2004, 137 (1): 185-189.View ArticlePubMedGoogle Scholar
- Pizarro-Cerda J, Cossart P: Bacterial adhesion and entry into host cells. Cell. 2006, 124 (4): 715-727.View ArticlePubMedGoogle Scholar
- Seepersaud R, Hanniffy SB, Mayne P, Sizer P, Le Page R, Wells JM: Characterization of a novel leucine-rich repeat protein antigen from group B streptococci that elicits protective immunity. Infect Immun. 2005, 73 (3): 1671-1683.PubMed CentralView ArticlePubMedGoogle Scholar
- Athman R, Philpott D: Innate immunity via Toll-like receptors and Nod proteins. Curr Opin Microbiol. 2004, 7 (1): 25-32.View ArticlePubMedGoogle Scholar
- McHale L, Tan X, Koehl P, Michelmore RW: Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006, 7 (4): 212-PubMed CentralView ArticlePubMedGoogle Scholar
- Lohman KL, Langer PJ, McMahon-Pratt D: Molecular cloning and characterization of the immunologically protective surface glycoprotein GP46/M-2 of Leishmania amazonensis. Proc Natl Acad Sci USA. 1990, 87 (21): 8393-8397.PubMed CentralView ArticlePubMedGoogle Scholar
- Myung KS, Beetham JK, Wilson ME, Donelson JE: Comparison of the post-transcriptional regulation of the mRNAs for the surface proteins PSA (GP46) and MSP (GP63) of Leishmania chagasi. J Biol Chem. 2002, 277 (19): 16489-16497.View ArticlePubMedGoogle Scholar
- Peacock CS, Seeger K, Harris D, Murphy L, Ruiz JC, Quail MA, Peters N, Adlem E, Tivey A, Aslett M, et al: Comparative genomic analysis of three Leishmania species that cause diverse human disease. Nat Genet. 2007, 39 (7): 839-847.PubMed CentralView ArticlePubMedGoogle Scholar
- Ivens AC, Peacock CS, Worthey EA, Murphy L, Aggarwal G, Berriman M, Sisk E, Rajandream MA, Adlem E, Aert R, et al: The genome of the kinetoplastid parasite, Leishmania major. Science. 2005, 309 (5733): 436-442.PubMed CentralView ArticlePubMedGoogle Scholar
- Kajava AV: Structural diversity of leucine-rich repeat proteins. J Mol Biol. 1998, 277 (3): 519-527.View ArticlePubMedGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73 (1): 152-160.PubMed CentralPubMedGoogle Scholar
- Anisimova M, Nielsen R, Yang Z: Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics. 2003, 164 (3): 1229-1236.PubMed CentralPubMedGoogle Scholar
- Kajava AV, Kobe B: Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information. Protein Sci. 2002, 11 (5): 1082-1090.PubMed CentralView ArticlePubMedGoogle Scholar
- Kahl LP, McMahon-Pratt D: Structural and antigenic characterization of a species- and promastigote-specific Leishmania mexicana amazonensis membrane protein. J Immunol. 1987, 138 (5): 1587-1595.PubMedGoogle Scholar
- Murray PJ, Spithill TW, Handman E: The PSA-2 glycoprotein complex of Leishmania major is a glycosylphosphatidylinositol-linked promastigote surface antigen. J Immunol. 1989, 143 (12): 4221-4226.PubMedGoogle Scholar
- Anisimova M, Bielawski J, Dunn K, Yang Z: Phylogenomic analysis of natural selection pressure in Streptococcus genomes. BMC Evol Biol. 2007, 7: 154-PubMed CentralView ArticlePubMedGoogle Scholar
- Orsi RH, Ripoll DR, Yeung M, Nightingale KK, Wiedmann M: Recombination and positive selection contribute to evolution of Listeria monocytogenes inlA. Microbiology. 2007, 153 (Pt 8): 2666-2678.View ArticlePubMedGoogle Scholar
- Lehmann P: Structure and evolution of plant disease resistance genes. J Appl Genet. 2002, 43 (4): 403-414.PubMedGoogle Scholar
- Sawyer SL, Wu LI, Emerman M, Malik HS: Positive selection of primate TRIM5alpha identifies a critical species-specific retroviral restriction domain. Proc Natl Acad Sci USA. 2005, 102 (8): 2832-2837.PubMed CentralView ArticlePubMedGoogle Scholar
- Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39 (12): 1461-1468.View ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882.PubMed CentralView ArticlePubMedGoogle Scholar
- Fankhauser N, Maser P: Identification of GPI anchor attachment signals by a Kohonen self-organizing map. Bioinformatics. 2005, 21 (9): 1846-1852.View ArticlePubMedGoogle Scholar
- Morgenstern B: DIALIGN 2: improvement of the segment-to-segment approach to multiple sequence alignment. Bioinformatics. 1999, 15 (3): 211-218.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704.View ArticlePubMedGoogle Scholar
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McLnerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163.View ArticlePubMedGoogle Scholar
- Chevenet F, Brun C, Banuls AL, Jacq B, Christen R: TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics. 2006, 7: 439-PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
- Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22 (4): 1107-1118.View ArticlePubMedGoogle Scholar
Copyright
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.