Poor Odors, Strength, and Persistence Give Their Rewards to Mutilla europaea Visiting Dangerous Wasp Nests

Social insect colonies are attractive for many arthropods. The rare velvet ant, Mutilla europaea, visit colonies of Polistes wasps, but to date it was unclear which resources it targeted therein. Our field observations and bioassays showed that velvet ants visit the colonies of the social wasp Polistes biglumis almost undisturbed and may feed on larval wasp saliva. Chemical insignificance and resistance are the characteristics that allow velvet ants to visit unharmed wasp colonies and gain such a nutritious reward.     

Vitellogenin Underwent Subfunctionalization to Acquire Caste and Behavioral Specific Expression in the Harvester Ant Pogonomyrmex barbatus

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.     

Dracula ant phylogeny as inferred by nuclear 28S rDNA sequences and implications for ant systematics (Hymenoptera: Formicidae: Amblyoponinae)

Ants are one of the most ecologically and numerically dominant families of organisms in almost every terrestrial habitat throughout the world, though they include only about 1% of all described insect species. The development of eusociality is thought to have been a driving force in the striking diversification and dominance of this group, yet we know little about the evolution of the major lineages of ants and have been unable to clearly determine their primitive characteristics. Ants within the subfamily Amblyoponinae are specialized arthropod predators, possess many anatomically and behaviorally primitive characters and have been proposed as a possible basal lineage within the ants. We investigate the phylogenetic relationships among the members of the subfamily, using nuclear 28S rDNA sequence data. Outgroups for the analysis include members of the poneromorph and leptanillomorph (Apomyrma, Leptanilla) ant subfamilies, as well as three wasp families. Parsimony, maximum likelihood, and Bayesian analyses provide strong support for the monophyly of a clade containing the two genera Apomyrma+Mystrium (100% bpp; 97% ML bs; and 97% MP bs), and moderate support for the monophyly of the Amblyoponinae as long as Apomyrma (Apomyrminae) is included (87% bpp; 57% ML bs; and 76% MP bs). Analyses did not recover evidence of monophyly of the Amblyopone genus, while the monophyly of the other genera in the subfamily is supported. Based on these results we provide a morphological diagnosis of the Amblyoponinae that includes Apomyrma. Among the outgroup taxa, Typhlomyrmex grouped consistently with Ectatomma, supporting the recent placement of Typhlomyrmex in the Ectatomminae. The results of this present study place the included ant subfamilies into roughly two clades with the basal placement of Leptanilla unclear. One clade contains all the Amblyoponinae (including Apomyrma), Ponerinae, and Proceratiinae (Poneroid clade). The other clade contains members from subfamilies Cerapachyinae, Dolichoderinae, Ectatomminae, Formicinae, Myrmeciinae, and Myrmicinae (Formicoid clade).     

A ten-year comparative study of the population dynamics and parasitoidism in the native paper wasp Polistes fuscatus and the invasive P. dominulus

The invasive success of Polistes dominulus in North America has been attributed to its greater productivity relative to native Polistes. Liberation from parasites and parasitoids are thought to be major factors contributing to the high productivity of P. dominulus. We analyzed historical records of colony relative abundance and productivity of P. dominulus and the sympatric, native Polistes fuscatus from 1995 to 2010 using historical data from our Michigan Polistes study site. We also analyzed evidence of parasitoids from 294 P. fuscatus and 507 P. dominulus archived combs from 2001 to 2010. Additionally, we examined field and laboratory colonies from outside of our study site for parasites and parasitoids in 2009 and 2010. We documented one parasite and three parasitoids exploiting Polistes in our Michigan study sites. Our historical records document that P. dominulus initially displaced P. fuscatus rapidly, then slowed, and finally the two populations stabilized. Furthermore, the historical pattern of decreasing displacement of P. fuscatus by P. dominulus corresponded temporally with a significant decline in the productivity and a significant increase in Dibrachys cavus infestation of P. dominulus. Our evidence indicates that the parasitoid, D. cavus, is a major factor in stabilizing the populations of the sympatric P. dominulus and P. fuscatus.     

A Chromosome-Level Assembly of Blunt Snout Bream (Megalobrama amblycephala) Genome Reveals an Expansion of Olfactory Receptor Genes in Freshwater Fish

The number of olfactory receptor genes (ORs), which are responsible for detecting diverse odor molecules varies extensively among mammals as a result of frequent gene gains and losses that contribute to olfactory specialization. However, how OR expansions/contractions in fish are influenced by habitat and feeding habit and which OR subfamilies are important in each ecological niche is unknown. Here, we report a major OR expansion in a freshwater herbivorous fish, Megalobrama amblycephala, using a highly contiguous, chromosome-level assembly. We evaluate the possible contribution of OR expansion to habitat and feeding specialization by comparing the OR repertoire in 28 phylogenetically and ecologically diverse teleosts. In total, we analyzed > 4,000 ORs including 3,253 intact, 122 truncated, and 913 pseudogenes. The number of intact ORs is highly variable ranging from 20 to 279. We estimate that the most recent common ancestor of Osteichthyes had 62 intact ORs, which declined in most lineages except the freshwater Otophysa clade that has a substantial expansion in subfamily β and ε ORs. Across teleosts, we found a strong association between duplications of β and ε ORs and freshwater habitat. Nearly, all ORs were expressed in the olfactory epithelium (OE) in three tested fish species. Specifically, all the expanded β and ε ORs were highly expressed in OE of M. amblycephala. Together, we provide molecular and functional evidence for how OR repertoires in fish have undergone gain and loss with respect to ecological factors and highlight the role of β and ε OR in freshwater adaptation.     

Functional Evolution of a Bark Beetle Odorant Receptor Clade Detecting Monoterpenoids of Different Ecological Origins

Insects detect odors using an array of odorant receptors (ORs), which may expand through gene duplication. How and which new functions may evolve among related ORs within a species remain poorly investigated. We addressed this question by functionally characterizing ORs from the Eurasian spruce bark beetle Ips typographus, in which physiological and behavioral responses to pheromones, volatiles from host and nonhost trees, and fungal symbionts are well described. In contrast, knowledge of OR function is restricted to two receptors detecting the pheromone compounds (S)-(–)-ipsenol (ItypOR46) and (R)-(–)-ipsdienol (ItypOR49). These receptors belong to an Ips-specific OR-lineage comprising seven ItypORs. To gain insight into the functional evolution of related ORs, we characterized the five remaining ORs in this clade using Xenopus oocytes. Two receptors responded primarily to the host tree monoterpenes (+)-3-carene (ItypOR25) and p-cymene (ItypOR27). Two receptors responded to oxygenated monoterpenoids produced in larger relative amounts by the beetle-associated fungi, with ItypOR23 specific for (+)-trans-(1R, 4S)-4-thujanol, and ItypOR29 responding to (+)-isopinocamphone and similar ketones. ItypOR28 responded to the pheromone E-myrcenol from the competitor Ips duplicatus. Overall, the OR responses match well with those of previously characterized olfactory sensory neuron classes except that neurons detecting E-myrcenol have not been identified. The characterized ORs are under strong purifying selection and demonstrate a shared functional property in that they all primarily respond to monoterpenoids. The variation in functional groups among OR ligands and their diverse ecological origins suggest that neofunctionalization has occurred early in the evolution of this OR-lineage following gene duplication.     

Ground nesting in the paper wasp <Emphasis Type="Italic">Polistes aurifer</Emphasis> (Hymenoptera,Vespidae)

Summary Polistes paper wasp species vary in their nest site selection, with some nesting in relatively exposed areas and others in protected cavities. Locating the nest sites of cavity dwelling species in natural habitat can prove difficult, and most behavioral studies on these species are therefore conducted using human-built structures. Since Hungerford and Williams (1912) and Rau (1929) noted the location of several P. fuscatus (Fabricius) nests in rodent burrows, there have been no published accounts of ground nesting Polistes species. I report the occurrence of a large, dense cluster of P. aurifer (Saussure) nests located within cracks in the dried soil of the Santa Monica Mountains of southern California, USA.Received 8 April 2003; revised 10 October 2003; accepted 31 October 2003.     

Predatory wasps learn to overcome the shelter defences of their larval prey

Larvae of Epargyreus clarus (Hesperiidae), the silver-spotted skipper, inhabit leaf-and-silk shelters that they construct on their leguminous host plants. In the field, Polistes spp. (Vespidae) wasps land on the shelters, quickly extracting and killing the larvae within. In marked contrast, wasps that emerge from field-collected colonies maintained in the laboratory visit and examine leaflets bearing sheltered caterpillars, but only rarely do they extract and kill the sheltered larvae. To examine whether learning is involved in the development of the ability of Polistes wasps to forage successfully on sheltered E. clarus larvae, we tested the responses of P. fuscatus and P. dominulus wasps to sheltered E. clarus larvae before and after their exposure to unsheltered larvae that were visible either on an opened host-leaf shelter (P. fuscatus and P. dominulus) or on a nonhost leaf in the absence of a shelter (P. fuscatus). After killing and processing an unsheltered larva that was visible on an opened leaf shelter, a majority of foragers subsequently extracted and killed larvae from closed shelters. Wasps that killed and processed an unsheltered larva on a nonhost leaf, on the other hand, generally did not later open shelters. Thus, it seems that experience with an exposed larva in the context of its shelter is necessary for a wasp to be able to prey on sheltered larvae. We conclude that the wasps must learn to associate the taste of the larva with shelter-related cues, such as presence of leaf damage and silk. In nature, this initial exposure may occur when the larva is visible in or near its shelter, perhaps when feeding or constructing a new shelter. Learning opportunities will thus depend on larval density. Our results show that invertebrate predators can learn to overcome their prey's defences, and are therefore able to make use of previously inaccessible prey.     

Diversity, Distribution, and Ancient Taxonomic Relationships Within the TIR and Non-TIR NBS-LRR Resistance Gene Subfamilies

Phylogenetic relationships among the NBS-LRR (nucleotide binding site–leucine-rich repeat) resistance gene homologues (RGHs) from 30 genera and nine families were evaluated relative to phylogenies for these taxa. More than 800 NBS-LRR RGHs were analyzed, primarily from Fabaceae, Brassicaceae, Poaceae, and Solanaceae species, but also from representatives of other angiosperm and gymnosperm families. Parsimony, maximum likelihood, and distance methods were used to classify these RGHs relative to previously observed gene subfamilies as well as within more closely related sequence clades. Grouping sequences using a distance cutoff of 250 PAM units (point accepted mutations per 100 residues) identified at least five ancient sequence clades with representatives from several plant families: the previously observed TIR gene subfamily and a minimum of four deep splits within the non-TIR gene subfamily. The deep splits in the non-TIR subfamily are also reflected in comparisons of amino acid substitution rates in various species and in ratios of nonsynonymous-to-synonymous nucleotide substitution rates (K _A/K _S values) in Arabidopsis thaliana. Lower K _A/K _S values in the TIR than the non-TIR sequences suggest greater functional constraints in the TIR subfamily. At least three of the five identified ancient clades appear to predate the angiosperm–gymnosperm radiation. Monocot sequences are absent from the TIR subfamily, as observed in previous studies. In both subfamilies, clades with sequences separated by approximately 150 PAM units are family but not genus specific, providing a rough measure of minimum dates for the first diversification event within these clades. Within any one clade, particular taxa may be dramatically over- or underrepresented, suggesting preferential expansions or losses of certain RGH types within particular taxa and suggesting that no one species will provide models for all major sequence types in other taxa. Received: 13 June 2001 / Accepted: 22 October 2001     

Phylogenetic analysis of the Argonaute protein family in platyhelminths

Argonaute proteins (AGOs) are mediators of gene silencing via recruitment of small regulatory RNAs to induce translational regression or degradation of targeted molecules. Platyhelminths have been reported to express microRNAs but the diversity of AGOs in the phylum has not been explored. Phylogenetic relationships of members of this protein family were studied using data from six platyhelminth genomes. Phylogenetic analysis showed that all cestode and trematode AGOs, along with some triclad planarian AGOs, were grouped into the Ago subfamily and its novel sister clade, here referred to as Cluster 1. These were very distant from Piwi and Class 3 subfamilies. By contrast, a number of planarian Piwi-like AGOs formed a novel sister clade to the Piwi subfamily. Extensive sequence searching revealed the presence of an additional locus for AGO2 in the cestode Echinococcus granulosus and exon expansion in this species and E. multilocularis. The current study suggests the absence of the Piwi subfamily and Class 3 AGOs in cestodes and trematodes and the Piwi-like AGO expansion in a free-living triclad planarian and the occurrence of exon expansion prior to or during the evolution of the most-recent common ancestor of the Echinococcus species studied.     

基因组水平蝙蝠嗅觉受体基因的分子进化

翼手目动物（俗称蝙蝠）的食性分化显著,不同食性的蝙蝠具有显著不同的嗅球大小。为了研究嗅觉是否影响了蝙蝠食性的进化,我们利用网上已公布的10种蝙蝠基因组的数据,通过同源比对的方法鉴定出所有的嗅觉受体基因,并进行嗅觉受体基因亚家族的分类,进而比较嗅觉受体基因亚家族的数目差异。结果显示,蝙蝠的嗅觉受体基因与其它哺乳动物一样,都可以分为13个单系起源的亚家族;在Yinpterochiroptera亚目中,OR1/3/7、OR2/13、OR5/8/9等3个嗅觉受体亚家族在食果蝙蝠中均发生了不同程度的扩张,基因数目显著地多于食虫蝙蝠,提示嗅觉在食果蝙蝠取食过程中具有重要的作用。因此,本研究在基因组水平上重现了蝙蝠嗅觉受体基因的进化历史,揭示了3个嗅觉受体基因亚家族的功能可能与食果蝙蝠的食性相关,突出了嗅觉对动物食性的重要作用.     

Odorant and Gustatory Receptors in the Tsetse Fly Glossina morsitans morsitans

Tsetse flies use olfactory and gustatory responses, through odorant and gustatory receptors (ORs and GRs), to interact with their environment. Glossina morsitans morsitans genome ORs and GRs were annotated using homologs of these genes in Drosophila melanogaster and an ab initio approach based on OR and GR specific motifs in G. m. morsitans gene models coupled to gene ontology (GO). Phylogenetic relationships among the ORs or GRs and the homologs were determined using Maximum Likelihood estimates. Relative expression levels among the G. m. morsitans ORs or GRs were established using RNA-seq data derived from adult female fly. Overall, 46 and 14 putative G. m. morsitans ORs and GRs respectively were recovered. These were reduced by 12 and 59 ORs and GRs respectively compared to D. melanogaster. Six of the ORs were homologous to a single D. melanogaster OR (DmOr67d) associated with mating deterrence in females. Sweet taste GRs, present in all the other Diptera, were not recovered in G. m. morsitans. The GRs associated with detection of CO₂ were conserved in G. m. morsitans relative to D. melanogaster. RNA-sequence data analysis revealed expression of GmmOR15 locus represented over 90% of expression profiles for the ORs. The G. m. morsitans ORs or GRs were phylogenetically closer to those in D. melanogaster than to other insects assessed. We found the chemoreceptor repertoire in G. m. morsitans smaller than other Diptera, and we postulate that this may be related to the restricted diet of blood-meal for both sexes of tsetse flies. However, the clade of some specific receptors has been expanded, indicative of their potential importance in chemoreception in the tsetse.     

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm—lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits—extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.     

The Molecular Phylogeny of a Nematode-Specific Clade of Heterotrimeric G-Protein α-Subunit Genes

In animal olfactory systems, odorant molecules are detected by olfactory receptors (ORs). ORs are part of the G-protein-coupled receptor (GPCR) superfamily. Heterotrimeric guanine nucleotide binding G-proteins (G-proteins) relay signals from GPCRs to intracellular effectors. G-proteins are comprised of three peptides. The G-protein α subunit confers functional specificity to G-proteins. Vertebrate and insect Gα-subunit genes are divided into four subfamilies based on functional and sequence attributes. The nematode Caenorhabditis elegans contains 21 Gα genes, 14 of which are exclusively expressed in sensory neurons. Most individual mammalian cells express multiple distinct GPCR gene products, however, individual mammalian and insect olfactory neurons express only one functional odorant OR. By contrast C. elegans expresses multiple ORs and multiple Gα subunits within each olfactory neuron. Here we show that, in addition to having at least one member of each of the four mammalian Gα gene classes, C. elegans and other nematodes also possess two lineage-specific Gα gene expansions, homologues of which are not found in any other organisms examined. We hypothesize that these novel nematode-specific Gα genes increase the functional complexity of individual chemosensory neurons, enabling them to integrate odor signals from the multiple distinct ORs expressed on their membranes. This neuronal gene expansion most likely occurred in nematodes to enable them to compensate for the small number of chemosensory cells and the limited emphasis on cephalization during nematode evolution. [Reviewing Editor: Dr. John Oakeshott] Damien M. O’Halloran and David A. Fitzpatrick contributed equally to this work.     

Genomic organization and evolution of olfactory receptors and trace amine-associated receptors in channel catfish,Ictalurus punctatus

BackgroundChannel catfish (Ictalurus punctatus) live in turbid waters with limited visibility to chase prey within a certain distance. This can be compensated through detecting specific water-soluble substances by the olfactory receptors (ORs) and trace amine associated receptors (TAARs) expressed on the olfactory epithelium.MethodsWe identified the OR and TAAR repertoires in channel catfish, and characterized the genomic organizations of these two gene families by data mining available genomic resources.ResultsA total of 47 putative OR genes and 36 putative TAAR genes were identified in the channel catfish genome, including 27 functional OR genes and 28 functional TAAR genes. Phylogenetic and orthogroup analyses were conducted to illustrate the evolutionary dynamics of the vertebrate ORs and TAARs. Collinear analysis revealed the presence of two conserved orthologous blocks that contain OR genes between the catfish genome and zebrafish genome. The complete loss of a conserved motif in fish OR family H may contribute to the divergence of family H from other families. The dN/dS analysis indicated that the highest degree of selection pressure was imposed on TAAR subfamily 14 among all fish ORs and TAARs.ConclusionsThe present study provides understanding of the evolutionary dynamics of the two gene families (OR and TAAR) associated with olfaction in channel catfish.General significanceThis is the first systematic study of ORs and TAARs in catfish, which could provide valuable genomic resources for further investigation of olfactory mechanisms in teleost fish.     

Displacement and replacement in real time: <Emphasis Type="Italic">Polistes dominula</Emphasis>’<Emphasis Type="Italic">s</Emphasis> impact on <Emphasis Type="Italic">P. fuscatus</Emphasis> in the northeastern U.S.

Two major challenges in studying the impacts of exotic invasive species on native species are identifying mechanisms of displacement and replacement and the lack of long-term population studies in these systems. A solution for the first is to study invasive and native congeners that occupy the same niche. A solution for the second is to study many populations for one year instead of one population for many years. We studied the invasion biology of the invasive European paper wasp Polistes dominula and its native congener the Northern paper wasp P. fuscatus, two species which compete for similar resources. We tracked the demography of the two wasps at sites in the northeastern United States. We found that the survival of P. dominula to the reproductive period in August was three times that of P. fuscatus, across all sites. The reproductive output of P. fuscatus declined in direct proportion to the percentage of P. dominula nests at the site. P. fuscatus nests at uninvaded sites had three times the nest cells of those at the most invaded sites. These findings suggest a positive feedback cycle in the establishment of P. dominula, in which the invasive wasp drives population declines in the native that in turn allow P. dominula to further establish. This system provides an example of a possible extinction vortex caused by competitive exclusion of a species by its invasive congener.     

Structure-Function Analysis of a Broad Specificity Populus trichocarpa Endo-β-glucanase Reveals an Evolutionary Link between Bacterial Licheninases and Plant XTH Gene Products

The large xyloglucan endotransglycosylase/hydrolase (XTH) gene family continues to be the focus of much attention in studies of plant cell wall morphogenesis due to the unique catalytic functions of the enzymes it encodes. The XTH gene products compose a subfamily of glycoside hydrolase family 16 (GH16), which also comprises a broad range of microbial endoglucanases and endogalactanases, as well as yeast cell wall chitin/β-glucan transglycosylases. Previous whole-family phylogenetic analyses have suggested that the closest relatives to the XTH gene products are the bacterial licheninases (EC 3.2.1.73), which specifically hydrolyze linear mixed linkage β(1→3)/β(1→4)-glucans. In addition to their specificity for the highly branched xyloglucan polysaccharide, XTH gene products are distinguished from the licheninases and other GH16 enzyme subfamilies by significant active site loop alterations and a large C-terminal extension. Given these differences, the molecular evolution of the XTH gene products in GH16 has remained enigmatic. Here, we present the biochemical and structural analysis of a unique, mixed function endoglucanase from black cottonwood (Populus trichocarpa), which reveals a small, newly recognized subfamily of GH16 members intermediate between the bacterial licheninases and plant XTH gene products. We postulate that this clade comprises an important link in the evolution of the large plant XTH gene families from a putative microbial ancestor. As such, this analysis provides new insights into the diversification of GH16 and further unites the apparently disparate members of this important family of proteins.     

Inventory of the cichlid olfactory receptor gene repertoires: identification of olfactory genes with more than one coding exon

Background

To help understand the molecular mechanisms underlying the remarkable phenotypic diversity displayed by cichlids, the genome sequences of O. niloticus, P. nyererei, H. burtoni, N. brichardi and M. zebra were recently determined. Here, we present the contents of the olfactory receptor (OR) repertoires in the genomes of these five fishes.

Results

We performed an exhaustive TBLASTN search of the five cichlid genomes to identify their OR repertoires as completely as possible. We used as bait a set of ORs described in the literature. The cichlid repertoires thereby extracted contained large numbers of complete genes (O. niloticus 158; H. burtoni 90; M. zebra 102; N. brichardi 69; P. nyererei 88), a small numbers of pseudogenes and many “edge genes” corresponding to incomplete genes located at the ends of contigs. A phylogenetic tree was constructed and showed these repertoires include a large number of families and subfamilies. It also allowed the identification of a large number of OR analogues between cichlids with very high amino-acid identity (≥99%). Nearly 9% of the full-length cichlid OR genes are composed of several coding exons. This is very unusual for vertebrate OR genes. Nevertheless, the evidence is strong, and includes the donor and acceptor splice junction sequences; also, the positions of these genes in the phylogenetic tree indicate that they constitute subfamilies well apart from non-OR G protein-coupled receptor families.

Conclusions

Cichlid OR repertoires are made up of a larger number of genes and fewer pseudogenes than those in other teleosts except zebrafish. These ORs share all identified properties common to all fish ORs; however, the large number of families and subfamilies, each containing few ORs implies that they have evolved more rapidly. This high level of OR diversity is consistent with the substantial phenotypic diversity that characterizes cichlids.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-586) contains supplementary material, which is available to authorized users.