Phylogeny, function, and evolution of the cupins, a structurally conserved, functionally diverse superfamily of proteins

The cupin superfamily is a group of functionally diverse proteins that are found in all three kingdoms of life, Archaea, Eubacteria, and Eukaryota. These proteins have a characteristic signature domain comprising two histidine- containing motifs separated by an intermotif region of variable length. This domain consists of six beta strands within a conserved beta barrel structure. Most cupins, such as microbial phosphomannose isomerases (PMIs), AraC- type transcriptional regulators, and cereal oxalate oxidases (OXOs), contain only a single domain, whereas others, such as seed storage proteins and oxalate decarboxylases (OXDCs), are bi-cupins with two pairs of motifs. Although some cupins have known functions and have been characterized at the biochemical level, the majority are known only from gene cloning or sequencing projects. In this study, phylogenetic analyses were conducted on the conserved domain to investigate the evolution and structure/function relationships of cupins, with an emphasis on single- domain plant germin-like proteins (GLPs). An unrooted phylogeny of cupins from a wide spectrum of evolutionary lineages identified three main clusters, microbial PMIs, OXDCs, and plant GLPs. The sister group to the plant GLPs in the global analysis was then used to root a phylogeny of all available plant GLPs. The resulting phylogeny contained three main clades, classifying the GLPs into distinct subfamilies. It is suggested that these subfamilies correlate with functional categories, one of which contains the bifunctional barley germin that has both OXO and superoxide dismutase (SOD) activity. It is proposed that GLPs function primarily as SODs, enzymes that protect plants from the effects of oxidative stress. Closer inspection of the DNA sequence encoding the intermotif region in plant GLPs showed global conservation of thymine in the second codon position, a character associated with hydrophobic residues. Since many of these proteins are multimeric and enzymatically inactive in their monomeric state, this conservation of hydrophobicity is thought to be associated with the need to maintain the various monomer- monomer interactions. The type of structure-based predictive analysis presented in this paper is an important approach for understanding gene function and evolution in an era when genomes from a wide range of organisms are being sequenced at a rapid rate.

Elusive relationships within order fabales: phylogenetic analyses using matK and rbcL sequence data

The order Fabales, including Leguminosae, Polygalaceae, Quillajaceae and Surianaceae, represents a novel hypothesis emerging from angiosperm molecular phylogenies. Despite good support for the order, molecular studies to date have suggested contradictory, poorly supported interfamilial relationships. Our reappraisal of relationships within Fabales addresses past taxon sampling deficiencies, and employs parsimony and Bayesian approaches using sequences from the plastid regions rbcL (166 spp.) and matK (78 spp.). Five alternative hypotheses for interfamilial relationships within Fabales were recovered. The Shimodaira-Hasegawa test found the likelihood of a resolved topology significantly higher than the one calculated for a polytomy, but did not favour any of the alternative hypotheses of relationship within Fabales. In the light of the morphological evidence available and the comparative behavior of rbcL and matK, the topology recovering Polygalaceae as sister to the rest of the order Fabales with Leguminosae more closely related to Quillajaceae + Surianaceae, is considered the most likely hypothesis of interfamilial relationships of the order. Dating of selected crown clades in the Fabales phylogeny using penalized likelihood suggests rapid radiation of the Leguminosae, Polygalaceae, and (Quillajaceae + Surianaceae) crown clades.

A lightweight supercomputing web portal for inferring phylogenetic trees

In this paper we describe a lightweight Web portal developed for running computational jobs on a IBM JS21 Bladecenter cluster, ThamesBlue, for inferring and analyzing evolutionary histories. We first discuss the need for leveraging HPC as a enabler for molecular phylogenetics research. We go on to describe how the portal is designed to interface with existing open-source software that is typical of a HPC resource configuration, and how by design this portal is generic enough to be portable to other similarly configured compute clusters, and for other applications.

Phylogenetic utility of more axillary growth (MAX4 ) - like genes: a case study in Digitalis / Isoplexis (Plantaginaceae)

We present the first assessment of phylogenetic utility of a potential novel low-copy nuclear gene region in flowering plants. A fragment of the MORE AXILLARY GROWTH 4 gene (MAX4, also known as RAMOSUS1 and DECREASED APICAL DOMINANCE1), predicted to span two introns, was isolated from members of Digitalis/Isoplexis. Phylogenetic analyses, under both maximum parsimony and Bayesian inference, were performed and revealed evidence of putative MAX4-like paralogues. The MAX4-like trees were compared with those obtained for Digitalis/Isoplexis using ITS and trnL-F, revealing a high degree of incongruence between these different DNA regions. Network analyses indicate complex patterns of evolution between the MAX4 sequences, which cannot be adequately represented on bifurcating trees. The incidence of paralogy restricts the use of MAX4 in phylogenetic inference within the study group, although MAX4 could potentially be used in combination with other DNA regions for resolving species relationships in cases where paralogues can be clearly identified.

Phylogenetic utility of MORE AXILLARY GROWTH4 (MAX4)-like genes: a case study in Digitalis/Isoplexis (Plantaginaceae)

We present the first assessment of phylogenetic utility of a potential novel low-copy nuclear gene region in flowering plants. A fragment of the MORE AXILLARY GROWTH 4 gene (MAX4, also known as RAMOSUS1 and DECREASED APICAL DOMINANCE1), predicted to span two introns, was isolated from members of Digitalis/Isoplexis. Phylogenetic analyses, under both maximum parsimony and Bayesian inference, were performed and revealed evidence of putative MAX4-like paralogues. The MAX4-like trees were compared with those obtained for Digitalis/Isoplexis using ITS and trnL-F, revealing a high degree of incongruence between these different DNA regions. Network analyses indicate complex patterns of evolution between the MAX4 sequences, which cannot be adequately represented on bifurcating trees. The incidence of paralogy restricts the use of MAX4 in phylogenetic inference within the study group, although MAX4 could potentially be used in combination with other DNA regions for resolving species relationships in cases where paralogues can be clearly identified.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Phylogenetic mixture models can reduce node-density artifacts

We investigate the performance of phylogenetic mixture models in reducing a well-known and pervasive artifact of phylogenetic inference known as the node-density effect, comparing them to partitioned analyses of the same data. The node-density effect refers to the tendency for the amount of evolutionary change in longer branches of phylogenies to be underestimated compared to that in regions of the tree where there are more nodes and thus branches are typically shorter. Mixture models allow more than one model of sequence evolution to describe the sites in an alignment without prior knowledge of the evolutionary processes that characterize the data or how they correspond to different sites. If multiple evolutionary patterns are common in sequence evolution, mixture models may be capable of reducing node-density effects by characterizing the evolutionary processes more accurately. In gene-sequence alignments simulated to have heterogeneous patterns of evolution, we find that mixture models can reduce node-density effects to negligible levels or remove them altogether, performing as well as partitioned analyses based on the known simulated patterns. The mixture models achieve this without knowledge of the patterns that generated the data and even in some cases without specifying the full or true model of sequence evolution known to underlie the data. The latter result is especially important in real applications, as the true model of evolution is seldom known. We find the same patterns of results for two real data sets with evidence of complex patterns of sequence evolution: mixture models substantially reduced node-density effects and returned better likelihoods compared to partitioning models specifically fitted to these data. We suggest that the presence of more than one pattern of evolution in the data is a common source of error in phylogenetic inference and that mixture models can often detect these patterns even without prior knowledge of their presence in the data. Routine use of mixture models alongside other approaches to phylogenetic inference may often reveal hidden or unexpected patterns of sequence evolution and can improve phylogenetic inference.

Floral meristem indeterminacy depends on flower position and is facilitated by acarpellate gynoecium development in Impatiens balsamina

Floral meristems are generally determinate. Termination of their activity varies with species, occurring after carpel or ovule development, depending on the placentation type. In terminal flowering Impatiens balsamina (cv. Dwarf Bush Flowered) some flowers exhibit meristem indeterminacy; they produce organs from the placenta after ovule development. Here we provide a detailed description of gynoecium development in this line and explore the basis of the indeterminate nature of some of its floral meristems. We find that the placenta is sometimes established without complete carpel fusion. Proliferative growth derives from meristematic remnants of the placenta and is more common in the terminal inflorescence. RNA in situ hybridization reveals that IbLFY (Impatiens LFY homologue) is expressed in all meristem states, even in proliferating meristems. Expression of IbAG in axillary flowers is as expected in the meristem, stamens and carpels but absent from the proliferating meristem. We conclude that I. balsamina has cauline placentation. Incomplete suppression of inflorescence identity in flowers of the terminal inflorescence leads to floral meristem proliferation after ovule development in this species.

Degradation pathway of bisphenol A: Does ipso substitution apply to phenols containing a quaternary alpha-carbon structure in the para position?

The degradation of bisphenol A and nonylphenol involves the unusual rearrangement of stable carboncarbon bonds. Some nonylphenol isomers and bisphenol A possess a quaternary alpha-carbon atom as a common structural feature. The degradation of nonylphenol in Sphingomonas sp. strain TTNP3 occurs via a type II ipso substitution with the presence of a quaternary alpha-carbon as a prerequisite. We report here a new degradation pathway of bisphenol A. Consequent to the hydroxylation at position C-4, according to a type 11 ipso substitution mechanism, the C-C bond between the phenolic moiety and the isopropyl group of bisphenol A is broken. Besides the formation of hydroquinone and 4-(2-hydroxypropan-2-yl) phenol as the main metabolites, further compounds resulting from molecular rearrangements consistent with a carbocationic intermediate were identified. Assays with resting cells or cell extracts of Sphingomonas sp. strain TTNP3 under an 18 02 atmosphere were performed. One atom of 180, was present in hydroquinone, resulting from the monooxygenation of bisphenol A and nonylphenol. The monooxygenase activity was dependent on both NADPH and flavin adenine dinucleotide. Various cytochrome P450 inhibitors had identical inhibition effects on the conversion of both xenobiotics. Using a mutant of Sphingomonas sp. strain TTNP3, which is defective for growth on nonylphenol, we demonstrated that the reaction is catalyzed by the same enzymatic system. In conclusion, the degradation of bisphenol A and nonylphenol is initiated by the same monooxygenase, which may also lead to ipso substitution in other xenobiotics containing phenol with a quaternary a-carbon.

Phylogenetic analysis of the pPT23A plasmid family of Pseudomonas syringae

The pPT23A plasmid family of Pseudomonas syringae contains members that contribute to the ecological and pathogenic fitness of their P. syringae hosts. In an effort to understand the evolution of these plasmids and their hosts, we undertook a comparative analysis of the phylogeny of plasmid genes and that of conserved chromosomal genes from P. syringae. In total, comparative sequence and phylogenetic analyses were done utilizing 47 pPT23A family plasmids (PFPs) from 16 pathovars belonging to six genomospecies. Our results showed that the plasmid replication gene (repA), the only gene currently known to be distributed among all the PFPs, had a phylogeny that was distinct from that of the P. syringae hosts of these plasmids and from those of other individual genes on PFPs. The phylogenies of two housekeeping chromosomal genes, those for DNA gyrase B subunit (gyrB) and primary sigma factor (rpoD), however, were strongly associated with genomospecies of P. syringae. Based on the results from this study, we conclude that the pPT23A plasmid family represents a dynamic genome that is mobile among P. syringae pathovars.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

A phylogenetic reappraisal of the Peltophorum group (Caesalpinieae: Leguminosae) based on the chloroplast trnL-F, rbcL and rps16 sequence data

The monophyly of the Peltophorum group, one of nine informal groups recognized by Polhill in the Caesalpinieae, was tested using sequence data from the trnL-F, rbcL, and rps16 regions of the chloroplast genome. Exemplars were included from all 16 genera of the Peltophorum group, and from 15 genera representing seven of the other eight informal groups in the tribe. The data were analyzed separately and in combined analyses using parsimony and Bayesian methods. The analysis method had little effect on the topology of well-supported relationships. The molecular data recovered a generally well-supported phylogeny with many intergeneric relationships resolved. Results show that the Peltophorum group as currently delimited is polyphyletic, but that eight genera plus one undescribed genus form a core Peltophorum group, which is referred to here as the Peltophorum group sensu stricto. These genera are Bussea, Conzattia, Colvillea, Delonix, Heteroflorum (inedit.), Lemuropisum, Parkinsonia, Peltophorum, and Schizolobium. The remaining eight genera of the Peltophorum group s.l. are distributed across the Caesalpinieae. Morphological support for the redelimited Peltophorum group and the other recovered clades was assessed, and no unique synapomorphy was found for the Peltophorum group s.s. A proposal for the reclassification of the Peltophorum group s.l. is presented.