Remnants of an Ancient Metabolism without Phosphate.

Phosphate is essential for all living systems, serving as a building block of genetic and metabolic machinery. However, it is unclear how phosphate could have assumed these central roles on primordial Earth, given its poor geochemical accessibility. We used systems biology approaches to explore the alternative hypothesis that a protometabolism could have emerged prior to the incorporation of phosphate. Surprisingly, we identified a cryptic phosphate-independent core metabolism producible from simple prebiotic compounds. This network is predicted to support the biosynthesis of a broad category of key biomolecules. Its enrichment for enzymes utilizing iron-sulfur clusters, and the fact that thermodynamic bottlenecks are more readily overcome by thioester rather than phosphate couplings, suggest that this network may constitute a "metabolic fossil" of an early phosphate-free nonenzymatic biochemistry. Our results corroborate and expand previous proposals that a putative thioester-based metabolism could have predated the incorporation of phosphate and an RNA-based genetic system. PAPERCLIP.

Pseudogenization of testis-specific Lfg5 predates human/Neanderthal divergence.

Recent reviews discussed the critical roles of apoptosis in human spermatogenesis and infertility. These reviews highlight the FasL-induced caspase cascade in apoptosis lending importance to our discovery of the pseudogene status of the Lfg5 gene in modern humans, Neanderthal and the Denisovan. This gene is a member of the ancient and highly conserved apoptosis Lifeguard family. This pseudogenization is the result of a premature stop codon at the 3'-end of exon 8 not found in any other ortholog. With the current exception of the domesticated bovine and buffalo, Lfg5's expression in mammals is testis-specific. A full analysis of this gene, its phylogenetic context and its recent hominin changes suggest its inactivation was likely under selection in human evolution.

Use of sustainable chemistry to produce an acyl amino acid surfactant.

Surfactants find wide commercial use as foaming agents, emulsifiers, and dispersants. Currently, surfactants are produced from petroleum, or from seed oils such as palm or coconut oil. Due to concerns with CO(2) emissions and the need to protect rainforests, there is a growing necessity to manufacture these chemicals using sustainable resources In this report, we describe the engineering of a native nonribosomal peptide synthetase pathway (i.e., surfactin synthetase), to generate a Bacillus strain that synthesizes a highly water-soluble acyl amino acid surfactant, rather than the water insoluble lipopeptide surfactin. This novel product has a lower CMC and higher water solubility than myristoyl glutamate, a commercial surfactant. This surfactant is produced by fermentation of cellulosic carbohydrate as feedstock. This method of surfactant production provides an approach to sustainable manufacturing of new surfactants.

The evolution of the cilium and the eukaryotic cell.

The Cilium, the Nucleus and the Mitochondrion are three important organelles whose evolutionary histories are intimately related to the evolution and origin of the eukaryotic cell. The cilium is involved in motility and sensory transduction. The cilium is only found in the eukaryotic cells. Here we show that eight gene duplications prior to the last common ancestor of all extant eukaryotes account for the expansion of the Heavy Chain Dynein family of motor proteins and the evolution of the complexity of the cilium. The ambiguities in the branching of the phylogenetic tree of the HC-Dyneins were resolved by creating well-defined subtrees and using them to create the full tree. Due to the intimate relationship between the nucleus, the division center, mitosis and the basal body/centriole, the evolution of the cilium can now be related to the evolution of mitosis. In addition, the analysis of the cilium rules out its endosymbiotic origin from a phagocytosis of a bacterium.

LFG: a candidate apoptosis regulatory gene family.

The expanding wealth of human, model and other organism's genomic data has allowed the identification of a distinct gene family of apoptotic related genes. Most of these genes are currently unannotated or have been subsumed under two questionably related gene families in the past. For example the transmembrane Bax inhibitor 1 (BI1) motif family has been reported to play a role in apoptosis and to consist of at least seven mammalian protein genes, GRINA, BI1, Lfg/FAIM2, Ghitm, RESC1/Tmbim1, GAAP/Tmbim4, and Tmbm1b. However, a detailed sequence and phylogenetic analysis shows that only five of these form a clear and unique protein family. This now provides information for understanding and investigating the biological roles of these proteins across a wide range of tissues in model organisms. The evolutionary relationships among these genes provide a powerful prospective for extrapolating to human conditions.

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes.

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps.

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

Diversity of WD-repeat proteins.

The WD-repeat-containing proteins form a very large family that is diverse in both its function and domain structure. Within all these proteins the WD-repeat domains are thought to have two common features: the domain folds into a beta propeller; and the domains form a platform without any catalytic activity on which multiple protein complexes assemble reversibly. The fact that these proteins play such key roles in the formation of protein-protein complexes in nearly all the major pathways and organelles unique to eukaryotic cells has two important implications. It supports both their ancient and proto eukaryotic origins and supports a likely association with many genetic diseases.

Origin of the Genetic Code Is Found at the Transition between a Thioester World of Peptides and the Phosphoester World of Polynucleotides.

The early metabolism arising in a Thioester world gave rise to amino acids and their simple peptides. The catalytic activity of these early simple peptides became instrumental in the transition from Thioester World to a Phosphate World. This transition involved the appearances of sugar phosphates, nucleotides, and polynucleotides. The coupling of the amino acids and peptides to nucleotides and polynucleotides is the origin for the genetic code. Many of the key steps in this transition are seen in in the catalytic cores of the nucleotidyltransferases, the class II tRNA synthetases (aaRSs) and the CCA adding enzyme. These catalytic cores are dominated by simple beta hairpin structures formed in the Thioester World. The code evolved from a proto-tRNA a tetramer XCCA interacting with a proto-aminoacyl-tRNA synthetase (aaRS) activating Glycine and Proline, the initial expanded code is found in the acceptor arm of the tRNA, the operational code. It is the coevolution of the tRNA with the aaRSs that is at the heart of the origin and evolution of the genetic code. There is also a close relationship between the accretion models of the evolving tRNA and that of the ribosome.

Transcription factor map alignment of promoter regions.

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels--to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human-mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Evolutionary changes in gene regulation from a comparative analysis of multiple Drosophila species.

Exploiting the ortholog/homolog information now available from the complete genomic sequences of twelve species of Drosophila, we have investigated the ability of regulatory site recognition methods to find regulatory changes for orthologs linked to chromosomal rearrangements. This has made use of the wealth of synteny information among these species. By comparing orthologs in multiple species, we found that the breakpoint of chromosomal rearrangements could have had an impact on regulatory changes of genes next to it with respect to the gene function and location. Extensions of our approach could be used to shed light on the role of gene regulation in the evolutionary adaptation to different environmental conditions.

Gene expansion in Trichomonas vaginalis: a case study on transmembrane cyclases.

The draft genome of Trichomonas vaginalis was recently published, but not much is known on why it has such a large genome. In part this size is due to many gene family expansions. For example we found over 100 members in the adenylyl cyclase family. About half are complete full length genes, and nearly half are initially confirmed to be pseudogenes, the remaining are either incomplete or the apparent result of assembly or sequencing problems. The family can be divided into two subgroups by sequence similarity. These can then be divided into functional and pseudo genes. Among all four of these sets the cyclase domain is very well conserved. We gave three possible hypotheses for that observation: a) Sequencing error or stop-codon read-through; b) Recency of duplication and mutation; c) The likelihood of functional pseudogene.

Constraining ribosomal RNA conformational space.

Despite the potential for many possible secondary-structure conformations, the native sequence of ribosomal RNA (rRNA) is able to find the correct and universally conserved core fold. This study reports a computational analysis investigating two mechanisms that appear to constrain rRNA secondary-structure conformational space: ribosomal proteins and rRNA sequence composition. The analysis was carried out by using rRNA-ribosomal protein interaction data for the Escherichia coli 16S rRNA and free energy minimization software for secondary-structure prediction. The results indicate that selection pressures on rRNA sequence composition and ribosomal protein-rRNA interaction play a key role in constraining the rRNA secondary structure to a single stable form.

Analysis of fracture healing by large-scale transcriptional profile identified temporal relationships between metalloproteinase and ADAMTS mRNA expression.

The aim of this study was to validate the use of transcriptional profiling as a means of characterizing the complex interactions of the thousands of genes that are expressed during fracture healing. Standard mid-diaphyseal tibia fractures were generated in C57/B6 murine tibiae and the transcriptional expression of approximately 13,000 genes was assessed. Three time points after fracture were assessed: day 3, representative of the inflammatory phase; day 10, representative of the peak of cartilage formation; and day 21, representative of the period of primary bone formation and coupled remodeling. A self-organizing mapping approach of the data revealed the temporal relationships between the expression of mRNAs for extracellular matrix proteins and the proteases that degrade the proteoglycan and collagenous matrices. A broad group of extracellular matrix protein mRNAs representative of basement membranes, blood vessels and cartilage all showed elevated expression over the first 21 days of fracture healing. The sorting of the data identified an orderly temporal expression of the metalloproteinases and ADAMTS during the progression of fracture healing with (MMP2/MMP14/TIMP2) and ADAMTS4 and 15 preceding the expression of (MMP9/MMP13). Based on their patterns of expression, relative to the known activities of the encoded proteolytic enzymes, our results suggest that the dissolution of cartilage protoeglycans proceeds before the underlying collagenous components of the matrix are removed. The exclusion of several mRNAs that are normally expressed by osteoclasts in the profiles of mRNAs from days 3 and 10 suggests that osteoclastic activity was largely absent during the early periods of cartilage tissue formation and that proteoglycan and specific collagenase activities, precedes or is prerequisite to later osteoclast infiltration into the remodeling tissues.

Techniques for multi-genome synteny analysis to overcome assembly limitations.

Genome scale synteny analysis, the analysis of relative gene-order conservation between species, can provide key insights into evolutionary chromosomal dynamics, rearrangement rates between species, and speciation analysis. With the rapid availability of multiple genomes, there is a need for efficient solutions to aid in comparative syntenic analysis. Current methods rely on homology assessment and multiple alignment based solutions to determine homologs of genetic markers between species and to infer syntenic relationships. One of the primary challenges facing multi-genome syntenic analysis is the uncertainty posed by genome assembly data with un-sequenced gaps and possible assembly errors. Currently, manual intervention is necessary to tune and correct the results of homology assessment and synteny inference. This paper presents a novel automated approach to overcome some of these limitations. It uses a graph based algorithm to infer sub-graphs denoting synteny chains with the objective of choosing the best locations for homologous elements, in the presence of paralogs, in order to maximize synteny. These synteny chains are expanded by merging sub-graphs based on various user defined thresholds for micro-syntenic scrambling. This approach comprehends and accommodates for contig and scaffold gaps in the assembly to determine homologous genetic elements that might either fall in unsequenced assembly gaps or lie on the edges of sequenced segments or on small fragments. Furthermore, it provides an automated solution for breakpoint analysis and a comparative study of chromosomal rearrangements between species. This approach was applied to a comparative study involving Drosophila.melanogaster and Drosophila.pseudoobscura genomes, as an example, and has been useful in analyzing inter-species syntenic relationships.

Probabilistic prediction of Saccharomyces cerevisiae mRNA 3'-processing sites.

We present a tool for the prediction of mRNA 3'-processing (cleavage and polyadenylation) sites in the yeast Saccharomyces cerevisiae, based on a discrete state-space model or hidden Markov model. Comparison of predicted sites with experimentally verified 3'-processing sites indicates good agreement. All predicted or known yeast genes were analyzed to find probable 3'-processing sites. Known alternative 3'-processing sites, both within the 3'-untranslated region and within the protein coding sequence were successfully identified, leading to the possibility of prediction of previously unknown alternative sites. The lack of an apparent 3'-processing site calls into question the validity of some predicted genes. This is specifically investigated for predicted genes with overlapping coding sequences.

Functional conservation between members of an ancient duplicated transcription factor family, LSF/Grainyhead.

The LSF/Grainyhead transcription factor family is involved in many important biological processes, including cell cycle, cell growth and development. In order to investigate the evolutionary conservation of these biological roles, we have characterized two new family members in Caenorhabditis elegans and Xenopus laevis. The C.elegans member, Ce-GRH-1, groups with the Grainyhead subfamily, while the X.laevis member, Xl-LSF, groups with the LSF subfamily. Ce-GRH-1 binds DNA in a sequence-specific manner identical to that of Drosophila melanogaster Grainyhead. In addition, Ce-GRH-1 binds to sequences upstream of the C.elegans gene encoding aromatic L-amino-acid decarboxylase and genes involved in post-embryonic development, mab-5 and dbl-1. All three C.elegans genes are homologs of D.melanogaster Grainyhead-regulated genes. RNA-mediated interference of Ce-grh-1 results in embryonic lethality in worms, accompanied by soft, defective cuticles. These phenotypes are strikingly similar to those observed previously in D.melanogaster grainyhead mutants, suggesting conservation of the developmental role of these family members over the course of evolution. Our phylogenetic analysis of the expanded LSF/GRH family (including other previously unrecognized proteins/ESTs) suggests that the structural and functional dichotomy of this family dates back more than 700 million years, i.e. to the time when the first multicellular organisms are thought to have arisen.

The evolution of Class II Aminoacyl-tRNA synthetases and the first code.

Class II Aminoacyl-tRNA synthetases are a set of very ancient multi domain proteins. The evolution of the catalytic domain of Class II synthetases can be reconstructed from three peptidyl-hairpins. Further evolution from this primordial catalytic core leads to a split of the Class II synthetases into two divisions potentially associated with the operational code. The earliest form of this code likely coded predominantly Glycine (Gly), Proline (Pro), Alanine (Ala) and "Lysine"/Aspartic acid (Lys/Asp). There is a paradox in these synthetases beginning with a hairpin structure before the Genetic Code existed. A resolution is found in the suggestion that the primordial Aminoacyl synthetases formed in a transition from a Thioester world to a Phosphate ester world.

Triage protein fold prediction.

We have constructed, in a completely automated fashion, a new structure template library for threading that represents 358 distinct SCOP folds where each model is mathematically represented as a Hidden Markov model (HMM). Because the large number of models in the library can potentially dilute the prediction measure, a new triage method for fold prediction is employed. In the first step of the triage method, the most probable structural class is predicted using a set of manually constructed, high-level, generalized structural HMMs that represent seven general protein structural classes: all-alpha, all-beta, alpha/beta, alpha+beta, irregular small metal-binding, transmembrane beta-barrel, and transmembrane alpha-helical. In the second step, only those fold models belonging to the determined structural class are selected for the final fold prediction. This triage method gave more predictions as well as more correct predictions compared with a simple prediction method that lacks the initial classification step. Two different schemes of assigning Bayesian model priors are presented and discussed.

Information-theoretic dissection of pairwise contact potentials.

Pairwise contact potentials have a long, successful history in protein structure prediction. They provide an easily-estimated representation of many attributes of protein structures, such as the hydrophobic effect. In order to improve on existing potentials, one should develop a clear understanding of precisely what information they convey. Here, using mutual information, we quantified the information in amino acid potentials, and the importance of hydropathy, charge, disulfide bonding, and burial. Sampling error in mutual information was controlled for by estimating how much information cannot be attributed to sampling bias. We found the information in amino acid contacts to be modest: 0.04 bits per contact. Of that, only 0.01 bits of information could not be attributed to hydropathy, charge, disulfide bonding, or burial.