Genome Synteny Research Articles

Finished Genome Sequence of Collimonas arenae Cal35.

We announce the finished genome sequence of soil forest isolate Collimonas arenae Cal35, which comprises a 5.6-Mbp chromosome and 41-kb plasmid. The Cal35 genome is the second one published for the bacterial genus Collimonas and represents the first opportunity for high-resolution comparison of genome content and synteny among collimonads.

Evolution of small prokaryotic genomes.

As revealed by genome sequencing, the biology of prokaryotes with reduced genomes is strikingly diverse. These include free-living prokaryotes with ∼800 genes as well as endosymbiotic bacteria with as few as ∼140 genes. Comparative genomics is revealing the evolutionary mechanisms that led to these small genomes. In the case of free-living prokaryotes, natural selection directly favored genome reduction, while in the case of endosymbiotic prokaryotes neutral processes played a more prominent role. However, new experimental data suggest that selective processes may be at operation as well for endosymbiotic prokaryotes at least during the first stages of genome reduction. Endosymbiotic prokaryotes have evolved diverse strategies for living with reduced gene sets inside a host-defined medium. These include utilization of host-encoded functions (some of them coded by genes acquired by gene transfer from the endosymbiont and/or other bacteria); metabolic complementation between co-symbionts; and forming consortiums with other bacteria within the host. Recent genome sequencing projects of intracellular mutualistic bacteria showed that previously believed universal evolutionary trends like reduced G+C content and conservation of genome synteny are not always present in highly reduced genomes. Finally, the simplified molecular machinery of some of these organisms with small genomes may be used to aid in the design of artificial minimal cells. Here we review recent genomic discoveries of the biology of prokaryotes endowed with small gene sets and discuss the evolutionary mechanisms that have been proposed to explain their peculiar nature.

Plant Comparative and Functional Genomics.

Plant Comparative and Functional Genomics.

Expression and evolution of hexamerins from the tobacco hornworm, Manduca sexta, and other Lepidoptera

Hexamerins are large hemolymph-proteins that accumulate during the late larval stages of insects. Hexamerins have emerged from hemocyanin, but have lost the ability to bind oxygen. Hexamerins are mainly considered as storage proteins for non-feeding stages, but may also have other functions, e.g. in cuticle formation, transport and immune response. The genome of the hornworm Manduca sexta harbors six hexamerin genes. Two of them code for arylphorins (Msex2.01690, Msex2.15504) and two genes correspond to a methionine-rich hexamerin (Msex2.10735) and a moderately methionine-rich hexamerin (Msex2.01694), respectively. Two other genes do not correspond to any known hexamerin and distantly resemble the arylphorins (Msex2.01691, Msex2.01693). Five of the six hexamerin genes are clustered within ∼45 kb on scaffold 00023, which shows conserved synteny in various lepidopteran genomes. The methionine-rich hexamerin gene is located at a distinct site. M. sexta and other Lepidoptera have lost the riboflavin-binding hexamerin. With the exception of Msex2.01691, which displays low mRNA levels throughout the life cycle, all hexamerins are most highly expressed during pre-wandering phase of the 5th larval instar of M. sexta, supporting their role as storage proteins. Notably, Msex2.01691 is most highly expressed in the brain, suggesting a divergent function. Phylogenetic analyses showed that hexamerin evolution basically follows insect systematics. Lepidoptera display an unparalleled diversity of hexamerins, which exceeds that of other hexapod orders. In contrast to previous analyses, the lepidopteran hexamerins were found monophyletic. Five distinct types of hexamerins have been identified in this order, which differ in terms of amino acid composition and evolutionary history: i. the arylphorins, which are rich in aromatic amino acids (∼20% phenylalanine and tyrosine), ii. the distantly related arylphorin-like hexamerins, iii. the methionine-rich hexamerins, iv. the moderately methionine rich hexamerins, and v. the riboflavin-binding hexamerins.

Genome-wide macrosynteny among Fusarium species in the Gibberella fujikuroi complex revealed by amplified fragment length polymorphisms.

The Gibberella fujikuroi complex includes many Fusarium species that cause significant losses in yield and quality of agricultural and forestry crops. Due to their economic importance, whole-genome sequence information has rapidly become available for species including Fusarium circinatum, Fusarium fujikuroi and Fusarium verticillioides, each of which represent one of the three main clades known in this complex. However, no previous studies have explored the genomic commonalities and differences among these fungi. In this study, a previously completed genetic linkage map for an interspecific cross between Fusarium temperatum and F. circinatum, together with genomic sequence data, was utilized to consider the level of synteny between the three Fusarium genomes. Regions that are homologous amongst the Fusarium genomes examined were identified using in silico and pyrosequenced amplified fragment length polymorphism (AFLP) fragment analyses. Homology was determined using BLAST analysis of the sequences, with 777 homologous regions aligned to F. fujikuroi and F. verticillioides. This also made it possible to assign the linkage groups from the interspecific cross to their corresponding chromosomes in F. verticillioides and F. fujikuroi, as well as to assign two previously unmapped supercontigs of F. verticillioides to probable chromosomal locations. We further found evidence of a reciprocal translocation between the distal ends of chromosome 8 and 11, which apparently originated before the divergence of F. circinatum and F. temperatum. Overall, a remarkable level of macrosynteny was observed among the three Fusarium genomes, when comparing AFLP fragments. This study not only demonstrates how in silico AFLPs can aid in the integration of a genetic linkage map to the physical genome, but it also highlights the benefits of using this tool to study genomic synteny and architecture.

Asexual propagation of a virulent clone complex in a human and feline outbreak of sporotrichosis.

Sporotrichosis is one of the most frequent subcutaneous fungal infections in humans and animals caused by members of the plant-associated, dimorphic genus Sporothrix. Three of the four medically important Sporothrix species found in Brazil have been considered asexual as no sexual stage has ever been reported in Sporothrix schenckii, Sporothrix brasiliensis, or Sporothrix globosa. We have identified the mating type (MAT) loci in the S. schenckii (strain 1099-18/ATCC MYA-4821) and S. brasiliensis (strain 5110/ATCC MYA-4823) genomes by using comparative genomic approaches to determine the mating type ratio in these pathogen populations. Our analysis revealed the presence of a MAT1-1 locus in S. schenckii while a MAT1-2 locus was found in S. brasiliensis representing genomic synteny to other Sordariomycetes. Furthermore, the components of the mitogen-activated protein kinase (MAPK)-pheromone pathway, pheromone processing enzymes, and meiotic regulators have also been identified in the two pathogens, suggesting the potential for sexual reproduction. The ratio of MAT1-1 to MAT1-2 was not significantly different from 1:1 for all three Sporothrix species, but the population of S. brasiliensis in the outbreaks originated from a single mating type. We also explored the population genetic structure of these pathogens using sequence data of two loci to improve our knowledge of the pattern of geographic distribution, genetic variation, and virulence phenotypes. Population genetics data showed significant population differentiation and clonality with a low level of haplotype diversity in S. brasiliensis isolates from different regions of sporotrichosis outbreaks in Brazil. In contrast, S. schenckii isolates demonstrated a high degree of genetic variability without significant geographic differentiation, indicating the presence of recombination. This study demonstrated that two species causing the same disease have contrasting reproductive strategies and genetic variability patterns.

High-Quality Draft Genome Sequence of Babesia divergens, the Etiological Agent of Cattle and Human Babesiosis.

Babesia divergens causes significant morbidity and mortality in cattle and splenectomized or immunocompromised individuals. Here, we present a 10.7-Mb high-quality draft genome of this parasite close to chromosome resolution that will enable comparative genome analyses and synteny studies among related parasites.

Lrp13 is a novel vertebrate lipoprotein receptor that binds vitellogenins in teleost fishes

Transcripts encoding a novel member of the lipoprotein receptor superfamily, termed LDL receptor-related protein (Lrp)13, were sequenced from striped bass (Morone saxatilis) and white perch (Morone americana) ovaries. Receptor proteins were purified from perch ovary membranes by protein-affinity chromatography employing an immobilized mixture of vitellogenins Aa and Ab. RT-PCR revealed lrp13 to be predominantly expressed in striped bass ovary, and in situ hybridization detected lrp13 transcripts in the ooplasm of early secondary growth oocytes. Quantitative RT-PCR confirmed peak lrp13 expression in the ovary during early secondary growth. Quantitative mass spectrometry revealed peak Lrp13 protein levels in striped bass ovary during late-vitellogenesis, and immunohistochemistry localized Lrp13 to the oolemma and zona radiata of vitellogenic oocytes. Previously unreported orthologs of lrp13 were identified in genome sequences of fishes, chicken (Gallus gallus), mouse (Mus musculus), and dog (Canis lupus familiaris). Zebrafish (Danio rerio) and Nile tilapia (Oreochromis niloticus) lrp13 loci are discrete and share genomic synteny. The Lrp13 appears to function as a vitellogenin receptor and may be an important mediator of yolk formation in fishes and other oviparous vertebrates. The presence of lrp13 orthologs in mammals suggests that this lipoprotein receptor is widely distributed among vertebrates, where it may generally play a role in lipoprotein metabolism.

Understanding the enormous diversity of bacteriophages: The tailed phages that infect the bacterial family Enterobacteriaceae

Bacteriophages are the predominant biological entity on the planet. The recent explosion of sequence information has made estimates of their diversity possible. We describe the genomic comparison of 337 fully sequenced tailed phages isolated on 18 genera and 31 species of bacteria in the Enterobacteriaceae. These phages were largely unambiguously grouped into 56 diverse clusters (32 lytic and 24 temperate) that have syntenic similarity over >50% of the genomes within each cluster, but substantially less sequence similarity between clusters. Most clusters naturally break into sets of more closely related subclusters, 78% of which are correlated with their host genera. The largest groups of related phages are superclusters united by genome synteny to lambda (81 phages) and T7 (51 phages). This study forms a robust framework for understanding diversity and evolutionary relationships of existing tailed phages, for relating newly discovered phages and for determining host/phage relationships.

An Atypical Orphan Carbohydrate-NRPS Genomic Island Encodes a Novel Lytic Transglycosylase

Microbial genome sequencing platforms have produced a deluge of orphan biosynthetic pathways suspected of biosynthesizing new small molecules with pharmacological relevance. Genome synteny analysis provides an assessment of genomic island content, which is enriched in natural product gene clusters. Here we identified an atypical orphan carbohydrate-nonribosomal peptide synthetase genomic island in Photorhabdus luminescens using genome synteny analysis. Heterologous expression of the pathway led to the characterization of five oligosaccharide metabolites with lysozyme inhibitory activities. The oligosaccharides harbor a 1,6-anhydro-β-D-N-acetyl-glucosamine moiety, a rare structural feature for natural products. Gene deletion analysis and biochemical reconstruction of oligosaccharide production led to the discovery that a hypothetical protein in the pathway is a lytic transglycosylase responsible for bicyclic sugar formation. The example presented here supports the notion that targeting select genomic islands with reduced reliance on known protein homologies could enhance the discovery of new metabolic chemistry and biology.

Genome sequencing of chimpanzee malaria parasites reveals possible pathways of adaptation to human hosts.

Plasmodium falciparum causes most human malaria deaths, having prehistorically evolved from parasites of African Great Apes. Here we explore the genomic basis of P. falciparum adaptation to human hosts by fully sequencing the genome of the closely related chimpanzee parasite species P. reichenowi, and obtaining partial sequence data from a more distantly related chimpanzee parasite (P. gaboni). The close relationship between P. reichenowi and P. falciparum is emphasized by almost complete conservation of genomic synteny, but against this strikingly conserved background we observe major differences at loci involved in erythrocyte invasion. The organization of most virulence-associated multigene families, including the hypervariable var genes, is broadly conserved, but P. falciparum has a smaller subset of rif and stevor genes whose products are expressed on the infected erythrocyte surface. Genome-wide analysis identifies other loci under recent positive selection, but a limited number of changes at the host–parasite interface may have mediated host switching.

Pseudomonas syringae CC1557: a highly virulent strain with an unusually small type III effector repertoire that includes a novel effector.

Both type III effector proteins and nonribosomal peptide toxins play important roles for Pseudomonas syringae pathogenicity in host plants, but whether and how these pathways interact to promote infection remains unclear. Genomic evidence from one clade of P. syringae suggests a tradeoff between the total number of type III effector proteins and presence of syringomycin, syringopeptin, and syringolin A toxins. Here, we report the complete genome sequence from P. syringae CC1557, which contains the lowest number of known type III effectors to date and has also acquired genes similar to sequences encoding syringomycin pathways from other strains. We demonstrate that this strain is pathogenic on Nicotiana benthamiana and that both the type III secretion system and a new type III effector, hopBJ1, contribute to pathogenicity. We further demonstrate that activity of HopBJ1 is dependent on residues structurally similar to the catalytic site of Escherichia coli CNF1 toxin. Taken together, our results provide additional support for a negative correlation between type III effector repertoires and the potential to produce syringomycin-like toxins while also highlighting how genomic synteny and bioinformatics can be used to identify and characterize novel virulence proteins.

Genomics and Proteomics Provide New Insight into the Commensal and Pathogenic Lifestyles of Bovine- and Human-Associated Staphylococcus epidermidis Strains.

The present study reports comparative genomics and proteomics of Staphylococcus epidermidis (SE) strains isolated from bovine intramammary infection (PM221) and human hosts (ATCC12228 and RP62A). Genome-level profiling and protein expression analyses revealed that the bovine strain and the mildly infectious ATCC12228 strain are highly similar. Their genomes share high sequence identity and synteny, and both were predicted to encode the commensal-associated fdr marker gene. In contrast, PM221 was judged to differ from the sepsis-associated virulent human RP62A strain on the basis of distinct protein expression patterns and overall lack of genome synteny. The 2D DIGE and phenotypic analyses suggest that PM221 and ATCC12228 coordinate the TCA cycle activity and the formation of small colony variants in a way that could result in increased viability. Pilot experimental infection studies indicated that although ATCC12228 was able to infect a bovine host, the PM221 strain caused more severe clinical signs. Further investigation revealed strain- and condition-specific differences among surface bound proteins with likely roles in adhesion, biofilm formation, and immunomodulatory functions. Thus, our findings revealed a close link between the bovine and commensal-type human strains and suggest that humans could act as a reservoir of bovine mastitis-causing SE strains.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes.

BackgroundThe bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes.DescriptionWe have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers.ConclusionsWith a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

Genome-Wide Identification, Evolution and Expression Analysis of the Grape (Vitis vinifera L.) Zinc Finger-Homeodomain Gene Family

Plant zinc finger-homeodomain (ZHD) genes encode a family of transcription factors that have been demonstrated to play an important role in the regulation of plant growth and development. In this study, we identified a total of 13 ZHD genes (VvZHD) in the grape genome that were further classified into at least seven groups. Genome synteny analysis revealed that a number of VvZHD genes were present in the corresponding syntenic blocks of Arabidopsis, indicating that they arose before the divergence of these two species. Gene expression analysis showed that the identified VvZHD genes displayed distinct spatiotemporal expression patterns, and were differentially regulated under various stress conditions and hormone treatments, suggesting that the grape VvZHDs might be also involved in plant response to a variety of biotic and abiotic insults. Our work provides insightful information and knowledge about the ZHD genes in grape, which provides a framework for further characterization of their roles in regulation of stress tolerance as well as other aspects of grape productivity.

Looking for the bird Kiss: evolutionary scenario in sauropsids

BackgroundThe neuropeptide Kiss and its receptor KissR are key-actors in the brain control of reproduction in mammals, where they are responsible for the stimulation of the activity of GnRH neurones. Investigation in other vertebrates revealed up to 3 Kiss and 4 KissR paralogs, originating from the two rounds of whole genome duplication in early vertebrates. In contrast, the absence of Kiss and KissR has been suggested in birds, as no homologs of these genes could be found in current genomic databases. This study aims at addressing the question of the existence, from an evolutionary perspective, of the Kisspeptin system in birds. It provides the first large-scale investigation of the Kisspeptin system in the sauropsid lineage, including ophidian, chelonian, crocodilian, and avian lineages.ResultsSauropsid Kiss and KissR genes were predicted from multiple genome and transcriptome databases by TBLASTN. Phylogenetic and syntenic analyses were performed to classify predicted sauropsid Kiss and KissR genes and to re-construct the evolutionary scenarios of both gene families across the sauropsid radiation.Genome search, phylogenetic and synteny analyses, demonstrated the presence of two Kiss genes (Kiss1 and Kiss2 types) and of two KissR genes (KissR1 and KissR4 types) in the sauropsid lineage. These four genes, also present in the mammalian lineage, would have been inherited from their common amniote ancestor. In contrast, synteny analyses supported that the other Kiss and KissR paralogs are missing in sauropsids as in mammals, indicating their absence in the amniote lineage. Among sauropsids, in the avian lineage, we demonstrated the existence of a Kiss2-like gene in three bird genomes. The divergence of these avian Kiss2-like sequences from those of other vertebrates, as well as their absence in the genomes of some other birds, revealed the processes of Kiss2 gene degeneration and loss in the avian lineage.ConclusionThese findings contribute to trace back the evolutionary history of the Kisspeptin system in amniotes and sauropsids, and provide the first molecular evidence of the existence and fate of a Kiss gene in birds.

Differential Locus Expansion Distinguishes Toxoplasmatinae Species and Closely Related Strains of Toxoplasma gondii

ABSTRACTToxoplasma gondii is a human obligate intracellular parasite that has infected over 20% of the world population and has a vast intermediate host range compared to those of its nearest relatives Neospora caninum and Hammondia hammondi. While these 3 species have highly syntenic genomes (80 to 99%), in this study we examined and compared species-specific structural variations, specifically at loci that have undergone local (i.e., tandem) duplication and expansion. To do so, we used genomic sequence coverage analysis to identify and curate T. gondii and N. caninum loci that have undergone duplication and expansion (expanded loci [ELs]). The 53 T. gondii ELs are significantly enriched for genes with predicted signal sequences and single-exon genes and genes that are developmentally regulated at the transcriptional level. We validated 24 T. gondii ELs using comparative genomic hybridization; these data suggested significant copy number variation at these loci. High-molecular-weight Southern blotting for 3 T. gondii ELs revealed that copy number varies across T. gondii lineages and also between members of the same clonal lineage. Using similar methods, we identified 64 N. caninum ELs which were significantly enriched genes belonging to the SAG-related surface (SRS) antigen family. Moreover, there is significantly less overlap (30%) between the expanded gene sets in T. gondii and N. caninum than would be predicted by overall genomic synteny (81%). Consistent with this finding, only 59% of queried T. gondii ELs are similarly duplicated/expanded in H. hammondi despite over 99% genomic synteny between these species.

Genomics and breeding of Brassicaceae crops

Brassicaceae crops include oil crops, vegetables, condiments, fodder crops, and ornamental plants. Most of these crops belong to or are closely related to the genus Brassica. The Brassicaceae also includes the most intensively studied model plant, Arabidopsis thaliana. In molecular genetic studies of Brassicaceae crops, information on A. thaliana genes is highly useful, and many genes important for plant breeding have been identified. Different species of Brassica have various numbers of chromosomes. The genome relationship of three monogenomic species and three digenomic species that are well known as U’s triangle was discovered in Japan, and relationships between different genera have also been studied (Mizushima 1980). In these studies, interspecific and intergeneric hybridizations have been performed, and some interspecific and intergeneric hybrids have been used as breeding materials. Studies on interspecific and intergeneric hybridization are reviewed by Kaneko and Bang in this special issue. Genome synteny relationships between different species in Brassica and closely related genera are complicated, as exemplified by the different numbers of chromosomes between species. The monogenomic species of Brassica and closely related genera have three copies, on average, of genes homologous to each A. thaliana gene, indicating whole genome triplication in the evolution of the ancestral species of Brassica and closely related genera. The recent development of next-generation sequencers has accelerated the determination of whole genome sequences from many species. A draft sequence of the Brassica rapa genome was published in 2011 (Wang et al. 2011), and whole genome sequences of other species of Brassicaceae crops will also be published soon. The progress with genome studies of Brassicaceae crops and genome relationships between them are described by Ashutosh et al. Many Brassicaceae crops have self-incompatibility, which is not common in other seed-propagated crops. Self-incompatibility is a commonly found trait in many angiosperm species, but self-compatible plants have been selected as crops. The genetics and molecular mechanisms of self-incompatibility have been most intensively studied using Brassica species, and this trait is also used as a tool for producing F1 hybrid seeds of Brassicaceae vegetables. Recent studies on Brassica self-incompatibility are introduced by Kitashiba and Nasrallah. On the other hand, the trait most widely used for F1 hybrid seed production in many crop species is cytoplasmic male sterility. Cytoplasmic male-sterile lines can be obtained by interspecific hybridization and repeated backcrossing with a paternal species, and this form of sterility is considered to be caused by incompatibility between the cytoplasm and nucleus. Cytoplasmic male sterility in Brassicaceae crops is discussed here by Yamagishi and Bhat. Brassicaceae crops generally contain glucosinolates, which are degraded to isothiocyanates by the endogenous enzyme myrosinase after disruption of cells. Isothiocyanates are pungent compounds, probably important for defense against herbivores, and they determine the taste or quality of vegetables, condiments, and feed. These compounds are reported to have anti-carcinogenic activity; therefore, Brassica vegetables are in the spotlight for medicinal uses. Studies on glucosinolates in Brassicaceae vegetables are presented by Ishida et al. Because of water deficit in farmland, probably due to global warming, drought tolerance has become increasingly important in rapeseed breeding. Drought is frequently accompanied by salinity; therefore, recent studies on rapeseed breeding for drought and salt tolerance are reviewed by Zhang et al. Furthermore, herbicide-tolerant cultivars developed by genetic engineering are an increasing percentage of rapeseed produced around the world. For the production of these genetically modified (GM) cultivars, pollen-mediated gene flow from GM cultivars to wild plants or landraces should be limited. Recent studies on the possibility of pollen-mediated gene flow from GM cultivars of rapeseed are described by Tsuda et al. After the Great East Japan Earthquake, the Rapeseed Project for Restoring Tsunami-Salt-Damaged Farmland was started by Tohoku University because of the high tolerance of Brassica crops to salinity. This project has been supported by Kureha Corporation and other companies. The cost of publication of these reviews is partially supported by this project.

Channel catfish hemoglobin genes: Identification, phylogenetic and syntenic analysis, and specific induction in response to heat stress

Hemoglobins transport oxygen from gill to inner organs in fish, and this process is affected by temperature, one of the major environmental factors for fish. The hemoglobin gene clusters have been well studied in humans and several model fish species, but remain largely unknown in catfish. Here, eight α- and six β-hemoglobin genes were identified and characterized in channel catfish. Genomic synteny analysis showed that these hemoglobin genes were separated into two unlinked clusters, the MN cluster containing six α- and six β-hemoglobin genes, and the LA cluster consisting of two α-hemoglobin genes. Channel catfish hemoglobin genes were ubiquitously expressed in all the 10 tested tissues from healthy fish, but exhibited higher expression level in spleen, head kidney, and trunk kidney. In response to heat stress, hemoglobin genes, especially MN Hbα4, MN Hbα5, MN Hbα6, MN Hbβ4, MN Hbβ5, MN Hbβ6, LA Hbα1, and LA Hbα2, presumably the embryonic hemoglobin genes, were drastically up-regulated in the gill and head kidney of heat-tolerant fishes, but not in these tissues of the heat-intolerant fish, suggesting the importance of the embryonic hemoglobin genes in coping with the low oxygen conditions under heat stress.

Evolution of the gastrin–cholecystokinin gene family revealed by synteny analysis

Gastrin (GAST) and cholecystokinin (CCK) are two structurally and functionally related peptide hormones that exert many functions, including regulation of gastric and pancreatic secretion, feeding behaviour and energy homeostasis. GAST and CCK genes are assumed to have diverged from a common ancestral gene, over 500 million years ago in the vertebrate lineage. However, although a large number of GAST and CCK-related sequences have been identified both in vertebrate and non-vertebrate species, the evolutionary history of the GAST/CCK family remains little understood. To address this issue, we used extensive genome synteny comparisons of vertebrate chromosomes, in particular to evaluate the impact of whole-genome duplications. In the present study, we confirm that the GAST/CCK family in vertebrates is composed of two paralogous genes, namely GAST and CCK, and even three in teleosts, namely GAST, CCK1 and CCK2. We also show that the GAST and CCK genes arose by duplications of a single ancestral gene through the 2R and that the two copies of the CCK gene found in teleosts have probably been generated through the 3R. Finally, our results suggest that the vertebrate ancestor possessed four members of the GAST/CCK family, of which two have likely been lost during evolution.