Singular value decomposition for genome-wide expression data processing and modeling

We describe the use of singular value decomposition in transforming genome-wide expression data from genes × arrays space to reduced diagonalized “eigengenes” × “eigenarrays” space, where the eigengenes (or eigenarrays) are unique orthonormal superpositions of the genes (or arrays). Normalizing the data by filtering out the eigengenes (and eigenarrays) that are inferred to represent noise or experimental artifacts enables meaningful comparison of the expression of different genes across different arrays in different experiments. Sorting the data according to the eigengenes and eigenarrays gives a global picture of the dynamics of gene expression, in which individual genes and arrays appear to be classified into groups of similar regulation and function, or similar cellular state and biological phenotype, respectively. After normalization and sorting, the significant eigengenes and eigenarrays can be associated with observed genome-wide effects of regulators, or with measured samples, in which these regulators are overactive or underactive, respectively.

The life sciences Global Image Database (GID)

Although a vast amount of life sciences data is generated in the form of images, most scientists still store images on extremely diverse and often incompatible storage media, without any type of metadata structure, and thus with no standard facility with which to conduct searches or analyses. Here we present a solution to unlock the value of scientific images. The Global Image Database (GID) is a web-based (http://www.g wer.ch/qv/gid/gid.htm) structured central repository for scientific annotated images. The GID was designed to manage images from a wide spectrum of imaging domains ranging from microscopy to automated screening. The annotations in the GID define the source experiment of the images by describing who the authors of the experiment are, when the images were created, the biological origin of the experimental sample and how the sample was processed for visualization. A collection of experimental imaging protocols provides details of the sample preparation, and labeling, or visualization procedures. In addition, the entries in the GID reference these imaging protocols with the probe sequences or antibody names used in labeling experiments. The GID annotations are searchable by field or globally. The query results are first shown as image thumbnail previews, enabling quick browsing prior to original-sized annotated image retrieval. The development of the GID continues, aiming at facilitating the management and exchange of image data in the scientific community, and at creating new query tools for mining image data.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

The observed global warming record: What does it tell us?

Global, near-surface temperature data sets and their derivations are discussed, and differences between the Jones and Intergovernmental Panel on Climate Change data sets are explained. Global-mean temperature changes are then interpreted in terms of anthropogenic forcing influences and natural variability. The inclusion of aerosol forcing improves the fit between modeled and observed changes but does not improve the agreement between the implied climate sensitivity value and the standard model-based range of 1.5–4.5°C equilibrium warming for a CO2 doubling. The implied sensitivity goes from below the model-based range of estimates to substantially above this range. The addition of a solar forcing effect further improves the fit and brings the best-fit sensitivity into the middle of the model-based range. Consistency is further improved when internally generated changes are considered. This consistency, however, hides many uncertainties that surround observed data/model comparisons. These uncertainties make it impossible currently to use observed global-scale temperature changes to narrow the uncertainty range in the climate sensitivity below that estimated directly from climate models.

Global environmental impacts of agricultural expansion: The need for sustainable and efficient practices

The recent intensification of agriculture, and the prospects of future intensification, will have major detrimental impacts on the nonagricultural terrestrial and aquatic ecosystems of the world. The doubling of agricultural food production during the past 35 years was associated with a 6.87-fold increase in nitrogen fertilization, a 3.48-fold increase in phosphorus fertilization, a 1.68-fold increase in the amount of irrigated cropland, and a 1.1-fold increase in land in cultivation. Based on a simple linear extension of past trends, the anticipated next doubling of global food production would be associated with approximately 3-fold increases in nitrogen and phosphorus fertilization rates, a doubling of the irrigated land area, and an 18% increase in cropland. These projected changes would have dramatic impacts on the diversity, composition, and functioning of the remaining natural ecosystems of the world, and on their ability to provide society with a variety of essential ecosystem services. The largest impacts would be on freshwater and marine ecosystems, which would be greatly eutrophied by high rates of nitrogen and phosphorus release from agricultural fields. Aquatic nutrient eutrophication can lead to loss of biodiversity, outbreaks of nuisance species, shifts in the structure of food chains, and impairment of fisheries. Because of aerial redistribution of various forms of nitrogen, agricultural intensification also would eutrophy many natural terrestrial ecosystems and contribute to atmospheric accumulation of greenhouse gases. These detrimental environmental impacts of agriculture can be minimized only if there is much more efficient use and recycling of nitrogen and phosphorus in agroecosystems.

Random circular permutation of genes and expressed polypeptide chains: application of the method to the catalytic chains of aspartate transcarbamoylase.

Recent studies on proteins whose N and C termini are in close proximity have demonstrated that folding of polypeptide chains and assembly of oligomers can be accomplished with circularly permuted chains. As yet no methodical study has been conducted to determine how extensively new termini can be introduced and where such termini cannot be tolerated. We have devised a procedure to generate random circular permutations of the catalytic chains of Escherichia coli aspartate transcarbamoylase (ATCase; EC 2.1.3.2) and to select clones that produce active or stable holoenzyme containing permuted chains. A tandem gene construct was made, based on the desired linkage between amino acid residues in the C- and N-terminal regions of the polypeptide chain, and this DNA was treated with a suitable restriction enzyme to yield a fragment containing the rearranged coding sequence for the chain. Circularization achieved with DNA ligase, followed by linearization at random with DNase I, and incorporation of the linearized, repaired, blunt-ended, rearranged genes into a suitable plasmid permitted the expression of randomly permuted polypeptide chains. The plasmid with appropriate stop codons also contained pyrI, the gene encoding the regulatory chain of ATCase. Colonies expressing detectable amounts of ATCase-like molecules containing permuted catalytic chains were identified by an immunoblot technique or by their ability to grow in the absence of pyrimidines in the growth medium. Sequencing of positive clones revealed a variety of novel circular permutations. Some had N and C termini within helices of the wild-type enzyme as well as deletions and insertions. Permutations were concentrated in the C-terminal domain and only few were detected in the N-terminal domain. The technique, which is adaptable generally to proteins whose N and C termini are near each other, can be of value in relating in vivo folding of nascent, growing polypeptide chains to in vitro renaturation of complete chains and determining the role of protein sequence in folding kinetics.