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.     

Proneural genes and the specification of neural cell types

Certain morphological, physiological and molecular characteristics are shared by all neurons. However, despite these similarities, neurons constitute the most diverse cell population of any organism. Recently, considerable attention has been focused on identifying the molecular mechanisms that underlie this cellular diversity. Parallel studies in Drosophila and vertebrates have revealed that proneural genes are key regulators of neurogenesis, coordinating the acquisition of a generic neuronal fate and of specific subtype identities that are appropriate for the location and time of neuronal generation. These studies reveal that, in spite of differences between invertebrate and vertebrate neural lineages, Drosophila and vertebrate proneural genes have remarkably similar roles.     

Do mosquitoes filter the access of Plasmodium cytochrome b lineages to an avian host?

Many parasites show fidelity to a set of hosts in ecological time but not evolutionary time and the determinants of this pattern are poorly understood. Malarial parasites use vertebrate hosts for the asexual stage of their life cycle but use Dipteran hosts for the sexual stage. Despite the potential evolutionary importance of Dipteran hosts, little is known of their role in determining a parasite's access to vertebrate hosts. Here, we use an avian malarial system in Panama to explore whether mosquitoes act as an access filter that limits the range of vertebrate hosts used by particular parasite lineages. We amplified and sequenced Plasmodium mitochondrial DNA (mtDNA) from Turdus grayi (clay-coloured robin) and from mosquitoes at the same study site. We trapped and identified to species 123 141 female mosquitoes and completed polymerase chain reaction (PCR) screening for Plasmodium parasites in 435 pools of 20 mosquitoes per pool (8700 individuals total) spanning the 11 most common mosquito species. Our primers amplified nine Plasmodium lineages, whose sequences differed by 1.72%–10.0%. Phylogenetic analyses revealed partial clustering of lineages that co-occurred in mosquito hosts. However PAN3 and PAN6, the two primary parasite lineages of T. grayi , exhibited sequence divergence of 8.59% and did not cluster in the phylogeny. We detected these two lineages exclusively in mosquitoes from different genera — PAN3 was found only in Culex (Melanoconion) ocossa , and PAN6 was found only in Aedeomyia squamipennis . Furthermore, each of these two parasite lineages co-occurred in mosquitoes with other Plasmodium lineages that were not found in the vertebrate host T. grayi . Together, this evidence suggests that parasite–mosquito associations do not restrict the access of parasites to birds but instead may actually facilitate the switching of vertebrate hosts that occurs over evolutionary time.     

Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani

The complete nucleotide sequence (16,650 bp) of the mitochondrial genome of the salamander Mertensiella luschani (Caudata, Amphibia) was determined. This molecule conforms to the consensus vertebrate mitochondrial gene order. However, it is characterized by a long non-coding intervening sequence with two 124-bp repeats between the tRNA(Thr) and tRNA(Pro) genes. The new sequence data were used to reconstruct a phylogeny of jawed vertebrates. Phylogenetic analyses of all mitochondrial protein-coding genes at the amino acid level recovered a robust vertebrate tree in which lungfishes are the closest living relatives of tetrapods, salamanders and frogs are grouped together to the exclusion of caecilians (the Batrachia hypothesis) in a monophyletic amphibian clade, turtles show diapsid affinities and are placed as sister group of crocodiles+birds, and the marsupials are grouped together with monotremes and basal to placental mammals. The deduced phylogeny was used to characterize the molecular evolution of vertebrate mitochondrial proteins. Amino acid frequencies were analyzed across the main lineages of jawed vertebrates, and leucine and cysteine were found to be the most and least abundant amino acids in mitochondrial proteins, respectively. Patterns of amino acid replacements were conserved among vertebrates. Overall, cartilaginous fishes showed the least variation in amino acid frequencies and replacements. Constancy of rates of evolution among the main lineages of jawed vertebrates was rejected.     

爬行动物MHC基因研究进展

主要组织相容性复合体(MHC)是有颌脊椎动物中发现的编码免疫球蛋白受体的高度多态的基因群,因其在免疫系统中的重要作用而备受关注。脊椎动物不同支系间MHC的结构和演化差异较大。尽管MHC基因特征在哺乳类、鸟类、两栖类和鱼类中已被较好地描述,但对爬行动物MHC的了解仍较少。鉴于爬行动物对于理解MHC基因的演化占据很重要的系统发育位置,研究其MHC具有重要意义。本文就近年来爬行动物MHC的分子结构、多态性维持机制、功能和主要应用的研究现状进行了系统地回顾和总结,并展望了其研究前景。     

Long-Lived Dichotomous Lineages of the Proteasome Subunit Beta Type 8 (PSMB8) Gene Surviving More than 500 Million Years as Alleles or Paralogs

On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans-species dimorphism of the proteasome subunit beta type 8 (PSMB8) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A(31)F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.     

Amphioxus (Branchiostoma floridae) has orthologs of vertebrate odorant receptors

Background  

A common feature of chemosensory systems is the involvement of G protein-coupled receptors (GPCRs) in the detection of environmental stimuli. Several lineages of GPCRs are involved in vertebrate olfaction, including trace amine-associated receptors, type 1 and 2 vomeronasal receptors and odorant receptors (ORs). Gene duplication and gene loss in different vertebrate lineages have lead to an enormous amount of variation in OR gene repertoire among species; some fish have fewer than 100 OR genes, while some mammals possess more than 1000. Fascinating features of the vertebrate olfactory system include allelic exclusion, where each olfactory neuron expresses only a single OR gene, and axonal guidance where neurons expressing the same receptor project axons to common glomerulae. By identifying homologous ORs in vertebrate and in non-vertebrate chordates, we hope to expose ancestral features of the chordate olfactory system that will help us to better understand the evolution of the receptors themselves and of the cellular components of the olfactory system.     

The amphioxus genome enlightens the evolution of the thyroid hormone signaling pathway

Thyroid hormones (THs) have pleiotropic effects on vertebrate development, with amphibian metamorphosis as the most spectacular example. However, developmental functions of THs in non-vertebrate chordates are largely hypothetical and even TH endogenous production has been poorly investigated. In order to get better insight into the evolution of the thyroid hormone signaling pathway in chordates, we have taken advantage of the recent release of the amphioxus genome. We found amphioxus homologous sequences to most of the genes encoding proteins involved in thyroid hormone signaling in vertebrates, except the fast-evolving thyroglobulin: sodium iodide symporter, thyroid peroxidase, deiodinases, thyroid hormone receptor, TBG, and CTHBP. As only some genes encoding proteins involved in TH synthesis regulation were retrieved (TRH, TSH receptor, and CRH receptor but not their corresponding receptors and ligands), there may be another mode of upstream regulation of TH synthesis in amphioxus. In accord with the notion that two whole genome duplications took place at the base of the vertebrate tree, one amphioxus gene often corresponded to several vertebrate homologs. However, some amphioxus specific duplications occurred, suggesting that several steps of the TH pathway were independently elaborated in the cephalochordate and vertebrate lineages. The present results therefore indicate that amphioxus is capable of producing THs. As several genes of the TH signaling pathway were also found in the sea urchin genome, we propose that the thyroid hormone signaling pathway is of ancestral origin in chordates, if not in deuterostomes, with specific elaborations in each lineage, including amphioxus.     

Comparative structure analysis of vertebrate ribonuclease P RNA.

Ribonuclease P cleaves 5'-precursor sequences from pre-tRNAs. All cellular RNase P holoenzymes contain homologous RNA elements; the eucaryal RNase P RNA, in contrast to the bacterial RNA, is catalytically inactive in the absence of the protein component(s). To understand the function of eucaryal RNase P RNA, knowledge of its structure is needed. Considerable effort has been devoted to comparative studies of the structure of this RNA from diverse organisms, including eucaryotes, primarily fungi, but also a limited set of vertebrates. The substantial differences in the sequences and structures of the vertebrate RNAs from those of other organisms have made it difficult to align the vertebrate sequences, thus limiting comparative studies. To expand our understanding of the structure of diverse RNase P RNAs, we have isolated by PCR and sequenced 13 partial RNase P RNA genes from 11 additional vertebrate taxa representing most extant major vertebrate lineages. Based on a recently proposed structure of the core elements of RNase P RNA, we aligned the sequences and propose a minimum consensus secondary structure for the vertebrate RNase P RNA.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

Hox gene duplication in fish

The study of Hox gene clusters continues to serve as a paradigm for those interested in vertebrate genome evolution. Recent exciting discoveries about Hox gene composition in fishes challenges conventional views about vertebrate Hox gene evolution, and has initiated lively debates concerning the evolutionary events making the divergence of the major vertebrate lineages. Comparative analyses indicate that Hox cluster duplications occurred in early vertebrate evolution, and again within the order Cypriniformes of teleost fish. Loss of Hox genes was more widespread than duplication during fish evolution.     

Patterns and consequences of vertebrate Emx gene duplications

SUMMARY We have cloned and analyzed two Emx genes from the lamprey Petromyzon marinus and our findings provide insight into the patterns and developmental consequences of gene duplications during early vertebrate evolution. The Emx gene family presents an excellent case for addressing these issues as gnathostome vertebrates possess two or three Emx paralogs that are highly pleiotropic, functioning in or being expressed during the development of several vertebrate synapomorphies. Lampreys are the most primitive extant vertebrates and characterization of their development and genomic organization is critical for understanding vertebrate origins. We identified two Emx genes from P. marinus and analyzed their phylogeny and their embryological expression relative to other chordate Emx genes. Our phylogenetic analysis shows that the two lamprey Emx genes group independently from the gnathostome Emx1, Emx2 , and Emx3 paralogy groups. Our expression analysis shows that the two lamprey Emx genes are expressed in distinct spatial and temporal patterns that together broadly encompass the combined sites of expression of all gnathostome Emx genes. Our data support a model wherein large-scale regulatory evolution of a single Emx gene occurred after the protochordate/vertebrate divergence, but before the vertebrate radiation. Both the lamprey and gnathostome lineages then underwent independent gene duplications followed by extensive paralog subfunctionalization. Emx subfunctionalization in the telencephalon is remarkably convergent and refines our understanding of lamprey forebrain patterning. We also identify lamprey-specific sites of expression that indicate either neofunctionalization in lampreys or sites-specific nonfunctionalization of all gnathostome Emx genes. Overall, we see only very limited correlation between Emx gene duplications and the acquisition of novel expression domains.     

Expression of three <Emphasis Type="Italic">spalt </Emphasis>(<Emphasis Type="Italic">sal</Emphasis>) gene homologues in zebrafish embryos

Three homologues of the Drosophilaregion-specific homeotic gene spalt (sal) have been isolated in zebrafish, sall1a, sall1b and sall3. Phylogenetic analysis of these genes against known salDNA sequences showed zebrafish sall1aand sall1b to be orthologous to other vertebrate sal-1 genes and zebrafish sall3to be orthologous to other vertebrate sal-3 genes, except Xenopus sall3. Phylogenetic reconstruction suggests that zebrafish sall1a and sall1bresulted from a gene duplication event occurring prior to the divergence of the ray-finned and lobe-finned fish lineages. Analysis of the expression pattern of the zebrafish sal genes shows that sall1a and sall3 share expression domains with both orthologous and non-orthologous vertebrate sal genes. Both are expressed in various regions of the CNS, including in primary motor neurons. Outside of the CNS, sall1a expression is observed in the otic vesicle (ear), heart and in a discrete region of the pronephric ducts. These analyses indicate that orthologies between zebrafish sal genes and other vertebrate sal genes do not imply equivalence of expression pattern and, therefore, that biological functions are not entirely conserved. However we suggest that, like other vertebrate sal genes, zebrafish sal genes have a role in neural development. Also, expression of zebrafish sall1a in the otic vesicle, heart sac and the pronephric ducts of zebrafish embryos is possibly consistent with some of the abnormalities seen in Sall1-deficient mice and in Townes-Brocks Syndrome, a human disorder which is caused by mutations in the human spalt gene SALL1.     

Microchromosomes Exhibit Distinct Features of Vertebrate Chromosome Structure and Function with Underappreciated Ramifications for Genome Evolution

Microchromosomes are common yet poorly understood components of many vertebrate genomes. Recent studies have revealed that microchromosomes contain a high density of genes and possess other distinct characteristics compared with macrochromosomes. Whether distinctive characteristics of microchromosomes extend to features of genome structure and organization, however, remains an open question. Here, we analyze Hi-C sequencing data from multiple vertebrate lineages and show that microchromosomes exhibit consistently high degrees of interchromosomal interaction (particularly with other microchromosomes), appear to be colocalized to a common central nuclear territory, and are comprised of a higher proportion of open chromatin than macrochromosomes. These findings highlight an unappreciated level of diversity in vertebrate genome structure and function, and raise important questions regarding the evolutionary origins and ramifications of microchromosomes and the genes that they house.     

Evolution of gliding in Southeast Asian geckos and other vertebrates is temporally congruent with dipterocarp forest development

Gliding morphologies occur in diverse vertebrate lineages in Southeast Asian rainforests, including three gecko genera, plus frogs, snakes, agamid lizards and squirrels. It has been hypothesized that repeated evolution of gliding is related to the dominance of Asian rainforest tree floras by dipterocarps. For dipterocarps to have influenced the evolution of gliding in Southeast Asian vertebrates, gliding lineages must have Eocene or later origins. However, divergence times are not known for most lineages. To investigate the temporal pattern of Asian gliding vertebrate evolution, we performed phylogenetic and molecular clock analyses. New sequence data for geckos incorporate exemplars of each gliding genus (Cosymbotus, Luperosaurus and Ptychozoon), whereas analyses of other vertebrate lineages use existing sequence data. Stem ages of most gliding vertebrates, including all geckos, cluster in the time period when dipterocarps came to dominate Asian tropical forests. These results demonstrate that a gliding/dipterocarp correlation is temporally viable, and caution against the assumption of early origins for apomorphic taxa.     

Adaptive evolution of the uncoupling protein 1 gene contributed to the acquisition of novel nonshivering thermogenesis in ancestral eutherian mammals

Homeotherms possess various physiological mechanisms to maintain their body temperature, thus allowing them to adapt to various environments. Under cold conditions, most eutherian mammals upregulate heat production in brown adipose tissue (BAT), and uncoupling protein (UCP) 1 is an essential factor in BAT thermogenesis. The evolutionary origin of UCP1 was believed to have been a specific event occurring in eutherian lineages. Recently, however, the UCP1 ortholog was found in fishes, which uncovers a more ancient origin of this gene than previously believed. Here we investigate the evolutionary process of UCP1 by comparative genomic approach. We found that UCP1 evolved rapidly by positive Darwinian selection in the common ancestor of eutherians, although this gene arose in the ancestral vertebrate, since the orthologous genes were shared among most of the vertebrate species. Adaptive evolution occurred after the divergence between eutherians and marsupials, which is consistent with the fact that BAT has been found only in eutherians. Our findings indicate that positive Darwinian selection acted on UCP1 contributed to the acquisition of an efficient mechanism for body temperature regulation in primitive eutherians. Phylogenetic reconstruction of UCP1 with two paralogs (UCP2 and UCP3) among vertebrate species revealed that the gene duplication events which produced these three genes occurred in the common ancestor of vertebrates much earlier than the emergence of eutherians. Thus, our data demonstrate that novel gene function can evolve without de novo gene duplication event.     

Zebrafish ribonucleases are bactericidal: implications for the origin of the vertebrate RNase A superfamily

Understanding the evolutionary origin of the ribonuclease (RNase) A superfamily is of great interest because the superfamily is the sole vertebrate-specific enzyme family known to date. Although mammalian RNases have a diverse array of biochemical and physiological functions, the original function of the superfamily at its birth is enigmatic. Such information may be obtained by studying basal lineages of the vertebrate phylogeny and is necessary for discerning how and why this superfamily originated. Here, we clone and characterize 3 RNase genes from the zebrafish, the most basal vertebrate examined for RNases. We report 1) that all the 3 zebrafish RNases are ribonucleolytically active, with one of them having an RNase activity comparable to that of bovine RNase A, the prototype of the superfamily; 2) that 2 zebrafish RNases have prominent expressions in adult liver and gut, whereas the 3rd is expressed in adult eye and heart; and 3) that all 3 RNases have antibacterial activities in vitro. These results, together with the presence of antibacterial and/or antiviral activities in multiple distantly related mammalian RNases, strongly suggest that the superfamily started as a host-defense mechanism in vertebrate evolution.