Nested gene fusions as markers of phylogenetic branchpoints in prokaryotes

The Operon as a Conundrum of Gene Dynamics and Biochemical Constraints: What We Have Learned from Histidine Biosynthesis

Operons represent one of the leading strategies of gene organization in prokaryotes, having a crucial influence on the regulation of gene expression and on bacterial chromosome organization. However, there is no consensus yet on why, how, and when operons are formed and conserved, and many different theories have been proposed. Histidine biosynthesis is a highly studied metabolic pathway, and many of the models suggested to explain operons origin and evolution can be applied to the histidine pathway, making this route an attractive model for the study of operon evolution. Indeed, the organization of his genes in operons can be due to a progressive clustering of biosynthetic genes during evolution, coupled with a horizontal transfer of these gene clusters. The necessity of physical interactions among the His enzymes could also have had a role in favoring gene closeness, of particular importance in extreme environmental conditions. In addition, the presence in this pathway of paralogous genes, heterodimeric enzymes and complex regulatory networks also support other operon evolution hypotheses. It is possible that histidine biosynthesis, and in general all bacterial operons, may result from a mixture of several models, being shaped by different forces and mechanisms during evolution.

Systematic identification and analysis of frequent gene fusion events in metabolic pathways

Background

Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur.

Results

Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/. A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent.

Conclusions

Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

Electronic supplementary material

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

Complex patterns of gene fission in the eukaryotic folate biosynthesis pathway

Shared derived genomic characters can be useful for polarizing phylogenetic relationships, for example, gene fusions have been used to identify deep-branching relationships in the eukaryotes. Here, we report the evolutionary analysis of a three-gene fusion of folB, folK, and folP, which encode enzymes that catalyze consecutive steps in de novo folate biosynthesis. The folK-folP fusion was found across the eukaryotes and a sparse collection of prokaryotes. This suggests an ancient derivation with a number of gene losses in the eukaryotes potentially as a consequence of adaptation to heterotrophic lifestyles. In contrast, the folB-folK-folP gene is specific to a mosaic collection of Amorphea taxa (a group encompassing: Amoebozoa, Apusomonadida, Breviatea, and Opisthokonta). Next, we investigated the stability of this character. We identified numerous gene losses and a total of nine gene fission events, either by break up of an open reading frame (four events identified) or loss of a component domain (five events identified). This indicates that this three gene fusion is highly labile. These data are consistent with a growing body of data indicating gene fission events occur at high relative rates. Accounting for these sources of homoplasy, our data suggest that the folB-folK-folP gene fusion was present in the last common ancestor of Amoebozoa and Opisthokonta but absent in the Metazoa including the human genome. Comparative genomic data of these genes provides an important resource for designing therapeutic strategies targeting the de novo folate biosynthesis pathway of a variety of eukaryotic pathogens such as Acanthamoeba castellanii.

Cohesion group approach for evolutionary analysis of aspartokinase, an enzyme that feeds a branched network of many biochemical pathways

Aspartokinase (Ask) exists within a variable network that supports the synthesis of 9 amino acids and a number of other important metabolites. Lysine, isoleucine, aromatic amino acids, and dipicolinate may arise from the ASK network or from alternative pathways. Ask proteins were subjected to cohesion group analysis, a methodology that sorts a given protein assemblage into groups in which evolutionary continuity is assured. Two subhomology divisions, ASK(alpha) and ASK(beta), have been recognized. The ASK(alpha) subhomology division is the most ancient, being widely distributed throughout the Archaea and Eukarya and in some Bacteria. Within an indel region of about 75 amino acids near the N terminus, ASK(beta) sequences differ from ASK(alpha) sequences by the possession of a proposed ancient deletion. ASK(beta) sequences are present in most Bacteria and usually exhibit an in-frame internal translational start site that can generate a small Ask subunit that is identical to the C-terminal portion of the larger subunit of a heterodimeric unit. Particularly novel are ask genes embedded in gene contexts that imply specialization for ectoine (osmotic agent) or aromatic amino acids. The cohesion group approach is well suited for the easy recognition of relatively recent lateral gene transfer (LGT) events, and many examples of these are described. Given the current density of genome representation for Proteobacteria, it is possible to reconstruct more ancient landmark LGT events. Thus, a plausible scenario in which the three well-studied and iconic Ask homologs of Escherichia coli are not within the vertical genealogy of Gammaproteobacteria, but rather originated via LGT from a Bacteroidetes donor, is supported.

Cohesion group approach for evolutionary analysis of TyrA, a protein family with wide-ranging substrate specificities

Many enzymes and other proteins are difficult subjects for bioinformatic analysis because they exhibit variant catalytic, structural, regulatory, and fusion mode features within a protein family whose sequences are not highly conserved. However, such features reflect dynamic and interesting scenarios of evolutionary importance. The value of experimental data obtained from individual organisms is instantly magnified to the extent that given features of the experimental organism can be projected upon related organisms. But how can one decide how far along the similarity scale it is reasonable to go before such inferences become doubtful? How can a credible picture of evolutionary events be deduced within the vertical trace of inheritance in combination with intervening events of lateral gene transfer (LGT)? We present a comprehensive analysis of a dehydrogenase protein family (TyrA) as a prototype example of how these goals can be accomplished through the use of cohesion group analysis. With this approach, the full collection of homologs is sorted into groups by a method that eliminates bias caused by an uneven representation of sequences from organisms whose phylogenetic spacing is not optimal. Each sufficiently populated cohesion group is phylogenetically coherent and defined by an overall congruence with a distinct section of the 16S rRNA gene tree. Exceptions that occasionally are found implicate a clearly defined LGT scenario whereby the recipient lineage is apparent and the donor lineage of the gene transferred is localized to those organisms that define the cohesion group. Systematic procedures to manage and organize otherwise overwhelming amounts of data are demonstrated.

The TyrA family of aromatic-pathway dehydrogenases in phylogenetic context

BACKGROUND

The TyrA protein family includes members that catalyze two dehydrogenase reactions in distinct pathways leading to L-tyrosine and a third reaction that is not part of tyrosine biosynthesis. Family members share a catalytic core region of about 30 kDa, where inhibitors operate competitively by acting as substrate mimics. This protein family typifies many that are challenging for bioinformatic analysis because of relatively modest sequence conservation and small size.

RESULTS

Phylogenetic relationships of TyrA domains were evaluated in the context of combinatorial patterns of specificity for the two substrates, as well as the presence or absence of a variety of fusions. An interactive tool is provided for prediction of substrate specificity. Interactive alignments for a suite of catalytic-core TyrA domains of differing specificity are also provided to facilitate phylogenetic analysis. tyrA membership in apparent operons (or supraoperons) was examined, and patterns of conserved synteny in relationship to organismal positions on the 16S rRNA tree were ascertained for members of the domain Bacteria. A number of aromatic-pathway genes (hisHb, aroF, aroQ) have fused with tyrA, and it must be more than coincidental that the free-standing counterparts of all of the latter fused genes exhibit a distinct trace of syntenic association.

CONCLUSION

We propose that the ancestral TyrA dehydrogenase had broad specificity for both the cyclohexadienyl and pyridine nucleotide substrates. Indeed, TyrA proteins of this type persist today, but it is also common to find instances of narrowed substrate specificities, as well as of acquisition via gene fusion of additional catalytic domains or regulatory domains. In some clades a qualitative change associated with either narrowed substrate specificity or gene fusion has produced an evolutionary "jump" in the vertical genealogy of TyrA homologs. The evolutionary history of gene organizations that include tyrA can be deduced in genome assemblages of sufficiently close relatives, the most fruitful opportunities currently being in the Proteobacteria. The evolution of TyrA proteins within the broader context of how their regulation evolved and to what extent TyrA co-evolved with other genes as common members of aromatic-pathway regulons is now feasible as an emerging topic of ongoing inquiry.

The origin and evolution of eucaryal HIS7 genes: from metabolon to bifunctional proteins?

The fifth step of histidine biosynthesis is catalysed by an imidazole glycerol-phosphate (IGP) synthase. In Archaea and Bacteria, the active form of IGP synthase is a stable 1:1 dimeric complex constituted by a glutamine amidotransferase (GAT) and a cyclase, the products of hisH and hisF. In Eucarya, the two activities are associated with a single bifunctional polypeptide encoded by HIS7. In this work, we report a comparative analysis of the amino acid sequence of all the available HisH, HisF and HIS7 proteins, which allowed depicting a likely evolutionary pathway leading to the present-day bifunctional HIS7 genes. According to the model that we propose, the bifunctional HIS7 gene is the outcome of a gene fusion event between two independent ancestral cistrons encoding an amidotransferase and a cyclase, respectively. The phylogenetic distribution of the eucaryal HIS7 genes and the analysis of all the available prokaryotic counterparts (hisH and hisF) revealed the absence of such fusions in prokaryotes, suggesting that the fusion event very likely occurred in an early stage of eucaryal evolution and was fixed in the nucleated cells. The biological significance of this gene fusion is also discussed.

Phylogenetic signal in nucleotide data from seed plants: implications for resolving the seed plant tree of life

Effects of taxonomic sampling and conflicting signal on the inference of seed plant trees supported in previous molecular analyses were explored using 13 single-locus data sets. Changing the number of taxa in single-locus analyses had limited effects on log likelihood differences between the gnepine (Gnetales plus Pinaceae) and gnetifer (Gnetales plus conifers) trees. Distinguishing among these trees also was little affected by the use of different substitution parameters. The 13-locus combined data set was partitioned into nine classes based on substitution rates. Sites evolving at intermediate rates had the best likelihood and parsimony scores on gnepine trees, and those evolving at the fastest rates had the best parsimony scores on Gnetales-sister trees (Gnetales plus other seed plants). When the fastest evolving sites were excluded from parsimony analyses, well-supported gnepine trees were inferred from the combined data and from each genomic partition. When all sites were included, Gnetales-sister trees were inferred from the combined data, whereas a different tree was inferred from each genomic partition. Maximum likelihood trees from the combined data and from each genomic partition were well-supported gnepine trees. A preliminary stratigraphic test highlights the poor fit of Gnetales-sister trees to the fossil data.

Ancient origin of the tryptophan operon and the dynamics of evolutionary change

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.

Evolution of the structure and chromosomal distribution of histidine biosynthetic genes

A database of more than 100 histidine biosynthetic genes from different organisms belonging to the three primary domains has been analyzed, including those found in the now completely sequenced genomes of Haemophilus influenzae, Mycoplasma genitalium, Synechocystis sp., Methanococcus jannaschii, and Saccharomyces cerevisiae. The ubiquity of his genes suggests that it is a highly conserved pathway that was probably already present in the last common ancestor of all extant life. The chromosomal distribution of the his genes shows that the enterobacterial histidine operon structure is not the only possible organization, and that there is a diversity of gene arrays for the his pathway. Analysis of the available sequences shows that gene fusions (like those involved in the origin of the Escherichia coli and Salmonella typhimurium hisIE and hisB gene structures) are not universal. In contrast, the elongation event that led to the extant hisA gene from two homologous ancestral modules, as well as the subsequent paralogous duplication that originated hisF, appear to be irreversible and are conserved in all known organisms. The available evidence supports the hypothesis that histidine biosynthesis was assembled by a gene recruitment process.

L-Arogenate Is a Chemoattractant Which Can Be Utilized as the Sole Source of Carbon and Nitrogen by Pseudomonas aeruginosa

L-Arogenate is a commonplace amino acid in nature in consideration of its role as a ubiquitous precursor of L-phenylalanine and/or L-tyrosine. However, the questions of whether it serves as a chemoattractant molecule and whether it can serve as a substrate for catabolism have never been studied. We found that Pseudomonas aeruginosa recognizes L-arogenate as a chemoattractant molecule which can be utilized as a source of both carbon and nitrogen. Mutants lacking expression of either cyclohexadienyl dehydratase or phenylalanine hydroxylase exhibited highly reduced growth rates when utilizing L-arogenate as a nitrogen source. Utilization of L-arogenate as a source of either carbon or nitrogen was dependent upon (sigma)(sup54), as revealed by the use of an rpoN null mutant. The evidence suggests that catabolism of L-arogenate proceeds via alternative pathways which converge at 4-hydroxyphenylpyruvate. In one pathway, prephenate formed in the periplasm by deamination of L-arogenate is converted to 4-hydroxyphenylpyruvate by cyclohexadienyl dehydrogenase. The second route depends upon the sequential action of periplasmic cyclohexadienyl dehydratase, phenylalanine hydroxylase, and aromatic aminotransferase.

Monofunctional chorismate mutase from Serratia rubidaea: a paradigm system for the three-isozyme gene family of enteric bacteria

Serratia rubidaea (ATCC 27614) typifies a substantial number of enteric bacteria which, unlike Escherichia coli, possess a monofunctional species of chorismate mutase (denoted CM-F). CM-F coexists with two additional species of chorismate mutase, each of the latter being one catalytic domain of a bifunctional protein. The two bifunctional proteins are utilized for phenylalanine (CM-P/prephenate dehydratase) and tyrosine (CM-T/cyclohexadienyl dehydrogenase) biosynthesis in all enteric bacteria. S. rubidaea was selected as the organism of choice for purification of CM-F because of the relatively abundant level of expression found for this enzyme. The monofunctional CM-F enzyme was purified about 1600-fold with a yield of about 16%. This is the first monofunctional chorismate mutase to be purified from any gram-negative prokaryote. The CM-F enzyme is a positively charged homodimer made up of 20-kDa subunits. It has a pH optimum of 5.5, exhibits a Km value of 0.33 mM for chorismate, and is sensitive to product inhibition by prephenate that is competitive with respect to chorismate. It is insensitive to feedback inhibition by any of the aromatic amino acids. Partial purification of the bifunctional P-protein and the bifunctional T-protein was also carried out in order to compare the properties of CM-F, CM-P, and CM-T in a common organism. The most striking differential properties of the three isozymes were those of pH optimum and degree of protection conferred by dithiothreitol.

Small single-copy region of plastid DNA in the non-photosynthetic angiosperm Epifagus virginiana contains only two genes. Differences among dicots, monocots and bryophytes in gene organization at a non-bioenergetic locus

We have determined the nucleotide sequence of a 7 kb (1 kb = 10(3) base-pairs) region that includes the entire small single-copy region (SSC) of the plastid genome of Epifagus virginiana, a non-photosynthetic, parasitic flowering plant. The SSC (4.8 kb) is considerably smaller than those of photosynthetic plants due to the complete deletion of all photosynthetic, chlororespiratory and ribosomal protein genes. This leaves only two genes: a protein gene of 1738 codons whose product is unlikely to be involved in bioenergetic processes and a leucine tRNA gene (trn(LUAG)). Both genes span junctions between the inverted repeat and the SSC, with the consequence that the terminal 20 base-pairs of the repeat is transcribed in both directions and functions both as the 3' end of the tRNA gene and as an internal segment of orf1738. We find that the region of tobacco plastid DNA homologous to Epifagus orf1738 contains a single open reading frame (ORF) of 1901 codons rather than the three ORFs of 1244, 273 and 228 codons originally reported. However, we confirm that the equivalent region of the bryophyte Marchantia contains two genes (1068 and 464 codons) corresponding to the N and C-terminal portions of the dicot protein. In contrast, rice plastid DNA contains a severely truncated pseudogene at this locus.