Microsatellite and single-nucleotide polymorphisms indicate recurrent transitions to asexuality in a microsporidian parasite.

Assessing the mode of reproduction of microparasites remains a difficult task because direct evidence for sexual processes is often absent and the biological covariates of sex and asex are poorly known. Species with geographically divergent modes of reproduction offer the possibility to explore some of these covariates, for example, the influence of life-history traits, mode of transmission and life-cycle complexity. Here, we present a phylogeographical study of a microsporidian parasite, which allows us to relate population genetic structure and mode of reproduction to its geographically diverged life histories. We show that in microsporidians from the genus Hamiltosporidium, that use the cladoceran Daphnia as host, an epidemic population structure has evolved, most probably since the last Ice Age. We partially sequenced three housekeeping genes (alpha tubulin, beta tubulin and hsp70) and genotyped seven microsatellite loci in 51 Hamiltosporidium isolates sampled within Europe and the Middle East. We found two phylogenetically related asexual parasite lines, one each from Fennoscandia and Israel, which share the unique ability of being transmitted both vertically and horizontally from Daphnia to Daphnia. The sexual forms cannot transmit horizontally among Daphnia, but presumably have a complex life cycle with a second host species. In spite of the similarities between the two asexual lineages, a clustering analysis based on microsatellite polymorphisms shows that asexual Fennoscandian parasites do not share ancestry with any other Hamiltosporidium that we have sampled. Moreover, allele sequence divergence at the hsp70 locus is twice as large in Fennoscandian than in Israeli parasites. Our results indicate that asexual reproduction evolved twice independently, first in Fennoscandian and more recently in the Israeli parasites. We conclude that the independent origin of asexuality in these two populations is associated with the altered parasite mode of transmission and the underlying dynamics of host populations.

Reappraising the theme of breeding systems in Echinococcus: is outcrossing a rare phenomenon?

Selfing has been considered the most common mode of reproduction in Echinococcus flatworms. However, population genetic studies on the asexual larval stage involving nuclear co-dominant markers have not always revealed significant heterozygote deficiencies--the expected outcome of a regularly and highly inbred population. In this study, we analysed the genetic structure of Echinococcus granulosus sensu lato populations from Southern Brazil during their adult (sexual) stage using 1 mitochondrial and 1 nuclear marker (cox 1 and mdh, respectively). We show that parasite genetic differentiation is largest among definitive hosts (domestic dogs) from different farms, suggesting that transmission is mostly maintained within a farm. Moreover, we show that heterozygote deficiencies are not significant, and we suggest that outbreeding is the most common mode of reproduction of the parasite in that region.

The EmsB tandemly repeated multilocus microsatellite: a new tool to investigate genetic diversity of Echinococcus granulosus sensu lato.

Cystic echinococcosis (CE) is a widespread and severe zoonotic disease caused by infection with the larval stage of the eucestode Echinococcus granulosus sensu lato. The polymorphism exhibited by nuclear and mitochondrial markers conventionally used for the genotyping of different parasite species and strains does not reach the level necessary for the identification of genetic variants linked to restricted geographical areas. EmsB is a tandemly repeated multilocus microsatellite that proved its usefulness for the study of genetic polymorphisms within the species E. multilocularis, the causative agent of alveolar echinococcosis. In the present study, EmsB was used to characterize E. granulosus sensu lato samples collected from different host species (sheep, cattle, dromedaries, dogs, and human patients) originating from six different countries (Algeria, Mauritania, Romania, Serbia, Brazil, and the People's Republic of China). The conventional mitochondrial cox1 and nad1 markers identified genotypes G1, G3, G5, G6, and G7, which are clustered into three groups corresponding to the species E. granulosus sensu stricto, E. ortleppi, and E. canadensis. With the same samples, EmsB provided a higher degree of genetic discrimination and identified variations that correlated with the relatively small-scale geographic origins of the samples. In addition, one of the Brazilian single hydatid cysts presented a hybrid genotypic profile that suggested genetic exchanges between E. granulosus sensu stricto and E. ortleppi. In summary, the EmsB microsatellite exhibits an interesting potential for the elaboration of a detailed map of the distribution of genetic variants and therefore for the determination and tracking of the source of CE.

Using mitochondrial and nuclear markers to evaluate the degree of genetic cohesion among Echinococcus populations.

Based on the distinctiveness of their mitochondrial haplotypes and other biological features, several recent publications have proposed that some Echinococcus granulosus strains should be regarded as separate species. However, the genetic cohesion of these species has not been extensively evaluated using nuclear markers. We assess the degree of polymorphism of the partial mitochondrial cox1 (366bp), the nuclear mdh (214bp) and EgAgB4 (281-283bp) genes of E. granulosus sensu lato isolates collected from areas where different strains occur sympatrically. Five distinct mitochondrial haplotypes were determined by direct sequencing (G1, G2, G5, G6 and G7). The mdh genotypes were first screened by SSCP: three alleles were identified (Md1-Md3), which were further confirmed by nucleotide sequencing. For EgAgB4, which was analysed by direct sequencing the PCR products, two groups of sequences were found: EgAgB4-1 and EgAgB4-2. No haplotype-specific mdh or EgAgB4 sequences occur. Nevertheless, alleles Md1 and Md2 and type 1 sequences of EgAgB4 showed a higher frequency within the group of haplotypes G1-G2, while allele Md3 and EgAgB4-2 are most frequent in the G5-G7 cluster. By AMOVA it is shown that 79% of the total genetic variability is found among haplotype groups. These findings are compatible with two not mutually exclusive evolutionary hypotheses: (a) that haplotypes share an ancestral polymorphism, or (b) that the reproductive isolation between parasites with distinct haplotypes is not complete, leading to gene introgression. The biologic and epidemiologic consequences of our findings are discussed.

Expression and diversity of Echinococcus multilocularis AgB genes in secondarily infected mice: evaluating the influence of T-cell immune selection on antigenic variation.

The T-cell-mediated immune response exhibits a crucial function in the control of the intrahepatic proliferation of Echinococcus multilocularis larvae in mice and humans, both being natural intermediate hosts of the parasite. Antigen B (AgB), a metabolized Echinococcus spp. lipoprotein, contributes to the modulation of the T-cell immune response, and distinct sites of the corresponding AgB1, AgB3 and AgB4 genes were shown to be under positive selection pressure. Since several AgB gene variants are present in a single Echinococcus metacestode, we used secondary E. multilocularis infections in BALB/c and in athymic nude mice (devoid of T-cell responses) to analyze the effect of the cellular immune response on the expression and diversity of EmAgB1-EmAgB4 genes. We demonstrated hereby that EmAgB transcripts were less abundant in nude mice during the early phase of infection (at one month post-infection), and that EmAgB2 is simultaneously down-regulated when compared to the other three genes. A negative relationship exists between the level of transcription and diversity of EmAgB genes. Moreover, no excess of non-synonymous substitutions was found among the distinct EmAgB alleles from a single host. Together, these results pointed to the effect of purifying selection, which seemed to eliminate the detrimental AgB variants generated during the development of the metacestode within the peritoneal cavity of its intermediate host.

Searching for antigen B genes and their adaptive sites in distinct strains and species of the helminth Echinococcus.

Twenty-seven PCR-derived antigen B (AgB) nucleotide sequences from four Echinococcus species (Echinococcus granulosus, Echinococcus multilocularis, Echinococcus oligarthrus and Echinococcus vogeli) were aligned with 78 already published sequences, to generate a maximum likelihood phylogeny of the AgB multigene family. The phylogenetic analysis confirms that the family is constituted by four groups of genes present in each one of the four species (AgB1, AgB2, AgB3 and AgB4), and suggests that it originated by ancient duplication events preceding speciation within the genus. AgB5 sequences, which had been formerly suggested to correspond to a putatively new AgB subunit, cluster with AgB3. Likelihood tests suggest that AgB gene evolution may have been driven by heterogeneous selection pressures acting on particular AgB1, AgB3 and AgB4 codons. No selection is detected in AgB2. We discuss implications of our findings in terms of AgB biology and its use as a diagnostic tool.

Contingent, non-neutral evolution in a multicellular parasite: natural selection and gene conversion in the Echinococcus granulosus antigen B gene family.

Recent studies have demonstrated that the Echinococcus granulosus antigen B (AgB) interferes with the intermediate hosts' immune response and is encoded by a multigene family. The number of members within the family is still uncertain, but there are several evidences of a large genetic variability. The E. granulosus AgB genomic sequences available in nucleotide databases can be grouped into four clades, corresponding to genes EgAgB1, EgAgB2, EgAgB3 and EgAgB4. In the present study, we use PCR amplifications followed by cloning and sequencing to evaluate the genetic variability for AgB isoforms. Two pairs of primers were independently used for PCR amplification. Both PCR reactions from each of three isolated protoscolex (larvae) were cloned in a plasmid vector and the plasmid inserts of 30 colonies from each cloning experiment were sequenced. Using phylogenetic tools, the 113 EgAgB clones are classified as follows: 25 are related to EgAgB1, 24 to EgAgB2, 9 to EgAgB3 and 39 to EgAgB4. The remaining 16 clones form a separate cluster, which we name EgAgB5, more closely related to EgAgB3 than to any of the other genes. Within each gene group, a number of variant sequences occur, which differ from one another by one or few nucleotides. One EgAgB3 clone has a premature stop codon (pseudogene) and an EgAgB2 clone lacks the region corresponding to the intron. The overall variation cannot be explained by differences among the asexual protoscoleces, or by experimental artifacts. Using Echinococcuss AgB genes from other species/strains as outgroups, neutrality is rejected for EgAgB2, and balancing selection is detected for EgAgB5, which also seems to be involved in gene conversion. We suggest that EgAgB1-EgAgB5 represent a family of contingency genes, that is, genes that are variably expressed, so that some but not others are expressed in each individual parasite. Contingency genes are common in parasitic protozoa and other microparasites, but the EgAgB family is the first set identified in a multicellular parasite.

Comparative analysis of two different subunits of antigen B from Echinococcus granulosus: gene sequences, expression in Escherichia coli and serological evaluation.

Two different Echinococcus granulosus antigen B subunits (AgB8/1 and AgB8/2) were characterized and the structure of the genes encoding these two proteins were compared. DNA sequences were expressed in Escherichia coli and the antigens' diagnostic value was then assessed. The genomic sequence of AgB8/1 has a 92 bp intron in the position corresponding to amino acid 16; the AgB8/2 genomic sequence presents a 68 bp intron in the position corresponding to amino acid 20. Both introns are located between the putative N-terminal hydrophobic sequence and the secreted peptide. A comparison between the AgB8/1 and AgB8/2 nucleotide sequences showed a 53.5% identity among exons and a 50% identity between introns. According to the molecular diversity analysis, the elapsed time since both genes shared a common ancestor would be around 4.2x10(7) years. When the native AgB and the two recombinant antigens (rAgB8/1 and rAgB8/2) were tested in an anti-IgG ELISA, the sensitivity of the native antigen B was 77.41% and its specificity was 81.9%, while rAgB8/1 showed 54.84% of sensitivity and 80.17% of specificity and rAg138/2 had an 83.87% sensitivity and a 98.28% specificity. Statistical analysis confirms that rAgB8/2 has a better performance than rAgB8/1 and native AgB in ELISA.

Livestock trade history, geography, and parasite strains: the mitochondrial genetic structure of Echinococcus granulosus in Argentina.

A sample of 114 isolates of Echinococcus granulosus (Cestoda: Taeniidae) collected from different host species and sites in Argentina has been sequenced for 391 bp from the mitochondrial cytochrome c oxidase subunit I gene to analyze genetic variability and population structure. Nine different haplotypes were identified, 5 of which correspond to already characterized strains. Analysis of molecular variance and nested clade analysis of the distribution of haplotypes among localities within 3 main geographic regions indicate that geographic differentiation accounts for the overall pattern of genetic variability in E. granulosus populations. Significant geographic differentiation is also present when the sheep strain alone is considered. Our results suggest that geographic patterns are not due to actual restricted gene flow between regions but are rather a consequence of past history, probably related to the time and origin of livestock introduction in Argentina.

Redundancy and recombination in the Echinococcus AgB multigene family: is there any similarity with protozoan contingency genes?

Numerous genetic variants of the Echinococcus antigen B (AgB) are encountered within a single metacestode. This could be a reflection of gene redundancy or the result of a somatic hypermutation process. We evaluate the complexity of the AgB multigene family by characterizing the upstream promoter regions of the 4 already known genes (EgAgB1-EgAgB4) and evaluating their redundancy in the genome of 3 Echinococcus species (E. granulosus, E. ortleppi and E. multilocularis) using PCR-based approaches. We have ascertained that the number of AgB gene copies is quite variable, both within and between species. The most repetitive gene seems to be AgB3, of which there are more than 110 copies in E. ortleppi. For E. granulosus, we have cloned and characterized 10 distinct upstream promoter regions of AgB3 from a single metacestode. Our sequences suggest that AgB1 and AgB3 are involved in gene conversion. These results are discussed in light of the role of gene redundancy and recombination in parasite evasion mechanisms of host immunity, which at present are known for protozoan organisms, but virtually unknown for multicellular parasites.

Genetic structure of natural populations of Dryas iulia (Lepidoptera: Nymphalidae) revealed by enzyme polymorphism and mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP).

Dryas iulia appears to have undergone a mode of evolution different from that of other members of its subfamily (Heliconiinae). While other species constitute highly subdivided and inbred populations, those of D. iulia are thought to be large and uniform. Analyzing six samples from Southern Brazil (state of Rio Grande do Sul) in relation to three enzyme systems (EST, LAP, and PGM) and their mtDNA RFLP patterns, we found that they are very similar at the molecular level. The F statistics for enzyme polymorphism data revealed that inbreeding makes a great contribution to the population homozygosity, since FIS equals 0.1322 and FST equals 0.0023. Since the chi-square test showed that FST is not significant, we conclude that all localities belong to the same population. The mtDNA differentiation was about 12 times greater than for nuclear genes; FST was equivalent to 0.0265. We suggest that this difference is due to a higher dispersal of males, in relation to females.

The Echinococcus granulosus antigen B shows a high degree of genetic variability.

Echinococcus granulosus larvae secret a polymeric lipoprotein known as antigen B (AgB) into the metacestode hydatid fluid. Three similar AgB subunits have been previously identified (AgB1, AgB2, and AgB3), and their respective genes isolated, but the actual number of genes encoding AgB subunits remains uncertain. In this study, we characterize the variability of genes encoding the AgB2 subunit, using PCR and RT-PCR followed by cloning and sequencing. We have analyzed 32 cDNA and 34 genomic sequences from a single metacestode, showing a high degree of sequence polymorphism. In addition, we have identified a possibly new AgB subunit, which we call AgB4. Additionally, we describe an AgB2 genomic clone lacking (i) a segment corresponding to the intron and (ii) a short, 45 bp sequence within exon II. The 45 bp segment encompasses the conserved splicing signals and corresponds to a highly conserved insect promoter motif.

Reduced genetic variability within coding and non-coding regions of the Echinococcus multilocularis genome.

Echinococcus multilocularis, a vulpine intestinal tapeworm, is the causative agent of alveolar echinococosis in humans, one of the most severe and lethal parasitic infections in man. To date, there is very little knowledge about the genetical polymorphism of this parasite. To assess sequence polymorphism, we analysed a sample of 33 E. multilocularis isolates from Europe, North America and Asia by PCR-SSCP followed by nucleotide sequencing. This assessment was performed comparatively to sheep, cattle and pig E. granulosus strains. Coding (nuclear antigen B and mitochondrial NADH dehydrogenase genes) and non-coding (introns of actin and homeobox-containing genes) regions of the parasite genome were chosen as targets. Since the estimated nucleotide diversity among genotypes of E. multilocularis were, in general, 10 times lower than among the recognized different strains of E. granulosus, we suggest that the conventional classification of the former species in 2 separated strains (European and North American) should be reviewed.

Selection, recombination and history in a parasitic flatworm (Echinococcus) inferred from nucleotide sequences.

Three species of flatworms from the genus Echinococcus (E. granulosus, E. multilocularis and E. vogeli) and four strains of E. granulosus (cattle, horse, pig and sheep strains) were analysed by the PCR-SSCP method followed by sequencing, using as targets two non-coding and two coding (one nuclear and one mitochondrial) genomic regions. The sequencing data was used to evaluate hypothesis about the parasite breeding system and the causes of genetic diversification. The calculated recombination parameters suggested that cross-fertilisation was rare in the history of the group. However, the relative rates of substitution in the coding sequences showed that positive selection (instead of purifying selection) drove the evolution of an elastase and neutrophil chemotaxis inhibitor gene (AgB/1). The phylogenetic analyses revealed several ambiguities, indicating that the taxonomic status of the E. granulosus horse strain should be revised.

Breeding systems in Echinococcus granulosus (Cestoda; Taeniidae): selfing or outcrossing?

We used the PCR-SSCP method followed by sequencing in order to assess the genetic variability of coding and noncoding parts of the genome of Echinococcus granulosus (Cestoda; Taeniidae) and to test whether or not the parasite populations are mainly self-fertilizing. For this, we analysed a sample of 110 E. granulosus metacestode isolates collected from different geographical regions (Southern Brazil, Europe and Australia) and from different intermediate hosts (ovine, bovine, human, macropod, swine and equine). Using appropriate controls, we were able to identify 4 strains in that sample (sheep, cattle, pig and horse strains). The high degree of genetic differentiation between strains, but not within, and the monomorphism found in most loci (EgAg4, EgActII, EgHbx2 and EgAg6-non-coding-EgAgB/1 and EgND1-coding) indicated that they are largely selfed. On the other hand, outcrossing was also shown to occur, since 5 potential hybrid genotypes between cattle and sheep strains were found in populations of Southern Brazil, but absent in other geographical areas. We suggest that both processes are adaptive. The article also reports, for the first time, the occurrence of the E. granulosus cattle strain in South America.

Isolation and characterization of microsatellites from the tapeworm Echinococcus granulosus.

The Echinococcus granulosus genome was searched for microsatellites using 8 different repeated oligonucleotides as probes (GT15, CT15, AT15, CG15, CAT10, CAA10, CGG10 and CATA10). Southern blot experiments revealed that DNA regions containing GT, CAA, CATA and CT repeats are the most frequent in the E. granulosus genome. AT and CG probes showed no hybridization signal. Two loci containing CA/GT (Egmsca1 and Egmsca2) and 1 locus containing GA/CT (Egmsga1) repeats were cloned and sequenced. The locus Egmsca1 was analysed in 73 isolates from Brazil and Argentina whose strains were previously characterized. Brazilian isolates from cattle strain and Argentinean isolates from camel strain were monomorphic and shared the allele (CA)7. Argentinean isolates of sheep and Tasmanian sheep strains shared 2 alleles [(CA)8 and (CA)10] with Brazilian isolates of sheep strain. The allele (CA)11 was found only in Brazilian isolates of sheep strain at a low frequency. The Brazilian and the Argentinean sheep strain populations were tested for the Hardy-Weinberg equilibrium, and only the former was in agreement with the expectations. No polymorphism was found among individual protoscoleces from a single hydatid cyst, validating the utilization of pooled protoscoleces from 1 cyst, grouped as an isolate, in population studies. This work describes for the first time the isolation and characterization of microsatellites from E. granulosus.

Selection, recombination and history in a parasitic flatworm (Echinococcus) inferred from nucleotide sequences

Three species of flatworms from the genus Echinococcus (E. granulosus, E. multilocularis and E. vogeli) and four strains of E. granulosus (cattle, horse, pig and sheep strains) were analysed by the PCR-SSCP method followed by sequencing, using as targets two non-coding and two coding (one nuclear and one mitochondrial) genomic regions. The sequencing data was used to evaluate hypothesis about the parasite breeding system and the causes of genetic diversification. The calculated recombination parameters suggested that cross-fertilisation was rare in the history of the group. However, the relative rates of substitution in the coding sequences showed that positive selection (instead of purifying selection) drove the evolution of an elastase and neutrophil chemotaxis inhibitor gene (AgB/1). The phylogenetic analysis revealed several ambiguities, indicationg that the taxonomic status of the E. granulosus horse strain should be revised.