Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets.

Reprogramming of a gene's expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100-500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

Quantitative proteomic analysis reveals concurrent RNA-protein interactions and identifies new RNA-binding proteins in Saccharomyces cerevisiae.

A growing body of evidence supports the existence of an extensive network of RNA-binding proteins (RBPs) whose combinatorial binding affects the post-transcriptional fate of every mRNA in the cell-yet we still do not have a complete understanding of which proteins bind to mRNA, which of these bind concurrently, and when and where in the cell they bind. We describe here a method to identify the proteins that bind to RNA concurrently with an RBP of interest, using quantitative mass spectrometry combined with RNase treatment of affinity-purified RNA-protein complexes. We applied this method to the known RBPs Pab1, Nab2, and Puf3. Our method significantly enriched for known RBPs and is a clear improvement upon previous approaches in yeast. Our data reveal that some reported protein-protein interactions may instead reflect simultaneous binding to shared RNA targets. We also discovered more than 100 candidate RBPs, and we independently confirmed that 77% (23/30) bind directly to RNA. The previously recognized functions of the confirmed novel RBPs were remarkably diverse, and we mapped the RNA-binding region of one of these proteins, the transcriptional coactivator Mbf1, to a region distinct from its DNA-binding domain. Our results also provided new insights into the roles of Nab2 and Puf3 in post-transcriptional regulation by identifying other RBPs that bind simultaneously to the same mRNAs. While existing methods can identify sets of RBPs that interact with common RNA targets, our approach can determine which of those interactions are concurrent-a crucial distinction for understanding post-transcriptional regulation.

Cell-type specific features of circular RNA expression.

Thousands of loci in the human and mouse genomes give rise to circular RNA transcripts; at many of these loci, the predominant RNA isoform is a circle. Using an improved computational approach for circular RNA identification, we found widespread circular RNA expression in Drosophila melanogaster and estimate that in humans, circular RNA may account for 1% as many molecules as poly(A) RNA. Analysis of data from the ENCODE consortium revealed that the repertoire of genes expressing circular RNA, the ratio of circular to linear transcripts for each gene, and even the pattern of splice isoforms of circular RNAs from each gene were cell-type specific. These results suggest that biogenesis of circular RNA is an integral, conserved, and regulated feature of the gene expression program.

The yeast Rab GTPase Ypt1 modulates unfolded protein response dynamics by regulating the stability of HAC1 RNA.

The unfolded protein response (UPR) is a conserved mechanism that mitigates accumulation of unfolded proteins in the ER. The yeast UPR is subject to intricate post-transcriptional regulation, involving recruitment of the RNA encoding the Hac1 transcription factor to the ER and its unconventional splicing. To investigate the mechanisms underlying regulation of the UPR, we screened the yeast proteome for proteins that specifically interact with HAC1 RNA. Protein microarray experiments revealed that HAC1 interacts specifically with small ras GTPases of the Ypt family. We characterized the interaction of HAC1 RNA with one of these proteins, the yeast Rab1 homolog Ypt1. We found that Ypt1 protein specifically associated in vivo with unspliced HAC1 RNA. This association was disrupted by conditions that impaired protein folding in the ER and induced the UPR. Also, the Ypt1-HAC1 interaction depended on IRE1 and ADA5, the two genes critical for UPR activation. Decreasing expression of the Ypt1 protein resulted in a reduced rate of HAC1 RNA decay, leading to significantly increased levels of both unspliced and spliced HAC1 RNA, and delayed attenuation of the UPR, when ER stress was relieved. Our findings establish that Ypt1 contributes to regulation of UPR signaling dynamics by promoting the decay of HAC1 RNA, suggesting a potential regulatory mechanism for linking vesicle trafficking to the UPR and ER homeostasis.

Defining the specificity of cotranslationally acting chaperones by systematic analysis of mRNAs associated with ribosome-nascent chain complexes.

Polypeptides exiting the ribosome must fold and assemble in the crowded environment of the cell. Chaperones and other protein homeostasis factors interact with newly translated polypeptides to facilitate their folding and correct localization. Despite the extensive efforts, little is known about the specificity of the chaperones and other factors that bind nascent polypeptides. To address this question we present an approach that systematically identifies cotranslational chaperone substrates through the mRNAs associated with ribosome-nascent chain-chaperone complexes. We here focused on two Saccharomyces cerevisiae chaperones: the Signal Recognition Particle (SRP), which acts cotranslationally to target proteins to the ER, and the Nascent chain Associated Complex (NAC), whose function has been elusive. Our results provide new insights into SRP selectivity and reveal that NAC is a general cotranslational chaperone. We found surprising differential substrate specificity for the three subunits of NAC, which appear to recognize distinct features within nascent chains. Our results also revealed a partial overlap between the sets of nascent polypeptides that interact with NAC and SRP, respectively, and showed that NAC modulates SRP specificity and fidelity in vivo. These findings give us new insight into the dynamic interplay of chaperones acting on nascent chains. The strategy we used should be generally applicable to mapping the specificity, interplay, and dynamics of the cotranslational protein homeostasis network.

ESRRA-C11orf20 is a recurrent gene fusion in serous ovarian carcinoma.

Every year, ovarian cancer kills approximately 14,000 women in the United States and more than 140,000 women worldwide. Most of these deaths are caused by tumors of the serous histological type, which is rarely diagnosed before it has disseminated. By deep paired-end sequencing of mRNA from serous ovarian cancers, followed by deep sequencing of the corresponding genomic region, we identified a recurrent fusion transcript. The fusion transcript joins the 5' exons of ESRRA, encoding a ligand-independent member of the nuclear-hormone receptor superfamily, to the 3' exons of C11orf20, a conserved but uncharacterized gene located immediately upstream of ESRRA in the reference genome. To estimate the prevalence of the fusion, we tested 67 cases of serous ovarian cancer by RT-PCR and sequencing and confirmed its presence in 10 of these. Targeted resequencing of the corresponding genomic region from two fusion-positive tumor samples identified a nearly clonal chromosomal rearrangement positioning ESRRA upstream of C11orf20 in one tumor, and evidence of local copy number variation in the ESRRA locus in the second tumor. We hypothesize that the recurrent novel fusion transcript may play a role in pathogenesis of a substantial fraction of serous ovarian cancers and could provide a molecular marker for detection of the cancer. Gene fusions involving adjacent or nearby genes can readily escape detection but may play important roles in the development and progression of cancer.

Identification of RNA recognition elements in the Saccharomyces cerevisiae transcriptome.

Post-transcriptional regulation of gene expression, including mRNA localization, translation and decay, is ubiquitous yet still largely unexplored. How is the post-transcriptional regulatory program of each mRNA encoded in its sequence? Hundreds of specific RNA-binding proteins (RBPs) appear to play roles in mediating the post-transcriptional regulatory program, akin to the roles of specific DNA-binding proteins in transcription. As a step toward decoding the regulatory programs encoded in each mRNA, we focused on specific mRNA-protein interactions. We computationally analyzed the sequences of Saccharomyces cerevisiae mRNAs bound in vivo by 29 specific RBPs, identifying eight novel candidate motifs and confirming or extending six earlier reported recognition elements. Biochemical selections for RNA sequences selectively recognized by 12 yeast RBPs yielded novel motifs bound by Pin4, Nsr1, Hrb1, Gbp2, Sgn1 and Mrn1, and recovered the known recognition elements for Puf3, She2, Vts1 and Whi3. Most of the RNA elements we uncovered were associated with coherent mRNA expression changes and were significantly conserved in related yeasts, supporting their functional importance and suggesting that the corresponding RNA-protein interactions are evolutionarily conserved.

Three-dimensional tracking of single mRNA particles in Saccharomyces cerevisiae using a double-helix point spread function.

Optical imaging of single biomolecules and complexes in living cells provides a useful window into cellular processes. However, the three-dimensional dynamics of most important biomolecules in living cells remains essentially uncharacterized. The precise subcellular localization of mRNA-protein complexes plays a critical role in the spatial and temporal control of gene expression, and a full understanding of the control of gene expression requires precise characterization of mRNA transport dynamics beyond the optical diffraction limit. In this paper, we describe three-dimensional tracking of single mRNA particles with 25-nm precision in the x and y dimensions and 50-nm precision in the z dimension in live budding yeast cells using a microscope with a double-helix point spread function. Two statistical methods to detect intermittently confined and directed transport were used to quantify the three-dimensional trajectories of mRNA for the first time, using ARG3 mRNA as a model. Measurements and analysis show that the dynamics of ARG3 mRNA molecules are mostly diffusive, although periods of non-Brownian confinement and directed transport are observed. The quantitative methods detailed in this paper can be broadly applied to the study of mRNA localization and the dynamics of diverse other biomolecules in a wide variety of cell types.

Concordant regulation of translation and mRNA abundance for hundreds of targets of a human microRNA.

MicroRNAs (miRNAs) regulate gene expression posttranscriptionally by interfering with a target mRNA's translation, stability, or both. We sought to dissect the respective contributions of translational inhibition and mRNA decay to microRNA regulation. We identified direct targets of a specific miRNA, miR-124, by virtue of their association with Argonaute proteins, core components of miRNA effector complexes, in response to miR-124 transfection in human tissue culture cells. In parallel, we assessed mRNA levels and obtained translation profiles using a novel global approach to analyze polysomes separated on sucrose gradients. Analysis of translation profiles for approximately 8,000 genes in these proliferative human cells revealed that basic features of translation are similar to those previously observed in rapidly growing Saccharomyces cerevisiae. For approximately 600 mRNAs specifically recruited to Argonaute proteins by miR-124, we found reductions in both the mRNA abundance and inferred translation rate spanning a large dynamic range. The changes in mRNA levels of these miR-124 targets were larger than the changes in translation, with average decreases of 35% and 12%, respectively. Further, there was no identifiable subgroup of mRNA targets for which the translational response was dominant. Both ribosome occupancy (the fraction of a given gene's transcripts associated with ribosomes) and ribosome density (the average number of ribosomes bound per unit length of coding sequence) were selectively reduced for hundreds of miR-124 targets by the presence of miR-124. Changes in protein abundance inferred from the observed changes in mRNA abundance and translation profiles closely matched changes directly determined by Western analysis for 11 of 12 proteins, suggesting that our assays captured most of miR-124-mediated regulation. These results suggest that miRNAs inhibit translation initiation or stimulate ribosome drop-off preferentially near the start site and are not consistent with inhibition of polypeptide elongation, or nascent polypeptide degradation contributing significantly to miRNA-mediated regulation in proliferating HEK293T cells. The observation of concordant changes in mRNA abundance and translational rate for hundreds of miR-124 targets is consistent with a functional link between these two regulatory outcomes of miRNA targeting, and the well-documented interrelationship between translation and mRNA decay.

Diverse RNA-binding proteins interact with functionally related sets of RNAs, suggesting an extensive regulatory system.

RNA-binding proteins (RBPs) have roles in the regulation of many post-transcriptional steps in gene expression, but relatively few RBPs have been systematically studied. We searched for the RNA targets of 40 proteins in the yeast Saccharomyces cerevisiae: a selective sample of the approximately 600 annotated and predicted RBPs, as well as several proteins not annotated as RBPs. At least 33 of these 40 proteins, including three of the four proteins that were not previously known or predicted to be RBPs, were reproducibly associated with specific sets of a few to several hundred RNAs. Remarkably, many of the RBPs we studied bound mRNAs whose protein products share identifiable functional or cytotopic features. We identified specific sequences or predicted structures significantly enriched in target mRNAs of 16 RBPs. These potential RNA-recognition elements were diverse in sequence, structure, and location: some were found predominantly in 3'-untranslated regions, others in 5'-untranslated regions, some in coding sequences, and many in two or more of these features. Although this study only examined a small fraction of the universe of yeast RBPs, 70% of the mRNA transcriptome had significant associations with at least one of these RBPs, and on average, each distinct yeast mRNA interacted with three of the RBPs, suggesting the potential for a rich, multidimensional network of regulation. These results strongly suggest that combinatorial binding of RBPs to specific recognition elements in mRNAs is a pervasive mechanism for multi-dimensional regulation of their post-transcriptional fate.

Autoantigen microarrays for multiplex characterization of autoantibody responses.

We constructed miniaturized autoantigen arrays to perform large-scale multiplex characterization of autoantibody responses directed against structurally diverse autoantigens, using submicroliter quantities of clinical samples. Autoantigen microarrays were produced by attaching hundreds of proteins, peptides and other biomolecules to the surface of derivatized glass slides using a robotic arrayer. Arrays were incubated with patient serum, and spectrally resolvable fluorescent labels were used to detect autoantibody binding to specific autoantigens on the array. We describe and characterize arrays containing the major autoantigens in eight distinct human autoimmune diseases, including systemic lupus erythematosus and rheumatoid arthritis. This represents the first report of application of such technology to multiple human disease sera, and will enable validated detection of antibodies recognizing autoantigens including proteins, peptides, enzyme complexes, ribonucleoprotein complexes, DNA and post-translationally modified antigens. Autoantigen microarrays represent a powerful tool to study the specificity and pathogenesis of autoantibody responses, and to identify and define relevant autoantigens in human autoimmune diseases.

Development of the human infant intestinal microbiota.

Almost immediately after a human being is born, so too is a new microbial ecosystem, one that resides in that person's gastrointestinal tract. Although it is a universal and integral part of human biology, the temporal progression of this process, the sources of the microbes that make up the ecosystem, how and why it varies from one infant to another, and how the composition of this ecosystem influences human physiology, development, and disease are still poorly understood. As a step toward systematically investigating these questions, we designed a microarray to detect and quantitate the small subunit ribosomal RNA (SSU rRNA) gene sequences of most currently recognized species and taxonomic groups of bacteria. We used this microarray, along with sequencing of cloned libraries of PCR-amplified SSU rDNA, to profile the microbial communities in an average of 26 stool samples each from 14 healthy, full-term human infants, including a pair of dizygotic twins, beginning with the first stool after birth and continuing at defined intervals throughout the first year of life. To investigate possible origins of the infant microbiota, we also profiled vaginal and milk samples from most of the mothers, and stool samples from all of the mothers, most of the fathers, and two siblings. The composition and temporal patterns of the microbial communities varied widely from baby to baby. Despite considerable temporal variation, the distinct features of each baby's microbial community were recognizable for intervals of weeks to months. The strikingly parallel temporal patterns of the twins suggested that incidental environmental exposures play a major role in determining the distinctive characteristics of the microbial community in each baby. By the end of the first year of life, the idiosyncratic microbial ecosystems in each baby, although still distinct, had converged toward a profile characteristic of the adult gastrointestinal tract.

The Stanford Tissue Microarray Database.

The Stanford Tissue Microarray Database (TMAD; http://tma.stanford.edu) is a public resource for disseminating annotated tissue images and associated expression data. Stanford University pathologists, researchers and their collaborators worldwide use TMAD for designing, viewing, scoring and analyzing their tissue microarrays. The use of tissue microarrays allows hundreds of human tissue cores to be simultaneously probed by antibodies to detect protein abundance (Immunohistochemistry; IHC), or by labeled nucleic acids (in situ hybridization; ISH) to detect transcript abundance. TMAD archives multi-wavelength fluorescence and bright-field images of tissue microarrays for scoring and analysis. As of July 2007, TMAD contained 205 161 images archiving 349 distinct probes on 1488 tissue microarray slides. Of these, 31 306 images for 68 probes on 125 slides have been released to the public. To date, 12 publications have been based on these raw public data. TMAD incorporates the NCI Thesaurus ontology for searching tissues in the cancer domain. Image processing researchers can extract images and scores for training and testing classification algorithms. The production server uses the Apache HTTP Server, Oracle Database and Perl application code. Source code is available to interested researchers under a no-cost license.

Gene expression programs of human smooth muscle cells: tissue-specific differentiation and prognostic significance in breast cancers.

Smooth muscle is present in a wide variety of anatomical locations, such as blood vessels, various visceral organs, and hair follicles. Contraction of smooth muscle is central to functions as diverse as peristalsis, urination, respiration, and the maintenance of vascular tone. Despite the varied physiological roles of smooth muscle cells (SMCs), we possess only a limited knowledge of the heterogeneity underlying their functional and anatomic specializations. As a step toward understanding the intrinsic differences between SMCs from different anatomical locations, we used DNA microarrays to profile global gene expression patterns in 36 SMC samples from various tissues after propagation under defined conditions in cell culture. Significant variations were found between the cells isolated from blood vessels, bronchi, and visceral organs. Furthermore, pervasive differences were noted within the visceral organ subgroups that appear to reflect the distinct molecular pathways essential for organogenesis as well as those involved in organ-specific contractile and physiological properties. Finally, we sought to understand how this diversity may contribute to SMC-involving pathology. We found that a gene expression signature of the responses of vascular SMCs to serum exposure is associated with a significantly poorer prognosis in human cancers, potentially linking vascular injury response to tumor progression.

Transcriptional modulation of genes encoding structural characteristics of differentiating enterocytes during development of a polarized epithelium in vitro.

Although there is considerable evidence implicating posttranslational mechanisms in the development of epithelial cell polarity, little is known about the patterns of gene expression and transcriptional regulation during this process. We characterized the temporal program of gene expression during cell-cell adhesion-initiated polarization of human Caco-2 cells in tissue culture, which develop structural and functional polarity similar to that of enterocytes in vivo. A distinctive switch in gene expression patterns occurred upon formation of cell-cell contacts between neighboring cells. Expression of genes involved in cell proliferation was down-regulated concomitant with induction of genes necessary for functional specialization of polarized epithelial cells. Transcriptional up-regulation of these latter genes correlated with formation of important structural and functional features in enterocyte differentiation and establishment of structural and functional cell polarity; components of the apical microvilli were induced as the brush border formed during polarization; as barrier function was established, expression of tight junction transmembrane proteins peaked; transcripts encoding components of the apical, but not the basal-lateral trafficking machinery were increased during polarization. Coordinated expression of genes encoding components of functional cell structures were often observed indicating temporal control of expression and assembly of multiprotein complexes.

Parallels between global transcriptional programs of polarizing Caco-2 intestinal epithelial cells in vitro and gene expression programs in normal colon and colon cancer.

Posttranslational mechanisms are implicated in the development of epithelial cell polarity, but little is known about the patterns of gene expression and transcriptional regulation during this process. We characterized temporal patterns of gene expression during cell-cell adhesion-initiated polarization of cultured human Caco-2 cells, which develop structural and functional polarity resembling enterocytes in vivo. A distinctive switch in gene expression patterns occurred upon formation of cell-cell contacts. Comparison to gene expression patterns in normal human colon and colon tumors revealed that the pattern in proliferating, nonpolarized Caco-2 cells paralleled patterns seen in human colon cancer in vivo, including expression of genes involved in cell proliferation. The pattern switched in polarized Caco-2 cells to one more closely resembling that in normal colon tissue, indicating that regulation of transcription underlying Caco-2 cell polarization is similar to that during enterocyte differentiation in vivo. Surprisingly, the temporal program of gene expression in polarizing Caco-2 cells involved changes in signaling pathways (e.g., Wnt, Hh, BMP, FGF) in patterns similar to those during migration and differentiation of intestinal epithelial cells in vivo, despite the absence of morphogen gradients and interactions with stromal cells characteristic of enterocyte differentiation in situ. The full data set is available at http://microarray-pubs.stanford.edu/CACO2.

The preclinical natural history of serous ovarian cancer: defining the target for early detection.

BACKGROUND: Ovarian cancer kills approximately 15,000 women in the United States every year, and more than 140,000 women worldwide. Most deaths from ovarian cancer are caused by tumors of the serous histological type, which are rarely diagnosed before the cancer has spread. Rational design of a potentially life-saving early detection and intervention strategy requires understanding the lesions we must detect in order to prevent lethal progression. Little is known about the natural history of lethal serous ovarian cancers before they become clinically apparent. We can learn about this occult period by studying the unsuspected serous cancers that are discovered in a small fraction of apparently healthy women who undergo prophylactic bilateral salpingo-oophorectomy (PBSO). METHODS AND FINDINGS: We developed models for the growth, progression, and detection of occult serous cancers on the basis of a comprehensive analysis of published data on serous cancers discovered by PBSO in BRCA1 mutation carriers. Our analysis yielded several critical insights into the early natural history of serous ovarian cancer. First, these cancers spend on average more than 4 y as in situ, stage I, or stage II cancers and approximately 1 y as stage III or IV cancers before they become clinically apparent. Second, for most of the occult period, serous cancers are less than 1 cm in diameter, and not visible on gross examination of the ovaries and Fallopian tubes. Third, the median diameter of a serous ovarian cancer when it progresses to an advanced stage (stage III or IV) is about 3 cm. Fourth, to achieve 50% sensitivity in detecting tumors before they advance to stage III, an annual screen would need to detect tumors of 1.3 cm in diameter; 80% detection sensitivity would require detecting tumors less than 0.4 cm in diameter. Fifth, to achieve a 50% reduction in serous ovarian cancer mortality with an annual screen, a test would need to detect tumors of 0.5 cm in diameter. CONCLUSIONS: Our analysis has formalized essential conditions for successful early detection of serous ovarian cancer. Although the window of opportunity for early detection of these cancers lasts for several years, developing a test sufficiently sensitive and specific to take advantage of that opportunity will be a challenge. We estimated that the tumors we would need to detect to achieve even 50% sensitivity are more than 200 times smaller than the clinically apparent serous cancers typically used to evaluate performance of candidate biomarkers; none of the biomarker assays reported to date comes close to the required level of performance. Overcoming the signal-to-noise problem inherent in detection of tiny tumors will likely require discovery of truly cancer-specific biomarkers or development of novel approaches beyond traditional blood protein biomarkers. While this study was limited to ovarian cancers of serous histological type and to those arising in BRCA1 mutation carriers specifically, we believe that the results are relevant to other hereditary serous cancers and to sporadic ovarian cancers. A similar approach could be applied to other cancers to aid in defining their early natural history and to guide rational design of an early detection strategy.

The Stanford Microarray Database: implementation of new analysis tools and open source release of software.

The Stanford Microarray Database (SMD; http://smd.stanford.edu/) is a research tool and archive that allows hundreds of researchers worldwide to store, annotate, analyze and share data generated by microarray technology. SMD supports most major microarray platforms, and is MIAME-supportive and can export or import MAGE-ML. The primary mission of SMD is to be a research tool that supports researchers from the point of data generation to data publication and dissemination, but it also provides unrestricted access to analysis tools and public data from 300 publications. In addition to supporting ongoing research, SMD makes its source code fully and freely available to others under an Open Source license, enabling other groups to create a local installation of SMD. In this article, we describe several data analysis tools implemented in SMD and we discuss features of our software release.

Genome-scale identification of membrane-associated human mRNAs.

The subcellular localization of proteins is critical to their biological roles. Moreover, whether a protein is membrane-bound, secreted, or intracellular affects the usefulness of, and the strategies for, using a protein as a diagnostic marker or a target for therapy. We employed a rapid and efficient experimental approach to classify thousands of human gene products as either "membrane-associated/secreted" (MS) or "cytosolic/nuclear" (CN). Using subcellular fractionation methods, we separated mRNAs associated with membranes from those associated with the soluble cytosolic fraction and analyzed these two pools by comparative hybridization to DNA microarrays. Analysis of 11 different human cell lines, representing lymphoid, myeloid, breast, ovarian, hepatic, colon, and prostate tissues, identified more than 5,000 previously uncharacterized MS and more than 6,400 putative CN genes at high confidence levels. The experimentally determined localizations correlated well with in silico predictions of signal peptides and transmembrane domains, but also significantly increased the number of human genes that could be cataloged as encoding either MS or CN proteins. Using gene expression data from a variety of primary human malignancies and normal tissues, we rationally identified hundreds of MS gene products that are significantly overexpressed in tumors compared to normal tissues and thus represent candidates for serum diagnostic tests or monoclonal antibody-based therapies. Finally, we used the catalog of CN gene products to generate sets of candidate markers of organ-specific tissue injury. The large-scale annotation of subcellular localization reported here will serve as a reference database and will aid in the rational design of diagnostic tests and molecular therapies for diverse diseases.

Anatomic demarcation by positional variation in fibroblast gene expression programs.

Fibroblasts are ubiquitous mesenchymal cells with many vital functions during development, tissue repair, and disease. Fibroblasts from different anatomic sites have distinct and characteristic gene expression patterns, but the principles that govern their molecular specialization are poorly understood. Spatial organization of cellular differentiation may be achieved by unique specification of each cell type; alternatively, organization may arise by cells interpreting their position along a coordinate system. Here we test these models by analyzing the genome-wide gene expression profiles of primary fibroblast populations from 43 unique anatomical sites spanning the human body. Large-scale differences in the gene expression programs were related to three anatomic divisions: anterior-posterior (rostral-caudal), proximal-distal, and dermal versus nondermal. A set of 337 genes that varied according to these positional divisions was able to group all 47 samples by their anatomic sites of origin. Genes involved in pattern formation, cell-cell signaling, and matrix remodeling were enriched among this minimal set of positional identifier genes. Many important features of the embryonic pattern of HOX gene expression were retained in fibroblasts and were confirmed both in vitro and in vivo. Together, these findings suggest that site-specific variations in fibroblast gene expression programs are not idiosyncratic but rather are systematically related to their positional identities relative to major anatomic axes.