A survey of innovation through duplication in the reduced genomes of twelve parasites

Essential Genes of the Parasitic Apicomplexa

Genome-scale mutagenesis screens for genes essential for apicomplexan parasite survival have been completed in three species: Plasmodium falciparum, the major human malaria parasite, Plasmodium berghei, a model rodent malaria parasite, and the more distantly related Toxoplasma gondii, the causative agent of toxoplasmosis. These three species share 2606 single-copy orthologs, 1500 of which have essentiality data in all three screens. In this review, we explore the overlap between these datasets to define the core essential genes of the phylum Apicomplexa. We further discuss the implications of these groundbreaking studies for understanding apicomplexan parasite biology, and we identify promising areas of focus for developing new pan-apicomplexan parasite interventions.

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or ‘dead’ enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

mRNA Polyadenylation Machineries in Intestinal Protozoan Parasites

In humans, mRNA polyadenylation involves the participation of about 20 factors in four main complexes that recognize specific RNA sequences. Notably, CFIm25, CPSF73, and PAP have essential roles for poly(A) site selection, mRNA cleavage, and adenosine residues polymerization. Besides the relevance of polyadenylation for gene expression, information is scarce in intestinal protozoan parasites that threaten human health. To better understand polyadenylation in Entamoeba histolytica, Giardia lamblia, and Cryptosporidium parvum, which represent leading causes of diarrhea worldwide, genomes were screened for orthologs of human factors. Results showed that Entamoeba histolytica and C. parvum have 16 and 12 proteins out of the 19 human proteins used as queries, respectively, while G. lamblia seems to have the smallest polyadenylation machinery with only six factors. Remarkably, CPSF30, CPSF73, CstF77, PABP2, and PAP, which were found in all parasites, could represent the core polyadenylation machinery. Multiple genes were detected for several proteins in Entamoeba, while gene redundancy is lower in Giardia and Cryptosporidium. Congruently with their relevance in the polyadenylation process, CPSF73 and PAP are present in all parasites, and CFIm25 is only missing in Giardia. They conserve the functional domains and predicted folding of human proteins, suggesting they may have the same roles in polyadenylation.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Environmental adaptation of Acanthamoeba castellanii and Entamoeba histolytica at genome level as seen by comparative genomic analysis

Amoebozoans are in many aspects interesting research objects, as they combine features of single-cell organisms with complex signaling and defense systems, comparable to multicellular organisms. Acanthamoeba castellanii is a cosmopolitan species and developed diverged feeding abilities and strong anti-bacterial resistance; Entamoeba histolytica is a parasitic amoeba, who underwent massive gene loss and its genome is almost twice smaller than that of A. castellanii. Nevertheless, both species prosper, demonstrating fitness to their specific environments. Here we compare transcriptomes of A. castellanii and E. histolytica with application of orthologs' search and gene ontology to learn how different life strategies influence genome evolution and restructuring of physiology. A. castellanii demonstrates great metabolic activity and plasticity, while E. histolytica reveals several interesting features in its translational machinery, cytoskeleton, antioxidant protection, and nutritional behavior. In addition, we suggest new features in E. histolytica physiology that may explain its successful colonization of human colon and may facilitate medical research.

Lineage-Specific and Highly Derived Gene Sequences Among Amoebozoa, Revealed by the Comparative Analysis of Transcriptomes from Twelve Amoebozoan Species

Amoebozoa represent a difficult group for traditional morphology-based taxonomy. Molecular approaches, such as gene sequencing and DNA barcoding have greatly enhanced our knowledge of the diversity of these organisms. However, metagenomic studies of Amoebozoa still did not provide as impressive results as they did among some other groups of protists. In environmental DNA surveys done on fragments of SSU rDNA gene and other traditional DNA barcodes, Amoebozoa genes normally constitute a minor part of the total gene diversity and represent only the most abundant lineages. A potential way to resolve this problem is the usage of DNA barcodes based on genes, which are unique or highly derived in this group of organisms. In the present study, we attempted to find such genes and gene families with a low level of paralogy, potentially appropriate as Amoebozoa-specific DNA barcodes. For this we re-assembled transcriptomes of 12 amoebozoan species available from the public databases and performed gene annotation and identification of orthologous genes. In our analysis Amoebozoa-specific and highly derived sequences formed 53,182 clusters of orthologs, containing from 2 to 299 proteins each. Some of these genes may be a potential target for DNA barcoding of Amoebozoa.

Comparative genomic analysis reveals occurrence of genetic recombination in virulent Cryptosporidium hominis subtypes and telomeric gene duplications in Cryptosporidium parvum

Background

Cryptosporidium hominis is a dominant species for human cryptosporidiosis. Within the species, IbA10G2 is the most virulent subtype responsible for all C. hominis–associated outbreaks in Europe and Australia, and is a dominant outbreak subtype in the United States. In recent yearsIaA28R4 is becoming a major new subtype in the United States. In this study, we sequenced the genomes of two field specimens from each of the two subtypes and conducted a comparative genomic analysis of the obtained sequences with those from the only fully sequenced Cryptosporidium parvum genome.

Results

Altogether, 8.59-9.05 Mb of Cryptosporidium sequences in 45–767 assembled contigs were obtained from the four specimens, representing 94.36-99.47% coverage of the expected genome. These genomes had complete synteny in gene organization and 96.86-97.0% and 99.72-99.83% nucleotide sequence similarities to the published genomes of C. parvum and C. hominis, respectively. Several major insertions and deletions were seen between C. hominis and C. parvum genomes, involving mostly members of multicopy gene families near telomeres. The four C. hominis genomes were highly similar to each other and divergent from the reference IaA25R3 genome in some highly polymorphic regions. Major sequence differences among the four specimens sequenced in this study were in the 5′ and 3′ ends of chromosome 6 and the gp60 region, largely the result of genetic recombination.

Conclusions

The sequence similarity among specimens of the two dominant outbreak subtypes and genetic recombination in chromosome 6, especially around the putative virulence determinant gp60 region, suggest that genetic recombination plays a potential role in the emergence of hyper-transmissible C. hominis subtypes. The high sequence conservation between C. parvum and C. hominis genomes and significant differences in copy numbers of MEDLE family secreted proteins and insulinase-like proteases indicate that telomeric gene duplications could potentially contribute to host expansion in C. parvum.

Electronic supplementary material

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