Evolutionary Signatures amongst Disease Genes Permit Novel Methods for Gene Prioritization and Construction of Informative Gene-Based Networks

Genes involved in the same function tend to have similar evolutionary histories, in that their rates of evolution covary over time. This coevolutionary signature, termed Evolutionary Rate Covariation (ERC), is calculated using only gene sequences from a set of closely related species and has demonstrated potential as a computational tool for inferring functional relationships between genes. To further define applications of ERC, we first established that roughly 55% of genetic diseases posses an ERC signature between their contributing genes. At a false discovery rate of 5% we report 40 such diseases including cancers, developmental disorders and mitochondrial diseases. Given these coevolutionary signatures between disease genes, we then assessed ERC''s ability to prioritize known disease genes out of a list of unrelated candidates. We found that in the presence of an ERC signature, the true disease gene is effectively prioritized to the top 6% of candidates on average. We then apply this strategy to a melanoma-associated region on chromosome 1 and identify MCL1 as a potential causative gene. Furthermore, to gain global insight into disease mechanisms, we used ERC to predict molecular connections between 310 nominally distinct diseases. The resulting “disease map” network associates several diseases with related pathogenic mechanisms and unveils many novel relationships between clinically distinct diseases, such as between Hirschsprung''s disease and melanoma. Taken together, these results demonstrate the utility of molecular evolution as a gene discovery platform and show that evolutionary signatures can be used to build informative gene-based networks.     

On the origins of Mendelian disease genes in man: the impact of gene duplication

Over 3,000 human diseases are known to be linked to heritable genetic variation, mapping to over 1,700 unique genes. Dating of the evolutionary age of these disease-associated genes has suggested that they have a tendency to be ancient, specifically coming into existence with early metazoa. The approach taken by past studies, however, assumes that the age of a disease is the same as the age of its common ancestor, ignoring the fundamental contribution of duplication events in the evolution of new genes and function. Here, we date both the common ancestor and the duplication history of known human disease-associated genes. We find that the majority of disease genes (80%) are genes that have been duplicated in their evolutionary history. Periods for which there are more disease-associated genes, for example, at the origins of bony vertebrates, are explained by the emergence of more genes at that time, and the majority of these are duplicates inferred to have arisen by whole-genome duplication. These relationships are similar for different disease types and the disease-associated gene's cellular function. This indicates that the emergence of duplication-associated diseases has been ongoing and approximately constant (relative to the retention of duplicate genes) throughout the evolution of life. This continued until approximately 390 Ma from which time relatively fewer novel genes came into existence on the human lineage, let alone disease genes. For single-copy genes associated with disease, we find that the numbers of disease genes decreases with recency. For the majority of duplicates, the disease-associated mutation is associated with just one of the duplicate copies. A universal explanation for heritable disease is, thus, that it is merely a by-product of the evolutionary process; the evolution of new genes (de novo or by duplication) results in the potential for new diseases to emerge.     

An ancient evolutionary origin of genes associated with human genetic diseases

Several thousand genes in the human genome have been linked to a heritable genetic disease. The majority of these appear to be nonessential genes (i.e., are not embryonically lethal when inactivated), and one could therefore speculate that they are late additions in the evolutionary lineage toward humans. Contrary to this expectation, we find that they are in fact significantly overrepresented among the genes that have emerged during the early evolution of the metazoa. Using a phylostratigraphic approach, we have studied the evolutionary emergence of such genes at 19 phylogenetic levels. The majority of disease genes was already present in the eukaryotic ancestor, and the second largest number has arisen around the time of evolution of multicellularity. Conversely, genes specific to the mammalian lineage are highly underrepresented. Hence, genes involved in genetic diseases are not simply a random subset of all genes in the genome but are biased toward ancient genes.     

Modification of gene duplicability during the evolution of protein interaction network

Duplications of genes encoding highly connected and essential proteins are selected against in several species but not in human, where duplicated genes encode highly connected proteins. To understand when and how gene duplicability changed in evolution, we compare gene and network properties in four species (Escherichia coli, yeast, fly, and human) that are representative of the increase in evolutionary complexity, defined as progressive growth in the number of genes, cells, and cell types. We find that the origin and conservation of a gene significantly correlates with the properties of the encoded protein in the protein-protein interaction network. All four species preserve a core of singleton and central hubs that originated early in evolution, are highly conserved, and accomplish basic biological functions. Another group of hubs appeared in metazoans and duplicated in vertebrates, mostly through vertebrate-specific whole genome duplication. Such recent and duplicated hubs are frequently targets of microRNAs and show tissue-selective expression, suggesting that these are alternative mechanisms to control their dosage. Our study shows how networks modified during evolution and contributes to explaining the occurrence of somatic genetic diseases, such as cancer, in terms of network perturbations.     

Conservation of regional gene expression in mouse and human brain

Many neurodegenerative diseases have a hallmark regional and cellular pathology. Gene expression analysis of healthy tissues may provide clues to the differences that distinguish resistant and sensitive tissues and cell types. Comparative analysis of gene expression in healthy mouse and human brain provides a framework to explore the ability of mice to model diseases of the human brain. It may also aid in understanding brain evolution and the basis for higher order cognitive abilities. Here we compare gene expression profiles of human motor cortex, caudate nucleus, and cerebellum to one another and identify genes that are more highly expressed in one region relative to another. We separately perform identical analysis on corresponding brain regions from mice. Within each species, we find that the different brain regions have distinctly different expression profiles. Contrasting between the two species shows that regionally enriched genes in one species are generally regionally enriched genes in the other species. Thus, even when considering thousands of genes, the expression ratios in two regions from one species are significantly correlated with expression ratios in the other species. Finally, genes whose expression is higher in one area of the brain relative to the other areas, in other words genes with patterned expression, tend to have greater conservation of nucleotide sequence than more widely expressed genes. Together these observations suggest that region-specific genes have been conserved in the mammalian brain at both the sequence and gene expression levels. Given the general similarity between patterns of gene expression in healthy human and mouse brains, we believe it is reasonable to expect a high degree of concordance between microarray phenotypes of human neurodegenerative diseases and their mouse models. Finally, these data on very divergent species provide context for studies in more closely related species that address questions such as the origins of cognitive differences.     

Adaptive evolution of the Populus tremula photoperiod pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light and genes in the photoperiodic pathway are known to be involved in controlling these processes. Here, we examine 25 of genes from the photoperiod pathway in Populus tremula (Salicaceae) for signatures of adaptive evolution. Overall, levels of synonymous polymorphism in the 25 genes are lower than at control loci selected randomly from the genome. This appears primarily to be caused by lower levels of synonymous polymorphism in genes associated with the circadian clock. Natural selection appears to play an important role in shaping protein evolution at several of the genes in the photoperiod pathways, which is highlighted by the fact that approximately 40% of the genes from the photoperiod pathway have estimates of selection on nonsynonymous polymorphisms that are significantly different from zero. A surprising observation we make is that circadian clock-associated genes appear to be over-represented among the genes showing elevated rates of protein evolution; seven genes are evolving under positive selection and all but one of these genes are involved in the circadian clock of Populus.     

Noncoding RNAs within the HOX gene network in tumor pathogenesis and progression

HOX genes are involved with normal development, cell identity, cell differentiation, cell metabolism, apoptosis, autophagy as well as with diseases such as tumor pathogenesis and progression. In particular, the genes belonging to HOX paralogous 13 seem to carry out a relevant role in both tumor development and disease progression. In recent years, several noncoding RNAs (ncRNA) sequences have been identified in HOX loci, including long noncoding RNA (lncRNA) and microRNA (miRNA), highly conserved during evolution. Many studies have shown that specific intergenic ncRNAs in HOX loci could directly modulate HOX genes expression in normal and pathological conditions. In the present review we attempt to describe the role of these ncRNAs, through the regulation of the HOX gene network, in normal cell biology, and, with particular emphasis, in diseases such as in cancer pathogenesis and progression.     

Colon cancer associated genes exhibit signatures of positive selection at functionally significant positions

ABSTRACT: BACKGROUND: Cancer, much like most human disease, is routinely studied by utilizing model organisms. Of these model organisms, mice are often dominant. However, our assumptions of functional equivalence fail to consider the opportunity for divergence conferred by ~180 Million Years (MY) of independent evolution between these species. For a given set of human disease related genes, it is therefore important to determine if functional equivalency has been retained between species. In this study we test the hypothesis that cancer associated genes have different patterns of substitution akin to adaptive evolution in different mammal lineages. RESULTS: Our analysis of the current literature and colon cancer databases identified 22 genes exhibiting colon cancer associated germline mutations. We identified orthologs for these 22 genes across a set of high coverage (>6X) vertebrate genomes. Analysis of these orthologous datasets revealed significant levels of positive selection. Evidence of lineage-specific positive selection was identified in 14 genes in both ancestral and extant lineages. Lineage-specific positive selection was detected in the ancestral Euarchontoglires and Hominidae lineages for STK11, in the ancestral primate lineage for CDH1, in the ancestral Murinae lineage for both SDHC and MSH6 genes and the ancestral Muridae lineage for TSC1. CONCLUSION: Identifying positive selection in the primate, Hominidae, Muridae and Murinae lineages suggest an ancestral functional shift in these genes between the rodent and primate lineages. Analyses such as this, combining evolutionary theory and predictions - along with medically relevant data, can thus provide us with important clues for modeling human diseases.     

Identification and mapping of 26 human testis mRNAs containing CAG/CTG repeats

Various types of pathologies, including neurodegenerative diseases, as well as different types of neoplasia, are related to genes exhibiting simple tandem repeat instabilities. In order to seek for new candidate genes for such disorders, we screened 4.10⁶ human testis cDNAs for CAG- and CTG-containing clones. Among 910 positive clones, we characterized 109 cDNAs corresponding to 26 independent mRNAs. Fourteen of these mRNAs represent new genes. The corresponding clones contain between 3 and 19 consecutive CAG or CTG triplets. We assigned 15 out of these 26 genes to 14 different human chromosomes. These genes represent new potential candidates for diseases associated with CAG or CTG repeat mutations. Received: 5 March 1998 / Accepted: 6 May 1998     

miRNAs

The evolution of complex animals such as insects and mammals is achieved with surprisingly few additions in protein coding genes. MicroRNAs (miRNAs), a class of non-coding RNAs, have emerged as important regulators of organogenesis in insects, fish and mammals. The microRNA repertoire of animals has expanded significantly during evolution especially in vertebrates, insects and nematodes, accompanying the appearance of complex body plans. MicroRNAs therefore have gained enormous interest in recent years. They are now regarded as key modulators of gene expression in many tissues during embryogenesis, in adult organisms and in disease processes. Therefore, these small RNA molecules have entered the center stage of molecular biology and are promising candidates not only for the regulation of key biological processes such as proliferation and apoptosis, but also for therapy of human diseases.     

Genome-wide Identification and Characterization of Fixed Human-Specific Regulatory Regions

Changes in gene regulatory networks are believed to have played an important role in the development of human-specific anatomy and behavior. We identified the human genome regions that show the typical chromatin marks of regulatory regions but cannot be aligned to other mammalian genomes. Most of these regions have become fixed in the human genome. Their regulatory targets are enriched in genes involved in neural processes, CNS development, and diseases such as autism, depression, and schizophrenia. Specific transposable elements contributing to the rewiring of the human regulatory network can be identified by the creation of human-specific regulatory regions. Our results confirm the relevance of regulatory evolution in the emergence of human traits and cognitive abilities and the importance of newly acquired genomic elements for such evolution.     

拟南芥花药绒毡层细胞中具有基因簇特征的基因进化和功能分析

在植物基因组中, 除了同源基因成簇现象外, 近年来还发现一些具有共表达特性的异源基因也能够以基因簇形式存在, 但这些异源基因簇的进化和生物学功能尚不清楚。花药发育和花粉形成是植物进化出的特有的生殖生物学过程, 同时产生了一些在花药绒毡层中特异表达和特定功能的基因簇基因。该研究通过筛选和分析花药绒毡层中基因簇基因的分子特性、表达调控、基因年龄和基因重复进化等信息, 探讨花药基因簇基因与植物开花功能进化之间的关系。结果表明, 在拟南芥(Arabidopsis thaliana)中共筛选到84个(13个基因簇)花药绒毡层特异高表达的基因簇基因, 它们主要产生于串联重复事件, 76%的基因出现在开花植物分化后的阶段, 主要参与生殖发育、花粉鞘组成和脂代谢等生物学过程。研究初步解析了拟南芥花药绒毡层中基因簇基因的基本特征、生物学功能和基因进化机制, 为深入揭示植物基因簇基因的遗传学功能奠定了基础。     

Genome sequence of the cat pathogen, Chlamydophila felis.

Chlamydophila felis (Chlamydia psittaci feline pneumonitis agent) is a worldwide spread pathogen for pneumonia and conjunctivitis in cats. Herein, we determined the entire genomic DNA sequence of the Japanese C. felis strain Fe/C-56 to understand the mechanism of diseases caused by this pathogen. The C. felis genome is composed of a circular 1,166,239 bp chromosome encoding 1005 protein-coding genes and a 7552 bp circular plasmid. Comparison of C. felis gene contents with other Chlamydia species shows that 795 genes are common in the family Chlamydiaceae species and 47 genes are specific to C. felis. Phylogenetic analysis of the common genes reveals that most of the orthologue sets exhibit a similar divergent pattern but 14 C. felis genes accumulate more mutations, implicating that these genes may be involved in the evolutional adaptation to the C. felis-specific niche. Gene distribution and orthologue analyses reveal that two distinctive regions, i.e. the plasticity zone and frequently gene-translocated regions (FGRs), may play important but different roles for chlamydial genome evolution. The genomic DNA sequence of C. felis provides information for comprehension of diseases and elucidation of the chlamydial evolution.     

Ancient and recent adaptive evolution of primate non-homologous end joining genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.     

Ancient gene transfer from algae to animals: Mechanisms and evolutionary significance

ABSTRACT: BACKGROUND: Horizontal gene transfer (HGT) is traditionally considered to be rare in multicellular eukaryotes such as animals. Recently, many genes of miscellaneous algal origins were discovered in choanoflagellates. Considering that choanoflagellates are the existing closest relatives of animals, we speculated that ancient HGT might have occurred in the unicellular ancestor of animals and affected the long-term evolution of animals. RESULTS: Through genome screening, phylogenetic and domain analyses, we identified 14 gene families, including 92 genes, in the tunicate Ciona intestinalis that are likely derived from miscellaneous photosynthetic eukaryotes. Almost all of these gene families are distributed in diverse animals, suggesting that they were mostly acquired by the common ancestor of animals. Their miscellaneous origins also suggest that these genes are not derived from a particular algal endosymbiont. In addition, most genes identified in our analyses are functionally related to molecule transport, cellular regulation and methylation signaling, suggesting that the acquisition of these genes might have facilitated the intercellular communication in the ancestral animal. CONCLUSIONS: Our findings provide additional evidence that algal genes in aplastidic eukaryotes are not exclusively derived from historical plastids and thus important for interpreting the evolution of eukaryotic photosynthesis. Most importantly, our data represent the first evidence that more anciently acquired genes might exist in animals and that ancient HGT events have played an important role in animal evolution.     

CRDB: database of chemosensory receptor gene families in vertebrate

Chemosensory receptors (CR) are crucial for animals to sense the environmental changes and survive on earth. The emergence of whole-genome sequences provides us an opportunity to identify the entire CR gene repertoires. To completely gain more insight into the evolution of CR genes in vertebrates, we identified the nearly all CR genes in 25 vertebrates using homology-based approaches. Among these CR gene repertoires, nearly half of them were identified for the first time in those previously uncharacterized species, such as the guinea pig, giant panda and elephant, etc. Consistent with previous findings, we found that the numbers of CR genes vary extensively among different species, suggesting an extreme form of 'birth-and-death' evolution. For the purpose of facilitating CR gene analysis, we constructed a database with the goals to provide a resource for CR genes annotation and a web tool for exploring their evolutionary patterns. Besides a search engine for the gene extraction from a specific chromosome region, an easy-to-use phylogenetic analysis tool was also provided to facilitate online phylogeny study of CR genes. Our work can provide a rigorous platform for further study on the evolution of CR genes in vertebrates.     

Genes "waiting" for recruitment by the adaptive immune system: the insights from amphioxus

In seeking evidence of the existence of adaptive immune system (AIS) in ancient chordate, cDNA clones of six libraries from a protochordate, the Chinese amphioxus, were sequenced. Although the key molecules such as TCR, MHC, Ig, and RAG in AIS have not been identified from our database, we demonstrated in this study the extensive molecular evidence for the presence of genes homologous to many genes that are involved in AIS directly or indirectly, including some of which may represent the putative precursors of vertebrate AIS-related genes. The comparative analyses of these genes in different model organisms revealed the different fates of these genes during evolution. Their gene expression pattern suggested that the primitive digestive system is the pivotal place of the origin and evolution of the AIS. Our studies support the general statement that AIS appears after the jawless/jawed vertebrate split. However our study further reveals the fact that AIS is in its twilight in amphioxus and the evolution of the molecules in amphioxus are waiting for recruitment by the emergence of AIS.