Occurrence of multiple, independent gene fusion events for the fifth and sixth enzymes of pyrimidine biosynthesis in different eukaryotic groups

Delineating transitions during the evolution of specialised peroxisomes: Glycosome formation in kinetoplastid and diplonemid protists

One peculiarity of protists belonging to classes Kinetoplastea and Diplonemea within the phylum Euglenozoa is compartmentalisation of most glycolytic enzymes within peroxisomes that are hence called glycosomes. This pathway is not sequestered in peroxisomes of the third Euglenozoan class, Euglenida. Previous analysis of well-studied kinetoplastids, the ‘TriTryps’ parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp., identified within glycosomes other metabolic processes usually not present in peroxisomes. In addition, trypanosomatid peroxins, i.e. proteins involved in biogenesis of these organelles, are divergent from human and yeast orthologues. In recent years, genomes, transcriptomes and proteomes for a variety of euglenozoans have become available. Here, we track the possible evolution of glycosomes by querying these databases, as well as the genome of Naegleria gruberi, a non-euglenozoan, which belongs to the same protist supergroup Discoba. We searched for orthologues of TriTryps proteins involved in glycosomal metabolism and biogenesis. Predicted cellular location(s) of each metabolic enzyme identified was inferred from presence or absence of peroxisomal-targeting signals. Combined with a survey of relevant literature, we refine extensively our previously postulated hypothesis about glycosome evolution. The data agree glycolysis was compartmentalised in a common ancestor of the kinetoplastids and diplonemids, yet additionally indicates most other processes found in glycosomes of extant trypanosomatids, but not in peroxisomes of other eukaryotes were either sequestered in this ancestor or shortly after separation of the two lineages. In contrast, peroxin divergence is evident in all euglenozoans. Following their gain of pathway complexity, subsequent evolution of peroxisome/glycosome function is complex. We hypothesize compartmentalisation in glycosomes of glycolytic enzymes, their cofactors and subsequently other metabolic enzymes provided selective advantage to kinetoplastids and diplonemids during their evolution in changing marine environments. We contend two specific properties derived from the ancestral peroxisomes were key: existence of nonselective pores for small solutes and the possibility of high turnover by pexophagy. Critically, such pores and pexophagy are characterised in extant trypanosomatids. Increasing amenability of free-living kinetoplastids and recently isolated diplonemids to experimental study means our hypothesis and interpretation of bioinformatic data are suited to experimental interrogation.

Anaerobic peroxisomes in Mastigamoeba balamuthi

The adaptation of eukaryotic cells to anaerobic conditions is reflected by substantial changes to mitochondrial metabolism and functional reduction. Hydrogenosomes belong among the most modified mitochondrial derivative and generate molecular hydrogen concomitant with ATP synthesis. The reduction of mitochondria is frequently associated with loss of peroxisomes, which compartmentalize pathways that generate reactive oxygen species (ROS) and thus protect against cellular damage. The biogenesis and function of peroxisomes are tightly coupled with mitochondria. These organelles share fission machinery components, oxidative metabolism pathways, ROS scavenging activities, and some metabolites. The loss of peroxisomes in eukaryotes with reduced mitochondria is thus not unexpected. Surprisingly, we identified peroxisomes in the anaerobic, hydrogenosome-bearing protist Mastigamoeba balamuthi We found a conserved set of peroxin (Pex) proteins that are required for protein import, peroxisomal growth, and division. Key membrane-associated Pexs (MbPex3, MbPex11, and MbPex14) were visualized in numerous vesicles distinct from hydrogenosomes, the endoplasmic reticulum (ER), and Golgi complex. Proteomic analysis of cellular fractions and prediction of peroxisomal targeting signals (PTS1/PTS2) identified 51 putative peroxisomal matrix proteins. Expression of selected proteins in Saccharomyces cerevisiae revealed specific targeting to peroxisomes. The matrix proteins identified included components of acyl-CoA and carbohydrate metabolism and pyrimidine and CoA biosynthesis, whereas no components related to either β-oxidation or catalase were present. In conclusion, we identified a subclass of peroxisomes, named "anaerobic" peroxisomes that shift the current paradigm and turn attention to the reductive evolution of peroxisomes in anaerobic organisms.

Construction of an oesophageal cancer-specific ceRNA network based on miRNA, lncRNA, and mRNA expression data

AIM

To explore the expression profiles of microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and mRNAs in oesophageal squamous cell carcinoma (ESCC) in order to construct an oesophageal cancer-specific competing endogenous RNA (ceRNA) network.

METHODS

In this work, the expression data of miRNAs, lncRNAs, and mRNAs in ESCC were obtained. An oesophageal cancer-specific ceRNA network was then constructed and investigated.

RESULTS

CeRNAs have the ability to reduce the targeting activity of miRNAs, leading to the de-repression of specific mRNAs with common miRNA response elements. CeRNA interactions have a critical effect in gene regulation and cancer development.

CONCLUSION

This study suggests a novel perspective on potential oesophageal cancer mechanisms as well as novel pathways for modulating ceRNA networks for treating cancers.

An exception among diatoms: unique organization of genes involved in isoprenoid biosynthesis in Rhizosolenia setigera CCMP 1694

The marine diatom Rhizosolenia setigera is unique among this group of microalgae given that it is only one of a handful of diatom species that can produce highly branched isoprenoid (HBI) hydrocarbons. In our efforts to determine distinguishing molecular characteristics in R. setigera CCMP 1694 that could help elucidate the underlying mechanisms for its ability to biosynthesize HBIs, we discovered the occurrence of independent genes encoding for two isopentenyl diphosphate isomerases (RsIDI1 and RsIDI2) and one squalene synthase (RsSQS), enzymes that catalyze non-consecutive steps in isoprenoid biosynthesis. These genes are peculiarly fused in all other genome-sequenced diatoms to date, making their organization in R. setigera CCMP 1694 a clear distinguishing molecular feature. Phylogenetic and sequence analysis of RsIDI1, RsIDI2, and RsSQS revealed that such an arrangement of individually transcribed genes involved in isoprenoid biosynthesis could have arisen through a secondary gene fission event. We further demonstrate that inhibition of squalene synthase (SQS) shifts the flux of exogenous isoprenoid precursors towards HBI biosynthesis suggesting the competition for isoprenoid substrates in the form of farnesyl diphosphate between the sterol and HBI biosynthetic pathways in this diatom.

Bifunctional activity of fused Plasmodium falciparum orotate phosphoribosyltransferase and orotidine 5'-monophosphate decarboxylase

Fusion of the last two enzymes in the pyrimidine biosynthetic pathway in the inversed order by having a COOH-terminal orotate phosphoribosyltransferase (OPRT) and an NH₂-terminal orotidine 5'-monophosphate decarboxylase (OMPDC), as OMPDC-OPRT, are described in many organisms. Here, we produced gene fusions of Plasmodium falciparum OMPDC-OPRT and expressed the bifunctional protein in Escherichia coli. The enzyme was purified to homogeneity using affinity and anion-exchange chromatography, exhibited enzymatic activities and functioned as a dimer. The activities, although unstable, were stabilized by its substrate and product during purification and long-term storage. Furthermore, the enzyme expressed a perfect catalytic efficiency (k_cat/K_m). The k_cat was selectively enhanced up to three orders of magnitude, while the K_m was not much affected and remained at low μM levels when compared to the monofunctional enzymes. The fusion of the two enzymes, creating a "super-enzyme" with perfect catalytic power and more flexibility, reflects cryptic relationship of enzymatic reactivities and metabolic functions on molecular evolution.

Orotidine Monophosphate Decarboxylase--A Fascinating Workhorse Enzyme with Therapeutic Potential

Orotidine 5'-monophosphate decarboxylase (ODCase) is known as one of the most proficient enzymes. The enzyme catalyzes the last reaction step of the de novo pyrimidine biosynthesis, the conversion from orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate. The enzyme is found in all three domains of life, Bacteria, Eukarya and Archaea. Multiple sequence alignment of 750 putative ODCase sequences resulted in five distinct groups. While the universally conserved DxKxxDx motif is present in all the groups, depending on the groups, several characteristic motifs and residues can be identified. Over 200 crystal structures of ODCases have been determined so far. The structures, together with biochemical assays and computational studies, elucidated that ODCase utilized both transition state stabilization and substrate distortion to accelerate the decarboxylation of its natural substrate. Stabilization of the vinyl anion intermediate by a conserved lysine residue at the catalytic site is considered the largest contributing factor to catalysis, while bending of the carboxyl group from the plane of the aromatic pyrimidine ring of OMP accounts for substrate distortion. A number of crystal structures of ODCases complexed with potential drug candidate molecules have also been determined, including with 6-iodo-uridine, a potential antimalarial agent.

A unique insertion of low complexity amino acid sequence underlies protein-protein interaction in human malaria parasite orotate phosphoribosyltransferase and orotidine 5'-monophosphate decarboxylase

OBJECTIVE

To investigate the multienzyme complex formation of human malaria parasite Plasmodium falciparum (P. falciparum) orotate phosphoribosyltransferase (OPRT) and orotidine 5'-monophosphate decarboxylase (OMPDC), the fifth and sixth enzyme of the de novo pyrimidine biosynthetic pathway. Previously, we have clearly established that the two enzymes in the malaria parasite exist physically as a heterotetrameric (OPRT)2(OMPDC)2 complex containing two subunits each of OPRT and OMPDC, and that the complex have catalytic kinetic advantages over the monofunctional enzyme.

METHODS

Both enzymes were cloned and expressed as recombinant proteins. The protein-protein interaction in the enzyme complex was identified using bifunctional chemical cross-linker, liquid chromatography-mass spectrometric analysis and homology modeling.

RESULTS

The unique insertions of low complexity region at the α 2 and α 5 helices of the parasite OMPDC, characterized by single amino acid repeat sequence which was not found in homologous proteins from other organisms, was located on the OPRT-OMPDC interface. The structural models for the protein-protein interaction of the heterotetrameric (OPRT)2(OMPDC)2 multienzyme complex were proposed.

CONCLUSIONS

Based on the proteomic data and structural modeling, it is surmised that the human malaria parasite low complexity region is responsible for the OPRT-OMPDC interaction. The structural complex of the parasite enzymes, thus, represents an efficient functional kinetic advantage, which in line with co-localization principles of evolutional origin, and allosteric control in protein-protein-interactions.

MosaicFinder: identification of fused gene families in sequence similarity networks

MOTIVATION

Gene fusion is an important evolutionary process. It can yield valuable information to infer the interactions and functions of proteins. Fused genes have been identified as non-transitive patterns of similarity in triplets of genes. To be computationally tractable, this approach usually imposes an a priori distinction between a dataset in which fused genes are searched for, and a dataset that may have provided genetic material for fusion. This reduces the 'genetic space' in which fusion can be discovered, as only a subset of triplets of genes is investigated. Moreover, this approach may have a high-false-positive rate, and it does not identify gene families descending from a common fusion event.

RESULTS

We represent similarities between sequences as a network. This leads to an efficient formulation of previous methods of fused gene identification, which we implemented in the Python program FusedTriplets. Furthermore, we propose a new characterization of families of fused genes, as clique minimal separators of the sequence similarity network. This well-studied graph topology provides a robust and fast method of detection, well suited for automatic analyses of big datasets. We implemented this method in the C++ program MosaicFinder, which additionally uses local alignments to discard false-positive candidates and indicates potential fusion points. The grouping into families will help distinguish sequencing or prediction errors from real biological fusions, and it will yield additional insight into the function and history of fused genes.

AVAILABILITY

FusedTriplets and MosaicFinder are published under the GPL license and are freely available with their source code at this address: http://sourceforge.net/projects/mosaicfinder.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phylogeny of parasitic parabasalia and free-living relatives inferred from conventional markers vs. Rpb1, a single-copy gene

BACKGROUND

Parabasalia are single-celled eukaryotes (protists) that are mainly comprised of endosymbionts of termites and wood roaches, intestinal commensals, human or veterinary parasites, and free-living species. Phylogenetic comparisons of parabasalids are typically based upon morphological characters and 18S ribosomal RNA gene sequence data (rDNA), while biochemical or molecular studies of parabasalids are limited to a few axenically cultivable parasites. These previous analyses and other studies based on PCR amplification of duplicated protein-coding genes are unable to fully resolve the evolutionary relationships of parabasalids. As a result, genetic studies of Parabasalia lag behind other organisms.

PRINCIPAL FINDINGS

Comparing parabasalid EF1α, α-tubulin, enolase and MDH protein-coding genes with information from the Trichomonas vaginalis genome reveals difficulty in resolving the history of species or isolates apart from duplicated genes. A conserved single-copy gene encodes the largest subunit of RNA polymerase II (Rpb1) in T. vaginalis and other eukaryotes. Here we directly sequenced Rpb1 degenerate PCR products from 10 parabasalid genera, including several T. vaginalis isolates and avian isolates, and compared these data by phylogenetic analyses. Rpb1 genes from parabasalids, diplomonads, Parabodo, Diplonema and Percolomonas were all intronless, unlike intron-rich homologs in Naegleria, Jakoba and Malawimonas.

CONCLUSIONS/SIGNIFICANCE

The phylogeny of Rpb1 from parasitic and free-living parabasalids, and conserved Rpb1 insertions, support Trichomonadea, Tritrichomonadea, and Hypotrichomonadea as monophyletic groups. These results are consistent with prior analyses of rDNA and GAPDH sequences and ultrastructural data. The Rpb1 phylogenetic tree also resolves species- and isolate-level relationships. These findings, together with the relative ease of Rpb1 isolation, make it an attractive tool for evaluating more extensive relationships within Parabasalia.

Compartmentalization of a glycolytic enzyme in Diplonema, a non-kinetoplastid euglenozoan

Glycosomes are peroxisome-related organelles containing glycolytic enzymes that have been found only in kinetoplastids. We show here that a glycolytic enzyme is compartmentalized in diplonemids, the sister group of kinetoplastids. We found that, similar to kinetoplastid aldolases, the fructose 1,6-bisphosphate aldolase of Diplonema papillatum possesses a type 2-peroxisomal targeting signal. Western blotting showed that this aldolase was present predominantly in the membrane/organellar fraction. Immunofluorescence analysis showed that this aldolase had a scattered distribution in the cytosol, suggesting its compartmentalization. In contrast, orotidine-5'-monophosphate decarboxylase, a non-glycolytic glycosomal enzyme in kinetoplastids, was shown to be a cytosolic enzyme in D. papillatum. Since euglenoids, the earliest diverging branch of Euglenozoa, do not possess glycolytic compartments, these findings suggest that the routing of glycolytic enzymes into peroxisomes may have occurred in a common ancestor of diplonemids and kinetoplastids, followed by diversification of these newly established organelles in each of these euglenozoan lineages.

Multiple horizontal gene transfer events and domain fusions have created novel regulatory and metabolic networks in the oomycete genome

Complex enzymes with multiple catalytic activities are hypothesized to have evolved from more primitive precursors. Global analysis of the Phytophthora sojae genome using conservative criteria for evaluation of complex proteins identified 273 novel multifunctional proteins that were also conserved in P. ramorum. Each of these proteins contains combinations of protein motifs that are not present in bacterial, plant, animal, or fungal genomes. A subset of these proteins were also identified in the two diatom genomes, but the majority of these proteins have formed after the split between diatoms and oomycetes. Documentation of multiple cases of domain fusions that are common to both oomycetes and diatom genomes lends additional support for the hypothesis that oomycetes and diatoms are monophyletic. Bifunctional proteins that catalyze two steps in a metabolic pathway can be used to infer the interaction of orthologous proteins that exist as separate entities in other genomes. We postulated that the novel multifunctional proteins of oomycetes could function as potential Rosetta Stones to identify interacting proteins of conserved metabolic and regulatory networks in other eukaryotic genomes. However ortholog analysis of each domain within our set of 273 multifunctional proteins against 39 sequenced bacterial and eukaryotic genomes, identified only 18 candidate Rosetta Stone proteins. Thus the majority of multifunctional proteins are not Rosetta Stones, but they may nonetheless be useful in identifying novel metabolic and regulatory networks in oomycetes. Phylogenetic analysis of all the enzymes in three pathways with one or more novel multifunctional proteins was conducted to determine the probable origins of individual enzymes. These analyses revealed multiple examples of horizontal transfer from both bacterial genomes and the photosynthetic endosymbiont in the ancestral genome of Stramenopiles. The complexity of the phylogenetic origins of these metabolic pathways and the paucity of Rosetta Stones relative to the total number of multifunctional proteins suggests that the proteome of oomycetes has few features in common with other Kingdoms.

Analysis of rare genomic changes does not support the unikont-bikont phylogeny and suggests cyanobacterial symbiosis as the point of primary radiation of eukaryotes

The deep phylogeny of eukaryotes is an important but extremely difficult problem of evolutionary biology. Five eukaryotic supergroups are relatively well established but the relationship between these supergroups remains elusive, and their divergence seems to best fit a “Big Bang” model. Attempts were made to root the tree of eukaryotes by using potential derived shared characters such as unique fusions of conserved genes. One popular model of eukaryotic evolution that emerged from this type of analysis is the unikont–bikont phylogeny: The unikont branch consists of Metazoa, Choanozoa, Fungi, and Amoebozoa, whereas bikonts include the rest of eukaryotes, namely, Plantae (green plants, Chlorophyta, and Rhodophyta), Chromalveolata, excavates, and Rhizaria. We reexamine the relationships between the eukaryotic supergroups using a genome-wide analysis of rare genomic changes (RGCs) associated with multiple, conserved amino acids (RGC_CAMs and RGC_CAs), to resolve trifurcations of major eukaryotic lineages. The results do not support the basal position of Chromalveolata with respect to Plantae and unikonts or the monophyly of the bikont group and appear to be best compatible with the monophyly of unikonts and Chromalveolata. Chromalveolata show a distinct, additional signal of affinity with Plantae, conceivably, owing to genes transferred from the secondary, red algal symbiont. Excavates are derived forms, with extremely long branches that complicate phylogenetic inference; nevertheless, the RGC analysis suggests that they are significantly more likely to cluster with the unikont–Chromalveolata assemblage than with the Plantae. Thus, the first split in eukaryotic evolution might lie between photosynthetic and nonphotosynthetic forms and so could have been triggered by the endosymbiosis between an ancestral unicellular eukaryote and a cyanobacterium that gave rise to the chloroplast.

EEF2 analysis challenges the monophyly of Archaeplastida and Chromalveolata

BACKGROUND

Classification of eukaryotes provides a fundamental phylogenetic framework for ecological, medical, and industrial research. In recent years eukaryotes have been classified into six major supergroups: Amoebozoa, Archaeplastida, Chromalveolata, Excavata, Opisthokonta, and Rhizaria. According to this supergroup classification, Archaeplastida and Chromalveolata each arose from a single plastid-generating endosymbiotic event involving a cyanobacterium (Archaeplastida) or red alga (Chromalveolata). Although the plastids within members of the Archaeplastida and Chromalveolata share some features, no nucleocytoplasmic synapomorphies supporting these supergroups are currently known.

METHODOLOGY/PRINCIPAL FINDINGS

This study was designed to test the validity of the Archaeplastida and Chromalveolata through the analysis of nucleus-encoded eukaryotic translation elongation factor 2 (EEF2) and cytosolic heat-shock protein of 70 kDa (HSP70) sequences generated from the glaucophyte Cyanophora paradoxa, the cryptophytes Goniomonas truncata and Guillardia theta, the katablepharid Leucocryptos marina, the rhizarian Thaumatomonas sp. and the green alga Mesostigma viride. The HSP70 phylogeny was largely unresolved except for certain well-established groups. In contrast, EEF2 phylogeny recovered many well-established eukaryotic groups and, most interestingly, revealed a well-supported clade composed of cryptophytes, katablepharids, haptophytes, rhodophytes, and Viridiplantae (green algae and land plants). This clade is further supported by the presence of a two amino acid signature within EEF2, which appears to have arisen from amino acid replacement before the common origin of these eukaryotic groups.

CONCLUSIONS/SIGNIFICANCE

Our EEF2 analysis strongly refutes the monophyly of the Archaeplastida and the Chromalveolata, adding to a growing body of evidence that limits the utility of these supergroups. In view of EEF2 phylogeny and other morphological evidence, we discuss the possibility of an alternative eukaryotic supergroup.

Evolutionary analysis of synteny and gene fusion for pyrimidine biosynthetic enzymes in Euglenozoa: an extraordinary gap between kinetoplastids and diplonemids

A unique feature of the genome architecture in the parasitic trypanosomatid protists is large-scale synteny. We addressed the evolutionary trait of synteny in the eukaryotic group, Euglenozoa, which consists of euglenoids (earliest branching), diplonemids, and kinetoplastids (trypanosomatids and bodonids). Synteny of the pyrimidine biosynthetic (pyr) gene cluster, which constitutes part of a large syntenic cluster in trypanosomatids and includes four separate genes (pyr1-pyr4) and one fused gene (pyr6/pyr5 fusion), was conserved in the bodonid, Parabodo caudatus. In the diplonemid, Diplonema papillatum, we identified pyr4 and pyr6 genes. Phylogenetic analyses of pyr4 and pyr6 showed the separate origin of each in kinetoplastids and euglenoids/diplonemids and suggested that kinetoplastids have acquired these genes via lateral gene transfer (LGT). Because replacement of genes by non-orthologs within the syntenic cluster is highly unlikely, we concluded that, after separation of the line leading to diplonemids, the syntenic pyr gene cluster was established in the common ancestor of kinetoplastids, preceded by their acquisition via LGT. Notably, we found that diplonemid pyr6 is a stand-alone gene, inconsistent with both euglenoid pyr5/pyr6 and kinetoplastid pyr6/pyr5 fusions. Our findings provide insights into the evolutionary gaps within Euglenozoa and the evolutionary trait of rearrangement of gene fusion in this lineage.