101
|
Sahoo D, Dill DL, Gentles AJ, Tibshirani R, Plevritis SK. Boolean implication networks derived from large scale, whole genome microarray datasets. Genome Biol 2008; 9:R157. [PMID: 18973690 PMCID: PMC2760884 DOI: 10.1186/gb-2008-9-10-r157] [Citation(s) in RCA: 80] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2008] [Revised: 09/06/2008] [Accepted: 10/30/2008] [Indexed: 11/23/2022] Open
Abstract
A method for analysis of microarray data is presented that extracts statistically significant Boolean implication relationships between pairs of genes. We describe a method for extracting Boolean implications (if-then relationships) in very large amounts of gene expression microarray data. A meta-analysis of data from thousands of microarrays for humans, mice, and fruit flies finds millions of implication relationships between genes that would be missed by other methods. These relationships capture gender differences, tissue differences, development, and differentiation. New relationships are discovered that are preserved across all three species.
Collapse
Affiliation(s)
- Debashis Sahoo
- Department of Computer Science, Stanford University, Stanford, CA 94305, USA
| | | | | | | | | |
Collapse
|
102
|
Abstract
It has been reported that increasingly microRNAs are associated with diseases. However, the patterns among the microRNA-disease associations remain largely unclear. In this study, in order to dissect the patterns of microRNA-disease associations, we performed a comprehensive analysis to the human microRNA-disease association data, which is manually collected from publications. We built a human microRNA associated disease network. Interestingly, microRNAs tend to show similar or different dysfunctional evidences for the similar or different disease clusters, respectively. A negative correlation between the tissue-specificity of a microRNA and the number of diseases it associated was uncovered. Furthermore, we observed an association between microRNA conservation and disease. Finally, we uncovered that microRNAs associated with the same disease tend to emerge as predefined microRNA groups. These findings can not only provide help in understanding the associations between microRNAs and human diseases but also suggest a new way to identify novel disease-associated microRNAs.
Collapse
|
103
|
Functional organization of the transcriptome in human brain. Nat Neurosci 2008; 11:1271-82. [PMID: 18849986 DOI: 10.1038/nn.2207] [Citation(s) in RCA: 558] [Impact Index Per Article: 34.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2008] [Accepted: 09/09/2008] [Indexed: 01/19/2023]
Abstract
The enormous complexity of the human brain ultimately derives from a finite set of molecular instructions encoded in the human genome. These instructions can be directly studied by exploring the organization of the brain's transcriptome through systematic analysis of gene coexpression relationships. We analyzed gene coexpression relationships in microarray data generated from specific human brain regions and identified modules of coexpressed genes that correspond to neurons, oligodendrocytes, astrocytes and microglia. These modules provide an initial description of the transcriptional programs that distinguish the major cell classes of the human brain and indicate that cell type-specific information can be obtained from whole brain tissue without isolating homogeneous populations of cells. Other modules corresponded to additional cell types, organelles, synaptic function, gender differences and the subventricular neurogenic niche. We found that subventricular zone astrocytes, which are thought to function as neural stem cells in adults, have a distinct gene expression pattern relative to protoplasmic astrocytes. Our findings provide a new foundation for neurogenetic inquiries by revealing a robust and previously unrecognized organization to the human brain transcriptome.
Collapse
|
104
|
Wolf MY, Wolf YI, Koonin EV. Comparable contributions of structural-functional constraints and expression level to the rate of protein sequence evolution. Biol Direct 2008; 3:40. [PMID: 18840284 PMCID: PMC2572155 DOI: 10.1186/1745-6150-3-40] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2008] [Accepted: 10/07/2008] [Indexed: 01/01/2023] Open
Abstract
Background Proteins show a broad range of evolutionary rates. Understanding the factors that are responsible for the characteristic rate of evolution of a given protein arguably is one of the major goals of evolutionary biology. A long-standing general assumption used to be that the evolution rate is, primarily, determined by the specific functional constraints that affect the given protein. These constrains were traditionally thought to depend both on the specific features of the protein's structure and its biological role. The advent of systems biology brought about new types of data, such as expression level and protein-protein interactions, and unexpectedly, a variety of correlations between protein evolution rate and these variables have been observed. The strongest connections by far were repeatedly seen between protein sequence evolution rate and the expression level of the respective gene. It has been hypothesized that this link is due to the selection for the robustness of the protein structure to mistranslation-induced misfolding that is particularly important for highly expressed proteins and is the dominant determinant of the sequence evolution rate. Results This work is an attempt to assess the relative contributions of protein domain structure and function, on the one hand, and expression level on the other hand, to the rate of sequence evolution. To this end, we performed a genome-wide analysis of the effect of the fusion of a pair of domains in multidomain proteins on the difference in the domain-specific evolutionary rates. The mistranslation-induced misfolding hypothesis would predict that, within multidomain proteins, fused domains, on average, should evolve at substantially closer rates than the same domains in different proteins because, within a mutlidomain protein, all domains are translated at the same rate. We performed a comprehensive comparison of the evolutionary rates of mammalian and plant protein domains that are either joined in multidomain proteins or contained in distinct proteins. Substantial homogenization of evolutionary rates in multidomain proteins was, indeed, observed in both animals and plants, although highly significant differences between domain-specific rates remained. The contributions of the translation rate, as determined by the effect of the fusion of a pair of domains within a multidomain protein, and intrinsic, domain-specific structural-functional constraints appear to be comparable in magnitude. Conclusion Fusion of domains in a multidomain protein results in substantial homogenization of the domain-specific evolutionary rates but significant differences between domain-specific evolution rates remain. Thus, the rate of translation and intrinsic structural-functional constraints both exert sizable and comparable effects on sequence evolution. Reviewers This article was reviewed by Sergei Maslov, Dennis Vitkup, Claus Wilke (nominated by Orly Alter), and Allan Drummond (nominated by Joel Bader). For the full reviews, please go to the Reviewers' Reports section.
Collapse
Affiliation(s)
- Maxim Y Wolf
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA.
| | | | | |
Collapse
|
105
|
An integrative approach to characterize disease-specific pathways and their coordination: a case study in cancer. BMC Genomics 2008; 9 Suppl 1:S12. [PMID: 18366601 PMCID: PMC2386054 DOI: 10.1186/1471-2164-9-s1-s12] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Background The most common application of microarray technology in disease research is to identify genes differentially expressed in disease versus normal tissues. However, it is known that, in complex diseases, phenotypes are determined not only by genes, but also by the underlying structure of genetic networks. Often, it is the interaction of many genes that causes phenotypic variations. Results In this work, using cancer as an example, we develop graph-based methods to integrate multiple microarray datasets to discover disease-related co-expression network modules. We propose an unsupervised method that take into account both co-expression dynamics and network topological information to simultaneously infer network modules and phenotype conditions in which they are activated or de-activated. Using our method, we have discovered network modules specific to cancer or subtypes of cancers. Many of these modules are consistent with or supported by their functional annotations or their previously known involvement in cancer. In particular, we identified a module that is predominately activated in breast cancer and is involved in tumor suppression. While individual components of this module have been suggested to be associated with tumor suppression, their coordinated function has never been elucidated. Here by adopting a network perspective, we have identified their interrelationships and, particularly, a hub gene PDGFRL that may play an important role in this tumor suppressor network. Conclusion Using a network-based approach, our method provides new insights into the complex cellular mechanisms that characterize cancer and cancer subtypes. By incorporating co-expression dynamics information, our approach can not only extract more functionally homogeneous modules than those based solely on network topology, but also reveal pathway coordination beyond co-expression.
Collapse
|
106
|
Yamashita R, Suzuki Y, Takeuchi N, Wakaguri H, Ueda T, Sugano S, Nakai K. Comprehensive detection of human terminal oligo-pyrimidine (TOP) genes and analysis of their characteristics. Nucleic Acids Res 2008; 36:3707-15. [PMID: 18480124 PMCID: PMC2441802 DOI: 10.1093/nar/gkn248] [Citation(s) in RCA: 92] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2007] [Revised: 03/25/2008] [Accepted: 04/17/2008] [Indexed: 12/03/2022] Open
Abstract
Although the knowledge accumulated on the transcriptional regulations of eukaryotes is significant, the knowledge on their translational regulations remains limited. Thus, we performed a comprehensive detection of terminal oligo-pyrimidine (TOP), which is one of the well-characterized cis-regulatory motifs for translational controls located immediately downstream of the transcriptional start sites of mRNAs. Utilizing our precise 5'-end information of the full-length cDNAs, we could screen 1645 candidate TOP genes by position specific matrix search. Among them, not only 75 out of 78 ribosomal protein genes but also eight previously identified non-ribosomal-protein TOP genes were included. We further experimentally validated the translational activities of 83 TOP candidate genes. Clear translational regulations exerted on the stimulation of 12-O-tetradecanoyl-1-phorbol-13-acetate for at least 41 of them was observed, indicating that there should be a few hundreds of human genes which are subjected to regulation at translation levels via TOPs. Our result suggests that TOP genes code not only formerly characterized ribosomal proteins and translation-related proteins but also a wider variety of proteins, such as lysosome-related proteins and metabolism-related proteins, playing pivotal roles in gene expression controls in the majority of cellular mRNAs.
Collapse
Affiliation(s)
- Riu Yamashita
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Yutaka Suzuki
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Nono Takeuchi
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Hiroyuki Wakaguri
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Takuya Ueda
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Sumio Sugano
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| | - Kenta Nakai
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo 108-8639, Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa, Chiba 277-8562 and Institute for Bioinformatics Research and Development (BIRD), Japan Science and Technology Agency (JST), 4-5-3 Chiyoda-ku, Tokyo, Japan
| |
Collapse
|
107
|
Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs. BMC Genomics 2008; 9:152. [PMID: 18384671 PMCID: PMC2335122 DOI: 10.1186/1471-2164-9-152] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2007] [Accepted: 04/02/2008] [Indexed: 12/30/2022] Open
Abstract
Background Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial. Results We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation. Conclusion Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.
Collapse
|
108
|
Ala U, Piro RM, Grassi E, Damasco C, Silengo L, Oti M, Provero P, Di Cunto F. Prediction of human disease genes by human-mouse conserved coexpression analysis. PLoS Comput Biol 2008; 4:e1000043. [PMID: 18369433 PMCID: PMC2268251 DOI: 10.1371/journal.pcbi.1000043] [Citation(s) in RCA: 104] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2007] [Accepted: 02/20/2008] [Indexed: 01/14/2023] Open
Abstract
BACKGROUND Even in the post-genomic era, the identification of candidate genes within loci associated with human genetic diseases is a very demanding task, because the critical region may typically contain hundreds of positional candidates. Since genes implicated in similar phenotypes tend to share very similar expression profiles, high throughput gene expression data may represent a very important resource to identify the best candidates for sequencing. However, so far, gene coexpression has not been used very successfully to prioritize positional candidates. METHODOLOGY/PRINCIPAL FINDINGS We show that it is possible to reliably identify disease-relevant relationships among genes from massive microarray datasets by concentrating only on genes sharing similar expression profiles in both human and mouse. Moreover, we show systematically that the integration of human-mouse conserved coexpression with a phenotype similarity map allows the efficient identification of disease genes in large genomic regions. Finally, using this approach on 850 OMIM loci characterized by an unknown molecular basis, we propose high-probability candidates for 81 genetic diseases. CONCLUSION Our results demonstrate that conserved coexpression, even at the human-mouse phylogenetic distance, represents a very strong criterion to predict disease-relevant relationships among human genes.
Collapse
Affiliation(s)
- Ugo Ala
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
| | - Rosario Michael Piro
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
| | - Elena Grassi
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
| | - Christian Damasco
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
| | - Lorenzo Silengo
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
| | - Martin Oti
- Department of Human Genetics and Centre for Molecular and Biomolecular Informatics, University Medical Centre Nijmegen, Nijmegen, The Netherlands
| | - Paolo Provero
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
- * E-mail: (PP); (FDC)
| | - Ferdinando Di Cunto
- Molecular Biotechnology Center, Department of Genetics, Biology and Biochemistry, University of Turin, Turin, Italy
- * E-mail: (PP); (FDC)
| |
Collapse
|
109
|
Kiesel J, Miller C, Abu-Amer Y, Aurora R. Systems level analysis of osteoclastogenesis reveals intrinsic and extrinsic regulatory interactions. Dev Dyn 2007; 236:2181-97. [PMID: 17584858 DOI: 10.1002/dvdy.21206] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Osteoclasts are bone-resorbing cells derived from the myeloid lineage that play a central role in bone remodeling and inflammatory bone erosion diseases. The receptor activator of NF-kappaB ligand (RANKL) produced by osteoblasts and activated immune cells initiates the development of osteoclasts in the bone marrow. Using time series gene expression data, the intrinsic processes and the extrinsic factors that control osteoclastogenesis were identified. The gene expression profiles display distinct commitment and differentiation phases. Analysis of the time course revealed several mechanistic details, including the complex role of cholesterol in osteoclast development. Epistatic interactions and the coordination between cellular processes that regulate development were inferred from the coexpression network. The coexpression network indicated that osteoclasts induce angiogenesis and recruit T-cells to the site of osteoclastogenesis early in the commitment phase. The resulting model provides an essential framework for a better understanding of the epigenetic program of osteoclastogenesis.
Collapse
Affiliation(s)
- Jennifer Kiesel
- Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, Saint Louis, Missouri, USA
| | | | | | | |
Collapse
|
110
|
Aurora R, Hihara Y, Singh AK, Pakrasi HB. A network of genes regulated by light in cyanobacteria. OMICS-A JOURNAL OF INTEGRATIVE BIOLOGY 2007; 11:166-85. [PMID: 17594236 DOI: 10.1089/omi.2007.4323] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Oxygenic photosynthetic organisms require light for their growth and development. However, exposure to high light is detrimental to them. Using time series microarray data from a model cyanobacterium, Synechocystis 6803 transferred from low to high light, we generated a gene co-expression network. The network has twelve sub-networks connected hierarchically, each consisting of an interconnected hub-and-spoke architecture. Within each sub-network, edges formed between genes that recapitulate known pathways. Analysis of the expression profiles shows that the cells undergo a phase transition 6-hours post-shift to high light, characterized by core sub-network. The core sub-network is enriched in proteins that (putatively) bind Fe-S clusters and proteins that mediate iron and sulfate homeostasis. At the center of this core is a sulfate permease, suggesting sulfate is rate limiting for cells grown in high light. To validate this novel finding, we demonstrate the limited ability of cell growth in sulfate-depleted medium in high light. This study highlights how understanding the organization of the networks can provide insights into the coordination of physiologic responses.
Collapse
Affiliation(s)
- Rajeev Aurora
- Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, St. Louis, Missouri 63104, USA.
| | | | | | | |
Collapse
|
111
|
Ucar D, Neuhaus I, Ross-MacDonald P, Tilford C, Parthasarathy S, Siemers N, Ji RR. Construction of a reference gene association network from multiple profiling data: application to data analysis. Bioinformatics 2007; 23:2716-24. [PMID: 17846039 DOI: 10.1093/bioinformatics/btm423] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
MOTIVATION Gene expression profiling is an important tool for gaining insight into biology. Novel strategies are required to analyze the growing archives of microarray data and extract useful information from them. One area of interest is in the construction of gene association networks from collections of profiling data. Various approaches have been proposed to construct gene networks using profiling data, and these networks have been used in functional inference as well as in data visualization. Here, we investigated a non-parametric approach to translate profiling data into a gene network. We explored the characteristics and utility of the resulting network and investigated the use of network information in analysis of variance models and hypothesis testing. RESULTS Our work is composed of two parts: gene network construction and partitioning and hypothesis testing using sub-networks as groups. In the first part, multiple independently collected microarray datasets from the Gene Expression Omnibus data repository were analyzed to identify probe pairs that are positively co-regulated across the samples. A co-expression network was constructed based on a reciprocal ranking criteria and a false discovery rate analysis. We named this network Reference Gene Association (RGA) network. Then, the network was partitioned into densely connected sub-networks of probes using a multilevel graph partitioning algorithm. In the second part, we proposed a new, MANOVA-based approach that can take individual probe expression values as input and perform hypothesis testing at the sub-network level. We applied this MANOVA methodology to two published studies and our analysis indicated that the methodology is both effective and sensitive for identifying transcriptional sub-networks or pathways that are perturbed across treatments.
Collapse
Affiliation(s)
- Duygu Ucar
- Computer Science and Engineering, The Ohio State University, Columbus, OH 43210, USA
| | | | | | | | | | | | | |
Collapse
|
112
|
Yeang CH, Darot JFJ, Noller HF, Haussler D. Detecting the coevolution of biosequences--an example of RNA interaction prediction. Mol Biol Evol 2007; 24:2119-31. [PMID: 17636042 DOI: 10.1093/molbev/msm142] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
A probabilistic graphical model is proposed in order to detect the coevolution between different sites in biological sequences. The model extends the continuous-time Markov process of sequence substitution for single nucleic or amino acids and imposes general constraints regarding simultaneous changes on the substitution rate matrix. Given a multiple sequence alignment for each molecule of interest and a phylogenetic tree, the model can predict potential interactions within or between nucleic acids and proteins. Initial validation of the model is carried out using tRNA and 16S rRNA sequence data. The model accurately identifies the secondary interactions of tRNA as well as several known tertiary interactions. In addition, results on 16S rRNA data indicate this general and simple coevolutionary model outperforms several other parametric and nonparametric methods in predicting secondary interactions. Furthermore, the majority of the putative predictions exhibit either direct contact or proximity of the nucleotide pairs in the 3-dimensional structure of the Thermus thermophilus ribosomal small subunit. The results on RNA data suggest a general model of coevolution might be applied to other types of interactions between protein, DNA, and RNA molecules.
Collapse
Affiliation(s)
- Chen-Hsiang Yeang
- Simons Center for Systems Biology, Institute for Advanced Study, Princeton, New Jersey, USA.
| | | | | | | |
Collapse
|
113
|
Carmel L, Rogozin IB, Wolf YI, Koonin EV. Evolutionarily conserved genes preferentially accumulate introns. Genes Dev 2007; 17:1045-50. [PMID: 17495009 PMCID: PMC1899115 DOI: 10.1101/gr.5978207] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2006] [Accepted: 03/28/2007] [Indexed: 11/24/2022]
Abstract
Introns that interrupt eukaryotic protein-coding sequences are generally thought to be nonfunctional. However, for reasons still poorly understood, positions of many introns are highly conserved in evolution. Previous reconstructions of intron gain and loss events during eukaryotic evolution used a variety of simplified evolutionary models that yielded contradicting conclusions and are not suited to reveal some of the key underlying processes. We combine a comprehensive probabilistic model and an extended data set, including 391 conserved genes from 19 eukaryotes, to uncover previously unnoticed aspects of intron evolution--in particular, to assign intron gain and loss rates to individual genes. The rates of intron gain and loss in a gene show moderate positive correlation. A gene's intron gain rate shows a highly significant negative correlation with the coding-sequence evolution rate; intron loss rate also significantly, but positively, correlates with the sequence evolution rate. Correlations of the opposite signs, albeit less significant ones, are observed between intron gain and loss rates and gene expression level. It is proposed that intron evolution includes a neutral component, which is manifest in the positive correlation between the gain and loss rates and a selection-driven component as reflected in the links between intron gain and loss and sequence evolution. The increased intron gain and decreased intron loss in evolutionarily conserved genes indicate that intron insertion often might be adaptive, whereas some of the intron losses might be deleterious. This apparent functional importance of introns is likely to be due, at least in part, to their multiple effects on gene expression.
Collapse
Affiliation(s)
- Liran Carmel
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894, USA
| | - Igor B. Rogozin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894, USA
| | - Yuri I. Wolf
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894, USA
| | - Eugene V. Koonin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894, USA
| |
Collapse
|
114
|
Gitelman I. Evolution of the vertebrate twist family and synfunctionalization: a mechanism for differential gene loss through merging of expression domains. Mol Biol Evol 2007; 24:1912-25. [PMID: 17567594 DOI: 10.1093/molbev/msm120] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
Twist genes are essential for embryonic development and are conserved from jellyfish to human. To study the vertebrate twist family and its evolution, the entire complement of twist genes was obtained for 9 representative species. Phylogenetic analysis showed that a single protochordate twist gene was duplicated at least twice before the teleost-tetrapod split to give rise to 3 ancestral genes, which were further duplicated or deleted, resulting in fluctuating number of twist paralogs in different vertebrate lineages. To find whether changes in gene copy number were associated with changes in gene function, embryonic expression patterns of twist orthologs were evaluated against the number of twist paralogs in different species. The results showed evidence for both neo- and subfunctionalization, and, in addition, for loss of an ancestral regulatory gene. For example, in Xenopus, twist2 was lost, but the twist1 paralog acquired, and therefore preserved, twist2 function. A general model is proposed to explain the data. In this process, termed synfunctionalization, one paralog acquires the expression domain(s) of another. The merging may lead to function shuffle. Alternatively, it may leave one paralog redundant and thus subject to deletion--while its function is retained by the surviving paralog(s). Synfunctionalization is a mechanism that, together with neo- and subfunctionalization, may work to establish equilibrium in the number of genes that regulate developmental processes; it may regulate the complexity of regulatory regions as well as gene copy number and therefore may play a role in evolution of gene function and the structure of genome.
Collapse
Affiliation(s)
- Inna Gitelman
- Department of Virology and Developmental Genetics, Faculty of Health Sciences, Ben Gurion University of the Negev, Beer Sheva, 84105, Israel.
| |
Collapse
|
115
|
Chen J, Blackwell TW, Fermin D, Menon R, Chen Y, Gao J, Lee AW, States DJ. Evolutionary-conserved gene expression response profiles across mammalian tissues. OMICS : A JOURNAL OF INTEGRATIVE BIOLOGY 2007; 11:96-115. [PMID: 17411398 DOI: 10.1089/omi.2006.0007] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Gene expression responses are complex and frequently involve the actions of many genes to effect coordinated patterns. We hypothesized these coordinated responses are evolutionarily conserved and used a comparison of human and mouse gene expression profiles to identify the most prominent conserved features across a set of normal mammalian tissues. Based on data from multiple studies across multiple tissues in human and mouse, 13 gene expression modes across multiple tissues were identified in each of these species using principal component analysis. Strikingly, 1-to-1 pairing of human and mouse modes was observed in 12 out of 13 modes obtained from the two species independently. These paired modes define evolutionarily conserved gene expression response modes (CGEMs). Notably, in this study we were able to extract biological responses that are not overwhelmed by laboratory-to-laboratory or species-to-species variation. Of the variation in our gene expression dataset, 84% can be explained using these CGEMs. Functional annotation was performed using Gene Ontology, pathway, and transcription factor binding site over representation. Our conclusion is that we found an unbiased way of obtaining conserved gene response modes that accounts for a considerable portion of gene expression variation in a given dataset, as well as validates the conservation of major gene expression response modes across the mammals.
Collapse
Affiliation(s)
- Ji Chen
- Bioinformatics Program, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | | | | | | | | | | | | |
Collapse
|
116
|
Oldham MC, Horvath S, Geschwind DH. Conservation and evolution of gene coexpression networks in human and chimpanzee brains. Proc Natl Acad Sci U S A 2006; 103:17973-8. [PMID: 17101986 PMCID: PMC1693857 DOI: 10.1073/pnas.0605938103] [Citation(s) in RCA: 416] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2006] [Indexed: 12/28/2022] Open
Abstract
Comparisons of gene expression between human and non-human primate brains have identified hundreds of differentially expressed genes, yet translating these lists into key functional distinctions between species has proved difficult. Here we provide a more integrated view of human brain evolution by examining the large-scale organization of gene coexpression networks in human and chimpanzee brains. We identify modules of coexpressed genes that correspond to discrete brain regions and quantify their conservation between the species. Module conservation in cerebral cortex is significantly weaker than module conservation in subcortical brain regions, revealing a striking gradient that parallels known evolutionary hierarchies. We introduce a method for identifying species-specific network connections and demonstrate how differential network connectivity can be used to identify key drivers of evolutionary change. By integrating our results with comparative genomic sequence data and estimates of protein sequence divergence rates, we confirm a number of network predictions and validate these findings. Our results provide insights into the molecular bases of primate brain organization and demonstrate the general utility of weighted gene coexpression network analysis.
Collapse
Affiliation(s)
- Michael C. Oldham
- *Interdepartmental Program for Neuroscience
- Neurology, Program in Neurogenetics, and
- Semel Institute, David Geffen School of Medicine, 710 Westwood Plaza, Los Angeles, CA 90095-6814
| | | | - Daniel H. Geschwind
- Human Genetics, and
- Neurology, Program in Neurogenetics, and
- Semel Institute, David Geffen School of Medicine, 710 Westwood Plaza, Los Angeles, CA 90095-6814
| |
Collapse
|
117
|
Donaldson IJ, Göttgens B. Evolution of candidate transcriptional regulatory motifs since the human-chimpanzee divergence. Genome Biol 2006; 7:R52. [PMID: 16808854 PMCID: PMC1779530 DOI: 10.1186/gb-2006-7-6-r52] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2006] [Revised: 06/01/2006] [Accepted: 06/09/2006] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Despite the recent completion of the chimpanzee genome project, few functionally significant sequence differences between humans and chimpanzees have thus far been identified. Alteration in transcriptional regulatory mechanisms represents an important platform for evolutionary change, suggesting that a significant proportion of functional human-chimpanzee sequence differences may affect regulatory elements. RESULTS To explore this hypothesis, we performed genome-wide identification of conserved candidate transcription-factor binding sites that have evolved since the divergence of humans and chimpanzees. Analysis of candidate transcription-factor binding sites conserved between mouse and chimpanzee yet absent in human indicated that loss of candidate transcription-factor binding sites in the human lineage was not random but instead correlated with the biologic functions of associated genes. CONCLUSION Our data support the notion that changes in transcriptional regulation have contributed to the recent evolution of humans. Moreover, genes associated with mutated candidate transcription-factor binding sites highlight potential pathways underlying human-chimpanzee divergence.
Collapse
Affiliation(s)
- Ian J Donaldson
- Department of Haematology, Cambridge Institute for Medical Research, University of Cambridge, Cambridge, CB2 2XY, UK
| | - Berthold Göttgens
- Department of Haematology, Cambridge Institute for Medical Research, University of Cambridge, Cambridge, CB2 2XY, UK
| |
Collapse
|
118
|
Tsaparas P, Mariño-Ramírez L, Bodenreider O, Koonin EV, Jordan IK. Global similarity and local divergence in human and mouse gene co-expression networks. BMC Evol Biol 2006; 6:70. [PMID: 16968540 PMCID: PMC1601971 DOI: 10.1186/1471-2148-6-70] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2006] [Accepted: 09/12/2006] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND A genome-wide comparative analysis of human and mouse gene expression patterns was performed in order to evaluate the evolutionary divergence of mammalian gene expression. Tissue-specific expression profiles were analyzed for 9,105 human-mouse orthologous gene pairs across 28 tissues. Expression profiles were resolved into species-specific coexpression networks, and the topological properties of the networks were compared between species. RESULTS At the global level, the topological properties of the human and mouse gene coexpression networks are, essentially, identical. For instance, both networks have topologies with small-world and scale-free properties as well as closely similar average node degrees, clustering coefficients, and path lengths. However, the human and mouse coexpression networks are highly divergent at the local level: only a small fraction (<10%) of coexpressed gene pair relationships are conserved between the two species. A series of controls for experimental and biological variance show that most of this divergence does not result from experimental noise. We further show that, while the expression divergence between species is genuinely rapid, expression does not evolve free from selective (functional) constraint. Indeed, the coexpression networks analyzed here are demonstrably functionally coherent as indicated by the functional similarity of coexpressed gene pairs, and this pattern is most pronounced in the conserved human-mouse intersection network. Numerous dense network clusters show evidence of dedicated functions, such as spermatogenesis and immune response, that are clearly consistent with the coherence of the expression patterns of their constituent gene members. CONCLUSION The dissonance between global versus local network divergence suggests that the interspecies similarity of the global network properties is of limited biological significance, at best, and that the biologically relevant aspects of the architectures of gene coexpression are specific and particular, rather than universal. Nevertheless, there is substantial evolutionary conservation of the local network structure which is compatible with the notion that gene coexpression networks are subject to purifying selection.
Collapse
Affiliation(s)
- Panayiotis Tsaparas
- Basic Research Unit, Helsinki Institute for Information Technology, University of Helsinki, Helsinki, Finland
| | - Leonardo Mariño-Ramírez
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA
| | - Olivier Bodenreider
- National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA
| | - Eugene V Koonin
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA
| | - I King Jordan
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD 20894, USA
- Present address: School of Biology, Georgia Institute of Technology, Atlanta, GA 30332, USA
| |
Collapse
|
119
|
Koonin EV, Wolf YI. Evolutionary systems biology: links between gene evolution and function. Curr Opin Biotechnol 2006; 17:481-7. [PMID: 16962765 DOI: 10.1016/j.copbio.2006.08.003] [Citation(s) in RCA: 76] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2006] [Revised: 08/10/2006] [Accepted: 08/30/2006] [Indexed: 11/25/2022]
Abstract
The recent accumulation of genome-wide data on various facets of gene expression, function and evolution stimulated the emergence of a new field, evolutionary systems biology. Many significant correlations were detected between variables that characterize the functioning of a gene, such as expression level, knockout effect, connectivity of genetic and protein-protein interaction networks, and variables that describe gene evolution, such as sequence evolution rate and propensity for gene loss. The first attempts on multidimensional analysis of genomic data yielded composite variables that describe the 'status' of a gene in the genomic community. However, it remains uncertain whether different functional variables affect gene evolution synergistically or there is a single, dominant factor. The number of translation events, linked to selection for translational robustness, was proposed as a candidate for such a major determinant of protein evolution. These developments show that, although the methodological basis of evolutionary systems biology is not yet fully solidified, this area of research is already starting to yield fundamental biological insights.
Collapse
Affiliation(s)
- Eugene V Koonin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA.
| | | |
Collapse
|
120
|
Abstract
It has been suggested that evolutionary changes in gene expression account for most phenotypic differences between species, in particular between humans and apes. What general rules can be described governing expression evolution? We find that a neutral model where negative selection and divergence time are the major factors is a useful null hypothesis for both transcriptome and genome evolution. Two tissues that stand out with regard to gene expression are the testes, where positive selection has exerted a substantial influence in both humans and chimpanzees, and the brain, where gene expression has changed less than in other organs but acceleration might have occurred in human ancestors.
Collapse
Affiliation(s)
- Philipp Khaitovich
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, D-04103 Leipzig, Germany
| | | | | | | |
Collapse
|
121
|
Hughes AL, Friedman R. Sharing of transcription factors after gene duplication in the yeast Saccharomyces cerevisiae. Genetica 2006; 129:301-8. [PMID: 16897462 DOI: 10.1007/s10709-006-0011-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2006] [Accepted: 04/20/2006] [Indexed: 10/24/2022]
Abstract
In a set of 190 duplicate gene pairs in yeast Saccharomyces cerevisiae, the sharing of transcription factors tended to decrease with increased divergence in coding sequence, at both synonymous and nonsynonymous sites. Our results showed a significantly higher sharing of transcription factors by duplicated gene pairs falling within duplicated genomic blocks than in other duplicated gene pairs; and genes in duplicated blocks also showed significantly greater conservation at the coding sequence level. In spite of the overall trends, there were certain gene pairs, both in duplicated blocks and in other genomic regions, which were highly divergent in coding sequence and yet had identical patterns of transcription factor binding. These results suggest that functional differentiation of genes after duplication is a multi-dimensional process, with different duplicate pairs differentiating in different ways.
Collapse
Affiliation(s)
- Austin L Hughes
- Department of Biological Sciences, University of South Carolina, Coker Life Sciences Bldg. 700 Sumter St., Columbia, SC 29208, USA.
| | | |
Collapse
|
122
|
Abstract
Recent genome analyses revealed intriguing correlations between variables characterizing the functioning of a gene, such as expression level (EL), connectivity of genetic and protein-protein interaction networks, and knockout effect, and variables describing gene evolution, such as sequence evolution rate (ER) and propensity for gene loss. Typically, variables within each of these classes are positively correlated, e.g. products of highly expressed genes also have a propensity to be involved in many protein-protein interactions, whereas variables between classes are negatively correlated, e.g. highly expressed genes, on average, evolve slower than weakly expressed genes. Here, we describe principal component (PC) analysis of seven genome-related variables and propose biological interpretations for the first three PCs. The first PC reflects a gene's 'importance', or the 'status' of a gene in the genomic community, with positive contributions from knockout lethality, EL, number of protein-protein interaction partners and the number of paralogues, and negative contributions from sequence ER and gene loss propensity. The next two PCs define a plane that seems to reflect the functional and evolutionary plasticity of a gene. Specifically, PC2 can be interpreted as a gene's 'adaptability' whereby genes with high adaptability readily duplicate, have many genetic interaction partners and tend to be non-essential. PC3 also might reflect the role of a gene in organismal adaptation albeit with a negative rather than a positive contribution of genetic interactions; we provisionally designate this PC 'reactivity'. The interpretation of PC2 and PC3 as measures of a gene's plasticity is compatible with the observation that genes with high values of these PCs tend to be expressed in a condition- or tissue-specific manner. Functional classes of genes substantially vary in status, adaptability and reactivity, with the highest status characteristic of the translation system and cytoskeletal proteins, highest adaptability seen in cellular processes and signalling genes, and top reactivity characteristic of metabolic enzymes.
Collapse
|
123
|
Abstract
Why do proteins evolve at different rates? Advances in systems biology and genomics have facilitated a move from studying individual proteins to characterizing global cellular factors. Systematic surveys indicate that protein evolution is not determined exclusively by selection on protein structure and function, but is also affected by the genomic position of the encoding genes, their expression patterns, their position in biological networks and possibly their robustness to mistranslation. Recent work has allowed insights into the relative importance of these factors. We discuss the status of a much-needed coherent view that integrates studies on protein evolution with biochemistry and functional and structural genomics.
Collapse
Affiliation(s)
- Csaba Pál
- European Molecular Biology Laboratory, Meyerhofstrasse 1, D-69012 Heidelberg, Germany
| | | | | |
Collapse
|
124
|
Feldmesser E, Olender T, Khen M, Yanai I, Ophir R, Lancet D. Widespread ectopic expression of olfactory receptor genes. BMC Genomics 2006; 7:121. [PMID: 16716209 PMCID: PMC1508154 DOI: 10.1186/1471-2164-7-121] [Citation(s) in RCA: 182] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2006] [Accepted: 05/22/2006] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Olfactory receptors (ORs) are the largest gene family in the human genome. Although they are expected to be expressed specifically in olfactory tissues, some ectopic expression has been reported, with special emphasis on sperm and testis. The present study systematically explores the expression patterns of OR genes in a large number of tissues and assesses the potential functional implication of such ectopic expression. RESULTS We analyzed the expression of hundreds of human and mouse OR transcripts, via EST and microarray data, in several dozens of human and mouse tissues. Different tissues had specific, relatively small OR gene subsets which had particularly high expression levels. In testis, average expression was not particularly high, and very few highly expressed genes were found, none corresponding to ORs previously implicated in sperm chemotaxis. Higher expression levels were more common for genes with a non-OR genomic neighbor. Importantly, no correlation in expression levels was detected for human-mouse orthologous pairs. Also, no significant difference in expression levels was seen between intact and pseudogenized ORs, except for the pseudogenes of subfamily 7E which has undergone a human-specific expansion. CONCLUSION The OR superfamily as a whole, show widespread, locus-dependent and heterogeneous expression, in agreement with a neutral or near neutral evolutionary model for transcription control. These results cannot reject the possibility that small OR subsets might play functional roles in different tissues, however considerable care should be exerted when offering a functional interpretation for ectopic OR expression based only on transcription information.
Collapse
Affiliation(s)
- Ester Feldmesser
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Tsviya Olender
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Miriam Khen
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Itai Yanai
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel
- Present address: Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA
| | - Ron Ophir
- Department of Biological Services, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Doron Lancet
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
125
|
Millan MJ. Multi-target strategies for the improved treatment of depressive states: Conceptual foundations and neuronal substrates, drug discovery and therapeutic application. Pharmacol Ther 2006; 110:135-370. [PMID: 16522330 DOI: 10.1016/j.pharmthera.2005.11.006] [Citation(s) in RCA: 389] [Impact Index Per Article: 21.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2005] [Accepted: 11/28/2005] [Indexed: 12/20/2022]
Abstract
Major depression is a debilitating and recurrent disorder with a substantial lifetime risk and a high social cost. Depressed patients generally display co-morbid symptoms, and depression frequently accompanies other serious disorders. Currently available drugs display limited efficacy and a pronounced delay to onset of action, and all provoke distressing side effects. Cloning of the human genome has fuelled expectations that symptomatic treatment may soon become more rapid and effective, and that depressive states may ultimately be "prevented" or "cured". In pursuing these objectives, in particular for genome-derived, non-monoaminergic targets, "specificity" of drug actions is often emphasized. That is, priority is afforded to agents that interact exclusively with a single site hypothesized as critically involved in the pathogenesis and/or control of depression. Certain highly selective drugs may prove effective, and they remain indispensable in the experimental (and clinical) evaluation of the significance of novel mechanisms. However, by analogy to other multifactorial disorders, "multi-target" agents may be better adapted to the improved treatment of depressive states. Support for this contention is garnered from a broad palette of observations, ranging from mechanisms of action of adjunctive drug combinations and electroconvulsive therapy to "network theory" analysis of the etiology and management of depressive states. The review also outlines opportunities to be exploited, and challenges to be addressed, in the discovery and characterization of drugs recognizing multiple targets. Finally, a diversity of multi-target strategies is proposed for the more efficacious and rapid control of core and co-morbid symptoms of depression, together with improved tolerance relative to currently available agents.
Collapse
Affiliation(s)
- Mark J Millan
- Institut de Recherches Servier, Centre de Recherches de Croissy, Psychopharmacology Department, 125, Chemin de Ronde, 78290-Croissy/Seine, France.
| |
Collapse
|
126
|
Sällström B, Arnaout RA, Davids W, Bjelkmar P, Andersson SGE. Protein evolutionary rates correlate with expression independently of synonymous substitutions in Helicobacter pylori. J Mol Evol 2006; 62:600-14. [PMID: 16586017 DOI: 10.1007/s00239-005-0104-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2005] [Accepted: 12/20/2005] [Indexed: 11/29/2022]
Abstract
In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites.
Collapse
Affiliation(s)
- Björn Sällström
- Program of Molecular Evolution, Department of Evolution, Genomics and Systematics, Evolutionary Biology Center, Uppsala University, 752 36 Uppsala, Sweden
| | | | | | | | | |
Collapse
|
127
|
Oldham MC, Geschwind DH. Comparative genomics: Grasping human transcriptome evolution: what does it all mean? Heredity (Edinb) 2006; 96:339-40. [PMID: 16552432 DOI: 10.1038/sj.hdy.6800807] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
|
128
|
Liu S, Zhang C, Zhou Y. Uneven size distribution of mammalian genes in the number of tissues expressed and in the number of co-expressed genes. Hum Mol Genet 2006; 15:1313-8. [PMID: 16537573 DOI: 10.1093/hmg/ddl051] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Tissue specificity, the traditional predictor of gene function, has recently been used to interpret the selective pressure associated with gene architecture. In this work, we examine gene structures and their relation to the number of tissues expressed and to the number of co-expressed genes, using a recent atlas of microarray-based mouse gene expression in 55 normal tissues. We define tissue specificity and expression-pattern specificity according to the number of tissues expressed and the number of co-expressed genes, respectively. We find that, consistent with previous findings, tissue non-specific (housekeeping) genes are short in all gene regions (coding regions, intron, 5' and 3' untranslated regions). However, in contrast to previous suggestion that tissue-specific genes are long, the genes that are the most tissue-specific (expressed only in one tissue) are also short. We further show that both expression-pattern-specific and non-specific genes are long in coding and non-coding regions. The origins for short tissue-specific genes and long expression-pattern-specific genes are not clear. Genes with highly non-specific expression patterns (i.e. genes with a large number of co-expressed genes) are composed of genes that spread all tissues but are overwhelmingly enriched in the central nervous system (e.g. brain). Thus, the large sizes of these genes are possibly related to the functional complexity and/or accelerated evolutions of the central nervous system.
Collapse
Affiliation(s)
- Song Liu
- Department of Physiology and Biophysics, Howard Hughes Medical Institute Center for Single Molecule Biophysics, State University of New York at Buffalo, 124 Sherman Hall, Buffalo, NY 14214, USA
| | | | | |
Collapse
|
129
|
Liao BY, Zhang J. Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolution. Mol Biol Evol 2006; 23:1119-28. [PMID: 16520335 DOI: 10.1093/molbev/msj119] [Citation(s) in RCA: 112] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Evolutionary rates provide important information about the pattern and mechanism of evolution. Although the rate of gene sequence evolution has been well studied, the rate of gene expression evolution is poorly understood. In particular, it is unclear whether the gene expression level and tissue specificity influence the divergence of expression profiles between orthologous genes. Here we address this question using a microarray data set comprising the expression signals of 10,607 pairs of orthologous human and mouse genes from over 60 tissues per species. We show that the level of gene expression and the degree of tissue specificity are generally conserved between the human and mouse orthologs. The rate of gene expression profile change during evolution is negatively correlated with the level of gene expression, measured by either the average or the highest level among all tissues examined. This is analogous to the observation that the rate of gene (or protein) sequence evolution is negatively correlated with the gene expression level. The impacts of the degree of tissue specificity on the evolutionary rate of gene sequence and that of expression profile, however, are opposite. Highly tissue-specific genes tend to evolve rapidly at the gene sequence level but slowly at the expression profile level. Thus, different forces and selective constraints must underlie the evolution of gene sequence and that of gene expression.
Collapse
Affiliation(s)
- Ben-Yang Liao
- Department of Ecology and Evolutionary Biology, University of Michigan, USA
| | | |
Collapse
|
130
|
Hao W, Golding GB. Asymmetrical Evolution of Cytochrome bd Subunits. J Mol Evol 2006; 62:132-42. [PMID: 16474982 DOI: 10.1007/s00239-005-0005-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2005] [Accepted: 09/18/2005] [Indexed: 11/27/2022]
Abstract
Functionally linked genes generally evolve at similar rates and the knowledge of this particular feature of genomic evolution has been used as the basis for the phylogenetic profiling method. We illustrate here an exception to this rule in the evolution of the cytochrome bd complex. This is a two-component oxidase complex, with the subunits I and II known to be widely present in bacteria. The subunits within the cytochrome bd complex are under the same evolutionary pressure and most likely behave in the same evolutionary manner. However, the sequence similarity of genes encoding subunit II varies considerably across species. Genes encoding subunit II evolve 1.2 times faster on most of the branches of their phylogeny than subunit I genes. Furthermore, the genes encoding subunit II in Oceanobacillus iheyensis, Bacillus halodurans, and Staphylococcus species do not have detectable homologues within E. coli due to their large divergence. Together, the two subunits of cytochrome bd reveal an interesting example of an asymmetric pattern of evolutionary change.
Collapse
Affiliation(s)
- Weilong Hao
- Department of Biology, McMaster University, Hamilton, Ontario, L8S 4K1, Canada
| | | |
Collapse
|
131
|
Chung WY, Albert R, Albert I, Nekrutenko A, Makova KD. Rapid and asymmetric divergence of duplicate genes in the human gene coexpression network. BMC Bioinformatics 2006; 7:46. [PMID: 16441884 PMCID: PMC1403810 DOI: 10.1186/1471-2105-7-46] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2005] [Accepted: 01/27/2006] [Indexed: 11/29/2022] Open
Abstract
Background While gene duplication is known to be one of the most common mechanisms of genome evolution, the fates of genes after duplication are still being debated. In particular, it is presently unknown whether most duplicate genes preserve (or subdivide) the functions of the parental gene or acquire new functions. One aspect of gene function, that is the expression profile in gene coexpression network, has been largely unexplored for duplicate genes. Results Here we build a human gene coexpression network using human tissue-specific microarray data and investigate the divergence of duplicate genes in it. The topology of this network is scale-free. Interestingly, our analysis indicates that duplicate genes rapidly lose shared coexpressed partners: after approximately 50 million years since duplication, the two duplicate genes in a pair have only slightly higher number of shared partners as compared with two random singletons. We also show that duplicate gene pairs quickly acquire new coexpressed partners: the average number of partners for a duplicate gene pair is significantly greater than that for a singleton (the latter number can be used as a proxy of the number of partners for a parental singleton gene before duplication). The divergence in gene expression between two duplicates in a pair occurs asymmetrically: one gene usually has more partners than the other one. The network is resilient to both random and degree-based in silico removal of either singletons or duplicate genes. In contrast, the network is especially vulnerable to the removal of highly connected genes when duplicate genes and singletons are considered together. Conclusion Duplicate genes rapidly diverge in their expression profiles in the network and play similar role in maintaining the network robustness as compared with singletons. Contact: kdm16@psu.edu Supplementary information: Please see additional files.
Collapse
Affiliation(s)
- Wen-Yu Chung
- Department of Computer Science and Engineering, Penn State University, University Park, PA, 16802, USA
- Department of Center for Comparative Genomics and Bioinformatics, Penn State University, University Park, PA, 16802, USA
| | - Reka Albert
- Department of Physics, Penn State University, University Park, PA, 16802, USA
| | - Istvan Albert
- Department of Huck Institute for Life Sciences, Penn State University, University Park, PA, 16802, USA
| | - Anton Nekrutenko
- Department of Biochemistry and Molecular Biology, Penn State University, University Park, PA, 16802, USA
- Department of Huck Institute for Life Sciences, Penn State University, University Park, PA, 16802, USA
- Department of Center for Comparative Genomics and Bioinformatics, Penn State University, University Park, PA, 16802, USA
| | - Kateryna D Makova
- Department of Biology, Penn State University, University Park, PA, 16802, USA
- Department of Center for Comparative Genomics and Bioinformatics, Penn State University, University Park, PA, 16802, USA
| |
Collapse
|
132
|
Barrett T, Edgar R. Gene expression omnibus: microarray data storage, submission, retrieval, and analysis. Methods Enzymol 2006; 411:352-69. [PMID: 16939800 PMCID: PMC1619900 DOI: 10.1016/s0076-6879(06)11019-8] [Citation(s) in RCA: 365] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The Gene Expression Omnibus (GEO) repository at the National Center for Biotechnology Information archives and freely distributes high-throughput molecular abundance data, predominantly gene expression data generated by DNA microarray technology. The database has a flexible design that can handle diverse styles of both unprocessed and processed data in a Minimum Information About a Microarray Experiment-supportive infrastructure that promotes fully annotated submissions. GEO currently stores about a billion individual gene expression measurements, derived from over 100 organisms, submitted by over 1500 laboratories, addressing a wide range of biological phenomena. To maximize the utility of these data, several user-friendly web-based interfaces and applications have been implemented that enable effective exploration, query, and visualization of these data at the level of individual genes or entire studies. This chapter describes how data are stored, submission procedures, and mechanisms for data retrieval and query. GEO is publicly accessible at http://www.ncbi.nlm.nih.gov/projects/geo/.
Collapse
Affiliation(s)
- Tanya Barrett
- National Center for Biotechnology Information, National Institutes of Health, Bethesda, MD, USA
| | | |
Collapse
|
133
|
|
134
|
Liao BY, Zhang J. Evolutionary Conservation of Expression Profiles Between Human and Mouse Orthologous Genes. Mol Biol Evol 2005; 23:530-40. [PMID: 16280543 DOI: 10.1093/molbev/msj054] [Citation(s) in RCA: 172] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Mouse models are often used to study human genes because it is believed that the expression and function are similar for the majority of orthologous genes between the two species. However, recent comparisons of microarray data from thousands of orthologous human and mouse genes suggested rapid evolution of gene expression profiles under minimal or no selective constraint. These findings appear to contradict non-array-based observations from many individual genes and imply the uselessness of mouse models for studying human genes. Because absolute levels of gene expression are not comparable between species when the data are generated by species-specific microarrays, use of relative mRNA abundance among tissues (RA) is preferred to that of absolute expression signals. We thus reanalyze human and mouse genome-wide gene expression data generated by oligonucleotide microarrays. We show that the mean correlation coefficient among expression profiles detected by different probe sets of the same gene is only 0.38 for humans and 0.28 for mice, indicating that current measures of expression divergence are flawed because the large estimation error (discrepancy in expression signal detected by different probe sets of the same gene) is mistakenly included in the between-species divergence. When this error is subtracted, 84% of human-mouse orthologous gene pairs show significantly lower expression divergence than that of random gene pairs. In contrast to a previous finding, but consistent with the common sense, expression profiles of orthologous tissues between species are more similar to each other than to those of nonorthologous tissues. Furthermore, the evolutionary rate of expression divergence and that of coding sequence divergence are found to be weakly, but significantly positively correlated, when RA and the Euclidean distance are used to measure expression-profile divergence. These results highlight the importance of proper consideration of various estimation errors in comparing the microarray data between species.
Collapse
Affiliation(s)
- Ben-Yang Liao
- Department of Ecology and Evolutionary Biology, University of Michigan, USA
| | | |
Collapse
|
135
|
Abstract
The notion of scale-freeness and its prevalence in both natural and artificial networks have recently attracted much attention. The concept of scale-freeness is enthusiastically applied to almost any conceivable network, usually with affirmative conclusions. Well-known scale-free examples include the internet, electric lines among power plants, the co-starring of movie actors, the co-authorship of researchers, food webs, and neural, protein-protein interactional, genetic, and metabolic networks. The purpose of this review is to clarify the relationship between scale-freeness and power-law distribution, and to assess critically the previous related works, especially on biological networks. In addition, I will focus on the close relationship between power-law distribution and lognormal distribution to show that power-law distribution is not a special characteristic of natural selection.
Collapse
Affiliation(s)
- Masanori Arita
- Department of Computational Biology, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwanoha 5-1-5 CB05, Kashiwa.
| |
Collapse
|
136
|
|
137
|
Khaitovich P, Hellmann I, Enard W, Nowick K, Leinweber M, Franz H, Weiss G, Lachmann M, Pääbo S. Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees. Science 2005; 309:1850-4. [PMID: 16141373 DOI: 10.1126/science.1108296] [Citation(s) in RCA: 419] [Impact Index Per Article: 22.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
The determination of the chimpanzee genome sequence provides a means to study both structural and functional aspects of the evolution of the human genome. Here we compare humans and chimpanzees with respect to differences in expression levels and protein-coding sequences for genes active in brain, heart, liver, kidney, and testis. We find that the patterns of differences in gene expression and gene sequences are markedly similar. In particular, there is a gradation of selective constraints among the tissues so that the brain shows the least differences between the species whereas liver shows the most. Furthermore, expression levels as well as amino acid sequences of genes active in more tissues have diverged less between the species than have genes active in fewer tissues. In general, these patterns are consistent with a model of neutral evolution with negative selection. However, for X-chromosomal genes expressed in testis, patterns suggestive of positive selection on sequence changes as well as expression changes are seen. Furthermore, although genes expressed in the brain have changed less than have genes expressed in other tissues, in agreement with previous work we find that genes active in brain have accumulated more changes on the human than on the chimpanzee lineage.
Collapse
MESH Headings
- Adult
- Aged
- Amino Acid Sequence
- Animals
- Base Sequence
- Child
- Chromosomes, Human, X/genetics
- Chromosomes, Mammalian/genetics
- Evolution, Molecular
- Female
- Gene Expression
- Gene Expression Profiling
- Gene Expression Regulation
- Genome
- Genome, Human
- Heart/physiology
- Humans
- Kidney/physiology
- Liver/physiology
- Male
- Middle Aged
- Models, Genetic
- Oligonucleotide Array Sequence Analysis
- Organ Specificity
- Pan troglodytes/genetics
- Prefrontal Cortex/physiology
- Promoter Regions, Genetic
- Proteins/genetics
- Selection, Genetic
- Sequence Analysis, DNA
- Species Specificity
- Testis/physiology
- Transcription, Genetic
- X Chromosome/genetics
Collapse
Affiliation(s)
- Philipp Khaitovich
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, D-04103 Leipzig, Germany
| | | | | | | | | | | | | | | | | |
Collapse
|
138
|
Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, Califano A. Reverse engineering of regulatory networks in human B cells. Nat Genet 2005; 37:382-90. [PMID: 15778709 DOI: 10.1038/ng1532] [Citation(s) in RCA: 917] [Impact Index Per Article: 48.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2004] [Accepted: 02/08/2005] [Indexed: 12/16/2022]
Abstract
Cellular phenotypes are determined by the differential activity of networks linking coregulated genes. Available methods for the reverse engineering of such networks from genome-wide expression profiles have been successful only in the analysis of lower eukaryotes with simple genomes. Using a new method called ARACNe (algorithm for the reconstruction of accurate cellular networks), we report the reconstruction of regulatory networks from expression profiles of human B cells. The results are suggestive a hierarchical, scale-free network, where a few highly interconnected genes (hubs) account for most of the interactions. Validation of the network against available data led to the identification of MYC as a major hub, which controls a network comprising known target genes as well as new ones, which were biochemically validated. The newly identified MYC targets include some major hubs. This approach can be generally useful for the analysis of normal and pathologic networks in mammalian cells.
Collapse
Affiliation(s)
- Katia Basso
- Institute for Cancer Genetics, 1300 St. Nicholas Avenue, Room 912, New York, New York 10032, USA
| | | | | | | | | | | |
Collapse
|
139
|
Fraser HB. Modularity and evolutionary constraint on proteins. Nat Genet 2005; 37:351-2. [PMID: 15750592 DOI: 10.1038/ng1530] [Citation(s) in RCA: 158] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2004] [Accepted: 02/03/2005] [Indexed: 11/09/2022]
Abstract
Modularity, which has been found in the functional and physical protein interaction networks of many organisms, has been postulated to affect both the mode and tempo of evolution. Here I show that in the yeast Saccharomyces cerevisiae, protein interaction hubs situated in single modules are highly constrained, whereas those connecting different modules are more plastic. This pattern of change could reflect a tendency for evolutionary innovations to occur by altering the proteins and interactions between rather than within modules, in a manner somewhat similar to the evolution of new proteins through the shuffling of conserved protein domains.
Collapse
Affiliation(s)
- Hunter B Fraser
- Department of Molecular and Cell Biology, University of California, Berkeley, California 94720, USA.
| |
Collapse
|
140
|
Abstract
This article discusses the most recent achievements in understanding the biological implications of the small-world and scale-free global topological properties of genetic, proteomic and metabolic networks. Most importantly, these networks are highly clustered and have small node-to-node distances. With their few very connected nodes, which are statistically unlikely to fail under random conditions, the proper functioning of these systems is maintained under external perturbations.
Collapse
Affiliation(s)
- Martin G Grigorov
- Nestlé Research Center, BioAnalytical Science, CH-1000 Lausanne 26, Switzerland.
| |
Collapse
|
141
|
Jordan IK, Mariño-Ramírez L, Koonin EV. Evolutionary significance of gene expression divergence. Gene 2005; 345:119-26. [PMID: 15716085 PMCID: PMC1859841 DOI: 10.1016/j.gene.2004.11.034] [Citation(s) in RCA: 104] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2004] [Revised: 11/08/2004] [Accepted: 11/15/2004] [Indexed: 11/25/2022]
Abstract
Recent large-scale studies of evolutionary changes in gene expression among mammalian species have led to the proposal that gene expression divergence may be neutral with respect to organismic fitness. Here, we employ a comparative analysis of mammalian gene sequence divergence and gene expression divergence to test the hypothesis that the evolution of gene expression is predominantly neutral. Two models of neutral gene expression evolution are considered: 1-purely neutral evolution (i.e., no selective constraint) of gene expression levels and patterns and 2-neutral evolution accompanied by selective constraint. With respect to purely neutral evolution, levels of change in gene expression between human-mouse orthologs are correlated with levels of gene sequence divergence that are determined largely by purifying selection. In contrast, evolutionary changes of tissue-specific gene expression profiles do not show such a correlation with sequence divergence. However, divergence of both gene expression levels and profiles are significantly lower for orthologous human-mouse gene pairs than for pairs of randomly chosen human and mouse genes. These data clearly point to the action of selective constraint on gene expression divergence and are inconsistent with the purely neutral model; however, there is likely to be a neutral component in evolution of gene expression, particularly, in tissues where the expression of a given gene is low and functionally irrelevant. The model of neutral evolution with selective constraint predicts a regular, clock-like accumulation of gene expression divergence. However, relative rate tests of the divergence among human-mouse-rat orthologous gene sets reveal clock-like evolution for gene sequence divergence, and to a lesser extent for gene expression level divergence, but not for the divergence of tissue-specific gene expression profiles. Taken together, these results indicate that gene expression divergence is subject to the effects of purifying selective constraint and suggest that it might also be substantially influenced by positive Darwinian selection.
Collapse
Affiliation(s)
- I King Jordan
- National Center for Biotechnology Information, National Institutes of Health 8600 Rockville Pike, Bldg 38A/Room 5N511-M, Bethesda, MD 20894, USA.
| | | | | |
Collapse
|