The caspase family in urochordates: distinct evolutionary fates in ascidians and larvaceans

BACKGROUND INFORMATION: Caspases are cysteine proteases that mediate apoptosis (programmed cell death) initiation and execution. Apoptosis is a conserved mechanism shared by all metazoans, although its physiological function and complexity show considerable taxon-dependent variations. To gain insight into the caspase repertoire of putative ancestors to vertebrates, we performed exhaustive genomic searches in urochordates, a sister taxon to vertebrates in which ascidians and appendicularians display chordate characters at early stages of their development. RESULTS: We identified the complete caspase families of two ascidians (Ciona intestinalis and C. savignyi) and one larvacean (Oikopleura dioica). We found in ascidian species an extremely high number of caspase genes (17 for C. intestinalis and 22 for C. savignyi), deriving from five founder gene orthologues to human pro-inflammatory, initiator and executioner caspases. Although considered to be sibling species, C. intestinalis and C. savignyi only share 11 orthologues, most of the additional genes resulting from recent mass duplications. A sharply contrasted picture was found in O. dioica, which displayed only three caspase genes deriving from a single founder gene distantly related to caspase 3/7. The difference between ascidian and larvacean caspase repertoires is discussed in the light of their developmental patterns and life cycles. CONCLUSIONS: The identification of caspase members in two ascidian species delineates five founder genes that bridge the gap between vertebrates and Ecdysozoa (arthropods and nematodes). Given the amazing diversity among urochordates, determination and comparison of the caspase repertoires in species from orders additional to Enterogona (ascidians) and Oikopleuridae might be highly informative on the evolution of caspase-dependent physiological processes.     

Evolution of ABC transporters by gene duplication and their role in human disease

The ATP-binding cassette (ABC) transporter genes represent the largest family of transporters and these genes are abundant in the genome of all vertebrates. Through analysis of the genome sequence databases we have characterized the full complement of ABC genes from several mammals and other vertebrates. Multiple gene duplication and deletion events were identified in ABC genes in different lineages indicating that the process of gene evolution is still ongoing. Gene duplication resulting in either gene birth or gene death plays a major role in the evolution of the vertebrate ABC genes. The understanding of this mechanism is important in the context of human health because these ABC genes are associated with human disease, involving nearly all organ systems of the body. In addition, ABC genes play an important role in the development of drug resistance in cancer cells. Future genetic, functional, and evolutionary studies of ABC transporters will provide important insight into human and animal biology.     

Evolutionary Origin of GnIH and NPFF in Chordates: Insights from Novel Amphioxus RFamide Peptides

Gonadotropin-inhibitory hormone (GnIH) is a newly identified hypothalamic neuropeptide that inhibits pituitary hormone secretion in vertebrates. GnIH has an LPXRFamide (X = L or Q) motif at the C-terminal in representative species of gnathostomes. On the other hand, neuropeptide FF (NPFF), a neuropeptide characterized as a pain-modulatory neuropeptide, in vertebrates has a PQRFamide motif similar to the C-terminal of GnIH, suggesting that GnIH and NPFF have diverged from a common ancestor. Because GnIH and NPFF belong to the RFamide peptide family in vertebrates, protochordate RFamide peptides may provide important insights into the evolutionary origin of GnIH and NPFF. In this study, we identified a novel gene encoding RFamide peptides and two genes of their putative receptors in the amphioxus Branchiostoma japonicum. Molecular phylogenetic analysis and synteny analysis indicated that these genes are closely related to the genes of GnIH and NPFF and their receptors of vertebrates. We further identified mature RFamide peptides and their receptors in protochordates. The identified amphioxus RFamide peptides inhibited forskolin induced cAMP signaling in the COS-7 cells with one of the identified amphioxus RFamide peptide receptors expressed. These results indicate that the identified protochordate RFamide peptide gene is a common ancestral form of GnIH and NPFF genes, suggesting that the origin of GnIH and NPFF may date back to the time of the emergence of early chordates. GnIH gene and NPFF gene may have diverged by whole-genome duplication in the course of vertebrate evolution.     

Similar numbers but different repertoires of olfactory receptor genes in humans and chimpanzees

Animals recognize their external world through the detection of tens of thousands of chemical odorants. Olfactory receptor (OR) genes encode proteins for detecting odorant molecules and form the largest multigene family in mammals. It is known that humans have fewer OR genes and a higher fraction of OR pseudogenes than mice or dogs. To investigate whether these features are human specific or common to all higher primates, we identified nearly complete sets of OR genes from the chimpanzee and macaque genomes and compared them with the human OR genes. In contrast to previous studies, here we show that the number of OR genes ( approximately 810) and the fraction of pseudogenes (51%) in chimpanzees are very similar to those in humans, though macaques have considerably fewer OR genes. The pseudogenization rates and the numbers of genes affected by positive selection are also similar between humans and chimpanzees. Moreover, the most recent common ancestor between humans and chimpanzees had a larger number of functional OR genes (>500) and a lower fraction of pseudogenes (41%) than its descendents, suggesting that the OR gene repertoires are in a phase of deterioration in both lineages. Interestingly, despite the close evolutionary relationship between the 2 species, approximately 25% of their functional gene repertoires are species specific due to massive gene losses. These findings suggest that the tempo of evolution of OR genes is similar between humans and chimpanzees, but the OR gene repertoires are quite different between them. This difference might be responsible for the species-specific ability of odor perception.     

Immune Related Genes Underpin the Evolution of Adaptive Immunity in Jawless Vertebrates

The study of immune related genes in lampreys and hagfish provides a unique perspective on the evolutionary genetic underpinnings of adaptive immunity and the evolution of vertebrate genomes. Separated from their jawed cousins at the stem of the vertebrate lineage, these jawless vertebrates have many of the gene families and gene regulatory networks associated with the defining morphological and physiological features of vertebrates. These include genes vital for innate immunity, inflammation, wound healing, protein degradation, and the development, signaling and trafficking of lymphocytes. Jawless vertebrates recognize antigen by using leucine-rich repeat (LRR) based variable lymphocyte receptors (VLRs), which are very different from the immunoglobulin (Ig) based T cell receptor (TCR) and B cell receptor (BCR) used for antigen recognition by jawed vertebrates. The somatically constructed VLR genes are expressed in monoallelic fashion by T-like and B-like lymphocytes. Jawless and jawed vertebrates thus share many of the genes that provide the molecular infrastructure and physiological context for adaptive immune responses, yet use entirely different genes and mechanisms of combinatorial assembly to generate diverse repertoires of antigen recognition receptors.     

The Dlx genes as clues to vertebrate genomics and craniofacial evolution

The group of Dlx genes belongs to the homeobox-containing superfamily, and its members are involved in various morphogenetic processes. In vertebrate genomes, Dlx genes exist as multiple paralogues generated by tandem duplication followed by whole genome duplications. In this review, we provide an overview of the Dlx gene phylogeny with an emphasis on the chordate lineage. Referring to the Dlx gene repertoire, we discuss the establishment and conservation of the nested expression patterns of the Dlx genes in craniofacial development. Despite the accumulating genomic sequence resources in diverse vertebrates, embryological analyses of Dlx gene expression and function remain limited in terms of species diversity. By supplementing our original analysis of shark embryos with previous data from other osteichthyans, such as mice and zebrafish, we support the previous speculation that the nested Dlx expression in the pharyngeal arch is likely a shared feature among all the extant jawed vertebrates. Here, we highlight several hitherto unaddressed issues regarding the evolution and function of Dlx genes, with special reference to the craniofacial development of vertebrates.     

A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori

A genome-wide analysis of innate immunity-related genes and gene families was conducted using the silkworm, Bombyx mori. We identified orthologs for a large number of genes involved in insect immunity that have been reported from Drosophila melanogaster (Diptera), Anopheles gambiae (Diptera), Apis mellifera (Hymenoptera) and Tribolium castaneum (Coleoptera). B. mori has a unique recognition gene and antimicrobial peptide genes that are not present in the Drosophila, Anopheles, Apis and Tribolium genomes, suggesting a lineage-specific gene evolution for lepidopteran insects. The comparative analysis of the insect immune repertoires indicated a dynamic and flexible gene expansion in recognition, modulation and effector mechanisms due to different selection pressures. Differential gene regulation by different bacterial species was found in PGRP and Serpin genes, suggesting that Bombyx has a highly selective gene regulation system depending on bacterial species.     

Gene Expression in Chicken Reveals Correlation with Structural Genomic Features and Conserved Patterns of Transcription in the Terrestrial Vertebrates

Background

The chicken is an important agricultural and avian-model species. A survey of gene expression in a range of different tissues will provide a benchmark for understanding expression levels under normal physiological conditions in birds. With expression data for birds being very scant, this benchmark is of particular interest for comparative expression analysis among various terrestrial vertebrates.

Methodology/Principal Findings

We carried out a gene expression survey in eight major chicken tissues using whole genome microarrays. A global picture of gene expression is presented for the eight tissues, and tissue specific as well as common gene expression were identified. A Gene Ontology (GO) term enrichment analysis showed that tissue-specific genes are enriched with GO terms reflecting the physiological functions of the specific tissue, and housekeeping genes are enriched with GO terms related to essential biological functions. Comparisons of structural genomic features between tissue-specific genes and housekeeping genes show that housekeeping genes are more compact. Specifically, coding sequence and particularly introns are shorter than genes that display more variation in expression between tissues, and in addition intergenic space was also shorter. Meanwhile, housekeeping genes are more likely to co-localize with other abundantly or highly expressed genes on the same chromosomal regions. Furthermore, comparisons of gene expression in a panel of five common tissues between birds, mammals and amphibians showed that the expression patterns across tissues are highly similar for orthologuous genes compared to random gene pairs within each pair-wise comparison, indicating a high degree of functional conservation in gene expression among terrestrial vertebrates.

Conclusions

The housekeeping genes identified in this study have shorter gene length, shorter coding sequence length, shorter introns, and shorter intergenic regions, there seems to be selection pressure on economy in genes with a wide tissue distribution, i.e. these genes are more compact. A comparative analysis showed that the expression patterns of orthologous genes are conserved in the terrestrial vertebrates during evolution.     

Comparative study of the <Emphasis Type="Italic">P2X</Emphasis> gene family in animals and plants

P2X receptors are ligand-gated ion channels that can bind with the adenosine triphosphate (ATP) and have diverse functional roles in neuropathic pain, inflammation, special sense, and so on. In this study, 180 putative P2X genes, including 176 members in 32 animal species and 4 members in 3 species of lower plants, were identified. These genes were divided into 13 groups, including 7 groups in vertebrates and 6 groups in invertebrates and lower plants, through phylogenetic analysis. Their gene organization and motif composition are conserved in most predicted P2X members, while group-specific features were also found. Moreover, synteny relationships of the putative P2X genes in vertebrates are conserved while simultaneously experiencing a series of gene insertion, inversion, and transposition. Recombination signals were detected in almost all of the vertebrates and invertebrates, suggesting that intragenic recombination may play a significant role in the evolution of P2X genes. Selection analysis also identified some positively selected sites that acted on the evolution of most of the predicted P2X proteins. The phenomenon of alternative splicing occurred commonly in the putative P2X genes of vertebrates. This article explored in depth the evolutional relationship among different subtypes of P2X genes in animal and plants and might serve as a solid foundation for deciphering their functions in further studies.     

Coevolution of immunoglobulin heavy- and light-chain variable-region gene families

The gene families encoding the immunoglobulin variable regions of heavy (VH) and light (VL) chains in vertebrates are composed of many genes. However, the gene number and the extent of diversity among VH and VL gene copies vary with species. To examine the causes of this variation and the evolutionary forces for these multigene families, we conducted a phylogenetic analysis of VH and VL genes from the species of amniotes. The results of our analysis showed that for each species, VH and VL genes have the same pattern of clustering in the trees, and, according to this clustering pattern, the species can be divided into two groups. In the first group of species (humans and mice), VH and VL genes were extensively intermingled with genes from other organisms; in the second group of species (chickens, rabbits, cattle, sheep, swine, and horses), the genes tended to form clusters within the same group of organisms. These results suggest that the VH and VL multigene families have evolved in the same fashion: they have undergone coordinated contraction and expansion of gene repertoires such that each group of organisms is characterized by a certain level of diversity of VH and VL genes. The extent of diversity among copies of VH and VL genes in each species is related to the mechanism of generation of antibody variety. In humans and mice, DNA rearrangement of immunoglobulin variable, diversity, and joining-segment genes is a main source of antibody diversity, whereas in chickens, rabbits, cattle, sheep, swine, and horses, somatic hypermutation and somatic gene conversion play important roles. The evolutionary pattern of VH and VL multigene families is consistent with the birth-and-death model of evolution, yet different levels of diversifying selection seem to operate in the VH and VL genes of these two groups of species.      

Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats

Pheromones are used by individuals of the same species to elicit behavioral or physiological changes, and they are perceived primarily by the vomeronasal organ (VNO) in terrestrial vertebrates. VNO pheromone receptors are encoded by the V1r and V2r gene superfamilies in mammals. A comparison of the V1r and V2r repertoires between closely related species can provide significant insights into the evolutionary genetic mechanisms responsible for species-specific pheromone communications. A total of 137 putatively functional V1r genes of 12 families were previously identified from the mouse genome. We report the identification of 95 putatively functional V1r genes from the draft rat genome sequence. These genes map primarily to four blocks in two chromosomes. The rat V1r genes can be phylogenetically grouped into 10 families, which are shared with mouse, and 2 new families, which are rat-specific. Even in many shared families, gene numbers differ between the two species, apparently due to frequent gene duplication and pseudogenization after the separation of the two species. Molecular dating suggests that most of the rat V1r families emerged before or during the radiation of mammalian orders, but many duplications within families occurred as recently as in the past 10 million years (MY). Our results show that the evolution of the V1r repertoire is characterized by exceptionally fast gene turnover via gains and losses of individual genes, suggesting rapid and substantial changes in pheromone communication between species.     

Evolution of genes involved in gamete interaction: evidence for positive selection, duplications and losses in vertebrates

Genes encoding proteins involved in sperm-egg interaction and fertilization exhibit a particularly fast evolution and may participate in prezygotic species isolation [1], [2]. Some of them (ZP3, ADAM1, ADAM2, ACR and CD9) have individually been shown to evolve under positive selection [3], [4], suggesting a role of positive Darwinian selection on sperm-egg interaction. However, the genes involved in this biological function have not been systematically and exhaustively studied with an evolutionary perspective, in particular across vertebrates with internal and external fertilization. Here we show that 33 genes among the 69 that have been experimentally shown to be involved in fertilization in at least one taxon in vertebrates are under positive selection. Moreover, we identified 17 pseudogenes and 39 genes that have at least one duplicate in one species. For 15 genes, we found neither positive selection, nor gene copies or pseudogenes. Genes of teleosts, especially genes involved in sperm-oolemma fusion, appear to be more frequently under positive selection than genes of birds and eutherians. In contrast, pseudogenization, gene loss and gene gain are more frequent in eutherians. Thus, each of the 19 studied vertebrate species exhibits a unique signature characterized by gene gain and loss, as well as position of amino acids under positive selection. Reflecting these clade-specific signatures, teleosts and eutherian mammals are recovered as clades in a parsimony analysis. Interestingly the same analysis places Xenopus apart from teleosts, with which it shares the primitive external fertilization, and locates it along with amniotes (which share internal fertilization), suggesting that external or internal environmental conditions of germ cell interaction may not be the unique factors that drive the evolution of fertilization genes. Our work should improve our understanding of the fertilization process and on the establishment of reproductive barriers, for example by offering new leads for experiments on genes identified as positively selected.     

Increased complexity of gene structure and base composition in vertebrates

How the structure and base composition of genes changed with the evolution of vertebrates remains a puzzling question. Here we analyzed 895 orthologous protein-coding genes in six multicellular animals: human, chicken, zebrafish, sea squirt, fruit fly, and worm. Our analyses reveal that many gene regions, particularly intron and 3' UTR, gradually expanded throughout the evolution of vertebrates from their invertebrate ancestors, and that the number of exons per gene increased. Studies based on all protein-coding genes in each genome provide consistent results.We also find that GC-content increased in many gene regions (especially 5' UTR) in the evolution of endotherms, except in coding-exons.Analysis of individual genomes shows that 3′ UTR demonstrated stronger length and CC-content correlation with intron than 5' UTR, and gene with large intron in all six species demonstrated relatively similar GC-content. Our data indicates a great increase in complexity in vertebrate genes and we propose that the requirement for morphological and functional changes is probably the driving force behind the evolution of structure and base composition complexity in multicellular animal genes.     

Population genomics of the immune evasion (var) genes of Plasmodium falciparum

Var genes encode the major surface antigen (PfEMP1) of the blood stages of the human malaria parasite Plasmodium falciparum. Differential expression of up to 60 diverse var genes in each parasite genome underlies immune evasion. We compared the diversity of the DBLalpha domain of var genes sampled from 30 parasite isolates from a malaria endemic area of Papua New Guinea (PNG) and 59 from widespread geographic origins (global). Overall, we obtained over 8,000 quality-controlled DBLalpha sequences. Within our sampling frame, the global population had a total of 895 distinct DBLalpha "types" and negligible overlap among repertoires. This indicated that var gene diversity on a global scale is so immense that many genomes would need to be sequenced to capture its true extent. In contrast, we found a much lower diversity in PNG of 185 DBLalpha types, with an average of approximately 7% overlap among repertoires. While we identify marked geographic structuring, nearly 40% of types identified in PNG were also found in samples from different countries showing a cosmopolitan distribution for much of the diversity. We also present evidence to suggest that recombination plays a key role in maintaining the unprecedented levels of polymorphism found in these immune evasion genes. This population genomic framework provides a cost effective molecular epidemiological tool to rapidly explore the geographic diversity of var genes.     

Evolution of different Y chromosomes in two medaka species, Oryzias dancena and O. latipes

Although the sex-determining gene DMY has been identified on the Y chromosome in the medaka (Oryzias latipes), this gene is absent in most Oryzias species, suggesting that closely related species have different sex-determining genes. Here, we investigated the sex-determination mechanism in O. dancena, which does not possess the DMY gene. Since heteromorphic sex chromosomes have not been reported in this species, a progeny test of sex-reversed individuals produced by hormone treatment was performed. Sex-reversed males yielded all-female progeny, indicating that O. dancena has an XX/XY sex-determination system. To uncover the cryptic sex chromosomes, sex-linked DNA markers were screened using expressed sequence tags (ESTs) established in O. latipes. Linkage analysis of isolated sex-linked ESTs showed a conserved synteny between the sex chromosomes in O. dancena and an autosome in O. latipes. Fluorescence in situ hybridization (FISH) analysis of these markers confirmed that sex chromosomes of these species are not homologous. These findings strongly suggest an independent origin of sex chromosomes in O. dancena and O. latipes. Further analysis of the sex-determining region in O. dancena should provide crucial insights into the evolution of sex-determination mechanisms in vertebrates.     

Characterization of a divergent non-classical MHC class I gene in sharks

Sharks are the most ancient group of vertebrates known to possess members of the major histocompatibility complex (MHC) gene family. For this reason, sharks provide a unique opportunity to gain insight into the evolution of the vertebrate immune system through comparative analysis. Two genes encoding proteins related to the MHC class I gene family were isolated from splenic cDNA derived from spiny dogfish shark ( Squalus acanthias). The genes have been designated MhcSqac-UAA*01 and MhcSqac-UAA*NC1. Comparative analysis demonstrates that the Sqac-UAA*01 protein sequence clusters with classical MHC class I of several shark species and has structural elements common to most classical MHC class I molecules. In contrast, Sqac-UAA*NC1 is highly divergent from all vertebrate classical MHC class I proteins, including the Sqac-UAA *01 sequence and those of other shark species. Although Sqac-UAA*NC1 is clearly related to the MHC class I gene family, no orthologous genes from other species were identified due to the high degree of sequence divergence. In fact, the Sqac NC1 protein sequence is the most divergent MHC class-I-like protein identified thus far in any shark species. This high degree of divergence is similar in magnitude to some of the MHC class-I-related genes found in mammals, such as MICA or CD1. These data support the existence of a class of highly divergent non-classical MHC class I genes in the most primitive vertebrates known to possess homologues of the MHC and other components of the adaptive immune system.