Statistical modeling of large microarray data sets to identify stimulus-response profiles

A statistical modeling approach is proposed for use in searching large microarray data sets for genes that have a transcriptional response to a stimulus. The approach is unrestricted with respect to the timing, magnitude or duration of the response, or the overall abundance of the transcript. The statistical model makes an accommodation for systematic heterogeneity in expression levels. Corresponding data analyses provide gene-specific information, and the approach provides a means for evaluating the statistical significance of such information. To illustrate this strategy we have derived a model to depict the profile expected for a periodically transcribed gene and used it to look for budding yeast transcripts that adhere to this profile. Using objective criteria, this method identifies 81% of the known periodic transcripts and 1,088 genes, which show significant periodicity in at least one of the three data sets analyzed. However, only one-quarter of these genes show significant oscillations in at least two data sets and can be classified as periodic with high confidence. The method provides estimates of the mean activation and deactivation times, induced and basal expression levels, and statistical measures of the precision of these estimates for each periodic transcript.

Evolutionary relationships among Rel domains indicate functional diversification by recombination

The recent sequencing of several complete genomes has made it possible to track the evolution of large gene families by their genomic structure. Following the large-scale association of exons encoding domains with well defined functions in invertebrates could be useful in predicting the function of complex multidomain proteins in mammals produced by accretion of domains. With this objective, we have determined the genomic structure of the 14 genes in invertebrates and vertebrates that contain rel domains. The sequence encoding the rel domain is defined by intronic boundaries and has been recombined with at least three structurally and functionally distinct genomic sequences to generate coding sequences for: (i) the rel/Dorsal/NFκB proteins that are retained in the cytoplasm by IkB-like proteins; (ii) the NFATc proteins that sense calcium signals and undergo cytoplasmic-to-nuclear translocation in response to dephosphorylation by calcineurin; and (iii) the TonEBP tonicity-responsive proteins. Remarkably, a single exon in each NFATc family member encodes the entire Ca2+/calcineurin sensing region, including nuclear import/export, calcineurin-binding, and substrate regions. The Rel/Dorsal proteins and the TonEBP proteins are present in Drosophila but not Caenorhabditis elegans. On the other hand, the calcium-responsive NFATc proteins are present only in vertebrates, suggesting that the NFATc family is dedicated to functions specific to vertebrates such as a recombinational immune response, cardiovascular development, and vertebrate-specific aspects of the development and function of the nervous system.

Mirror self-recognition in the bottlenose dolphin: A case of cognitive convergence

The ability to recognize oneself in a mirror is an exceedingly rare capacity in the animal kingdom. To date, only humans and great apes have shown convincing evidence of mirror self-recognition. Two dolphins were exposed to reflective surfaces, and both demonstrated responses consistent with the use of the mirror to investigate marked parts of the body. This ability to use a mirror to inspect parts of the body is a striking example of evolutionary convergence with great apes and humans.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Patterns of variance in stage-structured populations: Evolutionary predictions and ecological implications

Variability in population growth rate is thought to have negative consequences for organism fitness. Theory for matrix population models predicts that variance in population growth rate should be the sum of the variance in each matrix entry times the squared sensitivity term for that matrix entry. I analyzed the stage-specific demography of 30 field populations from 17 published studies for pattern between the variance of a demographic term and its contribution to population growth. There were no instances in which a matrix entry both was highly variable and had a large effect on population growth rate; instead, correlations between estimates of temporal variance in a term and contribution to population growth (sensitivity or elasticity) were overwhelmingly negative. In addition, survivorship or growth sensitivities or elasticities always exceeded those of fecundity, implying that the former two terms always contributed more to population growth rate. These results suggest that variable life history stages tend to contribute relatively little to population growth rates because natural selection may alter life histories to minimize stages with both high sensitivity and high variation.

Evolutionary genomics of Staphylococcus aureus: Insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic

An emerging theme in medical microbiology is that extensive variation exists in gene content among strains of many pathogenic bacterial species. However, this topic has not been investigated on a genome scale with strains recovered from patients with well-defined clinical conditions. Staphylococcus aureus is a major human pathogen and also causes economically important infections in cows and sheep. A DNA microarray representing >90% of the S. aureus genome was used to characterize genomic diversity, evolutionary relationships, and virulence gene distribution among 36 strains of divergent clonal lineages, including methicillin-resistant strains and organisms causing toxic shock syndrome. Genetic variation in S. aureus is very extensive, with ≈22% of the genome comprised of dispensable genetic material. Eighteen large regions of difference were identified, and 10 of these regions have genes that encode putative virulence factors or proteins mediating antibiotic resistance. We find that lateral gene transfer has played a fundamental role in the evolution of S. aureus. The mec gene has been horizontally transferred into distinct S. aureus chromosomal backgrounds at least five times, demonstrating that methicillin-resistant strains have evolved multiple independent times, rather than from a single ancestral strain. This finding resolves a long-standing controversy in S. aureus research. The epidemic of toxic shock syndrome that occurred in the 1970s was caused by a change in the host environment, rather than rapid geographic dissemination of a new hypervirulent strain. DNA microarray analysis of large samples of clinically characterized strains provides broad insights into evolution, pathogenesis, and disease emergence.

Comparison of aquatic food chains using nitrogen isotopes.

Recent studies have shown the utility of delta(15)N to model trophic structure and contaminant bioaccumulation in aquatic food webs. However, cross-system comparisons in delta(15)N can be complicated by differences in delta(15)N at the base of the food chain. Such baseline variation in delta(15)N is difficult to resolve using plankton because of the large temporal variability in the delta(15)N of small organisms that have fast nitrogen turnover. Comparisons using large primary consumers, which have stable tissue isotopic signatures because of their slower nitrogen turnover, show that delta(15)N increases markedly with the human population density in the lake watershed. This shift in delta(15)N likely reflects the high delta(15)N of human sewage. Correcting for this baseline variation in delta(15)N, we report that, contrary to expectations based on previous food-web analysis, the food chains leading up to fish varied by about only one trophic level among the 40 lakes studied. Our results also suggest that the delta(15)N signatures of nitrogen at the base of the food chain will provide a useful tool in the assessment of anthropogenic nutrient inputs.

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

Activation of Jak/STAT proteins involved in signal transduction pathway mediated by receptor for interleukin 2 in malignant T lymphocytes derived from cutaneous anaplastic large T-cell lymphoma and Sezary syndrome.

Signaling through the interleukin 2 receptor (IL-2R) involves phosphorylation of several proteins including Jak3, STAT5, and, in preactivated cells, STAT3. In the present study, we examined the functional status of the IL-2R-associated Jak/STAT pathway in malignant T lymphocytes from advanced skin-based lymphomas: anaplastic large T-cell lymphoma (ALCL) and Sezary syndrome (SzS). Proliferation of three ALCL cell lines (PB-1, 2A, and 2B) was partially inhibited by rapamycin, a blocker of some of the signals mediated by IL-2R, but not by cyclosporin A, FK-506, and prednisone, which suppress signals mediated by the T-cell receptor. All the cell lines expressed on their surface the high-affinity IL-2R (alpha, beta, and gamma c chains). They showed basal, constitutive phosphorylation, and coassociation of Jak3, STAT5, and STAT3. Weak basal phosphorylation of IL-2R gamma c was also detected. In regard to SzS, peripheral blood mononuclear cells from 10 of 14 patients showed basal phosphorylation of Jak3, accompanied by phosphorylation of STAT5 in 9 patients, and STAT3 in 4 patients. However, in vitro overnight culture of SzS cells without exogenous cytokines resulted in markedly decreased Jak3 and STAT5 phosphorylation, which could be reversed by stimulation with IL-2. This indicates that the basal phosphorylation of Jak3 and STAT5 in freshly isolated SzS cells is induced rather than constitutive. The basal activation of the Jak/STAT pathway involved in IL-2R signal transduction in ALCL and SzS cells reported here suggests that this pathway may play a role in the pathogenesis of cutaneous T-cell lymphomas, although the mechanism (induced versus constitutive) may vary between different lymphoma types.

Evolutionary relationships of the coelacanth, lungfishes, and tetrapods based on the 28S ribosomal RNA gene.

The origin of land vertebrates was one of the major transitions in the history of vertebrates. Yet, despite many studies that are based on either morphology or molecules, the phylogenetic relationships among tetrapods and the other two living groups of lobe-finned fishes, the coelacanth and the lungfishes, are still unresolved and debated. Knowledge of the relationships among these lineages, which originated back in the Devonian, has profound implications for the reconstruction of the evolutionary scenario of the conquest of land. We collected the largest molecular data set on this issue so far, about 3,500 base pairs from seven species of the large 28S nuclear ribosomal gene. All phylogenetic analyses (maximum parsimony, neighbor-joining, and maximum likelihood) point toward the hypothesis that lungfishes and coelacanths form a monophyletic group and are equally closely related to land vertebrates. This evolutionary hypothesis complicates the identification of morphological or physiological preadaptations that might have permitted the common ancestor of tetrapods to colonize land. This is because the reconstruction of its ancestral conditions would be hindered by the difficulty to separate uniquely derived characters from shared derived characters in the coelacanth/lungfish and tetrapod lineages. This molecular phylogeny aids in the reconstruction of morphological evolutionary steps by providing a framework; however, only paleontological evidence can determine the sequence of morphological acquisitions that allowed lobe-finned fishes to colonize land.

Antigenic peptides containing large PEG loops designed to extend out of the HLA-A2 binding site form stable complexes with class I major histocompatibility complex molecules.

Recognition of peptides bound to class I major histocompatibility complex (MHC) molecules by specific receptors on T cells regulates the development and activity of the cellular immune system. We have designed and synthesized de novo cyclic peptides that incorporate PEG in the ring structure for binding to class I MHC molecules. The large PEG loops are positioned to extend out of the peptide binding site, thus creating steric effects aimed at preventing the recognition of class I MHC complexes by T-cell receptors. Peptides were synthesized and cyclized on polymer support using high molecular weight symmetrical PEG dicarboxylic acids to link the side chains of lysine residues substituted at positions 4 and 8 in the sequence of the HLA-A2-restricted human T-lymphotrophic virus type I Tax peptide. Cyclic peptides promoted the in vitro folding and assembly of HLA-A2 complexes. Thermal denaturation studies using circular dichroism spectroscopy showed that these complexes are as stable as complexes formed with antigenic peptides.

Initiation process of earthquakes and its implications for seismic hazard reduction strategy.

For the average citizen and the public, "earthquake prediction" means "short-term prediction," a prediction of a specific earthquake on a relatively short time scale. Such prediction must specify the time, place, and magnitude of the earthquake in question with sufficiently high reliability. For this type of prediction, one must rely on some short-term precursors. Examinations of strain changes just before large earthquakes suggest that consistent detection of such precursory strain changes cannot be expected. Other precursory phenomena such as foreshocks and nonseismological anomalies do not occur consistently either. Thus, reliable short-term prediction would be very difficult. Although short-term predictions with large uncertainties could be useful for some areas if their social and economic environments can tolerate false alarms, such predictions would be impractical for most modern industrialized cities. A strategy for effective seismic hazard reduction is to take full advantage of the recent technical advancements in seismology, computers, and communication. In highly industrialized communities, rapid earthquake information is critically important for emergency services agencies, utilities, communications, financial companies, and media to make quick reports and damage estimates and to determine where emergency response is most needed. Long-term forecast, or prognosis, of earthquakes is important for development of realistic building codes, retrofitting existing structures, and land-use planning, but the distinction between short-term and long-term predictions needs to be clearly communicated to the public to avoid misunderstanding.

Sequence scanning: A method for rapid sequence acquisition from large-fragment DNA clones.

A strategy of "sequence scanning" is proposed for rapid acquisition of sequence from clones such as bacteriophage P1 clones, cosmids, or yeast artificial chromosomes. The approach makes use of a special vector, called LambdaScan, that reliably yields subclones with inserts in the size range 8-12 kb. A number of subclones, typically 96 or 192, are chosen at random, and the ends of the inserts are sequenced using vector-specific primers. Then long-range spectrum PCR is used to order and orient the clones. This combination of shotgun and directed sequencing results in a high-resolution physical map suitable for the identification of coding regions or for comparison of sequence organization among genomes. Computer simulations indicate that, for a target clone of 100 kb, the scanning of 192 subclones with sequencing reads as short as 350 bp results in an approximate ratio of 1:2:1 of regions of double-stranded sequence, single-stranded sequence, and gaps. Longer sequencing reads tip the ratio strongly toward increased double-stranded sequence.

Crystal structure of the large fragment of Thermus aquaticus DNA polymerase I at 2.5-A resolution: structural basis for thermostability.

The crystal structure of the large fragment of the Thermus aquaticus DNA polymerase (Klentaq1), determined at 2.5-A resolution, demonstrates a compact two-domain architecture. The C-terminal domain is identical in fold to the equivalent region of the Klenow fragment of Escherichia coli DNA polymerase I (Klenow pol I). Although the N-terminal domain of Klentaq1 differs greatly in sequence from its counterpart in Klenow pol I, it has clearly evolved from a common ancestor. The structure of Klentaq1 reveals the strategy utilized by this protein to maintain activity at high temperatures and provides the structural basis for future improvements of the enzyme.

Evolutionary consequences of food chain length in kelp forest communities.

Kelp forests are strongly influenced by macroinvertebrate grazing on fleshy macroalgae. In the North Pacific Ocean, sea otter predation on macroinvertebrates substantially reduces the intensity of herbivory on macroalgae. Temperate Australasia, in contrast, has no known predator of comparable influence. These ecological and biogeographic patterns led us to predict that (i) the intensity of herbivory should be greater in temperate Australasia than in the North Pacific Ocean; thus (ii) Australasian seaweeds have been under stronger selection to evolve chemical defenses and (iii) Australasian herbivores have been more strongly selected to tolerate these compounds. We tested these predictions first by measuring rates of algal tissue loss to herbivory at several locations in Australasian and North Pacific kelp forests. There were significant differences in grazing rates among sea otter-dominated locations in the North Pacific (0-2% day-1), Australasia (5-7% day-1), and a North Pacific location lacking sea otters (80% day-1). The expectations that chronically high rates of herbivory in Australasia have selected for high concentrations of defensive secondary metabolites (phlorotannins) in brown algae and increased tolerance of these defenses in the herbivores also were supported. Phlorotannin concentrations in kelps and fucoids from Australasia were, on average, 5-6 times higher than those in a comparable suite of North Pacific algae, confirming earlier findings. Furthermore, feeding rates of Australasian herbivores were largely unaffected by phlorotannins, regardless of the compounds' regional source. North Pacific herbivores, in contrast, were consistently deterred by phlorotannins from both Australasia and the North Pacific. These findings suggest that top-level consumers, acting through food chains of various lengths, can strongly influence the ecology and evolution of plantherbivore interactions.