Evolutionary conservation of expression profiles between human and mouse orthologous genes

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.     

Incongruent expression profiles between human and mouse orthologous genes suggest widespread neutral evolution of transcription control

Rapid rates of evolution can signify either a lack of selective constraint and the consequent accumulation of neutral alleles, or positive Darwinian selection driving the fixation of advantageous alleles. Based on a comparison of 1,350 orthologous gene pairs from human and mouse, we show that the evolution of gene expression profiles is so rapid that it is comparable to that of paralogous gene pairs or randomly paired genes. The expression divergence in the entire set of orthologous pairs neither strongly correlates with sequence divergence, nor focuses in any particular tissue. Moreover, comparing tissue expressions across the orthologous gene pairs, we observe that any human tissue is more similar to any other human tissue examined than to its corresponding mouse tissue. Collectively, these results indicate that, while some differences in expression profiles may be due to adaptive evolution, the levels of divergence are mostly compatible with a neutral mode of evolution, in which a mutation for ectopic expression may rise to fixation by random drift without significantly affecting the fitness. A disturbing corollary of these findings is that knowledge of where the gene is expressed may not carry information about its function.     

Evolutionary relationship between human and mouse immunoglobulin kappa light chain variable region genes

Similar to the Igh-V multigene family, the human or mouse Igk-V repertoirer is a distorted continuum of homologous genes that may be grouped into families displaying >80% nucleic acid sequence similarity among their members. systematic interspecies sequence comparisons reveal that most human Igk-V gene families exhibit clear homology to mouse Ogk-V families (sequence similarity >74%). A hypothetical phylogenetic tree of Igk-V genes predicts that a minimum of seven Igk-V genes/families predate mammalian radiation. In two cases, several interrelated mouse Igk-V families exhibit phylogenetic equidistance with just one human Igk-V family, implying a more pronounced divergence for the elevated number of Igk-V gene families in the mouse. Mouse-human Igk-V comaprisons, moreover, illustrate how expansion, contraction, and perhaps deletion of Igk-V gene families shape the Igk-V repertoire during mammalian evolution.     

Evolutionary conservation and network structure characterize genes of phenotypic relevance for mitosis in human

The impact of gene silencing on cellular phenotypes is difficult to establish due to the complexity of interactions in the associated biological processes and pathways. A recent genome-wide RNA knock-down study both identified and phenotypically characterized a set of important genes for the cell cycle in HeLa cells. Here, we combine a molecular interaction network analysis, based on physical and functional protein interactions, in conjunction with evolutionary information, to elucidate the common biological and topological properties of these key genes. Our results show that these genes tend to be conserved with their corresponding protein interactions across several species and are key constituents of the evolutionary conserved molecular interaction network. Moreover, a group of bistable network motifs is found to be conserved within this network, which are likely to influence the network stability and therefore the robustness of cellular functioning. They form a cluster, which displays functional homogeneity and is significantly enriched in genes phenotypically relevant for mitosis. Additional results reveal a relationship between specific cellular processes and the phenotypic outcomes induced by gene silencing. This study introduces new ideas regarding the relationship between genotype and phenotype in the context of the cell cycle. We show that the analysis of molecular interaction networks can result in the identification of genes relevant to cellular processes, which is a promising avenue for future research.     

Structural differences between the repertoires of mouse and human germline genes and their evolutionary implications

 Although human and mouse antibodies are similar when one considers their diversification strategies, they differ in the extent to which kappa and lambda light chains are present in their respective variable light chain repertoires. While the Igk-V germline genes are preponderant in mice (95% or more), they comprise only 60% in humans. This may account for differences in the structural repertoire encoded in the Igk-V germline genes of these species. However, this subject has not been properly investigated, partially because a systematic structural characterization of the mouse Igk-V germline genes has not been undertaken. In the present study we compiled all available information on mouse Igk-V germline genes to characterize their structural repertoire. As expected, comparison with the structural repertoire of human Igk-V germline genes indicates differences. The most interesting is that the mouse Igk-V germline gene repertoire is more diverse in structural terms than its human counterpart: the mouse encodes seven canonical structure classes (combination of canonical structures in L1 and L3). In contrast, the human encodes only four. Analysis of the evolutionary relationships of human and mouse Igk-V germline genes led us to propose that the difference reflects a strategy of mice to compensate for the small lambda chain contribution to the repertoire of their variable light chains. Received: 1 June 1997 / Revised: 6 October 1997     

Evolutionary and biochemical differences between human and monkey acidic mammalian chitinases

Acidic mammalian chitinase (AMCase), an enzyme implicated in the pathology of asthma, is capable of chitin cleavage at a low pH optimum. The corresponding gene (CHIA) can be found in genome databases of a variety of mammals, but the enzyme properties of only the human and mouse proteins were extensively studied. We wanted to compare enzymes of closely related species, such as humans and macaques. In our attempt to study macaque AMCase, we searched for CHIA-like genes in human and macaque genomes. We found that both genomes contain several additional CHIA-like sequences. In humans, CHIA-L1 (hCHIA-L1) is an apparent pseudogene and has the highest homology to CHIA. To determine which of the two genes is functional in monkeys, we assessed their tissue expression levels. In our experiments, CHIA-L1 expression was not detected in human stomach tissue, while CHIA was expressed at high levels. However, in the cynomolgus macaque stomach tissue, the expression pattern of these two genes was reversed: CHIA-L1 was expressed at high levels and CHIA was undetectable. We hypothesized that in macaques CHIA-L1 (mCHIA-L1), and not CHIA, is a gene encoding an acidic chitinase, and cloned it, using the sequence of human CHIA-L1 as a guide for the primer design. We named the new enzyme MACase (Macaca Acidic Chitinase) to emphasize its differences from AMCase. MACase shares a similar tissue expression pattern and pH optimum with human AMCase, but is 50 times more active in our enzymatic activity assay. DNA sequence of the mCHIA-L1 has higher percentage identity to the human pseudogene hCHIA-L1 (91.7%) than to hCHIA (84%). Our results suggest alternate evolutionary paths for human and monkey acidic chitinases.     

Evolutionary differences between alternative and constitutive protein-coding regions of alternatively spliced genes of Drosophila

A total of 790 Drosophila melanogaster genes that are alternatively spliced in a coding region and have orthologs in Drosophila pseudoobscura were studied. It proved that nucleotide substitutions are accumulated in alternative coding regions more rapidly than in constitutive coding regions. Moreover, the evolutionary patterns of alternative regions differing in insertion-deletion mechanisms (use of alternative promoters, splicing sites, or polyadenylation sites) differ significantly. The synonymous substitution rate in coding regions of genes varies more strongly than the nonsynonymous substitution rate. The patterns of substitutions in different classes of alternative regions of Drosophila melanogaster and mammals differ considerably.     

Comprehensive comparison between locations of orthologous genes on archaeal and bacterial genomes

MOTIVATION: Following an extensive search for orthologous genes between the complete genomes from archaea and bacteria, the spatial association of the orthologs has been investigated in terms of synteny, the conservation of the order of neighboring genes. However, the relationships between the relative locations of remote orthologs over entire genomes have not been shown. RESULTS: Comprehensive comparisons between the locations of orthologs on nineteen archaeal and bacterial genomes are presented by the location to location correspondence based on the gene-location distance. When the two genomes are rotated such that a pair of orthologs with the shortest distance is set in the same angle, a statistically significant number of orthologs maintain their relative locations between the genomes. Even by the short distances at the 5% significance level, the rotations are restricted within a narrow range, suggesting an intrinsic angle for realizing similar locations between the orthologs in each genome pair. Furthermore, the rotations in the restricted range agree with the replication origin and terminus sites for the analyzed genomes where such sites are known. The relationship between location-maintained orthologs and gene function is also discussed.     

Evolutionary modelling of feed forward loops in gene regulatory networks

Feed forward loops (FFLs) are gene regulatory network motifs. They exist in different types, defined by the signs of the effects of genes in the motif on one another. We examine 36 feed forward loops in Escherichia coli, using evolutionary simulations to predict the forms of FFL expected to evolve to generate the pattern of expression of the output gene. These predictions are tested using likelihood ratios, comparing likelihoods of the observed FFL structures with their likelihoods under null models. The very high likelihood ratios generated, of over 10(11), suggest that evolutionary simulation is a valuable component in the explanation of FFL structure.     

人与小鼠核糖体蛋白基因内含子中的转录调控位点分析

前期对酵母和果蝇核糖体蛋白(Ribosomal protein, RP)基因内含子序列中的寡核苷酸分析表明, 内含子中含有潜在的转录因子结合位点。为进一步发掘核糖体蛋白基因内含子参与转录调控的证据, 文章首先基于频率分析方法抽提出人和小鼠核糖体蛋白基因第一内含子中高频(Over-represented)出现的寡核苷酸片段 (亦称模体, Motif), 这些寡核苷酸中超过85%与已知的转录因子结合位点吻合, 是潜在的转录调控元件。对抽提出的寡核苷酸进行碱基组成分析, 发现95%以上的寡核苷酸富含碱基C和G, 而较少富含A和T。从寡核苷酸在内含子中的分布情况看, 它们相对靠近第一内含子的5′端, 即距离基因转录起始位点和上游区域较近。推测这些特征可能与基因转录调控有关。