Cluster analysis and display of genome-wide expression patterns

A system of cluster analysis for genome-wide expression data from DNA microarray hybridization is described that uses standard statistical algorithms to arrange genes according to similarity in pattern of gene expression. The output is displayed graphically, conveying the clustering and the underlying expression data simultaneously in a form intuitive for biologists. We have found in the budding yeast Saccharomyces cerevisiae that clustering gene expression data groups together efficiently genes of known similar function, and we find a similar tendency in human data. Thus patterns seen in genome-wide expression experiments can be interpreted as indications of the status of cellular processes. Also, coexpression of genes of known function with poorly characterized or novel genes may provide a simple means of gaining leads to the functions of many genes for which information is not available currently.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Proteome Analysis Database: online application of InterPro and CluSTr for the functional classification of proteins in whole genomes

The SWISS-PROT group at EBI has developed the Proteome Analysis Database utilising existing resources and providing comparative analysis of the predicted protein coding sequences of the complete genomes of bacteria, archaea and eukaryotes (http://www.ebi.ac.uk/proteome/). The two main projects used, InterPro and CluSTr, give a new perspective on families, domains and sites and cover 31–67% (InterPro statistics) of the proteins from each of the complete genomes. CluSTr covers the three complete eukaryotic genomes and the incomplete human genome data. The Proteome Analysis Database is accompanied by a program that has been designed to carry out InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database.

Analysis of endoplasmic reticulum trafficking signals by combinatorial screening in mammalian cells

To improve the accuracy of predicting membrane protein sorting signals, we developed a general methodology for defining trafficking signal consensus sequences in the environment of the living cell. Our approach uses retroviral gene transfer to create combinatorial expression libraries of trafficking signal variants in mammalian cells, flow cytometry to sort cells based on trafficking phenotype, and quantitative trafficking assays to measure the efficacy of individual signals. Using this strategy to analyze arginine- and lysine-based endoplasmic reticulum localization signals, we demonstrate that small changes in the local sequence context dramatically alter signal strength, generating a broad spectrum of trafficking phenotypes. Finally, using sequences from our screen, we found that the potency of di-lysine, but not di-arginine, mediated endoplasmic reticulum localization was correlated with the strength of interaction with α-COP.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

Analysis of gene expression during myc oncogene-induced lymphomagenesis in the bursa of Fabricius

The transcriptional effects of deregulated myc gene overexpression are implicated in tumorigenesis in a spectrum of experimental and naturally occurring neoplasms. In follicles of the chicken bursa of Fabricius, myc induction of B-cell neoplasia requires a target cell population present during early bursal development and progresses through preneoplastic transformed follicles to metastatic lymphomas. We developed a chicken immune system cDNA microarray to analyze broad changes in gene expression that occur during normal embryonic B-cell development and during myc-induced neoplastic transformation in the bursa. The number of mRNAs showing at least 3-fold change was greater during myc-induced lymphomagenesis than during normal development, and hierarchical cluster analysis of expression patterns revealed that levels of several hundred mRNAs varied in concert with levels of myc overexpression. A set of 41 mRNAs were most consistently elevated in myc-overexpressing preneoplastic and neoplastic cells, most involved in processes thought to be subject to regulation by Myc. The mRNAs for another cluster of genes were overexpressed in neoplasia independent of myc expression level, including a small subset with the expression signature of embryonic bursal lymphocytes. Overexpression of myc, and some of the genes overexpressed with myc, may be important for generation of preneoplastic transformed follicles. However, expression profiles of late metastatic tumors showed a large variation in concert with myc expression levels, and some showed minimal myc overexpression. Therefore, high-level myc overexpression may be more important in the early induction of these lymphomas than in maintenance of late-stage metastases.

Does Leaf Position within a Canopy Affect Acclimation of Photosynthesis to Elevated CO2?1 : Analysis of a Wheat Crop under Free-Air CO2 Enrichment

Previous studies of photosynthetic acclimation to elevated CO2 have focused on the most recently expanded, sunlit leaves in the canopy. We examined acclimation in a vertical profile of leaves through a canopy of wheat (Triticum aestivum L.). The crop was grown at an elevated CO2 partial pressure of 55 Pa within a replicated field experiment using free-air CO2 enrichment. Gas exchange was used to estimate in vivo carboxylation capacity and the maximum rate of ribulose-1,5-bisphosphate-limited photosynthesis. Net photosynthetic CO2 uptake was measured for leaves in situ within the canopy. Leaf contents of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco), light-harvesting-complex (LHC) proteins, and total N were determined. Elevated CO2 did not affect carboxylation capacity in the most recently expanded leaves but led to a decrease in lower, shaded leaves during grain development. Despite this acclimation, in situ photosynthetic CO2 uptake remained higher under elevated CO2. Acclimation at elevated CO2 was accompanied by decreases in both Rubisco and total leaf N contents and an increase in LHC content. Elevated CO2 led to a larger increase in LHC/Rubisco in lower canopy leaves than in the uppermost leaf. Acclimation of leaf photosynthesis to elevated CO2 therefore depended on both vertical position within the canopy and the developmental stage.

Analysis of genetic instability during mammary tumor progression using a novel selection-based assay for in vivo mutations in a bacteriophage lambda transgene target.

Genetic instability is thought to be responsible for the numerous genotypic changes that occur during neoplastic transformation and metastatic progression. To explore the role of genetic instability at the level of point mutations during mammary tumor development and malignant progression, we combined transgenic mouse models of mutagenesis detection and oncogenesis. Bitransgenic mice were generated that carried both a bacteriophage lambda transgene to assay mutagenesis and a polyomavirus middle T oncogene, mammary gland-targeted expression of which led to metastatic mammary adenocarcinomas. We developed a novel assay for the detection of mutations in the lambda transgene that selects for phage containing forward mutations only in the lambda cII gene, using an hfl- bacterial host. In addition to the relative ease of direct selection, the sensitivity of this assay for both spontaneous and chemically induced mutations was comparable to the widely used mutational target gene, lambda lacI, making the cII assay an attractive alternative for mutant phage recovery for any lambda-based mouse mutagenesis assay system. The frequencies of lambda cII- mutants were not significantly different in normal mammary epithelium, primary mammary adenocarcinomas, and pulmonary metastases. The cII mutational spectra in these tissues consisted mostly of G/C-->A/T transitions, a large fraction of which occurred at CpG dinucleotides. These data suggest that, in this middle T oncogene model of mammary tumor progression, a significant increase in mutagenesis is not required for tumor development or for metastatic progression.